*HEADING
 BASIC VERIFICATION - MPC TYPES SS BILINEAR, SSF BILINEAR AND SLIDER
 ELEMENTS: S8R AND C3D20, TEST with RIKS
*NODE
1, 0.,-.125,1.
3, 2.,-.125,1.
7, 0.,-.125,0.
9, 2.,-.125,0.
21, 0.,.125,1.
23, 2.,.125,1.
27, 0.,.125,0.
29, 2.,.125,0.
*NODE
101, 2.,0.,1.
105,10.,0.,1.
121, 2.,0.,0.
125,10.,0.,0.
*NGEN
1,3
7,9
21,23
27,29
1,7,3
3,9,3
21,27,3
23,29,3
1,21,10
3,23,10
7,27,10
9,29,10
*NGEN
101,105
121,125
101,121,10
105,125,10
111,115
*ELEMENT,TYPE=C3D20,ELSET=SOLID
1, 1,3,9,7,21,23,29,27,2,6,8,4,22,26,28,
   24,11,13,19,17
*SOLID SECTION,ELSET=SOLID,MATERIAL=MAT
*ELEMENT,TYPE=S8R
11, 101,103,123,121,102,113,122,111
*ELGEN,ELSET=SHELL
11,2,2,1
*SHELL SECTION,ELSET=SHELL,MATERIAL=MAT
.25,3
*MATERIAL,NAME=MAT
*ELASTIC
30.E6,.3
*NSET, NSET=SHELL_N1
111, 
*NSET,NSET=SHELL_N2
101, 
*NSET,NSET=INDP_N1
26, 
*NSET, NSET=INDP_N2
13, 
*TRANSFORM, NSET=SHELL_N1
** 30 ABOUT Y AXIS, THEN 60 ABOUT (TRANSF.) Z AXIS
0.4330127,0.866025,-0.25,-0.75,0.5,0.4330127
*TRANSFORM, NSET=SHELL_N2
** 30 ABOUT Y AXIS, THEN 45 ABOUT (TRANSF.) X AXIS
0.8660254,0.0,-0.5,0.35355339,0.7971,0.61237243
*TRANSFORM, NSET=INDP_N1
** 45 ABOUT Z AXIS, THEN 45 ABOUT (TRANSF) X AXIS
0.7071,0.7071,0.,-0.5,0.6,0.7071
*TRANSFORM, NSET=INDP_N2
0.7071,0.,-0.7071,0.,1.,0.
*MPC
SS BILINEAR, 101,3,13,23
SS BILINEAR, 121,9,19,29
SSF BILINEAR, 111,3,6,9,13,19,23,26,29
SLIDER, 13,3,23
SLIDER, 19,9,29
*NSET,NSET=SYM,GEN
7,9
17,17
27,29
122,125
*NSET,NSET=CLAMPED,GEN
1,7,3
11,17,3
21,27,3 
*NSET,NSET=EDGE
105,115,125
*BOUNDARY
CLAMPED,1,3
SYM,ZSYMM
**121,4,5
************************************
*STEP,NLGEOM
 Geometric nonlinear
*STATIC,RIKS
0.2,1.0,,,1.0
*CLOAD
105,2,-10.
115,2,-40.
125,2,-10.
*EL PRINT,FREQ=0
*NODE FILE,FREQ=999
 U,
 RF, 
*NODE PRINT,FREQ=999, GLO=YES
 U,
 RF, 
*END STEP