*HEADING
 CONTINUUM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
 elements with pore pressure
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 26.5: CPE4RP
**
*NODE,NSET=NALL
301,  15.0,  9.0
302,  16.0,  9.0
303,  16.0,  8.0
304,  15.0,  8.0
305,  17.0,  9.0
306,  17.0,  8.0
307,  17.0,  8.5
308,  15.0,  8.5
309,  15.0,  6.0
310,  16.0,  6.0
311,  16.0,  5.0
312,  15.0,  5.0
313,  17.0,  6.0
314,  17.0,  5.0
315,  17.0,  5.5
316,  15.0,  5.5
*NSET,NSET=LEFT
301,304,308,309,312,316
*NSET,NSET=MID
310,311
*NSET,NSET=RIGHT
313,314,315
*NSET,NSET=TIP
305,306,307,313,314,315
*NSET,NSET=REFER
301,304,305,306
*NSET,NSET=TEST
309,312,313,314
*ELEMENT,TYPE=CPE4RP,ELSET=TEST265
301, 301,304,306,305
302, 309,312,314,313
*ELSET,ELSET=REFER
301, 
*ELSET,ELSET=TEST
302, 
**
** TEST 27: CPE8P
**
*NODE,NSET=NALL
401,  15.0,  0.0
402,  16.0,  0.0
403,  16.0, -1.0
404,  15.0, -1.0
405,  17.0,  0.0
406,  17.0, -1.0
407,  17.0, -0.5
408,  15.0, -0.5
409,  15.0, -3.0
410,  16.0, -3.0
411,  16.0, -4.0
412,  15.0, -4.0
413,  17.0, -3.0
414,  17.0, -4.0
415,  17.0, -3.5
416,  15.0, -3.5
*NSET,NSET=LEFT
401,404,408,409,412,416
*NSET,NSET=MID
410,411
*NSET,NSET=RIGHT
413,414,415
*NSET,NSET=TIP
405,406,407,413,414,415
*NSET,NSET=REFER
401,404,405,406
*NSET,NSET=TEST
409,412,413,414
*ELEMENT,TYPE=CPE8P,ELSET=TEST27
401, 401,404,406,405,408,403,407,402
402, 409,412,414,413,416,411,415,410
*ELSET,ELSET=REFER
401, 
*ELSET,ELSET=TEST
402, 
**
** TEST 28: CPE8PH
**
*NODE,NSET=NALL
501,  15.0, -9.0
502,  16.0, -9.0
503,  16.0, -10.0
504,  15.0, -10.0
505,  17.0, -9.0
506,  17.0, -10.0
507,  17.0, -9.5
508,  15.0, -9.5
509,  15.0, -12.0
510,  16.0, -12.0
511,  16.0, -13.0
512,  15.0, -13.0
513,  17.0, -12.0
514,  17.0, -13.0
515,  17.0, -12.5
516,  15.0, -12.5
*NSET,NSET=LEFT
501,504,508,509,512,516
*NSET,NSET=MID
510,511
*NSET,NSET=RIGHT
513,514,515
*NSET,NSET=TIP
505,506,507,513,514,515
*NSET,NSET=REFER
501,504,505,506
*NSET,NSET=TEST
509,512,513,514
*ELEMENT,TYPE=CPE8PH,ELSET=TEST28
501, 501,504,506,505,508,503,507,502
502, 509,512,514,513,516,511,515,510
*ELSET,ELSET=REFER
501, 
*ELSET,ELSET=TEST
502, 
**
** TEST 29: CPE8RP
**
*NODE,NSET=NALL
601,  15.0, -18.0
602,  16.0, -18.0
603,  16.0, -19.0
604,  15.0, -19.0
605,  17.0, -18.0
606,  17.0, -19.0
607,  17.0, -18.5
608,  15.0, -18.5
609,  15.0, -21.0
610,  16.0, -21.0
611,  16.0, -22.0
612,  15.0, -22.0
613,  17.0, -21.0
614,  17.0, -22.0
615,  17.0, -21.5
616,  15.0, -21.5
*NSET,NSET=LEFT
601,604,608,609,612,616
*NSET,NSET=MID
610,611
*NSET,NSET=RIGHT
613,614,615
*NSET,NSET=TIP
605,606,607,613,614,615
*NSET,NSET=REFER
601,604,605,606
*NSET,NSET=TEST
609,612,613,614
*ELEMENT,TYPE=CPE8RP,ELSET=TEST29
601, 601,604,606,605,608,603,607,602
602, 609,612,614,613,616,611,615,610
*ELSET,ELSET=REFER
601, 
*ELSET,ELSET=TEST
602, 
**
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*SOLID SECTION,MATERIAL=ELAS,ELSET=EALL
1.0, 
*HOURGLASS STIFFNESS
1.E3
*MATERIAL,NAME=ELAS
*ELASTIC
1000., .25
*PERMEABILITY,SPECIFIC=1.0
 1.E-5, 
*EXPANSION
1.E-3, 
*DENSITY
1.E-3, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,70.
*INITIAL CONDITIONS,TYPE=RATIO
NALL,1.0
*INITIAL CONDITION,TYPE=PORE PRESSURE
REFER,10
*INITIAL CONDITION,TYPE=PORE PRESSURE
TEST,20
*AMPLITUDE,NAME=RAMP
0.,0.,  1.,1.
**
**
*STEP, UNSYMM=YES
 STEP #1: DISPLACE TEST SPECIMEN SO FINAL SHAPE IS
          IDENTICAL TO INITIAL SHAPE OF REFERENCE ELEMENT 
*SOILS
 1, 1.
*PRINT,MODEL CHANGE=YES
*BOUNDARY,OP=NEW
LEFT,1,2
TEST,8,8,10
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S, 
E, 
POR, 
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=999
E, S
*END STEP
**
**
*STEP
 STEP #2: REMOVE TEST ELEMENTS
*SOILS
 1., 1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*SOILS
 1., 1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST, 
*BOUNDARY,OP=NEW
LEFT,1,2
*END STEP
**
**
*STEP
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*SOILS
 0.1, 1.
*DLOAD
TEST,P4,-25
REFER,P4,-25
*END STEP
**
**
*STEP
 STEP #5: REMOVE TEST ELEMENT
*SOILS
 1,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*SOILS
1,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #7: APPLY THERMAL LOAD
*SOILS
.1,1.0
*BOUNDARY,OP=NEW,FIXED
TIP,1
LEFT,1,2
LEFT,8,8
*TEMPERATURE
NALL,100.
*EL PRINT,FREQ=999
S, 
POR, 
THE,TEMP
*EL FILE,FREQ=999
S, 
THE,TEMP
*END STEP
**