*HEADING
 COMPOSITE PLATE WITH EXPLICIT DAMAGE AND FAILURE MODEL
 1/4 SYMMETRY MODEL
** UNITS: IN, LBS
** -------------------------
** MODEL DEFINITION
** -------------------------
*NODE,INPUT=damagefailcomplate_node.inp,NSET=NALL
*ELEMENT,INPUT=damagefailcomplate_element.inp,TYPE=S4R,ELSET=PLATE
*ELEMENT, TYPE=MASS, ELSET=REF
999999, 2000
*NSET,NSET=XSYMMTRY
       1,       2,       3,       4,       5,       6,       7,       8,
       9,      10,      11,      12,      13,     190,     197,     204,
     211,     218,     225,     232,     239,     246,     253,     260,
     267,    1000
*NSET,NSET=YSYMMTRY
     170,     171,     172,     173,     174,     175,     176,     177,
     178,     179,     180,     181,     182,     287,     294,     301,
     308,     315,      322
*NSET,NSET=TOP
     267,     268,     269,     270,     271,     272,     273,     408,
     409,     410,     411,     412,      413
**
** MATERIAL: NONLINEAR SHEAR WITH BUILT-IN EXPLICIT FAILURE
**
** FV1: MATRIX COMPRESSIVE/TENSILE FAILURE
** FV2: FIBER-MATRIX SHEAR FAILURE
** FV3: SHEAR NONLINEARITY (DAMAGE) PRIOR TO FAILURE
** TOTAL OF 2^3 = 8 STATES
**
*MATERIAL,NAME=T300
*ELASTIC,TYPE=LAMINA,DEPENDENCIES=3
22.7E6,1.88E6,0.23,1.01E6,1.01E6,1.01E6,0.,0,
0,0
22.7E6,1.00E0,0.00,1.01E6,1.01E6,1.01E6,0.,1,
0,0
22.7E6,1.88E6,0.00,1.00E0,1.01E6,1.01E6,0.,0,
1,0
22.7E6,1.00E0,0.00,1.00E0,1.01E6,1.01E6,0.,1,
1,0
22.7E6,1.88E6,0.23,1.00E0,1.01E6,1.01E6,0.,0,
0,1
22.7E6,1.00E0,0.00,1.00E0,1.01E6,1.01E6,0.,1,
0,1
22.7E6,1.88E6,0.00,1.00E0,1.01E6,1.01E6,0.,0,
1,1
22.7E6,1.00E0,0.00,1.00E0,1.01E6,1.01E6,0.,1,
1,1
*DEPVAR
3, 
*USER DEFINED FIELD
**
** LOCAL ORIENTATION: 90 DEG. CW. WRT GLOBAL SYSTEM
**
*ORIENTATION,NAME=LOCAL
0.,1.,0.,-1.,0.,0.
3, 
*SHELL SECTION,COMPOSITE,ELSET=PLATE,ORIENTATION=LOCAL
0.06750,1,T300,+45.0
0.06750,1,T300,-45.0
*MASS, ELSET=REF
1.0E-5
**
** CONSTRAIN TOP FOR THE PURPOSES OF LOAD APPLICATION
**
*EQUATION
2, 
TOP,2,1.0,1000,2,-1.0
**
** CONSTRAIN OUTPUT NODE BETWEEN SURROUNDING NODES 
*MPC
LINEAR,2000,197,204
*ELSET,ELSET=EOUT
144, 
*NSET,NSET=OUTPUT
1000,2000,204,197
** -------------------
** ANALYSIS HISTORY
** -------------------
*STEP,INC=200,NLGEOM
*STATIC,DIRECT
0.05,1.0
*BOUNDARY
NALL,3,5,0.
XSYMMTRY,XSYMM
YSYMMTRY,YSYMM
1000,2,,-0.027
*RESTART,WRITE,OVERLAY
*NODE PRINT,NSET=OUTPUT
U2,RF2
*NODE FILE,NSET=OUTPUT
U,RF
*EL FILE,ELSET=EOUT
S,E
SDV,FV
*EL PRINT,FREQUENCY=5,ELSET=EOUT
S,E
*END STEP
** ---------------------
** USER-DEFINED FIELD
** ---------------------