基于LS-DYNA的复合材料内聚力失效仿真（双臂梁，隐式求解）

以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 

计算结果和理论的对比：