こんなのとかどうでしょう。
$\ell$が短いと局部座屈
$\ell\がある程度長いと全体座屈
ただし、sin1半端のグラフ（拘束側の回転が拘束されたような感じ）
座屈荷重は$\ellがある程度短いとオイラー（最小軸）よりも小さい。ある程度長いとオイラー（最小軸）とオイラー（中立軸$I_{y}$）の間に入る。ただし、間に入るけど、FEMの座屈荷重はオイラー（中立軸）に近い値。理由として考えられるのは、sin1半端に近いモードになることで、座屈荷重が大きくなる。座屈荷重が$\ell$が短いと小さいのは、全体座屈にならず局部座屈になるからであると考えられる。
| $\ell(m)$ | FEM(MN) | $オイラー$I_{y}$(MN) | オイラー$I_{v}$(MN) |
| 0.5 | 1384168 | 14231968 | 5855439 |
| 1 | 1122413 | 3557992 | 1463860 |
| 1.5 | 781544 | 1581330 | 650604 |
| 2 | 527298 | 889498 | 365965 |
| 2.5 | 370993 | 569279 | 234218 |
| 3 | 271934 | 395352 | 162651 |
| 4 | 162459 | 217566 | 91491 |
| 5 | 1047961 | 142320 | 58554 |
| 6 | 76354 | 98833 | 40662 |
| 7 | 57079 | 72162 | 29875 |
| 8 | 44389 | 55594 | 22872 |
| $\ell(m)$ | FEM(MN) | オイラー$I_{y}$(MN) | オイラー$I_{v}$(MN) |
| 0.5 | 1224567 | 6188872 | 3338744 |
| 1 | 754693 | 1547218 | 839686 |
| 1.5 | 446368 | 687652 | 373194 |
| 2 | 283395 | 386804 | 209922 |
| 2.5 | 192920 | 247554 | 134350 |
| 3 | 138721 | 171913 | 93298 |
| 4 | 81112 | 96701 | 52480 |
| 5 | 53017 | 61889 | 33587 |
| 6 | 37333 | 42978 | 23324 |
| 7 | 27717 | 31576 | 17136 |
| 8 | 21404 | 24175 | 13120 |
| L(m) | FEM(MN) | オイラー$I_{y}(MN) | オイラーI_{v}(MN) |
| 0.5 | 1320274 | 1423168 | 2171392 |
| 1 | 116584 | 355792 | 542848 |
| 2 | 905833 | 889498 | 135712 |
| 3 | 298460 | 395332 | 60316 |
| 4 | 169321 | 222375 | 33928 |
| 5 | 110431 | 142320 | 21714 |
| 6 | 77946 | 98833 | 15079 |
| 7 | 58047 | 72612 | 11079 |
| 8 | 44976 | 55594 | 8482 |
| $\ell$(m) | FEM | オイラー$I_{y}(MN)$ | オイラー$I_{v}(MN)$ |
| 0.5 | 1221219 | 6188872 | 3338744 |
| 1 | 851733 | 1547218 | 839686 |
| 2 | 496754 | 386804 | 209922 |
| 3 | 258737 | 171913 | 93298 |
| 4 | 86675 | 96701 | 52480 |
| 5 | 56973 | 61889 | 33587 |
| 6 | 40570 | 42978 | 23324 |
| 7 | 30508 | 31576 | 17136 |
| 8 | 23877 | 24175 | 13120 |
| L(m) | ccx座屈荷重(MN) | オイラーの理論値(x軸の断面二次)(MN) | 相対誤差(％) | 弱軸理論 |
| 1 | 1804685 | 355792 | 49.3 | 542848 |
| 2 | 904774 | 889498 | 1.7 | 135712 |
| 3 | 385280 | 395332 | 2.5 | 60316 |
| 4 | 217566 | 222375 | 2.2 | 33928 |
| 5 | 139465 | 142320 | 2.1 | 21714 |
| 6 | 96948 | 98833 | 1.9 | 15079 |
| 7 | 71276 | 72612 | 1.8 | 11079 |
| 8 | 54597 | 55594 | 1.8 | 8482 |
| L(m) | ccx座屈荷重(MN) | オイラーの理論値(x軸の断面二次)(MN) | 相対誤差(％) |
| 1 | 1476140 | 1309341 | 12.7 |
| 2 | 329270 | 327335 | 0.06 |
| 3 | 147350 | 145482 | 1.3 |
| 4 | 83102 | 81834 | 1.5 |
| 5 | 53225 | 52374 | 1.6 |
| 6 | 37010 | 36371 | 1.8 |
| 7 | 27203 | 26721 | 1.8 |
| 8 | 20833 | 20458 | 1.8 |
| L(m) | ccx座屈荷重(MN) | オイラーの理論値(x軸の断面二次)(MN) | 相対誤差(％) |
| 1 | 16374 | 16265 | 0.7 |
| 2 | 4113 | 4066 | 0.2 |
| 3 | 1830 | 1847 | -0.9 |
| 4 | 1030 | 1021 | 0.9 |
| 5 | 659 | 651 | 1.2 |
| 6 | 459 | 455 | 0.9 |
| 7 | 339 | 332 | 2.1 |
| 8 | 258 | 254 | 1.6 |
| L(m) | ccx座屈荷重(MN) | オイラー弱軸理論値(MN) | オイラーx軸理論値(MN) |
| 1 | 1543201 | 839686 | 3232690 |
| 2 | 485291 | 209922 | 808172 |
| 3 | 221807 | 93298 | 359187 |
| 4 | 125880 | 52480 | 202043 |
| 5 | 80895 | 33587 | 129308 |
| 6 | 56302 | 23324 | 89796 |
| 7 | 41421 | 17136 | 65973 |
| 8 | 31741 | 13120 | 50111 |
| L(m) | ccx座屈荷重(MN) | オイラー弱軸理論値(MN) | オイラーx軸理論値(MN) |
| 1 | 970933 | 1634643 | 9860722 |
| 2 | 682367 | 408661 | 2465181 |
| 3 | 361126 | 181627 | 1095635 |
| 4 | 215404 | 102165 | 616295 |
| 5 | 141431 | 65386 | 394429 |
| 6 | 99475 | 45407 | 273909 |
| 7 | 73623 | 33360 | 201239 |
| 8 | 56637 | 25541 | 154073 |
| L(m) | ccx座屈荷重(MN) | オイラー弱軸理論値(MN) | オイラーx軸理論値(MN) |
| 1 | 4019464 | 3273353 | 15899143 |
| 2 | 1728118 | 818338 | 3974785 |
| 3 | 827713 | 363706 | 1766571 |
| 4 | 476523 | 204585 | 993696 |
| 5 | 308130 | 130934 | 635966 |
| 6 | 215168 | 90926 | 441643 |
| 7 | 158609 | 66803 | 324472 |
| 8 | 121699 | 51146 | 248424 |
| L(m) | ccx座屈荷重(MN) | オイラー弱軸理論値(MN) | オイラーx軸理論値(MN) |
| 1 | 1611995 | 693300 | 2216120 |
| 2 | 405152 | 370450 | 554030 |
| 3 | 183168 | 77033 | 246235 |
| 4 | 103588 | 43331 | 138507 |
| 5 | 66470 | 27732 | 88644 |
| 6 | 46228 | 19258 | 61558 |
| 7 | 33933 | 14148 | 45226 |
| 8 | 26041 | 10832 | 34626 |
| L(m) | ccx座屈荷重(MN) | オイラー弱軸理論値(MN) | オイラーx軸理論値(MN) |
| 1 | 1481800 | 693300 | 2216210 |
| 2 | 409820 | 370450 | 554053 |
| 3 | 183090 | 77033 | 246245 |
| L(m) | ccx座屈荷重(MN) | オイラーの理論値(x軸の断面二次)(MN) | 相対誤差(％) |
| 1 | 16350 | 16265 | 0.52 |
| 2 | 4112 | 4066 | 1.1 |
| 3 | 1807 | 1847 | 2.1 |
| 4 | 1030 | 1021 | 0.7 |
| 5 | 659 | 651 | 1.2 |
shellsection OFFSET=0.5にすれば穴は空かない
http://www.str.ce.akita-u.ac.jp/~gotouhan/hatanaka/tyuukan.odp
ls /usr/local/bin/ccx*を実行して、パソコン内に古いバージョンが何々インストールされてるか調べてみてください。 他の人のパソコン(去年モード図を描いてた人の)に入ってるかも。
まあ、inpファイル内の座標にdatファイルの座屈モードを足して、それをfrdの形式に整形して出力する こともできるから、古いバージョンが利用できなくてもなんとかなるとは思います。 古いバージョンが利用できそうにない場合は、中間発表は直方体要素のモード図だけ見せて、 シェル要素のモード図は中間発表後に対処しましょう。
| x方向分割数 | y方向分割数 | ccxででた値 |
| 2 | 10 | 2186 |
| 2 | 30 | 2178 |
| 2 | 60 | 2189 |
| 2 | 80 | 2178 |
| 2 | 100 | 2164 |
| 2 | 1000 | 2164 |
| 4 | 4 | 2309 |
| x方向分割数 | z方向分割数 | ccxで出た座屈応力(MN) |
| 2 | 2 | 629.3 |
| 2 | 10 | 552.3 |
| 2 | 100 | 543.9 |
| 2 | 150 | 552.7 |
| 2 | 200 | 557.0 |
| 2 | 1000 | 538.5 |
| 分割数 | ccxで計算した座屈荷重(MN) |
| 1 | 計算しない |
| 2 | 7028 |
| 3 | 2813 |
| 4 | 1522 |
| 5 | 957 |
| 6 | 659 |
| 7 | 481 |
| 8 | 368 |
| 9 | 291 |
| 10 | 235 |
| 11 | 195 |
| 20 | 60 |
| 30 | 28 |
| 40 | 17 |
Pcr=2196MNであり、分割数でかなりの値の変動があった。(ちなみにtyoku.fでなくccxkataiza.fで分割させてもn=10が もっとも近くて同じような感じになった。)
モードは100
10MNでなければ計算してくれなかった。
バグ取りするときは、 まずは、なるべく簡単なモデル例えば、 7型の2辺をそれぞれ1要素、2辺の接続部分の角を1要素の計3要素16節点とかで 要素分割して、要素番号と節点番号の関係を手書きしたものと比べてみるとか して下さい。って、ことをほんとはやってもらった方がいいんだけど、 たぶん、以下の間違いでは?
n3=n1+(nfx+1)*(nfy+1)
n4=n3+1
n5=n1+(nfx+1)*(nfy+1)
直方体の手前左下節点n3と奥左上節点n5が同じ番号になっているのはおかしいです。
n3=n1+(nfx+1)
n4=n3+1
n5=n1+(nfx+1)*(nfy+1)
上記のI型のオリジナルのままに直すと7型になると思います。
1,-1.000000 -1.000000 -1.000000 2,-1.000000 1.000000 -1.000000 3,-1.000000 -1.000000 1.000000 4,-1.000000 1.000000 1.000000 5,0.000000 1.000000 -1.000000 6,-0.000000 -1.000000 -1.000000 7,0.000000 1.000000 1.000000 8,-0.000001 -1.000000 1.000000 上の情報は右側のobject 9,-0.000001 -1.000000 1.000000 10,0.000000 1.000000 1.000000 11,-0.000000 -1.000000 -1.000000 12,0.000000 1.000000 -1.000000 13,0.999999 -1.000001 1.000000 14,1.000000 0.999999 1.000000 15,1.000000 -1.000000 -1.000000 16,1.000000 1.000000 -1.000000 下の情報は左側のobjectこれは節点を共有(9,10,11,12)してますが別々のobjectを表記しています。
blenderのsplitという機能で節点共有でわけることができました
6-5 z / /| / 2-1 8 →x | |/ ↓ 3-4 y
20m×20m×0.1mの板のN分割したときにz軸方向に荷重1kNの等分布荷重をかけその際のZ方向の変位を分割数ごとで算出し、理論値と比較しました。
blenderをinpファイルにした。
| 日付 | 曜日 | 開始〜終了 | 作業時間(h) | 立合 | 作業内容 |
| 4/15 | 水曜日 | 10:00-12:00 | 2 | 後藤 | ログインの仕方 |
| 4/27 | 月曜日 | 12:00-15:00 | 3 | 後藤 | ブラインドタッチ、コマンドタッチの練習 |
| 5/1 | 金曜日 | 10:00-15:00 | 5 | 後藤 | 外国語文献購読翻訳、ブラインドタッチの練習 |
| 5/8 | 金曜日 | 10:00-16:00 | 6 | 後藤 | ブラインドタッチ、コマンドタッチの練習、翻訳 |
| 5/11 | 月曜日 | 12:00-16:00 | 4 | 後藤 | ブラインドタッチ、翻訳 |
| 5/12 | 火曜日 | 12:00-14:00 | 2 | 後藤 | ブラインドタッチ |
| 5/14 | 木曜日 | 10:00-14:00 | 2 | 後藤 | 課題1、2 |
| 5/18 | 月曜日 | 10:00-14:00 | 4 | 後藤 | 課題1、2、翻訳 |
| 5/19 | 火曜日 | 13:00-15:00 | 2 | 後藤 | 課題1、2 |
| 5/21 | 木曜日 | 12:00-16:00 | 4 | 後藤 | 課題3 |
| 5/22 | 金曜日 | 12:00-17:00 | 5 | 後藤 | 課題3 |
| 5/25 | 月曜日 | 13:00-16:00 | 3 | 後藤 | 課題3、翻訳 |
| 5/26 | 火曜日 | 13:00-15:00 | 2 | 後藤 | 課題3、翻訳 |
| 5/29 | 金曜日 | 12:00-16:00 | 4 | 後藤 | 課題4、翻訳 |
| 6/2 | 火曜日 | 10:00-16:00 | 6 | 後藤 | 課題4 |
| 6/5 | 金曜日 | 10:00-14:00 | 2 | 後藤 | 課題5 |
| 6/8 | 月曜日 | 10:00-15:00 | 5 | 後藤 | 課題5、翻訳 |
| 6/9 | 火曜日 | 11:00-15:00 | 4 | 後藤 | 課題5 |
| 6/12 | 金曜日 | 13:00-15:00 | 2 | 後藤 | 課題6 |
| 6/15 | 月曜日 | 12:00-16:00 | 4 | 後藤 | 課題6 |
| 6/16 | 火曜日 | 12:00-14:00 | 2 | 後藤 | 任意に分割できるプログラム作成 |
| 6/18 | 木曜日 | 11:00-13:30 | 2、5 | 後藤 | gnuplotの使い方について。 |
| 6/24 | 水曜日 | 14:30-18:00 | 3.5 | 後藤 | blenderの使い方 |
| 6/29 | 月曜日 | 11:00-17:00 | 6 | 後藤 | blenderの使い方をwikiに載せる。 |
| 7/1 | 水曜日 | 15:00-17:00 | 2 | 後藤 | blenderで作成した図を任意に分割する。 |
| 7/6 | 月曜日 | 10:00-16:00 | 6 | 後藤 | blenderをinp化 |
| 7/7 | 火曜日 | 11:00-14:00 | 3 | 後藤 | blenderファイルをinp化 |
| 7/9 | 木曜日 | 11:00-16:00 | 5 | 後藤 | 文献翻訳 |
| 7/13 | 月曜日 | 12:00-15:30 | 3.5 | 後藤 | blenderで作ったinpファイルの節点の法則性 |
| 7/30 | 木曜日 | 10:00-14:00 | 4 | 後藤 | blenderで切抜きを覚える |
| 8/3 | 月曜日 | 10:00-15:00 | 5 | 後藤 | blenderで切抜き |
| 8/10 | 月曜日 | 12:00-15:00 | 3 | 後藤 | 直方体をblenderで四面体分割する |
| 8/11 | 火曜日 | 12:00-16:00 | 3 | 後藤 | 直方体をblenderで四面体分割 |
| 8/17 | 月曜日 | 13:00-15:00 | 2 | 後藤 | 直方体を立体要素で分割 |
| 9/3 | 木曜日 | 12:00-15:00 | 3 | 後藤 | ccxで四面体分割を表示させる。 |
| 9/8 | 火曜日 | 14:30-16:00 | 2.5 | 後藤 | ccxで四面体要素で分割させる |
| 9/15 | 火曜日 | 8:30-18:00 | 8.5 | 後藤 | 木橋見学、木材加工所見学 |
| 9/16 | 水曜日 | 8:30-15:00 | 6.5 | 後藤 | 木橋をblenderで作成する |
| 9/18 | 金曜日 | 9:00-13:00 | 4 | 後藤 | blenderでモデル作成 |
| 9/24 | 木曜日 | 10:00-16:30 | 6.5 | 後藤 | Latexで夏休みの課題をまとめる |
| 9/25 | 金曜日 | 10:00-16:30 | 5.5 | 後藤 | Latexで夏休みの課題をまとめる |
| 9/29 | 火曜日 | 10:00-16:30 | 6.5 | 後藤 | ccxでC3D4で作ったモデルに荷重をかける |
| 9/30 | 水曜日 | 11:00-16:00 | 5 | 後藤 | ccxで四面体要素に応力をかける。 |
| 10/5 | 月曜日 | 10:00-16:00 | 6 | 後藤 | 夏休みの課題をまとめる |
| 10/6 | 火曜日 | 10:00-16:00 | 6 | 後藤 | 夏休みのまとめ |
| 10/9 | 金曜日 | 12:00-15:00 | 3 | 後藤 | tetgen調査 |
| 10/13 | 火曜日 | 9:00-16:00 | 7 | 後藤 | blenderで節点を入力して任意の図形をかく |
| 10/15 | 木曜日 | 12:00-13:00 | 1 | 後藤 | blenderで任意の図形を描く練習 |
| 10/19 | 月曜日 | 10:00-14:00 | 4 | 後藤 | blenderで節点共有を調べる |
| 10/20 | 火曜日 | 10:00-15:00 | 5 | 後藤 | blenderで違う材料諸元を用いたモデルを書く |
| 10/21 | 水曜日 | 13:00-16:00 | 3 | 後藤 | blenderで任意モデルの作成 |
| 10/22 | 木曜日 | 10:00-16:00 | 6 | 後藤 | 寸法計算 |
| 10/23 | 金曜日 | 10:30-16:00 | 5.5 | 後藤 | sagcadの練習 |
| 10/26 | 月曜日 | 10:30-16:00 | 5.5 | 後藤 | sagcadの練習 |
| 10/27 | 火曜日 | 9:00-16:00 | 7 | 後藤 | 設計の勉強 |
| 11/6 | 金曜日 | 8:50-16:00 | 7 | 後藤 | 部材力の基本を学習 |
| 11/9 | 月曜日 | 9:00-18:00 | 9 | 後藤 | 各々の部材力の影響線の算出 |
| 11/10 | 火曜日 | 10:30-17:00 | 6.5 | 後藤 | 各々の部材力の影響線の算出 |
| 11/16 | 月曜日 | 10:00-17:00 | 7 | 後藤 | 積分で断面2次をだすプログラムの作成 |
| 11/17 | 火曜日 | 11:00-15:00 | 4 | 後藤 | 山形綱の断面二次を出すプログラムの作成 |
| 11/18 | 水曜日 | 12:00-16:00 | 4 | 後藤 | 山形綱の断面二次半径を出すプログラムの作成 |
| 11/19 | 木曜日 | 12:00-16:00 | 4 | 後藤 | tetgenで分割 |
| 11/20 | 金曜日 | 9:00-17:00 | 8 | 後藤 | tetgenで山形綱を分割 |
| 11/24 | 火曜日 | 10:00-17:00 | 7 | 後藤 | 上に同じ |
| 11/25 | 水曜日 | 13:00-17:00 | 7 | 後藤 | L字綱のプログラム作成 |
| 11/30 | 月曜日 | 9:30-17:00 | 7.5 | 後藤 | L字綱に荷重をかける |
| 12/1 | 火曜日 | 12:00-16:00 | 4 | 後藤 | 影響線で荷重の算出 |
| 12/3 | 木曜日 | 10:30-16:00 | 6.5 | 後藤 | L字綱の境界条件の設定 |
| 12/4 | 金曜日 | 10:30-16:00 | 6.5 | 後藤 | L字綱に荷重をかけ、座屈、静的たわみを見る |
| 12/7 | 月曜日 | 10:30-16:00 | 6.5 | 後藤 | 手計算で算出したΔLと解析の比較 |
| 12/8 | 火曜日 | 9:00-17:00 | 8 | 後藤 | 静的変位を見る |
| 12/9 | 水曜日 | 12:00-16:00 | 4 | 後藤 | フックの法則との比較 |
| 12/10 | 木曜日 | 10:30-16:00 | 5.5 | 後藤 | 座屈モード |
| 12/11 | 金曜日 | 10:30-16:00 | 5.5 | 後藤 | 座屈モード |
| 12/16 | 水曜日 | 12:00-1:00 | 13 | 後藤 | shell座屈 |
| 12/17 | 木曜日 | 10:30-16:30 | 6 | 後藤 | shell座屈 |
| 12/18 | 金曜日 | 12:00-22:00 | 10 | 後藤 | 中間発表のまとめ |
| 12/19 | 土曜日 | 20:00-21:00 | 1 | 後藤 | 中間発表のグラフつくり |
| 12/20 | 日曜日 | 13:00-16:00 | 3 | 後藤 | 中間発表の練習 |
| 12/22 | 火曜日 | 12:00-17:00 | 5 | 後藤 | 山形鋼の節点番号をだすプログラムづくり |
| 12/25 | 金曜日 | 12:00-19:00 | 7 | 後藤 | 山形綱に荷重をかける。 |
| 12/27 | 日曜日 | 17:30-20:00 | 2.5 | 後藤 | 解析 |
| 1/6 | 水曜日 | 9:00-17:00 | 8 | 後藤 | 解析 |
| 1/7 | 木曜日 | 10:00-17:30 | 7.5 | 後藤 | 解析 |
| 1/8 | 金曜日 | 9:00-0:00 | 15 | 後藤 | 解析 |
| 1/9 | 土曜日 | 11:30-14:00 | 3 | 後藤 | 座屈モードの図をかく |
| 1/11 | 月曜日 | 11:00-16:00 | 5 | 後藤 | 支部原稿 |
| 1/12 | 火曜日 | 11:00-20:00 | 9 | 後藤 | 支部原稿 |
| 1/13 | 水曜日 | 11:00-17:00 | 6 | 後藤 | 支部原稿 |
| 1/19 | 火曜日 | 10:00-21:30 | 11.5 | 後藤 | 解析 |
| 1/20 | 水曜日 | 10:00-21:00 | 11. | 後藤 | 解析 |
| 1/21 | 木曜日 | 11:00-21:00 | 10 | 後藤 | 支部原稿 |
| 1/22 | 金曜日 | 11:00-24:00 | 13 | 後藤 | 支部原稿 |
| 1/25 | 月曜日 | 14:00-17:00 | 3 | 後藤 | 支部原稿と発表用スライド |
| 1/26 | 火曜日 | 14:00-17:00 | 3 | 後藤 | 発表用スライド |
| 1/27 | 水曜日 | 14:00-17:00 | 3 | 後藤 | 発表用スライド |
| 1/28 | 木曜日 | 14:00-17:00 | 3 | 後藤 | 発表用スライド |
| 1/29 | 金曜日 | 13:00-16:00 | 3 | 後藤 | 発表用スライド |
| 1/30 | 土曜日 | 13:00-17:00 | 4 | 後藤 | 卒論レジュメ作成 |
| 2/1 | 月曜日 | 13:00-18:00 | 5 | 後藤 | 卒論レジュメ作成 |
| 2/2 | 火曜日 | 10:00-19:00 | 9 | 後藤 | 発表練習、スライド訂正 |
| 2/3 | 水曜日 | 11:00-20:00 | 9 | 後藤 | 発表練習、スライド訂正 |
