*HEADING
C3D27RH ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
**
*NODE,NSET=NALL
**
**
*NSET,NSET=BOT
*NSET,NSET=TIP
*NSET,NSET=DEFBC
*NSET,NSET=MIDSIDEW
*NSET,NSET=MIDSIDE
*ELEMENT,TYPE=C3D27RH,ELSET=TEST1
**
**
**
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
**
**
**
*ELSET,ELSET=EALL
TEST,REFER
*SOLID SECTION,MATERIAL=ELAS,ELSET=EALL
*MATERIAL,NAME=ELAS
*ELASTIC
1.E10, .25
*EXPANSION
1.E-3,
*DENSITY
8.E3,
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,70.,0.
*AMPLITUDE,NAME=RAMP
0.,0.,1.,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,OP=NEW
BOT,PINNED
DEFBC,2,2,-0.03
DEFBC,1,1,0.0
DEFBC,3,3,0.0
MIDSIDE,2,2,-0.015
MIDSIDE,1,1,0.0
MIDSIDE,3,3,0.0
*NODE PRINT,FREQ=999
RF,
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
E,
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF,
*ELEMENT OUTPUT,POSITION=CENTROID
E,
S,
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT, NSET=NALL
RF,
*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
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST,
*BOUNDARY,OP=NEW
BOT,PINNED
*END STEP
**
**
*STEP,NLGEOM
 STEP #4: APPLY LOAD TO BOTH ELEMENTS
*STATIC
0.25,1.
*DLOAD
TEST, GRAV,150.,0.,1.,0.
REFER,GRAV,150.,0.,1.,0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: REMOVE TEST ELEMENT
*STATIC
.25,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
.25,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #7: APPLY THERMAL LOAD
*STATIC
.25,1.0
*BOUNDARY,OP=NEW,FIXED
BOT,1,3
TIP,1,3
MIDSIDEW,1,3
*TEMPERATURE
NALL,100.
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
THE,
E,
TEMP,
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S,
THE,
E,
TEMP,
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos
200, ,
*BOUNDARY,OP=NEW
BOT,1,3
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*OUTPUT,FIELD,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
 STEP #9: MODAL DYNAMIC
*MODAL DYNAMIC
.001,.01
*CLOAD
TIP,1,.1
*NODE PRINT,FREQ=999
RF,
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
E,
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF,
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT, NSET=NALL
RF,
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,POSITION=CENTROID
E,
S,
*END STEP
**
**
*STEP,PERTURBATION
 STEP #10: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
90.,400., 5
*CLOAD
TIP,1,0.1
*NODE PRINT,FREQ=999
RF,
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
E,
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF,
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT, NSET=NALL
RF,
*OUTPUT,FIELD,FREQ=0
*ELEMENT OUTPUT,POSITION=CENTROID
*END STEP
**
**
**
*STEP,PERTURBATION
 STEP #11: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
90.,400., 5
*BOUNDARY,OP=NEW
BOT,1,3
*CLOAD
TIP,1,0.1
*NODE PRINT,FREQ=999
RF,
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
E,
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF,
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT,NSET=NALL
RF,
*OUTPUT,FIELD, FREQ=0
*ELEMENT OUTPUT,POSITION=CENTROID
*END STEP
**
**
**
*STEP,PERTURBATION
 STEP #12: STEADY-STATE DYNAMICS,SUBSPACE PROJECTION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
90.,400., 5
*CLOAD
TIP,1,0.1
*NODE PRINT,FREQ=999
RF,
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
E,
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF,
*END STEP
**
**
**
*STEP,NLGEOM
 STEP #13: DYNAMICS
*DYNAMIC
.001, .01
*BOUNDARY,OP=NEW
BOT,PINNED
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
E,
*NODE PRINT,FREQ=999
RF,
*CLOAD,AMP=RAMP
TIP,1,.1
*END STEP