堅固山の卒論日誌 - *1/4モデル

目次

• 1/4モデル
• メモ
• 分け方（その１）←失敗モデル
• 分け方（その２）←成功？
• meshの要素について
• 三角孔について｛孔大きめ｝
• 失敗作
• ↑の作成方法
• 失敗モデル
• geometryでfuseさせたとき
• 成功例
• 三角孔ありモデル(スパン７ｍ)
• スパン７ｍ三角孔あり補剛材あり（全体）
• エラー
• 荷重を線載荷にしたときにたまに出てくる
• Netgen1D2D3Dでメッシュを作成したとき
• 飛び出しについて
• 1/4モデル座屈荷重
• スパン7m孔なし補剛材なしモデル
• スパン7m孔あり補剛材なしモデル(三角孔大きめ)
• スパン7m孔あり補剛材ありモデル
• スパン7m孔あり補剛材なしモデル
• スパン７ｍ、幅員を半分にしたモデル
• 板の座屈
• 板座屈の公式
• FEMと板座屈の比較
• オイラー座屈との比較
• 鋼板のオイラー座屈
• 目標
• 1/28-2/3
• 1/21-1/27
• 1/14-1/20
• 1/4-1/13
• 12/21-12/28
• 12/10-12/20
• 12/3-12/9
• 11/26-12/2
• Trahairの座屈公式(メモ)
• 11/19-11/25
• 11/12-11/18
• 11/5-11/11
• 10/29-11/4
• 10/22-10/28
• 10/15-10/21
• 10/8-10/14
• 9/30-10/7
• 夏休み課題
• 目標：オイラー座屈をsalomeで解析できるようにする。
• 9/26-9/29
• 9/19-9/25
• 9/12-9/18
• 9/5-9/11
• 8/29-9/4
• 8/22-8/28
• 8/15-8/21
• 8/8-8/14
• 8/1-8/7
• salomeで座屈解析
• 参考
• Asterメモ
• データの出力（ParaViSからCalcへ）
• ParaViS
• Calc
• gnuplotメモ
• 卒論日誌
• 課題
• 6/20課題
• 6/13課題
• 6/6課題
• 手計算との比較
• 5/30課題
• 5/23課題
• 四面体要素（面載荷）
• 立方体要素
• 5/16課題
• 条件

1/4モデル

メモ

分け方（その１）←失敗モデル

のように鋼材の座屈させたい部分に対応する木材部分を0.1mmだけ削る。

• 鋼材と木材の２材料のためpartitionが必要。図の左側（寸法2835mmのところ）の木材部分の鋼材と接している面をCreate Groupで選択しておく。
• 木材と鋼材（左側のみ）をFuseする。
• Partitionする。objectsがFuse、Tool objectsにさっきの面を選択（上下の木材２つ）
• 木材、鋼材それぞれをgroupを作っておく。図にもあるように間の面を選択しておく。（仮にaida1,aida2としておく）
• meshに移動したら、Partitionした方のmeshを最初に切る。（今回はNetgen1D2D3Dで切った）
• 切ったmeshを右クリック。Create groups from Geometryを選択し、あらかじめgroupしておいたものを選択してapply（ここでは木材(ki)、鋼材(kou)、間の面(aida1,aida2,yoko)、拘束線(kotei)としている。）
• 木材部分は一体化したいので、meshをいつも通りに切らない。（ここが重要）
• まず、木材（図の右側）を普通に Create meshで同じように切りたい形式（ここではNetgen1D2D3D）にしてmeshを切る。
computeせずにmeshを選択して、Create submeshを使う。submesh画面で、まずgeometryのところにくっつけたい面を選択（aida1 or aida2）。 AlgorithmをUsing Existing 2D Elementsを選択。 すると、Hypothesisの右の歯車が選択できるようになるのでそれを選択する。 Source Facesというのがでてくるので選択。その中のGroup facesをくっつけたいもう一方の面を選択する。このとき注意するのは、geometryで面を選択するのではなく、あらかじめmeshを切ってgroupを入れたものから選択をする（mesh1の中のaida1 or aida2。このときgroups of facesから選択することに注意）。 Applyしてmesh画面に戻ったら、computeする。そうすると、最初に切ったmeshに合わせてmeshを切ることが出来る。
• 鋼板も同様に行い、できたmeshをBuild compoundで選択して、Merge coincident nodes and elementsにチェックを入れて、Toleranceはそのまま？にしてApplyして完成。

commファイルは-ここ

