Download - Shaped Charge Penetrating Two Plates
chapter 41: Shaped Charge Penetrating Two Plates
41 Shaped Charge Penetrating Two Plates
Summary 657
Introduction 658
Solution Requirements 659
FEM Solution 660
Results 712
Pre- and Postprocess with SimXpert 663
Input File(s) 714
657CHAPTER 41
Shaped Charge Penetrating Two Plates
SummaryTitle Chapter 41: Shaped Charge Penetrating Two Plates
Features • Wall Boundary of Euler Mesh • Transient Initial Condition of Euler Region• Axis-symmetric Analysis• Structural multi material with shear strength and void
Geometry
Material properties • ExplosiveMilitary Compound B (See EOSIG in MD Nastran QRG)
• Copper
Density = 8960 kg/m3
Shear Modulus = 0.477E11 PaJohnson-Cook Yield ModelMinimum Pressure of Spallation = -2.5E10 Pa
• Steel
Density = 7830 kg/m3
Shear Modulus = 0.818E11 PaEquivalent Yield Stress = 1.4E9 PaMinimum Pressure of Spallation = -3.8E9 Pa
Analysis characteristics Transient explicit dynamic analysis (SOL 700)
Boundary conditions Wall Boundary on the part of Explosive Case
Element types Euler: 8-node solid element for explosive, void, steel, and copper
FE results 1. Snap Shots of Liner Collapse, Jet Formation and Plates Penetrated2. Velocity field of explosive gases, liner, and jet at 20 s
Steel Plates
VoidsCopper Plate
Explosive
MD Demonstration Problems
CHAPTER 41658
Introduction
Figure 41-1 Model
When a metal cone is explosively collapsed onto its axis, a high-velocity rod of molten metal, the jet, is ejected out of the open end of the cone. The cone is called a liner and is typically made of copper. The jet has a mass approximately 20 percent of the cone mass, and elongates rapidly due to its high velocity gradient. This molten rod is followed by the rest of the mass of the collapsed cone, the slug. Typical shaped charges have liner slope angles of less than 42 degrees ensuring the development of a jet; with jet velocities ranging from 3000 to 8000 m/s. A typical construction of a shaped charge is shown in Figure 41-2.
Figure 41-2 Typical Construction of Shaped Charge
659CHAPTER 41
Shaped Charge Penetrating Two Plates
An example simulation of shaped charge formation is carried out to demonstrate the ability of SOL 700 to perform such a simulation. A simplified axisymmetric model of explosives and a copper liner is created in a finite volume Euler mesh. Explosive are detonated starting from a point on the axis of symmetry at the end of the explosives. The simulation is carried out for 60 s after detonation of the explosives. The jet is formed and penetrates two thick plates. See Figure 41-3 for the model layout.
Figure 41-3 SOL 700 Model Setup
Typical shaped charges are axisymmetric. However, aiming at higher velocity, 3-D designs are targeted. 3-D simulation of shaped charge formation would be necessary to avoid excessive experimental work. SOL 700 has full abilities to perform such a 3-D simulation.
Solution RequirementsSOL 700 Model
The model is simplified as shown in Figure 41-3. The aluminum casting is replaced with a rigid body.
Detonation is assumed to start at a point on the axis at the rear end of the explosives. The liner shape is slightly simplified as shown in the figure. The retaining ring is assumed rigid and is modeled as a wall boundary for the Euler Mesh (BARRIER). SI units are used in this example.
MD Demonstration Problems
CHAPTER 41660
FEM SolutionA. Euler Mesh and Liner:
A triangular prismatic Finite Volume Euler mesh is used with head angle of 5 degrees as shown in Figure 41-4. A very fine mesh is used to accurately simulate the behavior of the extremely thin liner. The liner is placed in this Euler mesh. Symmetry conditions (closed boundary, default Euler boundary condition) are imposed on the two rectangular faces of the prism to create an axisymmetric behavior.
Figure 41-4 Euler Mesh
The liner material pressure – density relationship is modeled with EOSPOL model. The liner is made of copper and the constants are taken as follows:
Material yield strength is modeled with a Johnson-Cook yield model. The constants are taken as follows:
a1 1.43E11 N/m2
a2 0.839E11 N/m2
a3 2.16E9 N/m2
b1 0.0
b2 0.0
b3 0.0
A 1.2E8 N/m2
B 1.43E9 N/m2
C 0.0
n 0.5
661CHAPTER 41
Shaped Charge Penetrating Two Plates
Other liner material properties of liner are as follows:
In the input file:
MATDEUL 701 8960. 711 712 713 714EOSPOL, 711, 1.43+11, 0.839+11, 2.16+9SHREL,712,0.477E11$ Johnson-Cook$ A B n C m EPS0 CvYLDJC,713, 1.2E8, 1.43E9, 0.5, 0.0, 1.0, 1.0, 399.0,+$ TMELT TROOM+, 1356.0, 293.0$PMINC,714,-2.5E10
It is very easy to define the shape and position of the liner by using the method of geometrical regions when creating the initial conditions of the liner material.
CYLINDR, 1,, -0.5391, -0.56, 0., 2.0, 0.4147, 0.,++,0.2958CYLINDR, 2,, -0.5391, -0.56, 0., 2.0, 0.4147, 0.,++,0.2939CYLINDR, 3,, 0.2, 2.0406, 0., 0.2047, 2.0406, 0.,++,2.0019TICVAL,2,,DENSITY,8960.
B. Casting and Retaining Ring:
The casting is assumed to be rigid. It is modeled by the default Eulerian boundary condition (closed boundary). The retaining ring is also assumed to be rigid and is modeled by a barrier.
C. Plates:
Two thick plates are placed in this Euler mesh. Plate material is defined as steel:
MATDEUL 801 7830. 811 812 813 814EOSPOL, 811, 1.64E+11SHREL,812,0.818E11YLDVM,813,1.4E9
m 1.0
0 1.0
Tmelt 1356.0 K
Troom 293.0 K
Cv 399.0 J/kg
Density 8960 Kg/m3
Constant shear model 0.477E11 N/m2
Constant spallation model -2.5E10 N/m2
MD Demonstration Problems
CHAPTER 41662
PMINC,814,-3.8E9
The shapes and positions of the plates are defined by using the method of geometrical regions.
CYLINDR, 4,, 0.22, 2.0406, 0., 0.223, 2.0406, 0.,++,2.05CYLINDR, 5,, 0.27, 2.0406, 0., 0.273, 2.0406, 0.,++,2.05TICVAL,3,,DENSITY,7830.
D. Explosive:
The explosive is modeled by ignition and growth equation of state. The explosive is placed in this Euler mesh.
EOSIG,100,,,,,,,,++,,,,,,,,,++,,,,,99,,MCOMPB,SI
The explosive material is taken from the database that is build into SOL 700.
To initialize the whole Euler mesh, a TICEUL entry will be defined.
TICEUL1 1 1TICREG 1 1 ELEM 1 100 1 1.TICREG 2 1 CYLINDER1 701 2 2.TICREG 3 1 CYLINDER2 3.TICREG 4 1 CYLINDER3 701 2 4.TICREG 5 1 CYLINDER4 801 3 5.TICREG 6 1 CYLINDER5 801 3 6.$SET1 1 1 THRU 15342TICVAL,1,,DENSITY,1630.,SIE,4.29E6
663CHAPTER 41
Shaped Charge Penetrating Two Plates
Pre- and Postprocess with SimXpert
Create a New DatabaseEnter the MD Explicit Workspace.
a. Click MD Explicit
b. Click Save As
c. File name, enter CH41
d. Click Save
dc
b
a
MD Demonstration Problems
CHAPTER 41664
Change the Units
a. Tools: Options
b. Select Units Manager
c. Click Standard Units
d. Select the line with m, kg, s, ...
e. Click OK
f. Return to User Options screen and click OK
f
e
d
c
b
a
665CHAPTER 41
Shaped Charge Penetrating Two Plates
Import the Model Geometry
a. File: Import
b. Select Nastran
c. Look in: CHAPTER41
d. Select sch_model.bdf
e. Click Open
e
dc
b
a
MD Demonstration Problems
CHAPTER 41666
Create Explosive Material Compound BIgnition and Growth Equation of State
a. Click: EOS
b. Select [07] EOS Ignition
c. For Name: enter EOSIG_100
d. For MID, enter 100
e. For DBEXP, select MCOMPB
f. For UNITCNV, select SI
g. For ITRMAX, enter 99
h. Click Create
i. EOSIG_100 is added
ih
g
fed
c
b
a
667CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Explosive Material Compound B (continued)Shear Model Explosive
a. Click: Shear
b. Select Elastic Shear Model
c. For Name: enter SHREL_101
d. For MID, enter 101
e. For G, enter 3.E9
f. Click Create
g. SHREL_101 is added
g
f
edc
ba
MD Demonstration Problems
CHAPTER 41668
Create Explosive Material Compound B (continued)Yield Model Explosive
a. Click: Yield
b. Select Von Mises Yield
c. For Name: enter YLDVM_102
d. For MID, enter 102
e. For YIELD, enter 2.E8
f. Click Create
g. YLDVM_102 is added
g
f
edc
b
a
669CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Explosive Material Compound B (continued)Eulerian Material Explosive
a. Click: Eulerian
b. Select Eulerian Material
c. For Name: enter MATDEUL_103
d. For MID, enter 103
e. For RHO, enter 1630
f. Double click EID; select Select
g. For Entity Selection, select EOSIG_100; click OK
h. Double click SID; select Select
i. For Entity Selection, select SHREL_101; click OK
j. Double click YID; select Select
k. For Entity Selection, select YLDVM_102; click OK
l. Click Create
m. MATDEUL_103 is added
m
l
k
j
i
h
g
fedc
ba
MD Demonstration Problems
CHAPTER 41670
Create Material Copper LinerLinear Polynomial Equation of State
a. Click: EOS
b. Select [01] EOS Linear Polynomial
c. For Name: enter EOSPOL_701
d. For MID, enter 701
e. For A1, enter 1.43E11
f. For A2, enter 8.39E10
g. For A3, enter 2.16E9
h. Click Create
i. EOSPOL_701 is added
i
h
gfed
c
ba
671CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Material Copper Liner (continued)Shear Model Copper Liner
a. Click: Shear
b. Select Elastic Shear Model
c. For Name: enter SHREL_702
d. For MID, enter 702
e. For G, enter 4.77E10
f. Click Create
g. SHREL_702 is added
g
f
edc
ba
MD Demonstration Problems
CHAPTER 41672
Create Material Copper Liner (continued)Yield Model Copper Liner
a. Click: Yield
b. Select Johnson-Cook Yield
c. For Name: enter YLDJC_703
d. For MID, enter 703
e. For A, enter 1.2E8
f. For B, enter 1.43E9
g. For N, enter 0.5
h. For CP, enter 399
i. For TM, enter 1356
j. For TR, enter 293
k. Click Create
l. YLDJC_703 is added
l
kji
hgfedc
ba
673CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Material Copper Liner (continued)Spall Limit Copper Liner
a. Click: Spall
b. Select PMINC
c. For Name: enter PMINC_704
d. For MID, enter 704
e. For Value, enter -2.5E10
f. Click Create
g. PMINC_704 is added
fedc
ba
MD Demonstration Problems
CHAPTER 41674
Create Material Copper Liner (continued)Eulerian Material Copper Liner
a. Click: Eulerian
b. Select Eulerian Material
c. For Name: enter MATDEUL_705
d. For MID, enter 705
e. For RHO, enter 8960
f. Double click EID; select Select
g. For Entity Selection, select EOSPOL_701; click OK
h. Double click SID; select Select;
i. For Entity Selection, select SHREL_702 click OK
j. Double click YID; select Select
k. For Entity Selection, select YLJC_703; click OK
l. Double click PID; select Select
m. For Entity Selection, select PMINC_704; click OK
n. Click Create
o. PMINC_704 is added
on
m
l
k
j
i
h
g
fedc
ba
675CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Material Steel PlatesLinear Polynomial Equation of State
a. Click: EOS
b. Select [01] EOS Linear Polynomial
c. For Name: enter EOSPOL_801
d. For MID, enter 801
e. For A1, enter 1.64E11
f. Click Create
g. EOSPOL_801 is added
g
f
edc
ba
MD Demonstration Problems
CHAPTER 41676
Create Material Steel Plates (continued)Shear Model Steel Plates
a. Click: Shear
b. Select Elastic Shear Model
c. For Name: enter SHREL_802
d. For MID, enter 802
e. For G, enter 8.18E10
f. Click Create
g. SHREL_902 is added
g
fedc
ba
677CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Material Steel Plates (continued)Yield Model Steel Plates
a. Click: Yield
b. Select Von Mises Yield
c. For Name: enter YLDVM_803
d. For MID, enter 803
e. For A, enter 1.4E9
f. Click Create
g. YLDJC_803 is added
g
f
edc
b
a
MD Demonstration Problems
CHAPTER 41678
Create Material Steel Plates (continued)Spall Limit Steel Plates
a. Click: Spall
b. Select PMINC
c. For Name: enter PMINC_804
d. For MID, enter 804
e. For Value, enter -3.8E9
f. Click Create
g. PMINC_804 is added
g
f
edc
ba
679CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Material Steel Plates (continued)Eulerian Material Steel Plates
a. Click: Eulerian
b. Select Eulerian Material
c. For Name: enter MATDEUL_805
d. For MID, enter 805
e. For RHO, enter 7830
f. Double click EID; select Select
g. For Entity Selection, select EOSPOL_801; click OK
h. Double click SID; select Select;
i. For Entity Selection, select SHREL_802 click OK
j. Double click YID; select Select
k. For Entity Selection, select YLJC_803; click OK
l. Double click PID; select Select
m. For Entity Selection, select PMINC_804; click OK
n. Click Create
o. PMINC_804 is added
o
n
m
l
k
j
i
h
g
fedc
ba
MD Demonstration Problems
CHAPTER 41680
Create CylindersCylinder 1 defining outer surface of the liner
a. Click: Cylinder
b. Select XYZ
c. For XYZ Input: enter -0.5391 -0.56 0 2 0.4147 0; click OK
d. For Radius, enter 0.2958
e. Click Modify
f. Cylinder_1 is added
f
e
d
c
ba
681CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Cylinders (continued)Cylinder 2 defining inner surface of the liner
a. Click: Cylinder
b. Select XYZ
c. For XYZ Input: enter -0.5391 -0.56 0 2 0.4147 0; click OK
d. For Radius, enter 0.2939
e. Click Modify
f. Cylinder_2 is added
f
e
d
c
ba
MD Demonstration Problems
CHAPTER 41682
Create Cylinders (continued)Cylinder 3 defining the rear end of the liner
a. Click: Cylinder
b. Select XYZ
c. For XYZ Input: enter 0.2 2.0406 0 0.2047 2.0406 0; click OK
d. For Radius, enter 2.0019
e. Click Modify
f. Cylinder_3 is added
f
e
d
c
ba
683CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Cylinders (continued)Cylinder 4 defining the rear end of the liner
a. Click: Cylinder
b. Select XYZ
c. For XYZ Input: enter 0.22 2.0406 0 0.223 2.0406 0; click OK
d. For Radius, enter 2.05
e. Click Modify
f. Cylinder_4 is added
f
e
d
c
ba
MD Demonstration Problems
CHAPTER 41684
Create Cylinders (continued)Cylinder 5 defining the rear end of the liner
a. Click: Cylinder
b. Select XYZ
c. For XYZ Input: enter 0.27 2.0406 0 0.273 2.0406 0; click OK
d. For Radius, enter 2.05
e. Click Modify
f. Cylinder_5 is added
f
e
d
c
ba
685CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Cylinders (continued)Sphere 6 covering the entire model
a. Click: Sphere
b. Select XYZ
c. For XYZ Input: enter 0 0 0; click OK
d. For Radius, enter 1
e. Click Modify
f. Shpere_6 is added
f
e
d
c
ba
MD Demonstration Problems
CHAPTER 41686
Create Initial ValuesInitial values explosive
a. Click: TIC
b. Click TICVAL
c. For ID: enter 1
d. For Title, enter TICVAL_1
e. For Density, enter 1630
f. For SIE, enter 4.2E6
g. Click Modify
h. TICVAL_1 is added
h
g
fe
dc
b
a
687CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Initial Values (continued)Initial values copper liner
a. Click: TIC
b. Click TICVAL
c. For ID: enter 2
d. For Title, enter TICVAL_2
e. For Density, enter 8960
f. Click Modify
g. TICVAL_2 is added
g
fe
dc
b
a
MD Demonstration Problems
CHAPTER 41688
Create Initial Values (continued)Initial values steel plates
a. Click: TIC
b. Click TICVAL
c. For ID: enter 3
d. For Title, enter TICVAL_3
e. For Density, enter 7830
f. Click Modify
g. TICVAL_3 is added
g
fe
dc
b
a
689CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Initial RegionsInitial region explosive
a. Click: TIC
b. Click TICREG
c. For ID: enter 1
d. For Title, enter TICREGL_1
e. Double click VID; select Select
f. For Entity Selection, select Sphere_6; click OK
g. Double click MID; select Select
h. For Entity Selection, select MATDEUL_103 click OK
i. Double click TICID; select Select
j. For Entity Selection, select TICVAL_1; click OK
k. For Level, enter 1
l. Click Modify
m. TICREG_1 is added
m
lk
j
i
h
g
f
edc
b
a
MD Demonstration Problems
CHAPTER 41690
Create Initial Regions (continued)Initial region copper liner
a. Click: TIC
b. Click TICREG
c. For ID: enter 2
d. For Title, enter TICREGL_2
e. Double click VID; select Select
f. For Entity Selection, select Cylinder_1; click OK
g. Double click MID; select Select
h. For Entity Selection, select MATDEUL_705 click OK
i. Double click TICID; select Select
j. For Entity Selection, select TICVAL_2; click OK
k. For Level, enter 2
l. Click Modify
m. TICREG_2 is added
m
lk
j
i
h
g
f
edc
b
a
691CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Initial Regions (continued)Initial region of void
a. Click: TIC
b. Click TICREG
c. For ID: enter 3
d. For Title, enter TICREGL_3
e. Double click VID; select Select
f. For Entity Selection, select Cylinder_2; click OK
g. For Level, enter 3
h. Click Modify
i. TICREG_3 is added
i
h
g
f
edc
b
a
MD Demonstration Problems
CHAPTER 41692
Create Initial Regions (continued)Initial region copper liner
a. Click: TIC
b. Click TICREG
c. For ID: enter 4
d. For Title, enter TICREGL_4
e. Double click VID; select Select
f. For Entity Selection, select Cylinder_3; click OK
g. Double click MID; select Select
h. For Entity Selection, select MATDEUL_705 click OK
i. Double click TICID; select Select
j. For Entity Selection, select TICVAL_2; click OK
k. For Level, enter 4
l. Click Modify
m. TICREG_4 is added
m
lk
j
i
h
g
f
edc
b
a
693CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Initial Regions (continued)Initial region steel plate 1
a. Click: TIC
b. Click TICREG
c. For ID: enter 5
d. For Title, enter TICREGL_5
e. Double click VID; select Select
f. For Entity Selection, select Cylinder_4; click OK
g. Double click MID; select Select
h. For Entity Selection, select MATDEUL_805 click OK
i. Double click TICID; select Select
j. For Entity Selection, select TICVAL_3; click OK
k. For Level, enter 5
l. Click Modify
m. TICREG_5 is added
m
lk
j
i
h
g
f
edc
b
a
MD Demonstration Problems
CHAPTER 41694
Create Initial Regions (continued)Initial region steel plate 2
a. Click: TIC
b. Click TICREG
c. For ID: enter 6
d. For Title, enter TICREGL_6
e. Double click VID; select Select
f. For Entity Selection, select Cylinder_5; click OK
g. Double click MID; select Select
h. For Entity Selection, select MATDEUL_805 click OK
i. Double click TICID; select Select
j. For Entity Selection, select TICVAL_3; click OK
k. For Level, enter 6
l. Click Modify
m. TICREG_6 is added
m
lk
j
i
h
g
f
edc
b
a
695CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Initial Condition Euler
a. Click: TIC
b. Click TICEU1
c. For ID: enter 1
d. For Title, enter TICEUL_1
e. Click NREG; enter 6
f. Click Modify
g. Double click TSID1; select Select
h. For Entity Selection, select TICREG_1 click OK
i. Double click TSID2; select Select
j. For Entity Selection, select TICREG_2 click OK
k. Double click TSID3; select Select
l. For Entity Selection, select TICREG_3 click OK
m. Double click TSID4; select Select
n. For Entity Selection, select TICREG_4 click OK
o. Double click TSID5; select Select
p. For Entity Selection, select TICREG_5 click OK
q. Double click TSID6; select Select
r.;For Entity Selection, select TICREG_6 click OK
s. Click Modify
t. TICEUL_1 is added
t
s
rq
po nm
l
k
j
i
h
g
fe
dc
ba
MD Demonstration Problems
CHAPTER 41696
Create Eulerian Element Property
a. Under Materials and Properties in Properties, click: 3D
b. Click PEULER1
c. For Name: enter PEULER_1
d. For Type, select MMSTREN
e. Double click SID1; select Select
f. For Entity Selection, select TICEUL1_1; click OK
g. Click Create
h. In the Model Browser tree, right click PEULER1
i. Select Properties
j. In the Modify PEULER_1 Property window, click Change Region
k. In the Pick Window, select All
l. Click Done
m. Click Modify
m
lk
j
i
h
g
fedc
b
a
697CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Node Set SegmentsLocate the rear end of the copper liner
a. Zoom around Cylinder_3 area in window
b. Tools: Identify
c. From Pick window Identify Entities, select Nodes
d. Select nodes next to Cylinder_3 (Node 23593); in Pick window Identify Entities, click Exit
e. Assemble: Contact Set
e. Click: Node Set Segment
f. Select five (5) nodes next to Cylinder_3
g. In the Node Set Segment window, for Name:, enter BCSEG_1
h. In the Node Set Segment window, for Node Set:, enter 10
i. Click OK
j. BCSEG_1_1 is added
j
ih
g
f
e
d
c
b
a
MD Demonstration Problems
CHAPTER 41698
Create Barrier
a. LBC tab: Couple: Eulerian
b. Select Barrier
c. From Pick window CREATE BARRIER, select Nodes
d. Click Node 23593
e. Select Plane YZ; click OK
f. For ID: enter 1
g. For Name: enter Barrier_1
h. Double click BCID
i. Select BCSEG_1; click OK
j. Click DIR to unselect
k. Click Modify
l. Barrier_1 is added
l
k
j
i
hgf
ed
c
b
a
699CHAPTER 41
Shaped Charge Penetrating Two Plates
Define Values for PARAMsDefine DYINISTEP parameters
a. Job Parameter: PARAM
b. For Name: enter PARAM_1
c. For SID: enter 1
d. For N: enter DYINISTEP
e. For V1: enter 1.E-11
f. Click Create
g. Click Exit
h. PARAM_1 is added
h
gf
edcb
a
MD Demonstration Problems
CHAPTER 41700
Define Values for PARAMs (continued)Define Results Output Frequency
a. Job Parameter: DYPARAM
b. For Name: enter DYPARAM_1
c. For SID: enter 1
d. For F1: enter LSDYNA
e. For F2: enter BINARY
f. For F3: enter D3PLOT
g. For F4: enter 5.E-6
h. Click Create
i. DYPARAM_1 is added
i
h
gfedc
b
a
701CHAPTER 41
Shaped Charge Penetrating Two Plates
Define Values for PARAMs (continued)Define VELMAX parameter
a. Job Parameter: DYPARAM
b. For Name: DYPARAM_2
c. For SID: 2
d. For F1: enter VELMAX
e. For F2: enter 20.E3
f. Click Create
g. Click Exi
h. DYPARAM_2 is added
h
gf
edc
b
a
MD Demonstration Problems
CHAPTER 41702
Create a New Nastran Job
a. Model Browser: Right click over sch_model.bdf
b. Select Create new Nastran job
c. For Solver Input File, choose Chapter41/SESSION/Chapter41.bdf
d. SXLaunch: For File name: enter Chapter41.bdf
e. Click Save
f. Click OK
fed
c
ba
703CHAPTER 41
Shaped Charge Penetrating Two Plates
Define Load Cases and Export a Nastran Input File
a. Model Browser: Right click over Load Case Control
b. Select Properties
c. Select Subcase Nonlinear Static Parameters
d. For Ending Time, enter 60.E-6
e. For Number of Time Steps, enter 12
f. Click Apply
g. Model Browser: Right click over Displacement Output
h. Click Delete
i. Model Browser: Right click over Element Stress Output
j. Click Delete
k. Right click NewJob
l. Click Export
m. Click Run (optional)
ml
k
ji
hg
f
edc
ba
MD Demonstration Problems
CHAPTER 41704
Run MD Nastran Solver
a. Double click MD Nastran icon
b. Select Chapter41.bdf
c. Click Open
d. Click Run
d
c
ba
705CHAPTER 41
Shaped Charge Penetrating Two Plates
Access the MD Nastran ResultsTo access the results, the ARC file is attached.
a. Under File, select Attach Results
b. File path, select CHAPTER41.DYTR_EULER_0.ARC
c. Click Open
d. Click Apply
Note: If SimX can’t access the results, do the following:
File -> Save
File -> New
File –> Attach Results
Attach Options: BOTH
OK
e d
cb
a
MD Demonstration Problems
CHAPTER 41706
Create Fringe Plot
a. Results: Fringe
b. File path, select CHAPTER41.DYTR_EULER_0.ARC
c. For Result Cases, select Cycle 0, Time 0
d. For Result type, select DENSITY
e. Click Update
d
cb
a
707CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Fringe Plot (continued)Adjust Spectrum Colors
a. Results: Spectrum
b. Spectrum Manager: click Add
c. Spectrum: enter Spectrum_1
d. Click Update
e. Click Calculator
f. Click Colors
g. Click and drag colors from the table to the bar
h. Click Apply
i. Click OK
j. Check the colors
j
i h
g
f
e
dc
ba
MD Demonstration Problems
CHAPTER 41708
Create Fringe Plot (continued)
a. State plot property editor: click Fringe
b. Spectrum Manager: click Add
c. Spectrum range, Spectrum: enter Spectrum_1
d. Click Update; observe graphic
e. Click Plot Data
f. For Result cases, select Cycle 2993, Time 1.0; observe graphic
g. Repeat e. and f. for Time 2, 3, 4, 5, and 6; observe graphics on following page
g
ffe
d cc
b
a
709CHAPTER 41
Shaped Charge Penetrating Two Plates
Create Fringe Plot (continued)
Time = 0
Time = 1.E-5
Time = 2.E-5
Time = 3.E-5
Time = 4.E-5
Time = 5.E-5
Time = 6.E-5
MD Demonstration Problems
CHAPTER 41710
Animate Fringe PlotCreate a Fringe Plot with animation.
a. Results: Fringe
b. State plot property editor: click Fringe
c. For Result cases, select CHAPTER41.DYTR
d. Click Density
e. Select Animate
f. Click Update
f
ed
c
b
a
711CHAPTER 41
Shaped Charge Penetrating Two Plates
Animate Fringe Plot (continued)Create mpeg file
a. State plot property editor: click Animation
Record Attributes, select Movie Filename
c. SimXpert Results Animation File: File name, enter Animation
d. Click Save
e. Click Record Animation button
f. Click Play Animation button
g. Click Stop Animation button
gf
e
dc
b
a
MD Demonstration Problems
CHAPTER 41712
ResultsFigure 41-5 shows the initial position of the copper liner and two thick plates at 0s, snap shots of
liner collapse, jet formation and plates penetrated at 10 s, 20 s, 30 s, 40 s, 50 s and 60 s.
Figure 41-5 Initial Position of the Copper Liner and Two Thick Plates, Snap Shots of Liner Collapse, Jet Formation and Plates Penetrated (Courtesy – Postprocessing by CEI Ensight)
Figure 41-6 shows the velocity field of explosive gases, liner, and jet at 20 s. A jet velocity of about 6000 m/s is achieved
Figure 41-6 Velocity Field of Explosive Gases, Liner, and Jet
Abbreviated SOL 700 Input File
SOL 700,NLTRAN STOP=1CENDTITLE = SHAPED CHARGES TEST$ for QA purpose, run shorter time$ENDTIME = 1.E-5
713CHAPTER 41
Shaped Charge Penetrating Two Plates
IC = 1TSTEPNL=1$$$BEGIN BULKTSTEPNL 1 10 1.E-06 1PARAM*,DYINISTEP,1.E-11PARAM*,DYMINSTEP,1.E-13DYPARAM,VELMAX,20.0E+03DYPARAM,LSDYNA,BINARY,D3PLOT,1.E-5$INCLUDE model.bdfINCLUDE wall.dat$ EXPLOSIVE$MATDEUL 100 1630. 100 101 102$EOSIG,100,,MCOMPB,SI,,,,,++,,,,,,,,,++,,,,,,,,,++,,,,,99$SHREL,101,3.E9$YLDVM,102,2.E8$$ COPPER$MATDEUL 701 8960. 711 712 713 714EOSPOL, 711, 1.43+11, 0.839+11, 2.16+9SHREL,712,0.477E11$ Johnson-Cook $ A B n C m EPS0 CvYLDJC,713, 1.2E8, 1.43E9, 0.5, 0.0, 1.0, 1.0, 399.0,+$ TMELT TROOM+, 1356.0, 293.0 $PMINC,714,-2.5E10$$ STEEL$MATDEUL 801 7830. 811 812 813 814EOSPOL, 811, 1.64E+11SHREL,812,0.818E11YLDVM,813,1.4E9PMINC,814,-3.8E9$TICEUL1 1 1TICREG 1 1 ELEM 1 100 1 1.TICREG 2 1 CYLINDER1 701 2 2.TICREG 3 1 CYLINDER2 3.TICREG 4 1 CYLINDER3 701 2 4.TICREG 5 1 CYLINDER4 801 3 5.TICREG 6 1 CYLINDER5 801 3 6.
MD Demonstration Problems
CHAPTER 41714
$PEULER1, 1 ,, MMSTREN, 1SET1 1 1 THRU 15342CYLINDR, 1,, -0.5391, -0.56, 0., 2.0, 0.4147, 0.,++,0.2958 CYLINDR, 2,, -0.5391, -0.56, 0., 2.0, 0.4147, 0.,++,0.2939CYLINDR, 3,, 0.2, 2.0406, 0., 0.2047, 2.0406, 0.,++,2.0019CYLINDR, 4,, 0.22, 2.0406, 0., 0.223, 2.0406, 0.,++,2.05CYLINDR, 5,, 0.27, 2.0406, 0., 0.273, 2.0406, 0.,++,2.05$TICVAL,1,,DENSITY,1630.,SIE,4.29E6TICVAL,2,,DENSITY,8960.TICVAL,3,,DENSITY,7830.$BARRIER,1,2$ENDDATA
Input File(s)
File Description
nug_41.dat MD Nastran input file for wall boundary of Euler element