工藤の修論日誌 -

目次


ディレクトリ内の容量確認
fs

xfigマニュアル

xfig
後藤ちゃちゃ09/3/24: xfigのマニュアルをコピーして公開する場合は、 公開したURLをメールで作者に通知しないといけないのでは?

gnuplot

3Dカラー

折り畳み

 ★半径3cm
 ★高さ12cm固定
 ★厚さ0.01cm
 ★ヤング率69GPa
 ★共有/kudo/oritatami/oritata4、6、8、
 拘束条件
 上端xy拘束
 下端xyz拘束
画像

折り畳み20段の引張り時の座屈荷重(MN)

角数
468101214161820222426
02.94378E-023.38603E-024.40660E-023.55454E-02
0.011.72966E-031.87694E-032.22064E-033.15384E-03
0.19.30121E-031.07625E-029.91085E-031.31063E-02
0.21.94507E-022.58106E-023.24781E-024.24159E-02
0.33.51496E-022.96471E-024.47350E-028.54233E-02
0.45.27016E-023.43917E-024.73805E-021.08832E-01
0.57.28046E-024.14533E-024.88657E-021.14019E-01
0.69.39485E-025.17711E-025.22749E-021.14470E-01
0.71.10517E-016.86519E-025.84378E-021.13247E-01
0.81.01838E-011.02075E-017.30412E-021.11831E-012.25767E-01
0.99.05073E-021.05412E-011.23934E-011.31591E-012.13389E-01
0.958.43785E-029.61809E-021.29739E-011.96433E-012.66367E-014.33512E-01
0.991.00685E-011.15826E-011.56764E-012.67843E-014.45691E-015.56433E-015.81734E-015.85812E-016.36417E-016.84472E-017.56588E-018.58851E-01
1

折り畳み20段の圧縮時の座屈荷重(MN)

角数
468101214161820222426
0.015.51870E-045.74364E-048.49146E-04
0.12.28896E-032.90749E-032.16647E-032.25596E-03
0.25.76448E-039.31832E-038.82674E-037.92893E-03
0.31.16045E-022.04675E-022.13694E-021.88578E-02
0.41.96068E-023.26609E-023.85016E-023.50836E-02
0.52.80914E-024.20508E-025.59717E-025.47319E-02
0.63.55675E-024.90318E-026.11344E-027.43287E-02
0.74.18731E-025.55928E-026.64093E-028.78798E-02
0.84.71071E-026.33780E-027.26329E-028.85483E-029.02423E-02
0.94.98287E-026.89664E-027.97171E-029.23725E-021.02775E-01
0.954.60403E-026.83853E-028.04821E-028.99273E-021.02051E-011.09545E-01
0.991.93257E-024.71795E-026.55079E-027.77603E-028.67922E-029.49093E-021.03643E-011.11977E-011.07775E-011.04116E-011.06832E-011.11989E-01
12.77669E-035.49290E-039.22701E-031.39479E-021.96921E-022.65621E-023.46656E-024.46685E-025.58788E-026.78575E-028.05723E-029.05698E-02

折り畳み20段のばね定数

画像
 台形の中央
468101214161820222426
00.290.520.91.98
0.010.410.640.841.31
0.10.841.011.051.26
0.21.171.321.241.4
0.31.511.651.471.61
0.41.882.071.791.9
0.52.322.632.242.3
0.62.913.252.952.89
0.73.844.674.143.9
0.85.526.886.435.95.99
0.99.5312.2312.311.5411.13
0.9514.4919.2620.7620.5519.9219.820.48
0.9929.9238.1344.4649.4353.0355.5257.5560.0262.9266.6871.6181.12
1198.9204.36207.16208.87210.04210.92211.6212.04212.37212.79213.14213.45

台形の角

