*HEADING
FeFp: cgax4rh, uniaxial tension and perturbation
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=Cgax4rh,ELSET=ONE
1,1,2,3,4
*ELEMENT,TYPE=Cgax4rh,ELSET=two
2,11,12,13,14
*SOLID SECTION,ELSET=two,MATERIAL=RUB2
*MATERIAL,NAME=RUB2
*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
*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
*DENSITY
0.43
*amplitude,name=bc, time=total
0.0,0.0, 10.0,0.5
*elset,elset=eall
one,two
*STEP,NLGEOM,INC=2000
UNIAXIAL TENSION
*STATIC
1.0,10.,,1.0
*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,FREQ=100
*ELEMENT OUTPUT
S
*OUTPUT,FIELD,FREQ=100
*NODE OUTPUT
U
*END STEP
*step, perturbation
*static
*boundary
top, 2,2, 0.01
*output, field
*element output,elset=eall
s
*node output
rf
*el print
s,
ne,
peeq
*node print
u,
rf
*end step