/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	`Ec`	`RPCT`

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]$
`Ec`	（オプション）圧縮のヤング率 2 （実数）	$[Pa]$
`RPCT`	`P`_cと`P`_tで使用されるスケールファクター。 2 （実数）
`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の式）の場合にのみ使用できます。 3 （整数）
`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
Aluminum
#              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                     Epsilon_c             Fscalec
        38                            10                   1                                                  
        40                            40                 1.6                                                  
#funct_IDt                     Epsilon_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----|

▸例（圧縮のヤング率 - オプション）

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
## Material Law 
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW66/1
20MAT124_20degree                                                                                   
#              RHO_I
1.25000000000000E-09
#                  E                  Nu              C_hard               F_cut  F_smooth Iyld_rate
            210000.0               0.33                 0.0                 0.0         0         1
#                P_c                 P_t                  EC                RPCT
                 0.0                 0.0              70000.
#funct_IDc funct_IDt             Fscalec             Fscalet
        34        34              100.0              200.0
#          Epsilon_0                   c            Sigma_Y0        VP
                 0.0                 0.0                 0.0         0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
##HWCOLOR curves 24 9
/FUNCT/24
20deg_TENSION                                                                                       
#                  X                   Y
                 0.0         0.200535124
1.11148000000000E-04          0.23893938
2.47763000000000E-04         0.274602617
3.84248000000000E-04         0.308164207
5.81974000000000E-04         0.354639794
7.98180000000000E-04         0.395218545
1.03131900000000E-03          0.43242102
         0.001382833         0.473801266
         0.001747862         0.508073173
         0.002139804         0.539268776
         0.002704889         0.577068218
         0.003299031         0.610692299
         0.004078009         0.646563755
         0.005715537         0.702190108
                 1.0                0.75
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/34
20deg_COMPRESSION                                                                                   
#                  X                   Y
                 0.0         0.709520996
         0.002556758         0.768542475
         0.005112572         0.814793042
         0.007700089         0.850991269
         0.010231756         0.877769462
         0.012846771         0.896395146
         0.015380438         0.909432849
         0.020500175         0.924867455
         0.030979299         0.948100519
         0.073545677         1.001804424
         0.126630132         1.075808748
         0.245056814          1.17707424
                 1.0                 1.2
#enddata

これは等方性弾塑性則です。降伏応力の定義には、圧縮降伏応力および引張降伏応力と、圧縮と引張の両方の有効塑性ひずみが使用されます。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$
Ecが定義されている場合、ヤング率は次のように計算されます：
- P > -RPCT * Pの場合、ヤング率はE。_t
- P < -RPCT * Pの場合、ヤング率はEc。_c
- P_t < P < RPCT * PRPCTの場合、線形補間は、EとEcの間で行われます。_c
降伏応力は以下のように計算されます：
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

フォーマット

定義

▸例（アルミニウム）

▸例（圧縮のヤング率 - オプション）

コメント