*HEADING
 BEAM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 1: PIPE21
**
*NODE,NSET=NALL
*NSET,NSET=LEFT
*NSET,NSET=RLEFT
*NSET,NSET=TLEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=PIPE21,ELSET=TEST1
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
**
** TEST 2: PIPE22
**
*NODE,NSET=NALL
*NSET,NSET=LEFT
*NSET,NSET=RLEFT
*NSET,NSET=TLEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=PIPE22,ELSET=TEST2
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
**
** TEST 3: PIPE21H
**
*NODE,NSET=NALL
*NSET,NSET=LEFT
*NSET,NSET=RLEFT
*NSET,NSET=TLEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=PIPE21H,ELSET=TEST3
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
**
** TEST 4: PIPE22H
**
*NODE,NSET=NALL
*NSET,NSET=LEFT
*NSET,NSET=RLEFT
*NSET,NSET=TLEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=PIPE22H,ELSET=TEST4
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
**
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*BEAM SECTION,SECTION=PIPE,TEMPERATURE=GRADIENTS,MATERIAL=ELAS,ELSET=EALL,ROTARY INERTIA=ISO
1.0,.05
0.,0.,-1.
*NORMAL
TEST,LEFT, -1., 0.,0.
TEST,MID,  -1., 0.,0.
TEST,RIGHT,-1., 0.,0.
*AQUA
*WIND
**   rho,    Z0,     c_x,    c_y,   dx,  dy,  alpha
*AMPLITUDE, NAME=TEN
0.,10.,  1.,10.
*AMPLITUDE, NAME= A1_5
0., 1.5, 1., 1.5
*AMPLITUDE,NAME=ZERO
0.,0.,  1.,0.
*MATERIAL,NAME=ELAS
*ELASTIC
30.E9, .3
*EXPANSION
1.E-10,
*DENSITY
2.0,
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,70.,0.
**
**
*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,OP=NEW
LEFT,1,2
RLEFT,6
TLEFT,6,6,-.7853982
RIGHT,1,1,10.0
RIGHT,2,2,0.0
RIGHT,6,6,-.7853982
MID,1,1,5.0
MID,2,2,0.0
MID,6,6,-.7853982
*NODE PRINT,FREQ=999
RF,
*EL PRINT,POSITION=CENTROID,FREQ=999
S,E
*NODE FILE,FREQ=999
RF,
*EL FILE,POSITION=CENTROID,FREQ=999
E, S
*OUTPUT,HISTORY,VARIABLE=PRESELECTED,FREQUENCY=999
*OUTPUT,FIELD,FREQUENCY=999
*NODE OUTPUT
U,
*ELEMENT OUTPUT
E,
*END STEP
**
**
*STEP,NLGEOM
 STEP #2: REMOVE TEST ELEMENTS
*STATIC
1.,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*STATIC
0.25,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST,
*BOUNDARY,OP=NEW,TYPE=VELOCITY
LEFT,1,2
LEFT,6
*END STEP
**
**
*STEP,NLGEOM
 STEP #4: APPLY DLOAD TYPE P2
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST, P2,2.
REFER,P2,2.
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: APPLY DLOAD TYPE GRAV
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST, GRAV,150.,0.,1.,0.
REFER,GRAV,150.,0.,1.,0.
*OUTPUT,HISTORY,VARIABLE=PRESELECTED,FREQUENCY=999
*OUTPUT,FIELD,FREQUENCY=999
*NODE OUTPUT
U,RF
*ELEMENT OUTPUT
S,
*END STEP
**
**
*STEP,NLGEOM
 STEP #6: APPLY DLOAD TYPE PX
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST, PX, 10.
REFER,PX, 10.
*END STEP
**
**
*STEP,NLGEOM
 STEP #7: APPLY DLOAD TYPE PY
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST, PY, 10.
REFER,PY, 10.
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: APPLY DLOAD TYPE PY
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST, PY, 10.
REFER,PY, 10.
*END STEP
**
**
*STEP,NLGEOM
 STEP #9: APPLY DLOAD TYPE PI
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST, PI, 10., .8
REFER,PI, 10., .8
*OUTPUT,HISTORY,VARIABLE=PRESELECTED,FREQUENCY=999
*OUTPUT,FIELD,FREQUENCY=999
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
*END STEP
**
**
*STEP,NLGEOM
 STEP #10: APPLY DLOAD TYPE PE
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST, PE, 10., 1.2
REFER,PE, 10., 1.2
*END STEP
**
**
*STEP,NLGEOM
 STEP #11: APPLY DLOAD TYPE PB
