*HEADING
 BEAM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 1: B31
**
*NODE,NSET=NALL
1,  0.0, 1.
2,  1.0, 1.
3,  2.0, 1.
11, 0.0, 0.0, .0
12, 1.5, -.1, .1
13, 3.0, -.3, .2
*NSET,NSET=LEFT
1,11
*NSET,NSET=MID
12, 
*NSET,NSET=RIGHT
13, 
*NSET,NSET=TIP
3,13
*ELEMENT,TYPE=B31,ELSET=TEST1
1, 1,2
2, 2,3
11,11,12
12,12,13
*ELSET,ELSET=REFER
1,2
*ELSET,ELSET=TEST
11,12
*BEAM SECTION,SECTION=RECT,TEMPERATURE=GRADIENTS,MATERIAL=ELAS,ELSET=TEST1
1.,2.
0.,0.,-1.
**
** TEST 2: B32
**
*NODE,NSET=NALL
101, 5.0, 1.
102, 6.0, 1.
103, 7.0, 1.
111, 5.0, 0., .0
112, 6.5, -.1,.1
113, 8.0, -.3,.2
*NSET,NSET=LEFT
101,111
*NSET,NSET=MID
112, 
*NSET,NSET=RIGHT
113, 
*NSET,NSET=TIP
103,113
*ELEMENT,TYPE=B32,ELSET=TEST2
101,101,102,103
111,111,112,113
*ELSET,ELSET=REFER
101,102
*ELSET,ELSET=TEST
111,112
*BEAM SECTION,SECTION=RECT,TEMPERATURE=GRADIENTS,MATERIAL=ELAS,ELSET=TEST2
1.,2.
0.,0.,-1.
**
** TEST 3: B33
**
***NODE,NSET=NALL
**201, 10.0, 1.
**202, 11.0, 1.
**203, 12.0, 1.
**211, 10.0, 0., .0
**212, 11.5, -.1,.1
**213, 13.0, -.3,.2
***NSET,NSET=LEFT
**201,211
***NSET,NSET=MID
**212
***NSET,NSET=RIGHT
**213
***NSET,NSET=TIP
**203,213
***ELEMENT,TYPE=B33,ELSET=TEST3
**201,201,202
**202,202,203
**211,211,212
**212,212,213
***ELSET,ELSET=REFER
**201,202
***ELSET,ELSET=TEST
**211,212
***BEAM SECTION,SECTION=RECT,TEMPERATURE=GRADIENTS,MATERIAL=ELAS,ELSET=TEST3
**1.,2.
**0.,0.,-1.
**
** TEST 4: B34
**
***NODE,NSET=NALL
**301, 15.0, 1.
**302, 16.0, 1.
**303, 17.0, 1.
**311, 15.0,  .0, .0
**312, 16.5, -.1, .1
**313, 18.0, -.3, .2
***NSET,NSET=LEFT
**301,311
***NSET,NSET=MID
**312
***NSET,NSET=RIGHT
**313
***NSET,NSET=TIP
**303,313
***ELEMENT,TYPE=B34,ELSET=TEST4
**301, 301,302
**302, 302,303
**311, 311,312
**312, 312,313
***ELSET,ELSET=REFER
**301,302
***ELSET,ELSET=TEST
**311,312
***BEAM SECTION,SECTION=RECT,TEMPERATURE=GRADIENTS,MATERIAL=ELAS,ELSET=TEST4
**1.,2.
**0.,0.,-1.
**
** TEST 5: B31H
**
*NODE,NSET=NALL
401, 0., 4.
402, 1., 4.
403, 2., 4.
411, 0.0, 3.0, .0
412, 1.5, 2.9, .1
413, 3.0, 2.7, .2
*NSET,NSET=LEFT
401,411
*NSET,NSET=MID
412, 
*NSET,NSET=RIGHT
413, 
*NSET,NSET=TIP
403,413
*ELEMENT,TYPE=B31H,ELSET=TEST5
401, 401,402
402, 402,403
411, 411,412
412, 412,413
*ELSET,ELSET=REFER
401,402
*ELSET,ELSET=TEST
411,412
*BEAM SECTION,SECTION=RECT,TEMPERATURE=GRADIENTS,MATERIAL=ELAS,ELSET=TEST5
1.,2.
0.,0.,-1.
**
** TEST 6: B31OS
**
*NODE,NSET=NALL
501, 5.0,  4.0
502, 6.0,  4.0
503, 7.0,  4.0
511, 5.0,  3.0, .0
512, 6.5,  2.9, .1
513, 8.0,  2.7, .2
*NSET,NSET=LEFT
501,511
*NSET,NSET=MID
512, 
*NSET,NSET=RIGHT
513, 
*NSET,NSET=TIP
503,513
*ELEMENT,TYPE=B31OS,ELSET=TEST6
501, 501,502
502, 502,503
511, 511,512
512, 512,513
*ELSET,ELSET=REFER
501,502
*ELSET,ELSET=TEST
511,512
*BEAM SECTION,SECTION=I,TEMPERATURE=GRADIENTS,MATERIAL=ELAS,ELSET=TEST6
1.,2.,2.,2.,.5,.5,.5 
0.,0.,-1.
**
** TEST 7: B31OSH
**
*NODE,NSET=NALL
601, 10.0,  4.0
602, 11.0,  4.0
603, 12.0,  4.0
611, 10.0,  3.0, .0
612, 11.5,  2.9, .1
613, 13.0,  2.7, .2
*NSET,NSET=LEFT
601,611
*NSET,NSET=MID
612, 
*NSET,NSET=RIGHT
613, 
*NSET,NSET=TIP
603,613
*ELEMENT,TYPE=B31OSH,ELSET=TEST7
601, 601,602
602, 602,603
611, 611,612
612, 612,613
*ELSET,ELSET=REFER
601,602
*ELSET,ELSET=TEST
611,612
*BEAM SECTION,SECTION=I,TEMPERATURE=GRADIENTS,MATERIAL=ELAS,ELSET=TEST7
1.,2.,2.,2.,.5,.5,.5 
0.,0.,-1.
**
** TEST 8: B32H
**
*NODE,NSET=NALL
701, 0.0,  7.0
702, 1.0,  7.0
703, 2.0,  7.0
711, 0.0,  6.0, .0
712, 1.5,  5.9, .1
713, 3.0,  5.7, .2
*NSET,NSET=LEFT
701,711
*NSET,NSET=MID
712, 
*NSET,NSET=RIGHT
713, 
*NSET,NSET=TIP
703,713
*ELEMENT,TYPE=B32H,ELSET=TEST8
701, 701,702,703
711, 711,712,713
*ELSET,ELSET=REFER
701, 
*ELSET,ELSET=TEST
711, 
*BEAM SECTION,SECTION=RECT,TEMPERATURE=GRADIENTS,MATERIAL=ELAS,ELSET=TEST8
1.,2.
0.,0.,-1.
***
** TEST 9: B32OS
**
*NODE,NSET=NALL
801, 5.0,  7.0
802, 6.0,  7.0
803, 7.0,  7.0
811, 5.0,  6.0, .0
812, 6.5,  5.9, .1
813, 8.0,  5.7, .2
*NSET,NSET=LEFT
801,811
*NSET,NSET=MID
812, 
*NSET,NSET=RIGHT
813, 
*NSET,NSET=TIP
803,813
*ELEMENT,TYPE=B32OS,ELSET=TEST9
801, 801,802,803
811, 811,812,813
*ELSET,ELSET=REFER
801, 
*ELSET,ELSET=TEST
811, 
*BEAM SECTION,SECTION=I,TEMPERATURE=GRADIENTS,MATERIAL=ELAS,ELSET=TEST9
1.,2.,2.,2.,.5,.5,.5 
0.,0.,-1.
***
** TEST 10: B32OSH
**
*NODE,NSET=NALL
901, 5.0,  7.0
902, 6.0,  7.0
903, 7.0,  7.0
911, 5.0,  6.0, .0
912, 6.5,  5.9, .1
913, 8.0,  5.7, .2
*NSET,NSET=LEFT
901,911
*NSET,NSET=MID
912, 
*NSET,NSET=RIGHT
913, 
*NSET,NSET=TIP
903,913
*ELEMENT,TYPE=B32OSH,ELSET=TEST10
901, 901,902,903
911, 911,912,913
*ELSET,ELSET=REFER
901, 
*ELSET,ELSET=TEST
911, 
*BEAM SECTION,SECTION=I,TEMPERATURE=GRADIENTS,MATERIAL=ELAS,ELSET=TEST10
1.,2.,2.,2.,.5,.5,.5 
0.,0.,-1.
**
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*NORMAL
TEST,LEFT,0.,1.,0.
TEST,MID,.086,.9963,0.
TEST,RIGHT,.11723,.9931,0.
*MATERIAL,NAME=ELAS
*ELASTIC
1000., .25
*EXPANSION
1.E-3, 
*DENSITY
1.E-3, 
*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
LEFT,1,7
RIGHT,1,1,-1.0
RIGHT,2,2, 0.3
RIGHT,3,3,-0.2
RIGHT,4,5, 0.0
RIGHT,6,6,.11750075
MID,1,1,-0.5
MID,2,2, 0.1
MID,3,3,-0.1
MID,4,5, 0.0
MID,6,6,.086105945
*NODE PRINT,NSET=LEFT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*NODE FILE,NSET=LEFT,FREQ=999
RF, 
*EL FILE,POSITION=CENTROID,FREQ=999
E, S
*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
LEFT,1,7
*NODE PRINT,NSET=LEFT,FREQ=999
RF, 
*END STEP
**
**
*STEP,NLGEOM
 STEP #4: APPLY DLOAD TYPE P1
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST,P1,15.
REFER,P1,15.
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: APPLY DLOAD TYPE P2
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST,P2,15.
REFER,P2,15.
*END STEP
**
**
*STEP,NLGEOM
 STEP #6: APPLY DLOAD TYPE GRAV
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST, GRAV,150.,0.,1.,0.
REFER,GRAV,150.,0.,1.,0.
*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