468101214161820222426
05.631.060.530.45
0.011.820.860.560.54
0.11.621.180.960.99
0.21.751.421.181.25
0.321.721.421.49
0.42.312.131.741.79
0.52.742.72.192.2
0.63.393.332.92.79
0.74.474.814.093.78
0.86.67.196.45.775.8
0.912.3413.2512.511.4310.89
0.9520.3721.8821.6920.7119.7219.4119.92
0.9943.2646.0749.2352.2754.5856.0857.5759.5962.165.4769.978.44
1220.06218.82218.25217.95217.77217.65217.57217.41217.18217.19217.19217.19

折り畳みプログラムの最低限の周方向の長さを出す計算式


 最低の長さ=(orih/cos(oh)*m)+(b*m/2)

 orih=h/n*cos(psi)
 b=(ue-sita)+(syu/m)
 sita=orih/cos(oh)
 ue=orih*tan(oh)
 oh=pi/2.-th
 psi=asin(dh/(h/n))

折り畳み

 ★半径3cm
 ★高さ12cm固定
 ★厚さ0.01cm
 ★ヤング率69GPa
 ★共有/kudo/oritatami/oritata4、6、8、10
 拘束
 上端xy拘束
 下端z方向、一点のx方向と、その対称な点のy方向

折り畳み20段の座屈荷重(MN)

引っ張り

角数
468101214161820222426
02.93072E-023.37846E-024.40172E-023.55455E-02
0.011.56980E-031.77812E-032.12183E-033.07778E-03
0.19.29769E-031.06469E-029.88641E-031.30776E-02
0.21.93881E-022.91287E-023.24674E-024.23851E-02
0.33.51071E-026.35070E-025.92944E-028.54028E-02
0.45.26318E-021.08404E-018.13602E-021.09071E-01
0.56.92713E-021.36902E-018.96667E-021.15981E-01
0.67.24543E-021.55433E-019.68851E-021.18702E-01
0.77.98070E-021.70036E-011.08196E-011.25789E-01
0.89.47991E-021.86855E-011.26732E-011.38364E-01
0.91.24190E-012.56533E-011.61918E-011.62782E-01
0.951.50298E-012.11678E-011.79669E-011.72216E-012.20548E-013.12412E-014.39393E-01
0.991.37738E-011.88969E-012.15737E-012.34231E-012.91331E-013.95956E-014.76250E-015.36143E-015.02131E-015.50999E-015.32395E-016.67961E-01
1

圧縮

角数
468101214161820222426
02.93072E-023.37846E-024.40172E-023.56145E-02
0.015.33856E-045.50767E-048.17830E-04
0.12.28902E-032.89655E-032.16141E-032.25221E-03
0.25.76271E-039.29864E-038.82751E-037.92635E-03
0.31.16054E-022.04505E-022.13685E-021.88478E-02
0.41.96103E-023.26538E-023.85053E-023.50554E-02
0.52.80904E-024.17122E-025.66936E-025.21030E-02
0.63.55684E-024.93104E-026.50769E-026.08982E-02
0.74.18667E-025.63681E-027.06608E-027.22538E-02
0.84.71084E-026.31481E-027.54070E-028.99139E-027.25777E-02
0.94.98307E-026.89513E-027.96288E-029.23639E-021.00254E-01
0.954.60190E-026.83915E-028.04664E-028.99227E-021.01868E-011.06556E-019.30233E-02
0.991.93257E-024.71794E-026.55078E-027.77601E-028.67919E-029.49089E-021.03852E-011.00893E-019.92306E-029.35945E-021.00310E-011.03686E-01
12.66845E-034.24695E-036.29202E-039.06017E-031.25576E-021.67409E-022.16751E-022.76871E-023.42327E-024.14293E-024.94297E-025.78355E-02

折り畳み20段のばね定数(MN/m)

 ※圧縮、引張共に変わらず
