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


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.


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


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


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


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


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


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


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


Create Explosive Material Compound B (continued)Shear Model Explosive

a. Click: Shear

b. Select Elastic Shear Model

c. For Name: enter SHREL_101


e. For G, enter 3.E9

f. Click Create

g. SHREL_101 is added

g

f

edc

ba


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


e. For YIELD, enter 2.E8

f. Click Create

g. YLDVM_102 is added

g

f

edc

b

a

669CHAPTER 41


Create Explosive Material Compound B (continued)Eulerian Material Explosive

a. Click: Eulerian

b. Select Eulerian Material

c. For Name: enter MATDEUL_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


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


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


Create Material Copper Liner (continued)Shear Model Copper Liner

a. Click: Shear




e. For G, enter 4.77E10

f. Click Create


g

f

edc

ba


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


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


Create Material Copper Liner (continued)Spall Limit Copper Liner

a. Click: Spall

b. Select PMINC

c. For Name: enter PMINC_704


e. For Value, enter -2.5E10

f. Click Create

g. PMINC_704 is added

fedc

ba


CHAPTER 41674

Create Material Copper Liner (continued)Eulerian Material Copper Liner

a. Click: Eulerian






g. For Entity Selection, select EOSPOL_701; click OK

h. Double click SID; select Select;

i. For Entity Selection, select SHREL_702 click OK


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


Create Material Steel PlatesLinear Polynomial Equation of State

a. Click: EOS

b. Select [01] EOS Linear Polynomial

c. For Name: enter EOSPOL_801


e. For A1, enter 1.64E11

f. Click Create

g. EOSPOL_801 is added

g

f

edc

ba


CHAPTER 41676

Create Material Steel Plates (continued)Shear Model Steel Plates

a. Click: Shear




e. For G, enter 8.18E10

f. Click Create


g

fedc

ba

677CHAPTER 41


Create Material Steel Plates (continued)Yield Model Steel Plates

a. Click: Yield

b. Select Von Mises Yield

c. For Name: enter YLDVM_803


e. For A, enter 1.4E9

f. Click Create

g. YLDJC_803 is added

g

f

edc

b

a


CHAPTER 41678

Create Material Steel Plates (continued)Spall Limit Steel Plates

a. Click: Spall

b. Select PMINC

c. For Name: enter PMINC_804


e. For Value, enter -3.8E9

f. Click Create

g. PMINC_804 is added

g

f

edc

ba

679CHAPTER 41


Create Material Steel Plates (continued)Eulerian Material Steel Plates

a. Click: Eulerian






g. For Entity Selection, select EOSPOL_801; click OK

h. Double click SID; select Select;

i. For Entity Selection, select SHREL_802 click OK


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


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


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


e. Click Modify


f

e

d

c

ba


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


e. Click Modify


f

e

d

c

ba

683CHAPTER 41



a. Click: Cylinder

b. Select XYZ



e. Click Modify


f

e

d

c

ba


CHAPTER 41684


a. Click: Cylinder

b. Select XYZ



e. Click Modify


f

e

d

c

ba

685CHAPTER 41


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


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


Create Initial Values (continued)Initial values copper liner

a. Click: TIC

b. Click TICVAL

c. For ID: enter 2



f. Click Modify

g. TICVAL_2 is added

g

fe

dc

b

a


CHAPTER 41688

Create Initial Values (continued)Initial values steel plates

a. Click: TIC

b. Click TICVAL

c. For ID: enter 3



f. Click Modify

g. TICVAL_3 is added

g

fe

dc

b

a

689CHAPTER 41


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


CHAPTER 41690

Create Initial Regions (continued)Initial region copper liner

a. Click: TIC

b. Click TICREG

c. For ID: enter 2



f. For Entity Selection, select Cylinder_1; click OK






l. Click Modify


m

lk

j

i

h

g

f

edc

b

a

691CHAPTER 41


Create Initial Regions (continued)Initial region of void

a. Click: TIC

b. Click TICREG

c. For ID: enter 3




g. For Level, enter 3

h. Click Modify

i. TICREG_3 is added

i

h

g

f

edc

b

a


CHAPTER 41692

Create Initial Regions (continued)Initial region copper liner

a. Click: TIC

b. Click TICREG

c. For ID: enter 4









l. Click Modify


m

lk

j

i

h

g

f

edc

b

a

693CHAPTER 41


Create Initial Regions (continued)Initial region steel plate 1

a. Click: TIC

b. Click TICREG

c. For ID: enter 5









l. Click Modify


m

lk

j

i

h

g

f

edc

b

a


CHAPTER 41694

Create Initial Regions (continued)Initial region steel plate 2

a. Click: TIC

b. Click TICREG

c. For ID: enter 6









l. Click Modify


m

lk

j

i

h

g

f

edc

b

a

695CHAPTER 41


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


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


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


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


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


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


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


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


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


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


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


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


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


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


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


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


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


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


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.


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

http://www.mscsoftware.com/doc/nastran/mdug/input_files/ch041/nug_41.dat

shaped charge penetrating two plates

Documents