弾塑性材料

Johnson-Cook（/MAT/LAW2）

LAW2には、応力計算のための3つのパートがあります。

塑性ひずみの影響
ひずみ速度の影響
温度変化の影響

材料パラメータ

LAW2で材料パラメータを入力する方法は2つあります。

Iflag=0： Johnson-Cookパラメータ $a$ $a$ 、 $b$ $b$ 、 $n$ $n$ ：アクティブ
Iflag=1：降伏応力、UTS（公称応力）、またはUTSでのひずみによる、新しい簡素化された入力

Iflag= 0: (1) $σ = a + b \cdot ε_{p}^{n}$ $σ = a + b \cdot ε_{p}^{n}$; ここで、

$a$ $a$

材料試験で読み取られ、真応力に変換される可能性のある降伏応力です。

$b$ $b$ および $n$ $n$

材料パラメータ材料の応力-ひずみ曲線のフィッティングにより、これら2つのパラメータが求められます。試験から生成された応力-ひずみ曲線がない場合は、式 2 および式 3を使用して、 $b$ $b$ と $n$ $n$ を計算するためには2つの状態（2つの応力-ひずみポイント）が必要となります。最初の点はネッキングポイントで選ばれ（まず、 R_mの取得を試みる）、次に $b$ $b$ と $n$ $n$ が、その結果は選ばれた点によって異なるために、曲線のそれぞれの他の点に対して計算され、平均化されます。(2) $n = \frac{\ln (\frac{σ_{1} - a}{σ_{2} - a})}{\ln (\frac{ε_{1}}{ε_{2}})}$ $n = \frac{\ln (\frac{σ_{1} - a}{σ_{2} - a})}{\ln (\frac{ε_{1}}{ε_{2}})}$ (3) $b = \frac{σ_{1} - a}{ε_{1}^{n}}$ $b = \frac{σ_{1} - a}{ε_{1}^{n}}$

このテストの目的は引張り試験結果を用いて材料則のパラメーターを引き出す方法を提案することです。

Iflag = 1

この新しい入力では、ネッキングポイントでの降伏応力（ $σ_{y}$ $σ_{y}$ ), 引張り強さ（UTS）および工学ひずみ（ $ε_{U T S}$ $ε_{U T S}$ ）が必要です。この新しい入力により、Radiossは自動的に $a$ $a$ 、 $b$ $b$ および $n$ $n$ の等価値を計算します。

図 2. 引張試験

ひずみ速度

ひずみ速度は、引張または破壊における衝突パフォーマンスで、材料特性に大きな影響を与えます。Johnson-Cook理論では、降伏応力は直接ひずみ速度の影響を受け、次のように表されます：(4)

σ = (a + b \cdot ε_{p}^{n}) (1 + c \ln \frac{\dot{ε}}{{\dot{ε}}_{0}})

$σ = (a + b \cdot ε_{p}^{n}) (1 + c \ln \frac{\dot{ε}}{{\dot{ε}}_{0}})$

一般に、試験ひずみ速度が増加すると降伏応力は増加します。ひずみ速度係数

c

$c$ により、降伏応力の増加係数をスケーリングできます。

c

$c$ =0の場合、または

{\dot{ε}}_{0} = 10^{30}

${\dot{ε}}_{0} = 10^{30}$ あるいは

\dot{ε} \leq {\dot{ε}}_{0}

$\dot{ε} \leq {\dot{ε}}_{0}$ の場合、ひずみ速度の影響もまた定義されません。

温度変化

温度が上昇すると降伏応力は低下します。LAW2では、影響は

(1 - T^{* m})

$(1 - T^{* m})$ により考慮されます。(5)

σ = (a + b \cdot ε_{p}^{n}) (1 + c \ln \frac{\dot{ε}}{{\dot{ε}}_{0}}) (1 - T^{* m})

$σ = (a + b \cdot ε_{p}^{n}) (1 + c \ln \frac{\dot{ε}}{{\dot{ε}}_{0}}) (1 - T^{* m})$ ここで、(6)

