/PROP/TYPE14（SOLID）

ブロックフォーマットキーワードこのプロパティセットは、一般的なソリッドプロパティセットの定義に使用されます。

フォーマット

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)
/PROP/TYPE14/`prop_ID`/`unit_ID`または/PROP/SOLID/`prop_ID`/`unit_ID`
`prop_title`
`I`_solid	`I`_smstr	`I`_ale	`I`_cpre	`I`_tetra10	`I`_npts	`I`_tetra4	`I`_frame	`d`_n
`q`_a		`q`_b		`h`		$λ_{v}$		$μ_{v}$
$Δ t_{\min}$		`Vdef_min`		`Vdef_max`		`APS_max`		`COL_min`

Sol2SPHオプションをアクティブにするオプション行 11

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
`Ndir`	`sphpart_ID`

定義

フィールド	内容	SI単位の例
`prop_ID`	プロパティの識別子（整数、最大10桁）
`unit_ID`	単位識別子。（整数、最大10桁）
`prop_title`	プロパティのタイトル（文字、最大100文字）
`I`_solid	ソリッド要素定式化フラグ。 1 2 = 0 /DEF_SOLID内の値を使用。 = 1　/DEF_SOLIDが定義されていない場合のデフォルト標準的な8節点ソリッド要素、1積分点。直交変形モードと剛体変形モードの補正を伴う粘性アワグラス定式化（Belytschko）。 = 2 標準的な8節点ソリッド要素、1Gauss積分点。直交性を伴わない粘性アワグラス定式化（Hallquist）。 = 14 HA8ロックフリー8節点ソリッド要素、共回転、完全積分、さまざまな数のGauss点。 = 16 2次20節点ソリッド要素、完全積分、さまざまな数のGauss点。 = 17 2次8節点ソリッド要素、完全積分、固定222 Gauss積分点。 = 18 8節点ソリッド要素、共回転、完全積分、固定222 Gauss積分点、せん断ロックフリー、`I`_cpreおよび`I`_smstrのデフォルトは材料に基づきます。 = 24 HEPH8節点ソリッド要素、1Gauss積分点、共回転系定式化、物理的なアワグラス安定化。（整数）
`I`_smstr	微小ひずみ定式化フラグ。 4 = -1 要素タイプと材料則を基に最適な値（/DEF_SOLID参照）を自動的に定義 = 0 右記内の値を使用。 /DEF_SOLID = 1 時間=0以降の微小ひずみ = 2 Radioss Engine（/DT/BRICK/CST）で微小ひずみ定式化に切り替えが可能である完全幾何非線形 = 3 時間=0以降の簡易微小ひずみ定式化（実現象に即していない定式化）。 = 4 /DEF_SOLIDが定義されていない場合のデフォルト。完全幾何非線形性（/DT/BRICK/CSTによる影響はありません） =10 ラグランジュ型全ひずみ。 = 11 t = 0からの全微小ひずみ定式化 =12 Radioss Engine（/DT/BRICK/CST）で全微小ひずみ定式化に切り替えが可能なLagrangeタイプ全ひずみ（整数）
`I`_ale	フレームワーク 12 = 0（デフォルト） Lagrange = 1 ALE = 2 EULER
`I`_cpre	定圧定式化フラグ 5 `I`_solid = 14、17、18または24の場合にのみ有効 = -1 要素タイプと材料則を基に最適な値（/DEF_SOLID参照）を自動的に定義 = 0 `I`_solidの値に応じて使用される定式化 1 = 1 `I`_solid = 17の場合のデフォルト。体積ロックを阻止するための定圧定式化。 $ν = 0.5$ である非圧縮性材料に使用。 = 2 使用される定式化は塑性の関数。これにより、材料が圧縮性である際の弾性領域および材料が非圧縮性になる際の塑性領域の正しいモデリングが可能に。弾塑性材料則に対してのみ使用できます。 = 3 `I`_solid=14または24の場合のデフォルト。フォームのような圧縮材料のための低減のない圧力積分（整数）
`I`_tetra10	10節点四面体要素定式化フラグ 7 = 0 /DEF_SOLIDの値を使用 = 2 4つの積分点および、/DT1/BRICKをともなう/TETRA4と同じ時間ステップの2次式による/TETRA10定式化。 = 3 4つの積分点および/TETRA4と同じ時間ステップを伴う2次式による/TETRA10定式化（形状の悪い要素に対する安定性が低い）。 = 1000　/DEF_SOLIDが定義されていない場合のデフォルト 4つの積分点を伴う2次式による/TETRA10定式化（整数）
`I`_npts	積分点の数 6 `I`_solid =14、16に対してのみ有効（整数） = ijk （デフォルト=222）： 2 < i,j,k < 9（`I`_solid =14の場合） 2 < i,k ≤ 3、2 ≤ `j` < 9（`I`_solid = 16の場合）ここで、 `i` r方向の積分点の数 `j` s方向の積分点の数 `k` t方向の積分点の数
`I`_tetra4	4節点四面体要素定式化フラグ 7 = 0 /DEF_SOLIDの値を使用 = 1 節点ごとの6種類の自由度と4つの積分点を伴う2次式による/TETRA4定式化 = 3 節点圧力平均をともなう線形TETRA4。体積ロックを低減。 = 1000　/DEF_SOLIDが定義されていない場合のデフォルト１つの積分点を伴う線形/TETRA4定式化（整数）
`I`_frame	要素座標系定式化フラグ 8 2次元ソリッド要素、および3次元ソリッド要素（`I`_solid ==1、2または17）。`I`_solid= 14または24は常に共回転定式化を使用。 = -1 要素タイプと材料則を基に最適な値（/DEF_SOLID参照）を自動的に定義 = 0 /DEF_SOLID内の値を使用。 = 1　/DEF_SOLIDが定義されていない場合のデフォルト非共回転定式化。 = 2 共回転定式化を使用。大きい回転を伴うモデルに推奨されます。（整数）
`d`_n	安定化のための数値減衰 9 `I`_solid =24に対してのみ有効デフォルト = 0.1（実数）
`q`_a	2次体積粘性。デフォルト = 1.10（実数）デフォルト = 0.0（/MAT/LAW70の場合）
`q`_b	線形体積粘性。デフォルト = 0.05（実数）デフォルト = 0.0（/MAT/LAW70の場合）
`h`	アワグラス粘性係数。 `I`_solid =1、2に対してのみ有効デフォルト = 0.10（実数）、0.0 <`h` < 0.15である必要があります。
$λ_{v}$	数値Navier Stokes粘性 $λ_{v}$ デフォルト = 0.0（実数）
$μ_{v}$	数値Navier Stokes粘性 $μ_{v}$ デフォルト = 0.0（実数）
$Δ t_{\min}$	ソリッド要素に対する最小時間ステップ。 /DT/BRICK/CSTまたは/DT/BRICK/DELを使用する際にのみ有効デフォルト = 0.0（実数）	$[s]$
`Vdef_min`	ソリッド要素を削除するための最小体積比（V/Vo）。`Ndir` ≥ 1の場合、要素は削除され、SPH要素がアクティブになります。デフォルト = 0
`Vdef_max`	ソリッド要素を削除するための最大体積比（V/Vo）。`Ndir` ≥ 1の場合、要素は削除され、SPH要素がアクティブになります。 0と異なる場合のみ使用されます。デフォルト = 0
`ASP_max`	ソリッド要素を削除するための最大アスペクト比。`Ndir` ≥ 1の場合、要素は削除され、SPH要素がアクティブになります。 0と異なる場合のみ使用されます。デフォルト = 0
`COL_min`	ソリッド要素を削除するための最小つぶれ値。`Ndir` ≥ 1の場合、要素は削除され、SPH要素がアクティブになります。デフォルト = 0
`Ndir`	各ソリッド要素の方向ごとの粒子の数（整数） 11 = 1 各方向に1つの粒子 = 2 各方向に2つの粒子で、合計8つの粒子 = 3 各方向に3つの粒子で、合計27の粒子（整数）
`sphpart_ID`	Sol2SPHのSPHプロパティを記述するパート識別子（整数）

