陽解法動解析

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

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

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

非線形性ソース

幾何学的非線形性

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

材料の非線形性

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

拘束と接触の非線形性

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

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

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

追従荷重

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

適用される荷重は、FLLWERバルクおよびサブケースエントリやPARAM,FLLWERエントリを使用して、追従荷重として示すことができます。
注: 追従荷重は現在、DLOAD/TLOAD#を介して指定された荷重、すべての圧力荷重、FORCE1FORCE2MOMENT1、およびMOMENT2でサポートされています。

陽解法有限要素解析法

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

中央差分法

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

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 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCytaiaahg gadaahaaWcbeqaaiaad6gaaaGccqGH9aqpcaWHMbWaaSbaaSqaaiaa dwgaaeqaaOWaaeWaaeaacaWH1bWaaWbaaSqabeaacaWGUbaaaOGaai ilaiaadshadaahaaWcbeqaaiaad6gaaaaakiaawIcacaGLPaaacqGH RaWkcaWHMbWaaSbaaSqaaiaadsgaaeqaaOWaaeWaaeaacaWH2bWaaW baaSqabeaacaWGUbaaaOGaaiilaiaadshadaahaaWcbeqaaiaad6ga aaaakiaawIcacaGLPaaacqGHRaWkcaWHMbWaaSbaaSqaaiaadogaae qaaOWaaeWaaeaacaWH1bWaaWbaaSqabeaacaWGUbaaaOGaaiilaiaa dshadaahaaWcbeqaaiaad6gaaaaakiaawIcacaGLPaaacqGHRaWkca WHMbWaaSbaaSqaaiaadIgaaeqaaOWaaeWaaeaacaWH1bWaaWbaaSqa beaacaWGUbaaaOGaaiilaiaadshadaahaaWcbeqaaiaad6gaaaaaki aawIcacaGLPaaacqGHsislcaWHMbWaaSbaaSqaaiaadMgaaeqaaOWa aeWaaeaacaWH1bWaaWbaaSqabeaacaWGUbaaaOGaaiilaiaadshada ahaaWcbeqaaiaad6gaaaaakiaawIcacaGLPaaaaaa@685E@

ここで、
M MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCytaaaa@36CD@
集中質量マトリックス
f e MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOzamaaBa aaleaacaWGLbaabeaaaaa@37FC@ f d MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOzamaaBa aaleaacaWGLbaabeaaaaa@37FC@ f c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOzamaaBa aaleaacaWGLbaabeaaaaa@37FC@ f h MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOzamaaBa aaleaacaWGLbaabeaaaaa@37FC@ 、および f i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOzamaaBa aaleaacaWGLbaabeaaaaa@37FC@
は、それぞれ外力、減衰力、接触力、アワグラス力、要素内力のベクトルです。
a n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCyyamaaCa aaleqabaGaamOBaaaaaaa@3801@
加速度ベクトルは釣り合い方程式から直接計算されます。

a n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCyyamaaCa aaleqabaGaamOBaaaaaaa@3801@ から、速度と変位のベクトルは、次のように更新できます:

v n + 1 2 = v n 1 2 + 1 2 ( t n + 1 2 t n 1 2 ) a n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCODamaaCa aaleqabaGaamOBaiabgUcaRmaaleaameaacaaIXaaabaGaaGOmaaaa aaGccqGH9aqpcaWH2bWaaWbaaSqabeaacaWGUbGaeyOeI0YaaSqaaW qaaiaaigdaaeaacaaIYaaaaaaakiabgUcaRmaalaaabaGaaGymaaqa aiaaikdaaaWaaeWaaeaacaWG0bWaaWbaaSqabeaacaWGUbGaey4kaS YaaSqaaWqaaiaaigdaaeaacaaIYaaaaaaakiabgkHiTiaadshadaah aaWcbeqaaiaad6gacqGHsisldaWcbaadbaGaaGymaaqaaiaaikdaaa aaaaGccaGLOaGaayzkaaGaaCyyamaaCaaaleqabaGaamOBaaaaaaa@506C@
d n+1 = d n +( t n+1 t n ) v n+ 1 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCizamaaCa aaleqabaGaamOBaiabgUcaRiaaigdaaaGccqGH9aqpcaWHKbWaaWba aSqabeaacaWGUbaaaOGaey4kaSYaaeWaaeaacaWG0bWaaWbaaSqabe aacaWGUbGaey4kaSIaaGymaaaakiabgkHiTiaadshadaahaaWcbeqa aiaad6gaaaaakiaawIcacaGLPaaacaWH2bWaaWbaaSqabeaacaWGUb Gaey4kaSYaaSqaaWqaaiaaigdaaeaacaaIYaaaaaaaaaa@4A98@

ここで、
t n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiDamaaCa aaleqabaGaamOBaaaaaaa@3810@
現在時刻
t n + 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiDamaaCa aaleqabaGaamOBaiabgUcaRiaaigdaaaaaaa@39AD@
次の時刻

次の時間増分が定義されます:

Δ t n = t n + 1 t n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyiLdqKaam iDamaaCaaaleqabaGaamOBaaaakiabg2da9iaadshadaahaaWcbeqa aiaad6gacqGHRaWkcaaIXaaaaOGaeyOeI0IaamiDamaaCaaaleqaba GaamOBaaaaaaa@414D@
Δ t n 1 = t n t n 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyiLdqKaam iDamaaCaaaleqabaGaamOBaiabgkHiTiaaigdaaaGccqGH9aqpcaWG 0bWaaWbaaSqabeaacaWGUbaaaOGaeyOeI0IaamiDamaaCaaaleqaba GaamOBaiabgkHiTiaaigdaaaaaaa@4300@

したがって:

v n + 1 2 = v n 1 2 + 1 2 ( Δ t n 1 + Δ t n ) a n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCODamaaCa aaleqabaGaamOBaiabgUcaRmaaleaameaacaaIXaaabaGaaGOmaaaa aaGccqGH9aqpcaWH2bWaaWbaaSqabeaacaWGUbGaeyOeI0YaaSqaaW qaaiaaigdaaeaacaaIYaaaaaaakiabgUcaRmaalaaabaGaaGymaaqa aiaaikdaaaWaaeWaaeaacqGHuoarcaWG0bWaaWbaaSqabeaacaWGUb GaeyOeI0IaaGymaaaakiabgUcaRiabgs5aejaadshadaahaaWcbeqa aiaad6gaaaaakiaawIcacaGLPaaacaWHHbWaaWbaaSqabeaacaWGUb aaaaaa@4FE0@
d n+1 = d n +Δ t n v n+ 1 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCizamaaCa aaleqabaGaamOBaiabgUcaRiaaigdaaaGccqGH9aqpcaWHKbWaaWba aSqabeaacaWGUbaaaOGaey4kaSIaeyiLdqKaamiDamaaCaaaleqaba GaamOBaaaakiaahAhadaahaaWcbeqaaiaad6gacqGHRaWkdaWcbaad baGaaGymaaqaaiaaikdaaaaaaaaa@45C9@

臨界時間ステップ

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

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