計算結果がおかしいため失敗？ 座屈も分かりづらい。

分け方（その２）←成功？

ｍ単位で作成すると飛び出しが生じてしまう。mm単位で作るとうまくいく。

• 木材、鋼材を座屈させたい部分とそうでない部分とで分割して作成。
• 固定する線を作成し、Partitionをしておく。（後でも大丈夫？）
• 木材、鋼材をそれぞれFuseする。→木材が２つ、鋼材が１つの３つがある状態
• 木材の鋼材と接する面をcreate groupsで選択しておく。
• 木材、鋼材をFuseする。
• 先程選択しておいた面をPartitionで切る。あとはgroupを作っておく。

• meshを切る。3D Automatic Tetrahedralizationで切るとNetgen1D2D3Dで切るよりも要素を少なくできる。
• 2次要素に変換。あとは計算すれば完成。

追記

3D Automatic Tetrahedralizationで計算できたが、見直してみると細かい部分がうまく切れない可能性がある。切れない場合は、Netgen1D2D3Dを使って、submeshを駆使するとうまく切ることができた。

meshの要素について

三角孔について｛孔大きめ｝

図は適当に作ったので見にくいが、とりあえず寸法が高さ147.92mm、底辺143.66mmになるように円の入れる位置を調整した。 また、円に接するように接線を引いた。寸法は図にあるとおり。（手計算で算出）

細かいメモ

図の左下の角を原点としたとき、
-円の中心の座標
→左下(20,20)、左上(20,127.92)、右(123.66,20)

-円の接点の座標
→左下(20,0)(0,20)、左上(0,127.92)(34.43,141.77)、右(123.66,0)(138.09,33.85)



上記の点を組み合わせて三角孔の面を作る。

面の状態の三角孔ができたら立体にする。

面から立体にする方法は Create an extrusionを選択→画面がでてきたらExtrusionのところは一番右を選択→Baseに立体にしたい面を選択→DX,DY,DZに立体の寸法を入れる。（三角孔の場合は厚さ9mmなので厚くしたい方向に9を入れる。）→Apply and Closeで完成。

あとは適当に調整する。

↑のやり方ではmeshがうまく切ることが出来なかった。

手計算のため小数点数桁程度の誤差から、段差が生じてしまった。 そのため、新しいやり方で修正。

まず、元の三角形を作成（面）。寸法は上の図のとおり。 面なのでFillet2Dを用いる。

Operations→Fillet2D→Planar Faceに作った三角形を選択→Vertexesに頂点を選択→Radiusは今回は半径20で切りたいので20にする。→Apply and Close

３つの角に丸みを帯びさせることが出来る。

（ｋ２に三角孔あり鋼板をanaakikoubase2.brepで、三角孔部分をsankakufillet.brepで保存しています。）

失敗作

• 鋼板が座屈しない上になんか飛び出てきた。commファイルかpartitonに問題がある？

↑の作成方法

geometryにて

-鋼板(3.15m,0.009m,0.5m)、木材(3.15m,0.42m,0.12m)を作成。
-鋼板の座屈させたいところに対応する木材部分を0.1mmだけcutする。
-Fuseする。（前段階で鋼板と木材の境界面をgroupで選択しておく）
-二材料なので分けるためにPartitionを入れる。
-境界条件に合わせてgroupを作る。


mesh以降はいつも通りに行った。

• 拘束条件を少し変えたら別のところから飛び出してきた。

失敗モデル

geometryでfuseさせたとき

• geometryで鋼板（3150mm,9mm,500mm）と木材部分（2835mm,420mm,120mm　×２/315mm,419.9mm,120mm　×２）をfuseする。(cutを使わない)
• 鋼板と木材を分けるためにpartitonを入れる。
• meshでNetgen1D2D3Dで切る。←この時点でうまく切れていない感じがする。

計算してpostproで見てみると鋼板は座屈せず、飛び出しだけ起こる。

成功例

• 1次モード

• 2次モード

• 3次モード

• 4次モード

• 5次モード

荷重は大きいままなので、おそらく合っていない。→二次要素にしたら、値が大きく変わった。

間違えていた箇所を訂正↓

斉藤さんの座屈荷重にかなり近づいた！

三角孔ありモデル(スパン７ｍ)

• 1次モード

• 2次モード

• 3次モード

• 4次モード

• 5次モード

上は3D Automatic Tetrahedralizationでの結果。 Netgen1D2D3Dで計算したとき、座屈荷重が変わった。座屈変形も4次モードと5次モードで逆になった。

