後藤資料

井比の卒論日誌 - *座屈荷重

目次

座屈荷重

無次元化

1倍

周方向分割数468101214161820
40.05890.02980.01620.01100.008760.007710.007250.007010.00701
60.03710.03040.02410.02070.01860.01850.01940.01910.0187
80.02740.03130.02670.01800.01660.01610.01640.01660.0168
100.03800.03530.02270.01970.01720.01620.01470.01470.0150
120.03590.03200.03160.03280.03300.01710.01610.01580.0161
140.02940.02480.02560.02660.02700.02670.02400.02070.0188
160.02530.02060.02060.02140.02210.02250.02290.02010.0180
180.02130.01830.01860.01900.01940.01980.02000.02010.0178
200.01920.01590.01550.01540.01570.01600.01620.01640.0166

4倍

周方向分割数468101214161820
40.002190.002050.001880.002120.002360.002480.002510.002490.00249
60.003800.004190.004300.004600.005270.006030.006580.007000.00734
80.005700.005720.006640.006760.007510.008410.009390.001040.00114
100.01560.01010.007920.008690.009680.01080.01150.01270.0139
120.01990.01250.01070.01090.01210.01190.01340.01440.0149
140.01850.01650.01310.01210.01290.01430.01630.01590.0162
160.01750.01670.01550.01360.01380.01500.01540.01550.0157
180.01680.01590.01580.01550.01520.01540.01490.01470.0146
200.01620.01500.01460.01460.01480.01490.01460.01460.0146

16倍

周方向分割数468101214161820
40.878E-30.988E-30.001200.001420.001630.001790.001890.001960.00196
60.001980.001920.002230.002770.003510.004420.005380.006140.00673
80.004020.003270.003370.003870.004670.005760.007080.008090.00905
100.006790.005480.004870.005260.005870.006850.007830.008870.0100
120.007610.006750.006420.006610.007050.007560.008300.009280.0102
140.008560.007220.007390.007510.007790.008260.008890.009580.0102
160.009720.007760.007540.007900.008530.008960.009350.009780.0102
180.01040.008620.007890.007990.008370.008790.009240.009610.00990
200.01070.009390.008610.008490.008690.009010.009330.009670.00989

kakutyuu

後藤さんが手直ししたkakutyuuで円を作ろうとしたら形や変位もおかしくなったので元のプログラミングで円を作りました。

後藤ちゃちゃ12/12/11

r=pi*r/sin(th/2)/m !多角形の外周を半径rの円周と揃える場合
上の部分に/mを入れてみて下さい。

バネ定数

無次元化

周方向分割数468101214161820
412.918.420.220.921.221.421.521.521.6
68.4015.318.820.721.121.321.421.521.5
85.5912.817.617.920.921.221.321.421.5
103.9410.716.016.920.821.121.321.421.5
122.999.0915.116.320.721.121.221.321.4
142.428.1314.115.817.821.021.221.321.4
162.056.9013.210.617.520.921.121.321.4
181.786.1712.410.017.220.821.121.221.4
201.785.6511.69.6216.920.821.121.221.4

バネ定数4倍

周方向分割数468101214161820
41.387.5215.023.825.526.226.626.826.9
61.414.6511.420.523.625.326.126.526.7
80.9233.759.0914.621.624.025.426.026.5
100.7252.936.9912.219.622.624.425.526.1
120.5882.365.7110.318.021.323.524.925.7
140.4881.974.818.8513.420.022.624.225.3
160.4121.684.137.6912.018.821.723.524.8
180.3531.463.626.8510.917.720.822.724.3
200.3531.293.256.1710.016.619.922.023.9

バネ定数16倍

周方向分割数468101214161820
40.4393.3810.117.321.222.923.724.024.1
60.3001.946.1012.217.020.422.423.323.7
80.2481.394.248.7013.717.620.422.123.0
100.2091.083.186.7111.215.218.520.722.1
120.1740.8702.515.359.2613.316.819.421.1
140.1440.7142.054.417.7511.715.318.220.2
160.1200.5991.733.806.7110.414.117.119.4
180.09900.5101.503.315.929.3512.916.118.7
200.09900.4411.322.955.328.4712.015.217.9

