*HEADING
CONTINUUM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
WITH REBAR
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 39: CPE4
**
*NODE,NSET=NALL
1, -5.0, 0.0
2, -4.0, 0.0
3, -4.0, -1.0
4, -5.0, -1.0
5, -3.0, 0.0
6, -3.0, -1.0
7, -5.0, -3.0
8, -3.5, -3.1
9, -3.5, -4.1
10, -5.0, -4.0
11, -2.0, -3.2
12, -2.0, -4.2
*NSET,NSET=LEFT
1,4,7,10
*NSET,NSET=MID
8,9
*NSET,NSET=RIGHT
11,12
*NSET,NSET=TIP
5,6,11,12
*ELEMENT,TYPE=CPE4,ELSET=TEST39
1, 1,4,3,2
2, 2,3,6,5
3, 7,10,9,8
4, 8,9,12,11
*ELSET,ELSET=REFER
1,2
*ELSET,ELSET=TEST
3,4
**
** TEST 40: CPE8
**
*NODE,NSET=NALL
401, 5.0, 0.0
402, 6.0, 0.0
403, 6.0, -1.0
404, 5.0, -1.0
405, 7.0, 0.0
406, 7.0, -1.0
407, 7.0, -0.5
408, 5.0, -0.5
409, 5.0, -3.0
410, 6.5, -3.1
411, 6.5, -4.1
412, 5.0, -4.0
413, 8.0, -3.2
414, 8.0, -4.2
415, 8.0, -3.7
416, 5.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
*ELEMENT,TYPE=CPE8,ELSET=TEST40
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 41: CPS4
**
*NODE,NSET=NALL
1001, 15.0, 0.0
1002, 16.0, 0.0
1003, 16.0, -1.0
1004, 15.0, -1.0
1005, 17.0, 0.0
1006, 17.0, -1.0
1007, 15.0, -3.0
1008, 16.5, -3.1
1009, 16.5, -4.1
1010, 15.0, -4.0
1011, 18.0, -3.2
1012, 18.0, -4.2
*NSET,NSET=LEFT
1001,1004,1007,1010
*NSET,NSET=MID
1008,1009
*NSET,NSET=RIGHT
1011,1012
*NSET,NSET=TIP
1005,1006,1011,1012
*ELEMENT,TYPE=CPS4,ELSET=TEST41
1001, 1001,1004,1003,1002
1002, 1002,1003,1006,1005
1003, 1007,1010,1009,1008
1004, 1008,1009,1012,1011
*ELSET,ELSET=REFER
1001,1002
*ELSET,ELSET=TEST
1003,1004
**
** TEST 42: CPS8
**
*NODE,NSET=NALL
4001, 25.0, 0.0
4002, 26.0, 0.0
4003, 26.0, -1.0
4004, 25.0, -1.0
4005, 27.0, 0.0
4006, 27.0, -1.0
4007, 27.0, -0.5
4008, 25.0, -0.5
4009, 25.0, -3.0
4010, 26.5, -3.1
4011, 26.5, -4.1
4012, 25.0, -4.0
4013, 28.0, -3.2
4014, 28.0, -4.2
4015, 28.0, -3.7
4016, 25.0, -3.5
*NSET,NSET=LEFT
4001,4004,4008,4009,4012,4016
*NSET,NSET=MID
4010,4011
*NSET,NSET=RIGHT
4013,4014,4015
*NSET,NSET=TIP
4005,4006,4007,4013,4014,4015
*ELEMENT,TYPE=CPS8,ELSET=TEST42
4001, 4001,4004,4006,4005,4008,4003,4007,4002
4002, 4009,4012,4014,4013,4016,4011,4015,4010
*ELSET,ELSET=REFER
4001,
*ELSET,ELSET=TEST
4002,
**
** Axisymmetric Elements
**
*Node, Nset=NALL
6001, 35., 0.
6002, 36., 0.
6003, 36., -1.
6004, 35., -1.
6005, 37., 0.
6006, 37., -1.
6007, 35., -3.
6008, 36.5, -3.1
6009, 36.5, -4.1
6010, 35., -4.
6011, 38., -3.2
6012, 38., -4.2
10001, 45., 0.
10002, 46., 0.
10003, 46., -1.
10004, 45., -1.
10005, 47., 0.
10006, 47., -1.
10007, 47., -0.5
10008, 45., -0.5
10009, 45., -3.
10010, 46.5, -3.1
10011, 46.5, -4.1
10012, 45., -4.
10013, 48., -3.2
10014, 48., -4.2
10015, 48., -3.7
10016, 45., -3.5
11016, 35.4, -1.
11017, 35.4, 0.
11018, 36.4, -1.
11019, 36.4, 0.
12016, 45.8, -1.
12017, 45.8, -0.5
12018, 45.8, 0.
13016, 35.6, -4.04
13017, 35.6, -3.04
13018, 37.1, -4.14
13019, 37.1, -3.14
14016, 46.2, -4.08
14017, 46.2, -3.58
14018, 46.2, -3.08
*Element, type=CAX4
6001, 6001, 6004, 6003, 6002
6002, 6002, 6003, 6006, 6005
6003, 6007, 6010, 6009, 6008
6004, 6008, 6009, 6012, 6011
*Element, type=CAX8
10001, 10001, 10004, 10006, 10005, 10008, 10003, 10007, 10002
10002, 10009, 10012, 10014, 10013, 10016, 10011, 10015, 10010
**
*Elset,Elset=ONE
6001, 6002
*Elset,Elset=TWO
10001,
*Elset,Elset=THREE
6003, 6004
*Elset,Elset=FOUR
10002,
**
*ELEMENT,TYPE=SFMAX1,ELSET=SFMAX1-1
11002,11016,11017
11003,11018,11019
*SURFACE SECTION,ELSET=SFMAX1-1
*REBAR LAYER
REBAR1,1.000000,0.300000,0.000000,RBMAT1,-30.000000,1
*EMBEDDED ELEMENT, HOST ELSET=ONE
11002,
11003,
*ELEMENT,TYPE=SFMAX2,ELSET=SFMAX2-2
12002,12016,12017,12018
*SURFACE SECTION,ELSET=SFMAX2-2
*REBAR LAYER
REBAR1,1.000000,0.300000,0.000000,RBMAT1,-30.000000,1
*EMBEDDED ELEMENT, HOST ELSET=TWO
12002,
*ELEMENT,TYPE=SFMAX1,ELSET=SFMAX1-3
13002,13016,13017
13003,13018,13019
*SURFACE SECTION,ELSET=SFMAX1-3
*REBAR LAYER
REBAR2,1.000000,0.300000,0.000000,RBMAT2,-30.000000,1
*EMBEDDED ELEMENT, HOST ELSET=THREE
13002,
13003,
*ELEMENT,TYPE=SFMAX2,ELSET=SFMAX2-4
14002,14016,14017,14018
*SURFACE SECTION,ELSET=SFMAX2-4
*REBAR LAYER
REBAR2,1.000000,0.300000,0.000000,RBMAT2,-30.000000,1
*EMBEDDED ELEMENT, HOST ELSET=FOUR
14002,
**
**
*NSET,NSET=LEFT
6001,6004,6007,6010
*NSET,NSET=MID
6008,6009
*NSET,NSET=RIGHT
6011,6012
*NSET,NSET=TIP
6005,6006,6011,6012
*NSET,NSET=LEFT
10001,10004,10008,10009,10012,10016
*NSET,NSET=MID
10010,10011
*NSET,NSET=RIGHT
10013,10014,10015
*NSET,NSET=TIP
10005,10006,10007,10013,10014,10015
**
**
*Elset, elset=EALL_CAX
6001, 6002, 6003, 6004, 10001, 10002
*Solid Section, elset=EALL_CAX, material=ELAS
1.,
**
** REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*SOLID SECTION,MATERIAL=ELAS,ELSET=EALL
1.0,
*ELSET,ELSET=RTEST
TEST,6003,6004,10002,SFMAX1-3,SFMAX2-4
*MATERIAL,NAME=ELAS
*ELASTIC
10.e6, .25
*EXPANSION
1.E-6,
*DENSITY
1.E-4,
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,70.
*AMPLITUDE,NAME=RAMP
0.,0., 1.,1.
*REBAR,ELEMENT=CONTINUUM,MATERIAL=RBMAT1,NAME=REBAR1
REFER,1.0,0.3,-30,0.4,1
*REBAR,ELEMENT=CONTINUUM,MATERIAL=RBMAT2,NAME=REBAR2
TEST,1.0,0.3,-30,0.4,1
*MATERIAL,NAME=RBMAT1
*ELASTIC
30.E6,
*MATERIAL,NAME=RBMAT2
*ELASTIC
30.E6,
**
**
*FILE FORMAT,ZERO INCREMENT
*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
RIGHT,1,1,-1.0
RIGHT,2,2, 0.2
MID,1,1,-0.5
MID,2,2, 0.1
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S,E
*ELEMENT OUTPUT, REBAR
S,E
*NODE PRINT,FREQ=999
RF,
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE FILE,FREQ=999
RF,
*EL FILE,FREQ=999
E, S
*EL FILE, REBAR, FREQ=999
E, S
*END STEP
**
**
*STEP
STEP #2: REMOVE TEST ELEMENTS
*STATIC
1.,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
RTEST,
*END STEP
**
**
*STEP
STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*STATIC
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
RTEST,
*BOUNDARY,OP=NEW
LEFT,1,2
*EL PRINT,REBAR,FREQ=999
S,E
*END STEP
**
**
*STEP
STEP #4: APPLY LOAD TO BOTH ELEMENTS
*STATIC
0.1,1.
*DLOAD
TEST,P4,-0.25
REFER,P4,-0.25
6003,P4,-0.25
6004,P4,-0.25
6001,P4,-0.25
6002,P4,-0.25
10001,P4,-0.25
10002,P4,-0.25
*EL PRINT,REBAR,FREQ=999
S,E
*EL PRINT,FREQ=999
S,E
*END STEP
**
**
*STEP
STEP #5: REMOVE TEST ELEMENT
*STATIC
.1,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
RTEST,
*END STEP
**
**
*STEP
STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
.1,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
RTEST,
*EL PRINT,REBAR,FREQ=999
S,E
*EL PRINT,FREQ=999
S,E
*END STEP
**
**
*STEP
STEP #7: APPLY THERMAL LOAD
*STATIC
.1,1.0
*BOUNDARY,OP=NEW,FIXED
TIP,1
LEFT,1,2
*TEMPERATURE
NALL,100.
*EL PRINT,FREQ=999
S,THE,TEMP
*EL PRINT,REBAR,FREQ=999
S,THE,TEMP
*EL FILE,FREQ=999
S,THE,TEMP
*EL FILE,REBAR,FREQ=999
S,E
*ENERGY FILE,FREQ=999
*END STEP
**
**
*STEP
STEP #8: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos
32, ,
*BOUNDARY,OP=NEW
LEFT,1,2
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*MODAL FILE
*OUTPUT,FIELD,VAR=PRESELECT
*OUTPUT,HISTORY
*END STEP
**
**
*STEP,PERTURBATION
STEP #9: MODAL DYNAMIC
*MODAL DYNAMIC
.0004,.004
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF,
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE FILE,FREQ=999
RF,
*EL FILE,FREQ=999
E,S
*EL FILE,REBAR,FREQ=999
E,S
*END STEP
**
**
*STEP,PERTURBATION
STEP #10: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
1000., 10000., 5
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF,
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE FILE,FREQ=999
RF,
*EL FILE,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
STEP #11: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
1000., 10000., 5
*BOUNDARY,OP=NEW
LEFT,1,2
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF,
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE FILE,FREQ=999
RF,
*EL FILE,FREQ=0
*output,field
*elementoutput
inten,
*elementoutput,rebar
inten,
*END STEP
**
**
*STEP
STEP #12: DYNAMICS
*DYNAMIC,DIRECT
.0004, .004
*BOUNDARY,OP=NEW
LEFT,1,2
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE PRINT,FREQ=999
RF,
*CLOAD,AMP=RAMP
TIP,2,0.2
*ENERGY FILE,FREQ=999
*END STEP