こっちの方があっている気がする。

meshをさらに細かくしてみた。荷重のかけ方を線載荷に変更。

同じmeshの細かさで、荷重を面載荷に。

荷重を後輪荷重に変更。

スパン７ｍ三角孔あり補剛材あり（全体）

• １次モード

補剛材が先に座屈してしまう。

エラー

荷重を線載荷にしたときにたまに出てくる

Erreur utilisateur dans un calcul élémentaire :                                 
      Le matériau est nécessaire sur la maille : M78108                             
      - option de calcul élémentaire : RIGI_MECA                                    
      - type_element                 : MECA_TETRA4                                  
                                                                                    
    Conseils :                                                                      
      * Peut-etre avez-vous oublié de renseigner le mot clé CHAM_MATER dans la      
    commande courante.                                                              
      * Dans la commande AFFE_MATERIAU, avez-vous affecté un matériau sur la maille 
    incriminée ?                                                                    

→解決。面載荷で使っていたモデルを線載荷に変えた際、Create groupsで作ったグループがmeshで反映されていなかった。このエラーは恐らくグループがcommファイルに書いてあるグループと合っていないときに発生？同じ名前のグループが２つ以上あるときも発生する。

Netgen1D2D3Dでメッシュを作成したとき

Netgen1D2D3Dで作成したはずなのにalgorthmがNetgen2D3Dに変わってしまっている。

飛び出しについて

飛び出しがどこから出ているか見てみたところ、木材から飛び出していた。 試しに応力を確認してみたところ、

木材の削った部分の角に集中していた。

1/4モデル座屈荷重

スパン7m孔なし補剛材なしモデル

飛び出しが出てしまったモードがあった。そのモードは無視している。

スパン7m孔あり補剛材なしモデル(三角孔大きめ)

スパン7m孔あり補剛材ありモデル

スパン7m孔あり補剛材なしモデル

荷重を1/4$P$で計算

スパン７ｍ、幅員を半分にしたモデル

• 1次モード

• 2次モード

• 3次モード

• 4次モード

• 5次モード

板の座屈

板座屈の公式

\[ P=\frac{\pi^2 E t^3}{12b^2\{ 1-\nu^2 \}} \times m \{ \frac{mb}{a}+\frac{a}{mb} \}^2 \]

FEMと板座屈の比較

ヤング率E=2.06e11

境界条件は４辺単純支持

厚さを1mmずつ大きくしていく。

厚さ4mmあたりが限界？

オイラー座屈との比較

• オイラー座屈：\[ P=\frac{\pi^2*EI}{\ell^2} \]
• 直径1mm、高さ50mmの円柱供試体

• 正方形断面（1辺1mm、高さ50mm）

• 長方形断面（縦5mm、横10mm、高さ50mm）

長方形は誤差がやや大きめ？

鋼板のオイラー座屈

境界条件は両端単純支持条件 ヤング率E=2.06(GPa)

目標

1/28-2/3

• ハウトラスの場合を確かめてみる。
[状況]

斜材の座屈を確認したいがモード数を相当多くしないと確認できない。

また、モード数を多くしすぎると座屈変形が見れない。

モード数を50まで増やしたらようやく確認できた。

三角孔は実際のモデルのもので、斜材の座屈荷重は5440kNだった。

木材と鋼材を一体化しはじめたところの斜材が先に座屈した。

孔ありモデルで、三角孔寸法を5%ずつ大きくしたときの座屈荷重が比例に近い状態だったため 回帰直線を作ってみた。

縦軸は座屈荷重(kN)、横軸は三角孔なしの鋼板体積を１としたときの鋼板体積の倍率

孔ありの回帰直線を見ると、

孔なし鋼板の方が5kNほど孔ありモデルよりも強いという結果となった？

回帰直線メモ

y=62.0299x+239.6109

相関係数r=0.999595

x平均=0.8209080243

y平均=258.7284286

1/21-1/27

• 斜材部分の局部座屈について（拡大したときの斜材部分の寸法とか）
[状況]

斜材幅は約4.1228mmずつ、垂直部材幅は約5.2744mmずつ小さくなっていた。

斜材は座屈しないと思われる。（引張方向のため）

垂直部材では座屈が生じた。 幅が63.3533mmで座屈荷重約4000kNで座屈した。

垂直部材でも相当大きな荷重でないと座屈は起きづらいと考えられる。

