陽解法動解析

この新たに開発されたOptiStructの陽解法解析タイプ（ANALSIS=NLEXPL）は、OptiStructの陰解法解析と同じように、OptiStruct内部ソルバーとして開発されました。陽解法解析の入力データ（要素、材料、特性、荷重など）は陰解法解析と同じであり、出力データ構造も陰解法解析と同じです。

重要: この機能は、Altair Radiossとの統合による既存の陽解法動解析とは異なります（ANALYSIS=EXPDYN）。

このソリューションシーケンスでは、非線形陽解法有限要素解析が実行されます。非線形陽解法有限要素解析と非線形陰解法過渡解析の主な違いは、時間積分スキームです。非線形陽解法有限要素解析では、時間ステップは通常小さめであり、陰解法とは対照的に、陽解法ではマトリックスのアセンブリや反転は不要です。OptiStructの非線形陽解法ソリューションシーケンスは一般に、幾何学的大変位非線形性、材料非線形性、接触を含む、すべての主要な非線形フィーチャー（NLSTAT（LGDISP）など）をサポートしています。サブケースの継続は現時点ではサポートされていません。最適化も現在サポートされていません。

SMPとMPI（DDM）の並列処理は、OptiStructの非線形陽解法解析でサポートされています。

非線形性ソース

幾何学的非線形性

幾何学的非線形性を含む解析では、構造変形としての形状変化は構成方程式と釣り合い方程式の定式化に基づいて考慮されます。多くの工学的分野において、幾何学的非線形性に基づく大変位解析を使用する必要があります。金属成形、タイヤ解析、医療装置解析などの応用例がこれに該当します。

材料の非線形性

材料の非線形性には、現在の変形状態、変形の履歴、変形率、温度、圧力などに基づく材料の非線形挙動が関係します。

拘束と接触の非線形性

システム内の拘束の非線形性は、モデル内に運動学的拘束条件が存在する場合に発生します。モデルの運動学的自由度は、その動きに制約を課すことで拘束できます。OptiStructの陽解法では、MPCはサポートされていません。RBE2とTIEの接触については、拘束は運動学的に適用されます。

接触の場合は、拘束条件はペナルティ法によって適用されます。

Auto-Contactは、CONTACTバルクデータエントリでTYPEフィールドをAUTOに設定することで利用できます。

追従荷重

大変形が含まれる際、適用される荷重は構造の変形に依存することがあります。幾何学的には、適用される荷重（力または圧力）は、荷重のかかる位置においてモデルがどのように変形するかに基づき、その初期方向から変化し得ます。OptiStructでは、適用される荷重が追従荷重として扱われる場合、その荷重の向きまたは統合された大きさ、またはその両方は、解析全体を通して幾何形状の変化と共に更新されます。

適用される荷重は、FLLWERバルクおよびサブケースエントリやPARAM,FLLWERエントリを使用して、追従荷重として示すことができます。

注: 追従荷重は現在、DLOAD/TLOAD#を介して指定された荷重、すべての圧力荷重、FORCE1、FORCE2、MOMENT1、およびMOMENT2でサポートされています。

陽解法有限要素解析法

陽解法有限要素法では、時間-離散化方程式は、陽解法時間積分法を使用して解かれます。陽解法時間積分法は、中央差分スキームに基づいています。

中央差分法

中央差分法では、釣り合い方程式は次のようなかたちをとります：(1)

M a^{n} = f_{e} (u^{n}, t^{n}) + f_{d} (v^{n}, t^{n}) + f_{c} (u^{n}, t^{n}) + f_{h} (u^{n}, t^{n}) - f_{i} (u^{n}, t^{n})

ここで、

$M$: 集中質量マトリックス
$f_{e}$ 、 $f_{d}$ 、 $f_{c}$ 、 $f_{h}$ 、および $f_{i}$: は、それぞれ外力、減衰力、接触力、アワグラス力、要素内力のベクトルです。
$a^{n}$: 加速度ベクトルは釣り合い方程式から直接計算されます。

a^{n}

から、速度と変位のベクトルは、次のように更新できます：(2)

v^{n + \frac{1}{2}} = v^{n - \frac{1}{2}} + \frac{1}{2} (t^{n + \frac{1}{2}} - t^{n - \frac{1}{2}}) a^{n}

(3)

