*HEADING
Test embedded element with 3-D first order surface and solid elements
with model change
*PREPRINT,MODEL=YES,ECHO=YES,HISTORY=YES
*NODE
1,1,0,0
2,4,0,0
3,7,0,0
4,10,0,0
5,1,3,0
6,4,3,0
7,7,3,0
8,10,3,0
9,1,0,-3
10,4,0,-3
11,7,0,-3
12,10,0,-3
13,1,3,-3
14,4,3,-3
15,7,3,-3
16,10,3,-3
17,1,0,-2
18,3.996,0,-2
19,4,3,-2
20,0.98,3,-2
21,8,0,0
22,9,0,0
23,9,3,0
24,8,3,0
25,8,0,-3
26,9,0,-3
27,9,3,-3
28,8,3,-3
*ELEMENT, TYPE=C3D8, ELSET=HOSTSOLID3D
1, 1, 2, 10, 9, 5, 6, 14, 13
2, 2, 3, 11,10, 6, 7, 15, 14
3, 3, 4, 12,11, 7, 8, 16, 15 
*ELEMENT, TYPE=C3D8, ELSET=SOLID3D
4, 21, 22, 26, 25, 24, 23, 27, 28
*ELEMENT, TYPE=SFM3D4, ELSET=MEM3D
5, 17, 18, 19, 20
*SOLID SECTION,ELSET=HOSTSOLID3D,MATERIAL=STEEL
*MATERIAL,NAME=STEEL
*ELASTIC
2.1E11
*SURFACE SECTION,ELSET=MEM3D
*REBAR LAYER,ORI=REBAR
REB,0.3,0.5,0,MAT2,20.0,1
REB2,0.3,0.5,0,MAT2,60.0,1
*MATERIAL,NAME=STEEL2
*ELASTIC
5.0E11
*ORIENTATION, NAME=REBAR, DEFINITION=OFFSET TO NODES
2, 4
3
*MATERIAL,NAME=MAT2
*ELASTIC
 7.0E11,
*SOLID SECTION,ELSET=SOLID3D,MATERIAL=STEEL2
*EMBEDDED ELEMENT,HOST ELSET=HOSTSOLID3D,ROUND OFF TOLERANCE=0.005, EXTERIOR TOLERANCE=0.01
SOLID3D,MEM3D
*BOUNDARY
1,1,3
5,1,3
9,1,3
13,1,3
*RESTART,WRITE,FREQ=1
*STEP,NLGEOM
 STEP #1: DISPLACE TEST SPECIMEN SO FINAL SHAPE IS
          IDENTICAL TO INITIAL SHAPE OF REFERENCE ELEMENT 
**********************************************************
*STATIC
 0.25,1.
*PRINT,MODEL CHANGE=YES
*BOUNDARY
4,2,2,-0.1
8,2,2,-0.1
12,2,2,-0.1
16,2,2,-0.1
*ENERGY PRINT
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=MEM3D,REBAR
S
LE
*ELEMENT OUTPUT,ELSET=SOLID3D
S
LE
*ELEMENT OUTPUT,ELSET=HOSTSOLID3D
S
LE
*ELEMENT OUTPUT,ELSET=MEM3D,REBAR,POSITION=CENTRO
S
LE
*EL FILE,ELSET=MEM3D,REBAR
S
E
*EL FILE,ELSET=SOLID3D
S
E
*EL FILE,ELSET=HOSTSOLID3D
S
E
*EL PRINT,ELSET=MEM3D,REBAR
S
E
*EL PRINT,ELSET=MEM3D,REBAR,POSITION=CENTRO
S
E
*EL PRINT,ELSET=SOLID3D
S
E
*EL PRINT,ELSET=HOSTSOLID3D
S
E
*NODE PRINT
U
RF
*END STEP
*STEP,NLGEOM
 STEP #2: REMOVE TEST ELEMENTS
*STATIC
 1.,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
 MEM3D,
*END STEP
**
**
*STEP,NLGEOM
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*STATIC
 1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
 MEM3D, 
*END STEP
*STEP,nlgeom
*STATIC
0.2,1.0
*CLOAD
4,2,-1.0E6
8,2,-1.0E6
12,2,-1.0E6
16,2,-1.0E6
*ENERGY PRINT
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=MEM3D,REBAR
S
LE
*ELEMENT OUTPUT,ELSET=SOLID3D
S
LE
*ELEMENT OUTPUT,ELSET=HOSTSOLID3D
S
LE
*ELEMENT OUTPUT,ELSET=MEM3D,REBAR,POSITION=CENTRO
S
LE
*EL FILE,ELSET=MEM3D,REBAR
S
E
*EL FILE,ELSET=SOLID3D
S
E
*EL FILE,ELSET=HOSTSOLID3D
S
E
*EL PRINT,ELSET=MEM3D,REBAR
S
E
*EL PRINT,ELSET=MEM3D,REBAR,POSITION=CENTRO
S
E
*EL PRINT,ELSET=SOLID3D
S
E
*EL PRINT,ELSET=HOSTSOLID3D
S
E
*NODE PRINT
U
RF
*END STEP
*STEP,NLGEOM
 STEP #2: REMOVE TEST ELEMENTS
*STATIC
 1.,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
 MEM3D,
*END STEP
*STEP,NLGEOM,INC=100
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
 .025,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
 MEM3D,
*END STEP