例

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/2
unit for prop
#              MUNIT               LUNIT               TUNIT
                  kg                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/PROP/SOLID/1/2
Solid
#   Isolid    Ismstr      Iale     Icpre  Itetra10     Inpts   Itetra4    Iframe                  dn
        24        -1         0        -1         2         0         3        -1                   0
#                q_a                 q_b                   h            LAMBDA_V                MU_V
                   0                   0                   0                   0                   0
#         deltaT_min            vdef_min            vdef_max             ASP_max             COL_min
                   0                   0                   0                   0                   0
#     NdirSPHPART_ID
         0         0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

I_solid - ソリッド要素の定式化

ほとんどの場合、I_solid = 24（HEPH）六面体要素が計算コストと品質の間で最良の妥協点です。
要素定式化I_solid = 1、2および24は1 Gauss積分点を伴う低減積分要素、I_solid = 14、17および18は完全積分要素です。
I_solid = 24（HEPH）ソリッド要素は、I_shell = 24（QEPH）シェル要素によって使用されるアワグラス定式化と類似した物理的アワグラス定式化を使用します。このアワグラス定式化は、I_solid = 1または2で使用される粘性アワグラス定式化よりも良い結果を返します。

I_solid = 18、I_cpreおよびI_smstrのデフォルト値は材料に依存し、これあの推奨される値を使用：