要素の時間ステップ

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

    時間ステップサイズは次の条件を満たす必要があります:

    Δ t 2 ω max MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyiLdqKaam iDaiabgsMiJoaalaaabaGaaGOmaaqaaiabeM8a3naaBaaaleaaciGG TbGaaiyyaiaacIhaaeqaaaaaaaa@3FA5@

    ここで、 ω max MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyYdC3aaS baaSqaaiGac2gacaGGHbGaaiiEaaqabaaaaa@3AC4@ はシステムの最大固有振動数を表します。

    ソリッド要素の場合は、臨界時間ステップサイズは次のように計算されます:

    Δ t e = l e Q + ( Q 2 + c 2 ) 1 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyiLdqKaam iDamaaBaaaleaacaWGLbaabeaakiabg2da9maalaaabaGaamiBamaa BaaaleaacaWGLbaabeaaaOqaaiaadgfacqGHRaWkdaqadaqaaiaadg fadaahaaWcbeqaaiaaikdaaaGccqGHRaWkcaWGJbWaaWbaaSqabeaa caaIYaaaaaGccaGLOaGaayzkaaWaaWbaaSqabeaadaWccaqaaiaaig daaeaacaaIYaaaaaaaaaaaaa@461B@

    ここで、
    c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yaaaa@36DF@
    断熱音速
    Q MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yaaaa@36DF@
    体積粘性係数 C 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ および C 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@

    Q = C 1 c + C 0 l e max ( 0 , ε ˙ k k ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyuaiabg2 da9iaadoeadaWgaaWcbaGaaGymaaqabaGccaWGJbGaey4kaSIaam4q amaaBaaaleaacaaIWaaabeaakiaadYgadaWgaaWcbaGaamyzaaqaba GcciGGTbGaaiyyaiaacIhadaqadaqaaiaaicdacaGGSaGaeyOeI0Ia fqyTduMbaiaadaWgaaWcbaGaam4AaiaadUgaaeqaaaGccaGLOaGaay zkaaaaaa@4999@

    ここで、
    C 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ および C 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@
    体積粘性係数は、デフォルト値がそれぞれ1.5と0.06である無次元定数です。
    l e MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiBamaaBa aaleaacaWGLbaabeaaaaa@37FE@
    要素の特性長さ。
    8節点六面体
    l e = V e A e max MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiBamaaBa aaleaacaWGLbaabeaakiabg2da9maalaaabaGaamOvamaaBaaaleaa caWGLbaabeaaaOqaaiaadgeadaWgaaWcbaGaamyzamaaBaaameaaci GGTbGaaiyyaiaacIhaaeqaaaWcbeaaaaaaaa@4001@
    10節点四面体
    l e = 1 ( B i j B i j ) 1 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiBamaaBa aaleaacaWGLbaabeaakiabg2da9maalaaabaGaaGymaaqaamaabmaa baGaamOqamaaBaaaleaacaWGPbGaamOAaaqabaGccaWGcbWaaSbaaS qaaiaadMgacaWGQbaabeaaaOGaayjkaiaawMcaamaaCaaaleqabaWa aSqaaWqaaiaaigdaaeaacaaIYaaaaaaaaaaaaa@42D7@
    6節点五面体
    l e = 1 ( B i j B i j ) 1 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiBamaaBa aaleaacaWGLbaabeaakiabg2da9maalaaabaGaaGymaaqaamaabmaa baGaamOqamaaBaaaleaacaWGPbGaamOAaaqabaGccaWGcbWaaSbaaS qaaiaadMgacaWGQbaabeaaaOGaayjkaiaawMcaamaaCaaaleqabaWa aSqaaWqaaiaaigdaaeaacaaIYaaaaaaaaaaaaa@42D7@
    4節点四面体
    l e = 3 V e A e max MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiBamaaBa aaleaacaWGLbaabeaakiabg2da9maalaaabaGaaG4maiaadAfadaWg aaWcbaGaamyzaaqabaaakeaacaWGbbWaaSbaaSqaaiaadwgadaWgaa adbaGaciyBaiaacggacaGG4baabeaaaSqabaaaaaaa@40BE@
    ここで、
    B i j MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOqamaaBa aaleaacaWGPbGaamOAaaqabaaaaa@38C7@
    形状関数の対称勾配
    V e MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOvamaaBa aaleaacaWGLbaabeaaaaa@37E8@
    六面体要素の体積
    A e max MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyqamaaBa aaleaacaWGLbWaaSbaaWqaaiGac2gacaGGHbGaaiiEaaqabaaaleqa aaaa@3ADF@
    六面体要素の全6フェイスのうちの最大面積
  • シェル要素

    シェル要素の場合は、時間ステップサイズは次の式によって決定されます:

    Δ t = L c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqqHuoarcaWG0bGaeyypa0ZaaSaaaeaacaWGmbaabaGaam4qaaaa aaa@3DE5@

    ここで、 c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGJbaaaa@39BF@ は、次の式で計算される音速です:

    c = E ρ ( 1 ν 2 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGJbGaeyypa0ZaaOaaaeaadaWcaaqaaiaadweaaeaacqaHbpGC daqadaqaaiaaigdacqGHsislcqaH9oGBdaahaaWcbeqaaiaaikdaaa aakiaawIcacaGLPaaaaaaaleqaaaaa@4356@

    ここで、
    E MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGfbaaaa@39A1@
    ヤング率
    ρ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqaHbpGCaaa@3A97@
    密度
    ν MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqaH9oGBaaa@3A8F@
    ポアソン比
    L MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGfbaaaa@39A1@
    特性長。四角形要素については次のように計算されます:
    L = A max ( L 1 , L 2 , L 3 , L 4 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGmbGaeyypa0ZaaSaaaeaacaWGbbaabaGaciyBaiaacggacaGG 4bWaaeWaaeaacaWGmbWaaSbaaSqaaiaaigdaaeqaaOGaaiilaiaadY eadaWgaaWcbaGaaGOmaaqabaGccaGGSaGaamitamaaBaaaleaacaaI ZaaabeaakiaacYcacaWGmbWaaSbaaSqaaiaaisdaaeqaaaGccaGLOa Gaayzkaaaaaaaa@48FF@
    ここで、
    A MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGbbaaaa@399D@
    面積
    L 1 , L 2 , L 3 , L 4 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGmbWaaSbaaSqaaiaaigdaaeqaaOGaaiilaiaadYeadaWgaaWc baGaaGOmaaqabaGccaGGSaGaamitamaaBaaaleaacaaIZaaabeaaki aacYcacaWGmbWaaSbaaSqaaiaaisdaaeqaaaaa@41EB@
    三角形要素の辺の長さ。次のように計算されます:
    L = 2 A max ( L 1 , L 2 , L 3 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGmbGaeyypa0ZaaSaaaeaacaaIYaGaey4fIOIaamyqaaqaaiGa c2gacaGGHbGaaiiEamaabmaabaGaamitamaaBaaaleaacaaIXaaabe aakiaacYcacaWGmbWaaSbaaSqaaiaaikdaaeqaaOGaaiilaiaadYea daWgaaWcbaGaaG4maaqabaaakiaawIcacaGLPaaaaaaaaa@4835@
    ここで、
    A MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGbbaaaa@399D@
    面積
    L 1 , L 2 , L 3 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGmbWaaSbaaSqaaiaaigdaaeqaaOGaaiilaiaadYeadaWgaaWc baGaaGOmaaqabaGccaGGSaGaamitamaaBaaaleaacaaIZaaabeaaaa a@3F76@
    要素の辺の長さ
  • スプリング要素

    スプリング要素(集中スプリング質量系)では、臨界時間ステップサイズを計算するための波動伝播速度がありません。

    節点質量が m 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyBamaaBa aaleaacaaIXaaabeaaaaa@37D0@ および m 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyBamaaBa aaleaacaaIXaaabeaaaaa@37D0@ 、剛性が k MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4Aaaaa@36E7@ のスプリングの自由振動の固有値問題は次のように表されます:

    [ k k k k ] [ u 1 u 2 ] ω [ m 1 m 2 ] [ u 1 u 2 ] = [ 0 0 ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaafa qabeGacaaabaGaam4AaaqaaiabgkHiTiaadUgaaeaacqGHsislcaWG RbaabaGaam4AaaaaaiaawUfacaGLDbaadaWadaqaauaabeqaceaaae aacaWG1bWaaSbaaSqaaiaaigdaaeqaaaGcbaGaamyDamaaBaaaleaa caaIYaaabeaaaaaakiaawUfacaGLDbaacqGHsislcqaHjpWDdaWada qaauaabeqaciaaaeaacaWGTbWaaSbaaSqaaiaaigdaaeqaaaGcbaaa baaabaGaamyBamaaBaaaleaacaaIYaaabeaaaaaakiaawUfacaGLDb aadaWadaqaauaabeqaceaaaeaacaWG1bWaaSbaaSqaaiaaigdaaeqa aaGcbaGaamyDamaaBaaaleaacaaIYaaabeaaaaaakiaawUfacaGLDb aacqGH9aqpdaWadaqaauaabeqaceaaaeaacaaIWaaabaGaaGimaaaa aiaawUfacaGLDbaaaaa@563B@

    特性方程式の行列式はゼロに等しいため、次のように最大固有値を求めることができます:

    | k ω 2 m 1 k k k ω 2 m 2 | = 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaqWaaeaafa qabeGacaaabaGaam4AaiabgkHiTiabeM8a3naaCaaaleqabaGaaGOm aaaakiaad2gadaWgaaWcbaGaaGymaaqabaaakeaacqGHsislcaWGRb aabaGaeyOeI0Iaam4AaaqaaiaadUgacqGHsislcqaHjpWDdaahaaWc beqaaiaaikdaaaGccaWGTbWaaSbaaSqaaiaaikdaaeqaaaaaaOGaay 5bSlaawIa7aiabg2da9iaaicdaaaa@4BA4@
    ここで、 ω max 2 = k ( m 1 + m 2 ) m 1 m 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyYdC3aa0 baaSqaaiGac2gacaGGHbGaaiiEaaqaaiaaikdaaaGccqGH9aqpdaWc aaqaaiaadUgadaqadaqaaiaad2gadaWgaaWcbaGaaGymaaqabaGccq GHRaWkcaWGTbWaaSbaaSqaaiaaikdaaeqaaaGccaGLOaGaayzkaaaa baGaamyBamaaBaaaleaacaaIXaaabeaakiabgwSixlaad2gadaWgaa WcbaGaaGOmaaqabaaaaaaa@49CA@ です。

    トラス要素の臨界時間ステップに基づいて:

    Δt l c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam iDaiabgsMiJoaalaaabaGaamiBaaqaaiaadogaaaaaaa@3BF4@ かつ ω max = 2c l MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyYdC3aaS baaSqaaiGac2gacaGGHbGaaiiEaaqabaGccqGH9aqpdaWcaaqaaiaa ikdacaWGJbaabaGaamiBaaaaaaa@3E79@ とすると、次のように記述できます:

    Δt 2 ω max MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam iDaiabgsMiJoaalaaabaGaaGOmaaqaaiabeM8a3naaBaaaleaaciGG TbGaaiyyaiaacIhaaeqaaaaaaaa@3FA4@

    実際のモード質量の半分を使用してスプリング質量を近似すると、次のような解が得られます:

    Δt=2 m 1 m 2 m 1 + m 2 1 k MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam iDaiabg2da9iaaikdadaGcaaqaamaalaaabaGaamyBamaaBaaaleaa caaIXaaabeaakiaad2gadaWgaaWcbaGaaGOmaaqabaaakeaacaWGTb WaaSbaaSqaaiaaigdaaeqaaOGaey4kaSIaamyBamaaBaaaleaacaaI YaaabeaaaaGcdaWcaaqaaiaaigdaaeaacaWGRbaaaaWcbeaaaaa@446E@

    よって、節点質量の観点では、次のように臨界時間ステップを表現できます:

    Δ t e = 2 m n1 m n2 m n1 + m n2 1 k MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam iDamaaBaaaleaacaWGLbaabeaakiabg2da9maakaaabaWaaSaaaeaa caaIYaGaamyBamaaBaaaleaacaWGUbGaaGymaaqabaGccaWGTbWaaS baaSqaaiaad6gacaaIYaaabeaaaOqaaiaad2gadaWgaaWcbaGaamOB aiaaigdaaeqaaOGaey4kaSIaamyBamaaBaaaleaacaWGUbGaaGOmaa qabaaaaOWaaSaaaeaacaaIXaaabaGaam4AaaaaaSqabaaaaa@495A@

    ここでは、減衰は考慮されません。減衰が定義されている場合は、時間ステップが次のようにスケールされます:

    Δ t e = 2 m n1 m n2 m n1 + m n2 1 k ( 1 ξ 2 ξ ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam iDamaaBaaaleaacaWGLbaabeaakiabg2da9maakaaabaWaaSaaaeaa caaIYaGaamyBamaaBaaaleaacaWGUbGaaGymaaqabaGccaWGTbWaaS baaSqaaiaad6gacaaIYaaabeaaaOqaaiaad2gadaWgaaWcbaGaamOB aiaaigdaaeqaaOGaey4kaSIaamyBamaaBaaaleaacaWGUbGaaGOmaa qabaaaaOWaaSaaaeaacaaIXaaabaGaam4AaaaaaSqabaGcdaqadaqa amaakaaabaGaaGymaiabgkHiTiabe67a4naaCaaaleqabaGaaGOmaa aaaeqaaOGaeyOeI0IaeqOVdGhacaGLOaGaayzkaaaaaa@520B@

    ここで、
    m n 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyBamaaBa aaleaacaWGUbGaaGymaaqabaaaaa@38C3@ および m n 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyBamaaBa aaleaacaWGUbGaaGymaaqabaaaaa@38C3@
    節点質量。
    k MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4Aaaaa@36E7@
    対応する自由度の剛性。
    ξ= c c cr = c 2mω MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqOVdGNaey ypa0ZaaSaaaeaacaWGJbaabaGaam4yamaaBaaaleaacaWGJbGaamOC aaqabaaaaOGaeyypa0ZaaSaaaeaacaWGJbaabaGaaGOmaiaad2gacq aHjpWDaaaaaa@422E@
    c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4Aaaaa@36E7@
    減衰係数(CBUSH要素の場合は、PBUSHバルクデータエントリのBiフィールドを介して定義されます)。

節点時間ステップ

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

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

Δ t n = 2 m n k n MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqqHuoarcaWG0bWaaSbaaSqaaiaad6gaaeqaaOGaeyypa0ZaaOaa aeaadaWcaaqaaiaaikdacaWGTbWaaSbaaSqaaiaad6gaaeqaaaGcba Gaam4AamaaBaaaleaacaWGUbaabeaaaaaabeaaaaa@426B@

ここで、
m n MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGTbWaaSbaaSqaaiaad6gaaeqaaaaa@3AE8@
節点質量
k n MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGTbWaaSbaaSqaaiaad6gaaeqaaaaa@3AE8@
節点剛性(要素剛性から計算されます)

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

各要素について、臨界時間ステップ Δ t e MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqqHuoarcaWG0bWaaSbaaSqaaiaadwgaaeqaaaaa@3C4C@ が最初に計算され、各節点は同じ時間ステップ Δ t e MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqqHuoarcaWG0bWaaSbaaSqaaiaadwgaaeqaaaaa@3C4C@ を持つものと想定されて、その後各節点について、次の式から節点剛性を推定できます。

Δ t e = 2 m e i k e i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqqHuoarcaWG0bWaaSbaaSqaaiaadwgaaeqaaOGaeyypa0ZaaOaa aeaadaWcaaqaaiaaikdacaWGTbWaaSbaaSqaaiaadwgadaWgaaadba GaamyAaaqabaaaleqaaaGcbaGaam4AamaaBaaaleaacaWGLbWaaSba aWqaaiaadMgaaeqaaaWcbeaaaaaabeaaaaa@449C@

ここで、
i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaale aacaWGPbaaaa@39C6@
要素のi番目の節点
m e i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGTbWaaSbaaSqaaiaadwgadaWgaaadbaGaamyAaaqabaaaleqa aaaa@3C05@
i番目の節点の節点質量
k e i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGTbWaaSbaaSqaaiaadwgadaWgaaadbaGaamyAaaqabaaaleqa aaaa@3C05@
この要素のi番目の節点の節点剛性

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

k e i = 2 m e i Δ t e 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGRbWaaSbaaSqaaiaadwgadaWgaaadbaGaamyAaaqabaaaleqa aOGaeyypa0ZaaSaaaeaacaaIYaGaamyBamaaBaaaleaacaWGLbWaaS baaWqaaiaadMgaaeqaaaWcbeaaaOqaaiabfs5aejaadshadaqhaaWc baGaamyzaaqaaiaaikdaaaaaaaaa@4549@

最後の節点剛性は次のとおりです:

k n = e k e i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGRbWaaSbaaSqaaiaad6gaaeqaaOGaeyypa0ZaaabuaeaacaWG RbWaaSbaaSqaaiaadwgadaWgaaadbaGaamyAaaqabaaaleqaaaqaai aadwgaaeqaniabggHiLdaaaa@4224@

k n MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGRbWaaSbaaSqaaiaad6gaaeqaaaaa@3AE6@ を使用して、節点臨界時間ステップ Δ t n MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqqHuoarcaWG0bWaaSbaaSqaaiaad6gaaeqaaaaa@3C55@ を計算できます。

質量スケーリング

  • 要素マススケーリング

    スケーリングされた要素臨界時間ステップ(DTFACによってスケーリング)がDTMINを下回る場合は、要素質量をスケーリングすることで、 Δ t e MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyiLdqKaam iDamaaBaaaleaacaWGLbaabeaaaaa@396D@ を増大させることができます。これが可能となる理由は、要素時間ステップ式には、材料密度( ρ )に依存する音速の項( c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGJbaaaa@39BF@ )が含まれているからです。

  • 節点マススケーリング

    スケーリングされた節点臨界時間ステップ(DTFACによってスケーリング)がDTMINを下回る場合は、節点質量 m n MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGTbWaaSbaaSqaaiaad6gaaeqaaaaa@3AE8@ をスケーリングすることで、 Δ t n MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqqHuoarcaWG0bWaaSbaaSqaaiaad6gaaeqaaaaa@3C55@ を増大させることができます。

  • 質量スケーリング・コントロール

    次の陽解法動的解析サブケースの質量のスケーリングは、MSCALEサブケース情報入力で制御できます。MSCALEが定義されていない場合、質量のスケーリングは直前の陽解法動的解析サブケースから継続されます。

アワグラス制御

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

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=FULL) 通常の要素(ISOPE=URI 通常の要素(ISOPE=AURI 通常の要素(ISOPE=SRI
CHEXA

(1次)

アワグラス制御は不要です。 アワグラス制御はデフォルトでオンになります。 1
MAT1/MATS1の場合
アワグラスタイプ2
MATHEの場合
アワグラスタイプ4
アワグラス制御はデフォルトでオンになります。1 デフォルトは次のとおりです:
MAT1/MATS1の場合
アワグラスタイプ2
MATHEの場合
アワグラスタイプ4
アワグラス制御は不要です。
CPENTA

(1次)

NA NA NA NA

シェル要素の場合

シェル要素については、次の2タイプのアワグラス制御のみが用意されています:
  • タイプ1(FlanaganとBelytschko - 粘性形式)
  • タイプ2(FlanaganとBelytschko - 剛性形式)タイプ2は、CQUAD4MAT1/MATS1材料のデフォルトのアワグラスタイプとして選択されます。
アワグラス制御(シェル要素ベース)
要素 Belytschko-Tsay(ISOPE=BT 完全投影を使用したBelytschko-Wong-Chiang(ISOPE=BWC
CQUAD4 アワグラス制御はデフォルトでオンになります。 1
MAT1/MATS1の場合
アワグラスタイプ2
アワグラス制御はデフォルトでオンになります。 1
MAT1/MATS1の場合
アワグラスタイプ2

材料

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

ソリッドとシェル:Flanagan-Belytschko粘性形式

タイプ2

ソリッド: Puso拡張想定ひずみ剛性形式

シェル: Flanagan-Belytschko剛性形式

タイプ4

ソリッド:Reeseアワグラス制御

シェル:タイプ4はシェルでのサポートなし

MAT1/MAT2/MAT8/MATS1 使用可能 2 デフォルト 6 NA
MATHE NA 使用可能 2 デフォルト 6
MATVE NA 使用可能 2 デフォルト 6

アダプティブ動的緩和

動的緩和を使用すれば、陽解法動解析を使用して動的振動を回避することにより、静的または準静的問題を解析できます。陰解法解析と比較して、非線形性が高いケース(複雑な接触が多いなど)でより効率的で堅牢です。代表的な適用例として、電話構造の3点曲げシミュレーションやシートメタル成形におけるスプリングバックシミュレーションが挙げられます。

少なくとも1つの入力が必要な従来型の動的緩和と違って、OptiStructは、入力パラメータが必要ないDYRELエントリを介したアダプティブ動的緩和をサポートしています。減衰係数は、系の最も高い固有振動数に基づいて自動的に決定されます。

Material Failure Criterion

材料破壊基準を定義するには、MATFバルクデータエントリまたはMATS1バルクデータエントリ(損傷開始 / 進展基準の場合のみ)を使用します。材料の破壊は、荷重条件に大きく影響され、その結果として応力状態も大きな影響を受けます。したがって、使用可能ないくつかの基準では、応力三軸度の概念およびオプションでLodeパラメータを参照して、荷重条件(単軸引張、純せん断、平面ひずみなど)が記述されます。

塑性と応力状態に基づいて破壊基準を記述するには、値の応力三軸度( η )とLodeパラメータ( ξ )が必要です。シェルの場合のみ、応力三軸度が必要となります。

応力三軸度

応力三軸度( η )は、圧縮荷重と引張荷重を区別するために使用され、応力テンソルのトレースに依存しています。図 1に示すように、応力三軸度によって、静水圧軸上の応力状態の位置を決定できます。
1. 静水圧軸と偏差面上の応力状態の説明


これは次のように計算されます:
η = 1 3 σ x x + σ y y + σ z z σ V M MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH3oaAcq GH9aqpdaWcaaqaamaalaaabaGaaGymaaqaaiaaiodaaaWaaeWaceaa cqaHdpWCdaWgaaWcbaGaamiEaiaadIhaaeqaaOGaey4kaSIaeq4Wdm 3aaSbaaSqaaiaadMhacaWG5baabeaakiabgUcaRiabeo8aZnaaBaaa leaacaWG6bGaamOEaaqabaaakiaawIcacaGLPaaaaeaacqaHdpWCda WgaaWcbaGaamOvaiaad2eaaeqaaaaaaaa@4E78@
ここで、 σ V M MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaHdpWCda WgaaWcbaGaamOvaiaad2eaaeqaaaaa@3B00@ はフォンミーゼス相当応力です。
応力三軸度の値の変動範囲は、ソリッドの場合、 MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyOeI0Iaey OhIukaaa@3852@ + MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaey4kaSIaey OhIukaaa@3847@ であり(実際には-1~1の範囲に制限されています)、シェルの場合、-2/3~2/3です。
1. いくつかの一般的な応力状態の応力三軸度の値
荷重条件 ソリッド要素 シェル要素
拘束圧縮 -1
2軸圧縮 -2/3 -2/3
単軸圧縮 -1/3 -1/3
純せん断 0.0 0.0
単軸引張 1/3 1/3
平面ひずみ 0.5751 0.5751
2軸引張 2/3 2/3
拘束引張 1

Lode角

3D荷重条件を記述する上で、もう1つの重要な量は、次の式によって得られるLode角( θ )です:

cos 3 θ = 27 2 J 3 σ V M 3 MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaaciGGJbGaai 4BaiaacohadaqadaqaaiaaiodacqaH4oqCaiaawIcacaGLPaaacqGH 9aqpdaWcaaqaaiaaikdacaaI3aaabaGaaGOmaaaadaWcaaqaaiaadQ eadaWgaaWcbaGaaG4maaqabaaakeaacqaHdpWCdaqhaaWcbaGaamOv aiaad2eaaeaacaaIZaaaaaaaaaa@47AE@
ここで、 J 3 MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGkbWaaS baaSqaaiaaiodaaeqaaaaa@391C@ は第3偏差不変量です。

Lode角によって、偏差断面内の応力状態の位置が決定されます。その値は、0(引張の場合)と π / 3 MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaHapaCca GGVaGaaG4maaaa@3A91@ (圧縮の場合)の間になります。
2. Lode角の値に応じた偏差面上の応力状態位置


せん断ひずみと平面ひずみの条件では、 π / 6 MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaHapaCca GGVaGaaGOnaaaa@3A94@ というLode角の値が使用されます。

平面応力と仮定すると(シェル要素の場合)、Lode角と応力三軸度はリンクされているため、次のようにこれらの一方を使用して他方を復元できます:

cos 3 θ = 27 2 η η 2 1 3 MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaaciGGJbGaai 4BaiaacohadaqadaqaaiaaiodacqaH4oqCaiaawIcacaGLPaaacqGH 9aqpcqGHsisldaWcaaqaaiaaikdacaaI3aaabaGaaGOmaaaacqaH3o aAdaqadaqaaiabeE7aOnaaCaaaleqabaGaaGOmaaaakiabgkHiTmaa laaabaGaaGymaaqaaiaaiodaaaaacaGLOaGaayzkaaaaaa@4AB8@

ラジアンよりも正規化された値の方がはるかに扱いやすいため、Lode角は通常、次の式によって得られるLodeパラメータ ξ に置き換えられます:

ξ = 1 6 π θ MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH+oaEcq GH9aqpcaaIXaGaeyOeI0YaaSaaaeaacaaI2aaabaGaeqiWdahaaiab eI7aXbaa@4018@

Lodeパラメータの値は次のとおりです:
  • -1.0(圧縮状態)
  • 純せん断状態または平面ひずみ状態
  • 引張状態

サポートされている破壊基準

現在、陽解法動解析では、BIQUAD、TSTRN、表形式破壊基準、INIEVOの4つの破壊基準がサポートされています。
BIQUAD
BIQUAD基準は、応力三軸度に基づいた破壊基準であり、主に延性金属に使用されます。この二重の二次曲線形状は、応力三軸度 η に関する破壊時の塑性ひずみ ε p f MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTdu2aa0 baaSqaaiaadchaaeaacaWGMbaaaaaa@39A8@ の進展を表します(下図を参照)。
3. BIQUAD基準の応力三軸度に関する破壊時の塑性ひずみの進展


この基準には、5つのパラメータc1、c2、c3、c4、c5が必要であり、これらはそれぞれMATFバルクデータエントリ内のV1V2V3V4V5の値に対応しています。これら5つの値は、次の5種類の応力状態の破壊時塑性ひずみに対応しています:
  1. 単軸圧縮
  2. 純せん断
  3. 単軸引張
  4. 平面ひずみ
  5. 2軸引張
注: 高い応力三軸度における放物線の計算は、c4が常に最小値になるように行われます。
シェル要素の場合は、高いひずみ速度で発生しているひずみの局所化とネッキングは、正しく検出されない可能性があります。その理由は、板厚変化が純粋に数値的であるためです。したがって、同じサイズのソリッド要素と比べて、破壊が遅くなる可能性があります。このことを回避するために、シェルに対して追加の曲線を定義できます(下図の青い曲線を参照)。そのためには、INSTパラメータ(V6)を使用して、高い応力三軸度の放物線計算でc4を置き換えます。
4. シェルの高い応力三軸度における追加の破壊放物線(青色)


c1、c2、c3、c4、c5のパラメータすべてを特定するのに十分な実験データが得られない場合は、BIQUAD基準に対して材料セレクター入力も使用できます。上記のリストで選択されたキーワードMATERの値に応じて、次のようにc3の値に基づいてc1、c2、c4、c5のパラメータが自動的に計算されます。
c 1 = r 1 c 3 MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yaiaaig dacqGH9aqpcaWGYbGaaGymaiabgwSixlaadogacaaIZaaaaa@3E3E@ c 2 = r 2 c 3 MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yaiaaik dacqGH9aqpcaWGYbGaaGOmaiabgwSixlaadogacaaIZaaaaa@3E40@ c 3 = c 3 MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yaiaaio dacqGH9aqpcaWGJbGaaG4maaaa@3A44@ c 4 = r 4 c 3 MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yaiaais dacqGH9aqpcaWGYbGaaGinaiabgwSixlaadogacaaIZaaaaa@3E44@ c 5 = r 5 c 3 MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yaiaaiw dacqGH9aqpcaWGYbGaaGynaiabgwSixlaadogacaaIZaaaaa@3E46@
したがって、BIQUAD基準の材料入力を使用する際に、期待されるパラメータはc3の値のみです。ただし、c3の値が指定されていない場合は、c3のデフォルト値が自動的に設定されます。
2. MATERキーワードの自動パラメータ設定
キーワード c3(デフォルト) r1 r2 r4 r5
MILD 0.60 3.5 1.6 0.6 1.5
HSS 0.50 4.3 1.4 0.6 1.6
UHSS 0.12 5.2 3.1 0.8 3.5
AA5182 0.30 5.0 1.0 0.4 0.8
AA6082 0.17 7.8 3.5 0.6 2.8
PA6GF30 0.10 3.6 0.6 0.5 0.6
PP T40 0.11 10.0 2.7 0.6 0.7
時間ステップごとに、応力三軸度と放物線に従って、破壊時の塑性ひずみ ε p f ( η ) MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTdu2aa0 baaSqaaiaadchaaeaacaWGMbaaaOGaaiikaiabeE7aOjaacMcaaaa@3CB7@ が推測されます。これにより、応力状態の履歴を考慮して損傷変数の増大が可能になります:
D = t = 0 Δ ε p ε p f ( η ) MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiraiabg2 da9maaqahabaWaaSaaaeaacqGHuoarcqaH1oqzdaWgaaWcbaGaamiC aaqabaaakeaacqaH1oqzdaqhaaWcbaGaamiCaaqaaiaadAgaaaGcca GGOaGaeq4TdGMaaiykaaaaaSqaaiaadshacqGH9aqpcaaIWaaabaGa eyOhIukaniabggHiLdaaaa@493B@
TSTRN
TSTRN破壊基準はひずみに基づいた損傷モデルであり、完全連成であると想定されます(DAMAGEキーワードがアクティブ、かつ D C = 0 MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaWGdbaabeaakiabg2da9iaaicdaaaa@397B@ )。ただし、これを、破壊基準として使用することも、純粋な出力損傷変数として使用することもできます。ここでは、引張荷重条件において、ひずみの開始値と終了値の2点間の損傷変数の線形進展が考慮されます( η > 0 MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4TdGMaey Opa4JaaGimaaaa@3962@ ):
D = ε ε s t a r t ε e n d ε s t a r t MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiraiabg2 da9maalaaabaGaeqyTduMaeyOeI0IaeqyTdu2aaSbaaSqaaiaadoha caWG0bGaamyyaiaadkhacaWG0baabeaaaOqaaiabew7aLnaaBaaale aacaWGLbGaamOBaiaadsgaaeqaaOGaeyOeI0IaeqyTdu2aaSbaaSqa aiaadohacaWG0bGaamyyaiaadkhacaWG0baabeaaaaaaaa@4D35@
ε s t a r t MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTdu2aaS baaSqaaiaadohacaWG0bGaamyyaiaadkhacaWG0baabeaaaaa@3C8E@ ε e n d MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTdu2aaS baaSqaaiaadwgacaWGUbGaamizaaqabaaaaa@3A8D@ の値のペアがカード内で必要です。V1V2の値は、フォンミーゼス相当ひずみの開始値と終了値に対応しています。フォンミーゼス相当ひずみは次のように計算されます:
ε = 2 3 ε x x 2 ' + ε y y 2 ' + ε z z 2 ' + 2 ε x y 2 + 2 ε y z 2 + 2 ε z x 2 MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTduMaey ypa0ZaaOaaaeaadaWcaaqaaiaaikdaaeaacaaIZaaaamaabmaabaGa eqyTdu2aa0baaSqaaiaadIhacaWG4baabaGaaGOmaaaakmaaCaaale qabaGaai4jaaaakiabgUcaRiabew7aLnaaDaaaleaacaWG5bGaamyE aaqaaiaaikdaaaGcdaahaaWcbeqaaiaacEcaaaGccqGHRaWkcqaH1o qzdaqhaaWcbaGaamOEaiaadQhaaeaacaaIYaaaaOWaaWbaaSqabeaa caGGNaaaaOGaey4kaSIaaGOmaiabew7aLnaaDaaaleaacaWG4bGaam yEaaqaaiaaikdaaaGccqGHRaWkcaaIYaGaeqyTdu2aa0baaSqaaiaa dMhacaWG6baabaGaaGOmaaaakiabgUcaRiaaikdacqaH1oqzdaqhaa WcbaGaamOEaiaadIhaaeaacaaIYaaaaaGccaGLOaGaayzkaaaaleqa aaaa@6095@
ここで、 ε ' MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTdu2aaW baaSqabeaacaGGNaaaaaaa@3872@ は偏差ひずみテンソルです。
V3V4の値が指定されている場合、これらは最大主ひずみの開始値と終了値に対応しています。
ε = max ( ε 1 , ε 2 , ε 3 ) > 0 MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTduMaey ypa0JaciyBaiaacggacaGG4bGaaiikaiabew7aLnaaBaaaleaacaaI XaaabeaakiaacYcacqaH1oqzdaWgaaWcbaGaaGOmaaqabaGccaGGSa GaeqyTdu2aaSbaaSqaaiaaiodaaeqaaOGaaiykaiabg6da+iaaicda aaa@47BB@
注: V1/V2V3/V4の両方のペアが指定されている場合は、V3V4の値が常に優先されます。
表形式破壊基準
TAB破壊基準を使用すると、塑性ひずみに基づいた表形式基準を記述する上で、最大限の自由度が得られます。EPS_TIDによって定義されるTABLEMDエントリは、応力三軸度に関する、またソリッド要素の場合はオプションでLodeパラメータ ξ に関する、破壊時の塑性ひずみ ε p f MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTdu2aa0 baaSqaaiaadchaaeaacaWGMbaaaaaa@39A8@ の進展を示すマップを表します。(図 5を参照)。
5. 応力三軸度とLodeパラメータに関する破壊時の塑性ひずみの進展を表す表形式破壊基準のマップ


ソリッド要素の場合、考えられるすべての値のペア ( η , ξ ) MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikaiabeE 7aOjaacYcacqaH+oaEcaGGPaaaaa@3B6C@ を含むマップ全体が考慮されます。ただし、シェルの場合は、平面応力条件のために、応力三軸度とLodeパラメータはリンクされています。したがって、平面応力(図 5の青い線)のみが考慮されます。

V1の値は、マップ全体におけるすばやい増減を可能にするスケールファクターです。

損傷変数の進展は、V2の値で定義されたパラメータを使用した特定の式によって得られます:
D = t = 0 Δ ε p ε p f ( η , ξ ) n D 1 1 n MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiraiabg2 da9maaqahabaWaaSaaaeaacqGHuoarcqaH1oqzdaWgaaWcbaGaamiC aaqabaaakeaacqaH1oqzdaqhaaWcbaGaamiCaaqaaiaadAgaaaGcca GGOaGaeq4TdGMaaiilaiabe67a4jaacMcaaaaaleaacaWG0bGaeyyp a0JaaGimaaqaaiabg6HiLcqdcqGHris5aOGaeyyXICTaamOBaiabgw SixlaadseadaahaaWcbeqaamaabmaabaGaaGymaiabgkHiTmaalaaa baGaaGymaaqaaiaad6gaaaaacaGLOaGaayzkaaaaaaaa@5724@
したがって、それ自体の現行値も含め、損傷変数の進展では、応力状態の履歴だけでなく、損傷履歴も考慮されています。指数 n MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBaaaa@36E6@ を使用すると、塑性ひずみに関する損傷進展の形状を間接的に変更できます(図 6を参照)。 n MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBaaaa@36E6@ 指数パラメータが増大すると、図のように応力軟化効果が遅れる傾向にあります。
6. nパラメータが損傷と塑性ひずみの進展の関係に与える影響(左側)、およびnパラメータが単一要素の単軸引張挙動に与える影響(右側)




1行目のパラメータ(EPS_TIDV1V2)のみを定義するTAB基準を使用できます。この場合、他の使用可能な基準と同様に、DAMAGEを使用して要素の削除をアクティブにし、臨界損傷DCに対して応力軟化が低数値となる開始点を選択して、EXPにより応力軟化の形状を選択できます。

TAB基準を使用したもう1つの応力軟化アプローチは、ネッキング制御アプローチと呼ばれます。

この新しいアプローチを使用するには、2行目の最初の2つのパラメータINST_TIDV6を定義する必要があります。INST_TIDでは、応力三軸度およびオプションでLodeパラメータに関して、ネッキング不安定性、さらにはひずみの局所化が始まる、塑性ひずみ値( ε p l o c MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTdu2aa0 baaSqaaiaadchaaeaacaWGSbGaam4Baiaadogaaaaaaa@3B8A@ )の進展を示すマップを定義するTABLEMDエントリのIDを定義します。これは、主に高い応力三軸度において定義された不安定性限界曲線またはマップであり(図 3でBIQUAD基準について前述したものと同様です)、有効となるためには破壊曲線 / マップより低い値である必要があります。これは、ソリッドまたはシェルで使用できます。

このINST第2マップを使用すると、ネッキングトリガリング変数と呼ばれる新しい変数( f MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzaaaa@36DF@ )の進展を計算できます。この進展は、損傷変数の進展とよく似ています:

f = t = 0 Δ ε p ε p l o c ( η , ξ ) n f 1 1 n MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzaiabg2 da9maaqahabaWaaSaaaeaacqGHuoarcqaH1oqzdaWgaaWcbaGaamiC aaqabaaakeaacqaH1oqzdaqhaaWcbaGaamiCaaqaaiaadYgacaWGVb Gaam4yaaaakiaacIcacqaH3oaAcaGGSaGaeqOVdGNaaiykaaaaaSqa aiaadshacqGH9aqpcaaIWaaabaGaeyOhIukaniabggHiLdGccqGHfl Y1caWGUbGaeyyXICTaamOzamaaCaaaleqabaWaaeWaaeaacaaIXaGa eyOeI0YaaSaaaeaacaaIXaaabaGaamOBaaaaaiaawIcacaGLPaaaaa aaaa@594A@

この変数の値が1に達すると、応力軟化がトリガーされます(MATFバルクデータエントリのコメント12で定義されます)。ただし、MATFエントリで定数値 D C MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaWGdbaabeaaaaa@37B1@ を使用するのではなく、パラメータ D C MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaWGdbaabeaaaaa@37B1@ は積分点になります。したがって、 D C MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaWGdbaabeaaaaa@37B1@ は、要素応力状態の履歴に応じて、要素間で大きく異なる可能性があります。

このため、INST_TIDが使用されている場合、 D C MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaWGdbaabeaaaaa@37B1@ の値は、 f MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzaaaa@36DF@ の値が1以上になった瞬間に損傷変数 D MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiraaaa@36BD@ が取る値に一致します。すなわち、 D C MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaWGdbaabeaaaaa@37B1@ は、最初にネッキング基準に達したときの D MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiraaaa@36BD@ の値です。その後、 D C MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaWGdbaabeaaaaa@37B1@ は、シミュレーション終了まで変更されないままとなります。

D = Δ D f = Δ f D C = 1 while f < 1 D when f 1 σ = σ e f f 1 D D c 1 D c exp MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaacaWGeb Gaeyypa0Zaa8qaaeaacqGHuoarcaWGebaaleqabeqdcqGHRiI8aaGc baGaamOzaiabg2da9maapeaabaGaeyiLdqKaamOzaaWcbeqab0Gaey 4kIipaaOqaaiaadseadaWgaaWcbaGaam4qaaqabaGccqGH9aqpdaGa baqaauaabeqaceaaaeaafaqabeqadaaabaGaaGymaaqaaiaabEhaca qGObGaaeyAaiaabYgacaqGLbaabaGaamOzaiabgYda8iaaigdaaaaa baqbaeqabeWaaaqaaiaadseaaeaacaqG3bGaaeiAaiaabwgacaqGUb aabaGaamOzaiabgwMiZkaaigdaaaaaaaGaay5EaaaabaGaeq4WdmNa eyypa0Jaeq4Wdm3aaSbaaSqaaiaadwgacaWGMbGaamOzaaqabaGcda qadaqaaiaaigdacqGHsisldaqadaqaamaalaaabaGaamiraiabgkHi TiaadseadaWgaaWcbaGaam4yaaqabaaakeaacaaIXaGaeyOeI0Iaam iramaaBaaaleaacaWGJbaabeaaaaaakiaawIcacaGLPaaadaahaaWc beqaaiGacwgacaGG4bGaaiiCaaaaaOGaayjkaiaawMcaaaaaaa@6CE4@
D C MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaWGdbaabeaaaaa@37B1@ パラメータとは異なり、指数(EXP)は、すべての要素にわたって一定のパラメータです。

このネッキング制御アプローチでは、幅広い応力状態について予測精度を高めることができますが、特にネッキングが発生しやすい高い応力三軸度においては、不安定性マップを定義する必要があります。

最後に、パラメータV7V8の値は、以下で定義している要素サイズスケーリングの応力三軸度境界です。この値ペアが定義されている場合は、サイズのスケーリングは次の条件が満たされる場合にのみ行われます:

V 7 < η < V 8 MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOvaiaaiE dacqGH8aapcqaH3oaAcqGH8aapcaWGwbGaaGioaaaa@3CE1@

損傷の開始と進展(INIEVO)
INIEVO破壊基準は非常に具体的であり、DMGINIバルクデータエントリと、オプションでDMGEVOバルクデータエントリを使用して、破壊アプローチを定義することが可能になります。

DMGINIバルクデータエントリでは、DUCTILE基準のみを使用できます。DMGEVOバルクデータエントリでは、DISPとENERGYの進展のみを使用できます。

この基準は、次の2つの方法で定義できます:
  1. MATS1バルクデータエントリのDAMAGE継続行。この方法は、陰解法動解析と陽解法動解析の両方でサポートされています。
  2. MATFバルクデータエントリのCRI=INIEVO。この方法は、陽解法動解析でのみサポートされています。
注: INIEVOでは、ひずみ速度依存性と要素サイズ依存性は使用できません。

問題の設定

入力

  • アクティブ化:

    非線形陽解法サブケースは、ANALYSIS=NLEXPLによって指定できます。TTERMサブケースエントリは、終了時刻を定義するために必要です。また、対応するTSTEPEバルクデータエントリを指しているTSTEPEサブケースエントリを、非線形陽解法解析で使用することもできます。TSTEPEサブケースエントリが定義されていない場合は、ANALYSIS=NLEXPLTTERMと組み合わせて使用する必要があります。このサブケースエントリが定義されている場合は、TTERMTSTEPEを併用するだけで、陽解法非線形サブケースを指定できます。非線形陽解法解析は常に大変位解析です。

  • 初期条件:

    初期条件を定義するには、ICサブケースエントリをTICバルクデータエントリと組み合わせて使用します。初期温度フィールドは、TEMP(INIT)を使用して定義することができます。 TEMP(INIT)は、参照される温度フィールドを使用して、対応する MATS1エントリの初期材料データのTABLEMDエントリを検索します。

  • 荷重:

    荷重を定義するには、LOADDLOAD、およびTLOAD#バルクデータエントリを使用します。これらは、DLOADサブケースエントリを使用してサブケース内で参照される必要があります。LOADサブケースエントリまたはTLOAD#バルクエントリを介した参照では、荷重に対してFORCEFORCE1FORCE2MOMENTMOMENT1MOMENT2PLOAD2PLOAD4GRAVACCEL2、およびSPCDエントリのみがサポートされています。

  • 境界条件:

    境界条件は、対応するSPCサブケースエントリで参照されているSPCバルクデータを介して適用できます。MPCは現在サポートされていません。

  • サポートされている要素:
    ソリッド要素
    4節点CTETRA、10節点CTETRA、8節点CHEXA、6節点CPENTAの各要素がサポートされています。
    シェル要素
    CTRIA3CQUAD4がサポートされています。
    1次元要素
    CBUSHCBEAMおよびCBAR要素がサポートされています。
    現在、陽解法解析ではCBAR/CBEAM要素に対しては、Belytschko-Schwer Beamの定式化のみがサポートされています。
    質量要素
    CONM2がサポートされています。
    注:
    • シェル要素の要素またはプロパティのオフセットは、陽解法解析でサポートされています。
    • CBUSH要素の場合、PBUSH定義のMiフィールドは質量と慣性の計算に使用されます。詳細についてはリファレンスガイド内のPBUSHをご参照ください。
    • CBEAMCBAR要素の場合、
      • PBEAM/PBAR上の継続行は、陽解法解析ではサポートされません。
      • ピンフラグ(PAPB)は陽解法解析でサポートされています。
  • サポートされている材料:

    MAT1MAT2MAT8MATS1MATHEMATVE材料はサポートされています。MATVEエントリをMATHEエントリ下で定義する必要があります。

  • 積分スキーム:

    陽解法解析では、要素積分スキームは、PSOLIDPLSOLIDPLSOLIDPSHELLPCOMPPCOMPGPCOMPPエントリの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

NLDEBUG, CONT2TIENLDEBUG, RMNLMAT は、問題点を解決する目的で、特定の方法でモデルを単純化するために利用できます。

出力

一般的な出力エントリ(DISPLACEMENTVELOCITYおよび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要素力、複合材の応力、複合材のひずみ、複合破壊指数の等高線が出力されます。
MONVOLバルクデータエントリでモニターボリュームを定義すると、デフォルトで以下の出力結果が利用できます。圧力、温度、体積、面積、質量、内部エネルギー、質量流量、ベント面積、漏れ質量。
_expl.mvw
このセッションファイルは、対応する_expl.h3dファイルを自動的に読み込んで、_expl.h3dファイル内の結果出力をプロットできるようにします。
_s<ID>_e.expl
内部エネルギー、弾性接触エネルギー、塑性接触エネルギー、運動エネルギー、アワグラスエネルギー、および塑性散逸エネルギーの曲線が出力されます。
_expl_energy.mvw
このセッションファイルは、対応する_s<ID>_e.explファイルを自動的に読み込んで、さまざまなエネルギー出力をプロットできるようにします。
.out
陽解法では、.outファイルに、時間サイクル情報(PARAM,NOUTCYCに基づく)、現在時刻、現在の時間ステップ、最大ひずみエネルギー、この情報が出力される要素ID、運動エネルギー、接触仕事、合計エネルギー、最大貫通、この最大貫通に関連付けられた節点ID、最大垂直仕事、この最大垂直仕事に関連付けられた節点ID、質量変化率が含まれます。質量変化率は、マススケーリング後のスケーリングされた質量の変化に関する情報であり、これは(現在の質量 - 元の質量)/(元の質量)で算出されます。
_expl.cntf
メインのサーフェスの接触力出力結果を含んだASCIIファイルで、CONTF I/OオプションエントリでOPTIフォーマットが指定されている場合にアクティブになります。出力には、法線/接線力、大きさ、接触面積が含まれます。この出力は、陽解法時間ステップ毎に用意されます。
このファイルの出力間隔は、NLOUTエントリのNINTフィールドを使用して制御できます。
_TH.h5
陽解法動解析の時刻歴出力は、_TH.h5ファイル HDF5フォーマットファイルで入手できます。状況によっては、高い出力頻度で出力するための結果(エネルギーなど)のサブセットが必要です。ただし、NLOUTの出力頻度を上げると、すべての結果に影響してファイルサイズが大きくなり、望ましくありません。時刻歴出力は、このようなケースで役に立つ効率的なソリューションです。
サポートされている結果の詳細については、THISTバルクデータエントリをご参照ください。
3. 陽解法動解析の簡単なまとめ
非線形陽解法解析 サブケースまたは入出力 バルクデータ コメント
アクティブ化:
サブケースタイプ ANALYSIS=NLEXPL(オプション) NA TSTEPEが指定されていない場合は、ANALYSIS=NLEXPLは必須です。
非線形陽解法のアクティブ化 TTERM(必須)

TSTEPE(オプション)

TSTEPE(オプション) TSTEPEが指定されていない場合は、ANALYSIS=NLEXPLは必須です。
荷重:
節点荷重 LOADDLOAD サブケース内のLOADが使用されている場合:

FORCEFORCE1FORCE2MOMENTMOMENT1、およびMOMENT2

サブケース内のDLOADが使用されている場合:

TLOAD1またはTLOAD2

DLOADを使用して、複数のTLOADiデータを結合できます。

節点荷重では、TLOADiデータのEXCITEIDに設定できる値は、FORCEFORCE1FORCE2MOMENTMOMENT1、およびMOMENT2です。

このケースでTLOADiデータのTYPEフィールドは、0またはLOADに設定できます。
サーフェス荷重 LOADDLOAD サブケース内のLOADが使用されている場合:

PLOAD2およびPLOAD4

サブケース内のDLOADが使用されている場合:

TLOAD1またはTLOAD2

DLOADを使用して、複数のTLOADiデータを結合できます。

サーフェス荷重では、TLOADiデータのEXCITEIDは、PLOAD1およびPLOAD4に設定できます。

このケースでTLOADiデータのTYPEフィールドは、0またはLOADに設定できます。
ボディ荷重 LOADDLOAD サブケース内のLOADが使用されている場合:

GRAVおよびACCEL2

サブケース内のDLOADが使用されている場合:

TLOAD1またはTLOAD2

DLOADを使用して、複数のTLOADiデータを結合できます。

物体力荷重では、TLOADiデータのEXCITEIDは、GRAVおよびACCEL2に設定できます。

このケースでTLOADiデータのTYPEフィールドは、0またはLOADに設定できます。
強制変位、強制速度、強制加速度 LOADDLOAD サブケース内のLOADが使用されている場合:

SPCDまたはSPCDを使用した強制変位、強制速度、または強制加速度。

サブケース内のDLOADが使用されている場合:

TLOAD1またはTLOAD2

DLOADを使用して、複数のTLOADiデータを結合できます。

強制荷重では、TLOADiデータのEXCITEIDは、SPCまたはSPCDに設定できます。

TLOADiデータのTYPEフィールドは、以下の値に設定できます:
1またはDISP
強制変位の場合
2またはVELO
強制速度の場合
3またはACCE
強制加速度の場合
追従荷重 FLLWER FLLWER

PARAM,FLLWER

陰解法非線形解析と同様に、荷重として追従荷重を選択できます。

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

境界条件:
単点拘束 SPC SPC
初期条件:
初期変位 TIC IC
初速度 TIC IC
時間ステップコントロール:
基本的な時間コントロール TSTEPE TSTEPE TSTEPEエントリのTYPEフィールドを使用して、要素時間ステップコントロールまたは節点時間ステップコントロールのどちらかを選択します。

DTMINフィールドでは、最小時間ステップを定義できます。時間ステップがこれより小さい場合、節点 / 要素マススケーリングがアクティブ化されます。

DTFACフィールドでは、安定した時間増分のためのスケールファクターを定義できます。

質量要素:
質量要素のサポート CONM2がサポートされています。
構造要素:
サポートされている構造要素 NA 1次元要素:CBUSHCBEAMおよびCBARがサポートされています。
シェル要素
CTRIA3およびCQUAD4要素。
ソリッド要素
4節点CTETRA、10節点CTETRA、8節点CHEXA、6節点CPENTAの各要素。
積分スキーム NA PSOLIDPLSOLID、またはPSHELLISOPEフィールド。

PARAM,EXPISOP(パラメータはソリッド要素でのみサポートされています)。

ISOPEフィールドは、PARAM,EXPISOPで定義された設定を上書きします。

積分スキームの詳細については、ユーザーズガイド内の要素をご参照ください。

拘束:
剛体のサポート NA RBE2RBE3およびRBODYがサポートされています。
材料:
サポートされている材料 NA シェル:MAT1MAT2MAT8、および MATS1

ソリッド:MAT1MATS1MATVE、およびMATHE

MATS1の場合:陰解法でサポートされているMATS1上の材料に加えて、Johnson-Cookとクラッシャブルフォームの材料もサポートされています。

MATHEの場合:MATHEにリストされているすべての材料モデルは、ソリッド要素でサポートされます。

プロパティ:
サポートされているプロパティ NA PSHELL, PSOLID, PLSOLID, PCOMP, PCOMPG, PCOMPP, PCOMPLS PLYSTACKおよびDRAPEエントリをサポート。
接触:
サポートされている接触タイプ NA CONTACTおよびTIE N2SS2Sの接触離散化がサポートされています。

SMALLFINITE、およびCONSLIの接触がサポートされています。

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

陽解法でのTIEの場合:
1
運動学的TIEのみがサポートされています。すなわち、ペナルティに基づいた手法を使用する代わりに、運動条件が正確に拘束されます。
2
運動TIE内の階層はサポートされていません(すなわち、TIEのセカンダリ(旧称“スレーブ”)節点は別のTIE内のメイン(旧称“マスター”)節点になることはできません)。
3
過剰に拘束されたTIEは無視されます(.femファイル内の入力の順序に基づいて、そのようなケースの最初の拘束のみが保持されます)。
4
このような階層と過剰に拘束されたTIE節点はすべて、*_badtied.femファイル内の節点SET内に出力されます。
座標系:
サポートされているユーザー定義の座標系 NA CORD2RCORD1CCORD2CCORD1S、およびCORD2S
出力:
ASCIIの出力 NA PARAM,NOUTCYC 陽解法の時間サイクルサマリーおよび対応する情報(時間ステップ、エネルギー、最大貫通、質量変化率など)のみが.outファイルに出力されます。PARAM,NOUTCYCを使用して、.outファイルへのサマリー出力の頻度を選択できます。
バイナリファイル出力 DISPVELOCITYACCELERATIONSTRESSSTRAIN(塑性ひずみを含む)速度依存性塑性問題のひずみ速度、CBUSHFORCECSTRESSCSTRAINCFAILUREESE NA 結果は_expl.h3dファイルと_expl.mvwファイルのみに出力されます。
_expl.h3d
変位、回転、速度、加速度、応力、ひずみ、ひずみ速度(速度依存の塑性の場合)、CBUSH力、塑性ひずみ複合材応力、複合材ひずみ、複合材破壊指数の結果が出力されます。
_expl.mvw
対応する_expl.h3dファイルを自動的に読み込んで、_expl.h3dファイル内の結果出力をプロットできるようにします。
_s<ID>_e.expl
内部エネルギー、弾性接触エネルギー、塑性接触エネルギー、運動エネルギー、アワグラスエネルギー、および塑性散逸エネルギーの出力の曲線が格納されています。
_expl_energy.mvw
対応する_s<ID>_e.explファイルを自動的に読み込んで、さまざまなエネルギー出力をプロットできるようにします。

ESE出力は、.h3dフォーマットのみで、COMPグループとOCOMPグループのオプションに使用できます。

THISTは、_TH.h5ファイル内の特定の結果の時刻歴出力を生成するために使用できます。

MONVOLバルクデータエントリでモニター容積を定義すると、圧力、温度、体積、面積、質量、内部エネルギー、質量流量、ベント面積、漏れ質量の出力結果がデフォルトで利用できます。

出力制御 NLOUT, THIST NLOUT, THIST 陽解法解析ではNINTフィールドのみがサポートされています。

非線形陽解法解析では、NLADAPTエントリはサポートされていません。

その他:
大変位解析 NA NA 陽解法非線形解析は、デフォルトでは大変位非線形解析です。
アダプティブ動的緩和 DYREL NA
モニターボリューム NA MONVOL 膨張気体のハイブリッド入力により、単室の気体充填構造を定義します。
Material Failure Criterion NA MATF またはMATS1(損傷開始 / 進展基準の場合のみ)。
Mass scaling control MSCALE NA
アワグラス制御 HOURGLSHGIDフィールドは、PSOLID/PLSOLID/PSHELL上でこのカードを参照します)

PARAM,HOURGLS

デフォルトのアワグラス値は、PSOLID/PLSOLID/PSHELLエントリで参照されるHOURGLSエントリ、またはPARAM,HOURGLSによって上書きされます。

HOURGLSエントリを介したHGIDは、PARAM,HOURGLSを上書きします。

詳細については、アワグラス制御をご参照ください。

最適化:
最適化のサポート サポートされていません。 サポートされていません。 サポートされていません。
1 デフォルト値は、ユーザー定義のPARAM,HOURGLSまたはHOURGLSエントリ(プロパティエントリのHGIDで参照されます)によって上書きできます
2 ユーザーは、PARAM,HOURGLSまたはHOURGLSエントリ(プロパティエントリのHGIDで参照されます)を使用してアワグラス制御をオンにできます
3 ソリッド要素については、PSOLID/PLSOLIDエントリのISOPEフィールドを使用して、積分スキームを切り替えることができます
4 CTETRA(1次と2次)にはアワグラス制御は適用されません
5 シェル要素については、PSHELLエントリのISOPEフィールドを使用して、積分スキームを切り替えることができます
6 デフォルト値は、ユーザー定義のPARAM,HOURGLSまたはHOURGLSエントリ(プロパティエントリのHGIDで参照されます)によって上書きできます。MAT1/MATS1/MATHEについては、1次CHEXA要素の場合のみにデフォルト値が適用されます。CPENTA要素については、必要に応じてアワグラス制御をオンにする必要があります。
7 ここで列挙している一部の材料は、シェルに対してはサポートされていません(MATHEMATVEなど)。