| 日付 | 時間帯 | 作業時間 | 内容 | 立合 |
| 4/7 | 10:00〜11:00 | 1h | 卒論 | |
| 5/21 | 14:00〜16:30 | 2.5h | ゼミ課題（サロメ） | |
| 5/23 | 16:00〜18:30 | 2.5h | ゼミ課題（サロメ） | |
| 5/26 | 14:00〜17:00 | 3h | ゼミ課題（サロメ） | |
| 5月合計 | 8h | |||
| 6/1 | 16:00〜18:30 | 2.5h | ゼミ課題（サロメ） | |
| 6/2 | 13:00〜18:00 | 5h | ゼミ課題（サロメ） | |
| 6/3 | 15:00〜18:00 | 3h | ゼミ課題（サロメ） | |
| 6/7 | 15:00〜18:00 | 3h | ゼミ課題 | |
| 6/8 | 13:00〜20:00 | 5.5h | ゼミ課題 | |
| 6/13 | 13:30〜17:00 | 3.5h | ゼミ課題 | |
| 6/15 | 16:30〜20:00 | 3.5h | ゼミ課題 | |
| 6/18 | 14:30〜16:00 | 1.5h | ゼミ課題 | |
| 6/20 | 17:00〜18:00 | 1h | ゼミ課題 | |
| 6/21 | 16:00〜18:00 | 2h | ゼミ課題 | |
| 6/27 | 13:30〜16:00 | 2.5h | ゼミ課題、外論 | |
| 6/28 | 13:00〜16:00 | 3h | ゼミ課題 | |
| 6/29 | 13:00〜17:00 | 4h | ゼミ課題 | |
| 6月合計 | 40h(48h) | |||
| 7/4 | 15:00〜18:00 | 3h | ゼミ課題 | |
| 7/5 | 13:00〜15:30 | 2.5h | ゼミ課題 | |
| 7/10 | 13:00〜17:00 | 4h | ゼミ課題 | |
| 7/11 | 15:00〜18:00 | 3h | ゼミ課題 | |
| 7/12 | 13:00〜19:30 | 6.5h | ゼミ課題 | |
| 7/13 | 13:00〜16:30 | 3.5h | ゼミ課題 | |
| 7/19 | 13:00〜17:00 | 4h | ゼミ課題 | |
| 7/20 | 13:30〜15:30 | 2h | ゼミ課題 | |
| 7月合計 | 28.5h(76.5h) | |||
| 8/2 | 13:00〜17:30 | 4.5h | 中間 | |
| 8/17 | 14:30-17:30 | 3h | 中間 | |
| 8/20 | 15:30-18:00 | 2.5h | 中間 | |
| 8/23 | 14:00〜19:30 | 5.5h | 中間（エラーなおらない） | |
| 8/28 | 14:00〜19:30 | 5.5h | 中間 | |
| 8/31 | 15:00〜18:30 | 3.3h | 中間（エラー直った） | |
| 8月合計 | 24.5h(101) | |||
| 9/4 | 14:00〜20:30 | 6.5h | 中間 | |
| 9/6 | 13:00〜18:00 | 5h | 中間と実験手伝い | |
| 9/7 | 13:00〜16:00 | 3h | 中間 | |
| 9/12 | 14:30〜18:30 | 4h | 中間 | |
| 9/13 | 14:00〜18:00 | 4h | 中間 | |
| 9/14 | 13:00〜17:30 | 4.5h | vi復習 | |
| 9/15 | 13:00〜19:00 | 4.5h | データ編集 | |
| 9/16 | 14:30〜20:30 | 4h | 概要書く | |
| 9/19 | 15:00〜21:00 | 6h | 概要書く | |
| 9/20 | 11:00〜19:00 | 8h | 概要、スライド | |
| 9/21 | 15:00〜18:30 | 3.3h | 概要、スライド | |
| 9/26 | 14:00〜15:30 | 1.5h | 論文読む | |
| 9月合計 | 54.5h(155.5h) | |||
| 10/10 | 13:00~16:30 | 3.5h | 接触解析練習 | |
| 10/11 | 15:00~19:00 | 4h | 接触解析練習 | |
| 10/12 | 13:30~18:00 | 4.5ｈ | 接触解析練習と論文調べる | |
| 10/16 | 13:00~18:30 | 5.5h | 接触解析 | |
| 10/18 | 13:00~18:30 | 5.5h | 接触解析（エラー出まくる）と摩擦解析練習 | |
| 10/19 | 13:00~16:00 | 3h | 接触解析と摩擦練習 | |
| 10/23 | 13:00~17:30 | 4.5h | 接触解析 | |
| 10/24 | 13:30~17:30 | 4h | 摩擦解析練習 | |
| 10/25 | 14:30~18:30 | 4h | 摩擦解析練習 | |
| 10/26 | 14:00~19:30 | 5.5h | 摩擦解析練習 | |
| 10/31 | 13:30~17:00 | 3.3h | バネについて検証 | |
| 10月合計 | 47.5h(203h) | |||
| 11/1 | 14:30~18:30 | 4h | バネについて | |
| 11/2 | 13:00~17:00 | 4h | バネについて | |
| 11/7 | 15:00~18:00 | 3h | バネについて | |
| 11/8 | 13:00~17:00 | 4h | バネについて | |
| 11/9 | 13:00~16:30 | 3.5h | バネ | |
| 11/13 | 13:00~17:30 | 4.5h | 摩擦 | |
| 11/14 | 13:00~16:30 | 3.5h | 摩擦 | |
| 11/15 | 13:30~18:30 | 5h | 摩擦 | |
| 11/16 | 11:30~15:30 | 4h | 摩擦 | |
| 11/17 | 16:00~18:30 | 2.5h | 二面摩擦 | |
| 11/20 | 13:00~18:00 | 5h | 二面摩擦 | |
| 11/21 | 13:00~15:00 | 2h | 二面摩擦（収束しない） | |
| 11/22 | 13:00~18:00 | 5h | 二面摩擦 | |
| 11/23 | 17:00~21:30 | 4.5h | 二面摩擦 | |
| 11/27 | 17:00~18:00 | 1h | 一面摩擦(計算速度遅い) | |
| 11/28 | 13:00~17:00 | 4h | まだ計算終わらない | |
| 11/29 | 14:30~ | 3h | １面摩擦と２面摩擦 | |
| 11/30 | 13:00~20:30 | 7.5h | 1面摩擦と2面摩擦 | |
| 11月合計 | 86.5h(289.5h) | |||
| 12/5 | 14:00~ | 4h | salome | |
| 12/6 | 13:00~ | 4.5h | salome（収束しない＆結果安定せず） | |
| 12/7 | 13:00~ | 8h | salome | |
| 12/12 | 13:00~16:30 | 3.5h | salome(ボルト付き) | |
| 12/13 | 13:00~16:00 | 3h | salome | |
| 12/14 | 13:00~17:00 | 4h | salome | |
| 12/16 | 22:30~25:00 | 2.5h | salome | |
| 12/19 | 14:30~16:30 | 2h | salome | |
| 12/21 | 14:00~ | 4h | salomeと概要 | |
| 12/24 | 15:00~ | 5h | 概要 | |
| 12/25 | 16:30~24:00 | 7.5h | スライド | |
| 12月合計 | 48h(337.5h) | |||
| 1/6 | 17:00~18:00 | 1h | salome | |
| 1/9 | 13:00~17:00 | 4h | salome(ミーゼス応力とか) | |
| 1/10 | 11:30~14:00 | 2.5h | salome | |
| 1/11 | 17:00~20:00 | 3h | salome | |
| 1/12 | 14:00~15:30 | 1.5h | salome | |
| 1/15 | 15:00~18:00 | 3h | salome | |
| 1/16 | 13:00~17:00 | 4h | salome | |
| 1/17 | 15:00~18:30 | 3.5h | 支部概要 | |
| 1/18 | 15:00~ | 3h | 支部概要 | |
| 1/19 | 12:00~ | 3h | 支部概要 | |
| 1/22 | 13:00~18:30 | 5.5h | 支部概要 | |
| 1/23 | 13:30~18:30 | 5h | 支部概要 | |
| 1/24 | 13:00~ | 2h | 支部概要 | |
| 1/25 | 14:30~ | 2h | スライド | |
| 1/29 | 13:30~ | 3h | salome | |
| 1/30 | 13:00~ | 4h | salomeとスライド | |
| 1/31 | 16:00~18:30 | 2.5h | スライド | |
| 1月合計 | 51.5h(389h) | |||
| 2/1 | 14:00~19:00 | 5h | スライド | |
| 2/2 | 13:00~ | 2.5h | スライド | |
| 2/5 | 13:00~ | 4.5h | スライド・salome | |
| 2/6 | 12:00~ | 4h | 発表練習 | |
| 2/7 | 13:00~ | 3.5h | スライド | |
| 2/8 | 14:30~19:00 | 4.5h | スライド | |
| 2/9 | 13:00~16:30 | 3.5h | スライド | |
| 2/13 | 14:00~18:30 | 4.5h | スライド | |
| 2/14 | 14:00~22:00 | 5h | スライド | |
| 2/15 | 13:00~19:00 | 5h | スライド | |
| 2月合計 | 42h(431h) |
10mm×10mm断面の角材（ヤング率：6GPa, ポアソン比：0.3）の単純梁のたわみを求め、初等梁ティモシェンコ梁に対する相対誤差を求めた。 荷重は面荷重で100N. スパン長さは、100mmを分担。 理論値は初等梁＝0.4167mm、ティモシェンコ梁＝0.4297mm
<一次要素>
| length | 8 | 4 | 2 | 1 |
| たわみ(mm) | 0.3549 | 0.3549 | 0.4177 | 0.4177 |
| 相対誤差(理論値：初等梁) | 14.83% | 14.83% | 1.20% | 1.20% |
| 相対誤差(理論値：ティモシェンコ梁) | 17.41% | 17.41% | 4.19% | 4.19% |
<二次要素>
| length | 8 | 4 | 2 | 1 |
| たわみ(mm) | 0.4338 | 0.4338 | 0.4338 | 0.4338 |
| 相対誤差(理論値：初等梁) | 4.10% | 4.10% | 4.10% | 4.10% |
| 相対誤差(理論値：ティモシェンコ梁) | 0.95% | 0.95% | 0.95% | 0.95% |
<一次要素>
| length | 8 | 4 | 2 | 1 |
| たわみ(mm) | 0.4327 | 0.4699 | 0.5002 | 0.5361 |
| 相対誤差(理論値：初等梁) | 3.84% | 12.77% | 20.04% | 28.65% |
| 相対誤差(理論値：ティモシェンコ梁) | 12.00% | 4.43% | 1.735 | 9.03% |
<二次要素>
| length | 8 | 4 | 2 | 1 |
| たわみ(mm) | 0.5378 | 0.5485 | 0.5555 | 0.5711 |
| 相対誤差(理論値：初等梁) | 29.06% | 31.63% | 33.31% | 37.05% |
| 相対誤差(理論値：ティモシェンコ梁) | 9.38% | 1.16% | 1.30% | 1.61% |
サンドイッチ梁（縦に1mm幅の鋼材を入れた梁）を作ってたわみを見てみた。 木材のヤング率等は前回の課題と同じ。鋼材はE=206GPa、ν=0.3、G=79.2。 また、レンクスは１、メッシュは2次要素で切った。 たわみの理論値は0.09398mmとなった。
| たわみ(mm) | 相対誤差 | |
| 木材を等方性とした時 | 0.1119 | 19.07% |
| 木材を異方性とした時 | 0.1708 | 81.74% |
<木材を等方性としたとき>（E=6GPa、ν=0.4）
<木材を異方性としたとき>
誤差が大きくなってしまったのでせん断係数も考慮した式を使えばいいのかなと思いました。
海老さんから卒論のデータをもらい、同じようにサロメで結果を出せるかという課題だった。
同じように結果を出すことができた。Eficasのいじり方がよく分からないので今後勉強していこうと思った。
前回の課題のデータを使い、サロメと手計算それぞれで求めたひずみと応力の比較をした。 サロメでひずみが出せなかったので、応力からひずみを出した。
| 手計算 | サロメ | |
| 応力(N/m^2) | 3.6393 10^6(木) | 5.8131 10^6 |
| 97.438 10^6(鋼) | ||
| ひずみ | 4.73 10^-4 | 7.56 10^-4(木) |
| 0.282 10-4(鋼) |
6/22のサンドイッチ梁において、サロメと手計算それぞれで求めたひずみと応力の比較をした。サロメで見たのはＸ方向。
| 手計算 | サロメ | |
| ひずみ | 5.639 10^-4 | 5.2550 10^-4 |
| 応力(N/mm^2) | 116.1634(鋼)/3.3834(木) | 115.5825(鋼)/3.0491(木) |
前回の課題でのひずみの分布の仕方のグラフ。5.2〜5.8mmは鋼材部分の値でそれ以外は木材での値。鋼材のほうが若干ひずみが大きいことが分かる。
それから、サロメで求めたひずみにヤング率をかけて、サロメ・手計算で求めた応力と一致するかを確認した。
前回と同じサンドイッチ梁の木材部分において、この行列 を使って応力を求め、サロメで求めた応力と一致するか試した。
| 行列によって求めた応力 | サロメで求めた応力 | |
| 0mm | 2.8686 | 2.8688 |
| 2.5mm | 3.0521 | 3.0661 |
| 8mm | 3.0970 | 3.0967 |
小数点以下2桁あたりまで一致した。
近藤さんからもらったデータを使い、プレストレス木箱桁橋の荷重とひずみの関係を調べる。ひずみはサロメから求めた。荷重は対傾構ありなしそれぞれについて自重と積雪について計算した。以下の表のようになった。
| 荷重(MN/mm^2) | ひずみ(10^-4) | |
| 対傾構あり積雪あり | 0.0136759 | 4.630721 |
| 対傾構あり積雪なし | 0.0076566 | 2.592559 |
| 対傾構なし積雪あり | 0.0135336 | 4.589838 |
| 対傾構なし積雪なし | 0.0075144 | 2.548455 |
対傾構あり積雪ありの時のひずみの様子。
対傾構あり積雪なしの時のひずみの様子。
対傾構なし積雪ありの時のひずみの様子。
対傾構なし積雪なしの時のひずみの様子。
これらのひずみと荷重の関係を表にしてみた。
対傾構ありのほうが対傾構なしに比べてややひずみが大きかった
のエラーがなおらなくて半月以上たってしまったので、試しに適当に中央載荷でやってみたら （固定しすぎていたら3方向固定から2方向固定にするorバネは弱くないかorメッシュの切り方を変える、などでこのエラーは解決できた（2/7））
みたいになった。ネジ穴周りは
穴周りで色が変わるかと思ったけど、この場合は一色であった。エラーが解決できるようにがんばります。
やっとエラーが直った。メッシュの切り方を変えたら上手く行った。ネジ穴がない場合、4個の場合、6個の場合の軸力方向の応力を比べた。
ネジ穴2列と3列のときの添接板のネジ穴周りの応力を比べた。なお応力はほぼ対称の値だったので半円ぶんだけグラフにした。
穴3列のほうが応力が大きかった。$\sigma$=$E$$\epsilon$ より穴6つの方がひずみも大きくなると考えられる。
次回までに
今週できたことなど
次回までに
今週できたことなど
次回までに
今週出来たことなど
菊池さんのデータを応用して柔らかいゴムではなく弱いバネを使ったモデルを作った。このモデルで摩擦力は正しく発生しているか確認中である。
次回までに
今週出来たことなど
| バネ(x/y/z) | 変位(引っ張った方向) | 最大のひずみ | 最大の応力 |
| 0/0/0 | 1.439 | 15.6471 | 490 |
| 10/10/10 | 1.200 | 15.6331 | 500 |
| 10/0/0 | 1.1998 | 15.6334 | 490 |
| 50/50/50 | 1.199 | 15.6382 | 490 |
| 50/0/0 | 1.199 | 15.6394 | 490 |
| 100/100/100 | 1.200 | 15.6446 | 490 |
| 100/0/0 | 1.199 | 15.6471 | 490 |
今週できたことなど
| バネ定数(N/mm) | 静止摩擦係数 | 動き出す時の荷重(N)(手計算) | 動き出した時の荷重(N)(salome) | 変位(mm)(手計算) | 変位(mm)(salome) |
| 10 | 0.4 | 40 | 0.5 | 1 | 4.8077 |
| 100 | 0.4 | 40 | 0.5 | 0.1 | 0.4981 |
| 10 | 1 | 動かない | 0.5 | 動かない | 4.8077 |
| バネ無し | 0.4 | 40 | 0.5 | ？ | 125 |
次回までに
11/13 菊地さんのデータの時間を1000分割してみが、やはり最初から滑りだすし変位の仕方は線形である。摩擦の入れ方を確認したが間違いなどは見つけられなかった。解析の仕方を変えてみたが上手く収束しない。 11/15 ついに摩擦を取り入れた計算が上手く行った。マスター面とスレーブ面を逆にしていた。計算時間は時間を100分割で1時間くらい。 salomeの本P136に大きい面がマスターと書いてあった、たしかにそのとおりであった。 菊地さんのK2のデータのスレーブ面とマスター面を逆にすると、菊地さんの卒論日誌にあるような400Nになったときす滑りだす摩擦解析ができる。
| 静止摩擦係数 | 動き出す時の荷重(N)(手計算) | 動き出した時の荷重(N)(salome) | 変位(mm)(salome) |
| 0.4 | 40 | 40 | 149.494 |
| 0.6 | 60 | 60 | 99.4478 |
次回までに
今週出来たことなど
滑り出し方をグラフにしてみた。左のグラフが滑りだす瞬間までを拡大したグラフ。手計算では引っ張る力が165kN×2×0.4=132kN=185N/mm^2の時の滑りだすと出た。
次回までに
今週できたことなど
次回までに
今週できたことなど
次回までに
今週できたことなど
次回までにと今後やること
今週できたことなど
次回までに
今週できたことなど
ボルトを締めただけで鋼材が降伏しそうだ。同じような条件で高力ボルト継手の解析をしているこの資料（http://www.nilim.go.jp/lab/bcg/siryou/tnn/tnn0827pdf/ks082708.pdf）によると軸力はリラクセーションを考慮して7％くらい小さい軸力を与えている。
次回までに
今週できたことなど
「プレストレス木箱桁橋のボルト接合部の数値モデル化」
DEBUT();
kou=DEFI_MATERIAU(ELAS=_F(E=206000.0, #ヤング率（単位：MPa）
NU=0.3,),); #ポアソン比
MIAL=LIRE_MAILLAGE(FORMAT='MED',);
MIAL=MODI_MAILLAGE(reuse =MIAL,
MAILLAGE=MIAL,
ORIE_PEAU_3D=_F(GROUP_MA=('men2','men1','men3','men4',),),); ＃「DEFI_CONTACT」で使うような接触面はここに入れておく（やんなくてもあんまり問題はない気がする）
newmesh=CREA_MAILLAGE(MAILLAGE=MIAL, ＃バネについて
CREA_POI1=_F(NOM_GROUP_MA='bane',
GROUP_NO=('bane1','bane2','bane3',),),);#モデルを安定させるためのバネをここで定義。「GROUP_NO=」だとcommファイル上でグループ化みたいなことができる（？）
MODE=AFFE_MODELE(MAILLAGE=newmesh, ＃「MIAL」も選択できるけど新しく定義をした「newmesh」にする
AFFE=(_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='3D',),
_F(GROUP_MA='bane', ＃CREA_MAILLAGEで作った安定のためのバネのグループ名
PHENOMENE='MECANIQUE',
MODELISATION='DIS_T',),),);
MATE=AFFE_MATERIAU(MAILLAGE=newmesh, ＃材料の種類。材料の数だけAFFE=_Fが増える
AFFE=_F(TOUT='OUI',
MATER=kou,),);
sktdn=AFFE_CARA_ELEM(MODELE=MODE,
DISCRET=_F(CARA='K_T_D_N',
GROUP_MA='bane',
VALE=(0.1,0.1,0.1,),),); ＃バネ定数。VALE=(x方向、y方向、ｚ方向)
contact=DEFI_CONTACT(MODELE=MODE,
FORMULATION='DISCRETE',
FROTTEMENT='COULOMB', ＃摩擦係数をいれるならこれを入れる
ITER_GEOM_MAXI=100000, ＃
ZONE=(
_F(GROUP_MA_MAIT='men2', ＃マスター面
GROUP_MA_ESCL='men1', ＃スレーブ面 この2面を逆にすると正しい結果が得られない。
ALGO_CONT='PENALISATION', ＃ペナルティ方が収束しやすかった
E_N=206000.0, ＃ヤング率
COULOMB=0.4, ＃摩擦係数
ALGO_FROT='PENALISATION',
E_T=20600.0,), ＃ヤング率の10分の１の値？
_F(GROUP_MA_MAIT='men3',
GROUP_MA_ESCL='men4',
ALGO_CONT='PENALISATION',
E_N=206000.0,
COULOMB=0.4,
ALGO_FROT='PENALISATION',
E_T=20600.0,),),);
kotei=AFFE_CHAR_MECA(MODELE=MODE, ＃固定面。＃載荷、固定は「AFFE_CHAR_MECA」に入力する
DDL_IMPO=_F(GROUP_MA=('kotei1','kotei3',),
DX=0.0, ＃２方向固定にしてローラーヒンジのようにした
DZ=0.0,),);
osae=AFFE_CHAR_MECA(MODELE=MODE,
FORCE_FACE=(_F(GROUP_MA='joumen',
FY=20.24025,), ＃面載荷。165kN/載荷面の面積(mm^2)。単位はMPa
_F(GROUP_MA='kotei',
FY=-20.24025,),),);
tikara=AFFE_CHAR_MECA(MODELE=MODE,
FORCE_FACE=_F(GROUP_MA='saika', ＃真ん中の部材を引っ張る力
FZ=200.0,),);
tikara_f=DEFI_FONCTION( ＃折れ線グラフのイメージ
NOM_PARA='INST',
VALE=(0.0 ,0.0 , ＃VALE=(ｘ座標、ｙ座標、
10.0 ,0.0 ,
20.0 ,1.0 ,),); ＃10~20ｓで真ん中の部材を徐々に引っ張るようにした
osae_f=DEFI_FONCTION(
NOM_PARA='INST',
VALE=(0.0 ,0.0 ,
10.0 ,1.0 ,
20.0 ,1.0 ,),); ＃0~10ｓで徐々に添接板に荷重をかけ、10ｓ以降は与えた力を一定に保つようにした
kankaku=DEFI_LIST_REEL(DEBUT=0.0, ＃時間の分割数の設定
INTERVALLE=(_F(JUSQU_A=10.0,
NOMBRE=1,), ＃0~10ｓは1分割で力を与えた
_F(JUSQU_A=20.0,
NOMBRE=20,),),); ＃10~20ｓは20分割で力を与えた
resu=STAT_NON_LINE(MODELE=MODE,
CHAM_MATER=MATE,
CARA_ELEM=sktdn,
EXCIT=(_F(CHARGE=kotei,), ＃はじめからずっと固定するのでFONC_MULTはいらない
_F(CHARGE=osae,
FONC_MULT=osae_f,),
_F(CHARGE=tikara,
FONC_MULT=tikara_f,),),
CONTACT=contact,
COMPORTEMENT=_F(RELATION='ELAS',),
INCREMENT=_F(LIST_INST=kankaku,),
NEWTON=_F(REAC_ITER=1,),
CONVERGENCE=_F(RESI_GLOB_RELA=0.01, ＃messファイルのこの値の判定基準
ITER_GLOB_MAXI=300,), ＃収束計算を最大300回繰り返す。それで収束しなかったらそこで計算終了する
SOLVEUR=_F(SYME='OUI',),);
resu=CALC_CHAMP(reuse =resu,
MODELE=MODE,
CHAM_MATER=MATE,
RESULTAT=resu,
CONTRAINTE=('SIGM_ELNO','SIGM_NOEU',), ＃SIGMは応力
DEFORMATION=('EPSI_ELNO','EPSI_NOEU',), ＃EPSIはひずみ
CRITERES=('SIEQ_ELNO','SIEQ_NOEU',),); ＃SIEQは使わなかったので分からない
IMPR_RESU(FORMAT='MED',
UNITE=80,
RESU=_F(MAILLAGE=newmesh,
RESULTAT=resu,
NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','EPSI_NOEU','DEPL',),),); ＃DEPLは変位
FIN();
#CHECKSUM:b2f4a2740b714f58c842db0021057ff8 -:FIN CHECKSUM
・
・
・
contact=DEFI_CONTACT(MODELE=MODE,
FORMULATION='CONTINUE',
FROTTEMENT='COULOMB',
ALGO_RESO_GEOM='POINT_FIXE',
REAC_GEOM='AUTOMATIQUE',
ITER_GEOM_MAXI=100,
RESI_GEOM=0.1,
ALGO_RESO_CONT='POINT_FIXE',
ITER_CONT_MAXI=100,
ALGO_RESO_FROT='NEWTON', #「ALGO_RESO_GEOM=POINT_FICXE」,「ALGO_RESO_CONT=POINT_FICXE」,「ALGO_RESO_FROT='NEWTON」の組み合わせが収束しやすかった。
RESI_FROT=0.5,
ZONE=(
_F(GROUP_MA_MAIT='men1',
GROUP_MA_ESCL='boruto_m',
CONTACT_INIT='INTERPENETRE',
ALGO_CONT='PENALISATION',
COEF_PENA_CONT=206005.0,
COULOMB=0.4,
SEUIL_INIT=0.1, #これを入れると収束しやすくなる。結果が少し変わることがあるのであまり大きくしないほうが良いかも
ALGO_FROT='PENALISATION',
COEF_PENA_FROT=20600.0,),
