弾塑性材料

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$

材料パラメータ材料の応力-ひずみ曲線のフィッティング（例： Altair Composeスクリプト）により、これら2つのパラメータが求められます。

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

図 2. 引張試験

ひずみ速度

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

σ = (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})$ により考慮されます。(3)

σ = (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})$ ここで、(4)

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}$: 室温（単位はケルビン）。

T

$T$ は、以下で計算されます。(5)

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応力：(6) $\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}$
シェル要素：(7) $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異方性パラメータです。シェル要素で必要なHILLパラメータは、 $F$ $F$ 、 $G$ $G$ 、 $H$ $H$ 、および $N$ $N$ の4つのみです。

LAW78ではHILL基準は次のとおりです：(8) $φ (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}$ との比率です。(9)

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

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

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

$r_{α}$ $r_{α}$ は、直交異方性方向1に対してさまざまな角度で切断した多くの試料で測定できます。荷重の方向を直交異方性の方向1とした引張試験で測定した $r_{00}$ $r_{00}$ と同様。 $r_{90}$ $r_{90}$ 荷重の方向が直交異方性の方向1と直交する引張試験で測定した

ひずみ速度は、試料の厚み方向のひずみに対する試料の幅方向のひずみの比率です。

この場合のHILLパラメータは次のようになります：(10)

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

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

(11)

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

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

(12)

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

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

(13)

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基準は次のとおりです：(14)

σ_{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で使用する降伏応力比は次のとおりです：(15)

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によって自動的に計算されます。(16)

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}})$

(17)

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}})$

(18)

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}})$

(19)

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$ を取得します。