T^{*} = \frac{T - T_{r}}{T_{m e l t} - T_{r}}

$T^{*} = \frac{T - T_{r}}{T_{m e l t} - T_{r}}$

ここで、

$T_{m e l t}$ $T_{m e l t}$: 溶融温度（単位はケルビン）
$T_{r}$ $T_{r}$: 室温（単位はケルビン）

(7)

T = T_{i} + \frac{E_{int}}{ρ C_{p} (V o l u m e)}

$T = T_{i} + \frac{E_{int}}{ρ C_{p} (V o l u m e)}$

ここで、

$E_{int}$ $E_{int}$: 内部エネルギー; 内部エネルギーの変化は、Johnson-Cook則で降伏応力に影響を与えます。

硬化係数

金属は降伏するまで変形し、その後一般には硬化します（降伏応力は増加）。材料により硬化の様子は異なります（等方硬化、移動硬化など）。これは非常に重要な材料特性でもあります（スプリングバックの場合）。

LAW2では、オプションC_hard（硬化係数）を使用して、材料にどの硬化モデルを使用するかを記述します。この機能はLAW36、43、44、57、60、66、73、74でも使用できます。

C_hardの値は1～0です。等方モデルの場合はC_hard=0、移動Prager-Zieglerモデルの場合はC_hard=1、これら2つのモデルの間の硬化の場合は1と0の間となります。

C_hard= 0：等方性モデル: 1次元のケースでは、材料は降伏応力後に強化されます。前回の引張りの最大応力がそれに続く荷重での降伏となり、この新しい降伏応力はそれに続く引張りおよび圧縮での降伏応力と同じになります。

図 4.
C_hard= 1：運動学的Prager-Zieglerモデル: Bauschinger効果（引張りによる硬化の後、圧縮による軟化が発生し、圧縮での平均降伏が低下する）をモデル化するには、移動硬化を使用します。

図 5.

弾塑性区分線形材料（/MAT/LAW36）

LAW36では、さまざまなひずみ速度に対してさまざまな塑性応力-ひずみ曲線を直接定義できます。

大きなひずみ速度の塑性応力-ひずみ曲線は、必ず小さなひずみ速度の塑性応力-ひずみ曲線より上になります。

ヤング率

ヤング率は、オプションfct_ID_E、E_inf、およびC_Eを使用して、除荷時に更新（低減）できます。この機能の使用により、ハイテン鋼のスプリングバックの精度（除荷相時）が向上します。この機能は材料LAW43、LAW57、LAW60、LAW74およびLAW78でも使用できます。

fct_ID_Eを使用したヤング率の更新（fct_ID_E ≠ 0）：

図 7.
E_infおよびC_Eを使用したヤング率の更新（fct_ID_E = 0）：

図 8.

材料の挙動

fct_ID_pは、特定の材料における引張と圧縮の挙動の区別（圧力依存降伏）に使用されます。したがって、有効降伏応力は公称降伏応力に実際の圧力に対応する降伏係数を乗じることによって得られます。

HILL材料

Radioss材料則では、LAW32、LAW43、LAW72、LAW73、LAW74、LAW78およびLAW93はHILL基準を使用します。

HILL基準

代表的なHILL基準は：