デフォルト	材料則
`I`_cpre = 2	2、21、22、23、24、27、36、52、79、81、84
`I`_cpre = 3	12、14、15、25、28、50、53、68および $ν \leq 0.49$ の場合、1、13、16、33、34、35、38、40、41、70、および77 70, 77
`I`_cpre = 1	その他のすべての則および $ν \geq 0.49$ の場合、1、13、16、33、34、35、38、40、41、70、および77
`I`_smstr = 10	38、42、62、69、82、88、92、94、95、100および101
`I`_smstr = 11	70
`I`_smstr = 1	28
`I`_smstr = 2	すべてのその他の材料則

2次元ソリッド要素
以下の要素定式化は、/QUAD要素の使用時に2次元解析についてサポートされます。
- I_solid = 2、17
- I_smstr= 4
- I_frame = 1、2
- I_cpre = 1、2
自動設定オプションI_smstr = I_cpre = I_frame=-1を使用する場合、これらのオプションの値は、要素定式化、要素タイプ、材料に基づいて最適なオプションを使用して定義されます。あるいは、I_smstr = I_cpre = I_frame=-2と定義すると、このプロパティで定義されているこれらのオプションの値が、要素タイプと材料則に基づいた最適な値(/DEF_SOLID参照) で上書きされます。Radiossで定義された値を確認するには、Starter出力ファイルの“PART ELEMENT/MATERIAL PARAMETER REVIEW”セクションをご確認ください。

I_smstr - 微小ひずみ定式化フラグ

微小ひずみ定式化（I_smstr =1、3、11）、または微小ひずみ定式化（I_smstr =2、12）に切り替わる要素については、材料側で計算されるひずみと応力は工学ひずみと工学応力です。それ以外の場合は、真ひずみとCauchy応力になります。
I_smstr = 10、12は、全ひずみ定式化を使用する材料則1、38、42、62、69、82、88、92、94、100および101とのみ適合性があります。一般的に、左のCauchy-Greenひずみが使用されます。ユーザー則については、変形勾配テンソルと右ストレッチテンソルが使用できます。
I_smstr = 11は材料則70（フォーム）向けに開発され、全工学ひずみを使用している材料則（材料則1、38や70等）のみに適合しています。一般に、I_smstr =1の場合はより安定した結果を得ることができます。
I_smstr=12は/DT/BRICK/CSTと使用することで小さい時間ステップの要素を自動的にLagrange型全ひずみから全微小ひずみ（I_smstr = 11）に切り替えることができます。ただし、これらはひずみ定式化の変更中に応力の非連続性を多少生じる場合があります。
Radioss Engineオプション/DT/BRICK/CSTは、I_smstr = 2または12を使用するソリッドプロパティのみと機能します。

バージョン2017より、シミュレーション中に体積が負になるLagrange要素は、シミュレーションを継続できるよう、自動的にひずみ定式化を切り替えます。これが発生すると、WARNINGメッセージがEngine出力ファイルに出力されます。以下のオプションがサポートされています。

要素タイプおよび定式化	ひずみ定式化	負の体積の処理方法
/BRICK `I`_solid=1、2、14、17、18、24 /TETRA4、`I`_tetra4 = 1000 /TETRA10、`I`_tetra10 = 1000	完全に幾何学的非線形 `I`_smstr = 2、4	負の体積になる前のサイクルの要素形状を使用した微小ひずみに切り替え
Lagrange型全ひずみ。 `I`_smstr = 10、12	time=0.0における要素形状でのLagrange型全ひずみ

要素タイプおよび定式化

ひずみ定式化

負の体積の処理方法

/BRICK

I_solid=1、2、14、17、18、24

/TETRA4、I_tetra4 = 1000

/TETRA10、I_tetra10 = 1000

完全に幾何学的非線形

I_smstr = 2、4

負の体積になる前のサイクルの要素形状を使用した微小ひずみに切り替え

Lagrange型全ひずみ。