d^{n + 1} = d^{n} + (t^{n + 1} - t^{n}) v^{n + \frac{1}{2}}

ここで、

$t^{n}$: 現在時刻
$t^{n + 1}$: 次の時刻

次の時間増分が定義されます：(4)

Δ t^{n} = t^{n + 1} - t^{n}

(5)

Δ t^{n - 1} = t^{n} - t^{n - 1}

したがって：(6)

v^{n + \frac{1}{2}} = v^{n - \frac{1}{2}} + \frac{1}{2} (Δ t^{n - 1} + Δ t^{n}) a^{n}

(7)

d^{n + 1} = d^{n} + Δ t^{n} v^{n + \frac{1}{2}}

臨界時間ステップ

陰解法非線形過渡解析とは異なり、陽解法時間積分スキームは条件付きで安定しています。

陽解法解析は時間に合わせて前進します。各時間増分における時間ステップは、デフォルトで自動的に計算され（要素時間ステップがデフォルトです）、TSTEPEバルクデータエントリのTYPEフィールドを使用して、要素時間ステップと節点時間ステップの間で切り替えることができます。TSTEPEバルクデータエントリのDTMINフィールドを使用して、許容される最小節点時間増分を指定できます。陽解法動解析では、最も小さい上位10個の臨界時間ステップ（要素 / 節点）がデフォルトで.outファイルに出力されます。これは、PARAM, CRTELEMを使用して制御できます。

要素の時間ステップ

これは、非線形陽解法解析のデフォルトの時間ステップ制御タイプです。TSTEPEエントリのTYPEフィールドは、デフォルトでELEMに設定されます。

ソリッド要素
時間ステップサイズは次の条件を満たす必要があります：(8) $Δ t \leq \frac{2}{ω_{\max}}$

ここで、 $ω_{\max}$ はシステムの最大固有振動数を表します。

ソリッド要素の場合は、臨界時間ステップサイズは次のように計算されます：(9) $Δ t_{e} = \frac{l_{e}}{Q + {(Q^{2} + c^{2})}^{\frac{1}{2}}}$

ここで、

$c$

断熱音速

$Q$

体積粘性係数 $C_{0}$ および $C_{1}$

(10) $Q = C_{1} c + C_{0} l_{e} \max (0, - {\dot{ε}}_{k k})$

ここで、

$C_{0}$ および $C_{1}$

体積粘性係数は、デフォルト値がそれぞれ1.5と0.06である無次元定数です。

$l_{e}$

要素の特性長さ。

8節点六面体

(11) $l_{e} = \frac{V_{e}}{A_{e_{\max}}}$

10節点四面体

(12) $l_{e} = \frac{1}{{(B_{i j} B_{i j})}^{\frac{1}{2}}}$

6節点五面体

(13) $l_{e} = \frac{1}{{(B_{i j} B_{i j})}^{\frac{1}{2}}}$

4節点四面体

(14) $l_{e} = \frac{3 V_{e}}{A_{e_{\max}}}$

ここで、

$B_{i j}$

形状関数の対称勾配

$V_{e}$

六面体要素の体積

$A_{e_{\max}}$

六面体要素の全6フェイスのうちの最大面積
シェル要素
シェル要素の場合は、時間ステップサイズは次の式によって決定されます：(15) $Δ t = \frac{L}{c}$

ここで、 $c$ は、次の式で計算される音速です：(16) $c = \sqrt{\frac{E}{ρ (1 - ν^{2})}}$

ここで、

$E$

ヤング率

$ρ$

密度

$ν$

ポアソン比

$L$

特性長。四角形要素については次のように計算されます：

$L = \frac{A}{\max (L_{1}, L_{2}, L_{3}, L_{4})}$

ここで、

$A$

面積

$L_{1}, L_{2}, L_{3}, L_{4}$

三角形要素の辺の長さ。次のように計算されます：

$L = \frac{2 * A}{\max (L_{1}, L_{2}, L_{3})}$

ここで、

$A$

面積

$L_{1}, L_{2}, L_{3}$

要素の辺の長さ
要素マススケーリング
スケーリングされた要素臨界時間ステップ（DTFACによってスケーリング）がDTMINを下回る場合は、要素質量をスケーリングすることで、 $Δ t_{e}$ を増大させることができます。これが可能となる理由は、要素時間ステップ式には、材料密度（ $ρ$ ）に依存する音速の項（ $c$ ）が含まれているからです。