角数
468101214161820222426
00.290.520.91.98
0.010.410.640.831.26
0.10.810.9911.2
0.21.111.281.171.34
0.31.411.581.371.53
0.41.731.971.661.8
0.52.142.492.072.18
0.62.693.252.722.75
0.73.564.433.813.71
0.85.176.575.925.65
0.99.1111.8511.5111.12
0.9514.0918.9419.8719.98
0.9929.3637.6943.8648.5352.254.1856.3658.2661.3564.4769.678.62
1197.97203.68206.64208.45209.69210.61211.32211.91212.15212.58212.94213.26

* syutkotei.f

12/19

 シェルスプリクと
kaite.sh

 折り畳み
画像 画像

12/18

BUCKLE

inpファイルの中の
 *BUCKLE
 5,0.01
の0.01は要求する精度(default=0.01)って意味らしいけど、ここを
 *BUCKLE
 5,1
に書き換えれば荷重倍率が大きすぎても座屈荷重が出力される。これを座屈荷重が出力する荷重倍率に書き換えたものと比較しても、小数点以下6〜7桁は同じになる。

12/16

ダイヤカットに横荷重を載荷するプログラム

daiya4yka.f プログラムの中の荷重を書き換えて、コンパイルして実行する。横荷重の載荷節点番号と荷重が入力されたyokoka.inpが出来上がり、 (例) 自動でdaiya4yka.inpのCLOADの下に書き込まれる。。 daiya4yka.objはble6yka.f に通してください。

12/15

ダイヤカットに外圧を三角形の各頂点から載荷した時の座屈荷重

 載荷荷重=0.1KN
 ★2008/daiya
 モード図
画像

451015202530
45.1447E-021.9411E-039.2230E-047.9789E-047.7206E-047.6589E-04
52.6750E-02
101.4921E-014.2938E-03
151.3657E-011.6974E-03
201.1683E-011.0892E-03
258.8158E-029.0525E-04
307.3060E-028.3449E-04

12/10

ダイヤカットを引っ張ったときの座屈荷重(MN)

 載荷点を増やしたやつ
 載荷荷重=0.1KN
★共有bute/youso/dat

4567891011121314151617181920
44.4044E-023.0233E-022.2222E-031.7087E-021.7833E-02×××
51.0660E-024.2295E-022.8930E-022.1586E-021.9116E-022.2127E-02××
64.5968E-033.3143E-023.9433E-022.9927E-022.5177E-022.5891E-027.3997E-03×
73.2824E-031.3625E-024.9807E-023.8672E-023.1443E-023.2167E-023.5751E-02×
82.7439E-038.6218E-034.0195E-025.1631E-024.2115E-023.8159E-024.1260E-02×
92.6068E-036.9348E-032.4742E-025.9616E-025.6052E-024.8441E-024.8907E-02×
102.7675E-036.6193E-031.9459E-026.4907E-026.9662E-026.4534E-026.1090E-027.0953E-02
113.0800E-036.9001E-031.7640E-025.8376E-027.9629E-02×8.3974E-02
123.4805E-037.5155E-031.7276E-025.1766E-028.5779E-029.9775E-027.8259E-02×
134.0781E-038.7052E-031.8503E-024.9956E-029.2431E-021.0020E-011.0045E-01×
144.8896E-031.0450E-022.0793E-024.9310E-029.5603E-021.0393E-01××
156.0238E-031.3050E-022.4685E-025.3090E-021.0205E-019.9733E-02×
167.6525E-031.6778E-022.9938E-025.8961E-028.9817E-029.5842E-02×
179.9627E-032.1960E-023.7098E-026.6002E-028.5371E-028.9291E-02×
181.2809E-022.8204E-024.4910E-026.8561E-028.1032E-028.3595E-021.0722E-01
191.7550E-023.6697E-025.0584E-026.1520E-027.7846E-027.7908E-029.7021E-02
202.2750E-024.5614E-025.0755E-025.6773E-027.2208E-027.3112E-028.8470E-02

荷重載荷点を増やした ble6daiya2.f

12/4-2