先週のデータをグラフにしてみた。(縦軸：座屈荷重(kN)、横軸：鋼板体積($m^3$))

グラフを見る限り鋼板体積と座屈荷重は比例関係にあると思われる。 要素数がだいたい同じくらいになるようにしたが微妙に誤差があるため、きれいな直線にはなっていないと考えられる。

線載荷に変えてやり直し。

salomeのGeometryから、Operations→Transformation→Scale Transformで三角孔を拡大していたが、どうやら三角孔の面積の倍率を変えるのではなく寸法の倍率を変える機能だった。

今更ながら上部非一体、下部一体化モデルで解析してみた。

(１次モード)

解析結果は

1/14-1/20

• 先週の継続
• 三角孔を5%ずつ大きくする。
[状況]

先週に引き続き別の孔をあけてみた。

先週はワーレントラス風に孔をあけてみたが斜め部分が細くても座屈が起こりづらいことが分かったので、今回は斜めの部分を多くしてなるべく多めに孔をあけてみた。

モデルは、

計算結果は

のようになった。鋼板体積はワーレントラス風のものよりも約１$m^2$大きくなってしまったが座屈荷重は5kN程強くなった。

三角孔の開け方を変えてみた。 今までは三角孔の位置がちょうど良くなるように孔と孔の間の鋼板部分を調整していたが、三角孔の中心を求めて、その点から拡大するみたいな感じで三角孔を開けた。

三角孔の中心はGeometry画面から、Measures→Center of Massで中心点を作成した。

５％ずつ大きくして１次モードを比較してみる。

三角孔を大きめに作っていたため、三角孔を修正。 同様に５％ずつ大きくして１次モードを比較してみる。

25%のときmeshがうまく切れず、少しだけmeshの細かさを調節して計算できるようにしたが座屈挙動を見ることができなかった。たまに座屈挙動を見れない場合があるが原因はいまだ分からず。

→座屈挙動は前に出てきたエラーだったため改善して見ることができた。（モード数を増やすことで改善）

追記）

別の日に同じやり方でもう一度解き直してみたらなぜかわからないが計算ができた。 そのため、25%のところを修正。

1/4-1/13

• 先週の継続
[状況]

手計算で作ったモデルに欠陥があったため修正。 Filletを使ったほうが断然早いうえ、楽だった。

とりあえず改めて作りなおしたモデルでもう一度計算してみた。

元のモデルと比較すると1kN程低くなっているが問題はないと考えられる。

三角孔を20%拡大したモデルで計算してみたが、postproを見ることができなかった。 ちなみに計算結果は

となり、破壊試験の250kNを下回ってしまった。

今度は15%拡大したモデルで計算してみた。こちらは座屈モードが見れた。

となり、250kNを上回った。meshをもっと細かくすれば（今が要素数９万程度）座屈荷重が数kN低くなると考えられるので、三角孔の大きさは15%拡大したものが限界ではないかと思われる。

試しに別の形で孔をあけてみた。 元の三角孔を2倍にした形で孔をあけてみた。

のようになり、１次座屈モードは元のモデルと変わりはなかった。 計算結果は

座屈荷重は250kNを上回り、鋼板体積は元の三角孔の鋼板よりも約２$m^3$少なくできた。

孔の間の部分は細くなっているが、座屈は約1000kNで起こる計算結果となったため、大丈夫？

12/21-12/28

• 三角孔の大きさを大きくしてみる。
[状況]

とりあえず三角孔を10%ずつ大きくしてみた。

三角孔を30%拡大すると、破壊試験の250kNに耐えられなくなった。 25%拡大してみたところ、

とちょうど破壊試験の荷重と同じくらいの座屈荷重となった。

三角孔は25%の大きさまでは拡張できると考えられる。

三角孔との間の部分も拡大してたので、その部分も詰められればもう少し体積を減らすことが出来るかもしれない。

孔の数は元のモデルと同じにして解析してみた。（25%拡大） 計算はできたものの、座屈モードを見ることが出来なかった。 ちなみにエラーは、

  VERIFICATION A POSTERIORI DES MODES


   !--------------------------------------------------!
   ! <EXCEPTION> <ALGELINE5_23>                       !
   !                                                  !
   ! pour le concept  buckle                          !
   !                                                  !
   !   dans l'intervalle  ( -0.423861  ,  0.583837 )  !
   !   il y a th&#233;oriquement  11 charge(s) critique(s) !
   !   et on en a calcul&#233;  10                         !
   !--------------------------------------------------!

