*HEADING
 MIXED MODE FAILURE OF CONCRETE: FINE MESH, CPS4R
*** LOAD-CMSD CURVE OBTAINED FROM
*** LOAD: ADD RF2 FROM NODES 10005, 10007, 10009, 10011, 10013
*** CMSD: SUBSTRACT U2 OF NODE 10025 FROM THAT OF NODE 8000
*** UNITS: NEWTON, MILLIMETER
**
** COARSE AREA OF MESH
**
*NODE,NSET=COARSE
5001,    0.,     0.
5017, 366.5,     0.
5020,  580.,     0.
5032,  916.,     0.
5601,    0.,    82.
5617, 366.5,    82.
5620,  580.,    82.
5632,  916.,    82.
7401,    0.,   306.
7417, 366.5,   306.
7420,  580.,   306.
7432,  916.,   306.
*NGEN,NSET=BOTL
5001,5017,1
*NGEN,NSET=MIDL
5601,5617,1
*NGEN,NSET=TOPL
7401,7417,1
*NGEN,NSET=BOTR
5020,5032,1
*NGEN,NSET=MIDR
5620,5632,1
*NGEN,NSET=TOPR
7420,7432,1
*NFILL,NSET=COARSE
BOTL,MIDL,6,100
MIDL,TOPL,18,100
BOTR,MIDR,6,100
MIDR,TOPR,18,100
**
** CRACK REGION
**
*NODE,NSET=NALL
10001, 366.5,     0.
10057, 580.0,     0.
11201, 366.5,    82.
11257, 580.0,    82.
14801, 366.5,   306.
14857, 580.0,   306.
8000,  458.0,     0.
9200,  458.0,    82.
*NGEN,NSET=BOT
10001,10057,1
*NGEN,NSET=MID
11201,11257,1
*NGEN,NSET=TOP
14801,14857,1
*NGEN,NSET=FINE
8000,9200,100
*NFILL,NSET=FINE
BOT,MID,12,100
MID,TOP,36,100
*NSET,NSET=BOTFACE
FINE,COARSE
*NSET,NSET=NALL
BOTFACE,
**
*NODE,NSET=DUMMY
30000, 458.0,  306.
*BOUNDARY
DUMMY,1
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY
1000000,14841
2000000,10009
*ROTARY INERTIA, ELSET=ROTARY
1.E-7,1.E-7,1.E-7
*MPC
BEAM,14837,14841
BEAM,14839,14841
BEAM,14843,14841
BEAM,14845,14841
BEAM,10005,10009
BEAM,10007,10009
BEAM,10011,10009
BEAM,10013,10009
*EQUATION
3, 
 7401,2,1.0,   30000,2,-1.0, 30000,6,458. 
3, 
14841,2,1.0,   30000,2,-1.0, 30000,6,-61. 
2, 
7401,1,1.0,30000,1,-1.0
2, 
14841,1,1.0,30000,1,-1.0
*ELEMENT,TYPE=CPS4R,ELSET=COARSE
5001,  5001,5003,5203,5201
5020,  5020,5022,5222,5220
*ELGEN,ELSET=COARSE
5001,  8,2,2, 12,200,200
5020,  6,2,2, 12,200,200
**
*ELEMENT,TYPE=CPS4R,ELSET=FINE
10001, 10001,10003,10203,10201
10025,  8000,10027,10227, 8200
11125,  9000,11027,11227,11225
10027, 10027,10029,10229,10227
11225, 11225,11227,11427,11425
*ELGEN,ELSET=FINE
10001, 12,2,2,  24,200,200
10025, 5,200,200
10027, 15,2,2,  24,200,200
11225, 18,200,200
*ELSET,ELSET=EALL
COARSE,FINE
**
*NSET,NSET=NOUT
10009,20009, 5032,15032, 7401,17401,
14841,24841,30000,40000,10025,8000,20025,18000
10007,10011,10005,10013
*NSET,NSET=NOUT,GEN
10025,11225,200
20025,21225,200
 8000, 9200,200
18000,19200,200
*ELSET,ELSET=EOUT,GEN
10023,14423,200
10025,14425,200
**
*ELSET,ELSET=LCOARSE,GEN
5015,7415,200
*ELSET,ELSET=RCOARSE,GEN
5020,7420,200
*ELSET,ELSET=LFINE,GEN
10001,14601,200
*ELSET,ELSET=RFINE,GEN
10055,14655,200
**
*SOLID SECTION,ELSET=EALL,MATERIAL=CONC
156.,
**
** G_t  = .055 N/mm
** h    = 7.5 mm
** f_ct = 2.8 N/mm^2
** density = 2.4E-6 kg/mm^3
**
*MATERIAL,NAME=CONC
*ELASTIC
24800, .18
*BRITTLE CRACKING
2.8, 0.0
0.0, 0.005239
*BRITTLE SHEAR
1.0,  0.
.01,  2.910E-4
0.001,  52.39E-4
*DENSITY
2.4E-6,
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY
40000,30000
*ROTARY INERTIA,ELSET=ROTARY
.1,.1,.1
*AMPLITUDE,NAME=RAMP
0.,0.,  0.38,1.0
*SURFACE,TYPE=ELEMENT,NAME=RFINE
RFINE,S2
*SURFACE,TYPE=ELEMENT,NAME=LCOARSE
LCOARSE,S2
*SURFACE,TYPE=ELEMENT,NAME=LFINE
LFINE,S4
*SURFACE,TYPE=ELEMENT,NAME=RCOARSE
RCOARSE,S4
*TIE, NAME=TIE
LFINE, LCOARSE
RFINE, RCOARSE
**
*STEP
*DYNAMIC,EXPLICIT
,0.38
*BULK VISCOSITY
0.02, 0.4
*BOUNDARY
10009, 2
 5032, 2
*BOUNDARY,TYPE=VEL,AMP=RAMP
30000, 2,2,  -0.75
**
*OUTPUT,FIELD,VAR=PRESELECT,NUMBER INTERVAL=5
**
** Output requests for qa testing
**
*NSET,NSET=QA_TEST
NOUT,
*OUTPUT,FIELD,NUMBER INTERVAL=2,TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST
 U,
*OUTPUT,HISTORY,VAR=PRESELECT,TIME INTERVAL=0.095
*FILE OUTPUT,NUMBER INTERVAL=2,TIMEMARKS=YES
*NODE FILE, NSET=QA_TEST
 U,
**
*END STEP