*HEADING
 Quasi-static Three Point Bending 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=S4R, Elset=P1_DoubleChamberExtrusion
**
*Shell section, Elset=P1_DoubleChamberExtrusion, Material=AluminumAlloy
  2.5000e+00,     5
**
*Element, Type=S4R, Elset=P2_Punch
**
*Shell section, Elset=P2_Punch, Material=ElasticMat
 1.0000e+00,     3
**
*Element, Type=S4R, Elset=P3_LeftSupport
**
*Shell section, Elset=P3_LeftSupport, Material=ElasticMat
  1.0000e+00,     3
**
*Element, Type=S4R, Elset=P4_RightSupport
**
*Shell section, Elset=P4_RightSupport, Material=ElasticMat
  1.0000e+00,     3
**
** -----------------------  M A T E R I A L     D A T A  -----------------------
**
** Punch and Support
**
*Material, Name=ElasticMat
*Density
  7.820000e-09
*Elastic
  2.070000e+05,   3.000000e-01
**
** 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
*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
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
**
**
*NSET, NSET=P1_DoubleChamberExtrusion, ELSET=P1_DoubleChamberExtrusion
**
** ---------------------  R I G I D    B O D Y    D A T A  ---------------------
**
** Punch
**
*Elset, Elset=RBody_Punch
P2_Punch,
**
*Nset, Nset=RBody_Punch_RefNode
**
*Rigid body, Ref node=RBody_Punch_RefNode, Position=CENTER OF MASS, Elset=RBody_Punch
**
** Left Support
**
*Elset, Elset=RBody_LeftSupport
P3_LeftSupport,
**
*Nset, Nset=RBody_LeftSupport_RefNode
**
*Rigid body, Ref node=RBody_LeftSupport_RefNode, Position=CENTER OF MASS, Elset=RBody_LeftSupport
**
** Right Support
**
*Elset, Elset=RBody_RightSupport
P4_RightSupport,
**
*Nset, Nset=RBody_RightSupport_RefNode
**
*Rigid body, Ref node=RBody_RightSupport_RefNode, Position=CENTER OF MASS, Elset=RBody_RightSupport
**
** ----------------  C O N T A C T    S U R F A C E    D A T A  ----------------
**
** Contact Surfaces - Punch/DoubleChamberExtrusion
**
*Surface, Name=Contact_1_Master, Type=ELEMENT
 P2_Punch,
**
*Surface, Name=Contact_1_Slave, Type=NODE
 P1_DoubleChamberExtrusion,
**
** Contact Surfaces - LeftSupport/DoubleChamberExtrusion
**
*Surface, Name=Contact_2_Master, Type=ELEMENT
 P3_LeftSupport,
**
*Surface, Name=Contact_2_Slave, Type=ELEMENT
 P1_DoubleChamberExtrusion,
**
*Surface Interaction, Name=Contact_2_Prop
*Friction
 5.000000e-02,
**
** Contact Surfaces - RightSupport/DoubleChamberExtrusion
**
*Surface, Name=Contact_3_Master, Type=ELEMENT
 P4_RightSupport,
**
*Surface, Name=Contact_3_Slave, Type=ELEMENT
 P1_DoubleChamberExtrusion,
**
*Surface Interaction, Name=Contact_3_Prop
*Friction
 5.000000e-02,
**
** Self Contact
**
*Surface, Name=Contact_4_Self, Type=ELEMENT
 P1_DoubleChamberExtrusion,
**
*Surface Interaction, Name=Contact_4_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 (Punch)
**
*Boundary
 RBody_Punch_RefNode, 1, 1
 RBody_Punch_RefNode, 3, 3
 RBody_Punch_RefNode, 4, 4
 RBody_Punch_RefNode, 5, 5
 RBody_Punch_RefNode, 6, 6
*Nset, nset=QA_TEST
RBody_Punch_RefNode,
**
** Displacement Boundary Condition (Left Support)
**
*Boundary
 RBody_LeftSupport_RefNode, 1, 6
**
** Displacement Boundary Condition (Right Support)
**
*Boundary
 RBody_RightSupport_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
 , 7.200000e-02
**
*DIAGNOSTICS, CONTACT INITIAL OVERCLOSURE=DETAIL
**
** ------------  C O N T A C T    I N T E R A C T I O N    D A T A  ------------
**
*CONTACT PAIR, INTERACTION=CON01, CPSET=C_P_S, MECHANICAL CONSTRAINT=PENALTY
 Contact_1_Slave, Contact_1_Master
*SURFACE INTERACTION, NAME=CON01
*FRICTION
 0.05,
**
*Contact
**
*Contact inclusions
 Contact_2_Slave, Contact_2_Master
 Contact_3_Slave, Contact_3_Master
 Contact_4_Self, Contact_4_Self
**
*Contact property assignment
 Contact_2_Slave, Contact_2_Master, Contact_2_Prop
 Contact_3_Slave, Contact_3_Master, Contact_3_Prop
 Contact_4_Self, Contact_4_Self, Contact_4_Prop
**
*Surface property assignment, Property=THICKNESS
 Contact_2_Slave,    1.250000e+00, 1.
 Contact_2_Master,   1.250000e+00, 1.
 Contact_3_Slave,    1.250000e+00, 1.
 Contact_3_Master,   1.250000e+00, 1.
 Contact_4_Self,     2.500000e+00, 1.
**
*Contact formulation, Type=PURE MASTER-SLAVE
 Contact_2_Slave, Contact_2_Master, SLAVE
 Contact_3_Slave, Contact_3_Master, SLAVE
**
**
******************************************************************************
*AMPLITUDE, NAME=velocity
 .0, 0.0, 1.25000e-04, 79.972664, 2.50000e-04, 159.822006, 3.75000e-04, 239.424896,
 5.00e-04, 318.6586, 6.25e-04, 397.401062, 7.50e-04, 475.53067, 8.75e-04, 552.927002,
 0.001, 629.470764, 0.001125, 705.044128, 0.00125, 779.530212, 0.001375, 852.814331,
 0.0015, 924.783447, 0.001625, 995.326538, 0.00175, 1064.335083, 0.001875, 1131.702393,
 0.002, 1197.324707, 0.002125, 1261.100952, 0.00225, 1322.932495, 0.002375, 1382.724121,
 0.0025, 1440.383911, 0.002625, 1495.82251, 0.00275, 1548.954834, 0.002875, 1599.698975,
 0.003, 1647.976074, 0.003125, 1693.712402, 0.00325, 1736.836914, 0.003375, 1777.283691,
 0.0035, 1814.989624, 0.003625, 1849.897095, 0.00375, 1881.952148, 0.003875, 1911.105469,
*BOUNDARY, TYPE=VELOCITY, AMPLITUDE=velocity
 RBody_Punch_RefNode, 2, 2, -1
** ----------------------  O U T P U T     R E Q U E S T  ----------------------
**
** Default output requests
**
*Output, Field, TIME INTERVAL=0.001
*Node output
 U,
*Element output
 S, E, STATUS, DMICRT
*Output, Field, NUMBER INTERVAL = 4
*Node output, nset=QA_TEST
 U,
**
*OUTPUT, HISTORY, TIME INTERVAL=1.0e-04
*NODE OUTPUT,NSET=RBody_Punch_RefNode
 U2, RF2
*Energy output
 ALLKE, ALLIE, ALLSE, ALLPD, ALLVD, ALLAE, ALLWK, ETOTAL
**
*End step