ダイヤカットを引っ張ったときの座屈荷重(MN)

 載荷点は上端の要素の頂点だけ
 載荷荷重=0.1KN
 ★buckle/daiya/yousot/inpt/dat
画像

456789101112131415
44.41E-023.02E-022.18E-021.73E-021.77E-022.70E-02
51.07E-024.25E-022.88E-022.16E-021.91E-022.22E-02
64.60E-033.32E-023.96E-022.95E-022.52E-022.59E-027.41E-02
73.28E-031.36E-024.87E-023.88E-023.14E-023.04E-023.55E-02
82.74E-038.62E-034.06E-025.09E-024.23E-023.82E-024.14E-02
92.61E-036.93E-032.48E-024.57E-025.58E-024.85E-024.94E-02
102.77E-036.62E-031.95E-023.91E-025.02E-026.02E-026.12E-027.24E-02
113.08E-036.90E-031.76E-023.53E-024.36E-025.27E-026.39E-027.86E-02
123.48E-037.52E-031.73E-023.26E-024.06E-024.46E-025.71E-026.69E-02
134.08E-038.71E-031.85E-023.18E-023.76E-024.16E-025.15E-026.12E-02
144.89E-031.05E-022.08E-023.20E-023.78E-024.02E-024.83E-025.72E-026.33E-02
156.02E-031.31E-022.46E-023.35E-023.73E-023.96E-024.61E-025.40E-026.29E-02
167.65E-031.68E-022.99E-023.56E-023.95E-024.01E-024.56E-025.27E-025.88E-02
179.96E-032.19E-023.33E-023.79E-024.06E-024.10E-024.64E-025.00E-025.87E-02
181.28E-022.82E-023.75E-024.19E-024.52E-024.37E-024.83E-025.27E-025.83E-026.50E-02
191.73E-023.50E-024.17E-024.61E-024.74E-024.61E-025.04E-025.46E-026.01E-026.44E-02
202.36E-023.71E-024.50E-025.01E-025.30E-024.92E-025.33E-025.64E-026.03E-026.53E-02

12/4

ダイヤカットを圧縮させたときの座屈荷重(MN)
 ★buckle/yousoc/dat

画像

4567891011121314151617181920
41.9038E-031.7279E-031.6762E-031.7411E-032.0748E-032.3176E-032.5964E-032.6650E-032.7815E-032.7218E-032.6690E-032.5393E-032.4496E-032.3356E-032.2266E-032.1389E-032.0661E-03
52.1049E-032.4958E-032.4221E-032.4348E-032.6444E-032.9761E-033.3036E-033.5970E-033.8423E-033.9798E-034.0415E-034.0366E-034.0195E-033.9779E-033.9133E-033.8715E-033.8355E-03
62.8089E-033.3587E-033.0903E-033.1469E-033.2889E-033.5588E-033.8828E-034.2433E-034.6137E-034.9127E-035.1494E-035.3075E-035.4422E-035.5366E-035.5902E-035.6640E-035.7363E-03
73.2921E-033.5770E-033.8858E-034.0815E-033.8843E-034.1023E-034.3766E-034.7245E-035.1341E-035.5311E-035.9089E-036.2239E-036.5232E-036.7745E-036.9719E-037.1916E-037.4042E-03
84.2117E-034.1135E-034.2907E-034.5232E-034.7400E-035.0381E-034.9358E-035.2597E-035.6363E-036.0622E-036.4900E-036.9047E-037.3136E-037.6992E-038.0152E-038.3752E-038.7116E-03
95.9825E-035.3735E-035.2105E-035.2610E-035.4336E-035.6803E-035.9088E-036.2083E-036.3079E-036.6987E-037.1513E-037.5947E-038.0849E-038.5493E-038.9609E-039.2059E-039.3993E-03
101.0329E-029.0990E-039.7025E-035.9636E-036.0214E-036.1827E-036.4006E-036.6837E-037.0479E-037.4474E-037.9116E-038.0529E-038.4558E-038.6623E-038.8703E-039.0677E-039.2742E-03

