/MAT/LAW66

ブロックフォーマットキーワードこの材料則は、応力ひずみ（応力に対する塑性ひずみ）による加工硬化部について、ユーザー定義関数を使用して等方性引張 / 圧縮弾塑性材料則をモデル化します。この材料則は、圧縮と引張に対して定義できます。

フォーマット

(1)	(3)	(5)	(7)	(9)	(10)
/MAT/LAW66/`mat_ID`/`unit_ID`
`mat_title`
$ρ_{i}$
`E`	$υ$	`C`_hard	`F`_cut	`F`_smooth	`I`_{yld_rate}
`P`_c	`P`_t

I_{yld_rate} = 0、1、または2の場合にのみ読み出し

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
`fct_ID`_c	`fct_ID`_t	`Fscale`_c		`Fscale`_t
${\dot{ε}}_{0}$		`c`		$σ_{y 0}$		`VP`

I_{yld_rate} = 3の場合のみ読み込み

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
`fct_ID`_c	`fct_ID`_t	`Fscale`_c		`Fscale`_t
`Frate_ID`_c	`Frate_ID`_t	`Fscale_rate`_c		`Fscale_rate`_t

I_{yld_rate} = 4の場合のみ読み込み

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
`NFUNCC`	`NFUNCT`

右記のそれぞれに対して； NFUNCC

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
`fct_ID`_c		${\dot{ε}}_{i}^{c}$		`Fscale`_c

右記のそれぞれに対して； NFUNCT

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
`fct_ID`_t		${\dot{ε}}_{i}^{t}$		`Fscale`_t

定義

フィールド	内容	SI 単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	単位識別子。（整数、最大10桁）
`mat_title`	材料のタイトル（文字、最大100文字）
$ρ_{i}$	初期密度（実数）	$[\frac{kg}{m^{3}}]$
`E`	ヤング率（実数）	$[Pa]$
$υ$	ポアソン比。（実数）
`C`_hard	硬化係数。 = 0 硬化は完全等方性モデルです。 = 1 硬化は運動学的Prager-Zieglerモデルです。 = 0～1の値硬化は2つのモデル間で補間されます。（実数）
`F`_smooth	ひずみ速度スムージングオプションフラグ = 0（デフォルト）ひずみ速度を平滑化しません。 = 1 ひずみ速度スムージングはアクティブ。（整数）
`F`_cut	ひずみ速度フィルタリングのカットオフ周波数、 Appendix：フィルタリングデフォルト = 10³⁰（実数）	$[Hz]$
`I`_{yld_rate}	降伏応力に対する速度効果フラグ = 1（デフォルト） Cowper-Symondsを使用： $1 + {(\frac{\dot{ε}}{{\dot{ε}}_{0}})}^{\frac{1}{c}}$ = 2 次を使用： $1 + cLn (\frac{\dot{ε}}{{\dot{ε}}_{0}})$ = 3 2つの荷重曲線を使用して圧縮降伏応力（`fct_ID`_c）と引張降伏応力（`fct_ID`_t）をスケーリングします。 = 4 ひずみ速度の値が異なれば、異なる圧縮関数と引張関数を使用します。（整数）
`P`_c	圧縮圧力の制限デフォルト = 0（実数）	$[Pa]$
`P`_t	引張圧力の制限デフォルト = 0（実数）	$[Pa]$
`fct_ID`_c	圧縮降伏応力（整数）
`fct_ID`_t	引張降伏応力（整数）
`Fscale`_c	`fct_ID`_cの縦軸（応力）のスケールファクターデフォルト = 1.0（実数）	$[Pa]$
`Fscale`_t	`fct_ID`_tの縦軸（応力）のスケールファクターデフォルト = 1.0（実数）	$[Pa]$
`c`	ひずみ速度パラメータ（実数）
${\dot{ε}}_{0}$	参照ひずみ速度デフォルト = 1.0（実数）	$[\frac{1}{s}]$
$σ_{y 0}$	初期降伏応力デフォルト = 0（実数）	$[Pa]$
`VP`	ひずみ速度選択フラグ = 0 降伏応力に対するひずみ速度効果は全ひずみ速度に依存します。 = 1 降伏応力に対するひずみ速度効果は塑性ひずみ速度に依存します。この場合、ひずみ速度のフィルタリングはないので、`F`_smoothと`F`_cutは使用されません。`I`_{yld_rate} = 1（Cowper Symondsの式）の場合にのみ使用できます。 2 （整数）
`Frate_ID`_c	圧縮ひずみ速度効果関数の識別子（整数）
`Frate_ID`_t	引張ひずみ速度効果関数の識別子（整数）
`Fscale_rate`_c	`Frate_ID`_cの縦軸（応力）のスケールファクターデフォルト = 1.0（実数）	$[Pa]$
`Fscale_rate`_t	`Frate_ID`_tの縦軸（応力）のスケールファクターデフォルト = 1.0（実数）	$[Pa]$
`NFUNCC`	圧縮関数の数（整数）
`NFUNCT`	引張関数の数（整数）
${\dot{ε}}_{i}^{c}$	`i`番目の圧縮ひずみ速度`i` =1、`NFUNCC` （実数）	$[\frac{1}{s}]$
${\dot{ε}}_{i}^{t}$	`i`番目の引張ひずみ速度`i`=1、`NFUNCT` （実数）	$[\frac{1}{s}]$

