Image

Ansys LS-DYNA 2025 R1 Update

More Details
ansys transportation summit

Submit your abstract for the Ansys Transportation Summit 2025 including the European LS-DYNA Conference

More Details
Image

Download the presentations of the German LS-DYNA Forum 2024

More Details

ImageCommunity

LS-DYNA 2025 R1 (v251) Windows — User-Defined Material Package
<p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">Hi,</p><p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">I have Ansys LS-DYNA 2025 R1 (v251) installed on Windows, with solver version SMP R14.1.1 double precision.</p><p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">I am implementing a custom nonlinear constitutive model (resilient modulus k-theta) using a UMAT subroutine and need the Windows usermat package containing:</p><ul class="[li_&]:mb-0 [li_&]:mt-1 [li_&]:gap-1 [&:not(:last-child)_ul]:pb-1 [&:not(:last-child)_ol]:pb-1 list-disc flex flex-col gap-1 pl-8 mb-3"><li class="font-claude-response-body whitespace-normal break-words pl-2">dyn21.F / dyn21b.F<li class="font-claude-response-body whitespace-normal break-words pl-2">include files (.inc)<li class="font-claude-response-body whitespace-normal break-words pl-2">static libraries (.lib)<li class="font-claude-response-body whitespace-normal break-words pl-2">Windows build script (nmake makefile)<p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">Could you please advise on the correct procedure to obtain this package for my installation?</p><p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">Compiler environment: Intel oneAPI (ifort), Microsoft Visual C++ 2019.</p><p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">Thank you very much.</p><p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">Eloísa Castilho </p>
Avatar
emsantos@lnec.pt June 8, 2026
Bonded contact doesn’t hold.
<p>Hi everyone, I'm striking a rigid wall. LS-Dyna Workbench.<br>The impact is on the bolts. The bolt and nut have bonded contact. But when the impact occurs, the bolt slips through the nut.</p><p></p><p>I'm also getting penetration messages:</p><div id="paraDiv" class="text-color-primary" style="display: inline;"><p id="paragraphData" class="paragraphText text-color-primary" style="font-size: 12px; display: inline;">23098 initial penetrations were found for interface        763<br>     57 initial penetrations were found for interface        646<br>     70 initial penetrations were found for interface        649<br>     78 initial penetrations were found for interface        652<br>     32 initial penetrations were found for interface        655<br>      62 initial penetrations were found for interface        847<br>     27 initial penetrations were found for interface        658<br>   14564 initial penetrations were found for interface        712<br>      60 initial penetrations were found for interface        850<br>      50 initial penetrations were found for interface        865<br>      60 initial penetrations were found for interface        880<br>   11442 initial penetrations were found for interface        730<br>      60 initial penetrations were found for interface        841<br>   14557 initial penetrations were found for interface        718<br>      60 initial penetrations were found for interface        844<br>      48 initial penetrations were found for interface        859<br>   14514 initial penetrations were found for interface        724<br>      48 initial penetrations were found for interface        862<br>   14504 initial penetrations were found for interface        706<br>    122 initial penetrations were found for interface        673<br>     741 initial penetrations were found for interface        823<br>    107 initial penetrations were found for interface        676<br>     551 initial penetrations were found for interface        838</p></div><p>Help please.</p>
Avatar
javat33489 June 8, 2026
Define user material that calls existing LS-DYNA material subroutines
<p><p>I have read the user manual on user defined materials but it doens't address the following use case so I thought I would post here. </p><p>Is it possible define a user material that calls a predefined LS-DYNA material model? If so what would the syntax look like to accomplish this? </p><p>For example let's say I want to introduce a topology optimization variable as part of a user defined material model and I want to use a predefined elastoplatic model which already exists in LS-DYNA in combination with said TO variable. I am trying to take advantage of the material model that already exists rather than reimplementing it myself. </p><p>Any help or insight would be greatly appreciated. Thanks in advance. </p></p><p><p> A small update. I remembered being able to view the signatures of object files .o or .a files with command line tools like nm or objdump. So as a small first step I can at least view the function names available in precompiled object files. Which I believe correspond to the material subroutine names. However, this is not the complete function definition with variable names in the correct order but it is a start. It is also probably reasonable to assume that the variable ordering in the card definitions in the user manual are the same order as the variables in the function call. Although it will probably get a little tricky with cards that can be repeated and other such options. But it is a start...</p></p>
Avatar
mwhitte6@jhu.edu June 5, 2026
Sensor_Define_Element_Set comp direction not being recognized
<p>Hello, I am getting an error as shown below when trying to use sensor_define_element_set.    in the below parameter, if I set comp = PRIN1, below error shows up.   if I set it to XX, YY, ZZ, it doesn't</p><p>I have database_elout, database_history_solid and control_output all properly set.   If anyone can think of anything causing this, I would realy appreciate it.  </p><p> </p><p></p>
Avatar
Dennis Chen June 4, 2026
Explicit Dynamics not behaving as expected
<p>Hello all, I'm hoping to find some help here. Hopefully I am just misunderstanding the physics of what i going on.</p><p>So i set up an Explicit Dynamics sim of a rotating shaft coupling impacting a guard. The axis of the coupling is angled 45deg off the YZ plane and 30deg off the XY plane. Coupling is set to rotate around the coupling center axis at 2000rpm and translate 2000in/sec in the global Y direction. Or so i think. But the solution is showing the coupling rotating as expected, but translational movement is mostly in the X global direction with some translation in the intended Y global direction. </p><p>Attached are some pictures of the setup. I'm sure there will be questions and more information i have not thought to answer and provide, so please ask.</p><p>My question is why is the coupling translating mostly in X instead of Y as intended? Is it because i am misunderstanding the effects of angular momentum?</p><p></p><p></p><p></p><p></p><p> </p>
Avatar
david.mccul@gmail.com June 3, 2026
Cutter–Rock Cutting Model Stalls with Eroding Surface-to-Surface Contact
<p>I am modeling a rigid/inclined PDC cutter cutting into a rock block in LS-DYNA using an RHT-type rock material, MAT_ADD_EROSION, and CONTACT_ERODING_SURFACE_TO_SURFACE. The cutter motion is prescribed only in the positive x-direction, while the rock is constrained at the base and, in some trials, subjected to lateral confining pressure of 40 MPa. The model initializes without fatal errors in many runs, but it stalls very early or shortly after contact/erosion starts; CPU usage drops significantly, the time step becomes very small, and failed rock elements are either not visibly deleted or only a few elements erode before the simulation stops progressing.</p><p>I have tried changing erosion/contact parameters such as CON1/CON2, ISYM, EROSOP, IADJ, commenting out the confining pressure, and running in double precision, but the same stalling behavior persists. I suspect the issue may be related to the eroding contact definition, contact update after element deletion, conflicting contact/control keywords, or inappropriate RHT/erosion parameters. I would appreciate advice on the required/typical keyword setup for rock cutting with CONTACT_ERODING_SURFACE_TO_SURFACE, especially regarding contact cards, CONTROL_CONTACT, erosion deletion/update settings, and stable parameter ranges.</p>
Avatar
eliasedwarde@tamu.edu June 2, 2026
Hyperelastic material – spurious stress after rigid-to-deformable switch
<p class="ds-markdown-paragraph"><span class="">Solver Version:</span></p><p class="ds-markdown-paragraph"><span class="">Version : smp s R14                           <br>Revision: R14.1.1-8-g2ef819c24f <br>AnLicVer: 2025 R1 (20241030+dl-67-g73b144a) </span></p><p class="ds-markdown-paragraph"><span class="">Description:</span></p><p class="ds-markdown-paragraph"><span class="">I've encountered an issue where hyperelastic materials develop strange initial stresses (solver: SMP single precision, double precison solver will crash the simulation immediately) when a part is switched from rigid to deformable using </span>*DEFORMABLE_TO_RIGID_AUTOMATIC<span class="">. The same setup works perfectly with elastic materials.</span></p><p class="ds-markdown-paragraph"><span class="">Summary of behavior:</span></p><p class="ds-markdown-paragraph"><span class="">A part with hyperelastic material is initially rigid, then automatically switched to deformable at t=10.0</span></p><p class="ds-markdown-paragraph"><span class="">Immediately after the switch, the part shows spurious stresses</span></p><p class="ds-markdown-paragraph"><span class="">If I replace the hyperelastic material with an elastic material, the problem disappears</span></p><p class="ds-markdown-paragraph"><span class="">The same issue occurs whether using full restart with </span>*STRESS_INITIALIZATION<span class=""> to add a new hyperelastic part, or within a single continuous simulation using </span>*DEFORMABLE_TO_RIGID_AUTOMATIC</p><p class="ds-markdown-paragraph"><span class="">Attached test case:</span><br><span class="">A minimal working example demonstrating the issue. One single element uses </span>*DEFORMABLE_TO_RIGID_AUTOMATIC<span class=""> to switch from rigid to deformable at t=10.0 with hyperelastic rubber material. Simply swap the material to elastic (</span>*MAT_ELASTIC<span class="">) to see the problem disappear.</span></p><p class="ds-markdown-paragraph"><span class="">Questions:</span></p><p class="ds-markdown-paragraph"><span class="">Is this a known limitation/bug in R14?</span></p><p class="ds-markdown-paragraph"><span class="">Any insights would be greatly appreciated.</span></p><p class="ds-markdown-paragraph"><span class=""></span></p><p class="ds-markdown-paragraph"><span class="">Key file to reproduce:</span></p><div><div><span style="font-size: 12pt;">*KEYWORD</span></div><div><span style="font-size: 12pt;">*TITLE</span></div><div><span style="font-size: 12pt;">Hyperelastic part - deformable/rigid automatic switch</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">$  Units: mm, kg, ms, kN, GPa, kN-mm</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">*CONTROL_TERMINATION</span></div><div><span style="font-size: 12pt;">$   endtim    endcyc     dtmin    endeng    endmas</span></div><div><span style="font-size: 12pt;">      13.0         0       0.0       0.0       0.0</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">*DEFINE_CURVE</span></div><div><span style="font-size: 12pt;">$ LCID    SIDR    SFA     SFO     OFFA    OFFO    DATTYP  LCINT</span></div><div><span style="font-size: 12pt;">      1005         0       1.0       1.0       0.0       0.0         0         0</span></div><div><span style="font-size: 12pt;">$ time    timestep</span></div><div><span style="font-size: 12pt;">                 0.0                1e-2</span></div><div><span style="font-size: 12pt;">                 10.0               1e-2</span></div><div><span style="font-size: 12pt;">*CONTROL_TIMESTEP</span></div><div><span style="font-size: 12pt;">$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st</span></div><div><span style="font-size: 12pt;">      1e-2       1.0         0       0.0         0      1005         0         0</span></div><div><span style="font-size: 12pt;">*CONTROL_ENERGY</span></div><div><span style="font-size: 12pt;">$     hgen      rwen    slnten     rylen</span></div><div><span style="font-size: 12pt;">         2         2</span></div><div><span style="font-size: 12pt;">*CONTROL_OUTPUT</span></div><div><span style="font-size: 12pt;">$    npopt    neecho    nrefup    iaccop     opifs    ipnint    ikedit</span></div><div><span style="font-size: 12pt;">         1         3</span></div><div><span style="font-size: 12pt;">*CONTROL_SHELL</span></div><div><span style="font-size: 12pt;">$   wrpang    itrist     irnxx    istupd    theory       bwc     miter</span></div><div><span style="font-size: 12pt;">                                       1         2</span></div><div><span style="font-size: 12pt;">*DATABASE_BINARY_D3PLOT</span></div><div><span style="font-size: 12pt;">$       dt      lcdt</span></div><div><span style="font-size: 12pt;">      1.00</span></div><div><span style="font-size: 12pt;">*DATABASE_EXTENT_BINARY</span></div><div><span style="font-size: 12pt;">$    neiph     neips    maxint    strflg    sigflg    epsflg    rltflg    engflg</span></div><br><div><span style="font-size: 12pt;">$   cmpflg    ieverp    beamip</span></div><div><span style="font-size: 12pt;">                   1</span></div><div><span style="font-size: 12pt;">*DATABASE_BINARY_D3THDT</span></div><div><span style="font-size: 12pt;">$       dt      lcdt</span></div><div><span style="font-size: 12pt;">    999999</span></div><div><span style="font-size: 12pt;">*DATABASE_GLSTAT</span></div><div><span style="font-size: 12pt;">$       dt</span></div><div><span style="font-size: 12pt;">      0.10</span></div><div><span style="font-size: 12pt;">*DATABASE_MATSUM</span></div><div><span style="font-size: 12pt;">$       dt</span></div><div><span style="font-size: 12pt;">      0.10</span></div><div><span style="font-size: 12pt;">*DATABASE_NODOUT</span></div><div><span style="font-size: 12pt;">$       dt</span></div><div><span style="font-size: 12pt;">      0.10</span></div><div><span style="font-size: 12pt;">*DATABASE_RBDOUT</span></div><div><span style="font-size: 12pt;">$       dt</span></div><div><span style="font-size: 12pt;">      0.10</span></div><div><span style="font-size: 12pt;">*PART</span></div><div><span style="font-size: 12pt;">$      pid       sid       mid     eosid      hgid    adpopt</span></div><div><span style="font-size: 12pt;">part1</span></div><div><span style="font-size: 12pt;">         1         2         2</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">$  Aluminum</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">*MAT_PLASTIC_KINEMATIC</span></div><div><span style="font-size: 12pt;">$      mid        ro         e        pr      sigy      etan      beta</span></div><div><span style="font-size: 12pt;">         1   2.70e-6      68.9     0.330     0.286   0.00689</span></div><div><span style="font-size: 12pt;">$      src       srp        fs</span></div><br><div><span style="font-size: 12pt;">*SECTION_SHELL</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">$      sid    elform      shrf       nip     propt   qr/irid     icomp</span></div><div><span style="font-size: 12pt;">         1         2</span></div><div><span style="font-size: 12pt;">$       t1        t2        t3        t4      nloc</span></div><div><span style="font-size: 12pt;">       1.0       1.0       1.0       1.0</span></div><div><span style="font-size: 12pt;">*SECTION_SHELL</span></div><div><span style="font-size: 12pt;">$       ID    elform      shrf       nip     propt        qr     icomp     setyp</span></div><div><span style="font-size: 12pt;">         2         2  0.833333         3         0         0         0         0</span></div><div><span style="font-size: 12pt;">$       t1        t2        t3        t4      nloc     marea              unused</span></div><div><span style="font-size: 12pt;">       0.1       0.1       0.1       0.1         0         0                    </span></div><div><span style="font-size: 12pt;">*MAT_HYPERELASTIC_RUBBER</span></div><div><span style="font-size: 12pt;">$       ID        ro        pr         N        Nv         G      sigf       Ref</span></div><div><span style="font-size: 12pt;">         2 8.699E-07 0.4999172         0         0         0         0         0</span></div><div><span style="font-size: 12pt;">$      C10       C01       C11       C20       C02       C30             unused1</span></div><div><span style="font-size: 12pt;">   0.11e-3  0.028e-3         0         0         0         0                    </span></div><div><span style="font-size: 12pt;">*NODE</span></div><div><span style="font-size: 12pt;">$   node               x               y               z      tc      rc</span></div><div><span style="font-size: 12pt;">       1             0.0             0.0             0.0</span></div><div><span style="font-size: 12pt;">       2             1.0             0.0             0.0</span></div><div><span style="font-size: 12pt;">       3             1.0             1.0             0.0</span></div><div><span style="font-size: 12pt;">       4             0.0             1.0             0.0</span></div><div><span style="font-size: 12pt;">$$$$  Define Elements</span></div><div><span style="font-size: 12pt;">*ELEMENT_SHELL</span></div><div><span style="font-size: 12pt;">$    eid     pid      n1      n2      n3      n4</span></div><div><span style="font-size: 12pt;">       1       1       1       2       3       4</span></div><div><span style="font-size: 12pt;">*DEFORMABLE_TO_RIGID</span></div><div><span style="font-size: 12pt;">         1</span></div><div><span style="font-size: 12pt;">*DEFORMABLE_TO_RIGID_AUTOMATIC</span></div><div><span style="font-size: 12pt;">$#   SWSET      CODE     TIME1     TIME2     TIME3     ENTNO     RELSW    PAIRED</span></div><div><span style="font-size: 12pt;">         1         0      10.0      1.0e20     0.0         0         0         0</span></div><div><span style="font-size: 12pt;">$#    NRBF      NCSF       RWF     DTMAX       D2R       R2D    OFFSET</span></div><div><span style="font-size: 12pt;">         0         0         0       0.0         0         1       0.0</span></div><div><span style="font-size: 12pt;">$# Rigid to deformable card: one card because R2D = 1</span></div><div><span style="font-size: 12pt;">$#     PID     PTYPE</span></div><div><span style="font-size: 12pt;">         1      PART</span></div><div><span style="font-size: 12pt;">*END</span></div></div>
Avatar
maoy@umsystem.edu May 29, 2026
How to model a crack using MAT_ADD_EROSION in Ansys LS-Dyna?
I wish to use the MAT_ADD_EROSION model in LS-DYNA to capture material failure and erosion effectively. Are there any resources to help me get started? Thanks!    
Avatar
innovation customer May 29, 2026
Simulating a prestrained material
Is it possible to simulate a pre-strained material directly without any pre-simulation? We have conducted uniaxial prestraining on a material and would now like to run another simulation. The major and minor strains from our previous operation are available to us. Thanks!
Avatar
innovation customer May 29, 2026
Question Regarding Applicability of YLD2004-18p for Shell Elements
<p> Dear LS-DYNA Support Team, I hope you are doing well.</p><p>I would like to ask whether the yield locus model **YLD2004-18p** can be applied to **shell elements** in LS-DYNA.</p><p>Could you please let me know if this material model is supported for shell formulations, and if there are any specific limitations or requirements that need to be considered?</p><p>Thank you very much for your support.</p><p> Best regards,</p><p> Philipp</p>
Avatar
pneufeld@hs-koblenz.de May 29, 2026
LS-DYNA 2025 R1 (v251) Windows — User-Defined Material Package
<p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">Hi,</p><p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">I have Ansys LS-DYNA 2025 R1 (v251) installed on Windows, with solver version SMP R14.1.1 double precision.</p><p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">I am implementing a custom nonlinear constitutive model (resilient modulus k-theta) using a UMAT subroutine and need the Windows usermat package containing:</p><ul class="[li_&]:mb-0 [li_&]:mt-1 [li_&]:gap-1 [&:not(:last-child)_ul]:pb-1 [&:not(:last-child)_ol]:pb-1 list-disc flex flex-col gap-1 pl-8 mb-3"><li class="font-claude-response-body whitespace-normal break-words pl-2">dyn21.F / dyn21b.F<li class="font-claude-response-body whitespace-normal break-words pl-2">include files (.inc)<li class="font-claude-response-body whitespace-normal break-words pl-2">static libraries (.lib)<li class="font-claude-response-body whitespace-normal break-words pl-2">Windows build script (nmake makefile)<p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">Could you please advise on the correct procedure to obtain this package for my installation?</p><p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">Compiler environment: Intel oneAPI (ifort), Microsoft Visual C++ 2019.</p><p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">Thank you very much.</p><p class="font-claude-response-body break-words whitespace-normal leading-[1.7]">Eloísa Castilho </p>
Avatar
emsantos@lnec.pt June 8, 2026
Bonded contact doesn’t hold.
<p>Hi everyone, I'm striking a rigid wall. LS-Dyna Workbench.<br>The impact is on the bolts. The bolt and nut have bonded contact. But when the impact occurs, the bolt slips through the nut.</p><p></p><p>I'm also getting penetration messages:</p><div id="paraDiv" class="text-color-primary" style="display: inline;"><p id="paragraphData" class="paragraphText text-color-primary" style="font-size: 12px; display: inline;">23098 initial penetrations were found for interface        763<br>     57 initial penetrations were found for interface        646<br>     70 initial penetrations were found for interface        649<br>     78 initial penetrations were found for interface        652<br>     32 initial penetrations were found for interface        655<br>      62 initial penetrations were found for interface        847<br>     27 initial penetrations were found for interface        658<br>   14564 initial penetrations were found for interface        712<br>      60 initial penetrations were found for interface        850<br>      50 initial penetrations were found for interface        865<br>      60 initial penetrations were found for interface        880<br>   11442 initial penetrations were found for interface        730<br>      60 initial penetrations were found for interface        841<br>   14557 initial penetrations were found for interface        718<br>      60 initial penetrations were found for interface        844<br>      48 initial penetrations were found for interface        859<br>   14514 initial penetrations were found for interface        724<br>      48 initial penetrations were found for interface        862<br>   14504 initial penetrations were found for interface        706<br>    122 initial penetrations were found for interface        673<br>     741 initial penetrations were found for interface        823<br>    107 initial penetrations were found for interface        676<br>     551 initial penetrations were found for interface        838</p></div><p>Help please.</p>
Avatar
javat33489 June 8, 2026
Define user material that calls existing LS-DYNA material subroutines
<p><p>I have read the user manual on user defined materials but it doens't address the following use case so I thought I would post here. </p><p>Is it possible define a user material that calls a predefined LS-DYNA material model? If so what would the syntax look like to accomplish this? </p><p>For example let's say I want to introduce a topology optimization variable as part of a user defined material model and I want to use a predefined elastoplatic model which already exists in LS-DYNA in combination with said TO variable. I am trying to take advantage of the material model that already exists rather than reimplementing it myself. </p><p>Any help or insight would be greatly appreciated. Thanks in advance. </p></p><p><p> A small update. I remembered being able to view the signatures of object files .o or .a files with command line tools like nm or objdump. So as a small first step I can at least view the function names available in precompiled object files. Which I believe correspond to the material subroutine names. However, this is not the complete function definition with variable names in the correct order but it is a start. It is also probably reasonable to assume that the variable ordering in the card definitions in the user manual are the same order as the variables in the function call. Although it will probably get a little tricky with cards that can be repeated and other such options. But it is a start...</p></p>
Avatar
mwhitte6@jhu.edu June 5, 2026
Sensor_Define_Element_Set comp direction not being recognized
<p>Hello, I am getting an error as shown below when trying to use sensor_define_element_set.    in the below parameter, if I set comp = PRIN1, below error shows up.   if I set it to XX, YY, ZZ, it doesn't</p><p>I have database_elout, database_history_solid and control_output all properly set.   If anyone can think of anything causing this, I would realy appreciate it.  </p><p> </p><p></p>
Avatar
Dennis Chen June 4, 2026
Explicit Dynamics not behaving as expected
<p>Hello all, I'm hoping to find some help here. Hopefully I am just misunderstanding the physics of what i going on.</p><p>So i set up an Explicit Dynamics sim of a rotating shaft coupling impacting a guard. The axis of the coupling is angled 45deg off the YZ plane and 30deg off the XY plane. Coupling is set to rotate around the coupling center axis at 2000rpm and translate 2000in/sec in the global Y direction. Or so i think. But the solution is showing the coupling rotating as expected, but translational movement is mostly in the X global direction with some translation in the intended Y global direction. </p><p>Attached are some pictures of the setup. I'm sure there will be questions and more information i have not thought to answer and provide, so please ask.</p><p>My question is why is the coupling translating mostly in X instead of Y as intended? Is it because i am misunderstanding the effects of angular momentum?</p><p></p><p></p><p></p><p></p><p> </p>
Avatar
david.mccul@gmail.com June 3, 2026
Cutter–Rock Cutting Model Stalls with Eroding Surface-to-Surface Contact
<p>I am modeling a rigid/inclined PDC cutter cutting into a rock block in LS-DYNA using an RHT-type rock material, MAT_ADD_EROSION, and CONTACT_ERODING_SURFACE_TO_SURFACE. The cutter motion is prescribed only in the positive x-direction, while the rock is constrained at the base and, in some trials, subjected to lateral confining pressure of 40 MPa. The model initializes without fatal errors in many runs, but it stalls very early or shortly after contact/erosion starts; CPU usage drops significantly, the time step becomes very small, and failed rock elements are either not visibly deleted or only a few elements erode before the simulation stops progressing.</p><p>I have tried changing erosion/contact parameters such as CON1/CON2, ISYM, EROSOP, IADJ, commenting out the confining pressure, and running in double precision, but the same stalling behavior persists. I suspect the issue may be related to the eroding contact definition, contact update after element deletion, conflicting contact/control keywords, or inappropriate RHT/erosion parameters. I would appreciate advice on the required/typical keyword setup for rock cutting with CONTACT_ERODING_SURFACE_TO_SURFACE, especially regarding contact cards, CONTROL_CONTACT, erosion deletion/update settings, and stable parameter ranges.</p>
Avatar
eliasedwarde@tamu.edu June 2, 2026
Hyperelastic material – spurious stress after rigid-to-deformable switch
<p class="ds-markdown-paragraph"><span class="">Solver Version:</span></p><p class="ds-markdown-paragraph"><span class="">Version : smp s R14                           <br>Revision: R14.1.1-8-g2ef819c24f <br>AnLicVer: 2025 R1 (20241030+dl-67-g73b144a) </span></p><p class="ds-markdown-paragraph"><span class="">Description:</span></p><p class="ds-markdown-paragraph"><span class="">I've encountered an issue where hyperelastic materials develop strange initial stresses (solver: SMP single precision, double precison solver will crash the simulation immediately) when a part is switched from rigid to deformable using </span>*DEFORMABLE_TO_RIGID_AUTOMATIC<span class="">. The same setup works perfectly with elastic materials.</span></p><p class="ds-markdown-paragraph"><span class="">Summary of behavior:</span></p><p class="ds-markdown-paragraph"><span class="">A part with hyperelastic material is initially rigid, then automatically switched to deformable at t=10.0</span></p><p class="ds-markdown-paragraph"><span class="">Immediately after the switch, the part shows spurious stresses</span></p><p class="ds-markdown-paragraph"><span class="">If I replace the hyperelastic material with an elastic material, the problem disappears</span></p><p class="ds-markdown-paragraph"><span class="">The same issue occurs whether using full restart with </span>*STRESS_INITIALIZATION<span class=""> to add a new hyperelastic part, or within a single continuous simulation using </span>*DEFORMABLE_TO_RIGID_AUTOMATIC</p><p class="ds-markdown-paragraph"><span class="">Attached test case:</span><br><span class="">A minimal working example demonstrating the issue. One single element uses </span>*DEFORMABLE_TO_RIGID_AUTOMATIC<span class=""> to switch from rigid to deformable at t=10.0 with hyperelastic rubber material. Simply swap the material to elastic (</span>*MAT_ELASTIC<span class="">) to see the problem disappear.</span></p><p class="ds-markdown-paragraph"><span class="">Questions:</span></p><p class="ds-markdown-paragraph"><span class="">Is this a known limitation/bug in R14?</span></p><p class="ds-markdown-paragraph"><span class="">Any insights would be greatly appreciated.</span></p><p class="ds-markdown-paragraph"><span class=""></span></p><p class="ds-markdown-paragraph"><span class="">Key file to reproduce:</span></p><div><div><span style="font-size: 12pt;">*KEYWORD</span></div><div><span style="font-size: 12pt;">*TITLE</span></div><div><span style="font-size: 12pt;">Hyperelastic part - deformable/rigid automatic switch</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">$  Units: mm, kg, ms, kN, GPa, kN-mm</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">*CONTROL_TERMINATION</span></div><div><span style="font-size: 12pt;">$   endtim    endcyc     dtmin    endeng    endmas</span></div><div><span style="font-size: 12pt;">      13.0         0       0.0       0.0       0.0</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">*DEFINE_CURVE</span></div><div><span style="font-size: 12pt;">$ LCID    SIDR    SFA     SFO     OFFA    OFFO    DATTYP  LCINT</span></div><div><span style="font-size: 12pt;">      1005         0       1.0       1.0       0.0       0.0         0         0</span></div><div><span style="font-size: 12pt;">$ time    timestep</span></div><div><span style="font-size: 12pt;">                 0.0                1e-2</span></div><div><span style="font-size: 12pt;">                 10.0               1e-2</span></div><div><span style="font-size: 12pt;">*CONTROL_TIMESTEP</span></div><div><span style="font-size: 12pt;">$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st</span></div><div><span style="font-size: 12pt;">      1e-2       1.0         0       0.0         0      1005         0         0</span></div><div><span style="font-size: 12pt;">*CONTROL_ENERGY</span></div><div><span style="font-size: 12pt;">$     hgen      rwen    slnten     rylen</span></div><div><span style="font-size: 12pt;">         2         2</span></div><div><span style="font-size: 12pt;">*CONTROL_OUTPUT</span></div><div><span style="font-size: 12pt;">$    npopt    neecho    nrefup    iaccop     opifs    ipnint    ikedit</span></div><div><span style="font-size: 12pt;">         1         3</span></div><div><span style="font-size: 12pt;">*CONTROL_SHELL</span></div><div><span style="font-size: 12pt;">$   wrpang    itrist     irnxx    istupd    theory       bwc     miter</span></div><div><span style="font-size: 12pt;">                                       1         2</span></div><div><span style="font-size: 12pt;">*DATABASE_BINARY_D3PLOT</span></div><div><span style="font-size: 12pt;">$       dt      lcdt</span></div><div><span style="font-size: 12pt;">      1.00</span></div><div><span style="font-size: 12pt;">*DATABASE_EXTENT_BINARY</span></div><div><span style="font-size: 12pt;">$    neiph     neips    maxint    strflg    sigflg    epsflg    rltflg    engflg</span></div><br><div><span style="font-size: 12pt;">$   cmpflg    ieverp    beamip</span></div><div><span style="font-size: 12pt;">                   1</span></div><div><span style="font-size: 12pt;">*DATABASE_BINARY_D3THDT</span></div><div><span style="font-size: 12pt;">$       dt      lcdt</span></div><div><span style="font-size: 12pt;">    999999</span></div><div><span style="font-size: 12pt;">*DATABASE_GLSTAT</span></div><div><span style="font-size: 12pt;">$       dt</span></div><div><span style="font-size: 12pt;">      0.10</span></div><div><span style="font-size: 12pt;">*DATABASE_MATSUM</span></div><div><span style="font-size: 12pt;">$       dt</span></div><div><span style="font-size: 12pt;">      0.10</span></div><div><span style="font-size: 12pt;">*DATABASE_NODOUT</span></div><div><span style="font-size: 12pt;">$       dt</span></div><div><span style="font-size: 12pt;">      0.10</span></div><div><span style="font-size: 12pt;">*DATABASE_RBDOUT</span></div><div><span style="font-size: 12pt;">$       dt</span></div><div><span style="font-size: 12pt;">      0.10</span></div><div><span style="font-size: 12pt;">*PART</span></div><div><span style="font-size: 12pt;">$      pid       sid       mid     eosid      hgid    adpopt</span></div><div><span style="font-size: 12pt;">part1</span></div><div><span style="font-size: 12pt;">         1         2         2</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">$  Aluminum</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">*MAT_PLASTIC_KINEMATIC</span></div><div><span style="font-size: 12pt;">$      mid        ro         e        pr      sigy      etan      beta</span></div><div><span style="font-size: 12pt;">         1   2.70e-6      68.9     0.330     0.286   0.00689</span></div><div><span style="font-size: 12pt;">$      src       srp        fs</span></div><br><div><span style="font-size: 12pt;">*SECTION_SHELL</span></div><div><span style="font-size: 12pt;">$</span></div><div><span style="font-size: 12pt;">$      sid    elform      shrf       nip     propt   qr/irid     icomp</span></div><div><span style="font-size: 12pt;">         1         2</span></div><div><span style="font-size: 12pt;">$       t1        t2        t3        t4      nloc</span></div><div><span style="font-size: 12pt;">       1.0       1.0       1.0       1.0</span></div><div><span style="font-size: 12pt;">*SECTION_SHELL</span></div><div><span style="font-size: 12pt;">$       ID    elform      shrf       nip     propt        qr     icomp     setyp</span></div><div><span style="font-size: 12pt;">         2         2  0.833333         3         0         0         0         0</span></div><div><span style="font-size: 12pt;">$       t1        t2        t3        t4      nloc     marea              unused</span></div><div><span style="font-size: 12pt;">       0.1       0.1       0.1       0.1         0         0                    </span></div><div><span style="font-size: 12pt;">*MAT_HYPERELASTIC_RUBBER</span></div><div><span style="font-size: 12pt;">$       ID        ro        pr         N        Nv         G      sigf       Ref</span></div><div><span style="font-size: 12pt;">         2 8.699E-07 0.4999172         0         0         0         0         0</span></div><div><span style="font-size: 12pt;">$      C10       C01       C11       C20       C02       C30             unused1</span></div><div><span style="font-size: 12pt;">   0.11e-3  0.028e-3         0         0         0         0                    </span></div><div><span style="font-size: 12pt;">*NODE</span></div><div><span style="font-size: 12pt;">$   node               x               y               z      tc      rc</span></div><div><span style="font-size: 12pt;">       1             0.0             0.0             0.0</span></div><div><span style="font-size: 12pt;">       2             1.0             0.0             0.0</span></div><div><span style="font-size: 12pt;">       3             1.0             1.0             0.0</span></div><div><span style="font-size: 12pt;">       4             0.0             1.0             0.0</span></div><div><span style="font-size: 12pt;">$$$$  Define Elements</span></div><div><span style="font-size: 12pt;">*ELEMENT_SHELL</span></div><div><span style="font-size: 12pt;">$    eid     pid      n1      n2      n3      n4</span></div><div><span style="font-size: 12pt;">       1       1       1       2       3       4</span></div><div><span style="font-size: 12pt;">*DEFORMABLE_TO_RIGID</span></div><div><span style="font-size: 12pt;">         1</span></div><div><span style="font-size: 12pt;">*DEFORMABLE_TO_RIGID_AUTOMATIC</span></div><div><span style="font-size: 12pt;">$#   SWSET      CODE     TIME1     TIME2     TIME3     ENTNO     RELSW    PAIRED</span></div><div><span style="font-size: 12pt;">         1         0      10.0      1.0e20     0.0         0         0         0</span></div><div><span style="font-size: 12pt;">$#    NRBF      NCSF       RWF     DTMAX       D2R       R2D    OFFSET</span></div><div><span style="font-size: 12pt;">         0         0         0       0.0         0         1       0.0</span></div><div><span style="font-size: 12pt;">$# Rigid to deformable card: one card because R2D = 1</span></div><div><span style="font-size: 12pt;">$#     PID     PTYPE</span></div><div><span style="font-size: 12pt;">         1      PART</span></div><div><span style="font-size: 12pt;">*END</span></div></div>
Avatar
maoy@umsystem.edu May 29, 2026
How to model a crack using MAT_ADD_EROSION in Ansys LS-Dyna?
I wish to use the MAT_ADD_EROSION model in LS-DYNA to capture material failure and erosion effectively. Are there any resources to help me get started? Thanks!    
Avatar
innovation customer May 29, 2026
Simulating a prestrained material
Is it possible to simulate a pre-strained material directly without any pre-simulation? We have conducted uniaxial prestraining on a material and would now like to run another simulation. The major and minor strains from our previous operation are available to us. Thanks!
Avatar
innovation customer May 29, 2026
Question Regarding Applicability of YLD2004-18p for Shell Elements
<p> Dear LS-DYNA Support Team, I hope you are doing well.</p><p>I would like to ask whether the yield locus model **YLD2004-18p** can be applied to **shell elements** in LS-DYNA.</p><p>Could you please let me know if this material model is supported for shell formulations, and if there are any specific limitations or requirements that need to be considered?</p><p>Thank you very much for your support.</p><p> Best regards,</p><p> Philipp</p>
Avatar
pneufeld@hs-koblenz.de May 29, 2026

Still stuck or need help?

Submit a Support Request