*STATIC
0.25,1.
*CONTROLS,PARAMETERS=FIELD
,10.
*DLOAD,OP=NEW
TEST, PB,,2.25,1.25,1.9,2500.
REFER,PB,,2.25,1.25,1.9,2500.
*END STEP
**
**
*STEP,NLGEOM
 STEP #12: INTERMEDIATE STEP, REMOVE DRAG LOAD
*STATIC
1.,1.
*CONTROLS,PARAMETER=FIELD,FIELD=DISPLACEMENT
 , ,1.11
*CONTROLS,PARAMETER=FIELD,FIELD=ROTATION
 , ,12.8
*DLOAD,OP=NEW
*OUTPUT,HISTORY,VARIABLE=PRESELECTED,FREQUENCY=999
*OUTPUT,FIELD,FREQUENCY=999
*NODE OUTPUT
U,
*ELEMENT OUTPUT
E,
*END STEP
**
**
*STEP,NLGEOM
 STEP #13: APPLY DLOAD TYPE FDD
*DYNAMIC
.01,.01
*CONTROLS,RESET
*BOUNDARY,TYPE=VELOCITY,OP=NEW
LEFT,1,1, 10.0
LEFT,2,2,-20.0
LEFT,6,6
*DLOAD,OP=NEW
TEST,  FDD,   2.0, 7.5491, 1.2,   0.8,  TEN
REFER, FDD,   2.0, 7.5491, 1.2,   0.8,  TEN
*END STEP
**
**
*STEP,NLGEOM
 STEP #14: APPLY DLOAD TYPE FDT
*DYNAMIC
.01,.01
*BOUNDARY,TYPE=VELOCITY,OP=NEW
LEFT,1,1, 10.0
LEFT,2,2,-20.0
LEFT,6,6
*DLOAD,OP=NEW
TEST,  FDT,      ,  7.5491, 0.0381972,   0.8, 1.5, TEN
TEST,  FDT,      ,  7.5491, 0.0381972,   0.8,    , TEN
REFER, FDT,      ,  7.5491, 0.0381972,   0.8, 1.5, TEN
REFER, FDT,      ,  7.5491, 0.0381972,   0.8,    , TEN
*OUTPUT,HISTORY,VARIABLE=PRESELECTED,FREQUENCY=999
*OUTPUT,FIELD,FREQUENCY=999
*NODE OUTPUT
U,RF
*ELEMENT OUTPUT
S,
*END STEP
**
**
*STEP,NLGEOM
 STEP #15: APPLY DLOAD TYPE FI
*DYNAMIC
.01,.01
*BOUNDARY,TYPE=ACCELERATION,OP=NEW
LEFT,1,1, 2.0
LEFT,2,2,-6.0
LEFT,6
*DLOAD,OP=NEW
TEST,   FI,      ,  1.2, 0.0,   0.7
REFER,  FI,      ,  1.2, 0.0,   0.7
**TEST,   FI,      ,  7.4111, 0.0,   0.7
**REFER,  FI,      ,  7.4111, 0.0,   0.7
*END STEP
**
**
*STEP,NLGEOM
 STEP #16: INTERMEDIATE STEP, MOVE BEAM UPWARDS, CLEAR LOAD
*DYNAMIC
.01,.01
*BOUNDARY,TYPE=ACCELERATION,OP=NEW
LEFT,1,1, 2.0
LEFT,2,2,-6.0
LEFT,6,6
*DLOAD,OP=NEW
*END STEP
**
**
*STEP,NLGEOM
 STEP #17: APPLY DLOAD TYPE WDD
*DYNAMIC
.01,.01
*BOUNDARY,TYPE=VELOCITY,OP=NEW
LEFT,1,1
LEFT,2,2, 25.
LEFT,6
*DLOAD,OP=NEW
TEST, WDD,      , 7.5491, 1.2,   0.5,  A1_5, ZERO
REFER,WDD,      , 7.5491, 1.2,   0.5,  A1_5, ZERO
*END STEP
**
**
*STEP,NLGEOM
 STEP #18: APPLY DLOAD TYPE WD1 AND WD2
*DYNAMIC
.01,.01
*BOUNDARY,TYPE=VELOCITY,OP=NEW
LEFT,1,1,-20.0
LEFT,2
LEFT,6
*DLOAD,OP=NEW
TEST,  WD1,      , 7.5491, 1.2,       ,
TEST,  WD2,      , 7.5491, 1.2,       ,
REFER, WD1,      , 7.5491, 1.2,       ,
REFER, WD2,      , 7.5491, 1.2,       ,
*END STEP