/MAT/LAW24 (CONC)

ブロックフォーマットのキーワード この材料則は、コンクリートの脆性弾塑性挙動をモデル化します。降伏面は、Ottosen 3軸破壊面からの演繹で求められます。直交異方性損傷をモデル化し、亀裂が開いた状態と閉じた状態を扱うことができます。

オプションの組み込みモデルを使用すると、均質モデルで鉄筋を考慮できます。組み込みモデルを使用しない場合、通常は1D要素または3D要素で鉄筋がメッシュ化されます。

この材料則はソリッド要素のみと適合性があります。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MAT/LAW24/mat_ID/unit_IDまたは/MAT/CONC/mat_ID/unit_ID
mat_title
ρ i                
Ec ν Icap          
fc ft/fc fb/fc f2/fc s0/fc
Ht Dsup ε max        
ky ρ t ρ c Hbp    
α y α f vmax        
fk f0 Hv0 Eps0 hfac
E σ y Et        
α 1 α 2 α 3        

定義

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

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
mat_title 材料のタイトル

(文字、最大100文字)

 
ρ i 初期密度

(実数)

[ kg m 3 ]
Ec コンクリート弾性のヤング率

(実数)

[ Pa ]
ν ポアソン比

(実数)

 
Icap キャップ定式化フラグ 7
= 0(デフォルト)
改良されたキャップ定式化
= 1
元のキャップの定式化(Radioss 14.0以前で使用されていた初期定式化)。
= 2
更新された多軸および静水圧荷重の表現を取り入れたキャップモデル(推奨)。

(整数)

 
fc コンクリートの単軸圧縮強度

(実数)

[ Pa ]
ft/fc コンクリートの引張強度比

デフォルト = 0.10(実数)

 
fb/fc コンクリートの2軸強度比

デフォルト = 1.20(実数)

 
f2/fc コンクリートの拘束強度比

デフォルト = 4.00Icap= 0または1の場合)

デフォルト = 7.00Icap = 2の場合)

(実数)

 
s0/fc コンクリートの拘束応力比

デフォルト = 1.25(実数)

 
Ht コンクリートの引張接線係数

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

[ Pa ]
Dsup コンクリートの最大損傷

デフォルト = 0.99999(実数)

 
ε max コンクリートデータの合計破壊ひずみ

デフォルト = 1020(実数)

 
ky 硬化パラメータのコンクリート弾性初期値(第1パート)

デフォルト = 0.5(実数)

 
ρ t コンクリート弾性の破壊 / 塑性遷移圧力(第1パート)

デフォルト = 0.0(実数)

[ Pa ]
ρ c コンクリート弾性の比例降伏に遷移する圧力(第1パート)

デフォルト = -fc/3(実数)

[ Pa ]
Hbp コンクリート弾性ベースの塑性係数(第1パート)。

デフォルトはStarterによって計算されます(実数)。

[ Pa ]
α y コンクリート弾性の降伏時ダイタランシー係数(第2パート)

デフォルト = -0.2Icap= 0または2の場合)

デフォルト = 0.0Icap = 1の場合)

(実数)

 
α f コンクリート弾性の破壊時ダイタランシー係数(第2パート)

デフォルト = 0.0(実数)

 
vmax コンクリート弾性の最大体積圧縮(< 0)(第2パート)

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

 
fk 初期のキャップ開始 6

デフォルト = -fc/3(実数)

[ Pa ]
f0 初期のキャップ終了 6

デフォルト = -0.8 fcIcap = 0または1の場合)

デフォルト = 2 fcIcap = 2の場合)

(実数)

[ Pa ]
Hv0 初期3軸塑性係数

デフォルト = 0.2 Ec(実数)

[ Pa ]
Eps0 塑性硬化の参照塑性ひずみ(Icap = 2の場合のみ)

デフォルト = 0.02(実数)

 
hfac 塑性硬化のデフォルト減少係数(Icap = 2の場合のみ)。

デフォルト = 0.1(実数)

 
E 鋼材プロパティのヤング率

(実数)

[ Pa ]
σ y 降伏強さ

(実数)

[ Pa ]
Et 接線率

(実数)

[ Pa ]
α 1 方向1の補強で鉄筋断面が占める比率。

(実数)

 
α 2 方向2の補強で鉄筋断面が占める比率。

(実数)

 
α 3 方向3の補強で鉄筋断面が占める比率。

(実数)

 

例(コンクリート)

#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/CONC/1/1
concrete - air dry
#              RHO_I
               .0024  
#                E_c                  NU      Icap
               41200                  .2         0
#                 fc            ft_on_fc            fb_on_fc            f2_on_fc            s0_on_fc
                  44                   0                   0                   0                   0
#                H_t               D_sup             EPS_max
                   0                   0                   0
#                k_y                 r_t                 r_c                H_bp
                   0                   0                   0                   0
#            ALPHA_y             ALPHA_F               V_max
                -0.2                -0.1                   0
#                f_k                 f_0                H_v0                eps0               h_fac
                   0                   0                   0
#                  E             sigma_y                 E_t
                   0                   0                   0
#             ALPHA1              ALPHA2              ALPHA3
                   0                   0                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