→解決。モードの量を10から12に増やしたら見れるようになった。

250kNを下回ってしまったため、孔の数は変えず20%拡大で解析してみる。

三角孔をいじっていたら、寸法が微妙に違っていることが判明。（1,2mmくらい小さい）

そのため、自分で作りなおしてみた。おおまかな作り方はここ

今後はこっちをベースに作っていくことにする。

12/10-12/20

• 非一体化部分の鋼板上部のオイラー座屈、Trahair座屈を確認してみる。
• 中間発表準備
[状況]

鋼板上部のオイラー座屈 鋼板の非一体化部分(l=630mm,h=120mm,b=9mm)の両端固定条件

1/4モデルの鋼板上部(l=315mm,h=120mm,b=9mm)の片持ち条件

1/4モデルの鋼板上部(l=315mm,h=120mm,b=9mm)の横ねじれ座屈

補剛材ありモデルで解析してみたが、補剛材を全体につけるべきか、実際のモデルにあわせるべきなのかがわからない。全体につけたとき、現在使っているcommファイルの性質なのか、座屈変形が+-があり（+が鋼板下部で-が鋼板上部の座屈）、補剛材をつけると上側の座屈荷重のみがおおきくなるため、５次モードまで見ることができない（小さい順に表示される）。

12/3-12/9

• 孔ありモデルで解析してみる。
[状況]

いろいろ見直してたら、完全に自分のミスがあったことが判明した。 境界条件を間違えていた。

斉藤さんはスパン７ｍで座屈させていたので、今後は合わせていくことにする。

スパン７ｍの孔なし補剛材なしの1/4モデル

meshが粗めなので、もう少し細かくする。

3D Automastic Tetrahedralizationではうまく切れていなかったため、 Netgen1D2D3Dでmeshを切ってみた。

Netgen1D2D3Dを使うと、鋼板と木材の境界部分がやけに細かく切られてしまうため飛び出しが発生してしまっている気がする。

鋼材を座屈させたい部分とそうでない部分とで分けて(Fuse)作成していたが、一枚の鋼板にして計算させてみた。

3倍くらいmeshを細かくしてみたら、

5倍くらいmeshを細かくしてみたら、

11/26-12/2

• 横ねじれ座屈を確認してみる。
[状況]

Trahairの座屈公式とsalomeでの値とを比較してみる。

Trahairの座屈公式(メモ)

$P$: 座屈荷重(N) $E$: ヤング率(N/mm^2) $I_{y}$: Y軸回りの断面二次モーメント $G$: せん断弾性係数 $J$: ねじり定数 $I_{w}$: そりねじり定数

• b=500mm,t=9mm,l=3000mm　で、片持ち条件で横ねじれ座屈をさせてみる。

• l=315mmに変えて上と同条件で横ねじれ座屈させてみる。

↑スパンが短いとTrahairの座屈公式が使えないため、l=315mmは良くない。

11/19-11/25

• オイラー座屈で計算できるモデルで確かめてみる。
[状況]

• オイラー座屈で計算できるよう、座屈させたい鋼板の逆側から荷重をかける。

結果、座屈荷重は53.2979(MN)であった。 また、鋼板のみの場合は3.16431(MN)であった。

→二次要素に変えて計算してみたら、0.65355(MN)になった。

試しに鋼板の座屈させたい部分のみでやってみたら、0.278457(MN)になった。

→二次要素では、0.161215(MN)

線形要素と二次要素でかなりの差があった。もしかしたら、meshを改善すれば座屈荷重を近づけることができるかもしれない。

• 1Dアルゴリズムが入らないことについて

Using Existing 2D Elementsを使っていることが原因と考えたが違った。 試しに同じ大きさの立方体要素でmeshを切ってBuild Compoundさせてみたところ、同様に　'Global 1D algorithm is missing'　とエラーが出てきた。Build Compoundが原因である可能性が高い。

• 荷重がおおきくなってしまうことについて

