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

Contact doesn’t work
<p>Hi Everyone</p><p>I created a dome composed of many beams as shown in the figure 1 and 2, and made a cylinder imapct it. I applied *CONTACT_AUTOMATIC_BEAMS_TO_SURFACE_ID between the cylinder and each beam. As for beams, there is a little difference in height between 2 beams. I didn't add any contacts. I added *CONSTRAINED_GENERALIZED_WELD_SPOT or *CONSTRAINED_SPOTWELD between the nodes of 2 beams closest to the intersection. After I ran the simulation, the cylinder went through the dome rather than impacting it, as shown in figure 3. The beams appeared to have a little oscillation. It seems that there is still some interaction force between the beam and the cylinder. But the cylinder just directly went through the beams. Does anybody know the reason and how to fix it? </p><p>I attached my .k file : https://tuprd-my.sharepoint.com/:u:/g/personal/tut06916_temple_edu/IQAtEDBPZhKLRYw4Kr0EEa_AAeYXal-n6FkEF1ROSpG_67k?e=86nLOW</p><p></p><p></p><p></p>
Avatar
tut06916@temple.edu March 12, 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
Contact doesn’t work
<p>Hi Everyone</p><p>I created a dome composed of many beams as shown in the figure 1 and 2, and made a cylinder imapct it. I applied *CONTACT_AUTOMATIC_BEAMS_TO_SURFACE_ID between the cylinder and each beam. As for beams, there is a little difference in height between 2 beams. I didn't add any contacts. I added *CONSTRAINED_GENERALIZED_WELD_SPOT or *CONSTRAINED_SPOTWELD between the nodes of 2 beams closest to the intersection. After I ran the simulation, the cylinder went through the dome rather than impacting it, as shown in figure 3. The beams appeared to have a little oscillation. It seems that there is still some interaction force between the beam and the cylinder. But the cylinder just directly went through the beams. Does anybody know the reason and how to fix it? </p><p>I attached my .k file : https://tuprd-my.sharepoint.com/:u:/g/personal/tut06916_temple_edu/IQAtEDBPZhKLRYw4Kr0EEa_AAeYXal-n6FkEF1ROSpG_67k?e=86nLOW</p><p></p><p></p><p></p>
Avatar
tut06916@temple.edu March 12, 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

Still stuck or need help?

Submit a Support Request