節点時間ステップ

時間ステップの制御をデフォルトの要素時間ステップから節点時間ステップに切り替えるには、TSTEPEバルクエントリのTYPEフィールドをNODAに設定します。

節点時間ステップは次のように計算されます：(17)

Δ t_{n} = \sqrt{\frac{2 m_{n}}{k_{n}}}

ここで、

$m_{n}$: 節点質量
$k_{n}$: 節点剛性（要素剛性から計算されます）

節点剛性は次のように計算されます：

各要素について、臨界時間ステップ

Δ t_{e}

が最初に計算され、各節点は同じ時間ステップ

Δ t_{e}

を持つものと想定されて、その後各節点について、次の式から節点剛性を推定できます。(18)

Δ t_{e} = \sqrt{\frac{2 m_{e_{i}}}{k_{e_{i}}}}

ここで、

$i$: 要素のi番目の節点
$m_{e_{i}}$: i番目の節点の節点質量
$k_{e_{i}}$: この要素のi番目の節点の節点剛性

したがって、i番目の節点の節点剛性は次のとおりです：(19)

k_{e_{i}} = \frac{2 m_{e_{i}}}{Δ t_{e}^{2}}

最後の節点剛性は次のとおりです：(20)

k_{n} = \sum_{e} k_{e_{i}}

$k_{n}$ を使用して、節点臨界時間ステップ $Δ t_{n}$ を計算できます。

節点マススケーリング

スケーリングされた節点臨界時間ステップ（DTFACによってスケーリング）がDTMINを下回る場合は、節点質量 $m_{n}$ をスケーリングすることで、 $Δ t_{n}$ を増大させることができます。

アワグラス制御

アワグラス制御は、PARAM,HOURGLSエントリまたはHOURGLSエントリを使用してアクティブ化できます。これらのエントリは、アワグラス制御パラメータ（HGTYPとHGFAC）を調整するためのアクセスも可能にします。

HOURGLSエントリが入力された場合に、このエントリをアクティブにするには、対応するプロパティエントリのHGIDフィールドを通じてこのエントリを選択する必要があります。HGIDフィールドを介したHOURGLSエントリによって、PARAM,HOURGLSにより定義された設定が上書きされます。

ソリッド要素の場合

MAT1/MATS1材料を使用したソリッド要素については、次の2タイプのアワグラス制御が用意されています：

タイプ1（FlanaganとBelytschko、1981）は、粘性減衰を伴う望ましくないアワグラスモードに抵抗を与えます。
タイプ2（Puso、2000）は、拡張された想定ひずみ物理的安定化を使用して、粗いメッシュ精度で高い計算効率を提供します。タイプ2は、1次CHEXA要素のMAT1/MATS1材料のデフォルトのアワグラスタイプとして選択されます。

タイプ1とタイプ2のアワグラス制御の実装はよく似ていますが、アワグラス力の計算方法は異なります。

注: タイプ1と比べて、タイプ2は計算量が多い反面、アワグラスモードを排除する性能が優れています。タイプ2の唯一の短所は、大きな塑性変形を伴う曲げ問題で過度に剛性が高い応答をもたらす可能性があることです。

MATHEエントリの場合、デフォルトのアワグラス制御はタイプ4（Reese、2005）です。MATHEエントリにタイプ2を使用することもできます。

ソリッド要素の低減積分の場合（ISOPE=URI/AURI）、アワグラス制御はデフォルトでオンになります。

	アワグラス制御（ソリッド要素ベース）
要素	通常の要素（`ISOPE`=URI）	通常の要素（`ISOPE`=AURI）	通常の要素（`ISOPE`=SRI）
CHEXA （1次）	アワグラス制御はデフォルトでオンになります。 ¹ MAT1/MATS1の場合アワグラスタイプ2 MATHEの場合アワグラスタイプ4	アワグラス制御はデフォルトでオンになります。¹ デフォルトは次のとおりです： MAT1/MATS1の場合アワグラスタイプ2 MATHEの場合アワグラスタイプ4	アワグラス制御はデフォルトでオンになりません。 ²
CTETRA （2次）	アワグラス制御は不要です。	アワグラス制御は不要です。	アワグラス制御は不要です。
CPENTA （1次）	NA	NA	アワグラス制御はデフォルトでオンになりません。 ²
CTETRA（1次）	アワグラス制御は不要です。	アワグラス制御は不要です。	アワグラス制御は不要です。

