TNT model Validation

20 views
Skip to first unread message

Shubham Sharma

unread,
Apr 27, 2026, 8:59:41 PM (6 days ago) Apr 27
to LS-DYNA2
Hi all,
I am validating a 2D axisymmetric TNT model so it can be used to map 2D to 3D ale simulations.

 The expected overpressure at 400mm standoff should be around 3 MPa.
 The research paper that I am validating has a cylindrical TNT with a diameter-to-height ratio of 2.

When IVFG CONTTYP 4 was used to model a cylindrical TNT  in 2d axisymmetric model, it was not giving the correct blast wave propagation. Instead, CONTTYP 5 was used to define a rectangle in 2D with length and breadth proportional to the cylinder, but this simulation consistently overpredicts the pressure by at least a factor of 3, I am attaching .k file so you can have a look.

I need to understand why the cylindrical TNT made through IVFG does not give the desired wave propagation and also the rectangular charge's overprediction of blast overpressure. Also when spherical shape was used it gave overpressure values to some accuracy.

mapping_2d_rectangular.zip

l...@schwer.net

unread,
Apr 27, 2026, 10:19:33 PM (6 days ago) Apr 27
to Shubham Sharma, LS-DYNA2

Have you read this Remark 6 under the keyword *INITIAL_VOLUME_FRACTION_GEOMETRY:

 

“Two-dimensional geometry. If the ALE model is 2D (*SECTION_ALE2D instead of *SECTION_SOLID), CNTTYP = 4 defines a quadrilateral. In this case, the fields which, in the 3D case define a cone, are interpreted as the corner coordinates of a clockwise defined (inward normal) quadrilateral having the vertices: (X1, Y1), (X2, Y2), (X3, Y3), and (X4, Y4). The CNTTYP = 4 input fields X0, Y0, Z0, X1, Y1, Z1, R1, and R2 becomes X1, Y1, X2, Y2, X3, Y3, X4, and Y4 respectively. CNTTYP = 6 should be used to fill a circle.”

 

               --len

CONFIDENTIALITY NOTICE: The contents of this email message and any attachments are intended solely for the addressee(s) and may contain confidential and/or privileged information and may be legally protected from disclosure. If you are not the intended recipient of this message or their agent, or if this message has been addressed to you in error, please immediately alert the sender by reply email and then delete this message and any attachments. If you are not the intended recipient, you are hereby notified that any use, dissemination, copying, or storage of this message or its attachments is strictly prohibited.

--
You received this message because you are subscribed to the Google Groups "LS-DYNA2" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ls-dyna2+u...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/ls-dyna2/d88f91f3-7b71-4b2e-9d11-c91e75854840n%40googlegroups.com.

Shubham Sharma

unread,
Apr 28, 2026, 6:44:35 AM (5 days ago) Apr 28
to LS-DYNA2
I tested giving coordinates to IVFG according to the 2d ale, and it is giving the same result as given by CONTYP5 that is a rectangle, the issue I am facing is that for same weight and same scaled distance the spherical charge gives a peak overpressure value of 3 Mpa which is closer to the research paper's theoretical value of 2.88 Mpa, but the research paper uses cylindrical charge, and when I am using the cylindrical charge it is giving a value of 8.4 Mpa (for same weight and same scaled distance).
Reply all
Reply to author
Forward
0 new messages