*HEADING
 :  CGAX6M MODEL CHANGE VERIFICATION, HYPERELASTIC MATERIAL
*NODE,NSET=SHORT1
2,2.,0.5
202,1.,0.5
402,0.0,0.5
*NODE,NSET=SHORT2
10,2.,4.5
210,1.,4.5
410,0.0,4.5
*NFILL,NSET=SHORT
SHORT1,SHORT2,4,2
*NODE
1,2.
11,2.,5.
401,
411,,5.
*NGEN,NSET=RHS
1,11,1
*NGEN,NSET=LHS
401,411,1
*NFILL,NSET=NALL
RHS,LHS,4,100
*NSET, NSET=YFIX, GEN
1, 401, 100
*NSET, NSET=XFIX
1, 
*NSET, NSET=XFIX2
11, 
*NSET, NSET=YMOVE, GEN
11, 411, 100
*NSET, NSET=XMOVE, GEN
204, 208, 1
4, 8, 1
*NSET, NSET=XMOVE, GEN
104, 108, 1
304, 308, 1
404, 408, 1
*ELEMENT,TYPE=CGAX6M,ELSET=EALL
1  ,1,3,201, 2,102,101
101,201,3,203, 102,103,202
*ELGEN,ELSET=EALL
1,5,2,1,2,200,5
101,5,2,1,2,200,5
*ELSET, ELSET=EREMOV
2,3,4,7,8,9,
102,103,104,107,108,109
*SOLID SECTION,MATERIAL=RUBBER,ELSET=EALL
1.,
*MATERIAL,NAME=RUBBER
*HYPERELASTIC, N=1
56.00E4, 14.00E4, 1.43E-5
*CONDUCTIVITY
1., 
*BOUNDARY
YMOVE, 5, 5
YFIX, 2, 2
YFIX, 5, 5
*INITIAL CONDITIONS, TYPE=TEMPERATURE
NALL,0.0
**
*STEP,NLGEOM
 STEP 1: COMPRESS ELEMENTS
*STATIC
0.1,1.,0.1,0.1
*BOUNDARY
YMOVE, 2,2 ,-0.1
*NODE PRINT,FREQ=999
U, 
RF, 
*EL PRINT,FREQ=999
E, 
S, 
*NODE FILE,FREQ=999
U, RF
*EL FILE,FREQ=999
E, NE, S
*END STEP
**
*STEP,NLGEOM
STEP 2: REMOVE MIDDLE ELEMENTS
*STATIC
1.,1.
** CONTROL PARAMETERS NEEDED FOR ABAQUS QA PURPOSES ONLY
*CONTROLS,PARAMETERS=FIELD
1.e-6,,,1.,,1e-10
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
EREMOV,
*END STEP
**
*STEP,NLGEOM
 STEP 3: MOVE TOP EDGE DOWN
*STATIC
1.,1.
*CONTROLS,RESET
*BOUNDARY
XMOVE, 1,1, 0.0
*END STEP
**
*STEP,NLGEOM
 STEP 4: ADD ELEMENTS BACK IN
*STATIC
0.1,1.,0.1,0.1
*MODEL CHANGE,ADD,TYPE=ELEMENT
EREMOV,
*BOUNDARY,OP=NEW
YFIX, 2, 2
YFIX, 5, 5
YMOVE, 2,2 ,0.1
*END STEP