後藤資料

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

目次

座屈荷重

無次元化

１倍

周方向分割数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

４倍

周方向分割数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

１６倍

周方向分割数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が計算できなかった。
名前手順１手順２
Mesh1Assign a set of hypotheses2D AT
Mesh2Assign a set of hypothesesAlgorithmNetgen1D-2D
Mesh3Assign a set of hypotheses2D ATでsizeを下に２回クリック
Mesh4Assign a set of hypotheses2D ATでsizeを上に４回クリック
Mesh5Assign a set of hypotheses2D ATでsizeを上に５回クリック
Mesh6Assign a set of hypotheses2D ATでsizeを下に５回クリック

kadai2 mesh

名前手順１手順２
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で１つ１つ打ち込んでいく、そして打ち込んだら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/19木14:30~16:001.5後藤ブラインドタッチ
4/23月16:00~17:001後藤viの練習
4/26木14:30~16:001.5後藤viの練習
5/7月16:00~17:001後藤viの練習
5/10木14:30~15:301後藤salomeの練習
5/17木14:30~1515:301後藤salomeの練習
5/21月16:00~17:001後藤salomeの練習
6/4月16:00~17:001後藤salomeの練習
6/18月16:00~18:002後藤texの練習
6/20水14:30~16:302後藤cgxの練習
6/21木14:30~16:001.5後藤cgxの練習
6/25月16:00~17:001後藤texの練習
6/28木14:30~18:003.5後藤cgxの練習
7/2月16:00~17:001後藤cgxの練習
7/4水14:00~15:301.5後藤cgxの点載荷
7/5木14:30~16:001.5後藤gimpとオープンファイルでのポスター作成
7/9月16:00~17:301.5後藤cgxの練習
7/10火13:00~19:006後藤課題
7/11水12:00~19:007後藤課題
7/12木12:00~19:007後藤課題
7/13金12:00~19:007後藤課題
9/18火12:00~19:007後藤課題
9/19水12:00~19:007後藤課題
9/20木12:00~19:007後藤課題
9/21金12:00~19:007後藤課題
9/24月12:00~18:006後藤課題
9/25火12:00~18:006後藤課題
9/26水12:00~18:006後藤課題
9/27木12:00~18:006後藤課題
9/28金12:00~18:006後藤課題
10/3水11:00~22:0011後藤課題
10/4木11:00~17:006後藤salomeでのダイヤ作り
10/5金14:00~17:003後藤ダイヤ作り
10/9火13:00~18:005後藤ダイヤ作り、Mesh切り
10/10水13:30~17:304後藤cgx
10/11木15:00~18:003後藤freecad
10/12金13:00~18:005後藤cgx
10/16火13:00~18:005後藤変位からバネ定数を求める
10/17水13:00~18:005後藤freecad
10/18木13:30~17:304後藤freecad
10/19金14:00~17:003後藤freecad
10/29月13:00~18:005後藤モデル作り
10/30火13:30~17:304後藤cgx、ccxでの計算
10/31水13:30~19:306後藤cgx、ccx
11/1木16:00~19:003後藤cgx、ccx
11/5月13:00~18:005後藤ccx
11/6火13:00~18:005後藤ccx
11/8木15:00~17:002後藤ccx
11/9金12:00~17:005後藤ccx
11/12月13:00~17:004後藤salome
11/14水14:00~19:005後藤salome
11/15木15:00~18:003後藤ccx
11/20火13:00~18:005後藤ccx
11/21水13:00~17:004後藤ccx
11/22木15:00~18:003後藤ccx
11/27火14:00~18:004後藤ccx
11/28水14:00~18:004後藤平均の計算
12/3月14:00~17:003後藤平均の計算
12/4火14:00~19:305.5後藤バネ定数の計算、グラフ
12/5水13:00~18:005後藤バネ定数の計算、グラフ、ccx
12/6木15:00~18:003後藤円筒の計算
12/11火14:00~19:305.5後藤座屈荷重・座屈モード
12/14金14:00~19:005後藤ccx
12/17月14:00~18:004後藤プレゼン作成
12/18火10:00~20:0010後藤プレゼン作成
12/19水11:00~23:0012後藤プレゼン作成
12/20木9:00~13:004後藤プレゼン作成
1/10木14:00~18:004後藤座屈
1/11金11:00~18:007後藤座屈
1/15火12:00~18:006後藤座屈
1/16水11:00~18:007後藤座屈
1/17木15:00~18:003後藤座屈
1/18金9:00~17:308.5後藤卒論概要
1/21月9:00~24:0013後藤卒論概要
1/22火10:00~21:0011後藤卒論概要
1/23水10:00~22:0012後藤卒論概要
1/24木11:00~22:0011後藤卒論概要
1/25金12:00~20:008後藤卒論概要
1/28月10:00~22:0012後藤プレゼン作成
1/29火10:00~19:009後藤プレゼン作成
1/30水10:00~22:0012後藤プレゼン作成
1/31木10:00~22:0012後藤プレゼン作成
2/1金10:00~22:0012後藤プレゼン作成
2/2土10:00~22:0012後藤プレゼン作成
2/3日10:00~22:0012後藤プレゼン作成
2/4月10:00~22:0012後藤プレゼン作成
Powered by YukiWiki 2.1.3
Modified by GOTOU Humihiko.

*