Problem description: An explosive material is placed inside of a sphere. The explosive is lit, expands, and impacts the sphere. The sphere expands, reaches yield, seam lines fail (failure strain 0.2) and fragments fly apart. Due to use of high explosives time step must be scaled down for stability.
*CONTACT_SLIDING_ONLY *CONTROL_ENERGY *CONTROL_OUTPUT *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_BINARY_D3THDT *DATABASE_EXTENT_BINARY *DATABASE_GLSTAT *DATABASE_MATSUM *DATABASE_RCFORC *ELEMENT_SHELL *ELEMENT_SOLID *END *EOS_JWL *INITIAL_DETONATION *KEYWORD *MAT_HIGH_EXPLOSIVE_BURN *MAT_PLASTIC_KINEMATIC *NODE *PART *SECTION_SHELL *SECTION_SOLID *SET_SEGMENT *TITLE
*KEYWORD
*TITLE
an exploding sphere
$
$ John D. Reid 1/31/00 - original model obtained from K. Bui (LSTC) - 4/22/97
$
$ Units: unknown - possibly gm, cm, micros, 1E7N, Mbar, 1E7N-cm
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Control Ouput
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
*CONTROL_TERMINATION
$ endtim endcyc dtmin endeng endmas
50.0
$
*CONTROL_TIMESTEP
$ dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
0.30
$
*CONTROL_ENERGY
$ hgen rwen slnten rylen
2 2 2
$
*CONTROL_OUTPUT
$ npopt neecho nrefup iaccop opifs ipnint ikedit
1 3
$
$
*DATABASE_BINARY_D3PLOT
$ dt
1.0
$
$*DATABASE_EXTENT_BINARY
$ neiph neips maxint strflg sigflg epsflg rltflg engflg
$
$ cmpflg ieverp beamip
$ 1
$
*DATABASE_BINARY_D3THDT
$ dt
999999.9
$
*DATABASE_GLSTAT
$ dt
0.1
$
*DATABASE_MATSUM
$ dt
0.1
$
*DATABASE_RCFORC
$ dt
0.1
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Boundary and Initial Conditions
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
*INITIAL_DETONATION
$ pid x y z lt
2 0.00 0.00 0.00 0.00
$
$ pid = 2 Part ID of high explosive material to be lit.
$ x,y,z Coordinate of detonation point.
$ lt = 0 Lighting time.
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Define Parts and Materials
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
$$$ Part 1 shell: sphere
$
$$$ Part 2 solid: high explosive
$
*PART
$ pid sid mid eosid hgid grav adpopt
sphere
1 1 1
high explosive
2 2 2 2
seams
3 1 3
$
$
$$$$ Materials
$
*MAT_PLASTIC_KINEMATIC
$
$ mid ro e pr sigy etan beta
1 7.874 2.07 0.28 0.15 0.136 0
$
$ src srp fs
40 5
$
*MAT_PLASTIC_KINEMATIC
$
$ mid ro e pr sigy etan beta
3 7.874 2.07 0.28 0.15 0.136 0
$
$ src srp fs
40 5 0.2
$
$
*MAT_HIGH_EXPLOSIVE_BURN
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$ mid ro d pcj beta k g sigy
2 1.717 0.798 0.290
$
$ mid = 2 Material ID
$ ro = 1.717 mass density
$ d = 0.798 detonation velocity
$ pcj = 0.290 Chapman-Jouget pressure
$
$
$$$$ Sections
$
*SECTION_SHELL
$ sid elform shrf nip propt qr/irid icomp
1
$ t1 t2 t3 t4 nloc
0.50 0.50 0.50 0.50
$
*SECTION_SOLID
$ sid elform
2 1
$
$
$
*EOS_JWL
$ eosid A B r1 r2 omeg Eo Vo
2 5.2421999 0.07678 4.2000003 1.10 0.340 0.085 1.00
$
$ Defines the pressure in the solid elements as a function
$ of the relative volume and internal energy.
$
$ Eo = 0.085 initial internal energy
$ Vo = 1.000 initial relative volume
$ A, B, r1, r2, omeg = parameters of equation of state
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Contact
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ *CONTACT_SLIDING_ONLY - defined at end of model
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ Define Nodes and Elements
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$
*CONTACT_SLIDING_ONLY
$HMNAME GROUPS 1ContSlide_1
$HMCOLOR GROUPS 1 1
1 2 0 0
*SET_SEGMENT
1
*SET_SEGMENT
2
$
*END