シェル要素の場合

シェル要素については、次の2タイプのアワグラス制御のみが用意されています：

タイプ1（FlanaganとBelytschko - 粘性形式）
タイプ2（FlanaganとBelytschko - 剛性形式）タイプ2は、CQUAD4のMAT1/MATS1材料のデフォルトのアワグラスタイプとして選択されます。

	アワグラス制御（シェル要素ベース）
要素	Belytschko-Tsay（`ISOPE`=BT）	完全投影を使用したBelytschko-Wong-Chiang（`ISOPE`=BWC）
CQUAD4	アワグラス制御はデフォルトでオンになります。 ¹ MAT1/MATS1の場合アワグラスタイプ2	アワグラス制御はデフォルトでオンになります。 ¹ MAT1/MATS1の場合アワグラスタイプ2

材料

次の表では、サポートされている材料での各種アワグラス制御タイプとデフォルト値を示します。

	アワグラス制御（材料ベース）
材料	タイプ1 ソリッドとシェル： Flanagan-Belytschko粘性形式	タイプ2 ソリッド： Puso拡張想定ひずみ剛性形式シェル： Flanagan-Belytschko剛性形式	タイプ4 ソリッド： Reeseアワグラス制御シェル：タイプ4はシェルでのサポートなし
MAT1/MAT2/MAT8/MATS1	使用可能 ²	デフォルト ⁶	NA
MATHE	NA	使用可能 ²	デフォルト ⁶
MATVE	NA	使用可能 ²	デフォルト ⁶

問題の設定

入力

アクティブ化：
非線形陽解法サブケースは、ANALYSIS=NLEXPLによって指定できます。TTERMサブケースエントリは、終了時刻を定義するために必要です。また、対応するTSTEPEバルクデータエントリを指しているTSTEPEサブケースエントリを、非線形陽解法解析で使用することもできます。TSTEPEサブケースエントリが定義されていない場合は、ANALYSIS=NLEXPLをTTERMと組み合わせて使用する必要があります。このサブケースエントリが定義されている場合は、TTERMとTSTEPEを併用するだけで、陽解法非線形サブケースを指定できます。非線形陽解法解析は常に大変位解析です。
初期条件：
初期条件を定義するには、ICサブケースエントリをTICバルクデータエントリと組み合わせて使用します。
荷重：
荷重を定義するには、LOAD、DLOAD、およびTLOAD#バルクデータエントリを使用します。これらは、DLOADサブケースエントリを使用してサブケース内で参照される必要があります。LOADサブケースエントリまたはTLOAD#バルクエントリを介した参照では、荷重に対してFORCE、FORCE1、FORCE2、MOMENT、MOMENT1、MOMENT2、PLOAD2、PLOAD4、GRAV、ACCEL2、およびSPCDエントリのみがサポートされています。
境界条件：
境界条件は、対応するSPCサブケースエントリで参照されているSPCバルクデータを介して適用できます。MPCは現在サポートされていません。
サポートされている要素：

ソリッド要素

4節点CTETRA、10節点CTETRA、8節点CHEXA、6節点CPENTAの各要素がサポートされています。

シェル要素

CTRIA3とCQUAD4がサポートされています。

1次元要素

CBUSH、CBEAMおよびCBAR要素がサポートされています。

現在、陽解法解析ではCBAR/CBEAM要素に対しては、Belytschko-Schwer Beamの定式化のみがサポートされています。

質量要素

CONM2がサポートされています。
注:
- シェル要素の要素またはプロパティのオフセットは、陽解法解析でサポートされています。
- CBUSH要素の場合、PBUSH定義のMiフィールドは質量と慣性の計算に使用されます。詳細についてはリファレンスガイド内のPBUSHをご参照ください。
- CBEAM、CBAR要素の場合、
  - PBEAM/PBAR上の継続行は、陽解法解析ではサポートされません。
  - ピンフラグ（PAとPB）は陽解法解析でサポートされています。