3D等価HILL応力：(8) $\begin{array}{l} f = \sqrt{F {(σ_{y y} - σ_{z z})}^{2} + G {(σ_{z z} - σ_{x x})}^{2} + H {(σ_{x x} - σ_{y y})}^{2} + 2 L σ_{y z}^{2} + 2 M σ_{z x}^{2} + 2 N σ_{x y}^{2}} \\ = \sqrt{\underset{}{\underset{︸}{(G + H)}} σ_{x x}^{2} + \underset{}{\underset{︸}{(F + H)}} σ_{y y}^{2} + \underset{}{\underset{︸}{(F + G)}} σ_{z z}^{2} - \underset{}{\underset{︸}{2 H}} σ_{x x} σ_{y y} - \underset{}{\underset{︸}{2 F}} σ_{y y} σ_{z z} - \underset{}{\underset{︸}{2 G}} σ_{z z} σ_{x x} + \underset{}{\underset{︸}{ 2 L}} σ_{y z}^{2} + \underset{}{\underset{︸}{2 M}} σ_{z x}^{2} + \underset{}{\underset{︸}{2 N}} σ_{x y}^{2}} \end{array}$ $\begin{array}{l} f = \sqrt{F {(σ_{y y} - σ_{z z})}^{2} + G {(σ_{z z} - σ_{x x})}^{2} + H {(σ_{x x} - σ_{y y})}^{2} + 2 L σ_{y z}^{2} + 2 M σ_{z x}^{2} + 2 N σ_{x y}^{2}} \\ = \sqrt{\underset{}{\underset{︸}{(G + H)}} σ_{x x}^{2} + \underset{}{\underset{︸}{(F + H)}} σ_{y y}^{2} + \underset{}{\underset{︸}{(F + G)}} σ_{z z}^{2} - \underset{}{\underset{︸}{2 H}} σ_{x x} σ_{y y} - \underset{}{\underset{︸}{2 F}} σ_{y y} σ_{z z} - \underset{}{\underset{︸}{2 G}} σ_{z z} σ_{x x} + \underset{}{\underset{︸}{ 2 L}} σ_{y z}^{2} + \underset{}{\underset{︸}{2 M}} σ_{z x}^{2} + \underset{}{\underset{︸}{2 N}} σ_{x y}^{2}} \end{array}$
シェル要素：(9) $f = \sqrt{F σ_{y y}^{2} + G σ_{x x}^{2} + H {(σ_{x x} - σ_{y y})}^{2} + 2 N σ_{x y}^{2}} = \sqrt{\underset{}{\underset{︸}{(G + H)}} σ_{x x}^{2} + \underset{}{\underset{︸}{(F + H)}} σ_{y y}^{2} - \underset{}{\underset{︸}{2 H}} σ_{x x} σ_{y y} + \underset{}{\underset{︸}{2 N}} σ_{x y}^{2}}$ $f = \sqrt{F σ_{y y}^{2} + G σ_{x x}^{2} + H {(σ_{x x} - σ_{y y})}^{2} + 2 N σ_{x y}^{2}} = \sqrt{\underset{}{\underset{︸}{(G + H)}} σ_{x x}^{2} + \underset{}{\underset{︸}{(F + H)}} σ_{y y}^{2} - \underset{}{\underset{︸}{2 H}} σ_{x x} σ_{y y} + \underset{}{\underset{︸}{2 N}} σ_{x y}^{2}}$
ここで、 $F$ $F$ 、 $G$ $G$ 、 $H$ $H$ 、 $L$ $L$ 、 $M$ $M$ および $N$ $N$ は6つのHILL異方性パラメータ。シェル要素の場合、 $F$ $F$ 、 $G$ $G$ 、 $H$ $H$ および $N$ $N$ のみが、必要とされる4つのHILLパラメータです。

LAW78ではHILL基準は：(10) $φ (A) = \underset{G + H}{\underset{︸}{1}} \cdot A_{x x}^{2} - \underset{2 H}{\underset{︸}{\frac{2 r_{0}}{1 + r_{0}}}} A_{x x} A_{y y} + \underset{F + H}{\underset{︸}{\frac{r_{0} (1 + r_{90})}{r_{90} (1 + r_{0})}}} A_{y y}^{2} + \underset{2 N}{\underset{︸}{\frac{r_{0} + r_{90}}{r_{90} (1 + r_{0})} (2 r_{45} + 1)}} A_{x y}^{2}$ $φ (A) = \underset{G + H}{\underset{︸}{1}} \cdot A_{x x}^{2} - \underset{2 H}{\underset{︸}{\frac{2 r_{0}}{1 + r_{0}}}} A_{x x} A_{y y} + \underset{F + H}{\underset{︸}{\frac{r_{0} (1 + r_{90})}{r_{90} (1 + r_{0})}}} A_{y y}^{2} + \underset{2 N}{\underset{︸}{\frac{r_{0} + r_{90}}{r_{90} (1 + r_{0})} (2 r_{45} + 1)}} A_{x y}^{2}$
ランクフォードのパラメータを使ってHILLパラメータを決定する方法が2つあります。
- ひずみ比 $r_{00}, r_{45}, r_{90}$ $r_{00}, r_{45}, r_{90}$ （LAW32、LAW43、LAW72、LAW73）
- 降伏応力比 $R_{11}, R_{22}, R_{33}, R_{12}, R_{13}, R_{23}$ $R_{11}, R_{22}, R_{33}, R_{12}, R_{13}, R_{23}$ （LAW74、LAW93）

