![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
#author("2020-01-20T11:33:26+09:00","default:kouzouken","kouzouken") #contents *卒論タイトル案 [#k9de87f9] -CLT床版の数値モデル化 *作業メモ [#ue872716] **2018/1/29 [#j7a575f2] -異方性になると塑性を入れたらエラーが出て解けない。 -異方性は一先ず置いておき、片持ち梁と2層モデルについて解析を行った。片持ち梁は解析できたが、2層モデルについては、また異方性ではエラーが出て解析ができなかった。 **2018/1/27 [#x07d7194] 1/29までの目標 -これまで使用していたモデルからゴム板を外し、木材のみのモデルで計算することにした。このモデルについて、まずは弾性(等方性と異方性の両方)で解き、その後それに塑性を入れて解いていく。 →弾性では当然解け、塑性についても等方性は解けたので、次は異方性でやってみる。それも出来たら、今度はモデルを2層にして計算できるかやってみる。 -片持ち梁についても、梁を引っ張った時の弾塑性解析を行う。 **2018/1/26(後藤メモ) [#wf1c387a] -&link(ここの例題にある,http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wiki.cgi?Salome%2dMeca2017) 片持ち梁を曲げるだけの &link(katamoti0126p.hdf,http://www.str.ce.akita-u.ac.jp/~gotouhan/salome/hdf/katamoti0126p.hdf) hdfファイルに、AsterStudyのFunctions and ListsのDefine functionから、 cltdanso11.hdfとかと同様の以下の設定をしても、エラーが出ずに解ける。 --( elast_pl = DEFI_FONCTION( NOM_PARA='EPSI', NOM_RESU='SIGM', VALE=(0.5, 5000000.0, 1.0, 5000000.0) --) ということは、応力-ひずみ曲線の設定がまずいわけではなくて、 他のどこかの設定がまずいしわよせのせいで、 応力-ひずみ曲線が見つかりませんみたいなエラーを出している可能性もあります。 もう一度、解けるCLTの弾性問題のhdfに、Salome-Meca2017で、 応力-ひずみ曲線の設定をしてみましょう。 解ける片持ちばりとは何が違うのか。 あと、単位系は、N-mm系ではないですか。 降伏点が50MPaだとすると、50N/mm$^2$ではないですか。 **2018/1/25 [#y9bf52c8] -Salome-Meca2016で、去年の弾塑性の.commファイルを読み込んで、藤田さんと同じ解析ができるかどうかチェック。 →Salome-Meca2017と同じエラーが出た。 -もし、できるのであれば、その.commファイルの降伏点をスギ材の一般的な引張強度、圧縮強度に書き換えて計算。 --上記ができるのであれば、それを8層に修正できるかやってみる -めり込みを考慮して修正した実験値と上記を比較。 -めり込みを考慮しない場合も一緒にプロット。 **Salome-Meca2017 [#sc807ff3] 去年の弾塑性の.commファイルをSalome-Meca2017に読み込むと、降伏点の$(\varepsilon, \sigma)$の座標がData SettingsのFunctions and Listsに書き込まれている。 が、この設定で走らせると、 --( On ne trouve pas la courbe de traction (mot-clef TRACTION) dans le matériau ! ! fourni. --) traction(張力のことだが、応力のこと?)の曲線（たぶん、弾塑性の応力ひずみ曲線）が材料の設定に見つからないみたいなエラー。 弾塑性曲線の与え方（横軸、縦軸の設定、与える点の個数など）をいろいろ変えると、微妙にエラーが変わるが、解けない。 Salome-Meca2016でも解けなかったのでもう一度Salome-Meca2017で行ってみると、下のエラーに変わったが、解けなかった。 --( comportement TRACTION non trouve --) traction(張力か応力のこと?)の挙動が見つからないみたいな意味。 *卒論日誌 [#q4af9a7c] ,日付,時間帯,作業時間(h),内容,立会 ,4/7,10:00-11:00,1,ウィキの書き方,後藤さん ,4/17,16:00-18:00,2,vi, ,4/24,17:00-18:00,1,vi、fortran, ,4/27,17:00-18:00,1,vi、fortran, ,5/8,16:00-18:00,2,salome, ,5/25,16:00-23:00,7,salome, ,6/1,16:00-0:00,8,salome, ,6/2,0:00-3:00,3,salome, ,6/4,0:00-2:00,2,salome, ,6/12,16:00-21:00,5,salome, ,6/26,16:30-19:30,3,tex, ,7/3,17:00-19:00,2,tex, ,7/10,14:30-18:30,4,salome, ,7/18,15:00-17:30,2.5,salome, ,7/10,15:00-17:30,2.5,salome, ,7/27,16:30-18:30,2,salome, ,7/28,15:30-18:00,2.5,salome, ,7/29,17:30-22:30,5,salome, ,7/30,18:30-23:00,4.5,salome, ,9/5,14:30-16:30,2,salome, ,9/9,11:30-19:30,8,salome, ,9/18,12:30-22:30,10,salome, ,9/19,0:00-3:00,3,salome, ,9/20,14:00-21:00,7,salome, ,9/22,10:00-13:00,3,salome, ,9/22,15:30-18:30,3,salome, ,9/22,21:00-0:00,3,salome, ,9/23,0:00-2:30,2.5,salome, ,9/24,19:00-0:00,5,salome, ,9/25,0:00-3:00,3,salome, ,9/25,4:30-8:30,4,salome, ,9/29,16:00-20:00,4,salome, ,10/6,15:00-16:00,1,salome, ,10/20,9:00-13:00,4,salome, ,10/26,11:30-14:00,2.5,salome, ,10/27,12:00-13:00,1,salome, ,10/30,14:00-19:00,5,salome, ,11/3,23:30-2:30,3,salome, ,11/18,0:30-9:30,9,salome, ,11/25,0:30-5:30,5,salome, ,11/26,0:00-6:00,6,salome, ,12/21,16:30-22:30,5,salome, ,1/1,9:00-12:00,3,salome, ,1/16,15:30-3:30,12,salome, ,1/17,17:00-4:30,11.5,salome, ,1/19,13:00-0:00,11,salome, ,1/25,13:30-4:00,14.5,salome, ,1/26,5:00-7:30,2.5,salome, ,1/26,14:00-19:00,5,salome, ,1/27,1:00-6:00,5,salome, ,1/29,13:00-4:00,15,salome, ,2/1,10:00-2:00,16,salome, ,2/2,7:00-1:00,18,salome, ,2/4,8:00-2:00,18,salome, ,2/5,5:00-8:00,3,salome, ,2/5,8:30-5:30,21,salome, ,2/7,9:30-3:30,18,salome、tex, ,2/9,8:00-13:00,5,salome、tex, ,2/9,14:00-19:30,5.5,salome, ,2/10,7:00-3:00,18,salome, ,2/11,10:00-2:00,16,salome, ,2/15,9:00-2:30,15.5,salome、スライド作成, ,2/16,3:30-23:30,20,salome、スライド作成, ,合計,412.5, **卒論メモ(1/19) [#m2edc5ca] -salome2017で弾塑性だと等方性と異方性の両方でエラーが出てしまった。 --なので一旦弾塑性をやめて等方性と異方性で解き、それに手を加えて弾塑性で解けるかという方針に変更した。 →弾塑性では無い方は両方で解けたが、未だ弾塑性では解けていない。 **卒論メモ(1/16) [#u26e85c1] -古いsalomeで上手くいかなかった異方性で解く方法をsalome2017で解けるかやってみる。 --エラーが出て解けなかった。 -Analysis->STAT_NON_LINE->COMPORTEMENT->"+" -> Edit... --RELATIONとDEFORMATIONの部分に違いがあった →ここの違いがエラーの原因?かもしれない(確証はない)。 **卒論メモ(11/10) [#q6f7bba3] -去年の卒論の弾塑性をトレース --まずは等方性で解く。(11/13の12時までに) **卒論メモ(10/27) [#ff38f12f] -変位を細かくとる --グラフ化する →変位はとったがグラフ化はまだ。 -去年の卒論の弾塑性をトレース --まずは等方性で解く。新旧のsalomeで解く。 →まだ解けていない。commファイルの編集が上手く行かない。 --次に異方性で解く。 **卒論メモ(10/6) [#k08470b3] -作成したモデルの変位を細かくとる。 -ゴムを敷いて剛性を変える。 **後藤メモ(17/8/22)TeX関係 [#i6d938cf] -&link(sibuw.tar.gz,http://www.str.ce.akita-u.ac.jp/~gotou/tebiki/pdf/sibuw.tar.gz)を適当なディレクトリに保存。 -tar xvzf sibuw.tar.gzで展開。sibup2.texを書き換えながら、中間発表概要を作成。 -グラフはなるべくgnuplotで作成。 -gnuplotの図をTeXに貼り付けるやり方は、色々とあるが、 -&link(gnuplotの図に矢印や説明を加えてLaTeXに取り込む(PS/LaTeX編) ,http://www.str.ce.akita-u.ac.jp/~gotou/linux/vine.html#pstex)を標準的方法とする。 **後藤メモ(17/7/31)卒論ネタ [#pa190599] 後藤補足(17/8/3)以下の懸案を一つずつ片付けていきたいですが、まずは、藤田さんの論文と同じ弾性・弾塑性の計算結果が出せるようになったら、拘束面に実験と同じようなゴム板（程度のヤング率の物体）を挟んで解析結果がどう変わるか、あたりが手をつけやすいでしょう。 -（まずはこれを夏休み中にやる）実験では、支点部近くに載荷した場合、FEMのように支点側のたわみが小さくならない問題 --H形鋼、ゴム板を、実験と同じようにモデル化したらどうか（ゴムの変形の影響をチェック） --更に、H形鋼とゴム版を試験体より長めにとって解析したら上記と違いが出るか --下面中央部の変位を出力して、gnuplotで描ける --コムシートを敷いて同じ計算をする -能代の破壊試験の初期剛性が、山口大の静的載荷試験と合わない原因を調査 --山口大の初期剛性はFEMと合う。 --仮に破壊試験のたわみ計などの係数設定が間違っていたとして、 --山口大の剛性と合うように係数をかけ直したら、 --弾塑性解析の際に公称値の強度を降伏点に与えても、実験に近い破壊強度を示すのでは。 -6層モデルだと、引張側の降伏が起きやすくなるか（実物だと、8層では引張側の破断はしないが、6層だと破断する？） *5/15課題 [#r0ccc95b] 異方性材料にして計算する。 <1次要素> ,Length,1,2,4,8 ,たわみ(mm),409.3,387.4,355.9,252.3 ,相対誤差 (%) (理論値:初等梁=416.7mm),1.776,7.031,14.59,39.45 ,相対誤差 (%) (理論値:ティモシェンコ梁=416.8mm),1.325,6.731,13.59,37.85 <2次要素>Lengthが細かすぎて計算できなかった ,Length,1,2,4,8 ,たわみ(mm), ,相対誤差 (%) (理論値:初等梁=416.7mm), ,相対誤差 (%) (理論値:ティモシェンコ梁=mm), *5/8課題 [#m0fb400e] 10mm×10mm断面の角材（ヤング率：6GPa, ポアソン比：0.3）の単純梁のたわみと、初等梁($\frac{P\ell^{3}}{48EI}$)、ティモシェンコ梁($\frac{P\ell^{3}}{48EI}+\frac{P\ell}{4kGA}$)に対する相対誤差を求める。 ・荷重は面荷重で100N. ・スパン長さは1000mm. ・要素分割は、Lengthが8, 4, 2, 1 ・1次要素(linear)と2次要素(quadratic) <1次要素> ,Length,1,2,4,8 ,たわみ(mm),402.8,368.7,318.6,210.0 ,相対誤差 (%) (理論値:初等梁=416.7mm),3.336,11.52,23.54,49.60 ,相対誤差 (%) (理論値:ティモシェンコ梁=416.8mm),3.359,11.54,23.56,49.62 <2次要素> ,Length,1,2,4,8 ,たわみ(mm),Lengthが細かすぎて計算できなかった,416.8,416.8,416.8 ,相対誤差 (%) (理論値:初等梁=416.7mm),Lengthが細かすぎて計算できなかった,0.024,0.024,0.024 ,相対誤差 (%) (理論値:ティモシェンコ梁=416.8mm),Lengthが細かすぎて計算できなかった,0,0,0
