*HEADING
BRITTLE CRACKING MODEL, TENSION SHEAR TEST
***
*** 3D ELEMENTS
***
*NODE,NSET=BACK111
111,1.5,2.
112,2.5,2.
113,2.5,3.
114,1.5,3.
*NODE,NSET=FRONT111
115,1.5,2.,1.
116,2.5,2.,1.
117,2.5,3.,1.
118,1.5,3.,1.
*NSET,NSET=LEFT111
111,114,115,118
*NSET,NSET=RIGHT111
112,113,116,117
*NSET,NSET=BOT111
111,112,115,116
*NSET,NSET=TOP111
113,114,117,118
*ELEMENT,TYPE=C3D8R,ELSET=E111
111,111,112,113,114,115,116,117,118
*SOLID SECTION,ELSET=E111,MATE=CRACKE
*NODE,NSET=BACK211
211,1.5,4.
212,2.5,4.
213,2.5,5.
214,1.5,5.,
*NODE,NSET=FRONT211
215,1.5,4.,1.
216,2.5,4.,1.
217,2.5,5.,1.
218,1.5,5.,1.
*NSET,NSET=LEFT211
211,214,215,218
*NSET,NSET=RIGHT211
212,213,216,217
*NSET,NSET=BOT211
211,212,215,216
*NSET,NSET=TOP211
213,214,217,218
*ELEMENT,TYPE=C3D8R,ELSET=E211
211,211,212,213,214,215,216,217,218
*SOLID SECTION,ELSET=E211,MATE=CRACKT
*BOUNDARY
BOT111,2
BACK111,3
BOT211,2
BACK211,3
***
*** PLANE STRESS ELEMENTS
***
*NODE
121,3.,2.
122,4.,2.
123,4.,3.
124,3.,3.,
*NSET,NSET=LEFT121
121,124
*NSET,NSET=RIGHT121
122,123
*NSET,NSET=BOT121
121,122
*NSET,NSET=TOP121
123,124
*ELEMENT,TYPE=CPS4R,ELSET=E121
121,121,122,123,124
*SOLID SECTION,ELSET=E121,MATE=CRACKE
*NODE
221,3.,4.
222,4.,4.
223,4.,5.
224,3.,5.,
*NSET,NSET=LEFT221
221,224
*NSET,NSET=RIGHT221
222,223
*NSET,NSET=BOT221
221,222
*NSET,NSET=TOP221
223,224
*ELEMENT,TYPE=CPS4R,ELSET=E221
221,221,222,223,224
*SOLID SECTION,ELSET=E221,MATE=CRACKT
*BOUNDARY
BOT121,2
BOT221,2
***
*** PLANE STRAIN ELEMENTS
***
*NODE
131,4.5,2.
132,5.5,2.
133,5.5,3.
134,4.5,3.
*NSET,NSET=LEFT131
131,134
*NSET,NSET=RIGHT131
132,133
*NSET,NSET=BOT131
131,132
*NSET,NSET=TOP131
133,134
*ELEMENT,TYPE=CPE4R,ELSET=E131
131,131,132,133,134
*SOLID SECTION,ELSET=E131,MATE=CRACKE
*NODE
231,4.5,4.
232,5.5,4.
233,5.5,5.
234,4.5,5.,
*NSET,NSET=LEFT231
231,234
*NSET,NSET=RIGHT231
232,233
*NSET,NSET=BOT231
231,232
*NSET,NSET=TOP231
233,234
*ELEMENT,TYPE=CPE4R,ELSET=E231
231,231,232,233,234
*SOLID SECTION,ELSET=E231,MATE=CRACKT
*BOUNDARY
BOT131,2
BOT231,2
*MATERIAL,NAME=CRACKE
*DENSITY
2400., 
*ELASTIC
30.E9,.2
*BRITTLE CRACKING
3.E6,0.
0.,1.E-3
*** POWER LAW SHEAR RETENTION WITH P=2, EmaxCk=.001
*BRITTLE SHEAR,TYPE=POWER
1.E-3,2.
*MATERIAL,NAME=CRACKT
*DENSITY
2400., 
*ELASTIC
30.E9,.2
*BRITTLE CRACKING
3.E6,0.
0.,1.E-3
*** TABULAR SHEAR RETENTION FOR P=2, EmaxCk=.001
*BRITTLE SHEAR
1.,0.
.81,.1E-3
.64,.2E-3
.49,.3E-3
.36,.4E-3
.25,.5E-3
.16,.6E-3
.09,.7E-3
.04,.8E-3
.01,.9E-3
0.,1.E-3
*AMPLITUDE,NAME=RAMP1
0.,0.,0.001,1.
*AMPLITUDE,NAME=RAMP2
0.,0.,0.01,1.
*ELSET,ELSET=EALL,GEN
1,999,1
**
*STEP
*DYNAMIC,EXPLICIT
 ,0.001
*BOUNDARY,AMP=RAMP1
** 3D
111,1,,0.
114,1,,0.
115,1,,0.
118,1,,0.
112,1,,0.00011
113,1,,0.00011
116,1,,0.00011
117,1,,0.00011
211,1,,0.
214,1,,0.
215,1,,0.
218,1,,0.
212,1,,0.00011
213,1,,0.00011
216,1,,0.00011
217,1,,0.00011
** PS
121,1,,0.
124,1,,0.
122,1,,0.00011
123,1,,0.00011
221,1,,0.
224,1,,0.
222,1,,0.00011
223,1,,0.00011
** PE
131,1,,0.
134,1,,0.
132,1,,0.00011
133,1,,0.00011
231,1,,0.
234,1,,0.
232,1,,0.00011
233,1,,0.00011
**
*FILE OUTPUT,NUMBER INTERVAL=3,TIMEMARKS=YES
*EL FILE
 S, 
 LE, 
 CKE, 
 CKLE, 
 CKLS, 
 CKEMAG,
 CKSTAT,
 CRACK, 
*ELSET, ELSET=QA_TEST_ALLELEMENTS
EALL
*OUTPUT, FIELD, NUMBER INTERVAL=2, TIMEMARKS=YES
*ELEMENT OUTPUT, ELSET=QA_TEST_ALLELEMENTS
 S, 
 LE,
 CKE, 
 CKLE, 
 CKLS, 
 CKEMAG,
 CKSTAT, 
 *OUTPUT,HISTORY,TIME INTERVAL=0.015
*ENERGY OUTPUT, VAR=ALL
*END STEP
**
*STEP
*DYNAMIC,EXPLICIT
 ,0.01
*BOUNDARY,TYPE=VEL
** 3D
111,1,,0.
115,1,,0.
114,1,,0.2
118,1,,0.2
112,1,,0.1
116,1,,0.1
113,1,,0.3
117,1,,0.3
TOP111,2,,0.15
211,1,,0.
215,1,,0.
214,1,,0.2
218,1,,0.2
212,1,,0.1
216,1,,0.1
213,1,,0.3
217,1,,0.3
TOP211,2,,0.15
** PS
121,1,,0.
124,1,,0.2
122,1,,0.1
123,1,,0.3
TOP121,2,,0.15
221,1,,0.
224,1,,0.2
222,1,,0.1
223,1,,0.3
TOP221,2,,0.15
** PE
131,1,,0.
134,1,,0.2
132,1,,0.1
133,1,,0.3
TOP131,2,,0.15
231,1,,0.
234,1,,0.2
232,1,,0.1
233,1,,0.3
TOP231,2,,0.15
*END STEP