/FAIL/ALTER

破壊応力は、微小亀裂および亀裂伝播速度を定義するパラメータによって表されます。X-FEMでは、応力は亀裂方向に垂直で0に設定されます。

フォーマット

定義

例

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
#              MUNIT               LUNIT               TUNIT
                  Mg                  mm                   s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW36/72200103/1
Glass with linear hardening
#              RHO_I
             2.50E-9
#                  E                  NU           Eps_p_max               Eps_t               Eps_m
             70000.0                0.23                   0                   0                   0 
#  N_funct  F_smooth              C_hard               F_cut               Eps_f
         1         1                   0                1650                   0
#  fct_IDp              Fscale   fct_IDE                EInf                  CE
         0                   0         0                   0                   0
#  fct_ID1   fct_ID2   fct_ID3   fct_ID4   fct_ID5
 722001021
#           Fscale_1            Fscale_2            Fscale_3            Fscale_4            Fscale_5
                1000
#          Eps_dot_1           Eps_dot_2           Eps_dot_3           Eps_dot_4           Eps_dot_5
                   0 
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FAIL/ALTER/72200103/1
#              EXP_N                  V0                  VC   NCYCLES     IRATE     ISIDE      MODE
                16.0                 6.0             1520000         6         0         0         1
#            CR_FOIL              CR_AIR             CR_CORE             CR_EDGE    GRSH4N    GRSH3N
             0.00040             0.00100             0.00500                   0         0         0
#                KIC                 KTH                RLEN                TDEL      Iout
              23.717              7.9057                 1.0                   0         0
#              KRES1               KRES2
                   0                   0
#               ETA1               BETA1                TAU1            AREA_REF
                   0                   0                   0                   0
#               ETA2               BETA2                TAU2
                   0                   0                   0
#               SIG0             P_SCALE  P_SWITCH
                   0                   0         0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/722001021
Function for glass
#                  X                   Y
                 0.0               500.0
                 1.0               550.0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata

コメント

1. この破壊基準は、破壊の基準として最大応力を用いており、初期亀裂と亀裂進展速度によって決定される材料の強度に基づいて計算されます。modeスイッチフラグに応じて、近傍の要素間の異なる破壊進展モデルが使用され得ます。
2. Irate=0である際、指数移動平均フィルターが使用され、フィルタリングされた応力は：(1)
   $${\sigma }_f\left(t\right)=\alpha \sigma \left(t\right)+\left(1-\alpha \right)\sigma \left(t-\text{Δ}t\right)$$
   ここで、

   ${\sigma }_f=\text{filtered stress}$

   $\alpha =\frac{2}{N_{cycles}+1}$

3. この破壊モデルは低減積分シェル要素（I_shell =24およびI_sh3n =2を推奨）とのみ適合性があり、完全積分シェル要素とは適合性を有しません。また、使用されるシェルプロパティに制限はないものの、1層シェルモデルとのみ適合します。
4. グループgrsh4Nおよびgrsh3N内で定義された要素は、ウィンドシールドのエッジに沿っていなくてはならず、特定の破壊脆弱を受けます。
5. この破壊モデルは、重複節点を用いてポリビニルブチラール（PVB）ソリッド要素の層を挟むシェル要素に適用されます。アセンブリ全体がウィンドシールドをモデル化します。

   
   
   図 1. ウィンドシールドの有限要素モデル

   
   
   図 2. ウィンドシールドモデル - アセンブリ全体

6. この破壊モデルを用いたシェル要素は、法線が中間PVBとは別の方向を指すよう方向付けされていなくてはなりません。
7. シェル要素には、曲げを正しくモデル化するために、オフセットが適用される必要があります。これは、/PROP/TYPE51 I_pos=4を用いて行うことができます。
8. 破壊限界は、その位置および周りの要素の破壊状態に依存します。 ¹
9. fail_IDは、/STATE/BRICK/FAILと/INIBRI/FAILおよび/PERTURB/FAIL/BIQUADと共に使用されます。デフォルト値はありません。この行が空白の場合、/INIBRI/FAIL内の破壊モデル変数のために出力される値はありません（3次元ソリッドの場合は/STATE/BRICK/FAIL、シェルの場合は/STATE/SHELL/FAILで.staファイルに書き込まれます）。
10. Ch. Brokmann拡張²は、外部ガラスサーフェスのみの追加の破壊基準を定義します。これにより、機械的処理または化学的処理が原因で、ガラスサーフェスに初期応力が発生します。ガラスサーフェスの微小欠陥の統計的評価により、左切り捨てワイブル確率分布を使用して破壊の確率を定義できます。
    ここで、 $i$ = 1、2 （下部および上部のサーフェスの場合） (2)

    $$P(\sigma )=1-\exp \left[{{\displaystyle \left(\frac{{\tau }_i}{{\eta }_i}\right)}}^{{\beta }_i}-{{\displaystyle \left(\frac{\sigma }{{\eta }_i}\right)}}^{{\beta }_i}\right]$$
    切り捨て点 $\tau =0$ では、よく知られている2パラメータワイブル分布が得られます。Brokmannモデルは、ガラス内のランダムに配向された初期欠陥を計算し、それらをさまざまな長さと形状を持つすべての有限要素に配分します。亀裂の成長は、次の微分方程式で表すことができます： (3)

    $$da=V_0.{{\displaystyle \left(\frac{Y\sigma \sqrt{\pi a}}{{{\displaystyle K}}_{IC}}\right)}}^{Exp\_n}dt$$

    ここで、 $Y$ はワイブル分布を使用して求めた欠陥形状係数です。

    式 3を積分すると、応力速度に強く依存する実際の亀裂サイズが得られます。実際の応力拡大係数を計算して、破壊基準で使用することができます。

    Brokmannのモデルの面白さは、分布パラメータと破壊応力値によって、破壊の確率的な可能性を推定できることです。

    
    
    図 3.

    ガラスの欠陥をランダムに初期化して十分な回数のシミュレーションを実行した後、頭部損傷基準（HIC）の特定の値に到達する確率を推定する可能性を返します。

    
    
    図 4.

11. Ch. Brokmann基準を使用すると、フラグIrateは自動的に0に設定されます。この場合、指数平均を使用して、応力フィルタリング区間のサイクル数を定義する必要があります。