Marc/Mentatメモ
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
Marc/Mentatの操作方法等をメモしていく。
*情報統括センターでの起動方法 [#cc89e727]
windowsを検索でmarcと入力すると出てくる。情報統括センター...
-バージョン2016.0.0(2018/4/27現在)
*基本操作（バージョン2016.0.0） [#r6db65f7]
-平行移動（左クリックしながら）
-拡大・縮小（右クリックしながら）
-回転（ホイール押しながら）
*ブーリアン演算 [#l7de5c30]
「形状とメッシュ」 -> 「形状とメッシュ」 -> オペレーショ...
-結合
-取り去る
-交差
*円柱の作り方 [#f5983d80]
「形状とメッシュ」 -> 「形状とメッシュ」 -> ソリッドのプ...
-この通り入力するとZ座標に100で半径8の円柱が完成する。
*弾塑性 [#t01d475b]
-http://www.mscsoftware.com/training_videos/patran/Reverb...
**降伏後の傾き（完全弾塑性とかバイリニアとか）の設定（未...
PLASTICITY PROPERTIESのYIELD CRITERIONはVON MISESとして、
METHODでPOWER LAWを選ぶと、COEFFICIENT AとBが入力できるよ...
この辺に適当な数値を入れて、降伏後の傾きがどう変わるか調...
それでうまくいかなければ、この辺のオプションのキーワード
POWER LAWとかPIECEWISE LINEARとかとmarc mentat elasto-pla...
*座屈 [#l8a21d76]
-いったんSTATICで解析した後にBUCKLEで再び解析するとできる...
-SHELL要素で座屈荷重が得られなかったがANALYSIS OPTIONSを...
**倍率変更 [#ya31db49]
DEFORMATION SCALINGをMANUALにして自分で値を入力する
**座屈モード [#sa4b4845]
-BUCKLE MODESを変更すれば良い。初期は2となっているが5を入...
**梁で解く [#u58e3da9]
-手計算と同じ値が得られた。
***中身 [#x80a3ce8]
-GEOMTRIC
--TYPE
---PLANAR
---STRAIGHT BEAM
-MATERIAL
--GENERAL
---いじらなくても同じ結果が得られた
-JOBS
--INTIAL LOADS
---選択する。載荷方向を知るために必要かな。
--ANALYSIS OPTIONS
---LARGE STRAIN選択
---BUCKLE INCREMENTSをONにする。
---両方ともいじらなくても同じ結果
---PLANE STRAIN & PLANE STRESS & 2D のどれでもできた。
***座屈荷重を見る [#m58bc2b6]
-RESULTS画面のFac:○○○ ←これが手計算と一致
-その時のInc: 0:1 だった
-mentatもcalculix同様(*荷重)が座屈荷重かな
*終了番号の説明を載せる [#r5dcf0ff]
***13 [#q543e9e8]
-入力データにエラーがある。
***2004 [#p7364641]
-剛体変位が発生している。
-または全体剛性マトリクスが非正定マトリクスになっている．
-境界条件が正しいか確認する。
The determinant of the stiffness matrix becomes zero or
negative when indicated node has been reached during the
Gaussian elimination phase of the solution process. This
means that the stiffness matrix is non-positive definite.
If this happens at the start of the analysis, the condition
is usually caused by the existence of rigid body modes. It
may also be caused by incorrect material properties (e.g.
Poisson's ratio greater than 0.5; note that such situations
may arise through temperature dependence of properties). In
nonlinear cases, the structure may have buckled or reache...
plastic limit load. In rubber analysis, it may also be due
to the strain state being in a region where the input data
for the strain energy function is invalid. In contact
analysis with friction, lack of normal forces may result in
friction being absent. If desired, the program may be for...
to continue by use of the parameter PRINT or the model or
history definition option CONTROL. Either one of these
procedures may be used for restart. Whenever a non-positive
definite situation occurs one must exercise caution, as the
resultant numerical solution may be infeasible.
***3002 [#hffacc1f]
-収束公差内に収束しなかったときに出る。
Analysis has failed to converge to required convergence...
***3020 [#a2c152b6]
-ユーザーの介入で止めたときに出るエラー。こちらから何か止...
***3300 [#f1e0f70a]
-固有値抽出が最大許容反復回数内で収束できない。
-BUCKLEオプションで反復回数か収束判定許容差を大きくする。
***3301 [#u3bcb7a5]
-質量か初期応力剛性がゼロになっている。
-座屈解析では構造に荷重がすでに適用され、応力状態が存在し...
***3302 [#f57e99b8]
-解析においてランチョス法による固有値抽出の途中、非正定値...
-解析実行前に作用荷重を減らすか、BUCKLEを用いて逆べき乗法...
***3009 [#uf4b8e18]
-タイムステップ数が小さいので続行できないときに出るエラー
The time step size has become too small to continue the
analysis due to too many time step cut-backs in the current
increment.
It is not possible to obtain convergence with the specified
allowed maximum number of cut-backs.
***9999 [#y81a229b]
Marc内部でエラーが起きた場合に出る．お問い合わせしてと言...
*Mentatの画面が透明の時の対処方 [#ia09f71d]
デスクトップのシステム-設定-外観-視覚効果で効果なしを選ぶ...
*各種折り紙構造のstlファイル(ハニカムはまだわかんない) [#...
-S6のinpファイルをつくる。
-[[ここ:http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wik...
-それから[[こちら:http://www.str.ce.akita-u.ac.jp/cgi-bin...
-inpファイルを[[scp:http://www.str.ce.akita-u.ac.jp/cgi-b...
--(
$ssh -X hoge@hage.jp
--)
みたいにして、Marcサーバーにログイン-->
--(
$mentat2011
--)
でMarc起動。
-下段の「FILES」-->ウインドウ左の「IMPORT」-->「ABAQUS」-...
-できたstlをscpで持ってきて、FreeCADで開いて中身確認。
でもこのstlファイルはObjet Studioだｔ$\frac{1}{1000}$mmで...
*s6要素を要素数を4倍にしてs3要素に変換するプログラム [#yd...
-[[s6-s3.f90:https://www.str.ce.akita-u.ac.jp/~gotouhan/t...
*kisoken@へのファイルを送る方法とlocalに持ってくる方法 [#...
**kisokenに送る [#s5424aff]
-送りたいファイルがあるディレクトリに入る
-scp 送りたいファイル名 kisoken@quartet.gipc.akita-u.ac.j...
-パスワードを聞かれるので入力
**kisokenから取ってくる [#k498ede8]
-コピーしたいディレクトリに入る
-scp kisoken@quartet.gipc.akita-u.ac.jp:取ってきたいファ...
*今までもダイヤカットのinpファイルの読み方 [#id68e436]
(*SHELL SECTION,ELSET=Eall,MATERIAL=EL,OFFSET=0.)という部...
(*SHELL SECTION,ELSET=Eall,MATERIAL=EL)に変更
(*ELEMENT, TYPE=S6, ELSET=Eall)の部分を
(*ELEMENT, TYPE=S3R, ELSET=Eall)に変更すると読み込める
OFFSETはCOMPOSITEなしでは対応していないようだ
-上みたいに直しても読めないのもあった。細かいダイヤだと読...
*読み込んでからの要素タイプの変更 [#j698fdaa]
Main→Mesh generation → element type
*mentatへのインポートの仕方 [#s636ae2b]
FILES→IMPORT→ABAQUS等が選べる
*インポートしてからの手順 [#qc68e4e7]
*異方性材料を使ってる場合 [#j09347b7]
*CalculiXのinpファイルをインポートできるか？ [#rf3063ce]
-http://tech.groups.yahoo.com/group/marc-users/message/994j
*inpファイルを読むときの注意 [#t395bef1]
-inpファイルの中にd-3とかがあるとよめないので下のように変...
206d3→206E+003
*インポートできるものとできないの [#u639afa6]
-inpファイルは読み込める
-s6,s8要素は対応していない
Warning: Element Type S6 is not supported in this version
-S8R要素は読めるようだがcomposite optionがないとダメ
Error: MARC doesn't support offset of shell section witho...
-S3Rにすると読めるっぽい
-salomeのstlファイルも読み込める
-abaqus,acis,c-mold,dxf/dwg,i-deas,iges,marcinput,nastran...
*モデルの読み込み [#mc3f7314]
IGES,STEPファイルなら可能？
*直行異方性材料 [#ud495345]
nu12, nu23, nu31がどのポアソン比に対応しているかが不明。
詳細については、
[[構造メモ:https://www.str.ce.akita-u.ac.jp/~gotou/kouzou...
に書いた。
直行異方性の行列の表記は色んな流儀があるが、ここではまず...
以下のように書いてみる。
$$\left(\begin{array}{c} \varepsilon_{x}\\ \varepsilon_{y...
対称性が成り立つなら、以下も同じこと（こっちが主流かも）
$$\left(\begin{array}{c} \varepsilon_{x}\\ \varepsilon_{y...
-http://en.wikipedia.org/wiki/Orthotropic_material
-http://www.efunda.com/formulae/solid_mechanics/mat_mecha...
$x$軸方向に長い板を$x$軸方向に引っ張ってみる。
$\sigma_{y}=\sigma_{z}=0$とみなせるので、これを上の2式に...
$$\nu_{yx}=-\frac{E_{y}\varepsilon_{y}}{E_{x}\varepsilon_...
$$\nu_{zx}=-\frac{E_{z}\varepsilon_{z}}{E_{x}\varepsilon_...
同様に$y$軸方向に長い板を$y$軸方向に引っ張ってみる。
$\sigma_{x}=\sigma_{z}=0$とみなせるので、これを上の2式に...
$$\nu_{zy}=-\frac{E_{z}\varepsilon_{z}}{E_{y}\varepsilon_...
*例題 [#nf129d65]
$E_{x}=1$GPa, $E_{y}=10$GPa, $E_{z}=100$GPa,
として、対称性を満たす
$\nu_{xy}=0.03, \nu_{xz}=0.003, \nu_{yz}=0.03$,
$\nu_{yx}=\nu_{zx}=\nu_{zy}=0.3$
のポアソン比の材料で
$x$軸方向に長い板を直方体要素とかでモデル化して$x$軸方向...
（ひずみを直接出力するオプションがあるのかどうか知らない...
想定どおりのポアソン比が得られるかどうかをチェックする。
素直に
nu12$=\nu_{xy}=$0.03, nu23$=\nu_{yz}=0.003$, nu31$=\nu_{z...
とした場合はどうか（たぶん、これでいいような気が）。
nu12,nu23,nu31の組み合わせを変えてみてチェック。
[確認]
x=10 Y=1 Z=100の直方体要素をZ軸方向に引っ張る。
mentatにn12=0.016 n23=0.016 n31=0.4 Ex=0.28GPa Ey=0.28GPa...
解析後のそれぞれの軸方向の変化量からポアソン比を手計算に...
nzy=約0.4 nzx=約0.4 となった。直交異方性材料のポアソン比...
25n13=n31 n32=25n23 という関係が成り立つのでmentatに入力...
終了行:
Marc/Mentatの操作方法等をメモしていく。
*情報統括センターでの起動方法 [#cc89e727]
windowsを検索でmarcと入力すると出てくる。情報統括センター...
-バージョン2016.0.0(2018/4/27現在)
*基本操作（バージョン2016.0.0） [#r6db65f7]
-平行移動（左クリックしながら）
-拡大・縮小（右クリックしながら）
-回転（ホイール押しながら）
*ブーリアン演算 [#l7de5c30]
「形状とメッシュ」 -> 「形状とメッシュ」 -> オペレーショ...
-結合
-取り去る
-交差
*円柱の作り方 [#f5983d80]
「形状とメッシュ」 -> 「形状とメッシュ」 -> ソリッドのプ...
-この通り入力するとZ座標に100で半径8の円柱が完成する。
*弾塑性 [#t01d475b]
-http://www.mscsoftware.com/training_videos/patran/Reverb...
**降伏後の傾き（完全弾塑性とかバイリニアとか）の設定（未...
PLASTICITY PROPERTIESのYIELD CRITERIONはVON MISESとして、
METHODでPOWER LAWを選ぶと、COEFFICIENT AとBが入力できるよ...
この辺に適当な数値を入れて、降伏後の傾きがどう変わるか調...
それでうまくいかなければ、この辺のオプションのキーワード
POWER LAWとかPIECEWISE LINEARとかとmarc mentat elasto-pla...
*座屈 [#l8a21d76]
-いったんSTATICで解析した後にBUCKLEで再び解析するとできる...
-SHELL要素で座屈荷重が得られなかったがANALYSIS OPTIONSを...
**倍率変更 [#ya31db49]
DEFORMATION SCALINGをMANUALにして自分で値を入力する
**座屈モード [#sa4b4845]
-BUCKLE MODESを変更すれば良い。初期は2となっているが5を入...
**梁で解く [#u58e3da9]
-手計算と同じ値が得られた。
***中身 [#x80a3ce8]
-GEOMTRIC
--TYPE
---PLANAR
---STRAIGHT BEAM
-MATERIAL
--GENERAL
---いじらなくても同じ結果が得られた
-JOBS
--INTIAL LOADS
---選択する。載荷方向を知るために必要かな。
--ANALYSIS OPTIONS
---LARGE STRAIN選択
---BUCKLE INCREMENTSをONにする。
---両方ともいじらなくても同じ結果
---PLANE STRAIN & PLANE STRESS & 2D のどれでもできた。
***座屈荷重を見る [#m58bc2b6]
-RESULTS画面のFac:○○○ ←これが手計算と一致
-その時のInc: 0:1 だった
-mentatもcalculix同様(*荷重)が座屈荷重かな
*終了番号の説明を載せる [#r5dcf0ff]
***13 [#q543e9e8]
-入力データにエラーがある。
***2004 [#p7364641]
-剛体変位が発生している。
-または全体剛性マトリクスが非正定マトリクスになっている．
-境界条件が正しいか確認する。
The determinant of the stiffness matrix becomes zero or
negative when indicated node has been reached during the
Gaussian elimination phase of the solution process. This
means that the stiffness matrix is non-positive definite.
If this happens at the start of the analysis, the condition
is usually caused by the existence of rigid body modes. It
may also be caused by incorrect material properties (e.g.
Poisson's ratio greater than 0.5; note that such situations
may arise through temperature dependence of properties). In
nonlinear cases, the structure may have buckled or reache...
plastic limit load. In rubber analysis, it may also be due
to the strain state being in a region where the input data
for the strain energy function is invalid. In contact
analysis with friction, lack of normal forces may result in
friction being absent. If desired, the program may be for...
to continue by use of the parameter PRINT or the model or
history definition option CONTROL. Either one of these
procedures may be used for restart. Whenever a non-positive
definite situation occurs one must exercise caution, as the
resultant numerical solution may be infeasible.
***3002 [#hffacc1f]
-収束公差内に収束しなかったときに出る。
Analysis has failed to converge to required convergence...
***3020 [#a2c152b6]
-ユーザーの介入で止めたときに出るエラー。こちらから何か止...
***3300 [#f1e0f70a]
-固有値抽出が最大許容反復回数内で収束できない。
-BUCKLEオプションで反復回数か収束判定許容差を大きくする。
***3301 [#u3bcb7a5]
-質量か初期応力剛性がゼロになっている。
-座屈解析では構造に荷重がすでに適用され、応力状態が存在し...
***3302 [#f57e99b8]
-解析においてランチョス法による固有値抽出の途中、非正定値...
-解析実行前に作用荷重を減らすか、BUCKLEを用いて逆べき乗法...
***3009 [#uf4b8e18]
-タイムステップ数が小さいので続行できないときに出るエラー
The time step size has become too small to continue the
analysis due to too many time step cut-backs in the current
increment.
It is not possible to obtain convergence with the specified
allowed maximum number of cut-backs.
***9999 [#y81a229b]
Marc内部でエラーが起きた場合に出る．お問い合わせしてと言...
*Mentatの画面が透明の時の対処方 [#ia09f71d]
デスクトップのシステム-設定-外観-視覚効果で効果なしを選ぶ...
*各種折り紙構造のstlファイル(ハニカムはまだわかんない) [#...
-S6のinpファイルをつくる。
-[[ここ:http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wik...
-それから[[こちら:http://www.str.ce.akita-u.ac.jp/cgi-bin...
-inpファイルを[[scp:http://www.str.ce.akita-u.ac.jp/cgi-b...
--(
$ssh -X hoge@hage.jp
--)
みたいにして、Marcサーバーにログイン-->
--(
$mentat2011
--)
でMarc起動。
-下段の「FILES」-->ウインドウ左の「IMPORT」-->「ABAQUS」-...
-できたstlをscpで持ってきて、FreeCADで開いて中身確認。
でもこのstlファイルはObjet Studioだｔ$\frac{1}{1000}$mmで...
*s6要素を要素数を4倍にしてs3要素に変換するプログラム [#yd...
-[[s6-s3.f90:https://www.str.ce.akita-u.ac.jp/~gotouhan/t...
*kisoken@へのファイルを送る方法とlocalに持ってくる方法 [#...
**kisokenに送る [#s5424aff]
-送りたいファイルがあるディレクトリに入る
-scp 送りたいファイル名 kisoken@quartet.gipc.akita-u.ac.j...
-パスワードを聞かれるので入力
**kisokenから取ってくる [#k498ede8]
-コピーしたいディレクトリに入る
-scp kisoken@quartet.gipc.akita-u.ac.jp:取ってきたいファ...
*今までもダイヤカットのinpファイルの読み方 [#id68e436]
(*SHELL SECTION,ELSET=Eall,MATERIAL=EL,OFFSET=0.)という部...
(*SHELL SECTION,ELSET=Eall,MATERIAL=EL)に変更
(*ELEMENT, TYPE=S6, ELSET=Eall)の部分を
(*ELEMENT, TYPE=S3R, ELSET=Eall)に変更すると読み込める
OFFSETはCOMPOSITEなしでは対応していないようだ
-上みたいに直しても読めないのもあった。細かいダイヤだと読...
*読み込んでからの要素タイプの変更 [#j698fdaa]
Main→Mesh generation → element type
*mentatへのインポートの仕方 [#s636ae2b]
FILES→IMPORT→ABAQUS等が選べる
*インポートしてからの手順 [#qc68e4e7]
*異方性材料を使ってる場合 [#j09347b7]
*CalculiXのinpファイルをインポートできるか？ [#rf3063ce]
-http://tech.groups.yahoo.com/group/marc-users/message/994j
*inpファイルを読むときの注意 [#t395bef1]
-inpファイルの中にd-3とかがあるとよめないので下のように変...
206d3→206E+003
*インポートできるものとできないの [#u639afa6]
-inpファイルは読み込める
-s6,s8要素は対応していない
Warning: Element Type S6 is not supported in this version
-S8R要素は読めるようだがcomposite optionがないとダメ
Error: MARC doesn't support offset of shell section witho...
-S3Rにすると読めるっぽい
-salomeのstlファイルも読み込める
-abaqus,acis,c-mold,dxf/dwg,i-deas,iges,marcinput,nastran...
*モデルの読み込み [#mc3f7314]
IGES,STEPファイルなら可能？
*直行異方性材料 [#ud495345]
nu12, nu23, nu31がどのポアソン比に対応しているかが不明。
詳細については、
[[構造メモ:https://www.str.ce.akita-u.ac.jp/~gotou/kouzou...
に書いた。
直行異方性の行列の表記は色んな流儀があるが、ここではまず...
以下のように書いてみる。
$$\left(\begin{array}{c} \varepsilon_{x}\\ \varepsilon_{y...
対称性が成り立つなら、以下も同じこと（こっちが主流かも）
$$\left(\begin{array}{c} \varepsilon_{x}\\ \varepsilon_{y...
-http://en.wikipedia.org/wiki/Orthotropic_material
-http://www.efunda.com/formulae/solid_mechanics/mat_mecha...
$x$軸方向に長い板を$x$軸方向に引っ張ってみる。
$\sigma_{y}=\sigma_{z}=0$とみなせるので、これを上の2式に...
$$\nu_{yx}=-\frac{E_{y}\varepsilon_{y}}{E_{x}\varepsilon_...
$$\nu_{zx}=-\frac{E_{z}\varepsilon_{z}}{E_{x}\varepsilon_...
同様に$y$軸方向に長い板を$y$軸方向に引っ張ってみる。
$\sigma_{x}=\sigma_{z}=0$とみなせるので、これを上の2式に...
$$\nu_{zy}=-\frac{E_{z}\varepsilon_{z}}{E_{y}\varepsilon_...
*例題 [#nf129d65]
$E_{x}=1$GPa, $E_{y}=10$GPa, $E_{z}=100$GPa,
として、対称性を満たす
$\nu_{xy}=0.03, \nu_{xz}=0.003, \nu_{yz}=0.03$,
$\nu_{yx}=\nu_{zx}=\nu_{zy}=0.3$
のポアソン比の材料で
$x$軸方向に長い板を直方体要素とかでモデル化して$x$軸方向...
（ひずみを直接出力するオプションがあるのかどうか知らない...
想定どおりのポアソン比が得られるかどうかをチェックする。
素直に
nu12$=\nu_{xy}=$0.03, nu23$=\nu_{yz}=0.003$, nu31$=\nu_{z...
とした場合はどうか（たぶん、これでいいような気が）。
nu12,nu23,nu31の組み合わせを変えてみてチェック。
[確認]
x=10 Y=1 Z=100の直方体要素をZ軸方向に引っ張る。
mentatにn12=0.016 n23=0.016 n31=0.4 Ex=0.28GPa Ey=0.28GPa...
解析後のそれぞれの軸方向の変化量からポアソン比を手計算に...
nzy=約0.4 nzx=約0.4 となった。直交異方性材料のポアソン比...
25n13=n31 n32=25n23 という関係が成り立つのでmentatに入力...
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
