Workshop 12

Urban Blast

Introduction to ANSYS AUTODYN

1 © 2011 ANSYS, Inc. March 9, 2012
 Workshop Goal and Procedure 
Goal: 
Model a Explosive blast in an urban 
environment

Procedure:
Start AUTODYN standalone
Set up and solve a problem using 2D 
Axial Symmetry to model the 
initial spherical blast
Set up a 3D problem to model the 
urban area of the blast
Remap the results from the 2D 
model into the 3D model as initial 
conditions
Run the 3D solve

2 © 2011 ANSYS, Inc. March 9, 2012
 Quick Instructions

• Split into two problems
– 2D axisymmetric problem to model intial spherical 
blast to create remap file
– Remap file read in to provide initial conditions for the 
3D analysis
3 © 2011 ANSYS, Inc. March 9, 2012
 Quick Instructions

• Start AUTODYN and create a new 2D axially symmetric 
model with default units of mm, mg, ms
• Load in air and TNT from the material library
• Set up the 2D multi material Euler wedge predef as shown 
below

Air, internal energy

TNT, filled using an = 2.068x105
ellipse with x and y
axis of 1000 mm

i=1 i = 451
Min radius Max radius =
= 100 mm 45,000 mm

4 © 2011 ANSYS, Inc. March 9, 2012
 Quick Instructions

• Create a detonation point at origin (0.0,0.0)
• Set the wrap up criteria to a time of 50 ms
• Run the model to completion, which should take approximately 
3,500 cycles
• Create the .fil file by selecting write datafile under additional fill 
options and enter the filename urban_blast_wedge
• This creates a file urban_blast_wedge.fil in the working directory

5 © 2011 ANSYS, Inc. March 9, 2012
 Quick Instructions

• Create a new AUTODYN model file called urban blast with 
3D symmetry and default units (mm,mg,ms)
• Load air from the material library
• Create a flow_out Euler boundary condition called outflow 
with preferred material of air
• Create a 3D Euler, Ideal Gas part using manual definition 
with maximum i,j,k indices of 81,81,61
• Select fill block from within the zoning options and use the 
explicit file type option and select the file buildings.blk
• Load in the settings file urban_blast.set for an appropriate 
view of what’s been generated
• Switch on gauges in the plot menu to see the gauges that 
have been included
6 © 2011 ANSYS, Inc. March 9, 2012
 Quick Instructions

• Read in the previously written urban_blast_wedge.fil datafile to 
create the appropriate initial conditions of the explosive
• Select both materials to remap
• Select a remap origin of (300,000.0,225,000.0,0.0)
• Respond yes to all warnings as data is processed
• Apply the flow_out boundary conditions to all outer faces of the 
block apart from that of k=1

7 © 2011 ANSYS, Inc. March 9, 2012
 Quick Instructions

• Select boundaries in the plot 
section to ensure they have 
been correctly applied
• Set a wrap up criteria with a 
time limit of 1000ms
• Save the file every 25 ms
• Run the simulation and 
create some animations

8 © 2011 ANSYS, Inc. March 9, 2012
 Step 1 – Start AUTODYN standalone
1. Start AUTODYN standalone (Double click autodyn.exe)

9 © 2011 ANSYS, Inc. March 9, 2012
 Step 2 ‐ Initiate a new project for the Spherical Blast (2D)
2.a  Start a new project

2.b  Browse to select a project 
Folder 

2.c  Enter an Ident 
“spherical_blast” and a 
Heading

2.d  Select 2D Axial Symmetry

2.e  Keep the default set of units 
(recommended for 
accuracy)

2.f  Accept the data
10 © 2011 ANSYS, Inc. March 9, 2012
 Step 3 – Retrieve material data from library
3.a  Load material data

3.b  Select the following 
materials from the 
library:

AIR

TNT

3.c  Accept the data

<control>+<click> to select both materials

11 © 2011 ANSYS, Inc. March 9, 2012
 Step 4 – Generate an Euler Part (Wedge)