ひずみ比

ランクフォードパラメータ

r_{α}

$r_{α}$ は、面内の塑性ひずみと厚み方向の塑性ひずみ

ε_{33}

$ε_{33}$ との比率です。(11)

r_{α} = \frac{d ε_{α + π / 2}}{d ε_{33}}

$r_{α} = \frac{d ε_{α + π / 2}}{d ε_{33}}$

ここで、 $α$ $α$ は、直交異方性方向1に対して成す角度です。

$r_{α}$ $r_{α}$ は、直交方向1の異なる角度で切断した異なる試料で測定することができます。 $r_{00}$ $r_{00}$ は荷重方向が直交方向１に沿った引張試験から測定され、 $r_{90}$ $r_{90}$ は荷重が直交方向1に直交する引張試験から測定されます。

ひずみ比は、試料の幅方向のひずみと試料の厚さ方向のひずみとの比です。

この場合、HILLパラメータは：(12)

F = \frac{r_{00}}{r_{90} (r_{00} + 1)}

$F = \frac{r_{00}}{r_{90} (r_{00} + 1)}$

(13)

G = \frac{1}{(r_{00} + 1)}

$G = \frac{1}{(r_{00} + 1)}$

(14)

H = \frac{r_{00}}{(r_{00} + 1)}

$H = \frac{r_{00}}{(r_{00} + 1)}$

(15)

N = \frac{(1 + 2 r_{45}) (r_{00} + r_{90})}{2 r_{90} (r_{00} + 1)}

$N = \frac{(1 + 2 r_{45}) (r_{00} + r_{90})}{2 r_{90} (r_{00} + 1)}$

ここで、 $G + H = 1$ $G + H = 1$

LAW32、LAW43およびLAW73では、HILL基準は：(16)

σ_{e q} = \sqrt{A_{1} σ_{1}^{2} + A_{2} σ_{2}^{2} - A_{3} σ_{1} σ_{2} + A_{12} σ_{12}^{2}}

$σ_{e q} = \sqrt{A_{1} σ_{1}^{2} + A_{2} σ_{2}^{2} - A_{3} σ_{1} σ_{2} + A_{12} σ_{12}^{2}}$

$R = \frac{r_{00} + 2 r_{45} + r_{90}}{4}$ $R = \frac{r_{00} + 2 r_{45} + r_{90}}{4}$	$H = \frac{R}{1 + R}$ $H = \frac{R}{1 + R}$
$A_{1} = H (1 + \frac{1}{r_{00}})$ $A_{1} = H (1 + \frac{1}{r_{00}})$	$A_{2} = H (1 + \frac{1}{r_{90}})$ $A_{2} = H (1 + \frac{1}{r_{90}})$
$A_{3} = 2 H$ $A_{3} = 2 H$	$A_{12} = 2 H (r_{45} + 0.5) (\frac{1}{r_{00}} + \frac{1}{r_{90}})$ $A_{12} = 2 H (r_{45} + 0.5) (\frac{1}{r_{00}} + \frac{1}{r_{90}})$

これらはすべてランクフォードパラメータ（ひずみ比） $r_{00}, r_{45}, r_{90}$ $r_{00}, r_{45}, r_{90}$ を要求し、HILLパラメータ $A_{i}$ $A_{i}$ はRadiossによって自動的に計算されます。

降伏応力比

LAW93では、使用される降伏応力比は：(17)

R_{i j} = \frac{σ_{i j}}{σ_{0}}

