*HEADING
 Inertia Relief analysis for Automobile Axel
 Modeled as Beam elements
*PREPRINT,ECHO=YES
*NODE,NSET=AXEL
*NSET,NSET=NALL
AXEL
*NSET,NSET=FIXED
*ELEMENT,TYPE=B31,ELSET=AXEL1
*BEAM SECTION,SECTION=CIRC,MATERIAL=STEEL,ELSET=AXEL1
13.0,
*ELEMENT,TYPE=B31,ELSET=AXEL2
*BEAM SECTION,SECTION=CIRC,MATERIAL=STEEL,ELSET=AXEL2
15.0,
1.0, 0.0, 0.0
*ELEMENT,TYPE=B31,ELSET=AXEL3
*BEAM GENERAL SECTION,SECTION=GENERAL,ELSET=AXEL3,DENSITY=7.8e-9
391.0,17000.0,0.0,9400.0,22600.0
0.0,0.0,-1
2.1E+05,8.077E+04
*ELSET,ELSET=EALL
AXEL1,AXEL2,AXEL3
*BOUNDARY
FIXED,1,6
*MATERIAL,NAME=STEEL
*DENSITY
7.8E-9
*ELASTIC
2.1E+05,0.3
*STEP,NLGEOM,PERTURBATION
 INERTIA RELIEF ANALYSIS
*STATIC
*BOUNDARY
FIXED,1,6
*DLOAD
EALL,GRAV,9815.0,0.0,0.0,-1.0
*CLOAD,FOLLOWER
4,4,10
4,2,-10
5,6,10
5,1,-10
1,5,10
1,3,-10
*INERTIA RELIEF
*OUTPUT,FIELD,variable=preselect
*NODE PRINT,NSET=FIXED
RF
*NODE PRINT,NSET=FIXED
U
*NODE PRINT
CF
***EL PRINT,ELSET=EALL
**IRX
**IRA
**IRF
**IRRI
**IRMASS
*OUTPUT,HISTORY
*ELEMENT OUTPUT,ELSET=EALL
S
E
IRX
IRA
IRF
IRRI
IRMASS
*END STEP