*HEADING
 BASIC VERIFICATION - MPC TYPE=ELBOW with *TRANSFORM
  ELEMENT: ELBOW31, TEST TYPE:   TEST NUMBER:
**
** STEP 1: BENDING LOAD - LINEAR
** STEP 2: AXIAL LOAD - LINEAR
** STEP 3: BENDING LOAD - NLGEOM
** STEP 4: AXIAL LOAD - NLGEOM
**
** CASE: 1 (NO MPC)
*NODE, NSET=ALL
 1, 0.0E+0, 0.0E+0, 0.0E+0
 2, 10.0E+0, 0.0E+0, 0.0E+0
 4, 20.0E+0, 0.0E+0, 0.0E+0
**
*ELEMENT, TYPE=ELBOW31, ELSET=ALL1
 1,   1, 2
 2,   2, 4
**
*BEAM SECTION,  SECTION=ELBOW,  ELSET=ALL1,  MATERIAL=MAT
 5.0E+0, 0.25E+00
 0.E+00, 10.0E+01, 0.000E+00
 5, 20, 6
************************************************************
** CASE: 2 
*NODE, NSET=ALL
 101, 0.0E+0, 0.0E+0, 0.0E+0
 102, 10.0E+0, 0.0E+0, 0.0E+0
 103, 10.0E+0, 0.0E+0, 0.0E+0
 104, 20.0E+0, 0.0E+0, 0.0E+0
**
*ELEMENT, TYPE=ELBOW31, ELSET=SEC101
 101,   101, 102
*ELEMENT, TYPE=ELBOW31, ELSET=SEC102
 102,   103, 104
**
*BEAM SECTION,  SECTION=ELBOW,  ELSET=SEC101,  MATERIAL=MAT
 5.0E+0, 0.25E+00, 0.0E+0
 0.00E+00, 100.0E+00, 0.000E+00
 3, 16, 5
*BEAM SECTION,  SECTION=ELBOW,  ELSET=SEC102,  MATERIAL=MAT
 5.0E+0, 0.25E+00, 0.0E+0
 0.00E+00, 100.0E+00, 0.000E+00
 3, 16, 5
**
*MPC
ELBOW, 102, 103
************************************************************
** CASE: 3 
*NODE, NSET=ALL
 201, 0.0E+0, 0.0E+0, 0.0E+0
 202, 10.0E+0, 0.0E+0, 0.0E+0
 203, 10.0E+0, 0.0E+0, 0.0E+0
 204, 20.0E+0, 0.0E+0, 0.0E+0
**
*ELEMENT, TYPE=ELBOW31, ELSET=SEC201
 201,   201, 202
*ELEMENT, TYPE=ELBOW31, ELSET=SEC202
 202,   203, 204
**
*BEAM SECTION,  SECTION=ELBOW,  ELSET=SEC201,  MATERIAL=MAT
 5.0E+0, 0.25E+00, 0.0E+0
 0.00E+00, 100.0E+00, 0.000E+00
 5, 12, 4
**
*BEAM SECTION,  SECTION=ELBOW,  ELSET=SEC202,  MATERIAL=MAT
 5.0E+0, 0.25E+00, 0.0E+0
 0.00E+00, 0.000,   100.0E+00
 5, 12, 4
**
*MPC
ELBOW, 202, 203
***************************************************************
** CASE: 4
*NODE, NSET=ALL
 301, 0.0E+0,  0.0E+0, 0.0E+0
 302, 10.0E+0, 0.0E+0, 0.0E+0
 303, 10.0E+0, 0.0E+0, 0.0E+0
 304, 20.0E+0, 0.0E+0, 0.0E+0
**
*ELEMENT, TYPE=ELBOW31, ELSET=SEC301
 301,   301, 302
*ELEMENT, TYPE=ELBOW31, ELSET=SEC302
 302,   303, 304
**
*BEAM SECTION,  SECTION=ELBOW,  ELSET=SEC301,  MATERIAL=MAT
 5.0E+0, 0.25E+00, 0.0E+0
 0.00E+00, 100.0E+00, 0.000E+00
 5,  20, 6
**
*BEAM SECTION,  SECTION=ELBOW,  ELSET=SEC302,  MATERIAL=MAT
 5.0E+0, 0.25E+00, 0.0E+0
 0.00E+00, -100.0E+00, 0.000
 5, 20, 6
**
*MPC
ELBOW, 302, 303
**
************************************************************
*MATERIAL,  NAME=MAT
*ELASTIC
 30E+06,  0.3
**
*NSET, NSET=MPC_N
103,202
*TRANSFORM, NSET=MPC_N
** 45 about Y axis, then 60 about (transformed) Z axis
0.35355339,0.866025,-0.35355339,-0.6123724,0.5,0.6123724
**
** COMPLETELY FIX NODE 1, 101, 201
*BOUNDARY
 1,  1,6
 101,  1,6
 201,  1,6
 301, 1,6
 ALL, NODEFORM
**
*RESTART, WRITE, OVERLAY, FREQ=3
***********************************
*STEP, PERTURBATION
 STEP 1: LINEAR PURE BENDING APPLIED
**
*STATIC
**
*CLOAD
 4,   6,  1.0E+06
 104,  6,  1.0E+06
 204,  6,  1.0E+06
 304,  6,  1.0E+06
**
*EL PRINT, FREQ=10
S, 
SF, 
*NODE PRINT, FREQ=10, GLO=YES
U, 
RF, 
*EL FILE,FREQ=10
1,2,3,4,5
S, 
SF, 
*NODE FILE, FREQ=10
U, 
RF, 
*END STEP
****************************
*STEP, PERTURBATION
 STEP 2: LINEAR AXIAL LOAD APPLIED
**
*STATIC
**
*CLOAD, OP=NEW
 4,  1,  2.0E+6
 104,  1,  2.0E+6
 204,  1,  2.0E+6
 304,  1,  2.0E+6
**
**
*EL PRINT, FREQ=10
S, 
SF, 
*NODE PRINT, FREQ=10, GLO=YES
U, 
RF, 
*EL FILE,FREQ=10
1,2,3,4,5
S, 
SF, 
*NODE FILE, FREQ=10
U, 
RF, 
*END STEP
** 
** NON LINEAR ANALYSES
**
*****************************************************
*STEP, NLGEOM, INC=20
 STEP 3: NON LINEAR BENDING LOAD APPLIED
**
*STATIC
 0.1, 1.0
**
*CLOAD, OP=NEW
 4,  6,  1.0E+06
 104,  6,  1.0E+06
 204,  6,  1.0E+06
 304,  6,  1.0E+06
**
**
*EL PRINT, FREQ=10
S, 
SF, 
*NODE PRINT, FREQ=10, GLO=YES
U, 
RF, 
*EL FILE,FREQ=10
1,2,3,4,5
S, 
SF, 
*NODE FILE, FREQ=10
U, 
RF, 
*END STEP
***************************************************
*STEP, NLGEOM, INC=20
 STEP 4: NON LINEAR  AXIAL LOAD APPLIED
**
*STATIC
 0.1, 1.0
**
*CLOAD, OP=NEW
 4,  1,  2.0E+6
 104,  1,  2.0E+6
 204,  1,  2.0E+6
 304,  1,  2.0E+6
**
**
*EL PRINT, FREQ=10
S, 
SF, 
*NODE PRINT, FREQ=10, GLO=YES
U, 
RF, 
*EL FILE,FREQ=10
1,2,3,4,5
S, 
SF, 
*NODE FILE, FREQ=10
U, 
RF, 
*END STEP