$R_{i j} = \frac{σ_{i j}}{σ_{0}}$

降伏応力比

R_{i j}

$R_{i j}$ を得るには2つの荷重ケースでの降伏応力を測定する必要があります。

引張試験からの降伏応力 $σ_{11}, σ_{22}, σ_{33}$ $σ_{11}, σ_{22}, σ_{33}$
せん断試験からの降伏せん断応力 $σ_{12}, σ_{13}, σ_{23}$ $σ_{12}, σ_{13}, σ_{23}$

LAW93では、パラメータ入力が使用されている場合は初期応力パラメータ $σ_{y}$ $σ_{y}$ を基準降伏応力 $σ_{0}$ $σ_{0}$ とします。曲線入力を使用する場合は、曲線からの降伏応力を基準降伏応力 $σ_{0}$ $σ_{0}$ とします。

シェル用の4つのHILLパラメータがRadiossによって自動的に計算されます。(18)

F = \frac{1}{2} (\frac{1}{R_{22}^{2}} + \frac{1}{R_{33}^{2}} - \frac{1}{R_{11}^{2}})

$F = \frac{1}{2} (\frac{1}{R_{22}^{2}} + \frac{1}{R_{33}^{2}} - \frac{1}{R_{11}^{2}})$

(19)

G = \frac{1}{2} (\frac{1}{R_{33}^{2}} + \frac{1}{R_{11}^{2}} - \frac{1}{R_{22}^{2}})

$G = \frac{1}{2} (\frac{1}{R_{33}^{2}} + \frac{1}{R_{11}^{2}} - \frac{1}{R_{22}^{2}})$

(20)

H = \frac{1}{2} (\frac{1}{R_{22}^{2}} + \frac{1}{R_{11}^{2}} - \frac{1}{R_{33}^{2}})

$H = \frac{1}{2} (\frac{1}{R_{22}^{2}} + \frac{1}{R_{11}^{2}} - \frac{1}{R_{33}^{2}})$

(21)

N = \frac{3}{2 R_{12}^{2}}

$N = \frac{3}{2 R_{12}^{2}}$

LAW74では、降伏応力比

R_{i j}

$R_{i j}$ は降伏応力

σ_{11}, σ_{22}, σ_{33}

$σ_{11}, σ_{22}, σ_{33}$ および

σ_{12}, σ_{13}, σ_{23}

$σ_{12}, σ_{13}, σ_{23}$ 入力と直接使用され、ソリッド用の6つのHILLパラメータがRadiossによって自動的に計算されます。

$F = \frac{1}{2} (\frac{1}{σ_{22}^{2}} + \frac{1}{σ_{33}^{2}} - \frac{1}{σ_{11}^{2}})$ $F = \frac{1}{2} (\frac{1}{σ_{22}^{2}} + \frac{1}{σ_{33}^{2}} - \frac{1}{σ_{11}^{2}})$	$G = \frac{1}{2} (\frac{1}{σ_{22}^{2}} + \frac{1}{σ_{33}^{2}} - \frac{1}{σ_{11}^{2}})$ $G = \frac{1}{2} (\frac{1}{σ_{22}^{2}} + \frac{1}{σ_{33}^{2}} - \frac{1}{σ_{11}^{2}})$
$H = \frac{1}{2} (\frac{1}{σ_{22}^{2}} + \frac{1}{σ_{33}^{2}} - \frac{1}{σ_{11}^{2}})$ $H = \frac{1}{2} (\frac{1}{σ_{22}^{2}} + \frac{1}{σ_{33}^{2}} - \frac{1}{σ_{11}^{2}})$	$L = \frac{1}{2 σ_{23}^{2}}$ $L = \frac{1}{2 σ_{23}^{2}}$
$M = \frac{1}{2 σ_{31}^{2}}$ $M = \frac{1}{2 σ_{31}^{2}}$	$N = \frac{1}{2 σ_{12}^{2}}$ $N = \frac{1}{2 σ_{12}^{2}}$

シェル要素の場合、 $M = N$ $M = N$ 、 $L = N$ $L = N$ とします。