たわみδ(m)
| 分割数(nx.ny/2.nz/10) | nx*20*100 | 5*ny*100 | 5*20*nz |
| 1 | 3.82414E-04 | 4.07445E-04 | |
| 2 | 4.68018E-04 | 4.16297E-04 | 4.35141E-04 |
| 3 | 4.72480E-04 | 4.35032E-04 | 4.48660E-04 |
| 4 | 4.74400E-04 | 4.46627E-04 | 4.57106E-04 |
| 5 | 4.75477E-04 | 4.54794E-04 | 4.62870E-04 |
| 6 | 4.60947E-04 | 4.66992E-04 | |
| 7 | 4.65774E-04 | 4.70032E-04 | |
| 8 | 4.69658E-04 | 4.72327E-04 | |
| 9 | 4.72838E-04 | 4.74094E-04 | |
| 10 | 4.75477E-04 | 4.75477E-04 |
たわみδ(cm)
| 分割数(nx.ny/2.nz/10) | nx*20*100 | 5*ny*100 | 5*20*nz |
| 1 | 3.86670E-04 | 4.09447E-04 | |
| 2 | 4.79374E-04 | 4.20554E-04 | 4.37727E-04 |
| 3 | 4.79656E-04 | 4.39330E-04 | 4.51683E-04 |
| 4 | 4.79747E-04 | 4.50922E-04 | 4.60472E-04 |
| 5 | 4.79788E-04 | 4.59095E-04 | 4.66504E-04 |
| 6 | 4.65251E-04 | 4.70835E-04 | |
| 7 | 4.70080E-04 | 4.74037E-04 | |
| 8 | 4.73965E-04 | 4.76459E-04 | |
| 9 | 4.77147E-04 | 4.78326E-04 | |
| 10 | 4.79788E-04 | 4.79788E-04 |
たわみδ(m)
| 分割数(nx.ny/2.nz/10) | nx*20*100 | 5*ny*100 | 5*20*nz |
| 1 | 9.58186E-03 | 7.57345E-03 | |
| 2 | 9.75436E-03 | 9.70014E-03 | 9.09641E-03 |
| 3 | 9.76225E-03 | 9.73432E-03 | 9.45630E-03 |
| 4 | 9.76523E-03 | 9.74855E-03 | 9.59379E-03 |
| 5 | 9.76664E-03 | 9.75584E-03 | 9.66175E-03 |
| 6 | 9.76006E-03 | 9.70108E-03 | |
| 7 | 9.76272E-03 | 9.72638E-03 | |
| 8 | 9.76449E-03 | 9.74393E-03 | |
| 9 | 9.76574E-03 | 9.75679E-03 | |
| 10 | 9.76664E-03 | 9.76664E-03 |
たわみδ(cm)
| 分割数(nx.ny/2.nz/10) | nx*20*100 | 5*ny*100 | 5*20*nz |
| 1 | 9.59106E-03 | 7.57925E-03 | |
| 2 | 9.77531E-03 | 9.70910E-03 | 9.10378E-03 |
| 3 | 9.77544E-03 | 9.74325E-03 | 9.46412E-03 |
| 4 | 9.77549E-03 | 9.75745E-03 | 9.60185E-03 |
| 5 | 9.77551E-03 | 9.76473E-03 | 9.66999E-03 |
| 6 | 9.76894E-03 | 9.70947E-03 | |
| 7 | 9.77159E-03 | 9.73490E-03 | |
| 8 | 9.77336E-03 | 9.75257E-03 | |
| 9 | 9.77460E-03 | 9.76555E-03 | |
| 10 | 9.77551E-03 | 9.77551E-03 |
| 5点(h/l=10〜20) | 3点 | |
| 正方形 | 576MPa | 569.8MPa |
| 長方形 | 570MPa | 566.8MPa |
| k | 5点 | |
| 正方形 | 0.595 | 632MPa |
| 長方形 | 0.529 | 631MPa |
| 5点 | 3点 | |
| 正方形 | 10.23GPa | 10.055GPa |
| 長方形 | 10.22GPa | 10.055MPa |
| 5点 | 3点 | |
| 正方形 | 242MPa | 331MPa |
| 長方形 | 237MPa |
誤差がかなりあります。
というか5点の方が悪いようです。
線形回帰で
(h/L)^2が大きいほうに引っ張られて、
傾きが大きく出てしまいます。
で、１/E’を求め、ｙ軸に、
ｘ軸には(h/L)^2をプロットして線形回帰です。
今、手元にデータがないのですが、Eもマイナスとかになってました。
で、E（切片）、G（傾き）を逆算しました。
7PL^3/(769E'I)=7PL^3/(768EI)+73PL/(256kGA) でしたっけ?
768I/(7PL^3)を辺々かけて、
1/E'=1/E+768*73I/(7*256kGAL^2)
=1/E+219h^2/(7kGL^2)
とかそんな感じにならないすか(いま、子供に邪魔されながらなんで計算 間違ってるかも知れませんが)。
28はどこから出てきたのかな。
I=bh^3/12とおけばそうなりませんかね？？
この式はティモシェンコ梁のたわみと初等梁のたわみを等値した式な訳ですが、 FEMとティモシェンコとの相対誤差がじゅうぶん小さいなら、 ティモシェンコの式から逆算したkが0.6ぐらいに落ち着く範囲のL/hのところで k=0.6を使ってGを推定すれば、それなりにFEMのGと近い値が出そうな気がするのですが。
Eがマイナスってのもおかしいですし、
tyokuの結果は ここ にありますが。
L/ｈが小さいところでひっぱられるって感じはします。 それを除外してやってみてもいいですかね？？
荷重=0.0005MN
要素分割は5*20*100
サイズは0.02×0.02×0.4(cm)
| ccx | tyokum | 梁理論 | ティモシェンコ |
| 5.01677E-03 | 4.91779E-03 | 4.37500E-03 | 4.86384E-03 |
荷重=0.0005MN
要素分割は5*20*100
サイズは0.02×0.02×0.4(cm)
| ccx | tyokum | 梁理論 | ティモシェンコ |
| 2.23288E-4 | 2.19804E-4 | 2.18750E-04 | 2.23313E-04 |
ckhkata.fを改造してみる。
とりあえず走らせてみて
手計算と比較してみる。
大黒屋のパソコンでは。。。
nx=40,ny=40,nz=120での計算速度も無理らしい
nx=30,ny=30,nz=120での計算速度も無理らしい
nx=20,ny=30,nz=120での計算速度も無理らしい
nx=20,ny=15,nz=120での計算速度は20分くらい
nx=20,ny=18,nz=120での計算速度は25分くらい
nx=20,ny=19,nz=120での計算速度は30分くらい
nx=20,ny=20,nz=120での計算速度は40分くらい
というか、リブートしてすぐに計算したら、こないだできなかった
nx=20,ny=20,nz=120もOKでした
