#author("2022-10-07T15:12:43+09:00","default:kouzouken","kouzouken")
#contents
*Marcマニュアルの取得方法 [#p5d8b924]
-Marcを起動 --> ヘルプ(メニュータブ)--> でMarcのマニュアル(英語)を取得できる.~

・User's Guide~
・VolumeA: Theory and User Information~
・VolumeB: Element Library~
・VolumeC: Program Imput~
・VolumeD: User Subroutines~
・VolumeE: Demonstration Problems~
・検証ガイド~
・Pythonマニュアル~

#br

#br
#br

*線形解析 [#gf9921aa]
**解析手順 [#n3a56ee4]
1. 「形状とメッシュ」にて解析モデルを作成し,メッシュ分割を行う.~
2. 「形状特性」にて,要素のタイプ・断面形状・部材の主軸に関する設定を行う.~
3. 「材料特性」にて,材料のヤング率・ポアソン比等の設定を行う.~
4. 「境界条件」にて,モデルの拘束条件および載荷条件の設定を行う.~
5. 「ジョブ」にて,解析オプション・出力結果に関する設定を行い,解析を実行する.~
6.  結果ファイルを開き,解析結果を確認する.~


#br
#br

**要素の選択 [#qe96fdca]
***梁要素 [#l74cf3ff]
梁要素種類(vol.B_P41)~
要素番号~
9:3次元トラス要素~
64:アイソパラメトリック3節点トラス要素~
52:弾性梁(非線形も可)~
98:横せん断を有する弾性梁(非線形も可)~
14:閉断面梁~
25:3次元の閉断面梁~
78:そりの無い3次元薄肉梁~
76:そりの無い3次元薄肉梁~
79:そりを含む3次元薄肉梁~
77:そりを含む3次元薄肉梁~
51:ケーブルエレメント~
31:弾性局面パイプ/直線梁要素~


***シェル要素 [#zd02ae6f]


#br
#br
**座標系(ローカル・グローバル) [#f9f8e459]



**片持ち梁モデル(梁) [#e0a9687c]

***モデル作成 [#e2964c5d]

-長さ単位 (基本操作)~
長さを選ぶ(私はmm)~

-形状とメッシュ (基本操作)~
ポイント 追加 → 座標入力する(0 0 0 とか)(x y z の順,スペース開ける) → 節点 追加 (ポイントの上から節点を作る)→~
要素 追加 → 作った節点と節点を左クリック → 右クリック確定 → 線ができる~

-再分割 (オペレーション)~
分割 → 上からx,y,zの順に値を入力 → 要素 → 分割したいとこ左クリック → 右クリック確定 → 分割できる~

-スイープ (オペレーション)~
節点 → 全て → OK (必要のない節点を消す?みたいな)~

-再番号付け (基本操作)~
全ての形状/メッシュ → OK (節点の番号をつけ直す?みたいな)~


***形状設定 [#g2aa48c6]
-新規(構造)▼ (形状特性)~
3ーD ▶ → 中実断面梁(簡単な断面のとき使える) → 形状選ぶ(長方形,正方形とか) → 寸法入力 →~
X,Y,Zの値は変えたり変えなかったり(今のところよく分かっていない) → 要素 追加 →~
追加したいところ左クリック(全体選択でも可) → 右クリック確定 → 左側の要素の数増えてたらできてる~


***材料設定 [#bfe74861]
-新規(構造)▼ (材料特性)~
有限剛性領域 → スタンダード → ヤング率とポアソン比入力(n/mm^2で入力) → 要素 追加 →~
追加したいところ左クリック(全体選択でも可) → 右クリック確定 → OK~


***境界条件設定 [#j0b785bb]
-(固定)~
新規(構造)▼ (境界条件)~
変位指定 → 変位X,Y,Zと回転X,Y,Zにチェックを入れて0にする(固定端のとき) → 節点 追加 →~
固定端にする方の節点を左クリック → 右クリック → 矢印出てきたらできてる~

-(ローラー)~
新規(構造)▼ (境界条件)~
変位指定 → 変位X,Yにチェックを入れて0にする(ローラー端のとき) → 節点 追加 →~
ローラー端にする方の節点を左クリック → 右クリック → 矢印出てきたらできてる~

-(載荷)~
新規(構造)▼ (境界条件)~
集中荷重 → 力Yに値を入力 → 節点 追加 → 載荷するところの節点を左クリック → 右クリック → 矢印出てきたらできてる~


***解析 [#z06f102b]
-新規▼ (解析ジョブ)~
構造 → チェック → Job check complete になったらOK → 実行 → Marc実行(1) → 3004になったら解析できてる~


***結果見る [#w2500455]
-ポストファイルを開く(モデルプロット結果処理メニュー)~
→ 変形形状の形式を 変形図と元の形状図▼にする → 設定 → 自動にする → 変形図見れる~



**簡易モデル(梁) [#s592f8fa]
***モデル作成 [#nd903af8]

-長さ単位 (基本操作)~
長さを選ぶ(私はmm)~

-形状とメッシュ (基本操作)~
ポイント 追加(ケーブルのポイントも作っておくと後で楽) → 座標入力する → 節点 追加 → 要素 追加 →~
作った節点と節点を左クリック → 右クリック確定 → 線ができる~

-再分割 (オペレーション)~
分割 → 上からx,y,zの順に値を入力 → 要素 → 分割したいとこ左クリック → 右クリック確定 → 分割できる~

-スイープ (オペレーション) (ケーブルを作る前にスイープをしなければならない!!)~
節点 → 全て → OK ~

-再番号付け (基本操作)~
全ての形状/メッシュ → OK ~

-形状とメッシュ (基本操作)~
節点 追加 → ケーブルの始点と終点をクリックして節点を作る(部材と重なっている節点番号は後で使うので覚えておく) →~
要素 追加 → 部材と重なっていないケーブルをクリック → 重なっているケーブル番号を下に入力 → 線ができる~


***ケーブルと部材をつなげる [#i650c11a]
-リンク~
節点ダイイング▼ → 新規▶ → 梁ピン結合 → 上の 節点 に後から作った節点番号を入力 → 下の 節点 に重なってるとこクリック → OK~


***形状設定 [#qa7efa4d]
-(部材)~
新規(構造)▼ (形状特性)~
3ーD ▶ → 中実断面梁 → 形状選ぶ → 寸法入力 → X,Y,Zの値は変えたり変えなかったり → 要素 追加 →~
追加したいところ左クリック → 右クリック確定 → 左側の要素の数増えてたらできてる~

-(ケーブル)~
新規(構造)▼ (形状特性)~
3ーD ▶ → 中実断面梁 → 形状選ぶ(ケーブルは円形かな?) → 寸法入力(多分直径) → X,Y,Zの値は変えたり変えなかったり → 要素 追加 →~
追加したいところ左クリック → 右クリック確定 → 左側の要素の数増えてたらできてる~


***材料設定 [#wda89a73]
-(部材)~
新規(構造)▼ (材料特性)~
有限剛性領域 → スタンダード → ヤング率とポアソン比入力(n/mm^2で入力) → 要素 追加 →~
追加したいところ左クリック → 右クリック確定 → OK~

-(ケーブル)~
-新規(構造)▼ (材料特性)~
有限剛性領域 → スタンダード → ヤング率とポアソン比入力(n/mm^2で入力) → 熱膨張 → 値を入力(材料によって変わる) →~
要素 追加 → 追加したいところ左クリック(全体選択でも可) → 右クリック確定 → OK~


***境界条件設定 [#p069ea27]
-(部材固定)~
新規(構造)▼ (境界条件)~
変位指定 → 変位X,Y,Zと回転X,Y,Zにチェックを入れて0にする → 節点 追加 →~
固定端にする方の節点を左クリック → 右クリック → 矢印出てきたらできてる~

-(ケーブル固定)~
新規(構造)▼ (境界条件)~
変位指定 → 回転X以外にチェックを入れて0にする(x軸が軸方向のとき) → 節点 追加 →~
固定にする方の節点を左クリック → 右クリック → 矢印出てきたらできてる~

-(載荷)~
新規(構造)▼ (境界条件)~
集中荷重 → 力Yに値を入力 → 節点 追加 → 載荷するところの節点を左クリック → 右クリック → 矢印出てきたらできてる~

-(ケーブルに温度を入れる)(プレストレス入れるとき)
新規(状態変数)▼ (境界条件)~
節点温度 → 温度入力 → 節点 追加 → 固定した方のケーブルの節点を選ぶ → OK~


-解析ジョブ以降は片持ち梁と同じ




#br
#br
*非線形解析 [#c6264e85]
**弾塑性解析 [#c5058a35]
-材料の構成則はテーブルで入力する
-収束条件を緩くしすぎる(0.5くらい)と挙動がおかしくなる
-増分荷重はテーブルのtimeで入力する
-荷重1を加えて計算したものを10回ループさせる=固定荷重増分と,
最終的に荷重10を加えるけど,その増分の仕方をプログラムが自動的に
調整する方法=自動荷重増分がある.
-注意すべき点として、点と点を結ぶ線を作成するとき、線の方向によってモデルが異なってしまう。
(ローカル座標等の影響か)~
-シェル要素はサーフェス→要素へ変換するやり方だと要素タイプが薄肉シェルで作られる.~
ただ,その要素タイプでは回らないので要素タイプを厚肉シェルにしたら解析してくれた.(原因はわからないが)


#br
#br
#br
#br
*動的解析 [#v6c07612]





#br
#br
#br
#br

*メモ [#r4cbb759]
#br
***他ページ [#he8ae88f]

-[[Marc/mentat2018メモ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Marc/Mentat2018%E3%83%A1%E3%83%A2]]~


-[[Marc/mentatメモ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Marc/Mentat%E3%83%A1%E3%83%A2]]~


-[[Marc接触メモ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?MarcMentat-%E6%8E%A5%E8%A7%A6]]~


#br
#br



***梁層数 [#b3efdd10]
--解析job-->解析ジョブパラメータ-->シェル・梁層数~
シェルの厚さを通して積分に使用する層の数を設定する。




#br
#br

***ばね要素 [#y0384e5f]
--形状特性--3D・軸対称・平面-->リンク(バネ/ダッシュポット)~
自由度を一方向設定してその設定した方向の自由度を考慮したバネになる。~
並進と回転両方設定可能。剛性はEAと思われる。~

#br
#br

***RBE2 [#rb0b9c3d]
--リンク-->RBE2~
(RBE2とは、Nastranで定義されているRigid Body element type 2を意味する。)~
ノードとノードの自由度を自由に設定することができる。~
RBE2:剛体接続(節点の変位を制御)~
RBE3:柔体接続(力の分布を制御)~
弾塑性解析で計算すると自由度Y方向を固定しているはずが、変形図を見るとY方向固定されていなくなってしまう。~
→非線形解析時は使えない?
#br
#br

***自由度固定(全て・単一・A~B) [#o9db6985]
--リンク-->自由度固定~
ノードとノードのを自由度に設定することができる。~
ばね要素とも同時に使える。~
#br
#br


***荷重ケースと解析job [#xa0d37e0]
--解析の種類(機械的、熱伝導など)に応じて、メインメニューの解析パネルにある LOADCASE プロセッサでロード~
ケースを指定する。~
メインメニューのLOADCASEでは、負荷の増加(境界条件の選択)、ステップ数、自動ステップか固定ステップか、~
およびこの負荷ケースの制御を指定できる。~
解析jobは、解析プロセスの全体的な流れを制御するために使用される。これには、解析クラス、ロードケースの~
ロードケースの選択、解析オプション、要求される結果、初期荷重、接触制御、その他のパラメータが含まれる。~
また、要素タイプの指定、ジョブの整合性チェック、ジョブの実際の投入もこのプロセッサで行う。~


#br
#br

***薄肉構造を定義するための規約 [#ye9bd91d]
断面を一連の線となるように入力しなければならない。ある線の終点が次の線の始点となるように~
する必要があるということ。二重にして折り返す必要があるときは、厚さを0にする。~
各ブランチをセグメントに分割する必要がある。断面の応力点が分岐点となる。分割数は常に偶数で~
なければならない。応力がt/10?以下の距離で離れている場合、1つの点に統合される。~


#br
#br


***崩壊解析 [#hf07068d]
-【非正定】非正定値マトリックスを強制的に解くコマンド~
崩壊点を通過する場合にだけ利用する~

#br
#br

***鋼材の破断 [#s6f72685]
-相当塑性ひずみをしきい値として,ひずみがある値を超えた場合に,要素を消すあるいは,~
拘束を解くといった手法がある.~
また,相当塑性ひずみではなく,延性破壊式(例えば大矢根式)に着目した手法をとる方法もある~

#br
#br

***時刻歴解析 [#wd0a20be]

-複数荷重ステップを定義した場合、前のステップの評価結果(応力、ひずみ)の状態から次のステップの負荷が作用する

#br
#br


***節点座標の出力 [#yd85d713]
-PARTICLE TRACKINGという機能で,リメッシュ機能を使用しても,元の節点があった位置(連続体力学で言う物質点のイメージ?)の変形後の座標や,
その位置での応力の値などを追跡して,各タイムステップで出力する機能がある


#br
#br


***テーブル [#u9a14b3b]
-Table機能は時間増分法に固定増分法であるAuto Loadを選択すると,ひとつのLoadCase内で,変化するようなEdge Loadをかけることができるが,AutoStepなど,Marcが自動で次のタイムステップ~
を決める増分手法を選んだ場合,ひとつのLoadCase内で,荷重は線形的にしか変化できない制約がある.ひとつのLoadCase内で複雑に変化するような荷重を定義したい場合は,ユーザーサブルーチン~
を使用して記述するしかない



#br
#br



***ガウスポイントの結果 [#ce263d12]
- MentatでJOBS> JOB RESULTS> OUTPUT FILE> Full Element and Node Printを選択すると、ガウス点要素の結果を出力ファイルに直接出力することができる。
点要素の結果を直接出力ファイルに出力する。

#br
#br

***計算精度 [#kd3b4afa]
-特異性比が1e-6以下だと、計算結果の信頼度はない。そのため、再計算などする必要がある。



#br
#br

***Marcが立ち上がらない [#x7e01703]
ショートカットのアイコンをクリックしてMarcを立ち上げても、

 Security check failed!

 Check that the license server is running
 or contact MSC Software Support.

                                 OK
と表示され、起動しない。
コマンドプロンプトには
 License server is not at or above the minimum required version.
 Verify that your license server is at FLEXlm version 11.13.1.0 or above
と表示される。


▸時間を空けて再起動したら、問題なく立ち上げることができた。~
もしかするとサーバーが落ちていて、それが復旧したために起動できるようになったのかもしれない。

#br
#br


**エラーメッセージ [#q47d21dc]
<終了番号:3015>~
 Unable to reduce the time step below the minimum value
 allowed, and to still satisfy the user criteria in the AUTO
 STEP history definition option. Either change the criteria
 or reduce the minimum time step allowed. 

ー時間ステップを許容される最小値以下にすることができず、なおかつAUTO STEP履歴定義オプションのユーザー基準を満たすことができない。ー~
#br

解決策
-荷重ケース->複数条件_パラメータ->最小比における減衰の追加に☑を入れる
-荷重ケース->解法管理->最大反復回数を上げる


#br
#br


<終了番号:3002>~

 Analysis has failed to converge to required convergence
 tolerances. One of several error conditions has been
 detected and the run aborted. The output will specify
 additional messages.

ー解析が収束許容値に収束しなかったー~
#br

解決策
-荷重ケース->解法管理->最大反復回数を上げる


#br
#br

<終了番号:2004>~

 The determinant of the stiffness matrix becomes zero ornegative when indicated node has been reached during the
 Gaussian elimination phase of the solution process. Thismeans that the stiffness matrix is non-positive definite.
 If this happens at the start of the analysis, the conditionis 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). Innonlinear cases, the structure may have buckled or reached a
 plastic limit load. In rubber analysis, it may also be dueto the strain state being in a region where the input data
 for the strain energy function is invalid. In contactanalysis with friction, lack of normal forces may result in
 friction being absent. If desired, the program may be forcedto continue by use of the parameter PRINT or the model or
 history definition option CONTROL. Either one of theseprocedures may be used for restart. Whenever a non-positive
 definite situation occurs one must exercise caution, as theresultant numerical solution may be infeasible.

ー剛性マトリックスが非正定値であるー~
ー非線形の場合は、構造物が座屈したり、塑性限界荷重に達している可能性があるー~
ー線形解析の場合、剛体変位(部材の平行移動/回転移動)が発生している場合があるー~
ースイープで節点がマージ(融合)できていない可能性がある。ー~

 

#br

解決策
-荷重ケース->解法管理->非正定に☑を入れる
-材料定数の確認~
-シェル要素のときはサーフェスを細かく設定し、異なる部材が接点を共有できているか確認する。(線形解析)~
-梁要素のときは接点で要素が繋がっているか確認する。(線形解析)~
-スイープをやり直す。~
-変位指定できちんと固定されているかを確認する。~


#br
#br
<終了番号:13>~
 Data errors have been detected during data input.  Refer to
 output for location of error.  Likely causes are misspelled
 keywords, mistyped lines or invalid input options.

ー入力データがエラーであるー~

トップ   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS