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

Download the presentations of the German LS-DYNA Forum 2024

More Details

Community

Explicit Thermal Solver Crash with S-ALE Mesh – SETID & Phase Change
<div class="text-neutral-content"><div class="mb-sm mb-xs px-md xs:px-0 overflow-hidden" data-post-click-location="text-body"><div id="t3_1rr5r9z-post-rtjson-content" class="md text-14-scalable pb-2xs [--emote-size:20px]" dir="auto"><p>Hi,</p><p>I'm modeling hypervelocity impact (2.22 km/s tungsten carbide sphere on HSLA-100 steel) with coupled thermal-structural analysis using S-ALE (Structured ALE) mesh.</p><p>Initial Attempt & Issue:</p><p>I first tried the implicit thermal solver using:</p><p>- *CONTROL_SOLUTION (Thermal-Structural Analysis)</p><p>- *CONTROL_THERMAL_NONLINEAR</p><p>- *CONTROL_THERMAL_SOLVER</p><p>- *CONTROL_THERMAL_TIMESTEP</p><p>However, I observed physically unrealistic temperature rates (dT/dt on the order of billions of °K), which is physically impossible. I switched to the Explicit Thermal Solver (*CONTROL_EXPLICIT_THERMAL_* cards) as an alternative, but encountered a fatal initialization error.</p><p>Current Error:</p><p></p><p>forrtl: severe (164): Program Exception - integer divide by zero</p><p>Image PC Routine Line Source</p><p>lsdyna_mpp_dp_imp 00007FF731E6B4CC XPLCTH_INIT 1148 dyn20x.F</p><p>...</p><p></p><p>The crash occurs during initialization (XPLCTH_INIT), suggesting a division by zero in the thermal data setup.</p><p>Specific Questions:</p><p>Explicit Thermal Setup with S-ALE:</p><p>For cards like *DATABASE_ALE, *CONTROL_EXPLICIT_THERMAL_INITIAL, and *CONTROL_EXPLICIT_THERMAL_OUTPUT, the manual states SETID can be left blank for "all parts," but LS-PrePost/LS-DYNA requires an integer. Since S-ALE mesh elements are generated internally at runtime (not predefined in *ELEMENT_SHELL), I cannot reference them via traditional *SET_SHELL_LIST.</p><p>- How do I properly define SETID for S-ALE domains in these cards? Should I use *SET_MULTI_MATERIAL_GROUP instead?</p><p>Phase Change Modeling:</p><p>My material definition includes *MAT_THERMAL_ISOTROPIC_PHASE_CHANGE for the steel target (solid-to-liquid). Does the Explicit Thermal Solver support phase change energy (latent heat), or is this capability restricted to the implicit thermal solver? If supported, are there specific considerations for S-ALE multi-material groups?</p><p>Divide-by-Zero Cause:</p><p>The error trace points to XPLCTH_INIT (dyn20x.F:1148). This typically indicates zero thermal mass (density × specific heat) in a thermal part definition. Given that my vacuum/void material (*MAT_ALE_VACUUM) has near-zero thermal properties, could this be triggering the crash? Should vacuum regions be excluded from the thermal solver via *SET_PART_LIST, or is there a minimum non-zero thermal density requirement?</p><p>Model Context:</p><p>- 2D Axisymmetric S-ALE mesh (fine resolution)</p><p>- 3 ALE Multi-Material Groups: HSLA-100 Steel, Tungsten Carbide, Vacuum</p><p>- Thermal properties defined via *MAT_THERMAL_ISOTROPIC and *MAT_THERMAL_ISOTROPIC_PHASE_CHANGE</p><p>- Using MPP R14.1.1 on Windows</p><p>Keyword file attached for reference:</p><p><a class="relative pointer-events-auto a cursor-pointer underline" href="https://buffalo.box.com/s/y7rz3bqg4xsvuir6z94b50mcg2g3jp1s" target="_blank" rel="noopener">https://buffalo.box.com/s/y7rz3bqg4xsvuir6z94b50mcg2g3jp1s</p><p>Any guidance on proper explicit thermal initialization with S-ALE and phase change setup would be greatly appreciated.</p><p>Thanks,</p><p>Amrith</p></div></div></div><div> </div>
Avatar
Amrith Mariappan March 11, 2026
Difference between Autodyn standalone and Ansys/explicit dynamics plus Autodyn
<p>Hi all, I'm trying to learn simulating hypervelocity impact by myself, I notice that most of the journal papers choose Autodyn. And my question is what's diference between using Autodyn standalone and using Ansys/ explicit dynamics plus Autosdyn? Also, which way is easier to model layered composite material? </p>
Avatar
rileylau0913@gmail.com March 11, 2026
LS-DYNA error code 10125 – advice needed
<p></p><p>I got these error codes while trying to run an FEA simulation, and I have tried to go back in and define the section IDs, but it isn't working. Does anyone have advice to fix these?</p>
Avatar
willisbara0921@gmail.com March 11, 2026
Models repel each other
<p>Hi All, </p><p>I created 2 beams. After I applied cross section, there was an overlap between the two beams. I applied *contact_automatic_general between 2 beams. I found the beam would repel each other because there is model overlap. Does anybody know how to get rid of the repulsive force, even though there is overlap?</p><p></p><p></p><p></p><p></p>
Avatar
tut06916@temple.edu March 10, 2026
How to model laminated CFRP in Autodyn
<p>Hi all,</p><p>I am trying to simulate a hypervelocity impact (~5 km/s) on a CFRP laminate consisting of 16 plies with 15 resin-rich interlayers using the SPH method.</p><p>I already know the stacking sequence and fiber orientation of each ply. However, I am unsure how to properly construct the laminated structure in an SPH model.</p><p>In particular, I would like to understand:</p><p>• How to assign the orthotropic fiber orientation for each ply when using SPH particles.<br>• How to represent the interlaminar interfaces between plies (e.g., contact or cohesive behavior).<br>• How to build the layered laminate structure itself so that ply-level damage mechanisms such as fiber breakage, fiber–matrix debonding, and ply delamination can be captured.</p><p>What is the recommended approach to model the ply orientations, interply interfaces, and the laminate stacking structure when using SPH for CFRP impact simulations?</p>
Avatar
rileylau0913@gmail.com March 9, 2026
Stress Values Looks Weird on Johnson Cook Material Model
<p>Hello everyone</p><p>I'm working on a ballistic analysis on Ls-Dyna. Before I do my own simulations I'm trying to replicate an article right now. In that article a rigid ball smashes into a aluminum armor plate. Material model for aluminum is Johnson Cook and the values of it were given in article but when I run the simulation stress values looks weird. For instance X,Y and Z stress values all looks the same, there is NO shear stress whatsoever, there are only small dots of stress in von mises and tresca stress, effective plastic strain is non-existant and the bulge that appears when the ball penetrates through the armor seems too big. What could be causes of this? When I use simplified johnson cook model stress and strain values looks fine but when I switch to normal johnson cook everything looks wrong. Can somebody help please??? Pictures of my problem are in the link: https://imgur.com/a/nTRUaSL</p>
Avatar
yvzslm48@gmail.com March 8, 2026
RC blast door SALE simulation: deflection trend matches literature but peak disp
<p data-start="216" data-end="319">When simulating an RC blast door using <strong data-start="255" data-end="273">SALE (ALE–FSI) in <span class="whitespace-normal">LS-DYNA</span>,</p><p data-start="321" data-end="450">the <strong data-start="325" data-end="393">time-history trend of the door deflection matches the literature, but the <strong data-start="403" data-end="449">peak displacement is significantly smaller.</p><div class="TyagGW_tableContainer"><div class="group TyagGW_tableWrapper flex flex-col-reverse w-fit" tabindex="-1"><table class="w-fit min-w-(--thread-content-width)" data-start="452" data-end="532"><thead data-start="452" data-end="473"><tr data-start="452" data-end="473"><th class="" data-start="452" data-end="454" data-col-size="sm"> </th><th class="" data-start="454" data-end="473" data-col-size="sm">Peak deflection</th></tr></thead><tbody data-start="484" data-end="532"><tr data-start="484" data-end="508"><td data-start="484" data-end="497" data-col-size="sm">Literature</td><td data-start="497" data-end="508" data-col-size="sm">3.47 mm</td></tr><tr data-start="509" data-end="532"><td data-start="509" data-end="520" data-col-size="sm">My model</td><td data-start="520" data-end="532" data-col-size="sm">< 2.0 mm</td></tr></tbody></table></div></div><p data-start="534" data-end="608">So the <strong data-start="541" data-end="607">response shape is correct, but the amplitude is underestimated.</p><hr data-start="610" data-end="613"><h1 data-section-id="l0dx7p" data-start="615" data-end="630">Model summary</h1><h2 data-section-id="1fj62nz" data-start="632" data-end="643">Geometry</h2><ul data-start="644" data-end="794"><li data-section-id="19hxg5u" data-start="644" data-end="675"><p data-start="646" data-end="675">RC door thickness: <strong data-start="665" data-end="675">100 mm</p><li data-section-id="1g5rvbz" data-start="676" data-end="708"><p data-start="678" data-end="708">Door location: <strong data-start="693" data-end="708">z = 0–0.1 m</p><li data-section-id="d715vz" data-start="709" data-end="741"><p data-start="711" data-end="741">ALE domain: <strong data-start="723" data-end="741">z = -0.6–3.4 m</p><li data-section-id="28srp" data-start="742" data-end="794"><p data-start="744" data-end="794">Door <strong data-start="749" data-end="794">rear surface initially touching the frame</p><hr data-start="796" data-end="799"><h2 data-section-id="13lbijm" data-start="801" data-end="813">Explosive</h2><ul data-start="814" data-end="878"><li data-section-id="olji2y" data-start="814" data-end="832"><p data-start="816" data-end="832">TNT: <strong data-start="821" data-end="832">0.15 kg</p><li data-section-id="1s39v9c" data-start="833" data-end="878"><p data-start="835" data-end="878">Half model with <strong data-start="851" data-end="878">symmetry plane at x = 0</p><hr data-start="880" data-end="883"><h2 data-section-id="1g556vy" data-start="885" data-end="896">ALE mesh</h2><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">10 mm</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><hr data-start="916" data-end="919"><h2 data-section-id="1bvqh7m" data-start="921" data-end="932">Concrete</h2><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">MAT_072R3 (KCC)<br>f'c ≈ 19.7 MPa</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><hr data-start="977" data-end="980"><h2 data-section-id="ayu17e" data-start="982" data-end="998">Reinforcement</h2><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class="pointer-events-none absolute inset-x-4 top-12 bottom-4"><div class="pointer-events-none sticky z-40 shrink-0 z-1!"><div class="sticky bg-token-border-light"> </div></div></div><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">MAT_PLASTIC_KINEMATIC<br>E = 200 GPa<br>ETAN = 1e9 Pa</div></div></div></div></div></div></div></div></div></div></div></div><p data-start="1060" data-end="1075">Embedded using: CONSTRAINED_BEAM_IN_SOLID_PENALTY</p><hr data-start="1124" data-end="1127"><h2 data-section-id="qfdzp" data-start="1129" data-end="1146">Solid elements</h2><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">ELFORM = 1</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><hr data-start="1171" data-end="1174"><h2 data-section-id="142oigs" data-start="1176" data-end="1191">FSI coupling</h2><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class="pointer-events-none absolute inset-x-4 top-12 bottom-4"><div class="pointer-events-none sticky z-40 shrink-0 z-1!"><div class="sticky bg-token-border-light"> </div></div></div><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">ALE_STRUCTURED_FSI<br>LSTRSID = 1<br>(part set includes parts 1,2,3,4,5,11)<br>PFAC = 1<br>FLIP = 0</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><hr data-start="1293" data-end="1296"><h2 data-section-id="a2jb2d" data-start="1298" data-end="1308">Contact</h2><p data-start="1309" data-end="1332">Door–frame interaction:</p><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class="pointer-events-none absolute inset-x-4 top-12 bottom-4"><div class="pointer-events-none sticky z-40 shrink-0 z-1!"><div class="sticky bg-token-border-light"> </div></div></div><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">CONTACT_AUTOMATIC_SURFACE_TO_SURFACE</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><p data-start="1384" data-end="1466">Frame is <strong data-start="1393" data-end="1465">located behind the door and initially touching the door rear surface.</p><p data-start="1468" data-end="1477">Friction:</p><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">FS = 0.2<br>FD = 0.15</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><hr data-start="1512" data-end="1515"><h1 data-section-id="f1jprp" data-start="1517" data-end="1542">Tests already performed</h1><h3 data-section-id="1s8do4e" data-start="1544" data-end="1567">ALE backside domain</h3><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class="pointer-events-none absolute inset-x-4 top-12 bottom-4"><div class="pointer-events-none sticky z-40 shrink-0 z-1!"><div class="sticky bg-token-border-light"> </div></div></div><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">z = -0.6 → -0.3</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><p data-start="1596" data-end="1617">No noticeable effect.</p><hr data-start="1619" data-end="1622"><h3 data-section-id="1ixko64" data-start="1624" data-end="1658">Beam–solid coupling parameters</h3><p data-start="1659" data-end="1713">Adjusted penalty parameters → <strong data-start="1689" data-end="1713">no noticeable effect</p><hr data-start="1715" data-end="1718"><h3 data-section-id="19ok7fa" data-start="1720" data-end="1745">Frame-wall connection</h3><p data-start="1746" data-end="1801">Removing shared nodes → <strong data-start="1770" data-end="1801">response trend became worse</p><hr data-start="1803" data-end="1806"><h1 data-section-id="hh93l4" data-start="1808" data-end="1818">Question</h1><p data-start="1820" data-end="1859">In <strong data-start="1823" data-end="1852">ALE–FSI blast simulations, when:</p><ul data-start="1861" data-end="1993"><li data-section-id="laum6d" data-start="1861" data-end="1898"><p data-start="1863" data-end="1898">pressure history appears reasonable</p><li data-section-id="147b0ew" data-start="1899" data-end="1951"><p data-start="1901" data-end="1951">displacement time-history trend matches literature</p><li data-section-id="12d7n9f" data-start="1952" data-end="1993"><p data-start="1954" data-end="1993">but peak displacement is underestimated</p><p data-start="1995" data-end="2031">what are the <strong data-start="2008" data-end="2030">most common causes?</p><p data-start="2033" data-end="2068">In particular, I would like to ask:</p><p data-start="2072" data-end="2231">Could the <strong data-start="2082" data-end="2156">door–frame initial contact (rear surface initially touching the frame) act as an early support condition and reduce the door bending deformation?</p><p data-start="2233" data-end="2301">Any suggestions or similar experiences would be greatly appreciated.</p>
Avatar
n66134267@gs.ncku.edu.tw March 7, 2026
YOSHIDA UEMORI DEEP DRAWING MODELS
<p>Hello,</p><p>I would like to simulate the process of sheet metal stamping and then springback. In both processes, I would like to take into account the Bausinger effect, anisotropy , softening and workhardening stagnation stress reversals, and strain rate . The Yoshida Uemori model + Anisotropic model, e.g., Yld2000 + strain rate, would be suitable for this task. The problem is that LS Dyna does not have such a model . Therefore, I would like to ask if in the first simulation, i.e., the sheet metal forming simulation, I could use the mat_133 model, then transfer the simulation results to a new springback simulation and change the material model to mat_242, and then run a new springback simulation with the mat_242 material model? The mat_133 model has a backstress historyhttps://lsdyna.ansys.com/history-variables/</p><p> </p><p> </p>
Avatar
euklidespl@gmail.com March 6, 2026
Tied contact failed(DYNAMIC_RELAXATION in the total retart ) Need Advice
<p class="ds-markdown-paragraph">Hello everyone,</p><p class="ds-markdown-paragraph">I am encountering an issue with tied contacts during a total restart analysis for springback calculation and would appreciate some advice. (1)First analysis  A steel formwork falls  onto the rock,used shell element: and used TIED_SHELL_EDGE_TO_SURFACE_BEAM_OFFSET and TIED_SHELL_EDGE_TO_SURFACE_BEAM_OFFSET:<br> calculation in this step is normal,and tied contact works fine: (2)total start used CONTROL_DYNAMIC_RELAXATION keyword to calculate the springback but all tied contact was failed:<br> we could see that penetration occurred between the parts.  Has anyone encountered a similar issue? Any suggestions on what parameters I should check or adjust would be greatly appreciated.</p>
Avatar
383291376@qq.com March 6, 2026
Getting started with usermat
<p>Hello,</p><p>I requested the usermat version of LS-DYNA from my institution's software office, from the link found in here: https://ftp.lstc.com/anonymous/outgoing/support/FAQ/user_defined_materials.faq. However, after receiving the files, the file types were unexpected, and different from what the previous posts and user manual has mentioned. </p><p>They provided me with a folder named "ls-dyna....mpi" and inside were only a license folder, and 3 executable files named l2a, lstc_client, and ls-dyna...mpi. Are these the files I should be expecting, and if so, how do I use them?</p><p>Thanks </p>
Avatar
kshi84@gatech.edu March 5, 2026
Explicit Thermal Solver Crash with S-ALE Mesh – SETID & Phase Change
<div class="text-neutral-content"><div class="mb-sm mb-xs px-md xs:px-0 overflow-hidden" data-post-click-location="text-body"><div id="t3_1rr5r9z-post-rtjson-content" class="md text-14-scalable pb-2xs [--emote-size:20px]" dir="auto"><p>Hi,</p><p>I'm modeling hypervelocity impact (2.22 km/s tungsten carbide sphere on HSLA-100 steel) with coupled thermal-structural analysis using S-ALE (Structured ALE) mesh.</p><p>Initial Attempt & Issue:</p><p>I first tried the implicit thermal solver using:</p><p>- *CONTROL_SOLUTION (Thermal-Structural Analysis)</p><p>- *CONTROL_THERMAL_NONLINEAR</p><p>- *CONTROL_THERMAL_SOLVER</p><p>- *CONTROL_THERMAL_TIMESTEP</p><p>However, I observed physically unrealistic temperature rates (dT/dt on the order of billions of °K), which is physically impossible. I switched to the Explicit Thermal Solver (*CONTROL_EXPLICIT_THERMAL_* cards) as an alternative, but encountered a fatal initialization error.</p><p>Current Error:</p><p></p><p>forrtl: severe (164): Program Exception - integer divide by zero</p><p>Image PC Routine Line Source</p><p>lsdyna_mpp_dp_imp 00007FF731E6B4CC XPLCTH_INIT 1148 dyn20x.F</p><p>...</p><p></p><p>The crash occurs during initialization (XPLCTH_INIT), suggesting a division by zero in the thermal data setup.</p><p>Specific Questions:</p><p>Explicit Thermal Setup with S-ALE:</p><p>For cards like *DATABASE_ALE, *CONTROL_EXPLICIT_THERMAL_INITIAL, and *CONTROL_EXPLICIT_THERMAL_OUTPUT, the manual states SETID can be left blank for "all parts," but LS-PrePost/LS-DYNA requires an integer. Since S-ALE mesh elements are generated internally at runtime (not predefined in *ELEMENT_SHELL), I cannot reference them via traditional *SET_SHELL_LIST.</p><p>- How do I properly define SETID for S-ALE domains in these cards? Should I use *SET_MULTI_MATERIAL_GROUP instead?</p><p>Phase Change Modeling:</p><p>My material definition includes *MAT_THERMAL_ISOTROPIC_PHASE_CHANGE for the steel target (solid-to-liquid). Does the Explicit Thermal Solver support phase change energy (latent heat), or is this capability restricted to the implicit thermal solver? If supported, are there specific considerations for S-ALE multi-material groups?</p><p>Divide-by-Zero Cause:</p><p>The error trace points to XPLCTH_INIT (dyn20x.F:1148). This typically indicates zero thermal mass (density × specific heat) in a thermal part definition. Given that my vacuum/void material (*MAT_ALE_VACUUM) has near-zero thermal properties, could this be triggering the crash? Should vacuum regions be excluded from the thermal solver via *SET_PART_LIST, or is there a minimum non-zero thermal density requirement?</p><p>Model Context:</p><p>- 2D Axisymmetric S-ALE mesh (fine resolution)</p><p>- 3 ALE Multi-Material Groups: HSLA-100 Steel, Tungsten Carbide, Vacuum</p><p>- Thermal properties defined via *MAT_THERMAL_ISOTROPIC and *MAT_THERMAL_ISOTROPIC_PHASE_CHANGE</p><p>- Using MPP R14.1.1 on Windows</p><p>Keyword file attached for reference:</p><p><a class="relative pointer-events-auto a cursor-pointer underline" href="https://buffalo.box.com/s/y7rz3bqg4xsvuir6z94b50mcg2g3jp1s" target="_blank" rel="noopener">https://buffalo.box.com/s/y7rz3bqg4xsvuir6z94b50mcg2g3jp1s</p><p>Any guidance on proper explicit thermal initialization with S-ALE and phase change setup would be greatly appreciated.</p><p>Thanks,</p><p>Amrith</p></div></div></div><div> </div>
Avatar
Amrith Mariappan March 11, 2026
Difference between Autodyn standalone and Ansys/explicit dynamics plus Autodyn
<p>Hi all, I'm trying to learn simulating hypervelocity impact by myself, I notice that most of the journal papers choose Autodyn. And my question is what's diference between using Autodyn standalone and using Ansys/ explicit dynamics plus Autosdyn? Also, which way is easier to model layered composite material? </p>
Avatar
rileylau0913@gmail.com March 11, 2026
LS-DYNA error code 10125 – advice needed
<p></p><p>I got these error codes while trying to run an FEA simulation, and I have tried to go back in and define the section IDs, but it isn't working. Does anyone have advice to fix these?</p>
Avatar
willisbara0921@gmail.com March 11, 2026
Models repel each other
<p>Hi All, </p><p>I created 2 beams. After I applied cross section, there was an overlap between the two beams. I applied *contact_automatic_general between 2 beams. I found the beam would repel each other because there is model overlap. Does anybody know how to get rid of the repulsive force, even though there is overlap?</p><p></p><p></p><p></p><p></p>
Avatar
tut06916@temple.edu March 10, 2026
How to model laminated CFRP in Autodyn
<p>Hi all,</p><p>I am trying to simulate a hypervelocity impact (~5 km/s) on a CFRP laminate consisting of 16 plies with 15 resin-rich interlayers using the SPH method.</p><p>I already know the stacking sequence and fiber orientation of each ply. However, I am unsure how to properly construct the laminated structure in an SPH model.</p><p>In particular, I would like to understand:</p><p>• How to assign the orthotropic fiber orientation for each ply when using SPH particles.<br>• How to represent the interlaminar interfaces between plies (e.g., contact or cohesive behavior).<br>• How to build the layered laminate structure itself so that ply-level damage mechanisms such as fiber breakage, fiber–matrix debonding, and ply delamination can be captured.</p><p>What is the recommended approach to model the ply orientations, interply interfaces, and the laminate stacking structure when using SPH for CFRP impact simulations?</p>
Avatar
rileylau0913@gmail.com March 9, 2026
Stress Values Looks Weird on Johnson Cook Material Model
<p>Hello everyone</p><p>I'm working on a ballistic analysis on Ls-Dyna. Before I do my own simulations I'm trying to replicate an article right now. In that article a rigid ball smashes into a aluminum armor plate. Material model for aluminum is Johnson Cook and the values of it were given in article but when I run the simulation stress values looks weird. For instance X,Y and Z stress values all looks the same, there is NO shear stress whatsoever, there are only small dots of stress in von mises and tresca stress, effective plastic strain is non-existant and the bulge that appears when the ball penetrates through the armor seems too big. What could be causes of this? When I use simplified johnson cook model stress and strain values looks fine but when I switch to normal johnson cook everything looks wrong. Can somebody help please??? Pictures of my problem are in the link: https://imgur.com/a/nTRUaSL</p>
Avatar
yvzslm48@gmail.com March 8, 2026
RC blast door SALE simulation: deflection trend matches literature but peak disp
<p data-start="216" data-end="319">When simulating an RC blast door using <strong data-start="255" data-end="273">SALE (ALE–FSI) in <span class="whitespace-normal">LS-DYNA</span>,</p><p data-start="321" data-end="450">the <strong data-start="325" data-end="393">time-history trend of the door deflection matches the literature, but the <strong data-start="403" data-end="449">peak displacement is significantly smaller.</p><div class="TyagGW_tableContainer"><div class="group TyagGW_tableWrapper flex flex-col-reverse w-fit" tabindex="-1"><table class="w-fit min-w-(--thread-content-width)" data-start="452" data-end="532"><thead data-start="452" data-end="473"><tr data-start="452" data-end="473"><th class="" data-start="452" data-end="454" data-col-size="sm"> </th><th class="" data-start="454" data-end="473" data-col-size="sm">Peak deflection</th></tr></thead><tbody data-start="484" data-end="532"><tr data-start="484" data-end="508"><td data-start="484" data-end="497" data-col-size="sm">Literature</td><td data-start="497" data-end="508" data-col-size="sm">3.47 mm</td></tr><tr data-start="509" data-end="532"><td data-start="509" data-end="520" data-col-size="sm">My model</td><td data-start="520" data-end="532" data-col-size="sm">< 2.0 mm</td></tr></tbody></table></div></div><p data-start="534" data-end="608">So the <strong data-start="541" data-end="607">response shape is correct, but the amplitude is underestimated.</p><hr data-start="610" data-end="613"><h1 data-section-id="l0dx7p" data-start="615" data-end="630">Model summary</h1><h2 data-section-id="1fj62nz" data-start="632" data-end="643">Geometry</h2><ul data-start="644" data-end="794"><li data-section-id="19hxg5u" data-start="644" data-end="675"><p data-start="646" data-end="675">RC door thickness: <strong data-start="665" data-end="675">100 mm</p><li data-section-id="1g5rvbz" data-start="676" data-end="708"><p data-start="678" data-end="708">Door location: <strong data-start="693" data-end="708">z = 0–0.1 m</p><li data-section-id="d715vz" data-start="709" data-end="741"><p data-start="711" data-end="741">ALE domain: <strong data-start="723" data-end="741">z = -0.6–3.4 m</p><li data-section-id="28srp" data-start="742" data-end="794"><p data-start="744" data-end="794">Door <strong data-start="749" data-end="794">rear surface initially touching the frame</p><hr data-start="796" data-end="799"><h2 data-section-id="13lbijm" data-start="801" data-end="813">Explosive</h2><ul data-start="814" data-end="878"><li data-section-id="olji2y" data-start="814" data-end="832"><p data-start="816" data-end="832">TNT: <strong data-start="821" data-end="832">0.15 kg</p><li data-section-id="1s39v9c" data-start="833" data-end="878"><p data-start="835" data-end="878">Half model with <strong data-start="851" data-end="878">symmetry plane at x = 0</p><hr data-start="880" data-end="883"><h2 data-section-id="1g556vy" data-start="885" data-end="896">ALE mesh</h2><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">10 mm</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><hr data-start="916" data-end="919"><h2 data-section-id="1bvqh7m" data-start="921" data-end="932">Concrete</h2><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">MAT_072R3 (KCC)<br>f'c ≈ 19.7 MPa</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><hr data-start="977" data-end="980"><h2 data-section-id="ayu17e" data-start="982" data-end="998">Reinforcement</h2><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class="pointer-events-none absolute inset-x-4 top-12 bottom-4"><div class="pointer-events-none sticky z-40 shrink-0 z-1!"><div class="sticky bg-token-border-light"> </div></div></div><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">MAT_PLASTIC_KINEMATIC<br>E = 200 GPa<br>ETAN = 1e9 Pa</div></div></div></div></div></div></div></div></div></div></div></div><p data-start="1060" data-end="1075">Embedded using: CONSTRAINED_BEAM_IN_SOLID_PENALTY</p><hr data-start="1124" data-end="1127"><h2 data-section-id="qfdzp" data-start="1129" data-end="1146">Solid elements</h2><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">ELFORM = 1</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><hr data-start="1171" data-end="1174"><h2 data-section-id="142oigs" data-start="1176" data-end="1191">FSI coupling</h2><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class="pointer-events-none absolute inset-x-4 top-12 bottom-4"><div class="pointer-events-none sticky z-40 shrink-0 z-1!"><div class="sticky bg-token-border-light"> </div></div></div><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">ALE_STRUCTURED_FSI<br>LSTRSID = 1<br>(part set includes parts 1,2,3,4,5,11)<br>PFAC = 1<br>FLIP = 0</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><hr data-start="1293" data-end="1296"><h2 data-section-id="a2jb2d" data-start="1298" data-end="1308">Contact</h2><p data-start="1309" data-end="1332">Door–frame interaction:</p><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class="pointer-events-none absolute inset-x-4 top-12 bottom-4"><div class="pointer-events-none sticky z-40 shrink-0 z-1!"><div class="sticky bg-token-border-light"> </div></div></div><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">CONTACT_AUTOMATIC_SURFACE_TO_SURFACE</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><p data-start="1384" data-end="1466">Frame is <strong data-start="1393" data-end="1465">located behind the door and initially touching the door rear surface.</p><p data-start="1468" data-end="1477">Friction:</p><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">FS = 0.2<br>FD = 0.15</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><hr data-start="1512" data-end="1515"><h1 data-section-id="f1jprp" data-start="1517" data-end="1542">Tests already performed</h1><h3 data-section-id="1s8do4e" data-start="1544" data-end="1567">ALE backside domain</h3><div class="relative w-full mt-4 mb-1"><div class=""><div class="relative"><div class="h-full min-h-0 min-w-0"><div class="h-full min-h-0 min-w-0"><div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"><div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"><div class="pointer-events-none absolute inset-x-4 top-12 bottom-4"><div class="pointer-events-none sticky z-40 shrink-0 z-1!"><div class="sticky bg-token-border-light"> </div></div></div><div class=""><div class="relative z-0 flex max-w-full"><div id="code-block-viewer" class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy" dir="ltr"><div class="cm-scroller"><div class="cm-content q9tKkq_readonly">z = -0.6 → -0.3</div></div></div></div></div></div></div></div></div><div class=""><div class=""> </div></div></div></div></div><p data-start="1596" data-end="1617">No noticeable effect.</p><hr data-start="1619" data-end="1622"><h3 data-section-id="1ixko64" data-start="1624" data-end="1658">Beam–solid coupling parameters</h3><p data-start="1659" data-end="1713">Adjusted penalty parameters → <strong data-start="1689" data-end="1713">no noticeable effect</p><hr data-start="1715" data-end="1718"><h3 data-section-id="19ok7fa" data-start="1720" data-end="1745">Frame-wall connection</h3><p data-start="1746" data-end="1801">Removing shared nodes → <strong data-start="1770" data-end="1801">response trend became worse</p><hr data-start="1803" data-end="1806"><h1 data-section-id="hh93l4" data-start="1808" data-end="1818">Question</h1><p data-start="1820" data-end="1859">In <strong data-start="1823" data-end="1852">ALE–FSI blast simulations, when:</p><ul data-start="1861" data-end="1993"><li data-section-id="laum6d" data-start="1861" data-end="1898"><p data-start="1863" data-end="1898">pressure history appears reasonable</p><li data-section-id="147b0ew" data-start="1899" data-end="1951"><p data-start="1901" data-end="1951">displacement time-history trend matches literature</p><li data-section-id="12d7n9f" data-start="1952" data-end="1993"><p data-start="1954" data-end="1993">but peak displacement is underestimated</p><p data-start="1995" data-end="2031">what are the <strong data-start="2008" data-end="2030">most common causes?</p><p data-start="2033" data-end="2068">In particular, I would like to ask:</p><p data-start="2072" data-end="2231">Could the <strong data-start="2082" data-end="2156">door–frame initial contact (rear surface initially touching the frame) act as an early support condition and reduce the door bending deformation?</p><p data-start="2233" data-end="2301">Any suggestions or similar experiences would be greatly appreciated.</p>
Avatar
n66134267@gs.ncku.edu.tw March 7, 2026
YOSHIDA UEMORI DEEP DRAWING MODELS
<p>Hello,</p><p>I would like to simulate the process of sheet metal stamping and then springback. In both processes, I would like to take into account the Bausinger effect, anisotropy , softening and workhardening stagnation stress reversals, and strain rate . The Yoshida Uemori model + Anisotropic model, e.g., Yld2000 + strain rate, would be suitable for this task. The problem is that LS Dyna does not have such a model . Therefore, I would like to ask if in the first simulation, i.e., the sheet metal forming simulation, I could use the mat_133 model, then transfer the simulation results to a new springback simulation and change the material model to mat_242, and then run a new springback simulation with the mat_242 material model? The mat_133 model has a backstress historyhttps://lsdyna.ansys.com/history-variables/</p><p> </p><p> </p>
Avatar
euklidespl@gmail.com March 6, 2026
Tied contact failed(DYNAMIC_RELAXATION in the total retart ) Need Advice
<p class="ds-markdown-paragraph">Hello everyone,</p><p class="ds-markdown-paragraph">I am encountering an issue with tied contacts during a total restart analysis for springback calculation and would appreciate some advice. (1)First analysis  A steel formwork falls  onto the rock,used shell element: and used TIED_SHELL_EDGE_TO_SURFACE_BEAM_OFFSET and TIED_SHELL_EDGE_TO_SURFACE_BEAM_OFFSET:<br> calculation in this step is normal,and tied contact works fine: (2)total start used CONTROL_DYNAMIC_RELAXATION keyword to calculate the springback but all tied contact was failed:<br> we could see that penetration occurred between the parts.  Has anyone encountered a similar issue? Any suggestions on what parameters I should check or adjust would be greatly appreciated.</p>
Avatar
383291376@qq.com March 6, 2026
Getting started with usermat
<p>Hello,</p><p>I requested the usermat version of LS-DYNA from my institution's software office, from the link found in here: https://ftp.lstc.com/anonymous/outgoing/support/FAQ/user_defined_materials.faq. However, after receiving the files, the file types were unexpected, and different from what the previous posts and user manual has mentioned. </p><p>They provided me with a folder named "ls-dyna....mpi" and inside were only a license folder, and 3 executable files named l2a, lstc_client, and ls-dyna...mpi. Are these the files I should be expecting, and if so, how do I use them?</p><p>Thanks </p>
Avatar
kshi84@gatech.edu March 5, 2026

Still stuck or need help?

Submit a Support Request