*HEADING
 CONNECTOR PLASTICITY/DAMAGE BEHAVIOR IN LINEAR DYNAMICS PROCEDURES 
 /STANDARD VERIFICATION, 3D CONNECTORS
**
*NODE, NSET=NALL
 1,0.0,0.0,
 2,1.0,0.0,
 4,3.0,0.0,
*NSET, NSET=FIXED
 1,4
**
** CONNECTOR ELEMENTS
**
*ELEMENT, TYPE=CONN3D2, ELSET=ROT
 1,1,2
*ELEMENT, TYPE=CONN3D2, ELSET=ACTUATOR 
 3,2,4
*ELSET, ELSET=CONNECTORS
 ROT,ACTUATOR
**
** ROTARY INERTIA
**
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY
 1001,1
 1002,2
 1004,4
**
** SECTION DEFINITIONS
**
*ROTARY INERTIA, ELSET=ROTARY
 1.0,1.0,1.0
*CONNECTOR SECTION, ELSET=ROT, BEHAVIOR=CB_ROT
 CARDAN
 GLOBAL,
*CONNECTOR SECTION, ELSET=ACTUATOR
 CARDAN
 GLOBAL,
*ORIENTATION, NAME=GLOBAL
 1.0,0.0,0.0,0.0,1.0,0.0
3,0
**
** CONNECTOR BEHAVIOR DEFINITIONS
**
*CONNECTOR BEHAVIOR, NAME=CB_ROT
*CONNECTOR ELASTICITY, COMPONENT=6
100,
*CONNECTOR DAMPING, COMPONENT=6
1.0,
*CONNECTOR PLASTICITY, COMPONENT=6
*CONNECTOR HARDENING, TYPE=ISOTROPIC, DEFINITION=TABULAR
100,0.0
110,1.0
115,2.0
120,4.0
***CONNECTOR DAMAGE INITIATION, CRITERION=PLASTIC MOTION, COMPONENT=6
**1.0,
***CONNECTOR DAMAGE EVOLUTION, TYPE=MOTION, SOFTENING=LINEAR
**2.0,
**
** BOUNDARY CONSTRAINTS (BASE STATE)
**
*BOUNDARY
 NALL,1,5
 FIXED,6,6
**
**
*STEP
 STEP-1: OBTAIN NATURAL MODES OF THE 2-D.O.F. SYSTEM
**
*FREQUENCY,norm=mass
 1,,,-1000
**
*OUTPUT, FIELD
*NODE OUTPUT
 U,RF
*NODE PRINT
 UR3,RM3
*EL PRINT, ELSET=ROT
 CTM3,CEM3,CRM3
 CUR3,
 CSLST 
*EL PRINT, ELSET=ACTUATOR
 CTM3,CEM3,CRM3
 CUR3,
*NODE FILE
 U,RF
*EL FILE, ELSET=CONNECTORS
 CTF,CEF,CRF,CU
 CSLST,
*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
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,
 CSLST,
*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,
 CSLST 
*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
 CSLST,
*END STEP
**
**
*STEP, NLGEOM
 STEP-3: CHANGE THE BASE STATE - LOCK CONNECTOR 1
**
*STATIC
 0.05,1.0,,0.1
*CLOAD
 2,6,115.0
**
*OUTPUT, FIELD ,freq=100
*NODE OUTPUT
 U,RF
*ELEMENT OUTPUT
 NFORC,ELEN
*OUTPUT, HISTORY,freq=100
*ELEMENT OUTPUT, ELSET=CONNECTORS
 CTM3,CEM3,CVM3,CRM3
 CUR3,CPR3,CURP3,
 CSLST,
*NODE PRINT, FREQUENCY=100
 UR3,CM3,RM3
*EL PRINT, ELSET=ROT, FREQUENCY=100
 CTM3,CEM3,CVM3,CRM3
 CUR3,CPR3,CURP3,
 CSLST
*EL PRINT, ELSET=ACTUATOR, FREQUENCY=100
 CTM3,CEM3,CVM3,CRM3
 CUR3,CPR3,CURP3,
 *NODE FILE,freq=100
 U,RF
*EL FILE, ELSET=CONNECTORS,freq=100
 CTF,CEF,CVF,CRF,CP,CU
 CSLST,
*END STEP
**
** REPEAT ALL STEPS 
**
*STEP
 STEP-4: OBTAIN NATURAL MODES OF THE 2-D.O.F. SYSTEM
**
*FREQUENCY,norm=mass
 1,,,-1000
**
*OUTPUT, FIELD
*NODE OUTPUT
 U,RF
*NODE PRINT
 UR3,RM3
*EL PRINT, ELSET=ROT
 CTM3,CEM3,CRM3
 CUR3,
 CSLST 
*EL PRINT, ELSET=ACTUATOR
 CTM3,CEM3,CRM3
 CUR3,
*NODE FILE
 U,RF
*EL FILE, ELSET=CONNECTORS
 CTF,CEF,CRF,CU
 CSLST,
*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
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,
 CSLST,
*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
 CSLST 
*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
 CSLST,
*END STEP