Salome-Meca ビーム要素
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
Salome-Mecaのビーム要素（梁要素）に関するページを作成しま...
これから梁要素を用いた解析を扱うよって人は見ていただける...
また、何か分かったこととかあればこちらに書き込んで共有し...
*参考資料 [#mb640bf5]
[[code-asterによる解析 step-by-step:https://www.opencae.o...
[[結果の出力絡みのお話（CALC_CHAMP,IMPR_RESU）:https://ww...
[[構造要素の概要:https://www.opencae.or.jp/wp-content/upl...
[[code-asterの説明書:https://code-aster.org/doc/v15/en/in...
[[Beginning with code_aster:https://archives.framabook.or...
[[code-aster質問ページ:https://forum.code-aster.org/publi...
*ビーム要素(梁要素)とは [#z866693c]
ビーム要素（梁要素）とは点と点を線で結ぶ要素である.ソリッ...
例えば,橋梁モデルみたいな巨大なモデルで時刻歴応答解析をす...
そのため、モデルは全て梁要素、もしくは桁と塔をシェル要素...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
*単位系の注意点 [#k7f047ce]
Salome-Mecaは単位系について自分で設定しないといけません。~
モデル長mm想定なら...ヤング率(MPa),荷重(N),応力(MPa)~
m想定なら...ヤング率(Pa),荷重(N),応力(Pa)~
*気になること [#q683524a]
-非線形絡みの問題~
非線形には材料非線形、幾何学的非線形、境界非線形の3種類存...
今は弾性域内（線形で片付く範囲）でしか解析してない。一応...
→幾何学的非線形はcable要素でDeformation:GROT_GDEPを用いる...
-橋梁モデルを作れるのか~
簡易モデル（[[千代岡の卒論日誌参照:https://www.str.ce.aki...
→作成可.ケーブルに材料非線形を考慮しないなら、ケーブルをA...
-座屈解析はできるのか~
ビーム要素で座屈解析できるのだろうか？~
もし、試すなら、4年生の春課題で行なったモデルをビーム要素...
-固定条件について~
DDL_IMPOにBLOCAGEというコマンドがあり,これは選択した項目,...
しかし,例えば完全固定条件の1.DX=DY=DZ=DRX=DRY=DRZ=0と2.BL...
両者の違いは何だろうか？
-載荷条件について~
AFFE_CHAR_MECAで載荷条件の順番によって得られる結果が変わ...
[[事例:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/...
→AFFE_CHAR_MECAで使う載荷条件をそれぞれ1つずつ分けると改...
saika1 = AFFE_CHAR_MECA(FORCE_NODALE=_F(FY=-1175000.0,
GROUP_NO=('saika'...
MODELE=allmodel)
saika2 = AFFE_CHAR_MECA(FORCE_POUTRE=(_F(FY=-35250.0,
GROUP_MA=('keta'...
_F(FY=-122000.0,
GROUP_MA=('keta'...
_F(FY=-80000.0,
GROUP_MA=('CW', ...
_F(FY=-60000.0,
GROUP_MA=('syuto...
MODELE=allmodel)
saika4 = AFFE_CHAR_MECA(FORCE_NODALE=_F(FY=-9600000.0,
GROUP_NO=('CW_ten...
MODELE=allmodel)
追記：調査の結果、同じAFFE_CHAR_MRCAに2つ以上の載荷条件を...
[[参考資料:https://forum.code-aster.org/public/d/15289-re...
解決策を見つけたので,[[こちらを参考:https://www.str.ce.ak...
**ばねを用いたピン接合について [#fb9a4942]
Salomeで梁要素を用いてピン接合の再現をするにはばね要素を...
しかし、現状として解析が回るモデルと回らないモデルが混在...
（写真のような,横ばねが存在しているモデルとエラーが出る可...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↓メッセージ内容~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
現在、解決策を模索中であるが、検証をしたものをここにまと...
・バネ要素の両端（部材との共有節点）に自由度の設定をする。~
上の写真の場合、ばねの右端の共有節点DRXの自由度エラーが出...
**材料非線形に関する問題 [#h9737f8e]
***CREA_RESUの設定 [#vcda802b]
材料非線形を考慮する場合,CREA_RESU→COMPORTEMENTの設定をし...
[[説明書:https://code-aster.org/doc/v14/fr/man_u/u4/u4.44...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/CREA_R...
*モデルの作成・片持ち梁を例にした解析手法（編集中） [#p87...
100mm✕10mm✕10mmモデルを対象に説明します.（これまで先輩方...
モデルを片持ち梁として,変位を見てみましょう。
**Geometry [#ud86b1ec]
(1)ポイント(0,0,0)(100,0,0)を作成する.
(2)線（ポイントの右隣にあるアイコンをクリック）を作成する...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)(POU_D_EMを適用する場合のみ)~
ビーム断面を作成する.材料非線形を考慮した解析や複雑な断面...
これはAsterstudyでビームモデルに断面を適用するのに必要な...
長方形断面を作成する場合は画面上部にある長方形マーク（長...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
円形断面を作成する場合は長方形マークの右隣の楕円マーク（...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
方向はモデル軸方向をxにとっているため,OXYにする.そうしな...
（構造研スタイルと同じようにモデル軸方向をzにした場合,方...
モデルとしては,こんな感じになるはず.（作成したのは正方形...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**Mesh [#ed8a38b1]
(1)ビームモデル（line）のメッシュを切る.1D→Wire Discretis...
セグメントの数は分割数を示しているので,任意の数値を入れる.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(2)(POU_D_EMを適用する場合のみ)~
断面のメッシュは2D→NETGEN1D-2D,歯車マークからNETGEN2DPara...
最大・最小サイズの数値を一致させ,任意の数値を入れる.ただ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)メッシュを計算した後,グループを作成する.~
要素の型は,基本的にノード：点,Edge:線,Face:面の3種類を使...
手動で選択する場合,コンテンツにチェックを入れた後,グルー...
全てを選択するは,断面をグループ化するときに使う.
今回は,固定点と載荷点,線要素（ビーム）全て,（断面作成して...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**Asterstudy [#o2e2746d]
(数字-1)の箇所は断面を作成した場合(POU_D_EMを用いる場合)...
(1)LIRE_MAILAGEを選び,計算したメッシュ(ビームモデル)を選...
(2)AFFE_MODELEを選び,MAILAGEはビームモデルのものを選択す...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
梁要素で用いるAFFE_MODELEは[[こちら:https://www.str.ce.ak...
(2-1)DEFI_GEOM_FIBERを選択し,SECTION（断面定義）を選ぶ.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
TOUT_SECTorGROUP_MAで作成した断面を選ぶ。~
GROUP_FIBERは好きな名前をつけてよい。（ここでは,sectionと...
MAILAGE_SECTはLIRE_MAILAGEで適用した断面を選択。~
COOR_AXE_POUTREは0で。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)AFFE_CARA_ELEMを選択.ビーム要素において,ここでの設定は...
POUTREを選び,GROUP_MAをビームモデル,SECTIONをRECTANGLE（...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3-1)MULTIFIBERを選択.ここで,(2)で設定した断面定義をモデ...
(4)DEFI_MATERIAU→ELASで材料特性選択.鋼材を想定しているた...
(5)AFFE_MATERIAU→MODELEにチェック入れ,AFFE選択→(4)で作成...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(5-1){}+を選択→DEFI_COMPORを選択.選択すると,下画面のよう...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEFI_C...
MULTIFIBERを1itemにし,GROUP_FIBERを(2-1)で記入した名前に...
RELATIONの説明は...後でまとめます.(基本的にELAS(線形弾性)...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEFI_C...
(5-2)
(6)AFFE_CHAR_MECAを選択.固定条件と載荷条件は別々に分けて...
まず,固定条件はDDL_IMPOを選択し,GROUP_NO（グループ化[Grou...
DX=DY=DZ=DRX=DRY=DRZ=0と入れてもよいが,ENCASETR[完全固定]...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
載荷条件はFORCE_NODALEを選択.（ここで入れる値は[N],他のFO...
GROUP_NO（グループ化[Group]した点[Node]）に載荷点を入れ,...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(7)AnalysisはMECA_STATICSを選択（POU_D_EMには対応していな...
CHAM_MATERとCARA_ELEMにチェック入れ,EXCITにAFFE_CHAR_MECA...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(7-1)AnalysisはSTAT_NON_LINEを選択(POU_D_Eには対応してい...
(8)IMPR_RESUで.medファイルを作成,FORMATはMEDで！~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
RESUではRESULTATにAnalysisで設定したもの（MECA_STATICS or...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
なお,今回は変位のみを見ることが目的のためCALC_CHAMPは使わ...
**Paravis [#ed37e359]
Asterstudyで解析したmedファイルを開く.適用する前にGenerat...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
その後,Warp by Vectorを選択して変位図を見て,spread sheet ...
*つりあい解析（二段階解析）のやり方 [#gac275b5]
つりあい解析とは、事前に解析した構造解析の結果（1段階目の...
横軸時間をとる時刻歴解析やモデルを先に変形させてから力を...
大きなメリットとして,解析時間の短縮が挙げられる.そして,段...
・使用モデル~
卒論で扱ったケーブル1本の簡易モデルからケーブル(C1)をなく...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・解析条件~
荷重は死荷重とケーブルから得られた張力を掛けている.~
時刻歴応答解析 解析時間は115s,100sまで荷重等をゆっくり載...
109sまでを1段階目,109sから2段階目の解析を行う。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・設定法~
写真のようにStageを2種類用意する。~
chainanalysisが1段階目,chainanalysis-2が2段階目用のStage...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
1段階目 commファイルは[[こちら:http://www.str.ce.akita-u....
基本的に自分が行っている解析のaster-codeをいつもどおり組...
時刻歴解析の場合,区切る時間（今回は109s）までの解析をここ...
2段階目 commファイルは[[こちら:http://www.str.ce.akita-u....
AFFE_MODELEやDEFI_MATERIAUなど1段階目の解析と共通するもの...
Stage1で設定したものがそのまま使えます。
使用するaster-code~
絶対：Analysisに含まれるもの ,Output IMPR_RESU~
人によって変わる：Functions and Listsに含まれるもの （時...
ここでは時刻歴応答解析を行っているため,以下のように設定し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
1段階目の解析終了地点から解析を再開しないと,1段階目の結果...
そのため,解析時間は109sからスタートする。
Analysisでは,それぞれのaster-codeに入っているコマンド：ET...
（静的解析を行う場合,STAT_NON_LINEで。MECA_STATICSにはそ...
また、MODELEやCARA_ELEMといった選択式で入力するところはSt...
RESULTATに1段階目のAnalysisで用いたものを入れる。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
EVOL_NOLIに1段階目のAnalysisで用いたものを入れる。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
OUTPUTは1段階目と異なる.medの名前をつける。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
RESUにはNOM_CHAMの設定とTOUT_OREDEとTOUT_CMPをOUIにする。
・解析の回し方
写真のようにstage同士をつなげることでつりあい解析が可能と...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
初めて回す際、一気に2段階目まで解析を行おうとすると,2段階...
まず,1段階目の解析を回したいので,2段階目の解析に赤い四角...
解析が回ったら、1段階目の解析に青いチェックマーク（再利用...
・1段階目の解析結果が考慮さえているかの判断
Paravisで2段階目の解析で行った結果ファイルを開く。~
2段階目の解析結果のスタート地点の出力した値（DEPLとか）が...
つりあい解析を行った場合と行わない場合での結果を比較した...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
*その他のテクニック [#x988d31e]
ここでは,様々な解析手法を説明する.~
一部,千代岡の卒論日誌から引用し,編集したものを載せていま...
**部分腐食モデルの作り方 [#df1c52ea]
千代岡が使っていた部分腐食モデルの作成方法を載せる.~
部分腐食モデルは、腐食部分と健全部分を分けたモデルである...
2026年3月現在は、定着部の腐食を想定した位置に腐食部分を作...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/bubunc...
(1)いつもどおり、ビームモデルを作成する。（始点0,0,0、終...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/first....
(2)ポイントの作成→ポイントの真ん中のボタンを選択.作成した...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/second...
(3)(1)で作成したビームモデルと(2)で作成した点をパーティシ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/third....
(4)パーティションしたモデルをグループ化する。ここでは、腐...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/forth....
(5)mesh画面へ→メッシュ作成→1D:wire discretisationを選択→...
[余談；成功した試しはないが、今後サグ再現とかする場合は分...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/fifth....
(6)メッシュを切ると、写真のようになる。（写真は、Meshのグ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/sixth....
**要素番号や節点番号の探し方 [#w187d31e]
Meshで探すことができる。~
まず、作成した（計算した）メッシュを右クリックして,メッシ...
開くと,下の写真のような画面が出る。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
要素情報を選ぶと,3つの選択肢が出る。ノードと要素を選び,空...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
グループ化した要素についても、詳細情報が見れる。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**AFFE_CARA_ELEMについて（編集中） [#r8ad355d]
***AFFE_MODELEとAFFE_CARA_ELEMの対応コマンド [#n8e951b0]
ソリッド要素では,AFFE_MODELEを3Dと設定していたが,ビーム要...
AFFE_MODELEで設定したMODELISATIONはAFFE_CARA_ELEMのコマン...
対応しているものを以下にまとめたので、参考にしてください.
POUTRE：POU_D_E, POU_D_T, POU_D_TG, POU_D_T_GD, FLUI_STRU...
BARRE：BAR, CABLE_GAINE~
CABLE：CABLE, CABLE_POULIE~
COQUE：COQUE_AXIS, COQUE_C_PLAN, COQUE_D_PLAN, DKT, DST, ...
DISCRETE：DIS_T, DIS_TR~
DISCRETE_2D：2D_DIS_T, 2D_DIS_TR~
MULTI_FIBRE：POU_D_EM, POU_D_TGM~
→梁の種類によって、選択するものが異なる。自分がこれまで使...
・真っ直ぐな梁：POU_D_E(オイラー[Euler]梁), POU_D_T(ティ...
・マルチファイバーの梁(鉄筋コンクリートを想定したモデル作...
・棒(鉄筋とか断面が円のもの？) 2D_BARRE, BARRE 引張と圧...
・ケーブル CABLE 引張荷重 （引張だけ出せるが、斜張橋など...
・DIS_T,DIS_TR 離散要素（バネ要素） T[Translation]は変位,...
***断面特性の入力 [#e4c42938]
POUTRE→SECTION:GENERALEを選択した際,CARAの選択肢の意味を...
IZ GZと比較した主な幾何学的慣性モーメント~
IY GYに対する主慣性幾何モーメント~
AY GY方向のせん断係数~
AZ GZ方向のせん断係数~
EY ねじれ中心のオフセット（GYに対するCGの成分）~
EZ ねじれ中心のオフセット（GZによるCGの成分）~
JX ねじれ定数~
RY yに従って測定された外部ファイバーからの距離~
RZ zに従って測定された外部ファイバーからの距離~
RT トーションの有効光線~
JG 反りの定数（POU_D_TG、POU_D_TGM）~
IYR2 幾何剛性の計算に必要（POU_D_TGおよび POU_D_TGM)~
IZR2 幾何剛性の計算に必要（POU_D_TGおよびPOU_D_TGM)~
AI 梁内部の流体のバイパス部の表面~
***CABLEについて [#tb1796a6]
前にこの要素が橋梁モデルに使えるのか議論していたこと（未...
・SalomeのCABLEとは~
コンクリートのプレストレスなどを再現するのに使われるとの...
本研究室では,ケーブルに幾何学的非線形を入れた解析をする際...
・設定方法~
○AFFE_MODELEでの設定法~
MODELISATIONをCABLEに変更.
○AFFE_CARA_ELEMでの設定法~
N_INIT：初期張力の値を入力。しかし、橋梁のケーブルに引張...
→現状は,初期張力0と入力して,いつも使っているCREA_CHAMPを...
SECTION：断面積を値を入力~
FCX：DEFI_FONCTIONで作成した関数を入れる.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/setupC...
○DEFI_MATERIAUでの設定法~
CABLE→EC_SUR_E(直訳すると、圧縮弾性率と引張弾性率の寄与....
デフォルトで0.001.これを入れることで、圧縮をできる限り小...
***ピン接合の再現方法（DIS_TR[バネ要素]の設定法） [#v7995...
Marcでは.梁ピン接合という設定で再現できるが,Salomeではバ...
部材と部材の距離はできる限り短くする。(下の写真のモデルの...
DISCRETのK_TR_D_Lの設定画面~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
VALE[0],[1],[2]は変位に対する剛性(0,1,2の順にx,y,z方向を...
[3],[4],[5]は回転に対する剛性(3,4,5の順にDRX,DRY,DRZを示...
回転に対する剛性を緩め、変位に対する剛性を高めることでピ...
REPEREはGLOBAL（グローバル座標）でよい。
DISCRETのM_TR_D_Lの設定画面~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
VALE[0]〜[4]は質量行列に与える値です。質量行列を0にしない...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↓矢印説明書に書かれていた質量行列に関する内容を日本語訳し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**ビーム断面の適用方法 [#g75f9b28]
1D要素で弾塑性解析を始めとした材料非線形の解析を行う場合...
これらの要素にはビーム断面を作成してビームモデルに適用す...
ちなみに,MはMultifiber(マルチファイバー)という意味で、FEM...
[[ファイバー要素の説明はこちら:https://www.jepoc.or.jp/te...
(1)DEFI_GEOM_FIBER~
まず、このコマンドを使う前にビーム要素の断面として使う断...
SECTIONは断面定義、FIBREは断面の穴あき部分を定義する?(例...
ビーム断面が厚肉断面もしくは空白部分に何も入れる必要がな...
TOUT_SECTはビーム要素に作成した平面モデルを全て適用したい...
GROUP_FIBERは好きな名前をつけてよい。~
MAILAGE＿SECTはLIRE_MAILAGEで適用した平面モデルを選択。~
COOR_AXE_POUTREとANGLEは0でよい。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
(2)AFFE_CARA_ELEM~
GEOM_FIBREが適用されていることを確認して、MULTIFIBREの選...
(ただし、AFFE_MODELEで_*Mと書いてあるキーワードを設定する...
(3)DEFI_COMPOR~
ここでは断面の材料特性と状態(弾塑性を入れるとか)を定義す...
このコマンドは画面上の{a}+の横にある{}+(Show All Commands...
スクロールして探すなら、Otherの欄にあるので、それを選択し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
適用すると、下図のような画面が現れるが、GEOM_FIBREはDEFI_...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
MULTIFIBREを選択したら、GROUP_FIBREは(1)と同じ名前,mater...
最後に、DEFI_COMPORをAFFE_MATERIAUのAFFE_COMPORで適用する...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
**梁荷重の掛け方 [#vaf6c76b]
これまでFORCE_POUTREを使用していたが,載荷条件としてAFFE_C...
saika2 = AFFE_CHAR_MECA(FORCE_POUTRE=(_F(FY=-35250.0,
GROUP_MA=('keta'...
_F(FY=-122000.0,
GROUP_MA=('keta'...
_F(FY=-80000.0,
GROUP_MA=('CW', ...
_F(FY=-60000.0,
GROUP_MA=('syuto...
そのため、AFFE_CHAR_MECA1つにつき,載荷条件は1つしか与えて...
そこで、FORCE_NODALEを用いて梁荷重を再現する方法を見つけ...
1.Meshのグループ作成で,梁荷重として載荷したい要素のノード...
2.AFFE_CHAR_MECA→FORCE_NODALEを選択→載荷値は(分布荷重×要...
例：モデルに死荷重を掛けたい。（桁長：300m,2m間隔で要素分...
(50000×300)/150 = 100000N~
この値を入力する。
**動的応答解析のやり方 [#t4e8361f]
・注意点~
(1)DEFI_FONCTIONとAFFE_CHAR_MECAは、Pre Analysis・Analysi...
(2)DISCRETに質量マトリクスを用意しておくこと。DISCRETには...
この設定をしないと、ASSEMBLAGEでバネ要素に質量マトリクス...
[[詳しくはこちらを参照:https://code-aster.org/forum2/view...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)DYNA_VIBRAとCALC_CHAMPの名前は同じにすること。(同じ名...
(4)CALC_CHAMPのEXCITにDYNA_LINEで設定したCHARGEの分だけ入...
設定方法~
[1]Pre analysisのASSEMBLAGEを選択~
・ここでやること~
①NUME_DDLは{}押し、nddlと入力。~
②MATR_ASSEで質量マトリクスと剛性マトリクスを定義する。(質...
③CHARGEはAFFE_CHAR_MECAで適用した境界条件分だけ追加する。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
[2]DYNA_NON_LINE~
ここで弾塑性を入れた動的解析ができると思われるが、そのあ...
・最低限やること
①SCHEMA TEMPSでNEWMARK,FORMULATIONをDEPLACEMENT(変位),β=0...
ちなみに,[[説明書:http://dugi6514.odns.fr/doc/v15/en/man_...
更に詳しく調べると,これは直接積分法（正確には直接時間積分...
（Marcもニューマークβ法を用いて連鎖崩壊に関する解析をして...
[[参考サイト:http://jikosoftcom.blog25.fc2.com/blog-entry...
②INCREMENTはDEFI_LIST_INSTで設定したものを選ぶ。~
③EXCITにAFFE_CHAR_MECAで作成した数だけ設定する。載荷条件...
ビーム要素なら、CARA_ELEMにチェックを入れること。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**陰解法と陽解法について [#qe1cdcdd]
陰解法は時間変化を取り扱う解析で、現在の時間TからT+ΔTに進...
特徴として...
-非線形性が強い問題では収束しないことがあるため、静的や準...
-各step(時刻)で連立方程式を解く必要があるため、解析に時間...
陽解法は現在の時間Tをもとに、時刻T+ΔTの解を求める解法.~
特徴として...
-非線形性が強い問題を扱いやすいため、衝突や落下などの短い...
-時間増分ΔT（timestep）を[[クーラン条件:https://www.cyber...
→どうやら非常に細かいtimestepにしないといけないため、こち...
Salome-Mecaでは,DYNA_NON_LINE→SCHEMA_TEMPで設定可~
陰解法：NEWMARK,HHT~
陽解法：DIFF_CENT,TCHAMWA~
***各解法の説明 [#k9ed93bb]
-陰解法
・NEWMARK~
時間積分法の一種.平均加速度法の係数β=0.25,γ=0.50の組み合...
・HHT~
①修正平均加速度:組み合わせ「HHT,MODI_EQUI = 'NOT'」~
②Hilber-Hughes-Taylor法：組み合わせ「HHT,MODI_EQUI = 'OUI...
-陽解法
・DIFF_CENT~
・TCHAMWA~
**ある時間で力を0に落とす方法 [#s666447f]
FEMにおいて,要素をある時間で切断・消すことはできない。（...
そこで,力を0にしたい要素をあらかじめ消しておいて,元々要素...
ケーブル2本で桁を吊っている簡易モデルを例に説明する。解析...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
＜下準備＞~
(1)まず、破断させたい要素の軸力(EFGE_ELNO)をデータ化して...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(2)110sまで画面を移動させ,110.001sで=0を入力。そして,解析...
画像では「,」ついてますが、この時点では「,」要らないので...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)時間データを違うセルに貼り付け、その隣に,=","&並べた張...
「,」を入れないと、Salomeにデータを入れることができないた...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
ちなみに、0.005s掛けて張力を0に落としたい場合、破断初めの...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(4)作成したデータをテキストエディタに貼り付け、それを保存...
＜Salomeでの設定方法＞
(1)Meshの画面に行き,破断させたい要素を取り除いたcompound_...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(2)DEFI_FONCTION→VALEを選択。import tableを選択.（表に矢...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)import tableを選択すると,ファイルを選択する画面が出る...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(4)BC and Load → FORCE_NODALEを選択.~
ここで,元々要素があった向きに軸力を与える必要がある。元々...
要素をなくした部分を斜辺として,1つの三角形として表現する...
シンプルな方法として,三平方の定理を用いて斜辺を求め,底辺/...
今回は1:1:√2の三角形のため,1/√2の力の分解となる。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
そして,(2)で与えたDEFI_FONCTIONで与えた関数を与えることで...
**温度でのプレストレスの与え方（編集中） [#vac22495]
ここでは温度でのやり方を示す。（ひずみは後述のPRE_EPSIに...
(0)DEFI_FONCTIONで温度変化のfunctionを作成する。（MECA_ST...
NOM_PARA（横軸）をINST,NOM_RESU（縦軸）をTEMPに設定。~
PROL_はVALEで作成したグラフの範囲外（DROITEは右側,GAUCHE...
写真はEXCLUを使っているが,これは除外という意味。（VALEの...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
VALEの設定方法はSTAT_NON_LINE(静的解析)orDYNA_NON_LINEと...
写真は静的解析を想定したもので,0倍の変化時に温度0、1倍の...
時刻歴の場合だと、解析時間0秒のときに温度0℃、解析時間1秒...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
ちなみに...張力から温度を逆算する場合は、Δt = T/EAα を用...
(1)Post Processing → CREA_CHAMPを選択。~
TYPE_CHAMは「NOEU_TEMP_F」を選択する。（MECA_STATICの場合...
OPERATIONは「AFFE」,MODELEにチェックを入れる。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
AFFEを1itemにする.温度荷重を与えたい部材をGROUP_MAで選択.~
VALEはMECA_STATICSを用いる場合,ここで最終的に変化させたい...
VALE_FはMECA_STATICS以外で用い,(0)で作成した温度のfunctio...
NOM_CMPはTEMPと入力.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(2)Output → CREA_RESUを選択。（MECA_STATICSの場合は必要な...
この設定をしないと温度荷重を考慮して解析を回してくれない...
(3)AFFE_MATERIAU → AFFE_VARCを選択。~
**ひずみでのプレストレスの与え方[#c71fb860]
（追記）~
こちらは、プレストレスをひずみで与えるやり方です.~
簡易モデルでは、この方法でプレストレスを与えたが、PRE_EPS...
PRE EPSIは「設定したひずみになるように外力を与える」コマ...
それを回避するために,要素分割し,分割した1要素に以下の式で...
PRE_EPSIをプレストレスとして用いる際の式：ε = α T × 要素...
[[PRE_EPSIの説明書:http://dugi6514.odns.fr/doc/v15/en/man...
モデルをひずみを用いて変形させることができるコマンドであ...
1Dモデル(Multifiber要素も含む)ではEPX（梁のローカル軸に応...
(ky,kzはたわみの微分方程式から得られるdθ/dx = -M/EI の値...
X,Y,Zの軸の向きはグローバル座標(salome画面の色がついてる...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
単純梁モデル(10✕10✕120mm)にEPXX：0.01を与えると、以下の挙...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
線に色がついてるところが変形前のモデルの端であり、写真を...
・簡易モデルに用いた際に分かったこと~
現在、ケーブルにPRE_EPSIを入れており、縮む挙動(ケーブルが...
その理由について考えてみたが、要素を構成している節点にモ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑PRE_EPSIの意味 日本語に訳すと要素に変形の負荷を与えるみ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑簡単な例 PRE_EPSIで与えたひずみになるような外力をかけて...
**CALC_CHAMP [#t08db3b9]
***POU_D_E,POU_D_EMで計算できるもの（作成中） [#n061551b]
**IMPR_RESU [#m5434190]
***IMPR_RESUのGROUP_MA,GROUP_NO [#wafb540f]
IMPR_RESUのRESUの中に入っているGROUP_MA,GROUP_NOを使うこ...
これを使うメリットは容量の大きいモデルで解析を回す際に,見...
*影響線解析 [#h9310a98]
影響線：移動荷重に対してある点に生じる物理量の大きさを移...
影響線の定義から,まずは移動荷重を求め方が分からないと話に...
↓橋の上を走る電車を例に移動荷重を算出する方法を見つけた。...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
[[移動荷重について:https://forum.code-aster.org/public/d/...
[[移動荷重について_その2:https://forum.code-aster.org/pub...
AFFE_CHAR_MECA_F→FORCE_POUTREを使い、DYNA_NON_LINEなどの...
ちなみに...midasだと[[説明書の10章:https://in.midasit.co....
***軸力の見方 [#df51b2bf]
軸力はEFGE_ELNOの[N]で見る事ができることが分かった。Nは法...
見る方法だが、medファイルをParaViSで閲覧し、フィルター→me...
[[参考資料:http://dugi6514.odns.fr/doc/v15/en/man_u/u2/u2...
↑Salomeで出力される力やひずみなどについて詳しく書かれてい...
**断面定数の算出方法 [#g12212f4]
複雑な断面をビームモデルに適用したい場合,異なるhdfファイ...
ここでは,そのやり方を説明します.
・（橋梁の主塔のような内側に穴が空いている矩形断面の作成...
高さ2m,幅4m,厚さ0.025mの矩形断面を例に説明します。
○Geometry~
・平面を用いて断面を作る~
長方形のフェースを作成をクリック.ここで,高さor幅＋厚さの...
（big：高さ2m+厚さ0.025m,幅4m+0.025m, small：高さ2m−厚さ0...
断面を作成する際,座標の向きはX座標下向き,Y座標右向きで作...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
作成したら,Cutを用いて,厚さ0.025mの矩形断面を作成する.2種...
ここでは,bigをメインオブジェクト,smallをツールオブジェク...
後,下の写真では高度のオプションに2種類チェックを入れてし...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
Cutすると,以下の写真のようになる.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・外側の四辺、内側の四辺、外側の四辺上の角の点、内側の四...
外側の点と内側の点,外側の線と内側の線,計16個全てをグルー...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↓こんな感じになるかと.グループの名前は分かりやすいものに...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
○Mesh~
・メッシュを作成~
アルゴリズムはNETGEN1D-2D,詳細設定は歯車押して,適用すれば...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・ジオメトリグループの作成~
点はノード,線は要素に入れる.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・グループを作成~
Faceにチェックを入れ,全て選択します。を選ぶ.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・メッシュをmedファイルでエクスポート~
作成したMeshを右クリックして,エクスポート→MEDファイルを選...
○aster_study~
・LIRE_MAILAGEでメッシュ適用~
エクスポートしたMEDファイルを適用する.やり方は,…マークを...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・Model Definition > MACR_CARA_POUTRE >「GROUP_MA_BORD」...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
MACR_CARA_POUTREは1つの辺に対し,1つの点しか適用できないた...
面倒だと思うけど...一つずつ設定してください.~
また,設定の仕方ですが,一筆書きになるように作成してくださ...
そして,外側の下辺(lsoto3),左辺(lsoto4)と設定すると一筆書...
なお,「GROUP_NO」については黄色の線に対応するもの（今の状...
画像の例だと,soto1→soto2→soto3→soto4になります.（上辺の左...
↓「GROUP_MA_BORD」の設定画面~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↓「GROUP_NO」の設定画面~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・Model Definition > MACR_CARA_POUTRE にて内側も同様に~
外側の設定と同じようにやればOK.
・IMPR_RESU で FormatをRESULTATに、ファイル名を「kekka.re...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
○計算を回してresultファイルを見る。~
・ファイルに書かれているA(断面積),IY,IZ(断面2次モーメント...
なお、空白エリアが存在しない断面（密な断面）を作成する際...
また、.resuファイルはhdfファイルやFilesと同じところに保存...
○メインとなるhdfファイル（橋梁モデル等,主となるモデルがあ...
・断面形状の入れ方~
LIRE_MAILAGE→...をクリック→エクスポートした断面形状のmed...
*材料非線形 [#ad9f0c30]
**梁要素に適用できる（できそう？）弾塑性モデル [#y6458c84]
[[Astercode_U4.51.11:https://code-aster.org/doc/v15/en/ma...
[[Astercode_R5.03.09:https://code-aster.org/doc/v15/en/ma...
梁要素に適用できる,できそうな弾塑性モデルをまとめた。~
1Dモデルに対応している弾塑性モデルは以下のとおりである。~
-VMIS_ISOT_LINE 対称線形等方加工硬化を伴うフォンミーゼス
-VMIS_ISOT_TRAC 非特定の等方性加工硬化を伴うフォンミーゼス
-VMIS_CINE_LINE 対称線形キネマティック加工硬化を伴うフォ...
-ECRO_CINE_1D 対称線形キネマティック加工硬化を伴うフォン...
-GRILLE_ISOT_LINE 対称線形等方性加工硬化を伴うフォンミー...
-GRILLE_CINE_LINE 対称線形キネマティック加工硬化を伴うフ...
-PINTO_MENEGOTTO 鉄筋コンクリート鉄筋の挙動
-GRILLE_PINTO_MEN 鉄筋コンクリート配筋の挙動
-VMIS_ASYM_LINE 非対称線形加工硬化および復元によるフォン...
-VISC_IRRA_LOG,GRAN_IRRA_LOG 燃料集合体の粘塑性挙動： REF...
-MAZARS MAZARSのバージョン1Dの挙動.
-RELAX_ACIER プレストレストケーブルのリリーフをモデル化す...
説明は皆さんが使うことがある、可能性があるものを以下に記...
-VMIS_ISOT_LINE,VMIS_CINE_LINE~
適用可：BAR,POU_D_EM,POU_D_TGM~
フォンミーゼスの弾塑性挙動と線形加工硬化の関係。~
材料に必要なデータは、オペレータDEFI_MATERIAU [U4.43.01]...
VMIS_ISOT_LINEは大変形SIMO_MIEHEに対応しているとのこと。~
-VMIS_ISOT_TRAC~
適用可：BAR,POU_D_EM,POU_D_TGM~
フォンミーゼスの弾塑性挙動と非線形等方加工硬化の関係。~
曲線σ,εは、DEFI_MATERIAU [U4.43.01]のTRACTIONで提供される...
温度に応じていくつかのトラクションダイアグラムを定義する...
また、DEFI_MATERIAU演算子のELAS (_FO)は必ず必要。~
トラクションダイアグラムが提供された場合、挙動の関係に使...
大変形SIMO_MIEHEに対応しているとのこと。~
-VMIS_ASYM_LINE~
適用可：BAR~
説明：等方性加工硬化を伴うフォンミーゼスの弾塑性一軸等温...
制御装置や埋設ケーブルと地盤との相互作用をモデル化するこ...
材料に必要なデータは、DEFI_MATERIAU [U4.43.01]のECRO_ASYM...
*モーダル解析 [#r5d7cbd2]
モーダル解析とは~
対象物そのものが振動しやすい周波数とその振動形状を把握す...
構造物の固有振動数と固有モードが見れるとのこと.(よく皆が...
固有振動数とは？何故、固有振動数を求めるのか？~
→物体ごとがもつ周波数(1秒間に揺れる振動回数)、構造物の共...
固有モードとは？~
→それぞれの固有振動数においてどのように振動しているかの変...
片持梁を想定したモーダル解析コマンドファイル~
PreanalysisのASSEMBLAGEとAnalysisのCALC_MODEを使うとのこ...
mesh = LIRE_MAILLAGE(identifier='0:1',
FORMAT='MED',
UNITE=20)
mesh0 = LIRE_MAILLAGE(identifier='1:1',
FORMAT='MED',
UNITE=3)
model = AFFE_MODELE(identifier='2:1',
AFFE=_F(GROUP_MA=('beam', ),
MODELISATION='POU_D_EM',
PHENOMENE='MECANIQUE'),
MAILLAGE=mesh)
fiber = DEFI_GEOM_FIBRE(identifier='3:1',
SECTION=_F(COOR_AXE_POUTRE=(0.0, ...
GROUP_FIBRE='men',
GROUP_MA_SECT=('danmen...
MAILLAGE_SECT=mesh0))
elemprop = AFFE_CARA_ELEM(identifier='4:1',
GEOM_FIBRE=fiber,
MODELE=model,
MULTIFIBRE=_F(GROUP_FIBRE=('men...
GROUP_MA=('beam',...
PREC_AIRE=0.76,
PREC_INERTIE=0.95),
POUTRE=_F(CARA=('HY', 'HZ'),
GROUP_MA=('beam', ),
SECTION='RECTANGLE',
VALE=(10.0, 20.0)))
mater = DEFI_MATERIAU(identifier='5:1',
ELAS=_F(E=206000.0,
NU=0.3,
RHO=7.9e-09))
load = AFFE_CHAR_MECA(identifier='6:1',
DDL_IMPO=_F(DRX=0.0,
DRY=0.0,
DRZ=0.0,
DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_NO=('kotei', )),
MODELE=model)
behav = DEFI_COMPOR(identifier='7:1',
GEOM_FIBRE=fiber,
MATER_SECT=mater,
MULTIFIBRE=_F(GROUP_FIBRE=('men', ),
MATER=mater,
RELATION='ELAS'))
fieldmat = AFFE_MATERIAU(identifier='8:1',
AFFE=_F(GROUP_MA=('beam', ),
MATER=(mater, )),
AFFE_COMPOR=_F(COMPOR=behav,
GROUP_MA=('beam',...
MODELE=model)
ASSEMBLAGE(identifier='9:1',
CARA_ELEM=elemprop,
CHAM_MATER=fieldmat,
CHARGE=(load, ),
MATR_ASSE=(_F(MATRICE=CO('mass'),
OPTION='MASS_MECA'),
_F(MATRICE=CO('rigi'),
OPTION='RIGI_MECA')),
MODELE=model,
NUME_DDL=CO('nddl'))
modes = CALC_MODES(identifier='12:1',
CALC_FREQ=_F(NMAX_FREQ=10),
MATR_MASS=mass,
MATR_RIGI=rigi,
OPTION='PLUS_PETITE',
VERI_MODE=_F(STOP_ERREUR='OUI'))
IMPR_RESU(identifier='13:1',
FORMAT='MED',
RESU=_F(CARA_ELEM=elemprop,
NOM_CHAM=('DEPL', ),
RESULTAT=modes),
UNITE=4)
非線形性を考慮するには,DYNA_NON_LINEを使うらしい。要調査
結果比較~
,,鉛直1次モード,水平1次モード,
,結果,7.96Hz,16.35Hz,
,理論値,8.22Hz,16.45Hz,
理論値と近い値が出たため、ビーム要素におけるモーダル解析...
*エラーメッセージ一覧 [#m97c0c8e]
・このエラーメッセージは構造的に不安定なときに出る。~
これが出た場合は基本的にAsterstudyの設定を間違えているこ...
それでもダメなら,ファイバー要素を適用している断面のメッシ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・AFFE_CHAR_MECAの設定をモデルに適用できない場合に出る。...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・AFFE_MODELEとAFFE_CARA_ELEMの条件が一致しない場合に出る...
AFFE_MODELEで設定したMODELISATIONは、それぞれAFFE_CARA_EL...
(例) AFFE_MODELEでBARREを選択した場合→AFFE_CARA_ELEMではB...
AFFE_MODELEでPOU_D_Eを選択した場合→AFFE_CARA_ELEMで...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**収束エラーが出た場合の対処法（マルチファイバー要素を使...
収束エラーや行列絡みのエラーが出たときは、まずtimestepや...
ファイバー要素で作成した片持ち梁モデルで断面メッシュの細...
解析内容は片持ち梁の静的材料非線形解析である。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
***三角形メッシュの場合 [#m07ec30d]
載荷条件だが、与えた条件より値が小さくてもモデルが降伏し...
変えたのは断面の要素数のみで他の設定は変えていない。~
この条件では,要素数156では行列絡みのエラーが出たため解析...
ちなみに,解析結果だが,どの部分も設定した降伏応力を超える...
また,固定面を見ると,綺麗に引張,圧縮と分かれている.(どちら...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
解析結果を見て,分かったことがあったので、ここにまとめる。~
・ビーム要素(線要素)の解析において,一気に変形する挙動を起...
以下のアニメーションを見ると分かりやすいと思うが,大きな変...
実は,ちょうど解析エラーで止まった場所はこの大きな変化が起...
もしかすると,材料非線形を考慮した解析でエラーが出てしまう...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/ex1.gif
***四角形メッシュの場合 [#e01f2635]
しかし、三角形メッシュだと要素数が増やさないと解析が回ら...
そこで、三角形メッシュと同じ解析条件のまま、メッシュの形...
四角形メッシュの変え方はNETGEN1D-2Dを選択し,メッシュサイ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
四角形メッシュの要素数：100~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
この状態でも解析は回った。三角形メッシュでは解析が回らな...
応力をみると、こちらも設定した降伏応力を超えることはなか...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
***結果比較とまとめ [#g5773cf3]
とりあえず、先端変位と固定面の応力ひずみを確認した。
最終変位~
三角形:6.06mm 四角形：2.31mm~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
理論値は1.17mmで,結果を見ると,三角形メッシュの場合は1.30m...
グラフを見る限り,荷重が掛かり切ると,梁が大きくたわむこと...
応力ひずみ~
三角形メッシュの応力ひずみは以下の通り。中立軸データはグ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
四角形メッシュの応力ひずみは以下の通り。（要素数が少ない...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
グラフを見る限り、綺麗な三角形応力分布ができていると思わ...
まとめ~
どちらのメッシュの切り方でも応力ひずみ関係はほぼ同じ挙動...
これは,恐らく要素数の違いが原因だろう.なので,四角形メッシ...
ただ,たわみの理論値を少し超えた付近で線形関係から外れるた...
しかし、橋梁モデルのような大きいモデルを解析する際に三角...
以上から、大きめなモデルを解析する際は,要素数少なめで解析...
現在の閲覧人数: &online;~
総閲覧数:&counter;~
終了行:
#contents
Salome-Mecaのビーム要素（梁要素）に関するページを作成しま...
これから梁要素を用いた解析を扱うよって人は見ていただける...
また、何か分かったこととかあればこちらに書き込んで共有し...
*参考資料 [#mb640bf5]
[[code-asterによる解析 step-by-step:https://www.opencae.o...
[[結果の出力絡みのお話（CALC_CHAMP,IMPR_RESU）:https://ww...
[[構造要素の概要:https://www.opencae.or.jp/wp-content/upl...
[[code-asterの説明書:https://code-aster.org/doc/v15/en/in...
[[Beginning with code_aster:https://archives.framabook.or...
[[code-aster質問ページ:https://forum.code-aster.org/publi...
*ビーム要素(梁要素)とは [#z866693c]
ビーム要素（梁要素）とは点と点を線で結ぶ要素である.ソリッ...
例えば,橋梁モデルみたいな巨大なモデルで時刻歴応答解析をす...
そのため、モデルは全て梁要素、もしくは桁と塔をシェル要素...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
*単位系の注意点 [#k7f047ce]
Salome-Mecaは単位系について自分で設定しないといけません。~
モデル長mm想定なら...ヤング率(MPa),荷重(N),応力(MPa)~
m想定なら...ヤング率(Pa),荷重(N),応力(Pa)~
*気になること [#q683524a]
-非線形絡みの問題~
非線形には材料非線形、幾何学的非線形、境界非線形の3種類存...
今は弾性域内（線形で片付く範囲）でしか解析してない。一応...
→幾何学的非線形はcable要素でDeformation:GROT_GDEPを用いる...
-橋梁モデルを作れるのか~
簡易モデル（[[千代岡の卒論日誌参照:https://www.str.ce.aki...
→作成可.ケーブルに材料非線形を考慮しないなら、ケーブルをA...
-座屈解析はできるのか~
ビーム要素で座屈解析できるのだろうか？~
もし、試すなら、4年生の春課題で行なったモデルをビーム要素...
-固定条件について~
DDL_IMPOにBLOCAGEというコマンドがあり,これは選択した項目,...
しかし,例えば完全固定条件の1.DX=DY=DZ=DRX=DRY=DRZ=0と2.BL...
両者の違いは何だろうか？
-載荷条件について~
AFFE_CHAR_MECAで載荷条件の順番によって得られる結果が変わ...
[[事例:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/...
→AFFE_CHAR_MECAで使う載荷条件をそれぞれ1つずつ分けると改...
saika1 = AFFE_CHAR_MECA(FORCE_NODALE=_F(FY=-1175000.0,
GROUP_NO=('saika'...
MODELE=allmodel)
saika2 = AFFE_CHAR_MECA(FORCE_POUTRE=(_F(FY=-35250.0,
GROUP_MA=('keta'...
_F(FY=-122000.0,
GROUP_MA=('keta'...
_F(FY=-80000.0,
GROUP_MA=('CW', ...
_F(FY=-60000.0,
GROUP_MA=('syuto...
MODELE=allmodel)
saika4 = AFFE_CHAR_MECA(FORCE_NODALE=_F(FY=-9600000.0,
GROUP_NO=('CW_ten...
MODELE=allmodel)
追記：調査の結果、同じAFFE_CHAR_MRCAに2つ以上の載荷条件を...
[[参考資料:https://forum.code-aster.org/public/d/15289-re...
解決策を見つけたので,[[こちらを参考:https://www.str.ce.ak...
**ばねを用いたピン接合について [#fb9a4942]
Salomeで梁要素を用いてピン接合の再現をするにはばね要素を...
しかし、現状として解析が回るモデルと回らないモデルが混在...
（写真のような,横ばねが存在しているモデルとエラーが出る可...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↓メッセージ内容~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
現在、解決策を模索中であるが、検証をしたものをここにまと...
・バネ要素の両端（部材との共有節点）に自由度の設定をする。~
上の写真の場合、ばねの右端の共有節点DRXの自由度エラーが出...
**材料非線形に関する問題 [#h9737f8e]
***CREA_RESUの設定 [#vcda802b]
材料非線形を考慮する場合,CREA_RESU→COMPORTEMENTの設定をし...
[[説明書:https://code-aster.org/doc/v14/fr/man_u/u4/u4.44...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/CREA_R...
*モデルの作成・片持ち梁を例にした解析手法（編集中） [#p87...
100mm✕10mm✕10mmモデルを対象に説明します.（これまで先輩方...
モデルを片持ち梁として,変位を見てみましょう。
**Geometry [#ud86b1ec]
(1)ポイント(0,0,0)(100,0,0)を作成する.
(2)線（ポイントの右隣にあるアイコンをクリック）を作成する...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)(POU_D_EMを適用する場合のみ)~
ビーム断面を作成する.材料非線形を考慮した解析や複雑な断面...
これはAsterstudyでビームモデルに断面を適用するのに必要な...
長方形断面を作成する場合は画面上部にある長方形マーク（長...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
円形断面を作成する場合は長方形マークの右隣の楕円マーク（...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
方向はモデル軸方向をxにとっているため,OXYにする.そうしな...
（構造研スタイルと同じようにモデル軸方向をzにした場合,方...
モデルとしては,こんな感じになるはず.（作成したのは正方形...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**Mesh [#ed8a38b1]
(1)ビームモデル（line）のメッシュを切る.1D→Wire Discretis...
セグメントの数は分割数を示しているので,任意の数値を入れる.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(2)(POU_D_EMを適用する場合のみ)~
断面のメッシュは2D→NETGEN1D-2D,歯車マークからNETGEN2DPara...
最大・最小サイズの数値を一致させ,任意の数値を入れる.ただ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)メッシュを計算した後,グループを作成する.~
要素の型は,基本的にノード：点,Edge:線,Face:面の3種類を使...
手動で選択する場合,コンテンツにチェックを入れた後,グルー...
全てを選択するは,断面をグループ化するときに使う.
今回は,固定点と載荷点,線要素（ビーム）全て,（断面作成して...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**Asterstudy [#o2e2746d]
(数字-1)の箇所は断面を作成した場合(POU_D_EMを用いる場合)...
(1)LIRE_MAILAGEを選び,計算したメッシュ(ビームモデル)を選...
(2)AFFE_MODELEを選び,MAILAGEはビームモデルのものを選択す...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
梁要素で用いるAFFE_MODELEは[[こちら:https://www.str.ce.ak...
(2-1)DEFI_GEOM_FIBERを選択し,SECTION（断面定義）を選ぶ.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
TOUT_SECTorGROUP_MAで作成した断面を選ぶ。~
GROUP_FIBERは好きな名前をつけてよい。（ここでは,sectionと...
MAILAGE_SECTはLIRE_MAILAGEで適用した断面を選択。~
COOR_AXE_POUTREは0で。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)AFFE_CARA_ELEMを選択.ビーム要素において,ここでの設定は...
POUTREを選び,GROUP_MAをビームモデル,SECTIONをRECTANGLE（...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3-1)MULTIFIBERを選択.ここで,(2)で設定した断面定義をモデ...
(4)DEFI_MATERIAU→ELASで材料特性選択.鋼材を想定しているた...
(5)AFFE_MATERIAU→MODELEにチェック入れ,AFFE選択→(4)で作成...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(5-1){}+を選択→DEFI_COMPORを選択.選択すると,下画面のよう...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEFI_C...
MULTIFIBERを1itemにし,GROUP_FIBERを(2-1)で記入した名前に...
RELATIONの説明は...後でまとめます.(基本的にELAS(線形弾性)...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEFI_C...
(5-2)
(6)AFFE_CHAR_MECAを選択.固定条件と載荷条件は別々に分けて...
まず,固定条件はDDL_IMPOを選択し,GROUP_NO（グループ化[Grou...
DX=DY=DZ=DRX=DRY=DRZ=0と入れてもよいが,ENCASETR[完全固定]...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
載荷条件はFORCE_NODALEを選択.（ここで入れる値は[N],他のFO...
GROUP_NO（グループ化[Group]した点[Node]）に載荷点を入れ,...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(7)AnalysisはMECA_STATICSを選択（POU_D_EMには対応していな...
CHAM_MATERとCARA_ELEMにチェック入れ,EXCITにAFFE_CHAR_MECA...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(7-1)AnalysisはSTAT_NON_LINEを選択(POU_D_Eには対応してい...
(8)IMPR_RESUで.medファイルを作成,FORMATはMEDで！~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
RESUではRESULTATにAnalysisで設定したもの（MECA_STATICS or...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
なお,今回は変位のみを見ることが目的のためCALC_CHAMPは使わ...
**Paravis [#ed37e359]
Asterstudyで解析したmedファイルを開く.適用する前にGenerat...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
その後,Warp by Vectorを選択して変位図を見て,spread sheet ...
*つりあい解析（二段階解析）のやり方 [#gac275b5]
つりあい解析とは、事前に解析した構造解析の結果（1段階目の...
横軸時間をとる時刻歴解析やモデルを先に変形させてから力を...
大きなメリットとして,解析時間の短縮が挙げられる.そして,段...
・使用モデル~
卒論で扱ったケーブル1本の簡易モデルからケーブル(C1)をなく...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・解析条件~
荷重は死荷重とケーブルから得られた張力を掛けている.~
時刻歴応答解析 解析時間は115s,100sまで荷重等をゆっくり載...
109sまでを1段階目,109sから2段階目の解析を行う。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・設定法~
写真のようにStageを2種類用意する。~
chainanalysisが1段階目,chainanalysis-2が2段階目用のStage...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
1段階目 commファイルは[[こちら:http://www.str.ce.akita-u....
基本的に自分が行っている解析のaster-codeをいつもどおり組...
時刻歴解析の場合,区切る時間（今回は109s）までの解析をここ...
2段階目 commファイルは[[こちら:http://www.str.ce.akita-u....
AFFE_MODELEやDEFI_MATERIAUなど1段階目の解析と共通するもの...
Stage1で設定したものがそのまま使えます。
使用するaster-code~
絶対：Analysisに含まれるもの ,Output IMPR_RESU~
人によって変わる：Functions and Listsに含まれるもの （時...
ここでは時刻歴応答解析を行っているため,以下のように設定し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
1段階目の解析終了地点から解析を再開しないと,1段階目の結果...
そのため,解析時間は109sからスタートする。
Analysisでは,それぞれのaster-codeに入っているコマンド：ET...
（静的解析を行う場合,STAT_NON_LINEで。MECA_STATICSにはそ...
また、MODELEやCARA_ELEMといった選択式で入力するところはSt...
RESULTATに1段階目のAnalysisで用いたものを入れる。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
EVOL_NOLIに1段階目のAnalysisで用いたものを入れる。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
OUTPUTは1段階目と異なる.medの名前をつける。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
RESUにはNOM_CHAMの設定とTOUT_OREDEとTOUT_CMPをOUIにする。
・解析の回し方
写真のようにstage同士をつなげることでつりあい解析が可能と...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
初めて回す際、一気に2段階目まで解析を行おうとすると,2段階...
まず,1段階目の解析を回したいので,2段階目の解析に赤い四角...
解析が回ったら、1段階目の解析に青いチェックマーク（再利用...
・1段階目の解析結果が考慮さえているかの判断
Paravisで2段階目の解析で行った結果ファイルを開く。~
2段階目の解析結果のスタート地点の出力した値（DEPLとか）が...
つりあい解析を行った場合と行わない場合での結果を比較した...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
*その他のテクニック [#x988d31e]
ここでは,様々な解析手法を説明する.~
一部,千代岡の卒論日誌から引用し,編集したものを載せていま...
**部分腐食モデルの作り方 [#df1c52ea]
千代岡が使っていた部分腐食モデルの作成方法を載せる.~
部分腐食モデルは、腐食部分と健全部分を分けたモデルである...
2026年3月現在は、定着部の腐食を想定した位置に腐食部分を作...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/bubunc...
(1)いつもどおり、ビームモデルを作成する。（始点0,0,0、終...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/first....
(2)ポイントの作成→ポイントの真ん中のボタンを選択.作成した...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/second...
(3)(1)で作成したビームモデルと(2)で作成した点をパーティシ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/third....
(4)パーティションしたモデルをグループ化する。ここでは、腐...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/forth....
(5)mesh画面へ→メッシュ作成→1D:wire discretisationを選択→...
[余談；成功した試しはないが、今後サグ再現とかする場合は分...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/fifth....
(6)メッシュを切ると、写真のようになる。（写真は、Meshのグ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/sixth....
**要素番号や節点番号の探し方 [#w187d31e]
Meshで探すことができる。~
まず、作成した（計算した）メッシュを右クリックして,メッシ...
開くと,下の写真のような画面が出る。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
要素情報を選ぶと,3つの選択肢が出る。ノードと要素を選び,空...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
グループ化した要素についても、詳細情報が見れる。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**AFFE_CARA_ELEMについて（編集中） [#r8ad355d]
***AFFE_MODELEとAFFE_CARA_ELEMの対応コマンド [#n8e951b0]
ソリッド要素では,AFFE_MODELEを3Dと設定していたが,ビーム要...
AFFE_MODELEで設定したMODELISATIONはAFFE_CARA_ELEMのコマン...
対応しているものを以下にまとめたので、参考にしてください.
POUTRE：POU_D_E, POU_D_T, POU_D_TG, POU_D_T_GD, FLUI_STRU...
BARRE：BAR, CABLE_GAINE~
CABLE：CABLE, CABLE_POULIE~
COQUE：COQUE_AXIS, COQUE_C_PLAN, COQUE_D_PLAN, DKT, DST, ...
DISCRETE：DIS_T, DIS_TR~
DISCRETE_2D：2D_DIS_T, 2D_DIS_TR~
MULTI_FIBRE：POU_D_EM, POU_D_TGM~
→梁の種類によって、選択するものが異なる。自分がこれまで使...
・真っ直ぐな梁：POU_D_E(オイラー[Euler]梁), POU_D_T(ティ...
・マルチファイバーの梁(鉄筋コンクリートを想定したモデル作...
・棒(鉄筋とか断面が円のもの？) 2D_BARRE, BARRE 引張と圧...
・ケーブル CABLE 引張荷重 （引張だけ出せるが、斜張橋など...
・DIS_T,DIS_TR 離散要素（バネ要素） T[Translation]は変位,...
***断面特性の入力 [#e4c42938]
POUTRE→SECTION:GENERALEを選択した際,CARAの選択肢の意味を...
IZ GZと比較した主な幾何学的慣性モーメント~
IY GYに対する主慣性幾何モーメント~
AY GY方向のせん断係数~
AZ GZ方向のせん断係数~
EY ねじれ中心のオフセット（GYに対するCGの成分）~
EZ ねじれ中心のオフセット（GZによるCGの成分）~
JX ねじれ定数~
RY yに従って測定された外部ファイバーからの距離~
RZ zに従って測定された外部ファイバーからの距離~
RT トーションの有効光線~
JG 反りの定数（POU_D_TG、POU_D_TGM）~
IYR2 幾何剛性の計算に必要（POU_D_TGおよび POU_D_TGM)~
IZR2 幾何剛性の計算に必要（POU_D_TGおよびPOU_D_TGM)~
AI 梁内部の流体のバイパス部の表面~
***CABLEについて [#tb1796a6]
前にこの要素が橋梁モデルに使えるのか議論していたこと（未...
・SalomeのCABLEとは~
コンクリートのプレストレスなどを再現するのに使われるとの...
本研究室では,ケーブルに幾何学的非線形を入れた解析をする際...
・設定方法~
○AFFE_MODELEでの設定法~
MODELISATIONをCABLEに変更.
○AFFE_CARA_ELEMでの設定法~
N_INIT：初期張力の値を入力。しかし、橋梁のケーブルに引張...
→現状は,初期張力0と入力して,いつも使っているCREA_CHAMPを...
SECTION：断面積を値を入力~
FCX：DEFI_FONCTIONで作成した関数を入れる.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/setupC...
○DEFI_MATERIAUでの設定法~
CABLE→EC_SUR_E(直訳すると、圧縮弾性率と引張弾性率の寄与....
デフォルトで0.001.これを入れることで、圧縮をできる限り小...
***ピン接合の再現方法（DIS_TR[バネ要素]の設定法） [#v7995...
Marcでは.梁ピン接合という設定で再現できるが,Salomeではバ...
部材と部材の距離はできる限り短くする。(下の写真のモデルの...
DISCRETのK_TR_D_Lの設定画面~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
VALE[0],[1],[2]は変位に対する剛性(0,1,2の順にx,y,z方向を...
[3],[4],[5]は回転に対する剛性(3,4,5の順にDRX,DRY,DRZを示...
回転に対する剛性を緩め、変位に対する剛性を高めることでピ...
REPEREはGLOBAL（グローバル座標）でよい。
DISCRETのM_TR_D_Lの設定画面~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
VALE[0]〜[4]は質量行列に与える値です。質量行列を0にしない...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↓矢印説明書に書かれていた質量行列に関する内容を日本語訳し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**ビーム断面の適用方法 [#g75f9b28]
1D要素で弾塑性解析を始めとした材料非線形の解析を行う場合...
これらの要素にはビーム断面を作成してビームモデルに適用す...
ちなみに,MはMultifiber(マルチファイバー)という意味で、FEM...
[[ファイバー要素の説明はこちら:https://www.jepoc.or.jp/te...
(1)DEFI_GEOM_FIBER~
まず、このコマンドを使う前にビーム要素の断面として使う断...
SECTIONは断面定義、FIBREは断面の穴あき部分を定義する?(例...
ビーム断面が厚肉断面もしくは空白部分に何も入れる必要がな...
TOUT_SECTはビーム要素に作成した平面モデルを全て適用したい...
GROUP_FIBERは好きな名前をつけてよい。~
MAILAGE＿SECTはLIRE_MAILAGEで適用した平面モデルを選択。~
COOR_AXE_POUTREとANGLEは0でよい。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
(2)AFFE_CARA_ELEM~
GEOM_FIBREが適用されていることを確認して、MULTIFIBREの選...
(ただし、AFFE_MODELEで_*Mと書いてあるキーワードを設定する...
(3)DEFI_COMPOR~
ここでは断面の材料特性と状態(弾塑性を入れるとか)を定義す...
このコマンドは画面上の{a}+の横にある{}+(Show All Commands...
スクロールして探すなら、Otherの欄にあるので、それを選択し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
適用すると、下図のような画面が現れるが、GEOM_FIBREはDEFI_...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
MULTIFIBREを選択したら、GROUP_FIBREは(1)と同じ名前,mater...
最後に、DEFI_COMPORをAFFE_MATERIAUのAFFE_COMPORで適用する...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
**梁荷重の掛け方 [#vaf6c76b]
これまでFORCE_POUTREを使用していたが,載荷条件としてAFFE_C...
saika2 = AFFE_CHAR_MECA(FORCE_POUTRE=(_F(FY=-35250.0,
GROUP_MA=('keta'...
_F(FY=-122000.0,
GROUP_MA=('keta'...
_F(FY=-80000.0,
GROUP_MA=('CW', ...
_F(FY=-60000.0,
GROUP_MA=('syuto...
そのため、AFFE_CHAR_MECA1つにつき,載荷条件は1つしか与えて...
そこで、FORCE_NODALEを用いて梁荷重を再現する方法を見つけ...
1.Meshのグループ作成で,梁荷重として載荷したい要素のノード...
2.AFFE_CHAR_MECA→FORCE_NODALEを選択→載荷値は(分布荷重×要...
例：モデルに死荷重を掛けたい。（桁長：300m,2m間隔で要素分...
(50000×300)/150 = 100000N~
この値を入力する。
**動的応答解析のやり方 [#t4e8361f]
・注意点~
(1)DEFI_FONCTIONとAFFE_CHAR_MECAは、Pre Analysis・Analysi...
(2)DISCRETに質量マトリクスを用意しておくこと。DISCRETには...
この設定をしないと、ASSEMBLAGEでバネ要素に質量マトリクス...
[[詳しくはこちらを参照:https://code-aster.org/forum2/view...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)DYNA_VIBRAとCALC_CHAMPの名前は同じにすること。(同じ名...
(4)CALC_CHAMPのEXCITにDYNA_LINEで設定したCHARGEの分だけ入...
設定方法~
[1]Pre analysisのASSEMBLAGEを選択~
・ここでやること~
①NUME_DDLは{}押し、nddlと入力。~
②MATR_ASSEで質量マトリクスと剛性マトリクスを定義する。(質...
③CHARGEはAFFE_CHAR_MECAで適用した境界条件分だけ追加する。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
[2]DYNA_NON_LINE~
ここで弾塑性を入れた動的解析ができると思われるが、そのあ...
・最低限やること
①SCHEMA TEMPSでNEWMARK,FORMULATIONをDEPLACEMENT(変位),β=0...
ちなみに,[[説明書:http://dugi6514.odns.fr/doc/v15/en/man_...
更に詳しく調べると,これは直接積分法（正確には直接時間積分...
（Marcもニューマークβ法を用いて連鎖崩壊に関する解析をして...
[[参考サイト:http://jikosoftcom.blog25.fc2.com/blog-entry...
②INCREMENTはDEFI_LIST_INSTで設定したものを選ぶ。~
③EXCITにAFFE_CHAR_MECAで作成した数だけ設定する。載荷条件...
ビーム要素なら、CARA_ELEMにチェックを入れること。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**陰解法と陽解法について [#qe1cdcdd]
陰解法は時間変化を取り扱う解析で、現在の時間TからT+ΔTに進...
特徴として...
-非線形性が強い問題では収束しないことがあるため、静的や準...
-各step(時刻)で連立方程式を解く必要があるため、解析に時間...
陽解法は現在の時間Tをもとに、時刻T+ΔTの解を求める解法.~
特徴として...
-非線形性が強い問題を扱いやすいため、衝突や落下などの短い...
-時間増分ΔT（timestep）を[[クーラン条件:https://www.cyber...
→どうやら非常に細かいtimestepにしないといけないため、こち...
Salome-Mecaでは,DYNA_NON_LINE→SCHEMA_TEMPで設定可~
陰解法：NEWMARK,HHT~
陽解法：DIFF_CENT,TCHAMWA~
***各解法の説明 [#k9ed93bb]
-陰解法
・NEWMARK~
時間積分法の一種.平均加速度法の係数β=0.25,γ=0.50の組み合...
・HHT~
①修正平均加速度:組み合わせ「HHT,MODI_EQUI = 'NOT'」~
②Hilber-Hughes-Taylor法：組み合わせ「HHT,MODI_EQUI = 'OUI...
-陽解法
・DIFF_CENT~
・TCHAMWA~
**ある時間で力を0に落とす方法 [#s666447f]
FEMにおいて,要素をある時間で切断・消すことはできない。（...
そこで,力を0にしたい要素をあらかじめ消しておいて,元々要素...
ケーブル2本で桁を吊っている簡易モデルを例に説明する。解析...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
＜下準備＞~
(1)まず、破断させたい要素の軸力(EFGE_ELNO)をデータ化して...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(2)110sまで画面を移動させ,110.001sで=0を入力。そして,解析...
画像では「,」ついてますが、この時点では「,」要らないので...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)時間データを違うセルに貼り付け、その隣に,=","&並べた張...
「,」を入れないと、Salomeにデータを入れることができないた...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
ちなみに、0.005s掛けて張力を0に落としたい場合、破断初めの...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(4)作成したデータをテキストエディタに貼り付け、それを保存...
＜Salomeでの設定方法＞
(1)Meshの画面に行き,破断させたい要素を取り除いたcompound_...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(2)DEFI_FONCTION→VALEを選択。import tableを選択.（表に矢...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)import tableを選択すると,ファイルを選択する画面が出る...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(4)BC and Load → FORCE_NODALEを選択.~
ここで,元々要素があった向きに軸力を与える必要がある。元々...
要素をなくした部分を斜辺として,1つの三角形として表現する...
シンプルな方法として,三平方の定理を用いて斜辺を求め,底辺/...
今回は1:1:√2の三角形のため,1/√2の力の分解となる。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
そして,(2)で与えたDEFI_FONCTIONで与えた関数を与えることで...
**温度でのプレストレスの与え方（編集中） [#vac22495]
ここでは温度でのやり方を示す。（ひずみは後述のPRE_EPSIに...
(0)DEFI_FONCTIONで温度変化のfunctionを作成する。（MECA_ST...
NOM_PARA（横軸）をINST,NOM_RESU（縦軸）をTEMPに設定。~
PROL_はVALEで作成したグラフの範囲外（DROITEは右側,GAUCHE...
写真はEXCLUを使っているが,これは除外という意味。（VALEの...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
VALEの設定方法はSTAT_NON_LINE(静的解析)orDYNA_NON_LINEと...
写真は静的解析を想定したもので,0倍の変化時に温度0、1倍の...
時刻歴の場合だと、解析時間0秒のときに温度0℃、解析時間1秒...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
ちなみに...張力から温度を逆算する場合は、Δt = T/EAα を用...
(1)Post Processing → CREA_CHAMPを選択。~
TYPE_CHAMは「NOEU_TEMP_F」を選択する。（MECA_STATICの場合...
OPERATIONは「AFFE」,MODELEにチェックを入れる。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
AFFEを1itemにする.温度荷重を与えたい部材をGROUP_MAで選択.~
VALEはMECA_STATICSを用いる場合,ここで最終的に変化させたい...
VALE_FはMECA_STATICS以外で用い,(0)で作成した温度のfunctio...
NOM_CMPはTEMPと入力.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(2)Output → CREA_RESUを選択。（MECA_STATICSの場合は必要な...
この設定をしないと温度荷重を考慮して解析を回してくれない...
(3)AFFE_MATERIAU → AFFE_VARCを選択。~
**ひずみでのプレストレスの与え方[#c71fb860]
（追記）~
こちらは、プレストレスをひずみで与えるやり方です.~
簡易モデルでは、この方法でプレストレスを与えたが、PRE_EPS...
PRE EPSIは「設定したひずみになるように外力を与える」コマ...
それを回避するために,要素分割し,分割した1要素に以下の式で...
PRE_EPSIをプレストレスとして用いる際の式：ε = α T × 要素...
[[PRE_EPSIの説明書:http://dugi6514.odns.fr/doc/v15/en/man...
モデルをひずみを用いて変形させることができるコマンドであ...
1Dモデル(Multifiber要素も含む)ではEPX（梁のローカル軸に応...
(ky,kzはたわみの微分方程式から得られるdθ/dx = -M/EI の値...
X,Y,Zの軸の向きはグローバル座標(salome画面の色がついてる...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
単純梁モデル(10✕10✕120mm)にEPXX：0.01を与えると、以下の挙...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
線に色がついてるところが変形前のモデルの端であり、写真を...
・簡易モデルに用いた際に分かったこと~
現在、ケーブルにPRE_EPSIを入れており、縮む挙動(ケーブルが...
その理由について考えてみたが、要素を構成している節点にモ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑PRE_EPSIの意味 日本語に訳すと要素に変形の負荷を与えるみ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑簡単な例 PRE_EPSIで与えたひずみになるような外力をかけて...
**CALC_CHAMP [#t08db3b9]
***POU_D_E,POU_D_EMで計算できるもの（作成中） [#n061551b]
**IMPR_RESU [#m5434190]
***IMPR_RESUのGROUP_MA,GROUP_NO [#wafb540f]
IMPR_RESUのRESUの中に入っているGROUP_MA,GROUP_NOを使うこ...
これを使うメリットは容量の大きいモデルで解析を回す際に,見...
*影響線解析 [#h9310a98]
影響線：移動荷重に対してある点に生じる物理量の大きさを移...
影響線の定義から,まずは移動荷重を求め方が分からないと話に...
↓橋の上を走る電車を例に移動荷重を算出する方法を見つけた。...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
[[移動荷重について:https://forum.code-aster.org/public/d/...
[[移動荷重について_その2:https://forum.code-aster.org/pub...
AFFE_CHAR_MECA_F→FORCE_POUTREを使い、DYNA_NON_LINEなどの...
ちなみに...midasだと[[説明書の10章:https://in.midasit.co....
***軸力の見方 [#df51b2bf]
軸力はEFGE_ELNOの[N]で見る事ができることが分かった。Nは法...
見る方法だが、medファイルをParaViSで閲覧し、フィルター→me...
[[参考資料:http://dugi6514.odns.fr/doc/v15/en/man_u/u2/u2...
↑Salomeで出力される力やひずみなどについて詳しく書かれてい...
**断面定数の算出方法 [#g12212f4]
複雑な断面をビームモデルに適用したい場合,異なるhdfファイ...
ここでは,そのやり方を説明します.
・（橋梁の主塔のような内側に穴が空いている矩形断面の作成...
高さ2m,幅4m,厚さ0.025mの矩形断面を例に説明します。
○Geometry~
・平面を用いて断面を作る~
長方形のフェースを作成をクリック.ここで,高さor幅＋厚さの...
（big：高さ2m+厚さ0.025m,幅4m+0.025m, small：高さ2m−厚さ0...
断面を作成する際,座標の向きはX座標下向き,Y座標右向きで作...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
作成したら,Cutを用いて,厚さ0.025mの矩形断面を作成する.2種...
ここでは,bigをメインオブジェクト,smallをツールオブジェク...
後,下の写真では高度のオプションに2種類チェックを入れてし...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
Cutすると,以下の写真のようになる.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・外側の四辺、内側の四辺、外側の四辺上の角の点、内側の四...
外側の点と内側の点,外側の線と内側の線,計16個全てをグルー...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↓こんな感じになるかと.グループの名前は分かりやすいものに...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
○Mesh~
・メッシュを作成~
アルゴリズムはNETGEN1D-2D,詳細設定は歯車押して,適用すれば...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・ジオメトリグループの作成~
点はノード,線は要素に入れる.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・グループを作成~
Faceにチェックを入れ,全て選択します。を選ぶ.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・メッシュをmedファイルでエクスポート~
作成したMeshを右クリックして,エクスポート→MEDファイルを選...
○aster_study~
・LIRE_MAILAGEでメッシュ適用~
エクスポートしたMEDファイルを適用する.やり方は,…マークを...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・Model Definition > MACR_CARA_POUTRE >「GROUP_MA_BORD」...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
MACR_CARA_POUTREは1つの辺に対し,1つの点しか適用できないた...
面倒だと思うけど...一つずつ設定してください.~
また,設定の仕方ですが,一筆書きになるように作成してくださ...
そして,外側の下辺(lsoto3),左辺(lsoto4)と設定すると一筆書...
なお,「GROUP_NO」については黄色の線に対応するもの（今の状...
画像の例だと,soto1→soto2→soto3→soto4になります.（上辺の左...
↓「GROUP_MA_BORD」の設定画面~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↓「GROUP_NO」の設定画面~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・Model Definition > MACR_CARA_POUTRE にて内側も同様に~
外側の設定と同じようにやればOK.
・IMPR_RESU で FormatをRESULTATに、ファイル名を「kekka.re...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
○計算を回してresultファイルを見る。~
・ファイルに書かれているA(断面積),IY,IZ(断面2次モーメント...
なお、空白エリアが存在しない断面（密な断面）を作成する際...
また、.resuファイルはhdfファイルやFilesと同じところに保存...
○メインとなるhdfファイル（橋梁モデル等,主となるモデルがあ...
・断面形状の入れ方~
LIRE_MAILAGE→...をクリック→エクスポートした断面形状のmed...
*材料非線形 [#ad9f0c30]
**梁要素に適用できる（できそう？）弾塑性モデル [#y6458c84]
[[Astercode_U4.51.11:https://code-aster.org/doc/v15/en/ma...
[[Astercode_R5.03.09:https://code-aster.org/doc/v15/en/ma...
梁要素に適用できる,できそうな弾塑性モデルをまとめた。~
1Dモデルに対応している弾塑性モデルは以下のとおりである。~
-VMIS_ISOT_LINE 対称線形等方加工硬化を伴うフォンミーゼス
-VMIS_ISOT_TRAC 非特定の等方性加工硬化を伴うフォンミーゼス
-VMIS_CINE_LINE 対称線形キネマティック加工硬化を伴うフォ...
-ECRO_CINE_1D 対称線形キネマティック加工硬化を伴うフォン...
-GRILLE_ISOT_LINE 対称線形等方性加工硬化を伴うフォンミー...
-GRILLE_CINE_LINE 対称線形キネマティック加工硬化を伴うフ...
-PINTO_MENEGOTTO 鉄筋コンクリート鉄筋の挙動
-GRILLE_PINTO_MEN 鉄筋コンクリート配筋の挙動
-VMIS_ASYM_LINE 非対称線形加工硬化および復元によるフォン...
-VISC_IRRA_LOG,GRAN_IRRA_LOG 燃料集合体の粘塑性挙動： REF...
-MAZARS MAZARSのバージョン1Dの挙動.
-RELAX_ACIER プレストレストケーブルのリリーフをモデル化す...
説明は皆さんが使うことがある、可能性があるものを以下に記...
-VMIS_ISOT_LINE,VMIS_CINE_LINE~
適用可：BAR,POU_D_EM,POU_D_TGM~
フォンミーゼスの弾塑性挙動と線形加工硬化の関係。~
材料に必要なデータは、オペレータDEFI_MATERIAU [U4.43.01]...
VMIS_ISOT_LINEは大変形SIMO_MIEHEに対応しているとのこと。~
-VMIS_ISOT_TRAC~
適用可：BAR,POU_D_EM,POU_D_TGM~
フォンミーゼスの弾塑性挙動と非線形等方加工硬化の関係。~
曲線σ,εは、DEFI_MATERIAU [U4.43.01]のTRACTIONで提供される...
温度に応じていくつかのトラクションダイアグラムを定義する...
また、DEFI_MATERIAU演算子のELAS (_FO)は必ず必要。~
トラクションダイアグラムが提供された場合、挙動の関係に使...
大変形SIMO_MIEHEに対応しているとのこと。~
-VMIS_ASYM_LINE~
適用可：BAR~
説明：等方性加工硬化を伴うフォンミーゼスの弾塑性一軸等温...
制御装置や埋設ケーブルと地盤との相互作用をモデル化するこ...
材料に必要なデータは、DEFI_MATERIAU [U4.43.01]のECRO_ASYM...
*モーダル解析 [#r5d7cbd2]
モーダル解析とは~
対象物そのものが振動しやすい周波数とその振動形状を把握す...
構造物の固有振動数と固有モードが見れるとのこと.(よく皆が...
固有振動数とは？何故、固有振動数を求めるのか？~
→物体ごとがもつ周波数(1秒間に揺れる振動回数)、構造物の共...
固有モードとは？~
→それぞれの固有振動数においてどのように振動しているかの変...
片持梁を想定したモーダル解析コマンドファイル~
PreanalysisのASSEMBLAGEとAnalysisのCALC_MODEを使うとのこ...
mesh = LIRE_MAILLAGE(identifier='0:1',
FORMAT='MED',
UNITE=20)
mesh0 = LIRE_MAILLAGE(identifier='1:1',
FORMAT='MED',
UNITE=3)
model = AFFE_MODELE(identifier='2:1',
AFFE=_F(GROUP_MA=('beam', ),
MODELISATION='POU_D_EM',
PHENOMENE='MECANIQUE'),
MAILLAGE=mesh)
fiber = DEFI_GEOM_FIBRE(identifier='3:1',
SECTION=_F(COOR_AXE_POUTRE=(0.0, ...
GROUP_FIBRE='men',
GROUP_MA_SECT=('danmen...
MAILLAGE_SECT=mesh0))
elemprop = AFFE_CARA_ELEM(identifier='4:1',
GEOM_FIBRE=fiber,
MODELE=model,
MULTIFIBRE=_F(GROUP_FIBRE=('men...
GROUP_MA=('beam',...
PREC_AIRE=0.76,
PREC_INERTIE=0.95),
POUTRE=_F(CARA=('HY', 'HZ'),
GROUP_MA=('beam', ),
SECTION='RECTANGLE',
VALE=(10.0, 20.0)))
mater = DEFI_MATERIAU(identifier='5:1',
ELAS=_F(E=206000.0,
NU=0.3,
RHO=7.9e-09))
load = AFFE_CHAR_MECA(identifier='6:1',
DDL_IMPO=_F(DRX=0.0,
DRY=0.0,
DRZ=0.0,
DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_NO=('kotei', )),
MODELE=model)
behav = DEFI_COMPOR(identifier='7:1',
GEOM_FIBRE=fiber,
MATER_SECT=mater,
MULTIFIBRE=_F(GROUP_FIBRE=('men', ),
MATER=mater,
RELATION='ELAS'))
fieldmat = AFFE_MATERIAU(identifier='8:1',
AFFE=_F(GROUP_MA=('beam', ),
MATER=(mater, )),
AFFE_COMPOR=_F(COMPOR=behav,
GROUP_MA=('beam',...
MODELE=model)
ASSEMBLAGE(identifier='9:1',
CARA_ELEM=elemprop,
CHAM_MATER=fieldmat,
CHARGE=(load, ),
MATR_ASSE=(_F(MATRICE=CO('mass'),
OPTION='MASS_MECA'),
_F(MATRICE=CO('rigi'),
OPTION='RIGI_MECA')),
MODELE=model,
NUME_DDL=CO('nddl'))
modes = CALC_MODES(identifier='12:1',
CALC_FREQ=_F(NMAX_FREQ=10),
MATR_MASS=mass,
MATR_RIGI=rigi,
OPTION='PLUS_PETITE',
VERI_MODE=_F(STOP_ERREUR='OUI'))
IMPR_RESU(identifier='13:1',
FORMAT='MED',
RESU=_F(CARA_ELEM=elemprop,
NOM_CHAM=('DEPL', ),
RESULTAT=modes),
UNITE=4)
非線形性を考慮するには,DYNA_NON_LINEを使うらしい。要調査
結果比較~
,,鉛直1次モード,水平1次モード,
,結果,7.96Hz,16.35Hz,
,理論値,8.22Hz,16.45Hz,
理論値と近い値が出たため、ビーム要素におけるモーダル解析...
*エラーメッセージ一覧 [#m97c0c8e]
・このエラーメッセージは構造的に不安定なときに出る。~
これが出た場合は基本的にAsterstudyの設定を間違えているこ...
それでもダメなら,ファイバー要素を適用している断面のメッシ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・AFFE_CHAR_MECAの設定をモデルに適用できない場合に出る。...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・AFFE_MODELEとAFFE_CARA_ELEMの条件が一致しない場合に出る...
AFFE_MODELEで設定したMODELISATIONは、それぞれAFFE_CARA_EL...
(例) AFFE_MODELEでBARREを選択した場合→AFFE_CARA_ELEMではB...
AFFE_MODELEでPOU_D_Eを選択した場合→AFFE_CARA_ELEMで...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**収束エラーが出た場合の対処法（マルチファイバー要素を使...
収束エラーや行列絡みのエラーが出たときは、まずtimestepや...
ファイバー要素で作成した片持ち梁モデルで断面メッシュの細...
解析内容は片持ち梁の静的材料非線形解析である。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
***三角形メッシュの場合 [#m07ec30d]
載荷条件だが、与えた条件より値が小さくてもモデルが降伏し...
変えたのは断面の要素数のみで他の設定は変えていない。~
この条件では,要素数156では行列絡みのエラーが出たため解析...
ちなみに,解析結果だが,どの部分も設定した降伏応力を超える...
また,固定面を見ると,綺麗に引張,圧縮と分かれている.(どちら...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
解析結果を見て,分かったことがあったので、ここにまとめる。~
・ビーム要素(線要素)の解析において,一気に変形する挙動を起...
以下のアニメーションを見ると分かりやすいと思うが,大きな変...
実は,ちょうど解析エラーで止まった場所はこの大きな変化が起...
もしかすると,材料非線形を考慮した解析でエラーが出てしまう...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/ex1.gif
***四角形メッシュの場合 [#e01f2635]
しかし、三角形メッシュだと要素数が増やさないと解析が回ら...
そこで、三角形メッシュと同じ解析条件のまま、メッシュの形...
四角形メッシュの変え方はNETGEN1D-2Dを選択し,メッシュサイ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
四角形メッシュの要素数：100~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
この状態でも解析は回った。三角形メッシュでは解析が回らな...
応力をみると、こちらも設定した降伏応力を超えることはなか...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
***結果比較とまとめ [#g5773cf3]
とりあえず、先端変位と固定面の応力ひずみを確認した。
最終変位~
三角形:6.06mm 四角形：2.31mm~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
理論値は1.17mmで,結果を見ると,三角形メッシュの場合は1.30m...
グラフを見る限り,荷重が掛かり切ると,梁が大きくたわむこと...
応力ひずみ~
三角形メッシュの応力ひずみは以下の通り。中立軸データはグ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
四角形メッシュの応力ひずみは以下の通り。（要素数が少ない...
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/experi...
グラフを見る限り、綺麗な三角形応力分布ができていると思わ...
まとめ~
どちらのメッシュの切り方でも応力ひずみ関係はほぼ同じ挙動...
これは,恐らく要素数の違いが原因だろう.なので,四角形メッシ...
ただ,たわみの理論値を少し超えた付近で線形関係から外れるた...
しかし、橋梁モデルのような大きいモデルを解析する際に三角...
以上から、大きめなモデルを解析する際は,要素数少なめで解析...
現在の閲覧人数: &online;~
総閲覧数:&counter;~
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