サポートされている材料：
MAT1、MAT2、MAT8、MATS1、MATHE、MATVE材料はサポートされています。MATHE内のFOAM材料は、陽解法解析でサポートされていますが、非線形陰解法解析ではサポートされていません。MATVEエントリをMATHEエントリ下で定義する必要があります。
積分スキーム：
陽解法解析では、要素積分スキームは、PSOLID、PLSOLID、PLSOLID、PSHELL、PCOMP、PCOMPG、PCOMPPエントリのISOPEフィールドを使用するか、PARAM,EXPISOPを使用して変更できます。ISOPEフィールドの設定は、PARAM,EXPISOPの設定を上書きします。

例：

SUBCASE 10
   ANALYSIS=NLEXPL
   SPC = 1
   DLOAD = 2
   TSTEPE = 2
   NLOUT = 23
   IC = 12
   TTERM = 2.0
.
.
BEGIN BULK
TSTEPE,2,ELEM,0.8
NLOUT,23,NINT,12
IC,12,33,3,0.2
SPC,1,45,123,0.0
TLOAD1,2,3,,0,8
TABLED1,8
+,0.0,0.0,2.0,8.0,ENDT,ENDT

出力

一般的な出力エントリ（DISPLACEMENT、VELOCITYおよびACCELERATION）を使用して、非線形陽解法解析の対応する出力を要求できます。NLOUTサブケースおよびバルクデータエントリを使用して中間結果を要求できます。NINTパラメータのみがサポートされています。

NLOUTバルクデータエントリおよびNLOUTサブケース情報エントリを使用して、増分出力を制御できます。非線形陽解法解析では、NLOUTに対してNINTフィールドのみがサポートされています。NLADAPTエントリは非線形陽解法解析ではサポートされておらず、TSTEPEエントリ以外のTSTEP#エントリはサポートされていません。

現在は、Hyper3D（_expl.h3d）およびHyperGraphプレゼンテーション形式（_expl.mvw）ファイルのみがサポートされています。非線形陽解法解析の結果は通常の.h3dおよび.mvwファイルには出力されませんが、代わりにそれぞれ_expl.h3dファイル、_expl.mvwファイルに出力されます。

_expl.h3d: 変位、回転、速度、加速度、ひずみ、応力、塑性ひずみ、CBUSH要素力、複合材の応力、複合材のひずみ、複合破壊指数の等高線が出力されます。
_expl.mvw: ひずみエネルギー、弾性接触エネルギー、塑性接触エネルギー、運動エネルギー、アワグラスエネルギー、および塑性散逸エネルギーの曲線が出力されます。
.out: 陽解法では、.outファイルに、時間サイクル情報（PARAM,NOUTCYCに基づく）、現在時刻、現在の時間ステップ、最大ひずみエネルギー、この情報が出力される要素ID、運動エネルギー、接触仕事、合計エネルギー、最大貫通、この最大貫通に関連付けられた節点ID、最大垂直仕事、この最大垂直仕事に関連付けられた節点ID、質量変化率が含まれます。質量変化率は、マススケーリング後のスケーリングされた質量の変化に関する情報であり、これは（現在の質量 - 元の質量）/（元の質量）で算出されます。
_expl.cntf: メインのサーフェスの接触力出力結果を含んだASCIIファイルで、CONTF I/OオプションエントリでOPTIフォーマットが指定されている場合にアクティブになります。出力には、法線/接線力、大きさ、接触面積が含まれます。この出力は、陽解法時間ステップ毎に用意されます。; このファイルの出力間隔は、NLOUTエントリのNINTフィールドを使用して制御できます。

