![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
#author("2020-01-20T11:57:06+09:00","default:kouzouken","kouzouken") #contents *タッチタイピング記録 [#ue4321cc] ,日付,記録 ,4/20,2分19秒89 *ccx座屈荷重とオイラー座屈荷重 [#s1fb8e3f] 正方形断面、単純支持としてオイラー座屈荷重を求める。 詳しい寸法はf90ファイルの中にあります。 **正方形断面 [#m34a54fb] ***3点曲げ [#g22d1cd2] 細長比138.5640 Pcr=4.2357E-002 ccxでは単純支持は中立軸に線拘束、 載荷点は一面載荷にする。 --単純支持を線拘束でなく断面一面拘束してしまうとおかしい値になる -支点からはみ出すモデルがよいかどうか比較 ,,支点から少しはみ出す,支点からはみ出さない ,n次モード,座屈荷重（MN), ,1,0.4246390E-01, 0.4245872E-01 ,2,0.8387729E-01, 0.7136013E-01 ,3,0.1688250E+00, 0.1687385E+00 ,4,0.2461169E+00, 0.2113951E+00 ,5,0.3760653E+00, 0.3755901E+00 2次モードで座屈荷重が違う。。 --単純支持で端部を支点から少し飛び出すモデル &link(kouoz1.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/kouoz1.f90) --端部が支点から飛び出さないモデル &link(ifkou1.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/ifkou1.f90) **3点曲げの半解析 [#n503d033] 片方中立軸にローラーを線でいれる、 スパンを半分にした断面すべてでz方向拘束 -x方向拘束なし ,n次モード,座屈荷重（MN), ,1,0.4242539E-01 ,2,0.4246346E-01 ,3,0.3755969E+00 ,4,0.3760319E+00 ,5,0.1010655E+01 総荷重1MNを分配してccxで出てきたBUCKLINGの値をそのまま読めて、 そのまま読んだ値が3点曲げ解析の座屈荷重やオイラー座屈荷重と近い。 1次モードはかなり近い値になるが、2次モードからは違う値になる。 x方向の拘束なしでやってるけど入れられるところはあるのでしょうか。 &link(kouzahan2.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/kouzahan2.f90) **長方形断面 [#h8f2f7a7] --( b=3.2d-2!幅 h=2.0d-2!高さ --) くらいの長方形断面にしてみる。 細長比138.5640 オイラーの座屈荷重$Pcr=6.777E-002$ ,n次モード,3点曲げ座屈荷重(MN),半解析座屈荷重(MN) , 1, 0.6797019E-01 , 0.6796948E-01 , 2, 0.2702622E+00 , 0.1726180E+00 , 3, 0.3397393E+00 , 0.6020786E+00 , 4, 0.6021329E+00 , 0.1490727E+01 , 5, 0.9855231E+00 , 0.1621735E+01 やはり1次モードでは同じような値が得られるが、2次モードから違う値が出る。 座屈挙動を確認してみると3点曲げモデルでは1次,2次共にy方向に座屈がおこる。 半解析モデルでは、1次はy方向に座屈が起こっているが、 2次モードではx方向に座屈してしまっている。 ccxで座屈解析をして式の値と比較するときは、cgxで想定した方向に座屈しているか、全部のモードでちゃんと確認したほうがいいと思う。 *gnuplotで線形回帰 [#w0cf364d] 例えばkaiki.txtに、ひずみ-応力 のように並べて、 --( set term png set output 'hoge.png' f(x)=a*x+b fit f(x) 'tate' using 1:2 via a,b plot 'tate' pt 2,f(x) w l --) と書いて、コマンドラインで --( gnuplot<kaiki.txt --) yang.f90,syokignu2.f *鉄の引張 片持ち梁(弾性、弾塑性解析) [#r845dc6d] -鉄の種類:SS400 -降伏応力:235N/mm^2 －鉄の寸法 --長さ:1m --幅:10cm --高さ:10cm --ヤング率E=206GPa --ポアソン比ν=0.3 -断面積:100cm^2=10000mm^2 -分割数：nx=128,ny=10,nz=10 **calc [#wce20b6c] 節点荷重での解析(断面一様分布) ***弾性解析 [#g592ba05] ,応力(MPa),軸方向変位最大値(m),断面平均軸方向変位(m) ,200,1.019780e-03,9.71303E-04 ,235,1.198247E-03, ***弾塑性解析 [#ic1c66c2] --( print'("*DEFORMATION PLASTICITY")' print'("206.E3,0.3,235.,1.,0.2")' --) Exponent(n)を変えてみる。(式中の指数) DEFORMATION PLASTICITYについては&link(柴田さんの修論日誌,http://www.str.ce.akita-u.ac.jp/kako/j2011/shibata.html#i86)を参照。 datの中で各荷重段階にあわせてtimeごとに変位が出力される。 -n=1 ,応力(MPa),軸方向変位最大値(m) ,10,6.109180E-05 ,50,3.056134E-04 ,100,6.115850E-04 ,150,9.179159E-04 ,190,1.163240E-03 ,200,1.224603E-03 ,210,1.285988E-03 ,220,1.347380E-03 ,230,1.408793E-03 ,235,1.439502E-03 **mentat弾塑性解析 [#a5c45605] 面荷重での解析 -解析モデル ,応力(MPa),軸方向変位最大値(m),ひずみ,ヤング率(GPa) ,50,2.420E-4,2.420E-4,206.6 ,100,4.839E-4,4.839E-4,206.7 ,150,7.259E-4,7.259E-4,206.6 ,200,9.681E-4,9.681E-4,206.6 ,230,1.114E-3,1.114E-3,206.5 ,235,3.234E-3,3.234E-3, -応力とひずみの関係(縦軸:応力 横軸:ひずみ) http://www.str.ce.akita-u.ac.jp/~gotouhan/oyama/ss400.png *hariw.fの使い方 [#h2e4f18c] --( #倍率0のままで要素数、断面諸量などを適時変えて oir.fコンパイル実行 ./oir#1つのoiraa.dができる ./hariw<oiraa.d #1つのfort7とfort11とuvw.outができる ./modo2#fort7を読み込んでモード計算。modoが7個できる(ゴミの外乱) #oir2.fの倍率、弧長、弧長増分を適時変えてコンパイルして実行 ./oir2#oir001～oir010まで倍率を変えてファイルができる chmod 744 job ./job#oir001～010まで自動で読み込んで設定反復回数を終えると自動で止まる --) *3点曲げの半解析 [#idee4c37] -b=1cm,h=1cm,ell=50cmの鋼材の三点曲げ。荷重P=100N。 -荷重は支間中央部に線載荷する。半解析モデルでも線載荷とする。 -ヤング率206GPa,ポアソン比0.3 -C3D8要素、nx=ny=10,nz=600 **片方ヒンジ、片方ローラー [#e9fbe11e] ,初等梁理論のたわみ,ティモのたわみ ,1.5170E-003(m),1.5189E-003(m) ***3点曲げモデルたわみ [#wdd812a0] -拘束条件を簡略化 ヒンジ部分幅方向(x方向)に全てxyz=0 ローラー部分x方向に全てxy=0 支間中央線載荷 ,,,相対誤差(初等梁),相対誤差(ティモ) ,最大値(m),1.512610e-03,0.289%,0.413% ,中央載荷部断面平均(m),1.51235E-03,0.306%,0.430% --mentatで解く 最大値1.519(mm) -拘束条件を理想的に ヒンジ部分x方向の真ん中(b/2)だけxyz=0、他はyz=0 ローラー部分x方向の真ん中だけ(b/2)xy=0,他はy=0 ,最大値(m),中央載荷部断面平均(m) ,1.512610e-03,1.51235E-03 ***片持ちとして解く [#nbbe7c51] ell/2=25cm,荷重P/2=50N 自由端に線載荷 ,,,相対誤差(初等梁),相対誤差(ティモ) ,最大値(m),1.4955e-03,1.416%,1.540% ,自由端断面平均(m),1.49532E-03,1.429%,1.551% ***半解析 [#w40c5683] ell/2=25cm,荷重P/2=50N 片方ローラー、もう片方の断面をz方向拘束、断面の真ん中(x=b/2)だけx方向拘束 ,,,相対誤差(初等梁),相対誤差(ティモ) ,最大値,1.515650e-03,0.0883%,0.212% ,断面平均,1.51539E-03,0.106%,0.230% **両方ヒンジ [#w525d178] ***3点曲げモデルたわみ [#hd5e78d8] -拘束条件を簡略化 ヒンジ部分幅方向(x方向)に全てxyz=0 支間中央線載荷 ,最大値(m),中央載荷部断面平均(m) ,6.910970e-04,6.90936E-04 --mentatで解く 最大値6.978e-1 (mm) -拘束条件を理想的に ヒンジ部分x方向の真ん中(b/2)だけxyz=0、他はyz=0 ,最大値(m),中央載荷部断面平均(m) ,6.914600e-04,6.91299E-04 ***半解析 [#hdf97016] ell/2=25cm,P/2=50N 片方ヒンジ、もう片方は断面一面z=0で、b/2の部分だけxz=0 ,最大値(m),載荷部断面平均(m) ,6.938740e-04,6.93714E-04 *各種モデル最新版 [#kaf5a100] ***プレストレス木箱桁橋片持ち箱断面モデル [#k16aed87] &link(onsi14.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/model/onsi14.f90) ***プレストレス木箱桁橋片持ち補剛材ありモデル [#u77c5ef0] &link(onsih14.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/model/onsih14.f90) ***プレストレス木箱桁橋片持ち箱断面上部鋼板欠けモデル [#l46a3737] &link(onsikake.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/model/onsikake.f90) *調査用プログラム [#o5ad5246] **c3d6要素で直方体を作るプログラム [#h1afb6c6] ***一気にinpを作る [#dadc6f78] &link(ifc3d6.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/model/ccxtyousa/ifc3d6.f90) -1つのinp作る &link(c3d6t.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/model/ccxtyousa/c3d6t.f90) **c3d8要素で直方体を作るプログラム [#j98e10db] ***c3d8要素で一気にinpを作る [#l5b63f6e] &link(ifc3d8.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/model/ccxtyousa/ifc3d8.f90) *unv→inp変換プログラム [#td70de16] 大体プレストレス木箱桁橋用に特化してる **3点曲げ [#n3ceb722] &link(3tenc3d4unv.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/prog/3tenc3d4unv.f90) **半解析 [#mf20d7bc] &link(hanc3d4unv.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/prog/hanc3d4unv.f90) **d4d6d8unv.f90 [#j3b3fd1a] -c3d4,c3d6,c3d8要素に対応。 要素が混在してると、unv内でどこでどの要素が出てくるか、というのはモデルによって違ってくるので、要素の種類ごとに適当な名前の一時ファイルを作ってそこに要素の情報をいれといて、あとでそこから拾ってくるという方式にしてます。 &link(d4d6d8unv.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/d4d6d8unv.f90) -押し出しによるメッシュの作成の例を試したかったため作成 http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/salomeosidasi.png -参考 http://www.geocities.jp/penguinitis2002/study/Salome/extrusion/extrusion.html *計算用プログラム [#y007a4e9] -等方性材料片持ち梁の理論値計算確認用(&link(keitouhou.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/keitouhou.f90)) -異方性材料(木材)片持ち梁の理論値計算確認用(&link(keimoku.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/keimoku.f90)) *課題 [#m1bcaea2] **PCゼミ14/6/16 [#h1a9db26] **パソコンゼミ14/5/26 [#c8e7dd16] -mokuzai.f90 をコンパイルして、実行して出てくるテキストあるいはファイルに書きこませたテキストの、 --( ～～ *ELASTIC,TYPE=ORTHO ←ここから [#fc9b66c4] ～～ ～～ ～～ ←ここまで (要は、*SOLID SECTION の一行上まで) *SOLID SECTION,ELSET=Eall,MATERIAL=mokuzai [#pdfcc277] --) をコピーし、inpファイルの中身の --( *ELASTIC [#t183cb96] ヤング率、ポアソン比 --) という部分を消して、その部分に貼り付ける。 -片持ち梁の理論値計算確認用(&link(kei14.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/kei14.f90)) ***正方形断面の木材 片持ち梁 [#yd82ba1b] 12cm × 12cm 断面の木材、荷重100N,k=5/6 とする。 ヤング率7.000GPA,直行異方性の木材として片持ち梁で解く。 (&link(ifyousoz.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/ifyousoz.f90)) で長さを30cm～500cmまで10cmずつ変えてたわみの平均を取る。 C3D8要素とする。 要素分割はx=10,y=10,z=600とする。 -相対誤差と長さのグラフ http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/redsyoto-gretimo.png --y軸:相対誤差(%),x軸:長さ(cm) --赤線:初等梁理論のたわみとFEMの相対誤差,緑線:ティモシェンコ梁のたわみとFEMの相対誤差 -せん断項とell/hのグラフ http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/tawami-lh.png --y軸:せん断項/ティモシェンコ梁たわみ ,x軸:ell/h **14/5/12 [#zf1fd58d] -Salomeでモデル生成、C3D8でメッシュ作成、unvをinpに変換して解いてみる ヤング率7.694GPa ,ポアソン比0.4 幅80mm、高さ40mm、長さ800mm、荷重1MN -lengthは1つの要素の最大の長さ ,length,たわみmax(m) ,10,4.444930e+01 ,7,4.444540e+01 ,初等梁,せん断補正項,ティモ ,相対誤差(%) , 51.98856,0.10917598,52.09773,14.680 -このモデルではlengthを2以下にすると要素数が多くなりすぎて計算できない。 ***正方形断面の鋼材 [#e91c3673] 断面2cm×2cm、長さ40cm, ヤング率206.0GPA,ポアソン比0.3,荷重100N ,初等梁(m),ティモ(m) ,7.766990e-04,7.782136E-004 -まずccxプログラムと理論式で値を求める。 分割数はnx=10,ny=10 C3D8要素 ,nz,たわみmax(m),相対誤差(初等梁)(%),相対誤差(ティモ)(%),要素数 ,300,7.727840e-04,0.504060,6.97700E-01,30000 ,500,7.737010e-04,0.385996,5.79867E-01,50000 ,1500,7.741640e-04,3.26387E-001,5.20372E-001,150000 たわみの単位はメートル。 -salomeで作ったのを解く salomeでモデリング・unv生成、c3d4unvでinpに変換してccxで解く。 C3D4要素。 ,length,たわみmax(m),要素数,相対誤差(初等梁)(%),相対誤差(ティモ)(%) ,40.00,9.204240e-05, 128,8.81495E+01,8.81726E+01 ,20.00,4.252120e-04, 834,4.52540E+01,4.53605E+01 ,10.00,6.051090e-04, 3910,2.20922E+01,2.22438E+01 ,4.000,7.095070e-04, 21091,8.65097E+00,8.82876E+00 ,3.000,7.111990e-04, 21614,8.43313E+00,8.61134E+00 ,2.000,7.560170e-04,114338,2.66281E+00,2.85225E+00 ,1.800,7.553120e-04,109053,2.75358E+00,2.94284E+00 ,1.599,7.552970e-04,107468,2.75551E+00,2.94477E+00 ,1.300,7.596920e-04,201254,2.18965E+00,2.38001E+00 ,1.100,7.705140e-04,789275,7.96319E-01,9.89394E-01 ,1.000,7.706610e-04,817866,7.77393E-01,9.70505E-01 ,0.900,7.706280e-04,787666,7.81641E-01,9.74745E-01 ,0.770,7.706140e-04,789458,7.83444E-01,9.76544E-01 ,0.700,7.706920e-04,839416,7.73401E-01,9.66521E-01 ,0.600,7.724360e-04,2899514,5.48861E-01,7.42418E-01 ,0.500,ERROR,6119669,allocating memory ,corei5 メモリ16GBで発生 --現在のPCスペックでは、このモデルにおいてはlengthを0.600以下にするとメッシュ分割に時間がかかりすぎる(分割したところでccxで解ききれるかもわからない)ため、この辺りが限界。すべてC3D4でメッシュを切る場合、要素数が多くないと誤差が大きい。 *卒論日誌 [#d6e4e610] ,日付,作業時間,内容 ,4/11 , 3 ,ブライドタッチ ,4/15 , 3 ,vi操作 ,4/25 , 3 ,vi操作 ,5/13 , 4 ,fortran操作 ,5/23 , 4 ,fortran操作 ,5/30 , 4 ,ccx ,6/3 , 4 ,ccx ,6/10 , 3 ,ccxとSalome ,6/13 , 4 ,ccxとSalome ,6/14 , 4 ,Salome ,6/20 , 4 ,ccxとグラフ作成 ,6/24 , 4 ,ccxで箱型断面橋 の解析 ,6/28 , 4 ,箱型断面橋の解析、グラフ作成 ,7/1 , 4 ,ccxで木橋の解析 ,7/4 , 4 ,ccxで箱型断面の解析、グラフ作成 ,7/8 , 4 ,ccx ,7/10 , 9 ,継手あり破壊実験（能代) ,7/11 , 5 ,ccx ,7/12 , 5 ,gnuplot、グラフ ,9/27 , 5 ,断面の解析、tex ,10/3 , 5 ,モデルの解析プログラム作成、tex ,10/4 , 4 ,解析プログラム作成、グラフ作成 ,10/7 , 6 ,tex、解析プログラム ,10/8 , 4 ,tex ,10/9 , 4 ,tex、解析プログラム ,10/10 , 4 ,tex、解析プログラム ,10/17 , 5 ,解析プログラム ,10/18 , 4 ,解析プログラム ,10/22 , 4 ,解析プログラム ,10/23 , 5 ,解析プログラム ,10/25 , 5 ,解析プログラム補剛材あり ,10/29 , 5 ,解析プログラム補剛材あり ,10/31 , 5 ,解析プログラム補剛材あり ,11/1 , 5 ,解析プログラムの節点補剛材あり ,11/5 , 5 ,解析プログラムの要素補剛材あり ,11/6 , 5 ,解析プログラムの要素補剛材あり ,11/7 , 7 ,プレストレス鋼棒実験 ,11/13 , 6 ,解析プログラム補剛材なし作成 ,11/14 , 7 ,プレストレス鋼棒実験 ,11/22 , 6 ,解析プログラム修正・プログラム作成 ,11/25 , 6 ,プログラム、ccxによる解析 ,11/27 , 7 ,補剛材有り上部の鋼板無しプログラムの修正・節点 ,11/28 , 6 ,補剛材有り上部の鋼板無しプログラムの修正・要素 ,11/29 , 5 ,解析・データ整理 ,12/1 , 4 ,プログラムのヤング率変更・解析・グラフ ,12/2 , 7 ,プログラム・計算 ,12/3 , 5 ,プログラム ,12/4 , 6 ,プログラム製作 ,12/5 , 5 ,プログラム ,12/6 , 5 ,プログラム製作・データ整理 ,12/10 , 5 ,プログラム製作 ,12/13 , 5 ,グラフ作成・ccxによる解析・プログラム製作・修正 ,12/15 , 8 ,グラフ作成・データ整理・ifの活用(節点操作) ,12/16 , 5 ,グラフ・データ修正・プログラム ,12/17 , 8 ,データ整理・解析・グラフ修正・プログラム ,12/18 , 7 ,座屈プログラム製作修正・節点を増やすnaraxの活用 ,12/20 , 8 ,片持ち梁のグラフの間違いを修正、3点曲げと片持ちでデータ取得・グラフ整理 ,12/21 , 5 ,データ整理・プログラムで値の取得 ,12/22 , 10 ,中間発表用のスライド製作・データ整理 ,12/23 , 5 ,中間発表のためのデータ整理・グラフ更新 ,12/26 , 8 ,スライド更新・データ整理・発表練習 ,1/9 , 6 ,解析・今後の方針 ,1/10 , 11 ,穴開きプログラム節点 ,1/14 , 4 ,節点 ,1/16 , 12 ,プログラム ,1/21 , 7 ,プログラム・スライド更新 ,1/23 , 4 ,概要・スライド更新 ,1/24 , 13 ,プログラム・概要・スライド更新 ,1/26 , 7 ,卒論概要・スライド更新 ,1/27 , 12 ,卒論概要・スライド更新 ,1/29 , 13 ,卒論概要・スライド更新 ,1/31 , 8 ,発表練習・スライド更新 ,合計 ,412 *プレストレス木箱桁補剛材関連 [#k10ddb42] **FEMと穴開きモデルと2013/7/10、2013/2/13-14の実験との比較 [#c7196777] http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/FEM-jtest-tugite130710.png 穴開きモデルは木材部分の(saika)のたわみの平均値を取った。 **穴開きモデルとの比較 [#h20ddd71] ***補剛材あり鋼板あり穴開き [#h82d715a] http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/tugite130710kazyuu-tawami3ten-anaari.png -y軸:荷重(kN) x軸:たわみ(mm) 緑線-補剛材あり鋼板あり 青線-補剛材あり鋼板あり穴開きモデル 赤線-実験値(20130710) ***補剛材なし穴開き [#q7cf5898] http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/j-test1-hogounasianaari.png -y軸:荷重(kN) x軸:たわみ(mm) 緑線-補剛材なし穴なしモデル 青線-補剛材なし穴開きモデル 赤線-実験値(20130219.補剛材なし1回目) **12/20 ccxによる有限要素解析と2/17,18の実験との比較(3点曲げ) [#id7e1835] FEM鋼板有り、無しと2/17,18の実験結果 http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/FEM-zikken0218.png -y軸:荷重(kN) x軸:たわみ(mm) -FEMの方が固いので、座屈などの要素を取り入れて実験値に近づくかどうか -実際の実験で使われるような穴開きモデルを作ってccxで解析するなどが今後の課題 上のグラフに7月の実験結果も入れたグラフ http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/FEM-zikken0218-zikken07.png -y軸:荷重(kN) x軸:たわみ(mm) **12/17 [#e4c458cc] http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/j-test20130217-2.png -y軸:荷重(kN) x軸:たわみ(mm) --k2のdata/tugite130218内にある以下の補剛材なしの試験データとFEMの補剛材なしを比較する --( J-TEST-hs.xlsx 2/18 補剛材あり J-TEST-1.xlsx 2/19 1回目、補剛材なし J-TEST-2.xlsx 2/19 2回目、補剛材なし --) **データ類 [#sb531464] ***補剛材有り鋼板有りオンサイト3点曲げ座屈プログラム [#u5416da1] (&link(onsih3zakutu.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/onsih3zakutu.f90)) ***ccxで3点曲げで解いたもの [#bcf0ba81] ccxで荷重100kNを載荷した時のたわみ値を取得。ティモシェンコ理論値はkeionsite.f90類で取得。 http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/tugite130710kazyuu-tawami3ten.png -y軸:荷重(kN) x軸:たわみ(mm) 緑:補剛材有り鋼板有り 青:補剛材有り鋼板無し 紫:補剛材無し 赤:実験値 -補剛材あり鋼板有り座屈プログラム(&link(onsih3zakutu.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/onsih3zakutu.f90)) (&link(onsih3zakutu.dat,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/onsih3zakutu.dat)) ,,たわみ(m) ,,-5.1821E-03 ,,-5.8930E-03 ,平均,-5.5375E-03 -補剛材あり鋼板あり穴開き たわみ -9.85139050E-03 -補剛材なし穴開き -1.50161441E-02 平均点は : 2.09150003E-05 -補剛材あり鋼板有り(&link(onsih31220.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/onsih31220.f90)) ,たわみ(m) ,ティモシェンコ理論値(3点曲げ) ,6.82379538E-03 ,7.53422878067068914E-003 -補剛材あり鋼板無し(&link(onsihn31220.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/onsihn31220.f90)) ,たわみ(m) ,ティモシェンコ理論値(三点曲げ) ,8.35126638E-03 ,8.50147390384949769E-003 -補剛材無し(&link(onsi3ten1220.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/onsi3ten1220.f90)) ,たわみ(m),ティモシェンコ理論値(3点曲げ) ,8.65473226E-03,9.13542523676829851E-003 -ccxで解析したたわみ値とティモシェンコの式の理論値との比較(3点曲げ) http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/ccxtawami-timo3ten-2.png -y軸:荷重(kN) x軸:たわみ(mm) ***ccxで片持ち梁で解いたもの (上記の3点曲げに比べてellをell/2,荷重をp/2にしたもの) [#g4a472e5] http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/tugite130710kazyuu-tawamikata.png -y軸:荷重(kN) x軸:たわみ(mm) 緑:補剛材有り鋼板有り 青:補剛材有り鋼板無し 紫:補剛材無し 赤:実験値 -補剛材あり鋼板有り(&link(onsihogou1220.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/onsihogou1220.f90)) たわみ(mm) 6.87151263E-0 -補剛材あり鋼板無し(&link(onsihnasi1220.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/onsihnasi1220.f90)) たわみ(mm) 7.97978882E-0 -補剛材無し(&link(onsi21220.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/onsi21220.f90)) たわみ(mm) 8.60536564E-0 *モデル [#ac8f423b] **補剛材あり鋼板あり [#o9d69631] http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/Screenshot.png -剛性 118148936.16424280N・$m^2$ -σ木= 5031110.0908723157 Pa (&link(keionsited1.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/keionsited1.f90))にて取得 ,断面,f90孔なし,f90孔ありティモ,ティモ孔 ,補剛材あり,7.5342E-003,7.9970E-003 ,補剛材なし,,9.7931E-003 **補剛材あり鋼板なし [#vfa4117c] http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/Screenshot-1.png 剛性 104664411.84623867 N・$m^2$ σ木=6202551.3139656382 Pa (&link(keionsited2.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/keionsited2.f90))にて取得 **補剛材なし [#t1255234] http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/Screenshot3.png 剛性 96481417.91999998N・$m^2$ σ木=7774913.5449429713 Pa (&link(keionsited3.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/keionsited3.f90))にて取得 *目標 [#v16260fa] グラフの荷重を100kNまでにする **1/17まで [#w71157e1] nzhaji→ある nzhaji2→ない kの解決 **1/10まで [#pa48b6c8] 座屈プログラム完成 **12/27まで [#wa360eed] 有限要素解析での補剛材ありなしと2/17,18の補剛材ありなしをグラフで比較する 座屈プログラムで解析 データ整理 **12/20まで [#y0994c03] 座屈プログラム 過去の実験（2/18,19)で得たデータを補剛材有り、無しと比較しグラフを描く *パソコンゼミ [#p8921102] **7/29 [#l1670e1e] 課題:実験データからグラフを書く。※縦軸に載荷荷重をとるように -一例として下のグラフ -縦軸を載荷荷重、横軸を変位計の測定値に直してみました。 http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/ch7p.png **7/15 [#p8bb6892] 相対誤差とh2のグラフを描く。 ***木材と鋼板のモデル [#xafb036e] -&link(onsitekata3.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/onsitekata3.f90)で解析値を取得し、keionsite.f90を一部変更して理論値を取得。 -b1=1.2d-1!角材の幅 -b2=9.0d-3!鋼板の幅 -h1=9.0d-3!角材の高さ -h2=h2+5.00d-2!鋼板の高さ -h3=h2-h1*2!間の鋼板の高さ http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/moku-kou-graph.png -y軸:たわみ,x軸:h2[鋼板の高さ] 赤:理論値、緑:解析値 -相対誤差とh2のグラフ http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/mokukou-gosa-graph.png -y軸:相対誤差(%),x軸:h2(m) ***どっちも鋼板のモデル [#y95dee02] -onsitekata2.f90で解いたデータを解析値とし、hakotakasa.f90で理論値を取得。 -b1=1.2d-1!角材の幅 -b2=9.0d-3!鋼板の幅 -h1=1.2d-1!角材の高さ -h2=h2+5.00d-2!鋼板の高さ -h3=h2-h1*2!間の鋼板の高さ http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/kou-graph.png -y軸:たわみ,x軸:h2 赤:理論値, 緑:解析値 -相対誤差とh2のグラフ http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/kou-gosa-graph.png -y軸:相対誤差(%),x軸:h2(m) **7/1までの課題 [#q9f88bb7] 滝田さんのプログラムonsitekata2.f90を使い、たわみを解析してグラフ作成。 オンサイト木橋片持ち梁のモデルで、材料の寸法のh2(鋼板の高さ)を変化させてくれる。 要素分割はnz=500 ,h2,たわみ ,250,2.0820E-3 ,300,1.3374E-3 ,350,9.4549E-4 ,400,7.1509E-4 ,450,5.6755E-4 ,500,4.6677E-4 ,550,3.9444E-4 ,600,3.4045E-4 ,650,2.9887E-4 ,700,2.6602E-4 ,750,2.3951E-4 ,800,2.1773E-4 ,850,1.9956E-4 ,900,1.8419E-4 ,950,1.7106E-4 ,1000,1.5971E-4 http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/tawami.png -y軸:たわみ,x軸:h2 **6/24までの課題 [#sfe4cd8c] 私はnx1=2,nx2=3,ny1=6,ny2=8としてnzを変えてグラフを作成 ***応急橋片持ち [#i1dd6b13] ,nz,たわみ ,2,3.5076E-3 ,10,5.2700E-3 ,20,5.3724E-3 ,30,5.3948E-3 ,40,5.4037E-3 ,50,5.4083E-3 ,100,5.4154E-3 ,200,5.4177E-3 ,250,5.4180E-3 ,300,5.4182E-3 ,400,5.4183E-3 ,500,5.4184E-3 ,理論値,6.11052E-3 http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/ss.png -y軸:たわみ(E-3),x軸:nzの分割数 ***オンサイト木橋片持ち [#g53965f4] ,nz,たわみ ,2,1.9414E-3 ,10,5.4045E-3 ,20,5.7565E-3 ,30,5.8336E-3 ,40,5.8631E-3 ,50,5.8778E-3 ,60,5.8863E-3 ,100,5.9000E-3 ,200,5.9073E-3 ,500,5.9102E-3 ,1000,5.9108E-3 ,理論値,6.5039E-3 http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/ss2-2.png -y軸:たわみ(E-3),x軸:nzの分割数 **6/17までの課題 [#xf7c078d] 荷重は1MNとする。 6/10の課題での梁断面正方形 (b=0.040m ,h=0.040m,ell=0.200mのアクリル材)のモデルを考える。 たわみの理論値は0.0042967m=4.2967mmである。 ccx_2.3で出力されるdatファイルの中の大体一番大きそうな値をたわみ値にした。 ,Length,Nudes ,たわみ値(mm) ,0.030 ,156 ,3.0487 ,0.020 ,169 ,3.0505 ,0.015 ,660 ,3.7398 ,0.010 ,931 ,3.7085 ,0.008 ,1104 ,3.8722 ,0.006 ,4128,4.0899 ,0.005 ,4413,4.0920 ,0.004 ,5476,4.1287 ,0.003 ,20935,4.1846 ,0.002 ,34325,4.2114 ,0.001 ,181666,4.2321 -たわみとLengthのグラフ 見やすさを考慮しx軸は反転させている。右に行くほどLengthは小さい値になる。 http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/suuti1.png --y軸:たわみ値(mm),x軸:Length -たわみとNudesのグラフ 見づらかったので対数グラフにしました。 http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/nude0.png --y軸:たわみ値(mm),x軸:Nudes(要素数) ---- Lengthを小さくしてNudesを多くした方が理論値に近づくっぽい。 SalomeでLengthを入力するときは値が小さくなるなら2.0e-01とかのe-XXで入力する方が安心。 **6/10までの課題 [#b35ab7ce] b（幅）、h（高さ）、ell（長さ）とする。 また、荷重は1MNとする。 ***長方形（縦長） [#t3f9b8fe] b=0.020m ,h=0.040m,ell=0.200mのアクリル材 せん断補正係数によるたわみを計算したたわみ＝0.008593m http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/tatezk.png -nx=10,ny=10でnzを10～200まで変化させたグラフ -縦軸：たわみ(mm)、横軸：nz ***正方形 [#rea2f9cf] b=0.040m ,h=0.040m,ell=0.200mのアクリル材 せん断補正係数によるたわみを計算したたわみ＝0.0042967m http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/seihozk.png -nx=10,ny=10でnzを10～200まで変化させたグラフ -縦軸：たわみ(mm)、横軸：nz ***長方形(横長) [#cb7d032f] b=0.080m ,h=0.040m,ell=0.0200mのアクリル材 せん断補正係数によるたわみを計算したたわみ=0.0021483m http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/yokozk.png -nx=10,ny=10でnzを10～200まで変化させたグラフ -縦軸：たわみ(mm)、横軸：nz -ny,nzを増やしたほうがたわみ値は手計算の値に近づくのではないか。しかし、ny,nzを増やしていくと計算に時間がかかり、nx=10,ny=40,nz=200以上だと私のパソコンではものすごく時間がかかるようだったので今回はnx=10,ny=10としてnzを変化させていった。 *コマンド [#hf4c62fe] --( cp /media/Transcend/150416_1.CSV ~/zikken/ vi 150416_1.CSV ./syokignu2<150416_1.CSV>140416_1s.csv gnuplot --) 2015年度 以下は2014年度 --( cp /media/Transcend/140410_1.CSV ~/zikken/ vi 140410_1.CSV ./syokignu2<140410_1.CSV>140410_1s.csv gnuplot --)
