*HEADING
 Dynamic Axial Crushing of a Double Chambered Aluminum Extrusion
**
*PREPRINT, ECHO=YES, HISTORY=YES, CONTACT=YES, MODEL=YES
*RESTART, WRITE
**
*********************** M O D E L     D E F I N I T I O N  *********************
**
** --------------------  N O D A L    P O I N T    D A T A  --------------------
**
*Node, Nset=All_Nodes
**
** ------------------------  E L E M E N T     D A T A  ------------------------
**
*Element, Type=SFM3D4R, Elset=P2_BaseSupport
**
*Surface section, Elset=P2_BaseSupport, DENSITY=7.8200e-09
**
*Element, Type=SFM3D4R, Elset=P3_Impactor
**
*Surface section, Elset=P3_Impactor, DENSITY=7.8200e-09
**
*Element, Type=S3R, Elset=P1_DoubleChamberExtrusion
**
*Element, Type=S4R, Elset=P1_DoubleChamberExtrusion
**
*Shell section, Elset=P1_DoubleChamberExtrusion, Material=AluminumAlloy
  2.5000e+00,     5
**
** -----------------------  M A T E R I A L     D A T A  -----------------------
**
** Double Chamber Extrusion
**
*Material, Name=AluminumAlloy
*Density
2.70E-09
*Elastic
7.00E+04,	0.33
*plastic,hardening=isotropic,rate=0
3.1100E+02,	0.0000E+00
3.1664E+02,	3.4373E-03
3.2452E+02,	8.5932E-03
3.3520E+02,	1.6327E-02
3.4896E+02,	2.7928E-02
3.6552E+02,	4.5329E-02
3.8370E+02,	7.1431E-02
4.0165E+02,	1.1058E-01
*plastic,hardening=isotropic,rate=2.00e-03
3.1100E+02,	0.0000E+00
3.1664E+02,	3.4373E-03
3.2452E+02,	8.5932E-03
3.3520E+02,	1.6327E-02
3.4896E+02,	2.7928E-02
3.6552E+02,	4.5329E-02
3.8370E+02,	7.1431E-02
4.0165E+02,	1.1058E-01
*plastic,hardening=isotropic,rate=1.00e-02
3.1254E+02,	0.0000E+00
3.1817E+02,	3.4373E-03
3.2606E+02,	8.5932E-03
3.3673E+02,	1.6327E-02
3.5050E+02,	2.7928E-02
3.6706E+02,	4.5329E-02
3.8523E+02,	7.1431E-02
4.0318E+02,	1.1058E-01
*plastic,hardening=isotropic,rate=1.00e-01
3.1577E+02,	0.0000E+00
3.2140E+02,	3.4373E-03
3.2929E+02,	8.5932E-03
3.3996E+02,	1.6327E-02
3.5372E+02,	2.7928E-02
3.7028E+02,	4.5329E-02
3.8845E+02,	7.1431E-02
4.0639E+02,	1.1058E-01
*plastic,hardening=isotropic,rate=1.00e+00
3.2083E+02,	0.0000E+00
3.2646E+02,	3.4373E-03
3.3434E+02,	8.5932E-03
3.4500E+02,	1.6327E-02
3.5874E+02,	2.7928E-02
3.7526E+02,	4.5329E-02
3.9338E+02,	7.1431E-02
4.1123E+02,	1.1058E-01
*plastic,hardening=isotropic,rate=1.00e+01
3.2877E+02,	0.0000E+00
3.3433E+02,	3.4373E-03
3.4210E+02,	8.5932E-03
3.5258E+02,	1.6327E-02
3.6603E+02,	2.7928E-02
3.8208E+02,	4.5329E-02
3.9940E+02,	7.1431E-02
4.1597E+02,	1.1058E-01
*plastic,hardening=isotropic,rate=2.50e+01
3.3306E+02,	0.0000E+00
3.3846E+02,	3.4373E-03
3.4597E+02,	8.5932E-03
3.5604E+02,	1.6327E-02
3.6882E+02,	2.7928E-02
3.8376E+02,	4.5329E-02
3.9929E+02,	7.1431E-02
4.1299E+02,	1.1058E-01
*plastic,hardening=isotropic,rate=1.00e+02
3.4119E+02,	0.0000E+00
3.4639E+02,	3.4373E-03
3.5359E+02,	8.5932E-03
3.6316E+02,	1.6327E-02
3.7513E+02,	2.7928E-02
3.8877E+02,	4.5329E-02
4.0220E+02,	7.1431E-02
4.1262E+02,	1.1058E-01
*plastic,hardening=isotropic,rate=1.00e+03
3.4568E+02,	0.0000E+00
3.5083E+02,	3.4373E-03
3.5796E+02,	8.5932E-03
3.6741E+02,	1.6327E-02
3.7920E+02,	2.7928E-02
3.9254E+02,	4.5329E-02
4.0550E+02,	7.1431E-02
4.1519E+02,	1.1058E-01
*plastic,hardening=isotropic,rate=1.00e+04
3.4568E+02,	0.0000E+00
3.5083E+02,	3.4373E-03
3.5796E+02,	8.5932E-03
3.6741E+02,	1.6327E-02
3.7920E+02,	2.7928E-02
3.9254E+02,	4.5329E-02
4.0550E+02,	7.1431E-02
4.1519E+02,	1.1058E-01
*DAMAGE INITIATION, CRITERION=DUCTILE
3.3238E+01,	-3.3333E+00,	1.00E-04
3.3238E+01,	-3.3333E-01,	1.00E-04
2.3381E+01,	-2.6667E-01,	1.00E-04
1.6447E+01,	-2.0000E-01,	1.00E-04
1.1570E+01,	-1.3333E-01,	1.00E-04
8.1394E+00,	-6.6667E-02,	1.00E-04
5.7268E+00,	0.0000E+00,	1.00E-04
4.0303E+00,	6.6667E-02,	1.00E-04
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT
0.1
*DAMAGE INITIATION, CRITERION=SHEAR, ks=0.3
2.7610E-01,	-1.0000E+01,	1.00E-04
2.7610E-01,	1.4236E+00,	1.00E-04
2.6130E-01,	1.4625E+00,	1.00E-04
2.5300E-01,	1.5013E+00,	1.00E-04
2.5100E-01,	1.5401E+00,	1.00E-04
2.5510E-01,	1.5789E+00,	1.00E-04
2.6560E-01,	1.6177E+00,	1.00E-04
2.8250E-01,	1.6566E+00,	1.00E-04
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT
0.1
*DAMAGE INITIATION, CRITERION=MSFLD, DEFINITION=FLD, OMEGA=0.001
2.0370E-01,	-1.0185E-01
1.6590E-01,	-6.9125E-02
1.4548E-01,	-4.8493E-02
1.2454E-01,	-3.1136E-02
1.1234E-01,	-1.8723E-02
1.0013E-01,	-8.3438E-03
8.8206E-02,	0.0000E+00
9.2471E-02,	7.1132E-03
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT
0.1
**
**
** ---------------------  R I G I D    B O D Y    D A T A  --------------------
**
** Impactor
**
*Elset, Elset=RBody_Impactor
P3_Impactor,
**
*Nset, Nset=RBody_Impactor_RefNode
**
*Rigid body, Ref node=RBody_Impactor_RefNode, Position=CENTER OF MASS, Elset=RBody_Impactor
**
*ELEMENT, TYPE=MASS, ELSET=MASSE_500_kg
*MASS, ELSET=MASSE_500_kg
 0.500,
