*HEADING
 CONNECTOR FRICTION BEHAVIOR IN LINEAR DYNAMICS PROCEDURES
 /STANDARD VERIFICATION, 3D CONNECTORS
**
*NODE, NSET=NALL
*NSET, NSET=FIXED
*PARAMETER
 C=1.0
 K=100.0
**
** CONNECTOR ELEMENTS
**
*ELEMENT, TYPE=CONN3D2, ELSET=FR_ROT
*ELEMENT, TYPE=CONN3D2, ELSET=NOFR_ROT
*ELSET, ELSET=ROT
 FR_ROT,NOFR_ROT
*ELEMENT, TYPE=CONN3D2, ELSET=ACTUATOR
*ELSET, ELSET=CONNECTORS
 ROT,ACTUATOR
**
** ROTARY INERTIA
**
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY
**
** 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
**
** CONNECTOR BEHAVIOR DEFINITIONS
**
*CONNECTOR BEHAVIOR, NAME=CB_NOFR_ROT
*CONNECTOR ELASTICITY, COMPONENT=6
 <K>,
*CONNECTOR DAMPING, COMPONENT=6
 <C>,
**
*CONNECTOR BEHAVIOR, NAME=CB_FR_ROT
*CONNECTOR ELASTICITY, COMPONENT=6
 <K>,
*CONNECTOR DAMPING, COMPONENT=6
 <C>,
*CONNECTOR FRICTION, COMponent=6
*friction
0.1,
*CONNECTOR BEHAVIOR, NAME=ACTUATOR
*CONNECTOR ELASTICITY, COMPONENT=6
 <K>,
*CONNECTOR DAMPING, COMPONENT=6
 <C>,
**
**
** 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
 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
**
*STATIC
 0.0125,1.0,0.0125,0.0125
*CLOAD
 2,6,10.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
 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