例（アルミニウム）

#RADIOSS STARTER
/UNIT/1
unit for mat
                   g                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW66/1/1
Aluminium
#              RHO_I
               .0027
#                  E                  Nu              C_hard               F_cut  F_smooth Iyld_rate
               60400                 .33                   0                   0         0         4
#                P_c                 P_t
                 500                 600
#   NFUNCC    NFUNCT
         2         2
#funct_IDc                    Episilon_c             Fscalec
        38                            10                   1                                                  
        40                            40                 1.6                                                  
#funct_IDt                    Episilon_t             Fscalet
        38                            10                   1                                                  
        40                            40                 1.6                                                  
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  3. FUNCTIONS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/38
function_38
#                  X                   Y
                   0                  90                                                            
                 .08                 170                                                            
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/40
function_40
#                  X                   Y
                   0                  90                                                            
                 .08                 170                                                            
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

これは等方性弾塑性則です。降伏応力の定義には、圧縮降伏応力および引張降伏応力と、圧縮と引張の両方の有効塑性ひずみが使用されます。2つの圧力P_tまたはP_cを超えた場合、どちらを超えたかによって、引張降伏応力を使用するか、圧縮降伏応力を使用するかが決定されます。
圧力がこれら2つの値の間にある場合、降伏応力は次の式によって定義されます。

次の場合； $- P_{t} \leq P \leq P_{c}$ (1)
$\begin{array}{l} σ_{y} = α σ_{y}^{t} (ε_{p}) + (1 - α) σ_{y}^{c} (ε_{p}) \\ α = \frac{P_{c} - P}{P_{e} + P_{t}} \end{array}$

$P_{t} = P_{c} = 0$ の場合、または圧力が2つの値の範囲を超える場合、降伏応力は次の式によって定義されます：

$σ_{y} = σ_{y}^{t} (ε_{p})$ 次の場合； $P \leq 0$

$σ_{y} = σ_{y}^{c} (ε_{p})$ 次の場合； $P > 0$
降伏応力は以下のように計算されます：
VP= 1の場合：

$σ_{y} (ε_{p}, {\dot{ε}}_{p}) = σ_{y}^{s} (ε_{p}) + σ_{y 0} {(\frac{{\dot{ε}}_{p}}{{\dot{ε}}_{0}})}^{\frac{1}{c}}$ 次の場合； $σ_{y 0} > 0$

$σ_{y} (ε_{p}, {\dot{ε}}_{p}) = σ_{y}^{s} (ε_{p}) [1 + {(\frac{{\dot{ε}}_{p}}{{\dot{ε}}_{0}})}^{\frac{1}{c}}]$ 次の場合； $σ_{y 0} = 0$

VP= 0の場合：

$σ_{y} (ε_{p}, {\dot{ε}}_{p}) = σ_{y}^{s} (ε_{p})$ 次の場合； $σ_{y 0} > 0$

$σ_{y} (ε_{p}, {\dot{ε}}_{p}) = σ_{y}^{s} (ε_{p})$ 次の場合； $σ_{y 0} = 0$

ここで $σ_{y}^{t} (ε_{p})$ は静的降伏応力、 $σ_{y 0}$ は初期降伏応力です。
粘性効果を含めるには、この材料則と共に/VISC/PRONYを使用する必要があります。

/MAT/LAW66

フォーマット

定義

例（アルミニウム）

コメント