/MAT/LAW81

ブロックフォーマットキーワード この材料則はcapのあるDrücker-Prager降伏基準に基づいています。ここにはFosterの原理に基づくひずみ硬化capモデルがあります。塑性には等方硬化が含まれます。

破壊サーフェスは、圧力軸を中心に対称的な、標準の線形Drücker-Prager関係に制限されます。この材料則はLAG、ALE、EULERと適合性があります。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MAT/LAW81/mat_ID/unit_ID
mat_title
ρ i                
K0 G0 c0 Pb0  
ϕ ψ      
α Eps_max ε v 0 p    
fct_IDK fct_IDG fct_IDC fct_IDPb Isoft          
Kw n0 S0 U0    
Tol α v MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS baaSqaaiaadAhaaeqaaaaa@38BD@            

定義

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

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
mat_title 材料のタイトル

(文字、最大100文字)

 
ρ i 初期密度

(実数)

[ kg m 3 ]
K0 初期体積弾性率

(実数)

[ Pa ]
G0 初期せん断係数

(実数)

[ Pa ]
c0 初期材料粘着力

(実数)

[ Pa ]
Pb0 初期cap制限圧力

(実数)

[ Pa ]
ϕ 摩擦角

(実数)

[ deg ]
ψ 塑性流れ角度

(実数)

[ deg ]
α 比率:

α = p a p b

