*HEADING CONNECTOR FRICTION BEHAVIOR IN LINEAR DYNAMICS PROCEDURES /STANDARD VERIFICATION, 3D CONNECTORS ** *NODE, NSET=NALL 1,0.0,0.0, 2,1.0,0.0, 3,2.0,0.0, 4,3.0,0.0, *NSET, NSET=FIXED 1, *PARAMETER K=100.0 ** ** CONNECTOR ELEMENTS ** *ELEMENT, TYPE=CONN3D2, ELSET=FR_ROT 1,1,2 *ELEMENT, TYPE=CONN3D2, ELSET=NOFR_ROT 2,2,3 *ELSET, ELSET=ROT FR_ROT,NOFR_ROT *ELEMENT, TYPE=CONN3D2, ELSET=ACTUATOR 3,3,4 *ELSET, ELSET=CONNECTORS ROT,ACTUATOR ** ** ROTARY INERTIA ** *ELEMENT, TYPE=ROTARYI, ELSET=ROTARY 1001,1 1002,2 1003,3 1004,4 ** ** SECTION DEFINITIONS ** *ROTARY INERTIA, ELSET=ROTARY ,,1.0, *CONNECTOR SECTION, ELSET=FR_ROT, BEHAVIOR=CB_FR_ROT CARDAN GLOBAL, *CONNECTOR SECTION, ELSET=NOFR_ROT, BEHAVIOR=CB_NOFR_ROT CARDAN GLOBAL, *CONNECTOR SECTION, ELSET=ACTUATOR,behavior=actuator CARDAN GLOBAL, *ORIENTATION, NAME=GLOBAL 1.0,0.0,0.0,0.0,1.0,0.0 ** ** CONNECTOR BEHAVIOR DEFINITIONS ** *CONNECTOR BEHAVIOR, NAME=CB_NOFR_ROT *CONNECTOR ELASTICITY, COMPONENT=6 , ** *CONNECTOR BEHAVIOR, NAME=CB_FR_ROT *CONNECTOR ELASTICITY, COMPONENT=6 , *CONNECTOR FRICTION, COMponent=6, stick stiffness=1.e6 100. **1.0,0.0 **91.0,0.1 *friction 0.1, *CONNECTOR BEHAVIOR, NAME=ACTUATOR *CONNECTOR ELASTICITY, COMPONENT=6 , ** ** ** Boundary CONSTRAINTS (BASE STATE) ** *BOUNDARY NALL,1,3 NALL,4,5 FIXED,6,6 ** ** **------------------------- *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 UR3,RM3 *EL PRINT, ELSET=ROT CTM3,CEM3,CRM3 CUR3, *EL PRINT, ELSET=ACTUATOR CTM3,CEM3,CRM3 CUR3, *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,6,-1.0 *CONNECTOR LOAD, LOAD CASE=2 ACTUATOR,6,-1.0 ** *OUTPUT, FIELD *NODE OUTPUT U,RF *OUTPUT, HISTORY *ELEMENT OUTPUT, ELSET=CONNECTORS CTM3,CEM3,CVM3,CRM3 CUR3, *NODE PRINT, FREQUENCY=1 UR3,VR3,AR3,PUR3 CM3,RM3 *EL PRINT, ELSET=ROT, FREQUENCY=1 CTM3,CEM3,CVM3,CRM3 PHCTM3,PHCEM3,PHCVM3,PHCRM3 CUR3,PHCUR3 *EL PRINT, ELSET=ACTUATOR, FREQUENCY=1 CTM3,CEM3,CVM3,CRM3 PHCTM3,PHCEM3,PHCVM3,PHCRM3 CUR3,PHCUR3 *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, STICK ** *STATIC 0.0125,1.0,0.0125,0.0125 *CLOAD 2,6,5.0 ** *OUTPUT, FIELD *NODE OUTPUT U,RF *ELEMENT OUTPUT NFORC, *OUTPUT, HISTORY *ELEMENT OUTPUT, ELSET=CONNECTORS CTM3,CEM3,CVM3,CSM3,CRM3 CUR3,CPR3 *NODE PRINT, FREQUENCY=1 UR3,CM3,RM3 *EL PRINT, ELSET=ROT, FREQUENCY=1 CTM3,CEM3,CVM3,CRM3, CUR3,CPR3,CSM3,CNM3 *EL PRINT, ELSET=ACTUATOR, FREQUENCY=1 CTM3,CEM3,CVM3,CRM3 CUR3,CPR3 *NODE FILE U,RF *EL FILE, ELSET=CONNECTORS CTF,CEF,CVF,CSF,CRF,CP,CU *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 UR3,RM3 *EL PRINT, ELSET=ROT CTM3,CEM3,CRM3 CUR3, *EL PRINT, ELSET=ACTUATOR CTM3,CEM3,CRM3 CUR3, *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,6,1.0 *CONNECTOR LOAD, LOAD CASE=2 ACTUATOR,6,1.0 ** *OUTPUT, FIELD *NODE OUTPUT U,RF *OUTPUT, HISTORY *ELEMENT OUTPUT, ELSET=CONNECTORS CTM3,CEM3,CVM3,CRM3 CUR3, *NODE PRINT, FREQUENCY=1 UR3,VR3,AR3,PUR3 CM3,RM3 *EL PRINT, ELSET=ROT, FREQUENCY=1 CTM3,CEM3,CVM3,CRM3 PHCTM3,PHCEM3,PHCVM3,PHCRM3 CUR3,PHCUR3 *EL PRINT, ELSET=ACTUATOR, FREQUENCY=1 CTM3,CEM3,CVM3,CRM3 PHCTM3,PHCEM3,PHCVM3,PHCRM3 CUR3,PHCUR3 *NODE FILE U,RF *EL FILE, ELSET=CONNECTORS CTF,CEF,CVF,CRF PHCTF,PHCEF,PHCVF,PHCRF CU,PHCU *END STEP **------------------------- *STEP, NLGEOM STEP-6: CHANGE THE BASE STATE - ACTIVATE FRICTION, SLIP ** *STATIC 0.0125,1.0,0.0125,0.0125 *CLOAD 2,6,15.0 ** *OUTPUT, FIELD *NODE OUTPUT U,RF *ELEMENT OUTPUT NFORC, *OUTPUT, HISTORY *ELEMENT OUTPUT, ELSET=CONNECTORS CTM3,CEM3,CVM3,CSM3,CRM3 CUR3,CPR3 *NODE PRINT, FREQUENCY=1 UR3,CM3,RM3 *EL PRINT, ELSET=ROT, FREQUENCY=1 CTM3,CEM3,CVM3,CRM3, CUR3,CPR3,CSM3,CNM3 *EL PRINT, ELSET=ACTUATOR, FREQUENCY=1 CTM3,CEM3,CVM3,CRM3 CUR3,CPR3 *NODE FILE U,RF *EL FILE, ELSET=CONNECTORS CTF,CEF,CVF,CSF,CRF,CP,CU *END STEP ** ** REPEAT ALL STEPS ** **------------------------- *STEP STEP-7: OBTAIN NATURAL MODES OF THE 2-D.O.F. SYSTEM ** *FREQUENCY, EIGENSOLVER=SUBSPACE, SIM=NO 2, ** *OUTPUT, FIELD *NODE OUTPUT U,RF *NODE PRINT UR3,RM3 *EL PRINT, ELSET=ROT CTM3,CEM3,CRM3 CUR3, *EL PRINT, ELSET=ACTUATOR CTM3,CEM3,CRM3 CUR3, *NODE FILE U,RF *EL FILE, ELSET=CONNECTORS CTF,CEF,CRF,CU *END STEP ** ** **------------------------- *STEP STEP-8: 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,6,1.0 *CONNECTOR LOAD, LOAD CASE=2 ACTUATOR,6,1.0 ** *OUTPUT, FIELD *NODE OUTPUT U,RF *OUTPUT, HISTORY *ELEMENT OUTPUT, ELSET=CONNECTORS CTM3,CEM3,CVM3,CRM3 CUR3, *NODE PRINT, FREQUENCY=1 UR3,VR3,AR3,PUR3 CM3,RM3 *EL PRINT, ELSET=ROT, FREQUENCY=1 CTM3,CEM3,CVM3,CRM3 PHCTM3,PHCEM3,PHCVM3,PHCRM3 CUR3,PHCUR3 *EL PRINT, ELSET=ACTUATOR, FREQUENCY=1 CTM3,CEM3,CVM3,CRM3 PHCTM3,PHCEM3,PHCVM3,PHCRM3 CUR3,PHCUR3 *NODE FILE U,RF *EL FILE, ELSET=CONNECTORS CTF,CEF,CVF,CRF PHCTF,PHCEF,PHCVF,PHCRF CU,PHCU *END STEP