*HEADING
 MIXED MODE FAILURE OF CONCRETE: COARSE MESH, CPS4R
*** LOAD-CMSD CURVE OBTAINED FROM
*** LOAD: ADD RF2 FROM NODES 10005, 10009, 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,14845,14841
BEAM,10005,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,5005,5405,5401
5020,  5020,5024,5424,5420
*ELGEN,ELSET=COARSE
5001,  4,4,4, 6,400,400
5020,  3,4,4, 6,400,400
**
*ELEMENT,TYPE=CPS4R,ELSET=FINE
10001, 10001,10005,10405,10401
10025,  8000,10029,10429, 8400
10825,  8800,10829,11229,11225 
10029, 10029,10033,10433,10429
11225, 11225,11229,11629,11625
*ELGEN,ELSET=FINE
10001, 6,4,4,  12,400,400
10025, 2,400,500
10029, 7,4,4,  12,400,400
11225, 9,400,400
*ELSET,ELSET=EALL
COARSE,FINE
**
*NSET,NSET=NOUT
10009,20009,5032,15032, 7401,17401,
14841,24841,30000,40000,10025,8000,20025,18000,
10005,10013
*NSET,NSET=NOUT,GEN
10025,11225,400
20025,21225,400
 8000, 9200,400
18000,19200,400
*ELSET,ELSET=EOUT,GEN
10021,14421,400
10025,14425,400
**
*ELSET,ELSET=LCOARSE,GEN
5013,7413,400
*ELSET,ELSET=RCOARSE,GEN
5020,7420,400
*ELSET,ELSET=LFINE,GEN
10001,14401,400
*ELSET,ELSET=RFINE,GEN
10053,14453,400
**
*SOLID SECTION,ELSET=EALL,MATERIAL=CONC
156.,
**
** G_t  = .055 N/mm
** h    = 15 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.002619
*BRITTLE SHEAR
1.0,  0.
.01,  1.455E-4
0.0,  26.19E-4
*DENSITY
2.4E-6, 
*DAMPING, BETA=1.E-6
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY
30000,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,FILTER=ANTIALIASING
*FILE OUTPUT,NUMBER INTERVAL=2,TIMEMARKS=YES
*NODE FILE, NSET=QA_TEST
 U,
**
*END STEP