デフォルト = 0.5(実数)

 
Eps_max 最大ダイラタンシー( ρ ρ 0 1

デフォルト = -1020(実数)

 
ε v 0 p 体積塑性ひずみの初期値 3

(実数)

 
fct_IDK (オプション)体積弾性率のスケールファクター対塑性体積ひずみの関数識別子 4

(整数)

 
fct_IDG (オプション)せん断係数のスケールファクター対塑性体積ひずみの関数識別子

(整数)

 
fct_IDC (オプション)材料粘着力のスケールファクター対相当塑性ひずみの関数識別子

(整数)

 
fct_IDPb (オプション)cap制限圧力のスケールファクター対塑性体積ひずみの関数識別子

(整数)

 
Isoft cap軟化フラグ
= 0(デフォルト)
cap軟化が許可されます。
= 1
ε v p Pbが減少できないように設定されます。

(整数)

 
Kw 孔の体積弾性率(水)

(実数)

[ Pa ]
n0 初期空隙率

(実数)

 
S0 初期飽和度

(実数)

 
U0 初期孔圧力

(実数)

[ Pa ]
Tol capシフト粘性の許容値

デフォルト = 1.0E-4(実数)

 
α v MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS baaSqaaiaadAhaaeqaaaaa@38BD@ 粘性係数

デフォルト = 0.5(実数)

 

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                  kg                   m                   s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW81/1/1
LAW81
#              RHO_I
                1700
#                 K0                  G0                  c0                 PB0           
              2.83E9              1.31E9                   1                   1
#               BETA                 PSI
                  15                  10
#              ALPHA           EPS_p_max               EPS_0
                  .5                 .02                .002
#  Fct_IDK   Fct_IDG   Fct_IDc  Fct_IDPb    I_soft
         0         0         3         4         1
#                 Kw                  n0                  S0                  U0
              2.5E10                 0.1                0.99                 0.0
#                Tol             alpha_v
              0.0001                 0.5
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  3. FUNCTIONS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/3
Yield Hardening
#                  X                   Y
                   0                2000                                                            
                  .1             2002000                                                            
                   1             2002000                                                            
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/4
Cap Hardening
#                  X                   Y
                  -1                1000                                                            
                   0                1000                                                            
                .001               30000                                                            
               .0022               70000                                                            
               .0024               80000                                                            
                .004              100000                                                            
               .0056              200000                                                            
               .0078              800000                                                            
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

コメント

  1. 降伏サーフェスは、次のように定義されます:(1)
    F = q r c ( p ) ( p tan ϕ + c ) = 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOraiabg2 da9iaadghacqGHsislciGGYbWaaSbaaSqaaiaacogaaeqaaOWaaeWa aeaacaWGWbaacaGLOaGaayzkaaGaeyyXIC9aaeWaaeaacaWGWbGaci iDaiaacggacaGGUbGaeqy1dyMaey4kaSIaam4yaaGaayjkaiaawMca aiabg2da9iaaicdaaaa@4B28@

    ここで、

    r c ( p ) = 1 、右記の場合; p p a

    r c ( p ) = 1 ( p p a p b p a ) 2 、右記の場合; p a < p p b

    ここで、
    p
    圧力
    q
    フォンミーゼス応力
    c
    材料の粘着力
    P0
    圧力。ここで F p ( p 0 ) = 0
    pb
    cap制限圧力

    mat_law81
    図 1.

    この材料では、降伏サーフェスと破壊サーフェスは同じです。

  2. 塑性流れは、次に定義するように、関連付けのない潜在的な流れGによって制御されます:

    G = q p tan ψ = 0 右記の場合; p p a

    G = q tan ψ ( p ( p p a ) 2 2 ( p 0 p a ) ) = 0 、右記の場合; p a < p p 0

    G = F p > p 0 の場合、流れがcapで関連付けられるようになります。

  3. cap軟化が許可された場合、 ε v p は減少が可能になるので、それに続く曲線を適切な範囲で定義することをお奨めします。例えば、 ε v 0 p = 0 の場合は、負の値です。
  4. 体積弾性率、せん断係数、材料粘着力、cap制限圧力の各初期値は、それぞれの値ごとにスケールファクター曲線とした関数を定義することによってスケーリングできます。この関数を定義しない場合、その値は定数と見なされます。以下に例を挙げます。

    fct_IDK = 0の場合、 K = K 0

    fct_IDK ≠ 0の場合、 K = K 0 f K ( ε v p ) 、関数 f K は右記で定義; fct_IDK

  5. 初期体積弾性率と初期せん断係数は次のように計算できます。

    K 0 = E c 3 ( 1 2 ν ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbWaaS baaSqaaiaaicdaaeqaaOGaeyypa0ZaaSaaaeaacaWGfbWaaSbaaSqa aiaadogaaeqaaaGcbaGaaG4mamaabmaabaGaaGymaiabgkHiTiaaik dacqaH9oGBaiaawIcacaGLPaaaaaaaaa@4288@ G 0 = E c 2 ( 1 + ν ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGhbWaaS baaSqaaiaaicdaaeqaaOGaeyypa0ZaaSaaaeaacaWGfbWaaSbaaSqa aiaadogaaeqaaaGcbaGaaGOmamaabmaabaGaaGymaiabgUcaRiabe2 7aUbGaayjkaiaawMcaaaaaaaa@41BC@

    ここで、
    ν MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH9oGBaa a@3920@
    ポアソン比
    E c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGfbWaaS baaSqaaiaadogaaeqaaaaa@3946@
    コンクリートのヤング率
  6. 空隙率は、材料体積全体にボイドが占める体積の比率として定義します。(2)
    n =   V v o i d V t o t a l MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaamOBaiabg2da9iaacckadaWcaaWdaeaapeGaamOva8aadaWgaaWc baWdbiaadAhacaWGVbGaamyAaiaadsgaa8aabeaaaOqaa8qacaWGwb WdamaaBaaaleaapeGaamiDaiaad+gacaWG0bGaamyyaiaadYgaa8aa beaaaaaaaa@445F@
    弾性材料ではボイドの体積は変化しません。一方、塑性材料では、空隙率の変化を次のように定義します。(3)
    n=1( 1 n 0 ) e ε v p ε v0 p MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaqaaaaaaaaa Wdbiaad6gacqGH9aqpcaaIXaGaeyOeI0YaaeWaa8aabaWdbiaaigda cqGHsislcaWGUbWdamaaBaaaleaapeGaaGimaaWdaeqaaaGcpeGaay jkaiaawMcaaiaadwgapaWaaWbaaSqabeaapeGaeqyTdu2damaaDaaa meaapeGaamODaaWdaeaapeGaamiCaaaaliabgkHiTiabew7aL9aada qhaaadbaWdbiaadAhacaaIWaaapaqaa8qacaWGWbaaaaaaaaa@4A7F@
  7. 水で満たされた孔による効果:
    孔の初期状態は、初期空隙率n0、初期飽和度U0、初期孔圧力S0で定義し、次のように計算します。(4)
    n 0 =   V v o i d V t o t a l  and  S 0 =   V w a t e r V v o i d MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaamOBamaaBaaaleaacaaIWaaabeaakiabg2da9iaacckadaWcaaWd aeaapeGaamOva8aadaWgaaWcbaWdbiaadAhacaWGVbGaamyAaiaads gaa8aabeaaaOqaa8qacaWGwbWdamaaBaaaleaapeGaamiDaiaad+ga caWG0bGaamyyaiaadYgaa8aabeaaaaGcpeGaaeiiaiaabggacaqGUb GaaeizaiaabccacaWGtbWaaSbaaSqaaiaaicdaaeqaaOGaeyypa0Ja aiiOamaalaaapaqaa8qacaWGwbWdamaaBaaaleaapeGaam4Daiaadg gacaWG0bGaamyzaiaadkhaa8aabeaaaOqaa8qacaWGwbWdamaaBaaa leaapeGaamODaiaad+gacaWGPbGaamizaaWdaeqaaaaaaaa@5892@

    μ 0 >0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiVd02aaS baaSqaaiaaicdaaeqaaOGaeyOpa4JaaGimaaaa@3A5F@ であれば、S0に入力した値は使用されず、S0は再計算されます。

  8. 後処理では以下のユーザー変数を使用できます:

    USR1: 相当塑性ひずみEPSPD

    USR2: 体積塑性ひずみEPSPV

    USR3: 粘着力 c

    USR4: cap制限圧力 Pb

    USR5: 孔圧力 U

    USR6: 空隙率 n

    USR7: 飽和度 S

    USR8: capシフト u * MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaamyDa8aadaahaaWcbeqaa8qacaGGQaaaaaaa@3800@