! <A> <CALCULEL_26>                                                                                    !
   !                                                                                                      !
   ! Les commandes CALC_ELEM et CALC_NO ne doivent pas etre utilis&#233;es en "reuse" lorsque                  !
   ! l'on utilise l'un des mots cl&#233;s suivants :                                                           !
   !   * MODELE                                                                                           !
   !   * CHAM_MATER                                                                                       !
   !   * CARA_ELEM                                                                                        !
   !   * EXCIT                                                                                            !
   !   * GROUP_MA / MAILLE                                                                                !
   !                                                                                                      !
   ! Risques & conseils :                                                                                 !
   !   Si on utilise de tels mots cl&#233;s, la structure de donn&#233;es enrichie manque de coh&#233;rence.             !
   !   Par exemple, le champ de d&#233;placement DEPL a pu etre calcul&#233; avec un champ de mat&#233;riau,             !
   !   alors que le champ d'&#233;nergie cin&#233;tique ECIN_ELEM_DEPL est calcul&#233; avec un autre champ de mat&#233;riau. !
   !                                                                                                      !
   !                                                                                                      !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette                                       !
   ! alarme, vous pouvez obtenir des r&#233;sultats inattendus !                                               !
   !------------------------------------------------------------------------------------------------------!
 

 <A> <CALCULEL_26>！
   ！ ！
   ！ CALC_ELEM CALC_NOの注文とすると、「再利用」するために使用すべきではありません！
   ！ 1は、次のキーワードを使用しています。
   ！ *モデル！
   ！ * CHAM_MATER！
   ！ * CARA_ELEM！
   ！ * EXCIT！
   ！ * GROUP_MA / MESH！
   ！ ！
   ！リスク＆ヒント！
   ！このようなキーワードを使用している場合は、データ構造が矛盾して濃縮されました。 ！
   ！例えば、SHIFT変位フィールドは、材料の分野で計算することができます！
   ！運動エネルギーフィールドECIN_ELEM_DEPLは、別のフィールドの材料を用いて計算されます。 ！
   ！ ！
   ！ ！
   ！これはアラームです。あなたはこの意味を理解していない場合は！
   ！アラームは、あなたが予期しない結果を得ることができます！ ！
   ！------------------------------------------------- -------------------------------------------------- ---！

  

commファイルの中の、

result=CALC_ELEM(reuse =result,
                 MODELE=pmodel,
	　　　　 CHAM_MATER=matprops,
	         RESULTAT=result,
	         OPTION=('SIGM_ELNO','SIEQ_ELNO','SIEF_ELGA',),
		 EXCIT=_F(CHARGE=FLoad,),);

result=CALC_NO(reuse =result,
               RESULTAT=result,
	       OPTION=('SIGM_NOEU','SIEQ_NOEU','FORC_NODA',),);

reuseがエラーのでる原因？

• ここの書いてあるとおりにcommファイルを修正してみたが、計算することができなかった。

現状、1Dが入らないこととcommファイルのミス？の２つが解決できていない。

11/12-11/18

• Netgen1D2D3Dのアルゴリズムについて調べてみる。
• meshのミスを修正する。（境界条件、1Dが入らない状態を直す）
[状況]

• 境界条件を直したら、斉藤さんの座屈挙動と似たような座屈結果が得られた。しかし、荷重がまだ合わない。
Netgen1D2D3Dでやっているが、飛び出しがでないため、四面体でも六面体でも飛び出しをなくすことが可能？
• 六面体要素でやってみたところ、座屈変形はなぜか内側に変形してしまった。

しかし、座屈荷重は少しだけ変化して１次モードが770kNと、斉藤さんの結果の約10倍まで抑えることができた。

• compoundするとなぜか1Dがなくなってしまう。

このページ にNetgenについて

Netgenアルゴリズムは入力パラメータに厳密に従わない。 実際のメッシュは、必要以上に高密度になります。 そこにはいくつかの要因があります：

NETGENは、実際には、エッジの離散化のためのセグメントパラメータの数を使用しません。 このパラメータは、ローカル要素サイズ（指定された点でのサイズ）を定義するためにのみ使用されるため、近い辺のローカルサイズが相互に影響します。 NETGENは、エッジの曲率に応じて要素サイズを追加で制限します。 セグメントの局所的なサイズは、近い三角形のサイズに影響する。 NETGENによって生成される1Dメッシュの要素とそのサイズの順序は、Regular_1Dアルゴリズムで生成される1Dメッシュの要素の順序とそのサイズと異なるため、2Dと3Dのメッシュが異なります。

11/5-11/11

• calculixで解析して比較してみる。
[状況]

 calculixのやり方が分からず苦戦中。

10/29-11/4

• 1/4モデル解析
[状況]

座屈変形はうまく行っているように見えるが座屈荷重に差がありすぎる。 斉藤さんのデータでは78kNだが、自分は14000kNであきらかにおかしい。

10/22-10/28

