*HEADING
FeFp: cgax4r, uniaxial tension
Yeoh, compressible, mullins effect
Mullins test data for input
*preprint,model=yes,history=yes
*NODE,NSET=nall
1,
2,1.
3,1.,1.,
4,0.,1.,
11,2.,0.
12,3.,0.
13,3.,1.
14,2.,1.
*NSET,NSET=left
1,4,11,14
*NSET,NSET=top
3,4,13,14
*NSET,NSET=bot
1,2,11,12
*nset,nset=nout
3,4,13,14
*ELEMENT,TYPE=cax4r,ELSET=ONE
1,1,2,3,4
*ELEMENT,TYPE=cax4r,ELSET=two
2,11,12,13,14
*SOLID SECTION,ELSET=two,MATERIAL=RUB2
*MATERIAL,NAME=RUB2
*DENSITY
 1.e-3
*HYPERELASTIC,YEOH
191.322,4.9415,0.5346,5.E-06
*MULLINS EFFECT,TEST DATA INPUT
*INCLUDE,INPUT=mullins_calibrate_testdata.inp
*SOLID SECTION,ELSET=one,MATERIAL=RUB1
*MATERIAL,NAME=RUB1
*DENSITY
 1.e-3
*HYPERELASTIC,YEOH
191.322,4.9415,0.5346,5.E-06
*MULLINS EFFECT,TEST DATA INPUT
*INCLUDE,INPUT=mullins_calibrate_testdata.inp
*plastic
50.0,  0.0
200.0, 0.02
1600.0, 0.1
6000.0,0.25
65000.0,1.75
*amplitude,name=bc,definition=smooth step,time=total
0.0,0.0, 0.1,0.5
*elset,elset=eall
one,two
*STEP
UNIAXIAL TENSION
*dynamic, explicit
 ,0.1
*BOUNDARY
left,1,
bot,2,
*boundary,amp=bc
top,2,2,1.0
*output,history,freq=1
*element output, elset=eall
s,ne,peeq,eldmd
*node output,nset=nout
u2,rf2
*energy output
alldmd
*OUTPUT,FIELD
*ELEMENT OUTPUT
S
*NODE OUTPUT
U
*END STEP