編集(管理者用) | 編集 | 差分 | 新規作成 | 一覧 | RSS | FrontPage | 検索 | 更新履歴

工藤のCalculiX関係 -

目次

鋼板挿入を三本にしたときのモデル化

上下に1本ずつ鋼板挿入した場合のモデル化

★全解析・・・・・zgousei1.p(ver.1)
★載荷方法・・・・面積載荷
★断面・・・・・・6*12*108cm
★保存場所・・・・zgousei/youso

ny、nzの値を固定し、nxの鋼板の要素数(nxk)を増やした場合のy方向変位

zgouseip.f(ver.1)

★全解析・・・・・zgouseip.f(ver.1)
断面・・・・・・・6*12*108
★データを取る箇所
・変位・・・図心:nx*ny*nz-{(ny/2)*(nx+1)+(nx/2)} nx,nyは偶数
★保存場所・・・zgouseip/youso
＊ny、nzの値を固定し、nxの要素を増やした場合のy方向変位
nx*16*120
nxnxknynyky方向変位
1021640.0103657
1241640.0103762
1461640.0103834
1681640.0103886
18101640.0103927
20121640.0103959

＊nx、nzの値を固定し、nyの要素を増やした場合のy方向変位

12*ny*120
nxnxknynyky方向変位
1241240.0103191
1241460.0102419
1241680.00868550
12418100.00866173
12420120.00864131
12422140.00862435

ckzkata.f(ver.1)全節点載荷

★全解析・・・・・・ckzkata.f(ver.1)
★載荷方法・・・・全節点載荷(p,2p,4p)
★解析法・・・・・弾性解析
★材料定数・・・・E=200GPa,ν=0.33333
★断面・・・・・・・・2*2*20cm
★データを取る箇所
・変位・・・図心：nx*ny*nz-{(ny/2)*(nx+1)+(nx/2)} nx,nyは偶数
・応力・・・図心軸とx軸との交点：(nx/2)+1
★備考・・・・・・手計算で応力が238MPaになるように0.0016MNを載荷
      ・・・・・・応力では最大2.2%、たわみでは1.2%の誤差がある
★保存場所・・・・g3 2006/yasu/ckz
●y方向の変位(たわみ)
分割数nx*ny*nz/10nx*16*5012*ny*5014*ny*5012*16*nz
10.00117040
20.001573870.001529360.001529510.00147252
30.00154658
40.001584130.001572500.001572640.00157437
50.00158760
60.001586270.001581310.001581440.00159490
70.00159934
80.001587040.001584490.001584620.00160227
90.00160428
100.001587410.001585980.001586100.00160573
110.00160680
120.001587600.001586790.001586910.00160762
130.00160827
140.001587720.001587280.001587400.00160879
150.00160920
160.001587800.001587600.00158772
180.001587830.00158795
200.00158798
220.00158810
240.00158819
260.00158826
280.00158831
300.00158835
320.00158839
340.00158842

●z軸に掛かる応力

分割数nx*ny*nz/10nx*16*5012*ny*5014*ny*5012*16*nz
1170.090
2244.206245.569245.614221.356
3234.568
4242.148244.323244.359240.192
5242.953
6242.580245.317245.346244.752
7245.296
8242.791244.759244.783246.165
9246.672
10242.887244.154244.173246.909
11247.196
12242.953243.707243.722247.117
13247.614
14242.961243.119243.129247.682
15247.757
16242.981242.953242.961
18242.487242.493
20242.256
22242.004
24241.734
26241.757
28241.474
30241.484
32241.181
34241.193

ckzkata.f(ver.)中立軸載荷

★全解析・・・・・・ckzkata.f(ver.2)
★載荷方法・・・・中立軸載荷
★解析法・・・・・弾性解析
★材料定数・・・・E=200GPa,ν=0.33333
★断面・・・・・・・・2*2*20cm
★データを取る箇所
・変位・・・図心：nx*ny*nz-{(ny/2)*(nx+1)+(nx/2)} nx,nyは偶数
・応力・・・図心軸とx軸との交点：(nx/2)+1
★保存場所・・・・g3 2006/yasu/ckztyu

ー決定事項ー
要素数は12*16*120で決定。

●y方向の変位(たわみ)

要素数nx*ny*nz/10nx*16*5012*ny*5012*16*nz
10.00116823
20.001574080.001526610.00146982
30.00154378
40.001584310.001569700.00157156
50.00158480
60.001587450.001578510.00159212
70.00159658
80.001587230.001581690.00159953
90.00160156
100.001587590.001583180.00160303
110.00160412
120.001584800.001583990.00160494
130.00160502
140.001583940.001584480.00160614
150.00160657
160.001588000.00158480
17
180.001589040.00158503
19
200.001588080.00158519
220.00158530
240.00158539
260.00158546

●z軸にかかる応力

分割数nx*ny*nz/10nx*16*5012*ny*5012*16*nz
1170.457
2244.111245.029220.843
3234.041
4242.066243.762239.659
5242.417
6242.651244.758244.214
7244.758
8242.711244.207245.626
9246.132
10242.807243.608246.369
11246.656
12242.417243.165246.578
13247.074
14242.274242.581247.142
15247.217
16242.900242.417
17
18243.081241.955
19
20242.941241.726
22241.476
24241.208
26241.231

ckzkata.fの全節点載荷と中立軸載荷の比較(表を一つにまとめたもの)

●y方向の変位(たわみ)

載荷方法全節点中立軸全節点中立軸全節点中立軸
分割数nx*ny*nz/10nx*16*50nx*16*5012*ny*5012*ny*5012*16*nz12*16*nz
10.001170400.00116823
20.001573870.001574080.001529360.001526610.001472520.00146982
30.001546580.00154378
40.001584130.001584310.001572500.001569700.001574370.00157156
50.001587600.00158480
60.001586270.001587450.001581310.001578510.001594900.00159212
70.001599340.00159658
80.001587040.001587230.001584490.001581690.001602270.00159953
90.001604280.00160156
100.001587410.001587590.001585980.001583180.001605730.00160303
110.001606800.00160412
120.001587600.001584800.001586790.001583990.001607620.00160494
130.001608270.00160562
140.001587720.001583940.001587280.001584480.001608790.00160614
150.001609200.00160657
160.001587800.001588000.001587600.00158480
180.001589040.001587830.00158503
200.001588080.001587980.00158519
220.001588100.00158530
240.001588190.00158539
260.001588260.00158546
280.00158831
300.00158835
320.00158839
340.00158842

●z軸に掛かる応力

載荷方法全節点中立軸全節点中立軸全節点中立軸
分割数nx*ny*nz/10nx*16*50nx*16*5012*ny*5012*ny*5012*16*nz12*16*nz
1170.090170.457
2244.206244.111245.569245.029221.356220.843
3234.568234.041
4242.148242.066244.323243.762240.192239.659
5242.953242.417
6242.580242.651245.317244.758244.752244.214
7245.296244.758
8242.791242.711244.759244.207246.165245.626
9246.672246.132
10242.887242.807244.154243.608246.909246.369
11247.196246.656
12242.953242.417243.707243.165247.117246.578
13247.614247.074
14242.961242.274243.119242.581247.682247.142
15247.757247.217
16242.981242.900242.953242.417
18243.081242.487241.955
20242.941242.256241.726
22242.004241.476
24241.734241.208
26241.757241.231
28241.474
30241.484
32241.181
34241.193

図をwikiに貼り付ける練習

集成材の材料定数はこちら http://www.syuseizai.com/intensity.html

Powered by YukiWiki 2.1.2a
Modified by GOTOU Humihiko.

*