表 1. 非線形陽解法解析の簡単なまとめ
非線形陽解法解析	サブケースまたは入出力	バルクデータ	コメント
アクティブ化：
サブケースタイプ	ANALYSIS=`NLEXPL`（オプション）	NA	TSTEPEが指定されていない場合は、ANALYSIS=`NLEXPL`は必須です。
非線形陽解法のアクティブ化	TTERM（必須） TSTEPE（オプション）	TSTEPE（オプション）	TSTEPEが指定されていない場合は、ANALYSIS=`NLEXPL`は必須です。
荷重：
節点荷重	LOAD、DLOAD	サブケース内のLOADが使用されている場合： FORCE、FORCE1、FORCE2、MOMENT、MOMENT1、およびMOMENT2。サブケース内のDLOADが使用されている場合： TLOAD1またはTLOAD2。 DLOADを使用して、複数のTLOADiデータを結合できます。節点荷重では、TLOADiデータの`EXCITEID`に設定できる値は、FORCE、FORCE1、FORCE2、MOMENT、MOMENT1、およびMOMENT2です。	`TYPE` このケースでTLOADiデータのフィールドは、0またはLOADに設定できます。
サーフェス荷重	LOAD、DLOAD	サブケース内のLOADが使用されている場合： PLOAD2およびPLOAD4。サブケース内のDLOADが使用されている場合： TLOAD1またはTLOAD2。 DLOADを使用して、複数のTLOADiデータを結合できます。サーフェス荷重では、TLOADiデータの`EXCITEID`は、PLOAD1およびPLOAD4に設定できます。	`TYPE` このケースでTLOADiデータのフィールドは、0またはLOADに設定できます。
ボディ荷重	LOAD、DLOAD	サブケース内のLOADが使用されている場合： GRAVおよびACCEL2。サブケース内のDLOADが使用されている場合： TLOAD1またはTLOAD2。 DLOADを使用して、複数のTLOADiデータを結合できます。物体力荷重では、TLOADiデータの`EXCITEID`は、GRAVおよびACCEL2に設定できます。	`TYPE` このケースでTLOADiデータのフィールドは、0またはLOADに設定できます。
強制変位、強制速度、強制加速度	LOAD、DLOAD	サブケース内のLOADが使用されている場合： SPCDまたはSPCDを使用した強制変位、強制速度、または強制加速度。サブケース内のDLOADが使用されている場合： TLOAD1またはTLOAD2。 DLOADを使用して、複数のTLOADiデータを結合できます。強制荷重では、TLOADiデータの`EXCITEID`は、SPCまたはSPCDに設定できます。	`TYPE` TLOADiデータのフィールドは、以下の値に設定できます： 1またはDISP 強制変位の場合 2またはVELO 強制速度の場合 3またはACCE 強制加速度の場合
追従荷重	FLLWER	FLLWER PARAM,FLLWER	陰解法非線形解析と同様に、荷重として追従荷重を選択できます。追従荷重は現在、DLOAD/TLOAD#を介して指定された荷重、すべての圧力荷重、FORCE1、FORCE2、MOMENT1、およびMOMENT2でサポートされています。
境界条件：
単点拘束	SPC	SPC
初期条件：
初期変位	TIC	IC
初速度	TIC	IC
時間ステップコントロール：
基本的な時間コントロール	TSTEPE	TSTEPE	`TYPE` TSTEPEエントリのフィールドを使用して、要素時間ステップコントロールまたは節点時間ステップコントロールのどちらかを選択します。 `DTMIN` フィールドでは、最小時間ステップを定義できます。時間ステップがこれより小さい場合、節点 / 要素マススケーリングがアクティブ化されます。 `DTFAC` フィールドでは、安定した時間増分のためのスケールファクターを定義できます。
質量要素：
質量要素のサポート		CONM2がサポートされています。
構造要素：
サポートされている構造要素	NA	1次元要素： CBUSH、CBEAMおよびCBARがサポートされます。シェル要素 CTRIA3およびCQUAD4要素。ソリッド要素 4節点CTETRA、10節点CTETRA、8節点CHEXA、6節点CPENTAの各要素。
積分スキーム	NA	`ISOPE` PSOLID、PLSOLID、またはPSHELLのフィールド。 PARAM,EXPISOP（パラメータはソリッド要素でのみサポートされています）。	`ISOPE` フィールドは、PARAM,EXPISOPで定義された設定を上書きします。積分スキームの詳細については、ユーザーズガイド内の要素をご参照ください。
拘束：
剛体のサポート	NA	RBE2、RBE3およびRBODYがサポートされています。
材料：
サポートされている材料	NA	シェル： MAT1、MAT2、MAT8およびMATS1. ソリッド： MAT1、MATS1、MATVE、およびMATHE。	MATS1の場合：陰解法でサポートされているMATS1上の材料に加えて、Johnson-Cookとクラッシャブルフォームの材料もサポートされています。 MATHEの場合：陰解法でサポートされているすべての材料モデルは陽解法でもサポートされており、それに加えて、陽解法ではMATHEでフォーム材料がサポートされています。
プロパティ：
サポートされているプロパティ	NA	PSHELL、PSOLID、PLSOLID、PCOMP、PCOMPG、PCOMPP（PLY/STACKと共に）
接触：
サポートされている接触タイプ	NA	CONTACTおよびTIE	N2SとS2Sの接触離散化がサポートされています。 SMALL、FINITE、およびCONSLIの接触がサポートされています。 Auto-Contactは、CONTACTバルクデータエントリで`TYPE`フィールドをAUTOに設定することでサポートされます。陽解法でのTIEの場合： 1 運動学的TIEのみがサポートされています。すなわち、ペナルティに基づいた手法を使用する代わりに、運動条件が正確に拘束されます。 2 運動TIE内の階層はサポートされていません（すなわち、TIEのセカンダリ（旧称“スレーブ”）節点は別のTIE内のメイン（旧称“マスター”）節点になることはできません）。 3 過剰に拘束されたTIEは無視されます（.femファイル内の入力の順序に基づいて、そのようなケースの最初の拘束のみが保持されます）。 4 このような階層と過剰に拘束されたTIE節点はすべて、*_badtied.femファイル内の節点SET内に出力されます。
座標系：
サポートされているユーザー定義の座標系	NA	CORD2R、CORD1C、CORD2C、CORD1S、およびCORD2S
出力：
ASCIIの出力	NA	PARAM,NOUTCYC	陽解法の時間サイクルサマリーおよび対応する情報（時間ステップ、エネルギー、最大貫通、質量変化率など）のみが.outファイルに出力されます。PARAM,NOUTCYCを使用して、.outファイルへのサマリー出力の頻度を選択できます。
バイナリファイル出力	DISP、VELOCITY、ACCELERATION、STRESS、STRAIN（塑性ひずみを含む）CBUSH、FORCE、CSTRESS、CSTRAIN、CFAILURE	NA	結果は_expl.h3dファイルと_expl.mvwファイルのみに出力されます。 _expl.h3d 変位、回転、速度、加速度、応力、ひずみ、CBUSH力、塑性ひずみ複合材応力、複合材ひずみ、複合材破壊指数の結果が出力されます。 _expl.mvw ひずみエネルギー、弾性接触エネルギー、塑性接触エネルギー、運動エネルギー、アワグラスエネルギー、および塑性散逸エネルギーの曲線が出力されます。
出力制御	NLOUT	NLOUT	陽解法解析では`NINT`フィールドのみがサポートされています。非線形陽解法解析では、NLADAPTエントリはサポートされていません。
その他：
大変位解析	NA	NA	陽解法非線形解析は、デフォルトでは大変位非線形解析です。
アワグラス制御		HOURGLS（`HGID`フィールドは、PSOLID/PLSOLID/PSHELL上でこのカードを参照します） PARAM,HOURGLS	デフォルトのアワグラス値は、PSOLID/PLSOLID/PSHELLエントリで参照されるHOURGLSエントリ、またはPARAM,HOURGLSによって上書きされます。 `HGID` HOURGLSエントリを介したは、PARAM,HOURGLSを上書きします。詳細については、アワグラス制御をご参照ください。
最適化：
最適化のサポート	サポートされていません。	サポートされていません。		サポートされていません。

¹ デフォルト値は、ユーザー定義のPARAM,HOURGLSまたはHOURGLSエントリ（プロパティエントリのHGIDで参照されます）によって上書きできます

² ユーザーは、PARAMHOURGLSまたはHOURGLSエントリ（プロパティエントリのHGIDで参照されます）を使用してアワグラス制御をオンにできます

³ ソリッド要素については、PSOLID/PLSOLIDエントリのISOPEフィールドを使用して、積分スキームを切り替えることができます

⁴ CTRIA3要素にはアワグラス制御は不要です

⁵ シェル要素については、PSHELLエントリのISOPEフィールドを使用して、積分スキームを切り替えることができます

⁶ デフォルト値は、ユーザー定義のPARAM,HOURGLSまたはHOURGLSエントリ（プロパティエントリのHGIDで参照されます）によって上書きできます。MAT1/MATS1/MATHEについては、1次CHEXA要素の場合のみにデフォルト値が適用されます。CPENTA要素については、ユーザーは、必要に応じてアワグラス制御をオンにする必要があります。

⁷ ここで列挙している一部の材料は、シェルに対してはサポートされていません（MATHEやMATVEなど）。