• inpファイルからsalomeに導入？
[状況] Calculix→salome形式ファイルに直接変換するツールがないらしい。 Abaqus形式ファイルをNastran形式にフリーソフトで変換すれば Nastran形式ファイル→Gmsh→Universalファイルに出力できるみたい。

Abaqus形式ファイル、Nastran形式ファイルがよくわからないので調べ中

 Abaqus 入力(.inp)ファイル
    ノード座標、要素トポロジー、およびセクションの定義などのモデルデータが含まれる。等方性材料を使用している場合は、入力ファイルには材料特性データも含まれるらしい。

inpファイルはAbaqus形式ファイルらしい。

10/15-10/21

• 1/4モデルでいろいろ条件変えてやってみる
[状況]

拘束条件を少し変えてみたが失敗。

鋼板が座屈してくれないので、1/4モデルの鋼板だけでやってみた。

拘束条件がおかしいかもしれない

鋼板も一応座屈を見ることができたので、1/4モデルもできる？

鋼板は二次要素なので、二次要素でやるべき？

1/4モデルいろいろ試して見てるけど未だ成功せず。

10/8-10/14

• 板の座屈の理論値との比較（厚さを変えてみる）
• 斉藤さんのモデルを探す
[状況]

縦長比a/b=2で板座屈とFEMを比較してみた。 オンサイト木橋の1/4モデルで解析してみたが失敗。

9/30-10/7

• ４辺単純支持条件で板を座屈させてみる。
[状況]

１辺100mm、厚さ2mmの正方形板を４辺単純支持で座屈。

誤差はそこまで大きくないため座屈解析ができていると思われる。

meshを線形要素でやるとまったく理論値と合わない。二次要素が適している。 二次要素のnodesが6万か7万を超えると座屈解析ができなくなった。多少粗めのmesh分割にしないと解析できない。（メモリが足りない？みたいなエラーが出た）

画像：1次モード、2次モード

夏休み課題

目標：オイラー座屈をsalomeで解析できるようにする。

9/26-9/29

• 目標：発表準備

• 状況：概要は一応完成？

       meshを２次要素にすると板もオイラー座屈で確認できた。

9/19-9/25

• 目標：境界条件を変えてみる

• 状況：板の座屈をやってみたが、今までエラーがでなかったところにエラーが発生した。一端固定他端自由の条件では座屈ができた。

　　　板の座屈は、変形はうまくいっているように見えるが数値が合わない。

9/12-9/18

• 目標：先週の目標の継続

• 状況：meshを細かくしたところ理論値に近づけることができた。
　　　載荷面の中心に軸を持ってきたらうまくできた。ただ、誤差が5次モードだと大きいので、まだ修正箇所があるのかもしれない。

      両端単純支持でもできるようにしたい。（9/17完了）
      

9/5-9/11

• 目標：先週の目標の継続

• 状況：長さを変えて長方形断面部材を座屈させてみたが、誤差がやや大きいためどこかにミスがありそう。
　　　試しに断面を正方形にしたところ、座屈結果がでなかった。おそらくcommファイルが間違っているのだと思われる。（長方形断面も？）9/8

      正方形断面の問題はやはりcommファイルだった。しかし、理論値との差が大きすぎるため修正が必要。9/10
      自由端部分をx=0,y=0の拘束を加えてみたが、理論値との差は大きいままだった。

8/29-9/4

• 目標：長方形断面部材の座屈解析

• 状況：円柱の座屈でのオイラー座屈との誤差を求められた。
　　　長方形断面では、値が１桁も違うため、修正が必要。

8/22-8/28

• 目標：数値の読み取り方を確認。

• 状況：進展なし

8/15-8/21

• 目標：他の形での座屈解析

• 状況：幅120mm高さ9mmスパン315mmの板を単純支持にしてオイラー座屈を解析
　　　変形がおかしいため、commファイルに修正が必要だと思われる。もしくは変形の見方がよく分かっていないため、自分の見方が間違っているかもしれない。

8/8-8/14

• 目標：commファイルの修正、先週の目標の継続

• 状況：エラー修正（MATER_RIGIではなく、MATR_RIGIだった）
　　　修正後postproにより、座屈変形が見れるようになった。 　　　

8/1-8/7

• 目標：円柱供試体（直径1mm,高さ50mm）をsalomeで座屈解析する。

• 状況：commファイルを作ってみたが、
　　　『Mots cles inconnus :MATER_RIGI』 　　　とエラーが表示された。そのため、修正が必要。

salomeで座屈解析

参考

