Calculix 10本ノック: 8本目:Contact/Shell1
Calculixとは
フリーの解析ソフト。
- ccx ソルバー
- cgx プリポスト
環境
WSL
Xming
Calculix 2.11
gmsh
imagimagick
pip2
matplotlib
データの入手
公式サイトからのリンクで行けるGithubにあるCalculixのexample。
下記サイトのリンクから、自分のCalculixのバージョンにあったファイルを探す。
https://github.com/mkraska/CalculiX-Examples/releases
ディレクトリの移動。
cd ~/CalculiX-Examples-2.11/CalculiX-Examples-2.11/Contact/Shell1「~」は、サンプルを解凍した場所のアドレス。
内容
下のような8節点シェルで作ったモデルの接触解析。
ハニカム（ハチの巣、ハニーコム）構造のサンドイッチの最小構成みたいな形の上に半球。
半球に対して、z方向に-2の強制変位をかけている。
初期設定でうまく流れなかったが、contact pairのmaster,slaveを入れ替えたらうまくいった。
Too many cutbacks
初期.inpのまま解析実行したら、下記エラーで止まった。
too slow convergence; the increment size is decreased to 3.125000e-02
the increment is reattempted*ERROR: too many cutbacks
best solution and residuals are in the frd file
cutbacksの回数を増やしたら、止まらないんじゃないかと思い、cutbacksを増やす方法を検索。
controlをいじるといいらしい。 controlには引数が多すぎる。cutbackの回数以外は、初期設定にしておいた。
もしかすると、省略したら（何も書かずにコンマで区切ったら）よかったのかも。
下記*controlを*step(stepの外でも内でもいい？)に追加したら、このエラーは避けられたが…
** cutback 15
*CONTROLS, PARAMETERS=TIME INCREMENTATION
4,8,9,16,10,4,,15,,
0.25,0.5,0.75,0.85,,,1.5,結局別のエラーが出た。
下記は結局出てきたエラー。
too slow convergence; the increment size is decreased to 7.629395e-06
the increment is reattempted*ERROR: increment size smaller than minimum
best solution and residuals are in the frd file
これを解消する気にはならなかった。*controlは消して、別の方法を探った。
contact pairがおかしかった
下記の2行目で、slave,masterを決めていた。
*contact pair, interaction=tool, type=surface to surface
Stopneg,Sinddependent (slave), independent(master)の順番。
これを、下記のように逆さにしたら収束した。
*contact pair, interaction=tool, type=surface to surface
Sind, Stopnegちなみにsindは下記。
stopnegは下記。
結果
表示しているのはds 2 e 23で出せる、最悪主応力(worst Principal Stress)
READMEにある結果と最大値最小値一緒だし、同じ結果でしょ。多分。
