*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
*NGEN,NSET=BOTL
*NGEN,NSET=MIDL
*NGEN,NSET=TOPL
*NGEN,NSET=BOTR
*NGEN,NSET=MIDR
*NGEN,NSET=TOPR
*NFILL,NSET=COARSE
BOTL,MIDL,6,100
MIDL,TOPL,18,100
BOTR,MIDR,6,100
MIDR,TOPR,18,100
**
** CRACK REGION
**
*NODE,NSET=NALL
*NGEN,NSET=BOT
*NGEN,NSET=MID
*NGEN,NSET=TOP
*NGEN,NSET=FINE
*NFILL,NSET=FINE
BOT,MID,12,100
MID,TOP,36,100
*NSET,NSET=BOTFACE
FINE,COARSE
*NSET,NSET=NALL
BOTFACE,
**
*NODE,NSET=DUMMY
*BOUNDARY
DUMMY,1
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY
*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
*ELEMENT,TYPE=CPS4R,ELSET=COARSE
*ELGEN,ELSET=COARSE
**
*ELEMENT,TYPE=CPS4R,ELSET=FINE
*ELGEN,ELSET=FINE
*ELSET,ELSET=EALL
COARSE,FINE
**
*NSET,NSET=NOUT
*NSET,NSET=NOUT,GEN
*ELSET,ELSET=EOUT,GEN
**
*ELSET,ELSET=LCOARSE,GEN
*ELSET,ELSET=RCOARSE,GEN
*ELSET,ELSET=LFINE,GEN
*ELSET,ELSET=RFINE,GEN
**
*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
*BRITTLE SHEAR
*DENSITY
2.4E-6,
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY
*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.34
*BULK VISCOSITY
*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