4.c Enter “space” for 
the Part Name

4.d  Select the Euler, 
2D Multi‐
material Solver

4.e  Keep Part wizard
4.a Select
Parts On 4.f  Move to the next 
Navigation step
Bar

4.b Select
New in the
dialog panel
12 © 2011 ANSYS, Inc. March 9, 2012
 Step 5 – Define Geometry and Mesh for the Euler Part
5.a  Select the Wedge Predef

5.b  Enter “100.0” for the 
Minimum radius and 
“45000.0” for the Maximum 
radius

5.c  Move to the next step

5.d  Enter “450” for the number of 
Cells across the radius

5.e  Move to the next step

13 © 2011 ANSYS, Inc. March 9, 2012
 Step 6 – Fill the Euler Part with Air (1 bar)
6.a  DO NOT Check Fill with Initial 
Condition Set

6.b  Select AIR for the material 
(you fill the entire Euler mesh 
with AIR now, then refill next 
with TNT)

6.c  Keep the default density for 
air but enter “2.068e5” for the 
internal energy (this initializes 
the air to 1 bar pressure)

6.d  Accept the data

14 © 2011 ANSYS, Inc. March 9, 2012
 Step 7 – Fill  Euler mesh with TNT
7.a Select 7.d  Enter 
Parts on the “1000.0” for 
Navigation the X and Y 
Bar semi‐axis of 
the ellipse 
7.b Select Fill (i.e. a circle)

7.c Select the 7.e  Select TNT 

Ellipse fill for the 
option material

7.f  Accept the 
data

15 © 2011 ANSYS, Inc. March 9, 2012
 Step 8 – Define Detonation Point for the TNT
8.a  Select Detonation on the 
Navigation Bar
8.b  Select Point
8.c  Accept the data
(the defaults place a detonation 
point at the origin and the  
detonation starts at time 
zero)
8.d  Check Plot detonation 
points to see the placement 
of the point.

16 © 2011 ANSYS, Inc. March 9, 2012
 Step 9 – Set Solution Controls
9.a  Select Controls on the 
Navigation Bar

9.b  Enter a large number 
for the Cycle limit and 
“50.0” ms for the 
Time Limit 
(this is the time it takes for 
the blast wave to 
propagate to a radius 
of 45m)

17 © 2011 ANSYS, Inc. March 9, 2012
 Step 10 – Load a Plot Settings File (3D View)
10.a  Load the plot settings 
file “spherical_blast.set”

10.b Notice that a Contour
plot of Pressure has 
been specified

10.c  Notice that the loaded 
plot settings provide a 
3D view where the 
model is rotated through 
360 deg.

10.d  Extra credit: Try creating 
this view yourself 
without loading the 
settings file 

18 © 2011 ANSYS, Inc. March 9, 2012
 Step 11 – Solve the Problem
11.a  Save your Data During the solve, the view is
automatically updated to show the most
current solution
11.b  Solve the 
Notice that when the Solve is complete
problem
the pressure (blast) wave has almost
reached the outer radius of the wedge

19 © 2011 ANSYS, Inc. March 9, 2012
 Step 12 – Write the solution to a Remap File (.fil)
12.a Select Parts
on the
Navigation Bar
12.b Select Fill
12.c Select the
Datafile
12.d Select Write
Datafile
12.e Keep the
default
Filename
12.f Accept the
data

20 © 2011 ANSYS, Inc. March 9, 2012
 Step 13 ‐ Initiate a new project for the Urban Blast (3D)
13.a  Start a new project

13.b  Keep the same project 
Folder 

13.c  Enter an Ident
“urban_blast” and a 
Heading

13.d  Select 3D Symmetry

13.e  Keep the default set of 
units

13.f  Accept the data

21 © 2011 ANSYS, Inc. March 9, 2012
 Step 14 – Retrieve material data from library
14.a  Load material 
data

14.b  Select AIR from 
the library

14.c  Accept the data

22 © 2011 ANSYS, Inc. March 9, 2012
 Step 15 – Define an Outflow Boundary