• http://caelinux.org/wiki/index.php/Contrib:KeesWouters/beambuckling
• http://framabook.org/docs/Code_Aster/beginning_with_Code_Aster_JPAubry_20131206.pdf
• http://www.st.nagasaki-u.ac.jp/ken/matsuda/lecture/steel_st/ohp-7.pdf
• http://www.st.nagasaki-u.ac.jp/ken/matsuda/lecture/steel_st/ohp-8.pdf

Asterメモ

• コマンド
http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?plugin=attach&refer=SALOME-Meca%A4%CE%BB%C8%CD%D1%CB%A1%B2%F2%C0%E2&openfile=00-00.pdf

データの出力（ParaViSからCalcへ）

ParaViS

rmedファイルを開く。その後、見たい面だけ表示させる。（チェック外す）

Select Points Throughを選択、面を囲う。

右側にあるSplit Horizontalを選択。出てきた画面の一番下（Spreadsheet view）を選択。

見たい分の節点数が確認できたら、エクスポート（csvファイル）

Calc

必要な要素以外は削除しておく。

例えばxx方向の応力が見たい場合は、（RESU_SIGM_NOEU:0、Point1(y)、Point2(Z)）だけ残す。

順番を並び替える。'昇順に並び替え'を使うと便利。

（'昇順に並び替え'は'範囲を拡張'を選択）

区切る時は行間をあけておく。

できたら名前をつけて保存。ファイル形式はcsvファイルを選択。'フィルタ設定を編集する'にチェックをいれる。

フィールドの区切りをスペースに。テキストの区切りは削除。

あとは端末からgnuplotできる。

gnuplotメモ

set term png
set output 'gazou.png'
splot 'data' w l

• eog gazou.pngでgazou.pngを確認。
• gftpでstrのgotouhanにアクセス
• FTP→Preference→FTP→Passiveのチェックを外しておく
• 左側がローカル（自分のPC）、右側がサーバー
• サーバー側のpublic_htmlの中がウェブ上に公開されてる
• public_html/j2016/の中に各自のディレクトリをつくる（差し支えなければ、誰のディレクトリかわかるようなディレクトリ名で）
• さっき作ったPC上の画像ファイルを「→」をクリックして、サーバー側の自分のディレクトリにアップロード
• strサーバのpublic_htmlの下は、http://www.str.ce.akita-u.ac.jp/~gotouhan/と対応している。

gnuplot使い方 http://www.str.ce.akita-u.ac.jp/~gotou/linux/vine.html#gnuplot

gnuplotのグラフをTexに取り込む http://www.str.ce.akita-u.ac.jp/~gotou/linux/vine.html#gnutex

卒論日誌

課題

6/20課題

[固定端] （data1:固定端の応力分布、data4:手計算）

[固定端から1要素となり] （data2:固定端から1要素となりの応力分布、data5:手計算）

[梁の中央] （data3:梁の中央の応力分布、data6:手計算）

6/13課題

[固定端] （data8:固定端の応力分布、data5:手計算）

[固定端から1要素となり] （data9:固定端から1要素となりの応力分布、data6:手計算）

[梁の中央] （data10:梁の中央の応力分布、data7:手計算）

6/6課題

• 応力を縦軸(z)にして書き換える
• 断面(z=0)を3次元プロットに一緒にプロットする
• 直方体要素で、固定端から1要素となりの断面の応力分布は、どれくらい変わるか
• 梁の真ん中の断面の応力分布はどうか（以上の3箇所を比較）
• 手計算$\sigma=\frac{M}{I}y$とどれくらい合うか。手計算で求まる面を3次元プロットにあわせてプロット。

[固定端の応力分布] （data2:固定端の応力分布、data0:基準面）

[固定端から1要素となりの応力分布] （data3:固定端から1要素となりの応力分布、data0:基準面）

[梁の中央の応力分布] （data4:梁の中央の応力分布、data0:基準面）

[３箇所の応力分布の比較] （data2:固定端の応力分布、data3:固定端から1要素となりの応力分布、data4:梁の中央の応力分布)

手計算との比較

[固定端] （data2:固定端の応力分布、data5:手計算）

[固定端の一つとなり] （data3:固定端から1要素となりの応力分布、data6:手計算）

[梁の中央] （data4:梁の中央の応力分布、data7:手計算）

5/30課題

5/23課題

四面体要素（面載荷）

[2次要素]

立方体要素

[線形要素]

[2次要素]

5/16課題

Timoshenko梁: v=Pl^3/3EI + Pl/kGA

条件