小川の卒論日誌
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
*当面の課題 [#z290dd6a]
**20/5/11(後藤) [#q3a934a0]
-弾塑性解析ができるか確認(小川さん、佐藤さん、吉田さんで...
-[[SALOME-Mecaの使用法解説:http://opencae.gifu-nct.ac.jp/...
-固定端の引っ張り側の縁部の応力が、降伏点(105MPa)に達する...
-荷重（縦軸）、先端のたわみ（横軸）でグラフを描き、固定端...
-以上が確認できたら、片持ち梁の弾塑性解析のhdfファイル（...
**20/5/18 [#c33ced60]
-(x,y,z)=(10,20,100)のboxを作成し、メッシュサイズは１にし...
-メッシュを作成して保存し、サーバーを閉じようとしたら、固...
-強制終了して/home/salome_meca/appli_V2017/salomeで起動し...
-起動した時に端末にメッセージが表示された。
Warning: this type (Study,objref) already exists, it will...
found in 20.5 seconds
Start SALOME, elapsed time : 21.8 seconds
[#p5f7e7ba]
Warning: HYBRIDPLUGIN GUI resources are not found. [#v...
HYBRIDPLUGIN GUI will not be available. [#n38190e5]
[#jf596eb1]
[#z664fb54]
Warning: BLSURFPLUGIN GUI resources are not found. [#r59...
BLSURFPLUGIN GUI will not be available. [#h0f5c81d]
[#x05d1df1]
[#q3e1609c]
Warning: HexoticPLUGIN GUI resources are not found. [...
HexoticPLUGIN GUI will not be available. [#m001e348]
[#e1fbd5a5]
[#r94e0fd1]
Warning: NETGENPLUGIN GUI resources are not found. [#...
NETGENPLUGIN GUI will not be available. [#ac3e85e8]
[#tfc991fb]
[#i0dc08d4]
Warning: GHS3DPLUGIN GUI resources are not found. [#x...
GHS3DPLUGIN GUI will not be available. [#q2a11c06]
[#yba6c902]
**20/5/21 [#j953432e]
-AsterSutudyのDEFI_FONCTIONが見つからない
**20/5/28 [#qf12b966]
-6-5のCALC_ELEM、CALC_NOが見つからず、調べたところ、「バ...
-解析を実行したが、Outputのところに不備があり、赤色になっ...
-b_not_resueやb_noilなど、b_で始まるコマンドがほとんど見...
**20/6/5 [#vc356925]
-近藤さんのcommファイルやhttps://qiita.com/Jun_Tatsuno/it...
**20/6/12 [#v11bdfc0]
-BC and Loadのところに不備があり、直したら青色になった。...
**20/6/19 [#xc265596]
-ParaVisでTimeを動かし、変位を確認することができた。来週...
**20/6/26 [#l9bbef6e]
-変位をまとめてグラフを作成している最中です。
**20/6/30 [#l4c1f44f]
-縦軸が応力、横軸がひずみのグラフを作成した。グラフの傾き...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/0630.png
**20/7/10 [#m5d66eaa]
-functions and listsのlistrのDEFI_LIST_REELのPASを0.1から...
-Timeの10分割を変える方法が分からなかったので、PASを0.01...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/dan.png
**20/7/15 [#t5a1b6e3]
-根元の応力が105MPaになるときの荷重を計算したところ、荷重...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/da.png
**20/7/21 [#wb1b18d7]
-前回は応力とひずみの向きを考慮していなかったので、それを...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/danso.png
-上のグラフは先端のひずみで計算していたので、根元のひずみ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/dans.png
**20/7/24 [#f799d779]
-梁を軸方向に引っ張ったときの解析を行った。引張応力が105M...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/hippari.png
**20/8/7 [#m595aa64]
-今までは面に載荷していたが，集中荷重なので線載荷にしてみ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/sensaika....
**20/10/4 [#h52cbe94]
-Geometryで床版、主プレート、補助プレート、補助金具、支柱...
**20/10/7 [#t0d47e43]
-鋼材のメッシュサイズを決めるために支柱のみの解析をした．
-支柱を長さ600mmの円筒型の片持ち梁と仮定して理論値を計算...
-理論値計算と同じ条件で解析した．ポアソン比は0.3とした．...
,メッシュサイズ,変位
,6,0.000154
,7,0.00015381
,8,0.00015373
,9,0.00015379
,10,0.00015365
,11,0.00015372
,12,0.0001538
-メッシュサイズが1の時は他のサイズと同じ条件でも計算でき...
-この結果から，メッシュサイズが10の時に理論値に一番近くな...
**20/10/8 [#j78223ac]
-CLTを長さ2000mmの単純梁として手計算をした．
-荷重は1N，断面二次モーメントは281250000mm4とした．
-CLTの等級は異等級構成Mx60A(中村先生)
-CLTメモの等級ごとのヤング率の求め方の表3より，Mx60の5層5...
-計算したところ，変位は0.000423mmとなった．
**20/10/9 [#h3e3f621]
-鋼材のヤング率は正しくは206GPaだったので，その数値で理論...
**20/10/12 [#q89d4472]
-鋼材の正しいヤング率で解析した．結果は図のようになった．
,メッシュサイズ,変位
,6,0.000157
,7,0.0001568
,8,0.0001567
,9,0.0001568
,10,0.0001566
,11,0.0001567
,12,0.0001568
-メッシュサイズが10の時に理論値に一番近くなった．相対誤差...
-直交集成板の日本農林規格によると，「直交集成板を構成する...
-Mx 60は外層ラミナの強度等級がM60A、内層ラミナの強度等級...
-M60A,M30Aのヤング率の基準値は直交集成板の日本農林規格に...
**20/10/14 [#n55bd6e2]
-M60A,M30Aのヤング率をそれぞれ6.0GPa,3.0GPaとしてMx60のヤ...
-Shear Analogy Methodという方法で計算した．有効曲げ剛性(E...
-この有効曲げ剛性をCLT全体の断面二次モーメントで割ったと...
-有効せん断剛性も求めた．
-繊維方向のヤング率と各方向のせん断弾性係数の関係式を用い...
-求めた剛性を使ってCLTの変位を求めた．荷重は梁の中心に1N...
**20/10/16 [#f51c805e]
-メッシュサイズを更に増やして支柱の解析をした．結果は図と...
,メッシュサイズ,要素数,変位
,5,145813,0.000157025
,6,144052,0.0001569962
,7,38244,0.0001567935
,8,37009,0.0001567133
,9,36804,0.0001567753
,10,36116,0.000156634
,11,36035,0.0001567057
,12,35293,0.0001567914
,15,9570,0.0001533965
,20,9148,0.0001537212
,25,8747,0.0001540633
,30,2972,0.0001336225
,35,2364,0.0001379303
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ko.png
-このグラフから，要素数を増やすと0.000157当たりに収束して...
**20/10/19 [#n4f0e802]
-CLTの方の解析をした．メッシュサイズを20にしたところ，要...
-メモリの上限を超えているというメッセージが出た．
Error during the dynamic allocation. It was not possible...
! a memory zone length 1169 Mo, one exceeds the maxim...
! built-in with 1493 Mo and one occupies already 580 ...
! The last operation of memory deallocation made it p...
!-----------------------------------------------------...
**20/10/20 [#p7567c5d]
-要素数を少なくしたら，メッシュを作成することが出来なかっ...
-NgException:Stop meshing since too meny attemptsというエ...
-検索してみたがよく分からなかった．
https://www.salome-platform.org/forum/forum_10/156488534
-及川さんに聞いてみたところ，最小のメッシュの大きさが穴の...
**20/10/21 [#haab6b72]
-横2ミリ，縦1ミリの穴を横20ミリ，縦10ミリの大きさにしてみ...
-解析結果は図とグラフのようになった．
,メッシュサイズ,要素数,変位
,14,581948,0.0001126382
,15,443110,0.0001150793
,20,335445,0.0001133839
,25,255690,0.0001076753
,27,131296,0.0001036925
,30,83221,0.000097703671
,40,63683,0.000095167519
,50,54909,0.000093780481
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/syouban.png
-CLTを等方性としてたわみを計算しなおしたら，0.0001247878...
**20/10/28 [#i58662a3]
-支柱にブラケットをつけて計算し直した．
-ブラケットは一辺90ミリの立方体とした．等方性としたのでヤ...
-支柱の高さ588ミリの部分にグループを作成し，そこの変位を...
,メッシュサイズ,要素数,変位
,3,650732,0.0001525962
,4,253843,0.0001528185
,5,226615,0.000152853
,10,57390,0.0001527818
,20,13573,0.0001504562
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/kouz.png
**20/10/29 [#h7d7bf8f]
-プレートと床版が一体化していない状態で解析する必要がある...
-AsterStudyはMATE=AFFE_MATERIAUまでやった．
**20/10/30 [#ad3ee1e9]
-支柱のメッシュをもっと粗くして計算しなおした．結果は表と...
-グラフを見やすくするためにメッシュサイズ3のときの変位は...
,メッシュサイズ,要素数,変位
,3,650732,0.0001525962
,4,253843,0.0001528185
,5,226615,0.000152853
,10,57390,0.0001527818
,20,13573,0.0001504562
,30,5931,0.0001270705
,40,3982,0.0001377993
,50,3411,0.0001174439
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/kouz1.png
**20/11/2 [#n23f4312]
-接触解析の基本(1)の続き
-AsterStudyのCONTACTという項目が見つからない．下のサイト...
http://salome-meca.cocolog-nifty.com/blog/2012/03/60-1-2f...
-接触解析の基本(2)の荷重拘束
-CONTACTの代わりにDEFI-CONTACTを接触の設定として用いた．
-接触の部分はSalomeMecaの使用法解説の接触解析の基本(2)を...
-黄色マークで下のようなメッセージが出た．
<S> Exception user raised but not interceptee. !
! The bases are fermees. !
! Type of the exception: ArretCPUError !
! [('?', (), (), ())]
-このメッセージの意味はよく分からなかった． ...
-Salome接触解析メモを見たところ，部材ごとでメッシュを切ら...
**20/11/3 [#p73e8f5a]
-床版のメッシュを粗くして計算しなおした．結果は表とグラフ...
,メッシュサイズ,要素数,変位
,14,581948,0.0001126382
,15,443110,0.0001150793
,20,335445,0.0001133839
,25,255690,0.0001076753
,27,131296,0.0001036925
,30,83221,0.000097703671
,40,63683,0.000095167519
,50,54909,0.000093780481
,60,21413,0.000076680678
,70,18683,0.000073994244
,80,16084,0.000072255257
,100,14087,0.000070051691
,150,5070,0.000042707938
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/syouban1....
-計算し直す前よりは収束している様子を再現することが出来た．
-一応半解析の方もやってみた．床版を半分に切断し，切断面を...
-切断面にも穴を開けて載荷線を作るという方法でやってみた．...
<S> Exception user raised but not interceptee. ...
! The bases are fermees. ...
! Type of the exception: error ...
! ...
! Problème : la matrice est singulière ou presque sing...
! Lors de la factorisation de la matrice, on a renco...
! (pivot nul ou presque nul) à la ligne 37413 qui co...
! donné ci-dessus. ...
! Le nombre de décimales "perdues" pour cette ligne ...
! Ce nombre de décimales perdues est lié au degré de...
! Plus il est grand, plus le système est singulier. ...
! Quand on a perdu plus de 8 décimales, on estime qu...
! précision ...
! des nombres réels (qui ont 15 décimales environ). ...
! ...
! Risques et conseils : ...
! * Si la ligne correspond a un degré de liberté ph...
! probablement d'un mouvement ...
! de corps rigide mal bloqué. ...
! Vérifiez les conditions aux limites. ...
! Si vous faites du contact, il ne faut pas que l...
! par le contact. ...
! Vérifiez également les caractéristiques matéria...
! ...
! * Si la ligne correspond a un degré de liberté de...
! doute d'une condition ...
! limite redondante. ...
! En particulier, il se peut que la relation liné...
! des conditions de contact. ...
! Peut-être devriez vous exclure certains noeuds ...
! (mots clés SANS_NOEUD et SANS_GROUP_NO). ...
! ...
! * Si le solveur utilisé est LDLT ou MULT_FRONT, v...
! solveur MUMPS ...
! car celui-ci est le seul à pouvoir factoriser l...
! définies positives. ...
! ...
! * Parfois, en parallèle, le critère de détection ...
! trop pessimiste ! Il reste néanmoins souvent ...
! possible de faire passer le calcul complet en r...
! (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou ...
! en le débranchant (valeur du mot-clé NPREC=-1) ...
! sur moins de processeurs. ...
! ...
! * Il se peut aussi que ce phénomène soit tout à f...
! fissure passe ...
! très près d'un noeud. ...
! Si le nombre de décimales perdues n'est pas tro...
! vous pouvez relancer le calcul en augmentant le...
! autorisé : ...
! mot-clé NPREC du mot clé facteur SOLVEUR. ...
! Sinon, contactez l'équipe de développement.
-載荷線のxy方向の固定を解除してみたが，同じメッセージが出...
**20/11/6 [#p89342ae]
-ティモシェンコ梁で理論値の計算をした．
-せん断補正係数は，ポアソン比の0.3を代入して0.5306122449...
-せん断補正係数とせん断弾性係数(79GPa)を代入して計算した...
-支柱の長さを10倍ほどにして解析した．メッシュサイズを10,4...
-10のときに保存しないまま強制終了したのでそのときのhdfフ...
**20/11/8 [#ba52efdb]
-10倍だとできなかったので、とりあえず3倍でやってみること...
,メッシュサイズ,要素数,変位
,10,110765,0.00014746275
,20,27475,0.00014276444
,21,27419,0.000142449056
,23,19402,0.000142579938
,25,17792,0.0001425346
,30,10050,0.00012982425
,40,6912,0.000105412889
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/3sityuu.p...
**20/11/11 [#h10c4fe2]
-接触解析の基本(1)(摩擦なし)の続き
-MODI_MAILLAGE: メッシュの操作
-接触解析メモにMODI_MAILLAGEがなかったのでその項目を消し...
<ArretCPUError> <MECANONLINE9_1> ...
! ...
! Arrêt par manque de temps CPU pendant les itérati...
! - Temps moyen par itération de Newton : 81.619...
! - Temps restant : 130.51...
! La base globale est sauvegardée. Elle contient le...
! ...
! Conseil : ...
! - Augmentez le temps CPU. ...
!-------------------------------------
-このようなエラーメッセージが出た．CPU時間を増やすという...
-timeを15に戻し，STAT_NON_LINE(非線形解析の設定)を弾塑性...
**20/11/12 [#a92df122]
-今までやっていた接触解析のモデルはSTEPファイルを読みこん...
-荷重拘束の接触解析（摩擦あり）のペナルティ法
-Geometry：baceとtopの対称面はそれぞれでグループ化した．
-Mesh：メッシュ→アセンブリをビルドでコンパウンドを作成．G...
-AsterStudy：2.5次要素への変換は、CREA_MAILLAGEコマンドを...
**20/11/13 [#oa5b859b]
-床版(単純梁)の解析
-ふたつ目の固定線が全方向拘束されていたので，kotei2のy方...
-外す前とほとんど変化はなかった．
-理論値に一番近い値のときの相対誤差は-7.78%だった．
,メッシュサイズ,要素数,変位
,14,581948,0.0001126387
,15,443110,0.0001150798
,20,335445,0.000113384
,25,255690,0.000107676
,27,131296,0.0001036928
,30,83221,0.000097704057
,40,63683,0.000095167927
,50,54909,0.000093780519
,60,21413,0.000076680739
,70,18683,0.000073994344
,80,16084,0.000072255329
,100,14087,0.0000700518
,150,5070,0.000042708863
,200,3587,0.00003988833
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/1syouban....
-二次要素への変換方法
-作成された二次要素を右クリック→二次/一次要素への変換→二...
**20/11/16 [#w8a6ceff]
-二次要素への変換
-メッシュサイズ40の支柱モデルを二次要素に変換したら警告メ...
「ノードが （マジェンタ色で図示） 環境によって歪んだ要素...
-メッシュサイズ30のときも同じメッセージが出たが，25のとき...
-2次要素で解析した結果は表とグラフのようになった．(30と40...
,メッシュサイズ,要素数,変位
,10,110765,0.000148931
,20,27475,0.0001490144
,21,27419,0.0001489303
,23,19402,0.000148931
,25,17792,0.0001489472
,30,10050,0.00012982425
,40,6912,0.000105412889
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/2zisityuu...
-一番理論値に近い値の相対誤差は4.39%になった．
**20/11/17 [#q38e7f11]
-床版の方も二次要素で解析した．こちらでは支柱のようなエラ...
-しかし，要素数が大きすぎると計算時にエラーが起こった．
<S> Exception user raised but not interceptee. ...
! The bases are fermees. ...
! Type of the exception: error ...
! ...
! Error during the dynamic allocation. It was not pos...
! a memory zone length 41 Mo, one exceeds the maximum...
! built-in with 1499 Mo and one occupies already 1469...
! The last operation of memory deallocation made it p...
-メッシュサイズ25までは計算できたので，その結果を表とグラ...
,メッシュサイズ,要素数,変位
,200,3587,0.0001210055
,150,5070,0.0001209961
,100,14087,0.000121258
,80,16084,0.0001212511
,70,18683,0.0001212428
,60,21413,0.0001212307
,50,54909,0.0001212352
,40,63683,0.000121233
,30,83221,0.0001212595
,27,131296,0.000121299
,25,255690,0.0001213208
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/2zisyouba...
-グラフでは理論値と離れているように見えるが，一番近い値の...
-Geometryで作った物体を他のファイルに移動させる方法
-移動させたい物体を選択する→ファイル→エクスポート→STEPエ...
-ファイル→インポート→STEPインポート→保存したものを選んで...
**20/11/18 [#s2a5984b]
-軸直交方向のヤング率(床版)を求めた．Shear Analogy Method...
-軸方向のヤング率を4.8704GPaとしていたが，正しくは4.7488G...
-メッシュサイズ25で計算してみたら，より理論値に近くなった...
,メッシュサイズ,要素数,変位
,25,255690,0.0001244172
,27,131296,0.0001243948
,40,63683,0.0001243271
,60,21413,0.0001243248
,100,14087,0.0001243528
,20,3587,0.0001240943
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/2naosi.png
-理論値に一番近いところの相対誤差は-0.3%だった．遠いとこ...
**20/11/19 [#y9db4711]
-支柱が降伏するときの荷重
-I=2343194.144(mm4),σ=238(MPa),l=588(mm),y=57.15(mm)のと...
-16.6kN=1.69tfより，支柱が降伏するときの荷重は1.69tf
-先日の実験データのグラフ化
-荷重と変位の関係をLibreOfficeCalc内で作成した．変位はCLT...
-gnuplotでcsvファイルのグラフを作成する
-区切りをスペースからカンマに変更：set datafile separator...
**20/11/20 [#s2779479]
-実験の荷重と変位，ひずみの関係のグラフを作成
-荷重と変位の関係
-手計算では1.69tfで降伏するという結果になったが，グラフは...
-荷重とひずみの関係(支柱の根元)
-内側は荷重1.8tf当たりから急激にひずみが増加している．外...
-東側上面
-圧縮側にひずみが増加すると考えていたch11,12,13は引張側に...
-東側下面
-ch35〜38は圧縮側に増加していたが，しっかりとした線形では...
-真ん中
-少し線形から外れる点はあったが，特に問題はなかった．
-西側上面
-引張側は特に問題なかったが，圧縮側は0のままのものがあっ...
-西側下面
-ch49が0からほぼ動かない．ch47,48,50が不規則な動き方をし...
**20/11/25 [#cf74f64c]
-とりあえずLibreOfficeでグラフを作成した．荷重と変位，支...
-撮った写真はアプリでJPEG形式に変換した．
**20/11/26 [#l69ade34]
-gnuplotで荷重とひずみ(支柱根元)の関係のグラフを作成した...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/nem...
-CLTの縦ヤング率の手計算をしたら，150MPaになった．それを...
-確かめるため，解析も行った．繊維方向のヤング率，繊維直交...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/tat...
**20/11/27 [#la3d8e74]
-10倍の支柱の解析を行った．メッシュを作るとき，いつものや...
-しかし，メッシュサイズを変えても要素数が全く変わらない．...
**20/11/30 [#zcbef286]
-10倍の支柱の解析をした．Maxサイズを小さくすると要素数が...
,メッシュサイズ,要素数,変位
,Max100 min20,42192,1.3718466667
,M100 m10,67880,1.4224408333
,M100 m4.5,139661,1.4488306757
,M10 m4.5,357844,1.4893606757
,M8 m4.5,493323,1.4849333784
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/10s...
-全体モデル
-全部くっついている状態，CLTの切断面を固定，ブラケットに...
-ポアソン比をどう決めればいいかまだわからない．
**20/12/1 [#i663c898]
-写真をk2の/home/gakusei/zikken/cltにアップロードしました...
**20/12/2 [#ube1770d]
-ポアソン比の計算をした．構造関係メモの直交異方性材料のポ...
-1/25じゃなくて1/30にしたのは，CLTメモ内のShear Analogy M...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/poi...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/poi...
**20/12/3 [#h42dd35c]
-せん断弾性係数の計算
-全体モデルでの解析
-外桁の固定面のグループをメッシュに反映させようとしたらエ...
**20/12/4 [#s6814cf3]
-せん断弾性係数は、どの方向も繊維方向のヤング率の1/15
-エラーが起こる前からやり直してみたが、何も問題なく反映さ...
**20/12/7 [#y1dc24e8]
-弾性の全体モデルの解析を行った．荷重は降伏する直前の2.5t...
,全体モデル(弾性)
,荷重(t),変位(mm)
,0.5,1.329206
,1.0,2.658416
,1.5,3.98762
,2.0,5.31683
,2.5,6.646036
,実験値
,荷重(t),変位(mm)
,0.5,2.625
,1.0,5.5
,1.5,8.4
,2.0,11.55
,2.5,14.7
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/mod...
-以下のように修正した．
-切断面の縦方向を動けるようにした．
-外桁部分にゴムをつけた．
-ゴムとH鋼の接触面を全面固定した．
-ゴムのヤング率：0.002GPa,ポアソン比：0.49とした．
-とりあえず荷重2.5ｔで解析したところ，変位は9.66262375mm...
-次に，2次要素に直してみると，変位は12.7262142857となった．
-2次要素ならメッシュが粗くても精度は変わらないだろうと考...
-それでも相対誤差は-13.4%ほどある．
**20/12/8 [#qaf23da4]
-2次要素のメッシュサイズMax100,min3.2で0.5tから0.5刻みで...
-実験値は11/24のもの
,荷重(t),変位(mm)
,0.5,2.5460526316
,1,5.0921063158
,1.5,7.6381589474
,2,10.1842136842
,2.5,12.7302684211
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/rub...
**20/12/9 [#r8407b72]
-ゴムのヤング率を2MPaから100MPaに直した．ゴムが固くなった...
-実験でジャッキのヘッド側が回転し，斜めに荷重をかけていた...
-床版の真ん中の変位が実験と解析値でかなり差があったので理...
**20/12/16 [#qd837b19]
-接触解析ここを参考にやってみる．
https://qiita.com/Jun_Tatsuno/items/6159d914abdc75a91262
-変位制御で解析することが出来たが，WarpByVectorを見ること...
→GenerateVectorsを押したら出来た．
-荷重制御に変えてみたら，ニュートン法の最大繰り返し回数で...
**21/2/9 [#t3ae3af3]
-https://www.kochi-kensoku.org/q-a%E3%82%B3%E3%83%BC%E3%8...
-ガードレールの衝突荷重について
**21/2/11 [#cc6cfb40]
-http://www.pwrc.or.jp/thesis_shouroku/thesis_pdf/0406-P0...
-ガードレールの衝突荷重について２
**21/2/12 [#kb8827b2]
-http://www.jsce.or.jp/committee/steel/subcomit-active/hp...
-降伏点，安全率について
*卒論日誌 [#ff45fcb7]
,日付,時間帯,作業時間,内容,立会
,10/11,16:00-17:00,1h,ガイダンス,後藤さん
,10/18,14:30-16:30,2h,創造工房実習,後藤さん
,10/21,13:00-16:00,3h,課題,
,10/21,14:30-16:30,2h,創造工房実習,後藤さん
,10/28,13:30-16:00,2.5h,課題,
,11/1,14:30-16:00,1.5h,創造工房実習,後藤さん
,11/6,11:30-12:30,1h,課題,
,11/8,14:30-16:00,1.5h,創造工房実習,後藤さん
,11/11,14:00-15:00,1h,課題,
,11/13,12:00-13:30,1h,課題,
,11/15,14:30-16:00,1.5h,創造工房実習,後藤さん
,11/20,12:00-12:30,0.5h,課題,
,11/22,14:30-16:00,1.5h,創造工房実習,後藤さん
,11/27,12:00-14:00,2h,課題,
,12/6,14:30-16:00,1.5h,創造工房実習,後藤さん
,12/11,11:00-13:30,2.5h,課題,
,12/13,14:30-16:00,1.5h,創造工房実習,後藤さん
,12/18,12:00-14:00,2h,課題,
,12/19,14:30-16:30,2h,課題,
,12/20,14:30-16:00,1.5h,創造工房実習,
,1/10,14:30-16:00,1.5h,創造工房実習,後藤さん
,1/15,14:30-15:30,1h,創造工房実習,後藤さん
,1/24,14:30-16:00,1.5h,創造工房実習,後藤さん
,1/29,18:00-19:30,1.5h,課題,
,1/29,21:00-22:00,1h,課題,
,1/31,14:30-16:00,1.5h,創造工房実習,後藤さん
,1/31,16:00-17:30,1.5h,課題,
,2/1,17:00-17:30,0.5h,課題,
,2/3,14:00-14:30,0.5h,課題,
,2/3,16:30-19:00,2.5h,課題,
,2/3,20:00-22:00,2h,課題,
,2/4,10:30-11:30,1h,課題,
,2/4,13:00-16:00,3h,課題,
,2/5,12:00-17:00,5h,課題,
,2/6,14:30-15:00,0.5h,課題,
,2/7,14:30-15:00,0.5h,創造工房実習,後藤さん
,2/13,13:45-14:45,1h,課題,
,2/14,14:30-15:30,1h,創造工房実習,後藤さん
,2/18,9:30-10:00,0.5h,春休み課題,
,2/26,14:00-16:30,2.5h,春休み課題,
,2/27,15:00-16:30,1.5h,春休み課題,
,3/2,16:45-17:45,1h,春休み課題,
,3/5,16:30-17:30,1h,春休み課題,
,3/7,16:00-18:30,2.5h,春休み課題,
,3/9,16:00-17:30,1.5h,春休み課題,
,3/11,15:30-17:30,2h,春休み課題,
,3/16,14:00-18:30,4.5h,春休み課題,
,3/17,13:00-16:00,3h,春休み課題,
,3/18,13:30-17:30,4h,春休み課題,
,3/19,14:00-17:30,3.5h,春休み課題,
,3/23,16:00-18:30,2.5h,春休み課題,
,3/24,15:00-17:30,2.5h,春休み課題,
,3/25,15:00-16:30,1.5h,春休み課題,
,3/26,16:00-17:30,1.5h,春休み課題,
*10/18の課題 [#xef89c31]
,メッシュ長さ,要素数,先端変位,相対誤差,計算者
,0.5,394121,6.57938,-1.4%,佐藤
,1,72278,6.443,-3.4%,吉田
,2,11817,5.7382525,-13.9%,梅田
,4,2862,4.9428,-25.9%,青山
,8,897,4.0411725,-39.4%,小川
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/ogawa/kat...
*10/25の課題 [#k79ee126]
,メッシュ長さ,要素数,先端変位,相対誤差,計算者
,0.5,394121,6.57938,-1.4%,佐藤
,0.7,130916,6.4781,-2.9%,佐藤
,0.8,72101,6.437,-3.5%,青山
,0.9,71718,6.4314,-3.5%,吉田
,1,72278,6.443,-3.4%,吉田
,1.2,65575,6.443,-3.4%,小川
,1.4,41096,6.3162,-5.3%,青山
,1.5,23417,6.1206,-8.2%,梅田
,1.8,11758,5.7370,-13.98%,梅田
,2,11817,5.7382525,-13.9%,梅田
,4,2862,4.9428,-25.9%,青山
,8,897,4.0411725,-39.4%,小川
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/ogawa/kat...
考察
メッシュ長さ1の時の要素数が0.9、0.8の時よりも多くなってい...
*11/1の課題 [#w19e33c2]
,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.5,322687,0.428774714,2.9%,小川
,0.8,199968,0.41646,-0.050%,吉田
,1,62360,0.416513545,-0.037%,梅田
,1.5,50359,0.411953,-1.134%,青山
,2,20526,0.3706311667,-11.05%,佐藤
理論値 0.416667
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/ogawa/tan...
考察 メッシュを細かくするほど精度が高くなるわけではない...
*11/8の課題 [#n612ec71]
,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.5,322687,0.5214958571,6.07%,小川
,0.8,199968,0.5145437143,4.65%,吉田
,0.9,87953,0.5102753846,3.78%,梅田
,1,62360,0.505768727273,2.87%,青山
,1.1,83278,0.5033813636,2.38%,小川
,1.2,59515,0.5053267,2.78%,吉田
,1.3,57455,0.5038435556,2.47%,青山
,1.5,50359,0.50010225,1.71%,梅田
,1.8,28053,0.4926668571,0.20%,佐藤
,2,20526,0.4765718333,-3.08%,佐藤
ティモシェンコ梁 0.4916667
理論値 0.416667
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/ogawa/tan...
*11/15の課題 [#ee35932d]
わかったこと
・自動メッシュの3次元を選択までの流れ
わからなかったこと
・スケール係数の欄がどこにあるか
・固定面の節点の選択のしかた
・エンティティの追加の欄がどこにあるか
*11/22の課題 [#uf06888c]
,拡大係数,要素数,変位,相対誤差,計算者
,0.3,57431,6.422525,-3.7%,全員
,0.4,23349,6.26413,-6.0%,佐藤
,0.5,13741,6.03154,-9.5%,小川
,0.6,7257,5.91187,-11.3%,梅田
,0.7,4245,5.5716725,-16.4%,吉田
,0.8,3657,5.5734,-16.4%,青山
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/ogawa/kat...
*12/6の課題 [#e860c736]
,拡大係数,要素数,変位,相対誤差,計算者
,0.8,122515,6.91749,2.6%,吉田
,1,65490,6.8877,2.2%,梅田
,1.2,36693,6.8519225,1.6%,梅田
,1.4,23341,6.8117175,1.0%,吉田
,1.6,15371,6.6336,-1.6%,小川
,1.8,12480,6.74258,0.014%,佐藤
,2,8589,6.66499,-0.011%,梅田
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/ogawa/kat...
節点数で荷重を割るときに精度を上げていれば変位が下がらず...
有効数字
終了行:
#contents
*当面の課題 [#z290dd6a]
**20/5/11(後藤) [#q3a934a0]
-弾塑性解析ができるか確認(小川さん、佐藤さん、吉田さんで...
-[[SALOME-Mecaの使用法解説:http://opencae.gifu-nct.ac.jp/...
-固定端の引っ張り側の縁部の応力が、降伏点(105MPa)に達する...
-荷重（縦軸）、先端のたわみ（横軸）でグラフを描き、固定端...
-以上が確認できたら、片持ち梁の弾塑性解析のhdfファイル（...
**20/5/18 [#c33ced60]
-(x,y,z)=(10,20,100)のboxを作成し、メッシュサイズは１にし...
-メッシュを作成して保存し、サーバーを閉じようとしたら、固...
-強制終了して/home/salome_meca/appli_V2017/salomeで起動し...
-起動した時に端末にメッセージが表示された。
Warning: this type (Study,objref) already exists, it will...
found in 20.5 seconds
Start SALOME, elapsed time : 21.8 seconds
[#p5f7e7ba]
Warning: HYBRIDPLUGIN GUI resources are not found. [#v...
HYBRIDPLUGIN GUI will not be available. [#n38190e5]
[#jf596eb1]
[#z664fb54]
Warning: BLSURFPLUGIN GUI resources are not found. [#r59...
BLSURFPLUGIN GUI will not be available. [#h0f5c81d]
[#x05d1df1]
[#q3e1609c]
Warning: HexoticPLUGIN GUI resources are not found. [...
HexoticPLUGIN GUI will not be available. [#m001e348]
[#e1fbd5a5]
[#r94e0fd1]
Warning: NETGENPLUGIN GUI resources are not found. [#...
NETGENPLUGIN GUI will not be available. [#ac3e85e8]
[#tfc991fb]
[#i0dc08d4]
Warning: GHS3DPLUGIN GUI resources are not found. [#x...
GHS3DPLUGIN GUI will not be available. [#q2a11c06]
[#yba6c902]
**20/5/21 [#j953432e]
-AsterSutudyのDEFI_FONCTIONが見つからない
**20/5/28 [#qf12b966]
-6-5のCALC_ELEM、CALC_NOが見つからず、調べたところ、「バ...
-解析を実行したが、Outputのところに不備があり、赤色になっ...
-b_not_resueやb_noilなど、b_で始まるコマンドがほとんど見...
**20/6/5 [#vc356925]
-近藤さんのcommファイルやhttps://qiita.com/Jun_Tatsuno/it...
**20/6/12 [#v11bdfc0]
-BC and Loadのところに不備があり、直したら青色になった。...
**20/6/19 [#xc265596]
-ParaVisでTimeを動かし、変位を確認することができた。来週...
**20/6/26 [#l9bbef6e]
-変位をまとめてグラフを作成している最中です。
**20/6/30 [#l4c1f44f]
-縦軸が応力、横軸がひずみのグラフを作成した。グラフの傾き...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/0630.png
**20/7/10 [#m5d66eaa]
-functions and listsのlistrのDEFI_LIST_REELのPASを0.1から...
-Timeの10分割を変える方法が分からなかったので、PASを0.01...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/dan.png
**20/7/15 [#t5a1b6e3]
-根元の応力が105MPaになるときの荷重を計算したところ、荷重...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/da.png
**20/7/21 [#wb1b18d7]
-前回は応力とひずみの向きを考慮していなかったので、それを...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/danso.png
-上のグラフは先端のひずみで計算していたので、根元のひずみ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/dans.png
**20/7/24 [#f799d779]
-梁を軸方向に引っ張ったときの解析を行った。引張応力が105M...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/hippari.png
**20/8/7 [#m595aa64]
-今までは面に載荷していたが，集中荷重なので線載荷にしてみ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/sensaika....
**20/10/4 [#h52cbe94]
-Geometryで床版、主プレート、補助プレート、補助金具、支柱...
**20/10/7 [#t0d47e43]
-鋼材のメッシュサイズを決めるために支柱のみの解析をした．
-支柱を長さ600mmの円筒型の片持ち梁と仮定して理論値を計算...
-理論値計算と同じ条件で解析した．ポアソン比は0.3とした．...
,メッシュサイズ,変位
,6,0.000154
,7,0.00015381
,8,0.00015373
,9,0.00015379
,10,0.00015365
,11,0.00015372
,12,0.0001538
-メッシュサイズが1の時は他のサイズと同じ条件でも計算でき...
-この結果から，メッシュサイズが10の時に理論値に一番近くな...
**20/10/8 [#j78223ac]
-CLTを長さ2000mmの単純梁として手計算をした．
-荷重は1N，断面二次モーメントは281250000mm4とした．
-CLTの等級は異等級構成Mx60A(中村先生)
-CLTメモの等級ごとのヤング率の求め方の表3より，Mx60の5層5...
-計算したところ，変位は0.000423mmとなった．
**20/10/9 [#h3e3f621]
-鋼材のヤング率は正しくは206GPaだったので，その数値で理論...
**20/10/12 [#q89d4472]
-鋼材の正しいヤング率で解析した．結果は図のようになった．
,メッシュサイズ,変位
,6,0.000157
,7,0.0001568
,8,0.0001567
,9,0.0001568
,10,0.0001566
,11,0.0001567
,12,0.0001568
-メッシュサイズが10の時に理論値に一番近くなった．相対誤差...
-直交集成板の日本農林規格によると，「直交集成板を構成する...
-Mx 60は外層ラミナの強度等級がM60A、内層ラミナの強度等級...
-M60A,M30Aのヤング率の基準値は直交集成板の日本農林規格に...
**20/10/14 [#n55bd6e2]
-M60A,M30Aのヤング率をそれぞれ6.0GPa,3.0GPaとしてMx60のヤ...
-Shear Analogy Methodという方法で計算した．有効曲げ剛性(E...
-この有効曲げ剛性をCLT全体の断面二次モーメントで割ったと...
-有効せん断剛性も求めた．
-繊維方向のヤング率と各方向のせん断弾性係数の関係式を用い...
-求めた剛性を使ってCLTの変位を求めた．荷重は梁の中心に1N...
**20/10/16 [#f51c805e]
-メッシュサイズを更に増やして支柱の解析をした．結果は図と...
,メッシュサイズ,要素数,変位
,5,145813,0.000157025
,6,144052,0.0001569962
,7,38244,0.0001567935
,8,37009,0.0001567133
,9,36804,0.0001567753
,10,36116,0.000156634
,11,36035,0.0001567057
,12,35293,0.0001567914
,15,9570,0.0001533965
,20,9148,0.0001537212
,25,8747,0.0001540633
,30,2972,0.0001336225
,35,2364,0.0001379303
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ko.png
-このグラフから，要素数を増やすと0.000157当たりに収束して...
**20/10/19 [#n4f0e802]
-CLTの方の解析をした．メッシュサイズを20にしたところ，要...
-メモリの上限を超えているというメッセージが出た．
Error during the dynamic allocation. It was not possible...
! a memory zone length 1169 Mo, one exceeds the maxim...
! built-in with 1493 Mo and one occupies already 580 ...
! The last operation of memory deallocation made it p...
!-----------------------------------------------------...
**20/10/20 [#p7567c5d]
-要素数を少なくしたら，メッシュを作成することが出来なかっ...
-NgException:Stop meshing since too meny attemptsというエ...
-検索してみたがよく分からなかった．
https://www.salome-platform.org/forum/forum_10/156488534
-及川さんに聞いてみたところ，最小のメッシュの大きさが穴の...
**20/10/21 [#haab6b72]
-横2ミリ，縦1ミリの穴を横20ミリ，縦10ミリの大きさにしてみ...
-解析結果は図とグラフのようになった．
,メッシュサイズ,要素数,変位
,14,581948,0.0001126382
,15,443110,0.0001150793
,20,335445,0.0001133839
,25,255690,0.0001076753
,27,131296,0.0001036925
,30,83221,0.000097703671
,40,63683,0.000095167519
,50,54909,0.000093780481
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/syouban.png
-CLTを等方性としてたわみを計算しなおしたら，0.0001247878...
**20/10/28 [#i58662a3]
-支柱にブラケットをつけて計算し直した．
-ブラケットは一辺90ミリの立方体とした．等方性としたのでヤ...
-支柱の高さ588ミリの部分にグループを作成し，そこの変位を...
,メッシュサイズ,要素数,変位
,3,650732,0.0001525962
,4,253843,0.0001528185
,5,226615,0.000152853
,10,57390,0.0001527818
,20,13573,0.0001504562
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/kouz.png
**20/10/29 [#h7d7bf8f]
-プレートと床版が一体化していない状態で解析する必要がある...
-AsterStudyはMATE=AFFE_MATERIAUまでやった．
**20/10/30 [#ad3ee1e9]
-支柱のメッシュをもっと粗くして計算しなおした．結果は表と...
-グラフを見やすくするためにメッシュサイズ3のときの変位は...
,メッシュサイズ,要素数,変位
,3,650732,0.0001525962
,4,253843,0.0001528185
,5,226615,0.000152853
,10,57390,0.0001527818
,20,13573,0.0001504562
,30,5931,0.0001270705
,40,3982,0.0001377993
,50,3411,0.0001174439
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/kouz1.png
**20/11/2 [#n23f4312]
-接触解析の基本(1)の続き
-AsterStudyのCONTACTという項目が見つからない．下のサイト...
http://salome-meca.cocolog-nifty.com/blog/2012/03/60-1-2f...
-接触解析の基本(2)の荷重拘束
-CONTACTの代わりにDEFI-CONTACTを接触の設定として用いた．
-接触の部分はSalomeMecaの使用法解説の接触解析の基本(2)を...
-黄色マークで下のようなメッセージが出た．
<S> Exception user raised but not interceptee. !
! The bases are fermees. !
! Type of the exception: ArretCPUError !
! [('?', (), (), ())]
-このメッセージの意味はよく分からなかった． ...
-Salome接触解析メモを見たところ，部材ごとでメッシュを切ら...
**20/11/3 [#p73e8f5a]
-床版のメッシュを粗くして計算しなおした．結果は表とグラフ...
,メッシュサイズ,要素数,変位
,14,581948,0.0001126382
,15,443110,0.0001150793
,20,335445,0.0001133839
,25,255690,0.0001076753
,27,131296,0.0001036925
,30,83221,0.000097703671
,40,63683,0.000095167519
,50,54909,0.000093780481
,60,21413,0.000076680678
,70,18683,0.000073994244
,80,16084,0.000072255257
,100,14087,0.000070051691
,150,5070,0.000042707938
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/syouban1....
-計算し直す前よりは収束している様子を再現することが出来た．
-一応半解析の方もやってみた．床版を半分に切断し，切断面を...
-切断面にも穴を開けて載荷線を作るという方法でやってみた．...
<S> Exception user raised but not interceptee. ...
! The bases are fermees. ...
! Type of the exception: error ...
! ...
! Problème : la matrice est singulière ou presque sing...
! Lors de la factorisation de la matrice, on a renco...
! (pivot nul ou presque nul) à la ligne 37413 qui co...
! donné ci-dessus. ...
! Le nombre de décimales "perdues" pour cette ligne ...
! Ce nombre de décimales perdues est lié au degré de...
! Plus il est grand, plus le système est singulier. ...
! Quand on a perdu plus de 8 décimales, on estime qu...
! précision ...
! des nombres réels (qui ont 15 décimales environ). ...
! ...
! Risques et conseils : ...
! * Si la ligne correspond a un degré de liberté ph...
! probablement d'un mouvement ...
! de corps rigide mal bloqué. ...
! Vérifiez les conditions aux limites. ...
! Si vous faites du contact, il ne faut pas que l...
! par le contact. ...
! Vérifiez également les caractéristiques matéria...
! ...
! * Si la ligne correspond a un degré de liberté de...
! doute d'une condition ...
! limite redondante. ...
! En particulier, il se peut que la relation liné...
! des conditions de contact. ...
! Peut-être devriez vous exclure certains noeuds ...
! (mots clés SANS_NOEUD et SANS_GROUP_NO). ...
! ...
! * Si le solveur utilisé est LDLT ou MULT_FRONT, v...
! solveur MUMPS ...
! car celui-ci est le seul à pouvoir factoriser l...
! définies positives. ...
! ...
! * Parfois, en parallèle, le critère de détection ...
! trop pessimiste ! Il reste néanmoins souvent ...
! possible de faire passer le calcul complet en r...
! (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou ...
! en le débranchant (valeur du mot-clé NPREC=-1) ...
! sur moins de processeurs. ...
! ...
! * Il se peut aussi que ce phénomène soit tout à f...
! fissure passe ...
! très près d'un noeud. ...
! Si le nombre de décimales perdues n'est pas tro...
! vous pouvez relancer le calcul en augmentant le...
! autorisé : ...
! mot-clé NPREC du mot clé facteur SOLVEUR. ...
! Sinon, contactez l'équipe de développement.
-載荷線のxy方向の固定を解除してみたが，同じメッセージが出...
**20/11/6 [#p89342ae]
-ティモシェンコ梁で理論値の計算をした．
-せん断補正係数は，ポアソン比の0.3を代入して0.5306122449...
-せん断補正係数とせん断弾性係数(79GPa)を代入して計算した...
-支柱の長さを10倍ほどにして解析した．メッシュサイズを10,4...
-10のときに保存しないまま強制終了したのでそのときのhdfフ...
**20/11/8 [#ba52efdb]
-10倍だとできなかったので、とりあえず3倍でやってみること...
,メッシュサイズ,要素数,変位
,10,110765,0.00014746275
,20,27475,0.00014276444
,21,27419,0.000142449056
,23,19402,0.000142579938
,25,17792,0.0001425346
,30,10050,0.00012982425
,40,6912,0.000105412889
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/3sityuu.p...
**20/11/11 [#h10c4fe2]
-接触解析の基本(1)(摩擦なし)の続き
-MODI_MAILLAGE: メッシュの操作
-接触解析メモにMODI_MAILLAGEがなかったのでその項目を消し...
<ArretCPUError> <MECANONLINE9_1> ...
! ...
! Arrêt par manque de temps CPU pendant les itérati...
! - Temps moyen par itération de Newton : 81.619...
! - Temps restant : 130.51...
! La base globale est sauvegardée. Elle contient le...
! ...
! Conseil : ...
! - Augmentez le temps CPU. ...
!-------------------------------------
-このようなエラーメッセージが出た．CPU時間を増やすという...
-timeを15に戻し，STAT_NON_LINE(非線形解析の設定)を弾塑性...
**20/11/12 [#a92df122]
-今までやっていた接触解析のモデルはSTEPファイルを読みこん...
-荷重拘束の接触解析（摩擦あり）のペナルティ法
-Geometry：baceとtopの対称面はそれぞれでグループ化した．
-Mesh：メッシュ→アセンブリをビルドでコンパウンドを作成．G...
-AsterStudy：2.5次要素への変換は、CREA_MAILLAGEコマンドを...
**20/11/13 [#oa5b859b]
-床版(単純梁)の解析
-ふたつ目の固定線が全方向拘束されていたので，kotei2のy方...
-外す前とほとんど変化はなかった．
-理論値に一番近い値のときの相対誤差は-7.78%だった．
,メッシュサイズ,要素数,変位
,14,581948,0.0001126387
,15,443110,0.0001150798
,20,335445,0.000113384
,25,255690,0.000107676
,27,131296,0.0001036928
,30,83221,0.000097704057
,40,63683,0.000095167927
,50,54909,0.000093780519
,60,21413,0.000076680739
,70,18683,0.000073994344
,80,16084,0.000072255329
,100,14087,0.0000700518
,150,5070,0.000042708863
,200,3587,0.00003988833
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/1syouban....
-二次要素への変換方法
-作成された二次要素を右クリック→二次/一次要素への変換→二...
**20/11/16 [#w8a6ceff]
-二次要素への変換
-メッシュサイズ40の支柱モデルを二次要素に変換したら警告メ...
「ノードが （マジェンタ色で図示） 環境によって歪んだ要素...
-メッシュサイズ30のときも同じメッセージが出たが，25のとき...
-2次要素で解析した結果は表とグラフのようになった．(30と40...
,メッシュサイズ,要素数,変位
,10,110765,0.000148931
,20,27475,0.0001490144
,21,27419,0.0001489303
,23,19402,0.000148931
,25,17792,0.0001489472
,30,10050,0.00012982425
,40,6912,0.000105412889
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/2zisityuu...
-一番理論値に近い値の相対誤差は4.39%になった．
**20/11/17 [#q38e7f11]
-床版の方も二次要素で解析した．こちらでは支柱のようなエラ...
-しかし，要素数が大きすぎると計算時にエラーが起こった．
<S> Exception user raised but not interceptee. ...
! The bases are fermees. ...
! Type of the exception: error ...
! ...
! Error during the dynamic allocation. It was not pos...
! a memory zone length 41 Mo, one exceeds the maximum...
! built-in with 1499 Mo and one occupies already 1469...
! The last operation of memory deallocation made it p...
-メッシュサイズ25までは計算できたので，その結果を表とグラ...
,メッシュサイズ,要素数,変位
,200,3587,0.0001210055
,150,5070,0.0001209961
,100,14087,0.000121258
,80,16084,0.0001212511
,70,18683,0.0001212428
,60,21413,0.0001212307
,50,54909,0.0001212352
,40,63683,0.000121233
,30,83221,0.0001212595
,27,131296,0.000121299
,25,255690,0.0001213208
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/2zisyouba...
-グラフでは理論値と離れているように見えるが，一番近い値の...
-Geometryで作った物体を他のファイルに移動させる方法
-移動させたい物体を選択する→ファイル→エクスポート→STEPエ...
-ファイル→インポート→STEPインポート→保存したものを選んで...
**20/11/18 [#s2a5984b]
-軸直交方向のヤング率(床版)を求めた．Shear Analogy Method...
-軸方向のヤング率を4.8704GPaとしていたが，正しくは4.7488G...
-メッシュサイズ25で計算してみたら，より理論値に近くなった...
,メッシュサイズ,要素数,変位
,25,255690,0.0001244172
,27,131296,0.0001243948
,40,63683,0.0001243271
,60,21413,0.0001243248
,100,14087,0.0001243528
,20,3587,0.0001240943
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/2naosi.png
-理論値に一番近いところの相対誤差は-0.3%だった．遠いとこ...
**20/11/19 [#y9db4711]
-支柱が降伏するときの荷重
-I=2343194.144(mm4),σ=238(MPa),l=588(mm),y=57.15(mm)のと...
-16.6kN=1.69tfより，支柱が降伏するときの荷重は1.69tf
-先日の実験データのグラフ化
-荷重と変位の関係をLibreOfficeCalc内で作成した．変位はCLT...
-gnuplotでcsvファイルのグラフを作成する
-区切りをスペースからカンマに変更：set datafile separator...
**20/11/20 [#s2779479]
-実験の荷重と変位，ひずみの関係のグラフを作成
-荷重と変位の関係
-手計算では1.69tfで降伏するという結果になったが，グラフは...
-荷重とひずみの関係(支柱の根元)
-内側は荷重1.8tf当たりから急激にひずみが増加している．外...
-東側上面
-圧縮側にひずみが増加すると考えていたch11,12,13は引張側に...
-東側下面
-ch35〜38は圧縮側に増加していたが，しっかりとした線形では...
-真ん中
-少し線形から外れる点はあったが，特に問題はなかった．
-西側上面
-引張側は特に問題なかったが，圧縮側は0のままのものがあっ...
-西側下面
-ch49が0からほぼ動かない．ch47,48,50が不規則な動き方をし...
**20/11/25 [#cf74f64c]
-とりあえずLibreOfficeでグラフを作成した．荷重と変位，支...
-撮った写真はアプリでJPEG形式に変換した．
**20/11/26 [#l69ade34]
-gnuplotで荷重とひずみ(支柱根元)の関係のグラフを作成した...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/nem...
-CLTの縦ヤング率の手計算をしたら，150MPaになった．それを...
-確かめるため，解析も行った．繊維方向のヤング率，繊維直交...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/tat...
**20/11/27 [#la3d8e74]
-10倍の支柱の解析を行った．メッシュを作るとき，いつものや...
-しかし，メッシュサイズを変えても要素数が全く変わらない．...
**20/11/30 [#zcbef286]
-10倍の支柱の解析をした．Maxサイズを小さくすると要素数が...
,メッシュサイズ,要素数,変位
,Max100 min20,42192,1.3718466667
,M100 m10,67880,1.4224408333
,M100 m4.5,139661,1.4488306757
,M10 m4.5,357844,1.4893606757
,M8 m4.5,493323,1.4849333784
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/10s...
-全体モデル
-全部くっついている状態，CLTの切断面を固定，ブラケットに...
-ポアソン比をどう決めればいいかまだわからない．
**20/12/1 [#i663c898]
-写真をk2の/home/gakusei/zikken/cltにアップロードしました...
**20/12/2 [#ube1770d]
-ポアソン比の計算をした．構造関係メモの直交異方性材料のポ...
-1/25じゃなくて1/30にしたのは，CLTメモ内のShear Analogy M...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/poi...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/poi...
**20/12/3 [#h42dd35c]
-せん断弾性係数の計算
-全体モデルでの解析
-外桁の固定面のグループをメッシュに反映させようとしたらエ...
**20/12/4 [#s6814cf3]
-せん断弾性係数は、どの方向も繊維方向のヤング率の1/15
-エラーが起こる前からやり直してみたが、何も問題なく反映さ...
**20/12/7 [#y1dc24e8]
-弾性の全体モデルの解析を行った．荷重は降伏する直前の2.5t...
,全体モデル(弾性)
,荷重(t),変位(mm)
,0.5,1.329206
,1.0,2.658416
,1.5,3.98762
,2.0,5.31683
,2.5,6.646036
,実験値
,荷重(t),変位(mm)
,0.5,2.625
,1.0,5.5
,1.5,8.4
,2.0,11.55
,2.5,14.7
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/mod...
-以下のように修正した．
-切断面の縦方向を動けるようにした．
-外桁部分にゴムをつけた．
-ゴムとH鋼の接触面を全面固定した．
-ゴムのヤング率：0.002GPa,ポアソン比：0.49とした．
-とりあえず荷重2.5ｔで解析したところ，変位は9.66262375mm...
-次に，2次要素に直してみると，変位は12.7262142857となった．
-2次要素ならメッシュが粗くても精度は変わらないだろうと考...
-それでも相対誤差は-13.4%ほどある．
**20/12/8 [#qaf23da4]
-2次要素のメッシュサイズMax100,min3.2で0.5tから0.5刻みで...
-実験値は11/24のもの
,荷重(t),変位(mm)
,0.5,2.5460526316
,1,5.0921063158
,1.5,7.6381589474
,2,10.1842136842
,2.5,12.7302684211
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/rub...
**20/12/9 [#r8407b72]
-ゴムのヤング率を2MPaから100MPaに直した．ゴムが固くなった...
-実験でジャッキのヘッド側が回転し，斜めに荷重をかけていた...
-床版の真ん中の変位が実験と解析値でかなり差があったので理...
**20/12/16 [#qd837b19]
-接触解析ここを参考にやってみる．
https://qiita.com/Jun_Tatsuno/items/6159d914abdc75a91262
-変位制御で解析することが出来たが，WarpByVectorを見ること...
→GenerateVectorsを押したら出来た．
-荷重制御に変えてみたら，ニュートン法の最大繰り返し回数で...
**21/2/9 [#t3ae3af3]
-https://www.kochi-kensoku.org/q-a%E3%82%B3%E3%83%BC%E3%8...
-ガードレールの衝突荷重について
**21/2/11 [#cc6cfb40]
-http://www.pwrc.or.jp/thesis_shouroku/thesis_pdf/0406-P0...
-ガードレールの衝突荷重について２
**21/2/12 [#kb8827b2]
-http://www.jsce.or.jp/committee/steel/subcomit-active/hp...
-降伏点，安全率について
*卒論日誌 [#ff45fcb7]
,日付,時間帯,作業時間,内容,立会
,10/11,16:00-17:00,1h,ガイダンス,後藤さん
,10/18,14:30-16:30,2h,創造工房実習,後藤さん
,10/21,13:00-16:00,3h,課題,
,10/21,14:30-16:30,2h,創造工房実習,後藤さん
,10/28,13:30-16:00,2.5h,課題,
,11/1,14:30-16:00,1.5h,創造工房実習,後藤さん
,11/6,11:30-12:30,1h,課題,
,11/8,14:30-16:00,1.5h,創造工房実習,後藤さん
,11/11,14:00-15:00,1h,課題,
,11/13,12:00-13:30,1h,課題,
,11/15,14:30-16:00,1.5h,創造工房実習,後藤さん
,11/20,12:00-12:30,0.5h,課題,
,11/22,14:30-16:00,1.5h,創造工房実習,後藤さん
,11/27,12:00-14:00,2h,課題,
,12/6,14:30-16:00,1.5h,創造工房実習,後藤さん
,12/11,11:00-13:30,2.5h,課題,
,12/13,14:30-16:00,1.5h,創造工房実習,後藤さん
,12/18,12:00-14:00,2h,課題,
,12/19,14:30-16:30,2h,課題,
,12/20,14:30-16:00,1.5h,創造工房実習,
,1/10,14:30-16:00,1.5h,創造工房実習,後藤さん
,1/15,14:30-15:30,1h,創造工房実習,後藤さん
,1/24,14:30-16:00,1.5h,創造工房実習,後藤さん
,1/29,18:00-19:30,1.5h,課題,
,1/29,21:00-22:00,1h,課題,
,1/31,14:30-16:00,1.5h,創造工房実習,後藤さん
,1/31,16:00-17:30,1.5h,課題,
,2/1,17:00-17:30,0.5h,課題,
,2/3,14:00-14:30,0.5h,課題,
,2/3,16:30-19:00,2.5h,課題,
,2/3,20:00-22:00,2h,課題,
,2/4,10:30-11:30,1h,課題,
,2/4,13:00-16:00,3h,課題,
,2/5,12:00-17:00,5h,課題,
,2/6,14:30-15:00,0.5h,課題,
,2/7,14:30-15:00,0.5h,創造工房実習,後藤さん
,2/13,13:45-14:45,1h,課題,
,2/14,14:30-15:30,1h,創造工房実習,後藤さん
,2/18,9:30-10:00,0.5h,春休み課題,
,2/26,14:00-16:30,2.5h,春休み課題,
,2/27,15:00-16:30,1.5h,春休み課題,
,3/2,16:45-17:45,1h,春休み課題,
,3/5,16:30-17:30,1h,春休み課題,
,3/7,16:00-18:30,2.5h,春休み課題,
,3/9,16:00-17:30,1.5h,春休み課題,
,3/11,15:30-17:30,2h,春休み課題,
,3/16,14:00-18:30,4.5h,春休み課題,
,3/17,13:00-16:00,3h,春休み課題,
,3/18,13:30-17:30,4h,春休み課題,
,3/19,14:00-17:30,3.5h,春休み課題,
,3/23,16:00-18:30,2.5h,春休み課題,
,3/24,15:00-17:30,2.5h,春休み課題,
,3/25,15:00-16:30,1.5h,春休み課題,
,3/26,16:00-17:30,1.5h,春休み課題,
*10/18の課題 [#xef89c31]
,メッシュ長さ,要素数,先端変位,相対誤差,計算者
,0.5,394121,6.57938,-1.4%,佐藤
,1,72278,6.443,-3.4%,吉田
,2,11817,5.7382525,-13.9%,梅田
,4,2862,4.9428,-25.9%,青山
,8,897,4.0411725,-39.4%,小川
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/ogawa/kat...
*10/25の課題 [#k79ee126]
,メッシュ長さ,要素数,先端変位,相対誤差,計算者
,0.5,394121,6.57938,-1.4%,佐藤
,0.7,130916,6.4781,-2.9%,佐藤
,0.8,72101,6.437,-3.5%,青山
,0.9,71718,6.4314,-3.5%,吉田
,1,72278,6.443,-3.4%,吉田
,1.2,65575,6.443,-3.4%,小川
,1.4,41096,6.3162,-5.3%,青山
,1.5,23417,6.1206,-8.2%,梅田
,1.8,11758,5.7370,-13.98%,梅田
,2,11817,5.7382525,-13.9%,梅田
,4,2862,4.9428,-25.9%,青山
,8,897,4.0411725,-39.4%,小川
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/ogawa/kat...
考察
メッシュ長さ1の時の要素数が0.9、0.8の時よりも多くなってい...
*11/1の課題 [#w19e33c2]
,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.5,322687,0.428774714,2.9%,小川
,0.8,199968,0.41646,-0.050%,吉田
,1,62360,0.416513545,-0.037%,梅田
,1.5,50359,0.411953,-1.134%,青山
,2,20526,0.3706311667,-11.05%,佐藤
理論値 0.416667
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/ogawa/tan...
考察 メッシュを細かくするほど精度が高くなるわけではない...
*11/8の課題 [#n612ec71]
,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.5,322687,0.5214958571,6.07%,小川
,0.8,199968,0.5145437143,4.65%,吉田
,0.9,87953,0.5102753846,3.78%,梅田
,1,62360,0.505768727273,2.87%,青山
,1.1,83278,0.5033813636,2.38%,小川
,1.2,59515,0.5053267,2.78%,吉田
,1.3,57455,0.5038435556,2.47%,青山
,1.5,50359,0.50010225,1.71%,梅田
,1.8,28053,0.4926668571,0.20%,佐藤
,2,20526,0.4765718333,-3.08%,佐藤
ティモシェンコ梁 0.4916667
理論値 0.416667
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/ogawa/tan...
*11/15の課題 [#ee35932d]
わかったこと
・自動メッシュの3次元を選択までの流れ
わからなかったこと
・スケール係数の欄がどこにあるか
・固定面の節点の選択のしかた
・エンティティの追加の欄がどこにあるか
*11/22の課題 [#uf06888c]
,拡大係数,要素数,変位,相対誤差,計算者
,0.3,57431,6.422525,-3.7%,全員
,0.4,23349,6.26413,-6.0%,佐藤
,0.5,13741,6.03154,-9.5%,小川
,0.6,7257,5.91187,-11.3%,梅田
,0.7,4245,5.5716725,-16.4%,吉田
,0.8,3657,5.5734,-16.4%,青山
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/ogawa/kat...
*12/6の課題 [#e860c736]
,拡大係数,要素数,変位,相対誤差,計算者
,0.8,122515,6.91749,2.6%,吉田
,1,65490,6.8877,2.2%,梅田
,1.2,36693,6.8519225,1.6%,梅田
,1.4,23341,6.8117175,1.0%,吉田
,1.6,15371,6.6336,-1.6%,小川
,1.8,12480,6.74258,0.014%,佐藤
,2,8589,6.66499,-0.011%,梅田
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/ogawa/kat...
節点数で荷重を割るときに精度を上げていれば変位が下がらず...
有効数字
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