15.a Select Boundaries on
the Navigation Bar

15.b Select New

15.c Enter “Outflow” for

the Boundary Name

15.d Select Flow_Out for

the Type and keep Flow
out (Euler) for the Sub
option

15.e Keep AIR as the

Preferred Material

15.f Accept the data

23 © 2011 ANSYS, Inc. March 9, 2012
 Step 16 – Generate an Euler Part
16.c Enter “Environment” 
for the Part Name

16.d  Select the Euler, 
Ideal Gas Solver

16.e  Select Manual 
Definition

16.f  Enter (81,81,61) for 
16.a Select the maximum (I,J,K) 
Parts On indices (this is the 
Navigation range for the data  
Bar imported during the 
next step)

16.b Select 16.g  Accept the data

New in the
dialog panel
24 © 2011 ANSYS, Inc. March 9, 2012
 Step 17 – Import data for the Euler Part
17.a Select Zoning
17.c Select Explicit
for the File Type

17.d Browse and select 17.e Accept the data

17.b Select Fill Block
the file “buildings.blk”
25 © 2011 ANSYS, Inc. March 9, 2012
 Step 18 – Load a Plot Settings File
18.a  Load the plot 
settings file 
“urban_blast.set”

18.c  Check Gauges
under Additional 
Components

‐ This shows the 
location of three 
gauges that have 
been imported 
along with the 
model
When the Fill Block operation was performed, cells representing the buildings
were automatically defined as unused. All other cells were automatically filled
with air at room temperature and pressure.
26 © 2011 ANSYS, Inc. March 9, 2012
 Step 19 – Remap the Spherical Blast Results into 3D
19.a  Select Parts on the Navigation 
bar

19.b  Select (Fill) Datafile

19.c  Select Read Datafile

19.d  Browse and select the file 
“spherical_blast.fil”

19.e  Set the remap origine to 
(300000.0,225000.0,0.0)

19.f  Select AIR and TNT as materials 

to be remapped

19.g  Accept the data (respond yes to 
27 © 2011 ANSYS, Inc. March 9, 2012 all the warning messages)
 Step 20 ‐ Apply the Outflow Boundary
The quickest way to apply the Outflow boundary to the outer faces of the mesh is 
to apply it to the entire mesh and then clear it from all internal faces
20.a Select Parts on 20.b Select Boundary
the Navigation bar
20.c Select Block
20.d Accept the
defaults to apply the
Outflow boundary to
the entire mesh
20.e Select Clear
20.f Enter I=2-80, J=2-
80, K=1-60 to remove
the boundary from all
internal faces and
faces on the ground
20.g Accept the data
28 © 2011 ANSYS, Inc. March 9, 2012
 Step 21 – Set Solution and Output Controls
21.a  Select Controls on 
the Navigation Bar
21.b  Enter a large number 
for the Cycle limit and 
“1000.0” ms for the 
Time Limit 
(this is the time it 
takes for the blast 
wave to propagate 
through the urban 
environment)
21.c  Select Output on the 
Navigation Bar

21.d  Write Save files 
every 25.0 ms (set 
End time to a large 
number) 

29 © 2011 ANSYS, Inc. March 9, 2012
 Step 22 – Solve the Problem
22.a  Save your Data During the solve, the view is automatically
updated to show the most current solution
The calculation runs for about 25 minutes
22.b  Solve the 
problem Press Stop to stop the solve at any time (then
press Run if you want to continue)

30 © 2011 ANSYS, Inc. March 9, 2012
 Step 23 – View Pressure Histories at Gauges

23.a  After the Solve is 
complete, select History
on the Navigation Bar
23.b  Select Single Variable 
Plots
23.c  Select one or more 
gauges to view
23.d  Select Pressure as the Y 
variable
23.e  Accept the data
23.f  Try other options:

View different variables

Integrate Pressure to obtain 
Impulse

Select Scales

31 © 2011 ANSYS, Inc. March 9, 2012