4567891011121314151617181920
118.5098E-039.6599E-038.8314E-038.1633E-038.0335E-038.1662E-038.0594E-037.1799E-037.5274E-037.8833E-038.0675E-038.2752E-038.4526E-038.6128E-038.8327E-039.0434E-039.2507E-03
127.5746E-037.8853E-038.6984E-038.8752E-038.5336E-038.5183E-038.7598E-039.1050E-039.2680E-039.2609E-038.1095E-038.2565E-038.4517E-038.6215E-038.8225E-039.0121E-039.2083E-03
137.4008E-037.1651E-037.4817E-038.1566E-038.9482E-039.0580E-039.1591E-039.4649E-039.4659E-039.4281E-039.4422E-039.5139E-039.5890E-039.5356E-038.9817E-039.1454E-039.3106E-03
147.5923E-036.9518E-036.8711E-037.1737E-037.6565E-038.2367E-038.8196E-039.3463E-039.5760E-039.4790E-039.4369E-039.4634E-039.5014E-039.5754E-039.6946E-039.7881E-039.9083E-03
157.9506E-037.0884E-036.7259E-036.7705E-037.0299E-037.4261E-037.8926E-038.3597E-038.7823E-039.1495E-039.4734E-039.7321E-039.7215E-039.7576E-039.8643E-039.9232E-031.0008E-02
168.3392E-037.3791E-036.7854E-036.6171E-036.6910E-036.9180E-037.1868E-037.5706E-037.9755E-038.2775E-038.6610E-038.9112E-039.1766E-039.4198E-039.6780E-039.8880E-031.0043E-02
178.7529E-037.6917E-036.9546E-036.6331E-036.5215E-036.5847E-036.7728E-037.0277E-037.3106E-037.6074E-037.9095E-038.2073E-038.4660E-038.7134E-038.9678E-039.1695E-039.3820E-03
188.8081E-038.1426E-037.4904E-037.1600E-037.0013E-036.9455E-036.9960E-037.1210E-037.2611E-037.4286E-037.6199E-037.8089E-037.9747E-038.1641E-038.3573E-038.5066E-038.6672E-03
198.9405E-038.4546E-037.7465E-037.3504E-037.1313E-037.0112E-037.0020E-037.0675E-037.1644E-037.2903E-037.4478E-037.6204E-037.8011E-037.9376E-038.1323E-038.2558E-038.4011E-03
209.0507E-038.6090E-037.9794E-037.5468E-037.2615E-037.1081E-037.0381E-037.0291E-037.0731E-037.1623E-037.2770E-037.4203E-037.5647E-037.7083E-037.9078E-038.0311E-038.1890E-03

(12/2)ダイヤカット

比較モデル
半径3cm
高さ12cm
厚さ0.02cm
ヤング率69GPa

引張、圧縮共にばね定数は変わらない。(変位が同じ)

NLGEONを入力しても座屈荷重は変わらない

1125

氷結の実験

(11/20)

折りたたみ4段のの要素分割数を増やせるようにしたプログラム tubu4sub.f これもble6ori3.fで使えます。

(11/19)

折りたたみ4段のプログラムの高さ方向分割数を増やせるようにした tubu1sub.f これもble6ori3.f

(11/17)

こんな感じになりました画像

f

折りたたみ4段のプログラムをサブルーチンで整理した tubusub.f

(11/12)

創造工房で使う、長さとか与えてやると自動的にABDの長さが出るやつ。 sozokobo.f

(11/11)

座標出力うまくいったみたいです。稲荷くんのおかげ。。 rhen.f

(11/5)

ダイヤカット載荷、拘束を出力させるプログラム。両方保存しなおしてください。 ble6daiya.f daiya4bun.f

あとBUCKLEをいれても no convergence for the buckling factor; maybe no buckling occurs とエラーが出る。。荷重を変えても変わらず。 cgxでadd-displacementで描かせるとこんな感じ

