*HEADING CONNECTOR FRICTION BEHAVIOR IN LINEAR DYNAMICS PROCEDURES PERTURBATION/LOAD CASE ANALYSIS /STANDARD VERIFICATION, 2D CONNECTORS ** *NODE, NSET=NALL 1,0.0,0.0, 2,1.0,0.0, 3,2.0,0.0, 4,4.0,0.0, *NSET, NSET=FIXED 1,4 *PARAMETER K=100.0 ** ** CONNECTOR ELEMENTS ** *ELEMENT, TYPE=CONN2D2, ELSET=FR_AXIAL 1,1,2 *ELEMENT, TYPE=CONN2D2, ELSET=NOFR_AXIAL 2,2,3 *ELSET, ELSET=AXIAL FR_AXIAL,NOFR_AXIAL *ELEMENT, TYPE=CONN2D2, ELSET=ACTUATOR 3,3,4 *ELSET, ELSET=CONNECTORS AXIAL,ACTUATOR ** ** MASS ELEMENTS ** *ELEMENT, TYPE=MASS, ELSET=MASS 1001,1 1002,2 1003,3 1004,4 ** ** SECTION DEFINITIONS ** *MASS, ELSET=MASS 1.0, *CONNECTOR SECTION, ELSET=FR_AXIAL, BEHAVIOR=CB_FR_AXIAL AXIAL *CONNECTOR SECTION, ELSET=NOFR_AXIAL, BEHAVIOR=CB_NOFR_AXIAL AXIAL *CONNECTOR SECTION, ELSET=ACTUATOR AXIAL ** ** CONNECTOR BEHAVIOR DEFINITIONS ** *CONNECTOR BEHAVIOR, NAME=CB_NOFR_AXIAL *CONNECTOR ELASTICITY, COMPONENT=1 , ** *CONNECTOR BEHAVIOR, NAME=CB_FR_AXIAL *CONNECTOR ELASTICITY, COMPONENT=1 , *CONNECTOR FRICTION, COMPONENT=1, STICK STIFFNESS=1.e6 100. **1.0 ,0.0 **91.0,10.0 *friction 0.1, ** ** BOUNDARY CONSTRAINTS (BASE STATE) ** *BOUNDARY NALL,2,3 FIXED,1, ** ** **------------------------ *STEP STEP-1: OBTAIN NATURAL MODES OF THE 2-D.O.F. SYSTEM ** *FREQUENCY, EIGENSOLVER=SUBSPACE, SIM=NO 2, ** *OUTPUT, FIELD *NODE OUTPUT U,RF *NODE PRINT U1,RF1 *EL PRINT, ELSET=AXIAL CTF1,CEF1,CRF1 CU1, *EL PRINT, ELSET=ACTUATOR CTF1,CEF1,CRF1 CU1, *NODE FILE U,RF *EL FILE, ELSET=CONNECTORS CTF,CEF,CRF,CU *END STEP ** ** **------------------------ *STEP STEP-2: OBTAIN STEADY STATE DYNAMIC RESPONSE IN TERMS OF PHYSICAL DEGREES OF FREEDOM - DIRECT ** USE CONNECTOR LOAD ** *STEADY STATE DYNAMICS, DIRECT, INTERVAL=EIGENFREQUENCY 12, 60 ** *CONNECTOR LOAD, LOAD CASE=1 ACTUATOR,1,-1.0 *CONNECTOR LOAD, LOAD CASE=2 ACTUATOR,1,-1.0 ** *OUTPUT, FIELD *NODE OUTPUT U,RF *OUTPUT, HISTORY *ELEMENT OUTPUT, ELSET=CONNECTORS CTF1,CEF1,CVF1,CRF1 CU1, *NODE PRINT, FREQUENCY=1 U1,V1,A1,PU1 CF1,RF1 *EL PRINT, ELSET=AXIAL, FREQUENCY=1 CTF1,CEF1,CVF1,CRF1, PHCTF1,PHCEF1,PHCVF1,PHCRF1 CU1,PHCU1 *EL PRINT, ELSET=ACTUATOR, FREQUENCY=1 CTF1,CEF1,CVF1,CRF1 PHCTF1,PHCEF1,PHCVF1,PHCRF1 CU1,PHCU1 *NODE FILE U,RF *EL FILE, ELSET=CONNECTORS CTF,CEF,CVF,CRF PHCTF,PHCEF,PHCVF,PHCRF CU,PHCU *END STEP ** ** **------------------------ *STEP, NLGEOM STEP-3: CHANGE THE BASE STATE - ACTIVATE FRICTION ** *STATIC 0.05,1.0 *CLOAD 2,1,5. ** *OUTPUT, FIELD *NODE OUTPUT U,RF *ELEMENT OUTPUT NFORC,ELEN *OUTPUT, HISTORY *ELEMENT OUTPUT, ELSET=CONNECTORS CTF1,CEF1,CVF1,CSF1,CRF1 CU1,CP1 *NODE PRINT, FREQUENCY=1 U1,CF1,RF1 *EL PRINT, ELSET=AXIAL, FREQUENCY=1 CTF1,CEF1,CVF1,CSF1,CRF1, CU1,CP1 *EL PRINT, ELSET=ACTUATOR, FREQUENCY=1 CTF1,CEF1,CVF1,CRF1 CU1,CP1 *NODE FILE U,RF *EL FILE, ELSET=CONNECTORS CTF,CEF,CVF,CSF,CRF,CP,CU, NFORC, *END STEP ** ** REPEAT ALL STEPS ** **------------------------ *STEP STEP-4: OBTAIN NATURAL MODES OF THE 2-D.O.F. SYSTEM ** *FREQUENCY, EIGENSOLVER=SUBSPACE, SIM=NO 2, ** *OUTPUT, FIELD *NODE OUTPUT U,RF *NODE PRINT U1,RF1 *EL PRINT, ELSET=AXIAL CTF1,CEF1,CRF1 CU1, *EL PRINT, ELSET=ACTUATOR CTF1,CEF1,CRF1 CU1, *NODE FILE U,RF *EL FILE, ELSET=CONNECTORS CTF,CEF,CRF,CU *END STEP ** ** **------------------------ *STEP STEP-5: OBTAIN STEADY STATE DYNAMIC RESPONSE IN TERMS OF PHYSICAL DEGREES OF FREEDOM - DIRECT ** USE CONNECTOR LOAD ** *STEADY STATE DYNAMICS, DIRECT, INTERVAL=EIGENFREQUENCY 12, 60 ** *CONNECTOR LOAD, LOAD CASE=1 ACTUATOR,1,1.0 *CONNECTOR LOAD, LOAD CASE=2 ACTUATOR,1,1.0 ** *OUTPUT, FIELD *NODE OUTPUT U,RF *OUTPUT, HISTORY *ELEMENT OUTPUT, ELSET=CONNECTORS CTF1,CEF1,CVF1,CRF1 CU1, *NODE PRINT, FREQUENCY=1 U1,V1,A1,PU1 CF1,RF1 *EL PRINT, ELSET=AXIAL, FREQUENCY=1 CTF1,CEF1,CVF1,CRF1, PHCTF1,PHCEF1,PHCVF1,PHCRF1 CU1,PHCU1 *EL PRINT, ELSET=ACTUATOR, FREQUENCY=1 CTF1,CEF1,CVF1,CRF1 PHCTF1,PHCEF1,PHCVF1,PHCRF1 CU1,PHCU1 *NODE FILE U,RF *EL FILE, ELSET=CONNECTORS CTF,CEF,CVF,CRF PHCTF,PHCEF,PHCVF,PHCRF CU,PHCU *END STEP **------------------------ *STEP,PERTURBATION STEP 6: PERTURBATION ANALYSIS TO TEST LOAD CASES *STATIC *LOAD CASE,NAME=AXIAL-1 *CLOAD 2,1,10.0 *END LOAD CASE *LOAD CASE,NAME=AXIAL-2 *CLOAD 3,1,10.0 *END LOAD CASE *OUTPUT, FIELD *NODE OUTPUT U,RF *NODE PRINT, FREQUENCY=1 U1,V1,A1,PU1 CF1,RF1 *EL PRINT, ELSET=AXIAL, FREQUENCY=1 CTF1,CEF1,CVF1,CRF1, PHCTF1,PHCEF1,PHCVF1,PHCRF1 CU1,PHCU1 *EL PRINT, ELSET=ACTUATOR, FREQUENCY=1 CTF1,CEF1,CVF1,CRF1 PHCTF1,PHCEF1,PHCVF1,PHCRF1 CU1,PHCU1 *END STEP **------------------------------- *STEP, NLGEOM STEP-7: CHANGE THE BASE STATE - LOAD TO SLIP ** *STATIC 0.05,1.0 *CLOAD 2,1,15. ** *OUTPUT, FIELD *NODE OUTPUT U,RF *ELEMENT OUTPUT NFORC,ELEN *OUTPUT, HISTORY *ELEMENT OUTPUT, ELSET=CONNECTORS CTF1,CEF1,CVF1,CSF1,CRF1 CU1,CP1 *NODE PRINT, FREQUENCY=1 U1,CF1,RF1 *EL PRINT, ELSET=AXIAL, FREQUENCY=1 CTF1,CEF1,CVF1,CSF1,CRF1, CU1,CP1 *EL PRINT, ELSET=ACTUATOR, FREQUENCY=1 CTF1,CEF1,CVF1,CRF1 CU1,CP1 *NODE FILE U,RF *EL FILE, ELSET=CONNECTORS CTF,CEF,CVF,CSF,CRF,CP,CU, NFORC, *END STEP ** ** REPEAT ALL STEPS ** **------------------------ *STEP STEP-8: OBTAIN NATURAL MODES OF THE 2-D.O.F. SYSTEM ** *FREQUENCY, EIGENSOLVER=SUBSPACE, SIM=NO 2, ** *OUTPUT, FIELD *NODE OUTPUT U,RF *NODE PRINT U1,RF1 *EL PRINT, ELSET=AXIAL CTF1,CEF1,CRF1 CU1, *EL PRINT, ELSET=ACTUATOR CTF1,CEF1,CRF1 CU1, *NODE FILE U,RF *EL FILE, ELSET=CONNECTORS CTF,CEF,CRF,CU *END STEP ** **------------------------------- *STEP STEP-9: OBTAIN STEADY STATE DYNAMIC RESPONSE IN TERMS OF PHYSICAL DEGREES OF FREEDOM - DIRECT ** USE CONNECTOR LOAD ** *STEADY STATE DYNAMICS, DIRECT, INTERVAL=EIGENFREQUENCY 12, 60 ** *CONNECTOR LOAD, LOAD CASE=1 ACTUATOR,1,1.0 *CONNECTOR LOAD, LOAD CASE=2 ACTUATOR,1,1.0 ** *OUTPUT, FIELD *NODE OUTPUT U,RF *OUTPUT, HISTORY *ELEMENT OUTPUT, ELSET=CONNECTORS CTF1,CEF1,CVF1,CRF1 CU1, *NODE PRINT, FREQUENCY=1 U1,V1,A1,PU1 CF1,RF1 *EL PRINT, ELSET=AXIAL, FREQUENCY=1 CTF1,CEF1,CVF1,CRF1, PHCTF1,PHCEF1,PHCVF1,PHCRF1 CU1,PHCU1 *EL PRINT, ELSET=ACTUATOR, FREQUENCY=1 CTF1,CEF1,CVF1,CRF1 PHCTF1,PHCEF1,PHCVF1,PHCRF1 CU1,PHCU1 *NODE FILE U,RF *EL FILE, ELSET=CONNECTORS CTF,CEF,CVF,CRF PHCTF,PHCEF,PHCVF,PHCRF CU,PHCU *END STEP **----------------------------------------- *STEP,PERTURBATION STEP 10: PERTURBATION ANALYSIS TO TEST LOAD CASES *STATIC *LOAD CASE,NAME=AXIAL-11 *CLOAD 2,1,10.0 *END LOAD CASE *LOAD CASE,NAME=AXIAL-21 *CLOAD 3,1,10.0 *END LOAD CASE *OUTPUT, FIELD *NODE OUTPUT U,RF *NODE PRINT, FREQUENCY=1 U1,V1,A1,PU1 CF1,RF1 *EL PRINT, ELSET=AXIAL, FREQUENCY=1 CTF1,CEF1,CVF1,CRF1, PHCTF1,PHCEF1,PHCVF1,PHCRF1 CU1,PHCU1 *EL PRINT, ELSET=ACTUATOR, FREQUENCY=1 CTF1,CEF1,CVF1,CRF1 PHCTF1,PHCEF1,PHCVF1,PHCRF1 CU1,PHCU1 *END STEP **------------------------