*heading
 ANISOTROPIC HYPERELASTIC, kaliske with VISCOELASTICITY and MULLINS EFFECT
 Uniaxial test, displacement control
 Base problem for Explicit->Explicit Import
*preprint, model=yes
*node,nset=all
 11,
 12,1.
 13,1.,1.,
 14,0.,1.,
 15,0.,0.,1.
 16,1.,0.,1.
 17,1.,1.,1.
 18,0.,1.,1.
*nset,nset=face11
 11,12,13,14
*nset,nset=face12
 15,16,17,18
*nset,nset=face13
 11,12,15,16
*nset,nset=face14
 13,16,17
*nset,nset=face15
 13,14,17,18
*nset,nset=face16
 14,11,18,15
*nset, nset=ndisp1
 12
*element,type=c3d8r,elset=c3d
 1,11,12,13,14,15,16,17,18
**
*node,nset=all
 21,0.,2.,
 22,1.,2.,
 23,1.,3.,
 24,0.,3.,
*nset,nset=face21
 21,22
*nset,nset=face22
 23
*nset,nset=face23
 23,24
*nset,nset=face24
 21,24
*nset, nset=ndisp2
 22
*element,type=cpe4r,elset=cpe
 2,21,22,23,24
**
*node,nset=all
 31,0.,4.,
 32,1.,4.,
 33,1.,5.,
 34,0.,5.,
*nset,nset=face31
 31,32
*nset,nset=face32
 33
*nset,nset=face33
 33,34
*nset,nset=face34
 31,34
*nset, nset=ndisp3
 32
*element,type=cps4r,elset=cps
 3,31,32,33,34
**
*node,nset=all
 41,0.,6.,
 42,1.,6.,
 43,1.,7.,
 44,0.,7.,
*nset,nset=face41
 41,42
*nset,nset=face42
 43
*nset,nset=face43
 43,44
*nset,nset=face44
 41,44
*nset, nset=ndisp4
 42
*element,type=m3d4r,elset=m3d
 4,41,42,43,44
**
*node,nset=all
 51,0.,8.,
 52,1.,8.,
 53,1.,9.,
 54,0.,9.,
*nset,nset=face51
 51,52
*nset,nset=face52
 53
*nset,nset=face53
 53,54
*nset,nset=face54
 51,54
*nset, nset=ndisp5
 52
*nset, nset=QA_TEST
17,22,32,42,52
*element,type=s4r,elset=s4r
 5,51,52,53,54
**
*elset, elset=solid
 c3d, cpe, cps
*elset, elset=eall
 solid, m3d, s4r
*elset, elset = QA_TEST
eall,
**
*orientation,name=ori-1,local directions=2
 1.0,0.0,0.0,0.0,1.0,0.0
 3,0.0
 0.643055,0.76582,0.0
 0.643055,-0.76582,0.0
**
*solid section,elset=c3d,material=kaliske,orientation=ori-1
 1.0,
*solid section,elset=cpe,material=kaliske,orientation=ori-1
 1.0,
*solid section,elset=cps,material=kaliske,orientation=ori-1
 1.0,
*membrane section,elset=m3d,material=kaliske,orientation=ori-1
 1.0,
*shell section,elset=s4r,material=kaliske,orientation=ori-1
 1.0,
**
*parameter
 density=1.0e3
d= 1.e-11
a1= 0.4e9
a2= 0.0
a3= 0.0
b1= 0.2e9
b2= 0.0
b3= 0.0
c2= 0.0
c3= 1.0e9
c4= 0.0
c5= 0.0
c6= 0.0
d2= 0.3e9
d3= 0.0
d4= 0.0
d5= 0.0
d6= 0.0
e2= 0.5e9
e3= 0.0
e4= 0.0
e5= 0.0
e6= 0.0
f2= 0.0
f3= 1.0e9
f4= 0.0
f5= 0.0
f6= 0.0
g2= 0.7e8
g3= 0.0
g4= 0.0
g5= 0.0
g6= 0.0
**
*material,name=kaliske
*density
 <density>
*anisotropic hyperelastic,definition=kaliske,local directions=2
 <a1>,<a2>,<a3>,<b1>,<b2>,<b3>,<c2>,<c3>,
 <c4>,<c5>,<c6>,<d2>,<d3>,<d4>,<d5>,<d6>,
 <e2>,<e3>,<e4>,<e5>,<e6>,<f2>,<f3>,<f4>,
 <f5>,<f6>,<g2>,<g3>,<g4>,<g5>,<g6>,<d>
 
** 
*mullins effect
 2.0,,0.2
*viscoelastic, time=prony 
 0.4,   0.4,  1.43e-4
*equation
 2,
 face14,1,1,ndisp1,1,-1
 2,
 face22,1,1,ndisp2,1,-1
 2,
 face32,1,1,ndisp3,1,-1
 2,
 face42,1,1,ndisp4,1,-1
 2,
 face52,1,1,ndisp5,1,-1
*nset, nset=ndisp
 ndisp1, ndisp2, ndisp3, ndisp4, ndisp5
*boundary,type=velocity
 face11,3
 face13,2
 face16,1
**
 face21,2
 face24,1
**
 face31,2
 face34,1
**
 face41,2
 face44,1
 face41,3
 face43,3
**
 face51,2
 face54,1
 face51,3
 face53,3
**
*amplitude,name=bc
0.0,0.0,0.25,0.75,0.5,0.0,0.75,1.0,
**
*restart, write
**
*step
 uniaxial loading
*dynamic, explicit
 ,0.75
*boundary,amplitude=bc
 ndisp,1,1,0.3
*output,field, variable=preselect
*element output, elset=eall
 s, le
*node out, nset=all
 u, rf
*output,field, number=2
*element output, elset=QA_TEST
 s,
*node output, nset=QA_TEST
 u,
*output, history, variable=preselect
*element output, elset=eall
 s11, le11, le22, le33, ne11, ne22, ener
*node out, nset=ndisp
 rf1, u1
*end step