*HEADING C3D8 ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL *PREPRINT,MODEL=YES,HISTORY=NO ** ** *NODE,NSET=NALL 1,0,0,0 2,1,0,0 3,1,1,0 4,0,1,0 5,0,0,1 6,1,0,1 7,1,1,1 8,0,1,1 9,2,0,0 10,2,1,0 11,2,0,1 12,2,1,1 ** 21,3,0,0 22,4,0,0 23,4,1.03,0 24,3,1.03,0 25,3,0,1 26,4,0,1 27,4,1.03,1 28,3,1.03,1 29,5,0,0 30,5,1.03,0 31,5,0,1 32,5,1.03,1 ** 1032, 0., 1., 0.5 1033, 0., 0., 0.5 1034, 1., 0., 0.5 1035, 1., 1., 0.5 1036, 1., 1., 0.5 1037, 1., 0., 0.5 1038, 2., 0., 0.5 1039, 2., 1., 0.5 1040, 3., 1.03, 0.5 1041, 3., 0., 0.5 1042, 4., 0., 0.5 1043, 4., 1.03, 0.5 1044, 4., 1.03, 0.5 1045, 4., 0., 0.5 1046, 5., 0., 0.5 1047, 5., 1.03, 0.5 ** *NSET,NSET=BOT 1,2,9,5,6,11,21,22,29,25,26,31 *NSET,NSET=TIP 4,3,10,8,7,12,24,23,30,28,27,32 *NSET,NSET=DEFBC 24,23,30,28,27,32 *ELEMENT,TYPE=C3D8,ELSET=TEST1 1,1,2,3,4,5,6,7,8 2,2,9,10,3,6,11,12,7 *ELEMENT,TYPE=C3D8,ELSET=TEST2 3,21,22,23,24,25,26,27,28 4,22,29,30,23,26,31,32,27 ** *ELSET,ELSET=ONE 1 *ELSET,ELSET=TWO 2 *ELSET,ELSET=THREE 3 *ELSET,ELSET=FOUR 4 ** *ORIENTATION,NAME=REBAR1,DEFINITION=OFFSET TO NODES 2,4 3 *ELEMENT,TYPE=SFM3D4R,ELSET=SFM3D4R-1 1004,1032,1033,1034,1035 *SURFACE SECTION,ELSET=SFM3D4R-1 *REBAR LAYER,ORIENTATION=REBAR1 REB1A,0.500000,1.000000,0.000000,ELAS,59.261,1 *EMBEDDED ELEMENT,HOST ELSET=ONE 1004, *ELEMENT,TYPE=SFM3D4R,ELSET=SFM3D4R-2 1005,1036,1037,1038,1039 *SURFACE SECTION,ELSET=SFM3D4R-2 *REBAR LAYER,ORIENTATION=REBAR1 REB1B,0.500000,1.000000,0.000000,ELAS,59.261,1 *EMBEDDED ELEMENT,HOST ELSET=TWO 1005, *ELEMENT,TYPE=SFM3D4R,ELSET=SFM3D4R-3 1006,1040,1041,1042,1043 *SURFACE SECTION,ELSET=SFM3D4R-3 *REBAR LAYER,ORIENTATION=REBAR1 REB1C,0.500000,1.000000,0.000000,ELAS,59.261,1 *EMBEDDED ELEMENT,HOST ELSET=THREE 1006, *ELEMENT,TYPE=SFM3D4R,ELSET=SFM3D4R-4 1007,1044,1045,1046,1047 *SURFACE SECTION,ELSET=SFM3D4R-4 *REBAR LAYER,ORIENTATION=REBAR1 REB1D,0.500000,1.000000,0.000000,ELAS,59.261,1 *EMBEDDED ELEMENT,HOST ELSET=FOUR 1007, ** *ORIENTATION,NAME=CYL1,SYSTEM=C 0.,0., 0., 0.,0.,1. 3,0. *ORIENTATION,NAME=CYL2,SYSTEM=C 3.,0., 0., 3.,0.,1. 3,0. ** *ELSET,ELSET=REFER 1,2,1004,1005 *ELSET,ELSET=TEST 3,4,1006,1007 ** ** ** *ELSET,ELSET=EALL TEST,REFER,1004,1005,1006,1007 *SOLID SECTION,MATERIAL=ELAS,ELSET=TEST1,ORIENT=CYL1 *SOLID SECTION,MATERIAL=ELAS,ELSET=TEST2,ORIENT=CYL2 *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, *EL PRINT, REBAR S, E, RBANG, *OUTPUT,FIELD,FREQ=999 *ELEMENT OUTPUT S, *OUTPUT,FIELD,FREQ=999 *NODE OUTPUT RF, *ELEMENT OUTPUT,POSITION=CENTROID E, S, *ELEMENT OUTPUT, REBAR S, E, *OUTPUT,HISTORY,FREQ=999 *NODE OUTPUT, NSET=NALL RF, *ELEMENT OUTPUT MASS *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 *EL PRINT, REBAR S, E, *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, *EL PRINT, REBAR S, E, *OUTPUT,FIELD,FREQ=999 *ELEMENT OUTPUT S, THE, E, TEMP, *END STEP ** ** *STEP,NLGEOM STEP #8: FREQUENCY EXTRACTION *FREQUENCY, eigensolver=lanczos 50, , *BOUNDARY,OP=NEW BOT,1,3 *EL PRINT,FREQ=0 *NODE PRINT,FREQ=0 *OUTPUT,HISTORY,FREQ=0 *NODE OUTPUT,NSET=NALL *END STEP ** ** *STEP,PERTURBATION STEP #9: MODAL DYNAMIC *MODAL DYNAMIC .002,.02 *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 S, *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 50.,250., 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 50.,250., 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, *EL PRINT, REBAR S, E, *OUTPUT,HISTORY,FREQ=999 *NODE OUTPUT, NSET=NALL 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 50.,250., 5 *CLOAD TIP,1,0.1 *NODE PRINT,FREQ=999 RF, *EL PRINT,POSITION=CENTROID,FREQ=999 S, E, *EL PRINT, REBAR S, E, *OUTPUT,HISTORY,FREQ=999 *NODE OUTPUT, NSET=NALL RF, *OUTPUT,HISTORY,FREQ=999 *NODE OUTPUT,NSET=NALL RF, *END STEP ** ** *STEP,NLGEOM STEP #13: DYNAMICS *DYNAMIC .002, .02 *BOUNDARY,OP=NEW BOT,PINNED *EL PRINT,POSITION=CENTROID,FREQ=999 S, E, *EL PRINT, REBAR S, E, *NODE PRINT,FREQ=999 RF, *CLOAD,AMP=RAMP TIP,1,.1 *END STEP