*ELEMENT, TYPE=ROTARYI, ELSET=MASSE_500_kg_Rot_Inertia
*ROTARY INERTIA, ELSET=MASSE_500_kg_Rot_Inertia
 0.2088e+005,  0.1431e+005,  0.1430e+005
**
** Base Support
**
*Elset, Elset=RBody_BaseSupport
P2_BaseSupport,
**
*Nset, Nset=RBody_BaseSupport_RefNode
**
*Rigid body, Ref node=RBody_BaseSupport_RefNode, Position=CENTER OF MASS, Elset=RBody_BaseSupport
**
**
*NSET, NSET=P1_DoubleChamberExtrusion, ELSET=P1_DoubleChamberExtrusion
**
** ----------------  C O N T A C T    S U R F A C E    D A T A  ----------------
**
** Contact Impactor/DoubleChamberExtrusion
**
*SURFACE,TYPE=ELEMENT,NAME=Contact_1_Master
P3_Impactor,
**
*SURFACE,TYPE=NODE,NAME=Contact_Slave
P1_DoubleChamberExtrusion,
**
** Contact BaseSupport/DoubleChamberExtrusion
**
*SURFACE,TYPE=ELEMENT,NAME=Contact_2_Master
P2_BaseSupport,
**
** Self Contact
**
*Surface, Name=Contact_3_Self, Type=ELEMENT
P1_DoubleChamberExtrusion,
**
*Surface Interaction, Name=Contact_3_Prop
*Friction
  1.500000e-01,
