*HEADING
 COMPOSITE PLATE WITH EXPLICIT DAMAGE AND
 CPS4R FAILURE MODEL
** UNITS: IN, LBS
** -------------------------
** MODEL DEFINITION
** -------------------------
*PREPRINT,MODEL=YES
*NODE,INPUT=damagefailcomplate_node.inp,NSET=NALL
*ELEMENT,INPUT=damagefailcomplate_element.inp,
TYPE=CPS4R,ELSET=L1
***ELEMENT, TYPE=MASS, ELSET=REF
**999999,2000
*ELCOPY,ELEMENT SHIFT=1000,OLD SET=L1,NEW SET=L2,
SHIFT NODES=0
*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
*DENSITY
 1.0e-4
*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, 
1, EM, "Matrix failure index"
2, EFS, "Fiber-matrix shearing failure index"
3, DAMAGE, "Shear damage parameter"
*USER DEFINED FIELD, PROPERTIES=6
** yt  ,   yc  ,   sc  ,  xc    ,  g12  ,  alpha
14.86d3, 36.7d3, 15.5d3, 392.7d3, 1.01d6, 0.8d-14
**
** LOCAL ORIENTATIONS: P45 AT +45, N45 AT -45
**
*ORIENTATION,NAME=P45
0.,1.,0.,-1.,0.,0.
3,+45.0
*ORIENTATION,NAME=N45
0.,1.,0.,-1.,0.,0.
3,-45.0
*SOLID SECTION,ELSET=L1,ORIENTATION=P45,
MATERIAL=T300
0.06750, 
*SOLID SECTION,ELSET=L2,ORIENTATION=N45,
MATERIAL=T300
0.06750, 
***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
*NSET,NSET=OUTPUT
1000,2000,204,197
*ELSET,ELSET=EOUT,GENERATE
133,144
11133,11144
*NSET,NSET=NLOAD
1000
*NSET,NSET=QA_TEST
NLOAD
** -------------------
** ANALYSIS HISTORY
** -------------------
*amplitude,name=amp
 0.0,0.0,1.0e-3,1.0
*FILTER, NAME=SMOOTH, TYPE=BUTTERWORTH
 9.e3, 4
*STEP,nlgeom=yes
*DYNAMIC,EXPLICIT
 ,1.0e-3
*BOUNDARY
XSYMMTRY,XSYMM
YSYMMTRY,YSYMM
*BOUNDARY, amplitude=amp
1000,2,,-0.027
*OUTPUT, FIELD
*NODE OUTPUT
U,RF
*ELEMENT OUTPUT
S,LE
SDV,FV
*OUTPUT,HISTORY
*NODE OUTPUT, NSET=NLOAD
RF2, U2
******************************
** OUTPUT FOR QA PURPOSES
******************************
*OUTPUT, HISTORY, TIME INTERVAL=1.0e-5, FILTER=SMOOTH
*NODE OUTPUT, NSET=QA_TEST
 RF2, U2
******************************
*END STEP