以DCB復合材料雙懸臂梁實驗為研究對象，基于LS-DYNA隱式算法+內聚力單元（MAT138），給出了完整的k文件！

                                      裂紋尖端的應力云圖

                                     內聚界面的損傷演化

求解設置：

    內聚力單元采用了mat138 

    復合材料體系和幾何特征為：

    AS4/PEEK carbon fiber reinforced composite was simulated and compared with available beam

theory solutions. The specimen length, L, is given in Fig. 6 and 20.0 mm wide with two, 2h, 1.55

mm-thick arms, the latter providing a mode mixity of G II /G T = 43% for the FRMM models. The

initial delamination length is a 0 = 35 mm. The material properties of the AS4/PEEK specimen are

as follows: E 11 = 120 GPa, E 22 = E 33 = 11 GPa, ν 12 = ν 13 = 0.32, ν 23 = 0.45, G 12 = G 13 = 5.5 GPa,and G 23 = 3.7 GPa. The properties of the interface are given in Table 3. The average size of thecontinuum and cohesive elements in the direction of crack propagation is 0.5 mm

隱式求解設置：

*CONTROL_ENERGY

$#    hgen      rwen    slnten     rylen     

         2         2         2         2

*CONTROL_TIMESTEP

$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st

       0.0       0.9         0       0.0       0.0         0         0         0

$#  dt2msf   dt2mslc     imscl    unused    unused     rmscl     

       0.0         0         0                           0.0

$

*CONTROL_IMPLICIT_GENERAL

$#  imflag       dt0    imform     nsbs       igs     cnstn      form      

         1      0.01

*CONTROL_IMPLICIT_SOLUTION

$  nlsolvr    ilimit    maxref     dctol     ectol     rctol     lstol

         0         0         0      0.00      0.00         0         0

$    dnorm   divflag   inistif   nlprint

         1         0         0         0

$

*CONTROL_IMPLICIT_AUTO

$    iauto    iteopt    itewin     dtmin     dtmax

         1       50         20      1E-6      0.02

$

*CONTROL_IMPLICIT_STABILIZATION

$      ias     scale    tstart      tend

         1      0.00         0         0

$

*DATABASE_BINARY_D3PLOT

0.5000E-02

$

*CONTROL_TERMINATION

$$  ENDTIM    ENDCYC     DTMIN    ENDENG    ENDMAS     NOSOL

         1         0       0.0       0.0       0.0 

計算結果和理論的對比：