*HEADING
 FREQUENCY PROCEDURES - /STANDARD VERIFICATION, 2D CONNECTORS
**
*NODE
 1,0.0,0.0,
 11,0.0,10.0,
 100,-1.0,5.0,
*PARAMETER
 KAXIAL1=1.0E3
 KROT1=1.0E3
 U1=1.0E-3
 UR1=1.0E-3
 KAXIAL2=KAXIAL1/1000
 KROT2=KROT1/1000
 U2=100*U1
 UR2=100*UR1
 F1=KAXIAL1*U1
 M1=KROT1*UR1
 F2=F1+KAXIAL2*(U2-U1)
 M2=M1+KROT2*(UR2-UR1)
 mF1=-F1
 mF2=-F2
 mU1=-U1
 mU2=-U2
 mM1=-M1
 mM2=-M2
 mUR1=-UR1
 mUR2=-UR2
**
*NGEN,NSET=NCOL
 1,11,1
*ELEMENT, TYPE=B23, ELSET=COLUMN
 1,1,2
*ELGEN, ELSET=COLUMN
 1,10,1,1
**
** CONNECTOR ELEMENTS
**
*ELEMENT, TYPE=CONN2D2, ELSET=AXIAL
 1001,100,6
*ELEMENT, TYPE=CONN2D2, ELSET=BASE
 1002,,1
*ELEMENT, TYPE=CONN2D2, ELSET=ACTUATOR 
 1003,,11
*ELSET, ELSET=CART-ROT
 BASE, ACTUATOR
*ELSET, ELSET=CONNECTORS
 AXIAL, BASE, ACTUATOR
**
** SECTION DEFINITIONS
**
*BEAM SECTION, SECTION=RECT, ELSET=COLUMN, MATERIAL=STEEL
 0.1,0.1
 0.0,0.0,-1.0
*CONNECTOR SECTION, ELSET=AXIAL, BEHAVIOR= CB_AXIAL
 AXIAL
*CONNECTOR SECTION, ELSET=CART-ROT, BEHAVIOR= CB_ROT
 CARTESIAN, ROTATION
**
** MATERIAL AND CONNECTOR BEHAVIOR DEFINITIONS
**
*MATERIAL, NAME=STEEL
*ELASTIC
 200.0E9,0.3
*DENSITY
 7800.0,
**
*CONNECTOR BEHAVIOR, NAME= CB_AXIAL
*CONNECTOR ELASTICITY, COMPONENT=1, NONLINEAR
 <mF2>,<mU2>
 <mF1>,<mU1>
 <F1>,<U1>
 <F2>,<U2>
**
*CONNECTOR BEHAVIOR, NAME= CB_ROT
*CONNECTOR ELASTICITY, COMPONENT=6, NONLINEAR
 <mM2>,<mUR2>
 <mM1>,<mUR1>
 <M1>,<UR1>
 <M2>,<UR2>
**
** BOUNDARY CONSTRAINTS (BASE STATE)
**
*CONNECTOR MOTION
 ACTUATOR,1,
 BASE,1,
 BASE,2,
*BOUNDARY
 100,1,2
**
**
**
*STEP
 STEP-1: OBTAIN FIRST 5 NATURAL MODES
**
*FREQUENCY, EIGENSOLVER=SUBSPACE, SIM=NO
 5,
**
*OUTPUT, FIELD
*NODE OUTPUT
 U,RF
*ELEMENT OUTPUT, ELSET=COLUMN
 S,E 
*ELEMENT OUTPUT, ELSET=CONNECTORS
 CTF,CEF,CU 
*NODE PRINT
 U,
*EL PRINT, ELSET=CONNECTORS
 CTF
 CEF
 CRF
 CU
*NODE FILE
 U,RF
*EL FILE, ELSET=CONNECTORS
 CTF,CEF,CRF,CU 
*END STEP
**
** The following steps will test the connector behavior due to
** nonlinear elasticity.
**
*STEP, NLGEOM
 STEP-2: STATIC STEP (NLGEOM)
**
** Static step with NLGEOM to modify the base state
**
*STATIC
 0.1,1.0
*CONNECTOR LOAD, OP=NEW
 ACTUATOR,2,-1.0E5
*CLOAD
 11,6,1.0E4
**
*OUTPUT, FIELD
*NODE OUTPUT
 U,RF
*ELEMENT OUTPUT, ELSET=COLUMN
 S,E 
*ELEMENT OUTPUT, ELSET=CONNECTORS
 CTF,CEF,CU 
*OUTPUT,HISTORY,FREQ=1
*ELEMENT OUTPUT,ELSET=CONNECTORS
 CTF,CEF,CRF,CP,CU
*NODE PRINT
 U,
*EL PRINT, ELSET=CONNECTORS
 CTF
 CEF
 CRF
 CU
*NODE FILE
 U,RF
*EL FILE, ELSET=CONNECTORS
 CTF,CEF,CRF,CP,CU 
*END STEP
**
**
*STEP
 STEP-3: OBTAIN FIRST 5 NATURAL MODES
**
*FREQUENCY, EIGENSOLVER=SUBSPACE, SIM=NO
 5,
**
*OUTPUT, FIELD
*NODE OUTPUT
 U,RF
*ELEMENT OUTPUT, ELSET=COLUMN
 S,E 
*ELEMENT OUTPUT, ELSET=CONNECTORS
 CTF,CEF,CU
*NODE PRINT
 U,
*EL PRINT, ELSET=CONNECTORS
 CTF
 CEF
 CRF
 CU
*NODE FILE
 U,RF
*EL FILE, ELSET=CONNECTORS
 CTF,CEF,CRF,CU 
*END STEP