バネ定数64倍

周方向分割数468101214161820
40.2702.458.2014.518.921.122.122.522.7
60.1701.294.659.6214.618.320.521.622.2
80.1380.8673.046.9011.315.418.420.321.4
100.1170.6622.205.139.0113.016.418.920.4
120.09800.5291.703.987.3011.214.717.519.5
140.08200.4331.373.236.179.7113.316.318.6
160.06850.3641.152.758.5512.115.3
180.05750.3131.002.407.6311.1
200.05750.2750.8936.90

Mesh

やっていて思ったがこれはmeshの切り方じゃなくてsizeしか変えてないと思う。 そしてMesh3と6が計算できなかった。
名前手順1手順2
Mesh1Assign a set of hypotheses2D AT
Mesh2Assign a set of hypothesesAlgorithmNetgen1D-2D
Mesh3Assign a set of hypotheses2D ATでsizeを下に2回クリック
Mesh4Assign a set of hypotheses2D ATでsizeを上に4回クリック
Mesh5Assign a set of hypotheses2D ATでsizeを上に5回クリック
Mesh6Assign a set of hypotheses2D ATでsizeを下に5回クリック

kadai2 mesh

名前手順1手順2
Mesh1-5Assign a set of hypotheses2D AT
Mesh6-10NETGEN1D-2Dsimple parameters
Mesh11-15TriangleComposite、Maxsize
Mesh16-20TriangleComposite、no segment30

変位

daiya2の値と6がなぜかマイナスとなりおかしくなった。
名前変位
daiya5.00786502379924037E-004
daiya11.02497500000000013E-005
daiya2-202360692.80000001
daiya31.54191334634864098E-005
daiya45.09399999373272260E-006
daiya55.00169003316841520E-006
daiya620943766.449999999

バネ定数

名前k
daiya199685893.1
daiya19756335520
daiya36485448760
daiya419630938380
daiya519993242150

$\frac{kl}{EA}$

名前$\frac{kl}{EA}$
daiya0.02769030883
daiya11.352904501
daiya30.8993328286
daiya42.722209055
daiya52.772449476

エラー

salomeでのfuseをしていると
Boorean operation aborted non valid shape result
 と警告が出てfuseできない部分がある。

モデル ダイヤカット salome

daiya.fをコンパイルして、吐き出したobjファイルのx,y,z座標をsalomeのpointで1つ1つ打ち込んでいく、そして打ち込んだらlineで繋ぎ、faceにし、fuseにして形にした。

課題

円筒の作り方

円筒 円筒を作りsalome上でmeshを切る。 create mesh→Assign a set of hypothesesで3D Automatic Tetraを選んでApply→computeしてメッシュを切ると多分ダイヤカットではないが、それっぽい形に切ることができる。 それをunv fileにExportしてGmeshでマージしてそれをinp形式でsaveする。 inpfileから節点番号と要素番号をinp形式にまとめ、ccxs6.f90をコンパイルして実行する。 実行するときに
./ccxs6>hoge.inp
として吐き出したいinp形式のfileにまとめる。 そしてcgxという流れになるがsalomeでmeshを3Dで切るとcgxにしたときに形が欠けてしまう問題が起きる。 そこで2Dtetraの方でmeshを切って計算する。 ただサイズを適当に作ったモデルをccxで走らせるとうまく行ったが、同じようにモデルのサイズを決めてccxで計算させると
*ERROR in e_c3d: nonpositive jacobian
        determinant in element         190
とエラーが生じてしまう。

ダイヤ

nami84.step

リンク

**-Linux関連-
mesh
unix
**-tex関連-
コマンド
グラフィック
-calculix関連-
コマンド
試す
inpファイルの説明
**-その他-
単位変換
参考
プログラム
salome
salome2
CAE
vi
gnuplot
blender
freecad

