{"id":462391,"date":"2026-08-08T12:22:21","date_gmt":"2026-08-08T12:22:21","guid":{"rendered":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/implicit-simple-restart-fails-with-cohesive-elements-elform-19-19\/"},"modified":"2026-08-08T12:22:21","modified_gmt":"2026-08-08T12:22:21","slug":"implicit-simple-restart-fails-with-cohesive-elements-elform-19-19","status":"publish","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/implicit-simple-restart-fails-with-cohesive-elements-elform-19-19\/","title":{"rendered":"Implicit simple restart fails with cohesive elements (ELFORM 19\/-19)"},"content":{"rendered":"

Hello,<\/p>\n

I created a very simple LS-DYNA implicit model to test simple restart with cohesive elements.<\/p>\n

The model runs normally in a continuous implicit analysis. However, when I perform a simple restart using a d3dump<\/code> file, the analysis fails in the first implicit step after the restart.<\/p>\n

I tested both:<\/p>\n

\n
\n
\n
\n
\n
\n
\n
\n
 <\/div>\n
\n
\n
\n
\n
\n
ELFORM = 19\nELFORM = -19<\/code><\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
 <\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
and the same problem occurs with both cohesive element formulations.<\/p>\n
I have attached the minimal .k<\/code> test model, as well as the corresponding messag<\/code> files from the continuous analysis and the restart analysis.<\/p>\n
Is this a known limitation or bug related to cohesive solid elements in implicit simple restart? If simple restart is required for an implicit analysis, is there another cohesive element formulation that is recommended?<\/p>\n
Thank you.<\/p>\n
K file:<\/span>
$# LS-DYNA Keyword file created by LS-PrePost(R) 2025 R1 (v4.12.9) - 16Oct2025<\/code>
$# Created on Jul-1-2026 (16:49:32)<\/code>
*KEYWORD   <\/code>
*CONTROL_ACCURACY<\/code>
$#     osu       inn    pidosu      iacc     exacc    srtflg<\/code>
         1         4         0         1       2.0         0<\/code>
*CONTROL_BULK_VISCOSITY<\/code>
$#      q1        q2      type     btype    tstype<\/code>
       0.0       0.0         1         0         0<\/code>
*CONTROL_ENERGY<\/code>
$#    hgen      rwen    slnten     rylen     irgen     maten     drlen     disen<\/code>
         2         2         2         2         2         1         1         1<\/code>
*CONTROL_HOURGLASS<\/code>
$#     ihq        qh<\/code>
         6      0.01<\/code>
*CONTROL_IMPLICIT_AUTO<\/code>
$#   iauto    iteopt    itewin     dtmin     dtmax     dtexp     kfail    kcycle<\/code>
         1        11         51.00000E-7     0.001       0.0         0         0<\/code>
*CONTROL_IMPLICIT_DYNAMICS<\/code>
$#   imass     gamma      beta    tdybir    tdydth    tdybur     irate     alpha<\/code>
         0       0.5      0.25       0.01.00000E281.00000E28         0       0.0<\/code>
*CONTROL_IMPLICIT_GENERAL<\/code>
$#  imflag       dt0    imform      nsbs       igs     cnstn      form    zero_v<\/code>
         1     0.001         2         1         0         0         0         0<\/code>
*CONTROL_IMPLICIT_SOLUTION<\/code>
$   NSOLVR    ILIMIT    MAXREF     DCTOL     ECTOL     RCTOL     LSTOL<\/code>
$#  nsolvr    ilimit    maxref     dctol     ectol     rctol     lstol    abstol<\/code>
        12         1        301.0000E-101.0000E-101.00000E-5       0.9       0.0<\/code>
$    DNORM    DIVERG     ISTIF   NLPRINT    NLNORM<\/code>
$#   dnorm    diverg     istif   nlprint    nlnorm   d3itctl     cpchk<\/code>
        -1         1         1         3         2         0         0<\/code>
$    DMTOL     EMTOL     RMTOL               NTTOL     NRTOL     RTTOL     RRTOL<\/code>
$#   dmtol     emtol     rmtol         -     nttol     nrtol     rttol     rrtol<\/code>
1.0000E-101.0000E-101.00000E-4                 0.0       0.0       0.0       0.0<\/code>
$   ARCCTL    ARCDIR    ARCLEN    ARCMTH    ARCDMP<\/code>
$#  arcctl    arcdir    arclen    arcmth    arcdmp    arcpsi    arcalf    arctim<\/code>
         0         0       0.0         1         2       0.0       0.0       0.0<\/code>
$#   lsmtd     lsdir      irad      srad      awgt      sred    kssize<\/code>
         6         2       0.0       0.0       0.0       0.0        10<\/code>
*CONTROL_IMPLICIT_SOLVER<\/code>
$#  lsolvr    lprint     negev     order      drcm    drcprm   autospc   autotol<\/code>
         7         1         2         0         4       0.0         1       0.0<\/code>
$#  lcpack    mtxdmp    iparm1    rparm1    rparm2         -    rparm4    rparm5<\/code>
         2         0       5001.00000E-9     0.001               240.0       0.0<\/code>
$#  emxdmp    rdcmem    absmem   isingle   iblropt         -      ispd    memlvl<\/code>
         0      0.95       0.0         0         0                   0       0.0<\/code>
*CONTROL_TERMINATION<\/code>
$#  endtim    endcyc     dtmin    endeng    endmas     nosol<\/code>
 3.260E+01         0       0.0       0.0       0.0         0<\/code>
*CONTROL_TIMESTEP<\/code>
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st<\/code>
       0.0      1e-5         0       0.0       0.0         0         0         0<\/code>
$#  dt2msf   dt2mslc     imscl    unused    unused     rmscl     emscl      ihdo<\/code>
       0.0         0         0                           0.0       0.0         0<\/code>
*DATABASE_EXTENT_BINARY<\/code>
$#   neiph     neips    maxint    strflg    sigflg    epsflg    rltflg    engflg<\/code>
        80         3         3         1         1         1         1         1<\/code>
$#  cmpflg    ieverp    beamip     dcomp      shge     stssz    n3thdt   ialemat<\/code>
         0         0         4         1         1         1         2         0<\/code>
$# nintsld   pkp_sen      sclp     hydro     msscl     therm    intout    nodout<\/code>
         8         0       1.0         0         0         0STRESS    STRESS<\/code>
$#    dtdt    resplt      <\/code>
         0         0         3<\/code>
*DATABASE_BNDOUT<\/code>
$#      dt    binary      lcur     ioopt   option1   option2   option3   option4<\/code>
      0.01         2         0         1         0         0         0         0<\/code>
*DATABASE_GLSTAT<\/code>
$#      dt    binary      lcur     ioopt<\/code>
      0.01         2         0         1<\/code>
*DATABASE_MATSUM<\/code>
$#      dt    binary      lcur     ioopt<\/code>
     0.001         2         0         1<\/code>
*DATABASE_SPCFORC<\/code>
$#      dt    binary      lcur     ioopt<\/code>
      0.01         2         0         1<\/code>
*DATABASE_BINARY_D3PLOT<\/code>
$#      dt      lcdt      beam     npltc    psetid<\/code>
     0.001         0         0         0         0<\/code>
$#   ioopt      rate    cutoff    window      type      pset<\/code>
         0       0.0       0.0       0.0         0         0<\/code>
*DATABASE_BINARY_D3DUMP<\/code>
$#      dt      lcdt      beam     npltc    psetid<\/code>
         1         0         0         0         0<\/code>
*LOAD_SEGMENT<\/code>
     10001     0.5E6         0        13        14        16        15<\/code>
*DEFINE_CURVE<\/code>
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp     lcint<\/code>
     10001         0       1.0       1.0       0.0       0.0         0         0<\/code>
$#                a1                  o1<\/code>
  0.000000000000E+00  0.000000000000E+00<\/code>
  5.000000000000E-01  1.000000000000E+00<\/code>
  6.000000000000E-01  1.000000000000E+00<\/code>
  3.260000000000E+01  1.000000000000E+00<\/code>
*BOUNDARY_PRESCRIBED_MOTION_SET<\/code>
$#    nsid       dof       vad      lcid        sf       vid     death     birth<\/code>
         2         1         2     10000       1.0         0       0.0       0.0<\/code>
*DEFINE_CURVE<\/code>
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp     lcint<\/code>
     10000         0       1.0       1.0       0.0       0.0         0         0<\/code>
$#                a1                  o1<\/code>
  0.000000000000E+00  0.000000000000E+00<\/code>
  6.000000000000E-01  0.000000000000E+00<\/code>
  1.133333333333E+00  8.000000000000E-03<\/code>
  2.200000000000E+00 -8.000000000000E-03<\/code>
  2.733333333333E+00  0.000000000000E+00<\/code>
  3.266666666667E+00  8.000000000000E-03<\/code>
  4.333333333333E+00 -8.000000000000E-03<\/code>
  4.866666666667E+00  0.000000000000E+00<\/code>
  5.400000000000E+00  8.000000000000E-03<\/code>
  6.466666666667E+00 -8.000000000000E-03<\/code>
  7.000000000000E+00  0.000000000000E+00<\/code>
  7.533333333333E+00  8.000000000000E-03<\/code>
  8.600000000000E+00 -8.000000000000E-03<\/code>
  9.133333333333E+00  0.000000000000E+00<\/code>
  9.666666666667E+00  8.000000000000E-03<\/code>
  1.073333333333E+01 -8.000000000000E-03<\/code>
  1.126666666667E+01  0.000000000000E+00<\/code>
  1.180000000000E+01  8.000000000000E-03<\/code>
  1.286666666667E+01 -8.000000000000E-03<\/code>
  1.340000000000E+01  0.000000000000E+00<\/code>
  1.393333333333E+01  8.000000000000E-03<\/code>
  1.500000000000E+01 -8.000000000000E-03<\/code>
  1.553333333333E+01  0.000000000000E+00<\/code>
  1.606666666667E+01  8.000000000000E-03<\/code>
  1.713333333333E+01 -8.000000000000E-03<\/code>
  1.766666666667E+01  0.000000000000E+00<\/code>
  1.820000000000E+01  8.000000000000E-03<\/code>
  1.926666666667E+01 -8.000000000000E-03<\/code>
  1.980000000000E+01  0.000000000000E+00<\/code>
  2.033333333333E+01  8.000000000000E-03<\/code>
  2.140000000000E+01 -8.000000000000E-03<\/code>
  2.193333333333E+01  0.000000000000E+00<\/code>
  2.246666666667E+01  8.000000000000E-03<\/code>
  2.353333333333E+01 -8.000000000000E-03<\/code>
  2.406666666667E+01  0.000000000000E+00<\/code>
  2.460000000000E+01  8.000000000000E-03<\/code>
  2.566666666667E+01 -8.000000000000E-03<\/code>
  2.620000000000E+01  0.000000000000E+00<\/code>
  2.673333333333E+01  8.000000000000E-03<\/code>
  2.780000000000E+01 -8.000000000000E-03<\/code>
  2.833333333333E+01  0.000000000000E+00<\/code>
  2.886666666667E+01  8.000000000000E-03<\/code>
  2.993333333333E+01 -8.000000000000E-03<\/code>
  3.046666666667E+01  0.000000000000E+00<\/code>
  3.100000000000E+01  8.000000000000E-03<\/code>
  3.206666666667E+01 -8.000000000000E-03<\/code>
  3.260000000000E+01  0.000000000000E+00<\/code>
*SET_NODE_LIST<\/code>
$#     sid       da1       da2       da3       da4    solver       its         -<\/code>
         2       0.0       0.0       0.0       0.0MECH               1          <\/code>
$#    nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8<\/code>
         9        10        11        12        13        14        15        16<\/code>
*BOUNDARY_SPC_SET<\/code>
$#    nsid       cid      dofx      dofy      dofz     dofrx     dofry     dofrz<\/code>
         2         0         0         1         0<\/code>
*BOUNDARY_SPC_SET<\/code>
$#    nsid       cid      dofx      dofy      dofz     dofrx     dofry     dofrz<\/code>
         1         0         1         1         1<\/code>
*SET_NODE_LIST_TITLE<\/code>
NODESET(SPC) 1<\/code>
$#     sid       da1       da2       da3       da4    solver       its         -<\/code>
         1       0.0       0.0       0.0       0.0MECH               1          <\/code>
$#    nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8<\/code>
         1         2         3         4         5         6         7         8<\/code>
*MAT_ELASTIC_TITLE<\/code>
ROCK_OUTTER<\/code>
$#     mid        ro         e        pr        da        db  not used<\/code>
         1      2700 1.700E+10      0.23         0         0         0<\/code>
*MAT_COHESIVE_ELASTIC_TITLE<\/code>
COHESIVE<\/code>
         2 1.000E-06         1         0 2.400E+09 3.000E+10 1.000E+20 1.000E+20<\/code>
*MAT_ELASTIC_TITLE<\/code>
ROCK_OUTTER<\/code>
$#     mid        ro         e        pr        da        db  not used<\/code>
         3      2700 1.700E+10      0.23         0         0         0<\/code>
*SECTION_SOLID<\/code>
$#    secid    elform<\/code>
         1         2<\/code>
*SECTION_SOLID<\/code>
$#    secid    elform<\/code>
         2        19<\/code>
*PART<\/code>
$#                                                                         title<\/code>
boxsolid<\/code>
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid<\/code>
         1         1         1         0         0         0         0         0<\/code>
*PART<\/code>
$#                                                                         title<\/code>
boxsolid<\/code>
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid<\/code>
         2         2         2         0         0         0         0         0<\/code>
*PART<\/code>
$#                                                                         title<\/code>
boxsolid<\/code>
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid<\/code>
         3         1         3         0         0         0         0         0<\/code>
*NODE<\/code>
$#   nid               x               y               z      tc      rc<\/code>
       1             0.0             0.0             0.0       0       0<\/code>
       2            0.15             0.0             0.0       0       0<\/code>
       3             0.0            0.15             0.0       0       0<\/code>
       4            0.15            0.15             0.0       0       0<\/code>
       5             0.0             0.0           0.075       0       0<\/code>
       6            0.15             0.0           0.075       0       0<\/code>
       7             0.0            0.15           0.075       0       0<\/code>
       8            0.15            0.15           0.075       0       0<\/code>
       9             0.0             0.0           0.075       0       0<\/code>
      10            0.15             0.0           0.075       0       0<\/code>
      11             0.0            0.15           0.075       0       0<\/code>
      12            0.15            0.15           0.075       0       0<\/code>
      13             0.0             0.0            0.15       0       0<\/code>
      14            0.15             0.0            0.15       0       0<\/code>
      15             0.0            0.15            0.15       0       0<\/code>
      16            0.15            0.15            0.15       0       0<\/code>
*ELEMENT_SOLID<\/code>
$#   eid     pid      n1      n2      n3      n4      n5      n6      n7      n8<\/code>
       1       1       1       2       4       3       5       6       8       7<\/code>
       2       2       5       6       8       7       9      10      12      11<\/code>
       3       3       9      10      12      11      13      14      16      15<\/code>
*END<\/code><\/p>\n
messag01 file  after simple restart:<\/span><\/p>\n
read d3dump1337 file<\/p>\n
BEGIN implicit statics step 1337 t=1.3370E+00<\/p>\n
Linear step computed
Initial translational energy = NaN
Initial total energy         = NaN
Initial residual norm        = NaN
Translational nodes norm     = NaN<\/p>\n","protected":false},"template":"","class_list":["post-462391","topic","type-topic","status-publish","hentry","topic-tag-cohesive","topic-tag-implicit","topic-tag-restart-analysis-1"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t