(11/1)

段数を増やしてもue図のBの長さが一定になる。 プログラムの座標の出力を見直してみます。。

後藤メモ(08/10/31)

プログラムが4段ぶんしかないとしても、5段の場合の座標は計算できまそうですね。 例えば、創造工房の円筒折り紙構造の高さが19.2cmだとして、 初期高さをこれの14の4.8cmで5段に折りたい場合、 これの4段ぶんの高さは、4.8×45=3.84cmだから、 プログラム上は、折る前の紙の高さを19.2×45=15.36cmにしておいて、 初期高さを3.84cmで4段の折り紙構造にすれば、 台形の上底の長さを出せるのでは?

(10/29)ナオシカケ

ble6ori2.f iitubure.f

(10/21)

拘束は今まで通りで、さらに載荷節点をxy方向拘束。座屈モード1の時のcgxの図 画像

ble6dkata.f

ダイヤカット片側固定(一点xz拘束、その点に対象な点yz拘束、他z拘束)で軸方向に載荷した .objファイルをcalculix用に書き換えるプログラム ble6dkata.f。 円柱、ダイヤカットの1分割(daiya.f)の周方向分割数が偶数のときは0、奇数は1、ダイヤの4分割(daiya4bun.f)の時は2を20行目ぐらいにあるdaiya=のところに各自書き換えて。

(10/10)

おりたたみのプログラムを走らせるとバージョンによってエラーがでる1.6も7と同様 ccx1_5 ccx1_7

後藤ちゃちゃ(08/10/10)

稲荷の卒論日誌参照。

リンク

円筒 公式

(10/8)-1

ダイヤカットの要素分割を増やすプログラム daiya4bun.f daiya4bun.inp daiya4bun.obj 画像

今は拘束、載荷節点がinpファイルに出力するように書き換え中。。

(10-7)16:00

最下段の節点座標がうまくいかず。 クリック クリック

(10/6)16:00

段数は少ないけど、いちおう節点の座標出力はうまくいった。これから段数をふやしていくます。 クリック

(10/2)14:00

inpファイルの出力は出来た。2、4、6の節点座標が間違ったみたい。 クリック

(9/30)14:00

ダイヤカットの分割数を増やしているところ。 それらしくはなってきたが形がおかしい。inpファイルも出力せず クリック

(9/25)正午

折りたたみ4段のプログラムの修正したやつtubure_k42.f tubure_k42.inp

(9/25)

折りたたみ4段のプログラムtubure_k41.f tubure_k41.inp どこか間違えたので続きはまた明日。。

後藤ちゃちゃ(08/9/25)

プログラム名が私のやつと同じだとのちのち混同してしまうので、 工藤の4段だったら、tubure_k4.fとか適宜、名前を変更して下さい。 自分でプログラムを修正していくときも、 過去のバージョンは残しながら修正してください。 でないと、どこを直した途端におかしくなったのか/うまくいくようになったのかが分からなくなってしまう。

修論関係

蛇腹PET

座屈 座屈2

リンク

-Linux関連-
UNIXコマンド一覧

-Tex関連-
texコマンド一覧
texグラフィック系

-CalculiX関連-
ccxコマンド
CalculiXを試す
cgxでモデル作成(マニュアル)
fortran基礎
inpファイルの説明

-OpenOffice_-
オープンオフィスの使い方

-その他-
単位変換

後藤ちゃちゃ07/11/8

ちょっと苦労しましたが、I型梁データ作成用の試作品 ccxkataiza.fを作りました。 cgxで描かせてみる限り、メッシュ切りはうまくいっているのではないかと思います。 拘束条件だの載荷条件のところは知りませんが。 まずは、先端のたわみとか出して、手計算と合うかどうか確認しましょう。 というか、4年生二人にやらせてください。

前期時間割

時限
1.2建設構造工学
3.4構造力学特論材料設計学
5.6土質工学特論
7.8構造実験(TA)マーケティング論