日付曜日開始〜終了作業時間立合作業内容
4/1914:30~16:001.5後藤ブラインドタッチ
4/2316:00~17:001後藤viの練習
4/2614:30~16:001.5後藤viの練習
5/716:00~17:001後藤viの練習
5/1014:30~15:301後藤salomeの練習
5/1714:30~1515:301後藤salomeの練習
5/2116:00~17:001後藤salomeの練習
6/416:00~17:001後藤salomeの練習
6/1816:00~18:002後藤texの練習
6/2014:30~16:302後藤cgxの練習
6/2114:30~16:001.5後藤cgxの練習
6/2516:00~17:001後藤texの練習
6/2814:30~18:003.5後藤cgxの練習
7/216:00~17:001後藤cgxの練習
7/414:00~15:301.5後藤cgxの点載荷
7/514:30~16:001.5後藤gimpとオープンファイルでのポスター作成
7/916:00~17:301.5後藤cgxの練習
7/1013:00~19:006後藤課題
7/1112:00~19:007後藤課題
7/1212:00~19:007後藤課題
7/1312:00~19:007後藤課題
9/1812:00~19:007後藤課題
9/1912:00~19:007後藤課題
9/2012:00~19:007後藤課題
9/2112:00~19:007後藤課題
9/2412:00~18:006後藤課題
9/2512:00~18:006後藤課題
9/2612:00~18:006後藤課題
9/2712:00~18:006後藤課題
9/2812:00~18:006後藤課題
10/311:00~22:0011後藤課題
10/411:00~17:006後藤salomeでのダイヤ作り
10/514:00~17:003後藤ダイヤ作り
10/913:00~18:005後藤ダイヤ作り、Mesh切り
10/1013:30~17:304後藤cgx
10/1115:00~18:003後藤freecad
10/1213:00~18:005後藤cgx
10/1613:00~18:005後藤変位からバネ定数を求める
10/1713:00~18:005後藤freecad
10/1813:30~17:304後藤freecad
10/1914:00~17:003後藤freecad
10/2913:00~18:005後藤モデル作り
10/3013:30~17:304後藤cgx、ccxでの計算
10/3113:30~19:306後藤cgx、ccx
11/116:00~19:003後藤cgx、ccx
11/513:00~18:005後藤ccx
11/613:00~18:005後藤ccx
11/815:00~17:002後藤ccx
11/912:00~17:005後藤ccx
11/1213:00~17:004後藤salome
11/1414:00~19:005後藤salome
11/1515:00~18:003後藤ccx
11/2013:00~18:005後藤ccx
11/2113:00~17:004後藤ccx
11/2215:00~18:003後藤ccx
11/2714:00~18:004後藤ccx
11/2814:00~18:004後藤平均の計算
12/314:00~17:003後藤平均の計算
12/414:00~19:305.5後藤バネ定数の計算、グラフ
12/513:00~18:005後藤バネ定数の計算、グラフ、ccx
12/615:00~18:003後藤円筒の計算
12/1114:00~19:305.5後藤座屈荷重・座屈モード
12/1414:00~19:005後藤ccx
12/1714:00~18:004後藤プレゼン作成
12/1810:00~20:0010後藤プレゼン作成
12/1911:00~23:0012後藤プレゼン作成
12/209:00~13:004後藤プレゼン作成
1/1014:00~18:004後藤座屈
1/1111:00~18:007後藤座屈
1/1512:00~18:006後藤座屈
1/1611:00~18:007後藤座屈
1/1715:00~18:003後藤座屈
1/189:00~17:308.5後藤卒論概要
1/219:00~24:0013後藤卒論概要
1/2210:00~21:0011後藤卒論概要
1/2310:00~22:0012後藤卒論概要
1/2411:00~22:0011後藤卒論概要
1/2512:00~20:008後藤卒論概要
1/2810:00~22:0012後藤プレゼン作成
1/2910:00~19:009後藤プレゼン作成
1/3010:00~22:0012後藤プレゼン作成
1/3110:00~22:0012後藤プレゼン作成
2/110:00~22:0012後藤プレゼン作成
2/210:00~22:0012後藤プレゼン作成
2/310:00~22:0012後藤プレゼン作成
2/410:00~22:0012後藤プレゼン作成