<< Calculix 10本ノック: 9本目: exampleの実行 Linear/Mesh1 | TOP | Calculix 20本ノック:11本目: exampleの実行 Contact/Hertz2D >>

20200302

Calculix 10本ノック: 10本目: exampleの実行 Kasten

Calculixとは

フリーの解析ソフト。

  • ccx　ソルバー

  • cgx　プリポスト


自分の環境

  • WSL

  • Xming

  • Calculix 2.11

  • imagimagick


データの入手

公式サイトのリンクから行けるGithubのページCalculixのexample。

下記サイトのリンクから、使用しているCalculixのバージョンに合ったファイルを探す。

https://github.com/mkraska/CalculiX-Examples/releases



ノック10本目: exampleの実行 Kasten

Kasten

Kastenは、ドイツ語で"箱"。

大変形、塑性、不整のある座屈（小さい内圧）、シェル要素とreadmeの書き出しにあった。Buckling with imperfectionは不整のある座屈と翻訳するらしい。

箱が座屈するモデル。 S8 (8-節点 四角形シェル要素) 


 cd ~/CalculiX-Examples-2.11/CalculiX-Examples-2.11/Kasten

「~」は、サンプルを解凍したアドレス。


exampleには存在しないけれど

1/4モデルを作成して、

  1. 材料非線形解析なしの場合

  1. 非線形解析なしの場合。

を用意した。


1/4モデルの作成

1/4モデル作成のため、pre.fbdを下のように書き直した。

 pnt p1 0 0 0
 pnt p2 50 0 0
 pnt p3 0 50 0
 line l1 p1 p2 12
 line l2 p1 p3 12
 ​
 seta lower all
 swep all upper tra 0 0 150 20
 flip A002
 ​
 seta miryz l L004
 seta mirxz l L005
 ​
 plot la all
 plus pa all
 plus sa all
 rot y
 rot d 10
 rot r 10
 frame
 hcpy png
 #sys mv hcpy_1.png geo.png
 ​
 div all auto 3
 elty all qu8
 mesh all
 send all abq
 ​
 send all abq pres -0.01
 send lower abq nam
 send upper abq nam
 send mirxz abq spc 246
 send miryz abq spc 156

4bunno1.png

1/4モデルの為に、Kasten.inpの*boundary部分を下記のように書き直した。

 *boundary
 Nlower,1,6
 Nupper,1,6
 *include,input=mirxz_246.bou
 *include,input=miryz_156.bou

1/4モデルの為に、post.fbdをミラーコピーするように下記のように書き直した。

 copy all base1 mir x
 move base1 tra 100 0 0
 copy all base2 mir y
 move base2 tra 0 100 0

1/4モデルの解析結果

1.png


非線形解析なしの場合

非線形がない場合を見るために、Kasten.inpのStepのnlgeomを削った。

2.png


材料非線形がない場合

材料非線形がない場合を見るために、Kasten.inpの*plasticと、その下3行を削った。

あと、上記だけだと収束しなくなったので、*controlを追加して、収束条件を甘くしてしまった。

 *CONTROLS,PARAMETERS=FIELD
 5.0e-1,1e-1,1e-2,,2e-8,1.e-5,1.e-3,1.e-8

3.png


10本完成してしまった。次は25本くらいを目標に挙げておく。
やはり、example頼りになりそう。
【Calculixの最新記事】
posted by yuchan at 07:00 | Comment(1) | Calculix