コメント

  1. この材料則は、 ρ i , E c , ν MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdi3aaS baaSqaaiaadMgaaeqaaOGaaiilaiaadweadaWgaaWcbaGaam4yaaqa baGccaGGSaGaeqyVd4gaaa@3DDA@ および f c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaamOza8aadaWgaaWcbaWdbiaadogaa8aabeaaaaa@3839@ の4つのパラメータのみを指定して使用できます。 f t ,   f b ,   f 2 ,   s 0   MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaamOza8aadaWgaaWcbaWdbiaadshaa8aabeaak8qacaGGSaGaaeii aiaadAgapaWaaSbaaSqaa8qacaWGIbaapaqabaGcpeGaaiilaiaabc cacaWGMbWdamaaBaaaleaapeGaaGOmaaWdaeqaaOWdbiaacYcacaqG GaGaam4Ca8aadaWgaaWcbaWdbiaaicdaa8aabeaak8qacaGGGcaaaa@4408@ の各値は f c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaamOza8aadaWgaaWcbaWdbiaadogaa8aabeaaaaa@3839@ の比率で入力し、一般的なコンクリートの破壊に基づきます。
  2. この損傷モデルの各パラメータは次のとおりです。

    tensile_damage
    図 1. LAW24損傷モデルの応力ひずみ曲線. 破壊サーフェスと降伏サーフェスの変化

    f c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaamOza8aadaWgaaWcbaWdbiaadogaa8aabeaaaaa@3839@ は単軸圧縮強度です。

  3. 降伏エンベロープは、スケールファクター k ( σ m , k 0 ) を使用して破壊エンベロープから導かれます。
    (1)
    f = r k ( σ m , k 0 ) r f = 0
  4. スケールファクター k ( σ m , k 0 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaci4Aamaabm aabaGaeq4Wdm3aaSbaaSqaaiaad2gaaeqaaOGaaiilaiaadUgadaWg aaWcbaGaaGimaaqabaaakiaawIcacaGLPaaaaaa@3DEB@ は平均応力(圧力)の関数であり、 k 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4AamaaBa aaleaacaaIWaaabeaaaaa@37CC@ は硬化効果をモデル化しています。これは、弾性領域がその限界に達した状態の面に関連しています。引張では、 k ( σ m , k 0 ) = 1 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaci4Aamaabm aabaGaeq4Wdm3aaSbaaSqaaiaad2gaaeqaaOGaaiilaiaadUgadaWg aaWcbaGaaGimaaqabaaakiaawIcacaGLPaaacqGH9aqpcaaIXaaaaa@3FAC@ および降伏エンベロープと破壊エンベロープを重ねて表示しています。圧縮では、スケールファクター k ( σ m , k 0 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaci4Aamaabm aabaGaeq4Wdm3aaSbaaSqaaiaad2gaaeqaaOGaaiilaiaadUgadaWg aaWcbaGaaGimaaqabaaakiaawIcacaGLPaaaaaa@3DEB@ Icapの値(0、1、または2)によって異なります。


    図 2. スケールファクター k ( σ m , k 0 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaci4Aamaabm aabaGaeq4Wdm3aaSbaaSqaaiaad2gaaeqaaOGaaiilaiaadUgadaWg aaWcbaGaaGimaaqabaaakiaawIcacaGLPaaaaaa@3DEB@ によって破壊面から降伏面を形成
    • Icap=0または1の場合:


      図 3.
    • Icap=2(キャップあり)の場合:


      図 4.

    r = 2 J 2 = 2 3 σ V M σ m = I 1 3 は平均応力(圧力)、 I 1 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysamaaBa aaleaacaaIXaaabeaaaaa@37AB@ J 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOsamaaBa aaleaacaaIYaaabeaaaaa@37AD@ は、それぞれ第1応力不変量と第2応力不変量です。

  5. α y = α f = 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS baaSqaaiaadMhaaeqaaOGaeyypa0JaeqySde2aaSbaaSqaaiaadAga aeqaaOGaeyypa0JaaGimaaaa@3E4F@ の場合は圧縮が発生しないので、キャップ効果はありません。ダイタランシーパラメータ α y MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS baaSqaaiaadMhaaeqaaaaa@38BF@ α f MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS baaSqaaiaadMhaaeqaaaaa@38BF@ のデフォルト値は、Icap = 02の場合についてRadioss 2017で変更されています。これらのパラメータは負の値にする必要があり、推奨値は、それぞれ、-0.2と-0.1です。デフォルトのキャップ定式化はRadioss 2017で更新されています。Radioss 14.0以前で使用されていた元のキャップ定式化を使用するには、Icap=1を設定します。3軸荷重と静水圧荷重の場合は、Icap =2とした定式化の方が正確です。キャップ硬化は圧縮係数の関数です。せん断硬化係数は遷移領域で減少し、より良い安定性を確実にします。Icap =2の設定は、/ANIM/BRICK/EPSPまたは/H3D/SOLID/EPSPを使用して塑性ひずみを出力する場合に使用する必要があります。
  6. 組み込みの鉄筋モデルはオプションです。このモデルでは、硬化材料モデルで弾塑性材料を使用しています。鉄筋の断面積の比率を定義する際は、各要素が均質な挙動を示すことが前提となっています。この要素は、代表要素体積(REV)に該当する程度に十分に大きい必要があります。この均質モデルは、ほとんどの場合、大規模構造と粗いメッシュに使用します。それ以外の場合、鉄筋のモデル化には、トラス要素、スプリング要素、ビーム要素を使用できるほか、ソリッド要素も使用できます。鉄筋断面の比率 α 1 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS baaSqaaiaaigdaaeqaaaaa@387C@ α 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS baaSqaaiaaigdaaeqaaaaa@387C@ ,、および α 3 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS baaSqaaiaaigdaaeqaaaaa@387C@ を指定する必要があります。


    図 5.
  7. 直交異方性ソリッドプロパティ/PROP/TYPE6で鉄筋の方向を定義する必要があります。それ以外の場合は、局所要素座標r、s、tがそれぞれ方向1、方向2、方向3となります。ただし、Iframe = 2でIsolid = 1または2の場合を除きます。この場合の直交異方性の方向1、2、3は、time = 0における局所共回転座標r、s、tで定義されます。
  8. 軸対称解析では、 θ の方向が方向3になります。
  9. この材料則は、10節点四面体要素に適用できます。