Loading [MathJax]/jax/output/CommonHTML/jax.js

/MAT/LAW93 (ORTH_HILL)または(CONVERSE)

ブロックフォーマットキーワード この材料則はHill塑性を擁する直交異方性弾性挙動材料を記述し、シェルおよびソリッド要素に適用可能です。(/BRICK/TETRA4および/TETRA10)。

シェル要素の場合はプロパティセット/PROP/TYPE11/PROP/TYPE17/PROP/TYPE51/PROP/PCOMPP、ソリッド要素の場合は/PROP/TYPE6と使用されます。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MAT/LAW93/mat_ID/unit_IDまたは/MAT/ORTH_HILL/mat_ID/unit_IDまたは/MAT/CONVERSE/mat_ID/unit_ID/
mat_title
ρi                
E11 E22 E33 G12 ν12
G13 G23 ν13 ν23  
Nrate VP Fcut            
Nrate>0 の場合の降伏の曲線入力、行毎にNrate 塑性関数を定義:
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
fct_IDi   Fscalei ˙εi        
降伏のパラメータ入力:
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
σy QR1 CR1 QR2 CR2
HILL基準の降伏応力比:
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
R11 R22 R12    
R33 R13 R23    

定義

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

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
mat_title 材料のタイトル

(文字、最大100文字)

 
ρi 初期密度

(実数)

[kgm3]
E11 方向11のヤング率

(実数)

[Pa]
E22 方向22のヤング率

(実数)

[Pa]
E33 方向33のヤング率

(実数)

[Pa]
G12 方向12におけるせん断係数

(実数)

[Pa]
G13 方向13におけるせん断係数

(実数)

[Pa]
G23 方向23におけるせん断係数

(実数)

[Pa]
ν12 ポアソン比12

(実数)

 
ν13 ポアソン比13

(実数)

 
ν23 ポアソン比23

(実数)

 
Nrate 降伏関数の数  
VP ひずみ速度選択フラグ
= 1
降伏応力に対するひずみ速度効果は塑性ひずみ速度に依存します。
= 2(デフォルト)
降伏に対するひずみ速度効果は全ひずみ速度に依存します。
= 3
降伏に対するひずみ速度効果はひずみ速度の偏差に依存します。

(整数)

 
Fcut ひずみ速度フィルタリングのカットオフ周波数。

デフォルト = 1.0 x 104(実数)

[Hz]
fct_IDi 塑性曲線のi番目の関数の識別子(i=1、Nrate

(整数)

 
Fscalei i番目の関数のスケールファクター(i=1、Nrate

デフォルト = 1.0(実数)

[Pa]
˙εi i番目の関数のひずみ速度(i=1、Nrate

(実数)

[1s]
σy 初期降伏応力

デフォルト = 1E30 (実数)

[Pa]
QR1 硬化のパラメータ

デフォルト = 0.0(実数)

[Pa]
CR1 硬化のパラメータ

デフォルト = 0.0(実数)

 
QR2 硬化のパラメータ

デフォルト = 0.0(実数)

[Pa]
CR2 硬化のパラメータ

デフォルト = 0.0(実数)

 
R11 方向11の降伏応力比

デフォルト = 1.0(実数)

 
R22 方向22の降伏応力比

デフォルト = 1.0(実数)

 
R33 方向33の降伏応力比

デフォルト = 1.0(実数)

 
R12 方向12の降伏応力比

デフォルト = 1.0(実数)

 
R13 方向13の降伏応力比

デフォルト = 1.0(実数)

 
R23 方向23の降伏応力比

デフォルト = 1.0(実数)

 

例(曲線入力)

コメント

  1. 降伏応力は、直交異方性フレームの相当応力と比較されます。ソリッド要素の場合、この相当応力は次のように定義されます:(1)
    σeq=F(σ222σ233)+G(σ233σ211)+H(σ211σ222)+2Lσ223+2Mσ231+2Nσ212
    ここで、
    F=12(1R222+1R2331R211)
    G=12(1R233+1R2111R222)
    H=12(1R222+1R2111R233)
    L=32R223
    M=32R231
    N=32R212
    Rii=σiiFσF
    法線方向
    Rij=3σijFσF
    せん断方向
    σijF
    方向の降伏応力 ij
    σF
    Voce硬化の合計で定義したり、表形式にすることができるグローバル流動応力(以下を参照)。
    平面応力条件下では、シェル要素の場合、相当降伏応力は次のようになります:(2)
    σeq=(G+H)σ211+(F+H)σ2222Hσ11σ22+2Nσ212
  2. 降伏関数 Φ は、Hillの相当応力 σeq を流動応力 σF と次のように比較します:(3)
    Φ=σeqσF
    流動応力 σF を定義する方法には異なる2つの方法があります:パラメータ入力または曲線入力です。
    • パラメータ入力の場合、流応力は、初期降伏応力と二重Voce硬化を用いて次のように定義されます。(4)
      σF=σ0Y+R(εp)

      ここで、 R(εp)=2iQRi(1eCRiεp)

    • 曲線入力の場合、パラメータ入力の値は無視されます。
      降伏は、ひずみ速度効果を考慮し、応力対塑性ひずみ曲線を用いて定義できます。応力対ひずみ曲線が定義されている場合、これが硬化を定義するデフォルトの方法となります。
      1. ˙ε˙εn の場合、降伏は fn fn1 の間で補間されます。
      2. ˙ε˙ε1 の場合、関数 f1 が使用されます。
      3. ˙εmax を超えた場合、降伏は外挿されます。


        図 1.
  3. 表形式の流動応力の場合、ひずみ速度 ˙ε の計算はフラグVPの値に依存します。
    • VP= 1の場合、塑性ひずみ速度が使用されます
    • VP= 2の場合、全ひずみ速度が使用されます
    • VP= 3の場合、全ひずみ速度が使用されます

    すべての場合において、ひずみ速度の計算にはフィルタリングが含まれます。VP = 1の場合、カットオフ周波数が自動的に設定されます。ただし、VP = 1または3の場合、カットオフ周波数Fcutを入力できます。それ以外の場合は、デフォルト値が設定されます。