*HEADING
 FREQUENCY PROCEDURES - /STANDARD VERIFICATION, 3D CONNECTORS
 TEST 9-DIGIT LABELS
**
*PREPRINT, ECHO=YES,HIST=YES,MODEL=YES
*NODE
900000001,0.0,0.0
900000011,0.0,10.0
900000100,-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
900000001,900000011,1
*ELEMENT, TYPE=B33, ELSET=COLUMN
900000001,900000001,900000002
*ELGEN, ELSET=COLUMN
900000001,10,1,1
**
** CONNECTOR ELEMENTS
**
*ELEMENT, TYPE=CONN3D2, ELSET=AXIAL
900001001,900000100,900000006
*ELEMENT, TYPE=CONN3D2, ELSET=BASE
900001002,,900000001
*ELEMENT, TYPE=CONN3D2, ELSET=ACTUATOR 
900001003,,900000011
*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
 CARTESIAN
 GLOBAL,
*CONNECTOR SECTION, ELSET=CART-ROT, BEHAVIOR= CB_ROT
 CARTESIAN, CARDAN
 GLOBAL,
*ORIENTATION, NAME=GLOBAL
 1.0,0.0,0.0,0.0,1.0,0.0
**
** 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,
 ACTUATOR,4,
 ACTUATOR,5,
 BASE,1,
 BASE,2,
 BASE,3,
 BASE,4,
 BASE,5,
*BOUNDARY
900000100,1,6
 NCOL,3,5
**
**
**
*STEP
 STEP-1: OBTAIN FIRST 5 NATURAL MODES
**
*FREQUENCY, EIGENSOLVER=SUBSPACE, SIM=NO
 5,
**
*OUTPUT, FIELD
*NODE OUTPUT
 UT,RF
*ELEMENT OUTPUT, ELSET=COLUMN
 S,E 
*ELEMENT OUTPUT, ELSET=CONNECTORS
 CTF,CEF,CU
*NODE PRINT
 U1,U2,UR3
*EL PRINT, 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
900000011,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
 U1,U2,UR3
*EL PRINT, ELSET=CONNECTORS
 CTF
 CEF
 CRF
 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
 U1,U2,UR3
*EL PRINT, ELSET=CONNECTORS
 CTF
 CEF
 CRF
 CU
*END STEP