**
** -------------  B O U N D A R Y    C O N D I T I O N    D A T A  -------------
**
** Displacement Boundary Condition (Impactor)
**
*Boundary
RBody_Impactor_RefNode, 2, 2
RBody_Impactor_RefNode, 3, 3
RBody_Impactor_RefNode, 4, 4
RBody_Impactor_RefNode, 5, 5
RBody_Impactor_RefNode, 6, 6
**
** Displacement Boundary Condition (Base Support)
**
*Boundary
RBody_BaseSupport_RefNode, 1, 6
**
**
****************** H I S T O R Y     D E F I N I T I O N ************************
**
*Step, NLGEOM=YES
**
*Dynamic, Explicit, Element by element
,0.0501
**
*DIAGNOSTICS, CONTACT INITIAL OVERCLOSURE=DETAIL
**
*NSET, NSET=Output_Refnodes
 RBody_Impactor_RefNode, RBody_BaseSupport_RefNode
**
** ------------  C O N T A C T    I N T E R A C T I O N    D A T A  -----------
**
**Contact Impactor/DoubleChamberExtrusion
*CONTACT PAIR, INTERACTION=CON01, CPSET=C_P_I, MECHANICAL CONSTRAINT=PENALTY
 Contact_Slave, Contact_1_Master
*SURFACE INTERACTION, NAME=CON01
*FRICTION
 0.15,
**
**Contact BaseSupport/DoubleChamberExtrusion
*CONTACT PAIR, INTERACTION=CON02, CPSET=C_P_S, MECHANICAL CONSTRAINT=PENALTY
 Contact_Slave, Contact_2_Master
*SURFACE INTERACTION, NAME=CON02
*FRICTION
 0.15,
**
**Self Contact
**
*Contact
**
*Contact inclusions
 Contact_3_Self, Contact_3_Self
**
*Contact property assignment
 Contact_3_Self, Contact_3_Self, Contact_3_Prop
**
*Surface property assignment, Property=THICKNESS
 Contact_3_Self,     2.500000e+00, 1.
*Contact controls assignment, nodal erosion=yes
**
******************************************************************************
*INITIAL CONDITIONS, TYPE=VELOCITY
 RBody_Impactor_RefNode, 1, 10000.0
**
**
** ----------------------  O U T P U T     R E Q U E S T  ----------------------
**
** Default output requests
**
*Output, Field, Number interval= 25
*Node output
 U,
*Element output
 S, E, STATUS, DMICRT
**
*OUTPUT, HISTORY, TIME INTERVAL=1.0e-04
*NODE OUTPUT,NSET=RBody_Impactor_RefNode
U1,
*NODE OUTPUT,NSET=RBody_BaseSupport_RefNode
RF1,
*ENERGY OUTPUT
ALLKE, ALLIE, ALLSE, ALLPD, ALLVD, ALLAE, ALLWK, ETOTAL
*End step