I_smstr = 10、12

time=0.0における要素形状でのLagrange型全ひずみ

I_cpre - 定圧定式化フラグ
- I_cpre=1は、非圧縮性または疑似非圧縮性材料で体積ロックを阻止するために使用されます。この場合、応力テンソルは球状および偏差部分に分解されます。続いて圧力が一定に留まるよう、低減積分が球状部分に使用されます。
- I_cpre=2 は、弾塑性材料則に対してのみ使用できます。体積ロックを阻止するために、ポアソン係数を伴う追加の項がひずみに加えられます。材料内で尚も弾性で、そのため圧縮性である際は、ポアソン係数は小さく、材料が塑性となりそのため非圧縮性になると、ポアソン係数の項は増大して体積ロックを阻止します。詳細については、Radioss理論マニュアルをご参照ください。
I_npts - 積分点の数
- I_solid = 14および16について、推奨される値はI_npts = 222です。
テトラ要素
- I_solidフラグは４節点（/TETRA4）または１０節点（/TETRA10）四面体とは使用されません。
- ALEフレームワークによる4節点四面体は、I_tetra4 = 1000またはI_tetra4 = 3のみと互換性があります。
- I_tetra10 = 2を用いた2次四面体要素（/TETRA10）は、動的縮退を使用して、同じ元の積分スキームと線形四面体（/DT1/BRICKを使用）の時間ステップを持つ形状関数を使用します。メソッドの実装は独自のものです。
- I_tetra10=2の非互換性リストは次のとおりです。
  - AMS
  - 陰解法
  - Tetra10の中間節点が、/RBE3、/MPC、または/RLINKの従属節点として定義されている場合。
  - Tetra10の中間節点が、/RBODY、/RBE2のセカンダリ節点として単独で定義されている場合。
  - Tetra10の中間節点が、/BCS、/IMPACC、/IMPVELまたは/IMPDISで単独で定義されている場合。
  - 運動学的タイプの/INTER/TYPE2の場合：コーナー節点の1つがセカンダリ節点として定義されている場合、Spot_flag < 25および=30。
    （I_tetra10=2とともにインターフェースTYPE2のSpot_flag=27または28を使用することをお勧めします。）
I_frame - 要素座標系定式化フラグ
- 共回転定式化が使用されている場合（I_frame = 2）、応力テンソルは共回転座標系で計算されます。大きい回転が含まれる場合、この定式化はより精確ですが、計算コストは多くかかります。せん断変形が重要である弾性問題または粘弾性問題の場合に適用が推奨されます。
d - 数値減衰およびn - アワグラス粘性係数
- 数値減衰d_nは、I_solid = 24（HEPH）ソリッド要素のアワグラス応力計算に使用されます。数値減衰からのエネルギーは、時刻歴内部エネルギー出力内に含まれます。
- d_n=hであるI_solid = 24とI_solid = 1または2の間で結果を比較すると、数値減衰はI_solid = 24のほうがI_solid =1または2の際より $(2 / 3) \times 10^{- 3}$ 倍小さくなります。
- 数値Navier Stokes 粘性モデルは全ての材料側で利用可能です。粘性応力の出力はユーザー材料則とI_solid =1にのみ有効です（時刻歴出力では粘性応力が応力に加えられます）。
ポスト処理用出力
- ソリッド要素応力のポスト処理について、アニメーションは/ANIM/BRICK/TENSを、プロットファイルは/TH/BRICKをご参照ください。
- プロットファイルとアニメーションファイルでは、応力テンソルは共回転フレームに付加されます。
Solid to SPHプロパティ（Sol2SPH）
- Sol2SPHの使用時、1つのソリッドが接触、材料破壊基準または時間ステップ基準のせいで削除されると、ソリッド要素はSPH粒子に変換されます。
- アクティブになるSPH粒子の数は、パラメータNdirに依存します。粒子のプロパティはsphpart_IDパート番号を使用して計算されます。
- このオプションSol2SPHは、I_solid = 1、2または24、I_frame = 1または2とのみ互換性があります。
- インターフェース/INTER/TYPE7のオプション、I_del= -1または-2は、設定する必要があります。
- ソリッドとSPHパートについて両方の材料則を同じ種類にする必要があります。
フラグI_aleは、/ALE/MATでALEフレームワークを有効にする、または /EULER/MATでEULERフレームワークを有効にする従来の方法に取って代わるものです。

/PROP/TYPE14（SOLID）

フォーマット

定義

例

コメント