日付 | 時間帯 | 内容 | 立会 |
4/10 | 3 | パソコン操作 | |
4/14 | 3 | パソコン操作 | 斉藤 |
4/16 | 4 | パソコン操作 | 斉藤 |
4/17 | 4 | パソコン操作 | |
4/20 | 5 | vi操作 | 斉藤 |
4/21 | 5 | vi操作 | |
4/24 | 5 | fortran | |
4/28 | 5 | fortran | |
5/8 | 5 | fortran | |
5/10 | 5 | fortran | |
5/15 | 5 | salome | |
5/17 | 10 | salomeとccx解析 | 斉藤 |
5/20 | 8 | 構造力学 | |
5/22 | 5 | 構造力学 | |
5/25 | 6 | 構造力学 | |
6/1 | 8 | salomeとccx解析・木材直交異方性を解く。グラフ作成 | 斉藤 |
6/2 | 5 | gnuplotでグラフ作成・Latex |
6/3 | 6 | Latex・グラフ取り込み。 |
6/4 | 5 | 構造力学 | |
6/5 | 10 | Latex 課題 | |
6/6 | 7 | 構造力学 | |
6/9 | 6 | salomeとccx解析・2種類の材料を解く | 斉藤 |
6/12 | 5 | salomeとccx解析 | |
6/26 | 5 | salomeとccx解析 | |
6/30 | 10 | texレポート作成 | |
7/7 | 6 | salomeとccx解析、木杭のモデルを解く。 | 斉藤 |
9/3 | 12 | salomeとccx解析 | |
9/4 | 10 | salomeとccx解析 | |
9/5 | 8 | salomeとccx解析 | |
10/17 | 8 | ccx | |
10/19 | 5 | salome | |
10/20 | 8 | salome | |
10/24 | 10 | salome | |
10/25 | 11 | salomeとccx解析 | |
10/26 | 10 | salomeとccx解析 | |
10/27 | 8 | salomeとccx解析 | |
11/2 | 10 | salomeとccx解析 | |
11/5 | 9 | salomeとccx解析 | |
11/7 | 10 | salomeとccx解析 | |
11/9 | 8 | salomeとccx解析 | |
11/12 | 8 | salomeとccx解析 | |
12/17 | 10 | 中間 | |
1219 | 12 | 中間 | |
12/20 | 12 | 中間 | |
12/21 | 10 | 中間 | |
1/9 | 10 | 解析 | |
1/10 | 12 | 解析 | |
1/12 | 10 | 解析 | |
1/13 | 11 | 解析 | |
1/15 | 6 | 解析 | |
1/18 | 8 | 解析 | |
1/19 | 10 | tex | |
1/20 | 8 | tex | |
1/21 | 6 | 卒論概要 | |
1/22 | 6 | 卒論概要 | |
1/24 | 8 | 卒論概要 | |
1/25 | 8 | 卒論概要 | |
1/27 | 12 | 概要、卒論スライド | |
1/28 | 5 | 卒論スライド | |
1/29 | 5 | 卒論スライド | |
1/30 | 5 | 卒論スライド |
合計 | 460時間 |
・C3D4でメッシュ作成、ヤング率206e9Pa ,ポアソン比0.3 x50mm、y70mm、z1000mm、荷重500N
length | たわみ(max)[m] | 要素数 | 相対誤差[%] |
10 | 5.39236e-4 | 25341 | 4.79% |
5 | 5.58475e-4 | 146633 | 1.35% |
4 | 5.5862e-4 | 146141 | 1.32% |
3 | 1005503 | ||
2 | 5.642e-4 | 1053439 | 0.34% |
・たわみ理論値(初等梁)5.6611e-4[m]
lengthの値が小さくなるにつれてたわみ(max)の値が理論値に近づいていった。 相対誤差も同じように少なくなっていった。 よってlengthの値が小さいければ小さいほど理論値に近づくことがわかった。 ただ、lengthの値が小さいとmesh分割するのも細かくなるのでパソコンの容量的にとても 時間がかかるかerrorがおきて出来なかった。
・C3D8でメッシュ作成、ヤング率6GPa ,ポアソン比0.3 x50mm、y50mm、z100mm、荷重100N
初等梁理論値(\( v \)=\( \frac{Pl^3}{3EI} \))(単位 m) | ティモシェンコ梁理論値(\( v \)=\( \frac{Pl^3}{3EI}+\frac{Pl}{kGA} \))(単位 m) |
1.06667E-005 | 1.27467E-005 |
Number of Segments | element | たわみ最大値(m) | 相対誤差(初等梁)(%) | 相対誤差(ティモ)(%) |
5 | 125 | 1.1458E-5 | 7.418 | 10.11 |
10 | 1000 | 1.22179E-5 | 14.542 | 4.149 |
15 | 3375 | 1.23803E-5 | 16.065 | 2.874 |
20 | 8000 | 1.24379E-5 | 16.605 | 2.423 |
25 | 15625 | 1.2476E-5 | 16.844 | 2.222 |
30 | 27000 | 1.24762E-5 | 16.964 | 2.122 |
縦軸は相対誤差(%)、横軸はelementの値である グラフの緑線は初等梁、赤線はティモシェンコ梁である
初等梁理論値(\( v \)=\( \frac{Pl^3}{3EI} \))(単位 m) | ティモシェンコ梁理論値(\( v \)=\( \frac{Pl^3}{3EI}+\frac{Pl}{kGA} \))(単位 m) |
1.4062E-003 | 1.5187E-003 |
Number of Segments | たわみ最大値(m) | 相対誤差(初等梁)(%) | 相対誤差(ティモ)(%) |
5 | 1.31441e-3 | 6.53 | 13.5 |
10 | 1.46266e-3 | 4.01 | 3.69 |
15 | 1.4942e-3 | 6.254 | 1.616 |
20 | 1.5056e-3 | 7.1 | 0.86 |
25 | 1.51098e-3 | 7.45 | 0.51 |
30 | 1.5139e-3 | 7.66 | 0.32 |
40 | 1.5168e-3 | 7.862 | 0.128 |
50 | 1.5180e-03 |
縦軸は相対誤差(%)、横軸はNumber of Segmentsの値である グラフの赤線は初等梁、緑線はティモシェンコ梁である Number of Segmentsの値が増えると要素数が増える。
v=5.2358E-6
測定値 | 相対誤差 |
たわみmax=3.43999E-6 | 34.299% |
たわみ平均=1.7551E-6 | 66.479% |
x=120mm,y=50mm,z=150mmの梁を同一の寸法で、
初等梁理論値v=4.368932E-7 ccx測定値v=4.03814E-7 相対誤差 ~%
木材初等梁理論値v=1.1300853E-5 ティモシェンコ梁理論値v=1.695128E-5 ccx測定値v=1.04452E-5 相対誤差(初等梁) 相対誤差(ティモ)
理論値(\( v \)=\( \frac{Pl^3}{3(E木I木+E鋼I鋼)} \))(単位m)=2.208859E-6
要素数 | たわみ平均(m) | 相対誤差 |
57162 | 2.21546E-6 | 0.229 |
7957 | 2.25053E-6 | 1.887 |
3214 | 2.1255E-6 | 3.774 |
2222 | 2.0906E-6 | 5.354 |
&link(wdunv.f90,http://str.ce.akita-u.ac.jp/~gotouhan/j2014/r_watanabe/wdunv.f90)のプログラムの中の値を変えて解く。
木橋希望で