*HEADING
BRITTLE  MODEL, CRACK COMBINATIONS, UNLOADING/RELOADING
***
*** AXISYMMETRIC ELEMENTS
***
*NODE
1,
2,1.,
3,1.,1.,
4,0.,1.,
*NSET,NSET=LEFT1
1,4
*NSET,NSET=RIGHT1
2,3
*NSET,NSET=BOT1
1,2
*NSET,NSET=TOP1
3,4
*ELEMENT,TYPE=CAX4R,ELSET=E1
1,1,2,3,4
*SOLID SECTION,ELSET=E1,MATE=CRACKE
*NODE
101,,2.
102,1.,2.
103,1.,3.
104,0.,3.,
*NSET,NSET=LEFT101
101,104
*NSET,NSET=RIGHT101
102,103
*NSET,NSET=BOT101
101,102
*NSET,NSET=TOP101
103,104
*ELEMENT,TYPE=CAX4R,ELSET=E101
101,101,102,103,104
*SOLID SECTION,ELSET=E101,MATE=CRACKE
*NODE
201,,4.
202,1.,4.
203,1.,5.
204,0.,5.,
*NSET,NSET=LEFT201
201,204
*NSET,NSET=RIGHT201
202,203
*NSET,NSET=BOT201
201,202
*NSET,NSET=TOP201
203,204
*ELEMENT,TYPE=CAX4R,ELSET=E201
201,201,202,203,204
*SOLID SECTION,ELSET=E201,MATE=CRACKE
*NODE
301,,6.
302,1.,6.
303,1.,7.
304,0.,7.
*NSET,NSET=LEFT301
301,304
*NSET,NSET=RIGHT301
302,303
*NSET,NSET=BOT301
301,302
*NSET,NSET=TOP301
303,304
*ELEMENT,TYPE=CAX4R,ELSET=E301
301,301,302,303,304
*SOLID SECTION,ELSET=E301,MATE=CRACKD
*NODE
401,,8.
402,1.,8.
403,1.,9.
404,0.,9.,
*NSET,NSET=LEFT401
401,404
*NSET,NSET=RIGHT401
402,403
*NSET,NSET=BOT401
401,402
*NSET,NSET=TOP401
403,404
*ELEMENT,TYPE=CAX4R,ELSET=E401
401,401,402,403,404
*SOLID SECTION,ELSET=E401,MATE=CRACKG
*BOUNDARY
LEFT1,1
BOT1,2
LEFT101,1
BOT101,2
LEFT201,1
BOT201,2
LEFT301,1
BOT301,2
LEFT401,1
BOT401,2
***
*** 3D ELEMENTS
***
*NODE,NSET=BACK11
11,1.5
12,2.5,
13,2.5,1.,
14,1.5,1.,
*NODE,NSET=FRONT11
15,1.5,,1.
16,2.5,,1.
17,2.5,1.,1.
18,1.5,1.,1.
*NSET,NSET=LEFT11
11,14,15,18
*NSET,NSET=RIGHT11
12,13,16,17
*NSET,NSET=BOT11
11,12,15,16
*NSET,NSET=TOP11
13,14,17,18
*ELEMENT,TYPE=C3D8R,ELSET=E11
11,11,12,13,14,15,16,17,18
*SOLID SECTION,ELSET=E11,MATE=CRACKE
*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=CRACKE
*NODE,NSET=BACK311
311,1.5,6.
312,2.5,6.
313,2.5,7.
314,1.5,7.
*NODE,NSET=FRONT311
315,1.5,6.,1.
316,2.5,6.,1.
317,2.5,7.,1.
318,1.5,7.,1.
*NSET,NSET=LEFT311
311,314,315,318
*NSET,NSET=RIGHT311
312,313,316,317
*NSET,NSET=BOT311
311,312,315,316
*NSET,NSET=TOP311
313,314,317,318
*ELEMENT,TYPE=C3D8R,ELSET=E311
311,311,312,313,314,315,316,317,318
*SOLID SECTION,ELSET=E311,MATE=CRACKD
*NODE,NSET=BACK411
411,1.5,8.
412,2.5,8.
413,2.5,9.
414,1.5,9.,
*NODE,NSET=FRONT411
415,1.5,8.,1.
416,2.5,8.,1.
417,2.5,9.,1.
418,1.5,9.,1.
*NSET,NSET=LEFT411
411,414,415,418
*NSET,NSET=RIGHT411
412,413,416,417
*NSET,NSET=BOT411
411,412,415,416
*NSET,NSET=TOP411
413,414,417,418
*ELEMENT,TYPE=C3D8R,ELSET=E411
411,411,412,413,414,415,416,417,418
*SOLID SECTION,ELSET=E411,MATE=CRACKG
*BOUNDARY
LEFT11,1
BOT11,2
BACK11,3
LEFT111,1
BOT111,2
BACK111,3
LEFT211,1
BOT211,2
BACK211,3
LEFT311,1
BOT311,2
BACK311,3
LEFT411,1
BOT411,2
BACK411,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=CRACKE
*NODE
321,3.,6.
322,4.,6.
323,4.,7.
324,3.,7.
*NSET,NSET=LEFT321
321,324
*NSET,NSET=RIGHT321
322,323
*NSET,NSET=BOT321
321,322
*NSET,NSET=TOP321
323,324
*ELEMENT,TYPE=CPS4R,ELSET=E321
321,321,322,323,324
*SOLID SECTION,ELSET=E321,MATE=CRACKD
*NODE
421,3.,8.
422,4.,8.
423,4.,9.
424,3.,9.,
*NSET,NSET=LEFT421
421,424
*NSET,NSET=RIGHT421
422,423
*NSET,NSET=BOT421
421,422
*NSET,NSET=TOP421
423,424
*ELEMENT,TYPE=CPS4R,ELSET=E421
421,421,422,423,424
*SOLID SECTION,ELSET=E421,MATE=CRACKG
*BOUNDARY
LEFT121,1
BOT121,2
LEFT221,1
BOT221,2
LEFT321,1
BOT321,2
LEFT421,1
BOT421,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=CRACKE
*NODE
331,4.5,6.
332,5.5,6.
333,5.5,7.
334,4.5,7.
*NSET,NSET=LEFT331
331,334
*NSET,NSET=RIGHT331
332,333
*NSET,NSET=BOT331
331,332
*NSET,NSET=TOP331
333,334
*ELEMENT,TYPE=CPE4R,ELSET=E331
331,331,332,333,334
*SOLID SECTION,ELSET=E331,MATE=CRACKD
*NODE
431,4.5,8.
432,5.5,8.
433,5.5,9.
434,4.5,9.,
*NSET,NSET=LEFT431
431,434
*NSET,NSET=RIGHT431
432,433
*NSET,NSET=BOT431
431,432
*NSET,NSET=TOP431
433,434
*ELEMENT,TYPE=CPE4R,ELSET=E431
431,431,432,433,434
*SOLID SECTION,ELSET=E431,MATE=CRACKG
*BOUNDARY
LEFT131,1
BOT131,2
LEFT231,1
BOT231,2
LEFT331,1
BOT331,2
LEFT431,1
BOT431,2
***
*** BEAM ELEMENTS
***
*NODE,NSET=BEAMN
241,6.,4.
244,6.,5.,
*NSET,NSET=BOT241
241, 
*NSET,NSET=TOP241
244, 
*ELEMENT,TYPE=B21,ELSET=E241
241,241,244
*BEAM SECTION,SECTION=RECT,ELSET=E241,MATE=CRACKE
1.,1.
*NODE,NSET=BEAMN
341,6.,6.
344,6.,7.
*NSET,NSET=BOT341
341, 
*NSET,NSET=TOP341
344, 
*ELEMENT,TYPE=B21,ELSET=E341
341,341,344
*BEAM SECTION,SECTION=RECT,ELSET=E341,MATE=CRACKD
1.,1.
*NODE,NSET=BEAMN
441,6.,8.
444,6.,9.,
*NSET,NSET=BOT441
441, 
*NSET,NSET=TOP441
444, 
*ELEMENT,TYPE=B21,ELSET=E441
441,441,444
*BEAM SECTION,SECTION=RECT,ELSET=E441,MATE=CRACKG
1.,1.
*BOUNDARY
BEAMN,1
BEAMN,6
BOT241,2
BOT341,2
BOT441,2
*MATERIAL,NAME=CRACKE
*DENSITY
2400., 
*ELASTIC
30.E9,.2
*BRITTLE CRACKING
3.E6,0.
1.E6,.5E-3
0.,1.E-3
*BRITTLE SHEAR
1.,0.
.25,.5E-3
0.,1.E-3
*MATERIAL,NAME=CRACKD
*DENSITY
2400., 
*ELASTIC
30.E9,.2
*BRITTLE CRACKING,TYPE=DISP
3.E6,0.
1.E6,.5E-3
0.,1.E-3
*BRITTLE SHEAR
1.,0.
.25,.5E-3
0.,1.E-3
*MATERIAL,NAME=CRACKG
*DENSITY
2400., 
*ELASTIC
30.E9,.2
*BRITTLE CRACKING,TYPE=GFI
3.E6,1500.
*BRITTLE SHEAR
1.,0.
.25,.5E-3
0.,1.E-3
*AMPLITUDE,NAME=UNRELOAD
0.,0.,0.0025,0.0005,.0055,-0.0001,.015,.0018
*ELSET,ELSET=EALL,GEN
1,999,1
**
*STEP
*DYNAMIC,EXPLICIT
 ,0.015
*BOUNDARY,AMP=UNRELOAD
** AXI
RIGHT1,1,,1.
TOP1,2,,.75
RIGHT101,1,,1.
TOP201,2,,1.
TOP301,2,,1.
TOP401,2,,1.
** 3D
RIGHT11,1,,1.
TOP11,2,,.75
FRONT11,3,,1.
RIGHT111,1,,1.
TOP111,2,,.75
FRONT111,3,,.0
TOP211,2,,1.
TOP311,2,,1.
TOP411,2,,1.
** PS
RIGHT121,1,,1.
TOP121,2,,.75
TOP221,2,,1.
TOP321,2,,1.
TOP421,2,,1.
** PE
RIGHT131,1,,1.
TOP131,2,,.75
TOP231,2,,1.
TOP331,2,,1.
TOP431,2,,1.
** BEAM
TOP241,2,,1.
TOP341,2,,1.
TOP441,2,,1.
**
*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