*HEADING
Plane Stress FeFp
Elements: cps4r, cps3, cps6m, cps8, s4r, s3r, sc8r, m3d4r
Biaxial tension
Marlow model, incompressible, mullins effect
*NODE,NSET=one
1,
2,1.
3,1.,1.,
4,0.,1.,
**
11,2.
12,3.
13,3.,1.
14,2.,1.
**
21,4.0,0.0
22,5.0,0.0
23,5.0,1.0
24,4.0,1.0
25,4.5,0.0
26,5.0,0.5
27,4.5,1.0
28,4.0,0.5
29,4.5,0.5
**
31,6.0,0.0
32,7.0,0.0
33,7.0,1.0
34,6.0,1.0
35,6.5,0.0
36,7.0,0.5
37,6.5,1.0
38,6.0,0.5
**
41,0.0,3.0
42,1.0,3.0
43,1.0,4.0
44,0.0,4.0
**
51,2.0,3.0
52,3.0,3.0
53,3.0,4.0
54,2.0,4.0
**
61,4.0,3.0,0.0
62,5.0,3.0,0.0
63,5.0,4.0,0.0
64,4.0,4.0,0.0
65,4.0,3.0,1.0
66,5.0,3.0,1.0
67,5.0,4.0,1.0
68,4.0,4.0,1.0
**
71,6.0,3.0
72,7.0,3.0
73,7.0,4.0
74,6.0,4.0
**
*NSET, NSET=SHELL-NODE
41,42,43,44
51,52,53,54
61,62,63,64
*NSET, NSET=SC8R-FRONT
65,66,67,68
*NSET, NSET=BOT
1,2
11,12
21,22,25
31,32,35
41,42
51,52
61,62,65,66
71,72
*NSET, NSET=LEFT
1,4
11,14
21,24,28
31,34,38
41,44
51,54
61,64,65,68
71,74
*NSET, NSET=TOP
3,4
13,14
43,44
53,54
**63,64,67,68
73,74
*NSET, NSET=TOP-1-CPS6M
23,24
*NSET, NSET=TOP-1-8NODE
33,34
*NSET, NSET=TOP-2-CPS6M
27
*NSET, NSET=TOP-2-8NODE
37
*NSET, NSET=TOP-3DEL
63,64,67,68
*NSET, NSET=RIGHT-3DEL
62,63,66,67
*NSET, NSET=RIGHT
2,3
12,13
42,43
52,53
72,73
*NSET, NSET=RIGHT-1-CPS6M
22,23
*NSET, NSET=RIGHT-1-8NODE
32,33
*NSET, NSET=RIGHT-2-CPS6M
26
*NSET, NSET=RIGHT-2-8NODE
36
*ELEMENT,TYPE=CPS4R,ELSET=SOLID
1,1,2,3,4
*ELEMENT,TYPE=CPS3,ELSET=SOLID
2,11,12,13
3,11,13,14
*ELEMENT,TYPE=CPS6M,ELSET=SOLID
21,21,22,23,25,26,29
22,21,23,24,29,27,28
*ELEMENT,TYPE=CPS8,ELSET=SOLID
31,31,32,33,34,35,36,37,38
*ELEMENT,TYPE=S4R,ELSET=SHELL
41,41,42,43,44
*ELEMENT,TYPE=S3R,ELSET=SHELL
51,51,52,53
52,51,53,54
*ELEMENT,TYPE=SC8R,ELSET=SHELL
61,61,62,63,64,65,66,67,68
*ELEMENT,TYPE=M3D4R,ELSET=MEMBRANE
71,71,72,73,74
*ELSET, ELSET=ALLE
 SOLID
 SHELL
*SOLID SECTION,ELSET=SOLID,MATERIAL=mat-129
1.
*SHELL SECTION,ELSET=SHELL,MATERIAL=mat-129
1.,5
*MEMBRANE SECTION,ELSET=MEMBRANE,MATERIAL=mat-129
1.,5
*MATERIAL,NAME=mat-129
*HYPERELASTIC,marlow,poisson=0.5
*include, input=heplmu_matprops_bi.inp
*amplitude,name=nominal_stress
0.0,0.0, 10.,2.44444, 20.0,0.0, 40.,3.959985,
60.,0.0, 90.0,6.21102
*STEP,NLGEOM,INC=2000
*STATIC
1.0,90.,,1.0
*BOUNDARY
BOT,2
LEFT,1
SHELL-NODE,3,6
SC8R-FRONT,4,6
*cload, amplitude=nominal_stress
RIGHT,1,0.5
TOP,  2,0.5
TOP-3DEL,2,0.25
RIGHT-3DEL,1,0.25
** for 2nd-order elements
RIGHT-1-CPS6M,1,0.25
RIGHT-2-CPS6M,1,0.50
RIGHT-1-8NODE,1,0.1666666667
RIGHT-2-8NODE,1,0.6666666666
TOP-1-CPS6M,  2,0.25
TOP-2-CPS6M,  2,0.50
TOP-1-8NODE,  2,0.1666666667
TOP-2-8NODE,  2,0.6666666666
*output,history,freq=2
*element output, elset=ALLE
 S22
*output,field,freq=1000
*node output
 u,rf
*element output
 S,LE,PE
 PEEQ
*END STEP