/MAT/LAW73

ブロックフォーマットキーワード この材料則はThermal Hill直交異方性材料を記述します。これはシェル要素にのみ適用可能です。

この材料則は、降伏応力がひずみ速度と塑性ひずみだけでなく温度(ユーザーテーブルで定義)にも依存する点で、/MAT/LAW43(HILL_TAB)とは異なります。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MAT/LAW73/mat_ID/unit_ID
mat_title
ρi                
E v            
fct_IDE   Einf CE 空白
r00 r45 r90 Chard Iyield0  
εpmax εt εm        
Tab_ID   σscale ε˙scale        
Ti Cp            

定義

フィールド 内容 SI単位の例
mat_ID 材料識別子

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
mat_title 材料のタイトル

(文字、最大100文字)

 
ρi 初期密度

(実数)

[kgm3]
E 初期ヤング率

(実数)

[Pa]
v ポアソン比

(実数)

 
fct_IDE ヤング率のスケールファクターの関数ID(ヤング率が塑性ひずみの関数である場合)
= 0(デフォルト)
この場合はヤングの進展はEinfCEに応じて決まります。

(整数)

 
Einf 無限の塑性ひずみに対する飽和ヤング率

(実数)

[Pa]
CE ヤング率進展のパラメータ

(実数)

 
r00 0°のランクフォードパラメータ 3

デフォルト = 1.0(実数)

 
r45 45°のランクフォードパラメータ

デフォルト = 1.0(実数)

 
r90 90°のランクフォードパラメータ

デフォルト = 1.0(実数)

 
Chard 硬化係数。
= 0
硬化は完全等方性モデルです。
= 1
硬化は運動学的Prager-Zieglerモデルです。
= 01の値
硬化は2つのモデル間で補間されます。

(実数)

 
Iyield0 降伏応力フラグ
= 0
平均降伏応力入力
= 1
直交方向1の降伏応力

(整数)

 
εpmax 破壊塑性ひずみ。

デフォルト = 1.0x1030(実数)

 
εt 応力が減少し始める引張破壊ひずみ。

デフォルト = 1.0x1030(実数)

 
εm 要素内の応力が0に設定される最大引張破壊ひずみ

デフォルト = 2.0x1030(実数)

 
Tab_ID 降伏応力を定義するためのテーブル識別子

(整数)

 
σscale 降伏応力のスケールファクター

デフォルトは1.0に設定されます(実数)

[Pa]
ε˙scale ひずみ速度のスケールファクター

デフォルトは1.0に設定されます(実数)

[1s]
Ti 初期温度。

デフォルトは293Kに設定されます(実数)

[K]
Cp 単位質量あたりの比熱

(実数)

[JkgK] MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeaaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaadaWcaaqaaiaabQeaaeaacaqGRbGaae4zaiabgwSixlaabUeaaaaacaGLBbGaayzxaaaaaa@3DB3@

例(鋼材)

コメント

  1. この材料則は、プロパティセット/PROP/TYPE9 (SH_ORTH)または/PROP/TYPE10 (SH_COMP)と共に使用する必要があります。
  2. 降伏応力はユーザー関数で定義し、次の等価応力と比較されます:
    (1)
    σeq=A1σ12+A2σ22-A3σ1σ2+A12σ122
  3. ランクフォードパラメータの角度は、直交異方性方向1に対して定義されます。
    R=r00+2r45+r904 H=R1+R
    A1=H(1+1r00) A2=H(1+1r90)
    A3=2H A12=2H(r45+0.5)(1r00+1r90)
    r00=A32A1A3 r45=12(A12A1+A2A31)
    r90=A32A2A3  
    ランクフォードパラメータ rα は、面内の塑性ひずみと厚み方向の塑性ひずみ ε33 との比率です。(2)
    rα=dεα+π/2dε33

    ここで、α MathType@MTEF@5@5@+=feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@3792@ は、直交異方性方向1に対して成す角度です。

    このランクフォードパラメータ rα は、角度α MathType@MTEF@5@5@+=feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@3792@ で実施する簡易な引張試験によって特定できます。

    Rの大きな値は、成形性に優れていることを意味します。

  4. εp (塑性ひずみ)が1つの積分点で εpmax に到達すると、対応するシェル要素が削除されます。
  5. 最大主ひずみが ε1>εt の場合、応力は次の関係式に従って減少します: (3)
    σ=σ(εmε1εmεt)
  6. ε1>εm の場合、応力は0に減少します(ただし、要素は削除されません)。
  7. この材料則は、常に反復投影を使用して可塑性を計算します(プロパティセットのIplasは無視されます)。
  8. この材料則は、可塑性のための全体定式化では使用できません(シェルプロパティでN=0は使用できません)。
  9. 降伏応力定義用のテーブルは、塑性ひずみ、ひずみ速度、温度をそれぞれ表すパラメータ (εp,ε˙,T) を含む3次元テーブルである必要があります。テーブルの値は、降伏応力値です。
  10. εm1pεpεmp かつ ε˙n1ε˙ε˙n かつ Tq1TTq である場合、降伏は、 (εip,ε˙j,Tk),i=m1,m;j=n1,n;k=q1,q に対応するテーブルの8つの値の間で線形補間されます。
  11. (εp,ε˙,T) がテーブルの範囲に収まらない場合、降伏応力は線形外挿によって得られます。したがって、ひずみ速度ゼロに対応する静的曲線をテーブルに入力する必要があります( ε˙=0 のエントリをテーブル定義に含める必要があります)。
    /HEAT/MATオプションが材料識別子に関連付けられていない場合、断熱条件が仮定されて温度は次のように計算されます:(4)
    Τ=Ti+EintρCp(Volume)
    ここで、
    Eint
    Radiossによって計算される内部エネルギー
    σ
    Volumeは現在の密度と体積
    Cp
    単位質量当たりの熱容量

    上記とは別の方法を指定する場合、熱伝導のための有限要素定式化が必要になります(/HEAT/MATオプションでIform =1を設定します)。この場合、/HEAT/MATオプションに入力された初期温度と比熱が使用されます。

  12. 降伏応力が直交方向1で取得された場合は、Iyield0 =1を定義します。そうでない場合は、Iyield0=0を定義します。