*HEADING
 TIRE AIR COUPLED ACOUSTIC MODEL
 3D FULL TIRE MODEL
 STEP 0: TRANSFER BASE STATE FROM tire_acoustic_rev
 STEP 1: ESTABLISH EQUILIBRIUM
 STEP 2: CHANGE FRICTION
 STEP 3: FREQUENCY EXTRACTION
 STEP 4: STEADY STATE DYNAMICS WITH VERTICAL ROAD
 EXCITATION
 STEP 5: STEADY STATE DYNAMICS WITH HORIZONTAL
 ROAD EXCITATION
 UNITS: KG, M
*RESTART,WRITE,FREQ=100
*SYMMETRIC MODEL GENERATION,REFLECT=LINE,
ELEMENT=100000,NODE=100000
 0.0, 0.0, 0.0,    0.0, 0.0, 1.0
*SYMMETRIC RESULTS TRANSFER
*NSET,NSET=AIRCROWN
211,230,231
*NSET,NSET=AIRFOOT
15211,15230,15231
*NSET,NSET=AIRSIDE
3211,3230,3231
*NSET,NSET=SIDE
26012,26015,26018
********************************************
*STEP,INC=100,NLGEOM=YES
 1:  FOOTPRINT (LOAD CONTROLLED)
*STATIC
 1.0, 1.0
*BOUNDARY,OP=NEW
 TIRE-RIM,  1, 6
 ROAD, 1, 2
 ROAD, 4, 6
*DSLOAD,OP=NEW
 TIRE-LINER, P, 200.E3
*CLOAD, OP=NEW
 ROAD, 3, 3300.
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT
 U,
*NODE OUTPUT,NSET=AIR
 POR,
*ELEMENT OUTPUT,VAR=PRE
*OUTPUT,HISTORY,FREQ=1,VAR=PRE
*END STEP
********************************************
*STEP,INC=100,NLGEOM=YES
 2:  CHANGE FRICTION
*STATIC
 1.0, 1.0
*CHANGE FRICTION,INTERACTION=SRIGID
*FRICTION
 0.1,
*BOUNDARY,OP=NEW
 TIRE-RIM,  1, 6
 ROAD, 2, 2
 ROAD, 4, 6
*END STEP
*******************************
*STEP,NLGEOM
 3: EXTRACT EIGENVALUES
*FREQUENCY,EIGENSOLVER=LANCZOS,PROPERTY=230.
 12, 200.0, 260.0
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT
 U,  
*NODE OUTPUT,NSET=AIR
 POR,
*ELEMENT OUTPUT,VAR=PRE
*OUTPUT,HISTORY,FREQ=1,VAR=PRE
*END STEP
*******************************
*STEP,INC=50,NLGEOM
 4: RESPONSE TO VERTICAL EXCITATION
*STEADY STATE DYNAMICS, DIRECT,INTER=EIGEN
 225.0, 235.0, 3
** 200.0, 260.0, 7 -> to produce manual figures
*CLOAD
 ROAD, 3, 1.
*NODE PRINT,FREQ=0
*EL PRINT,FREQ=0
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT
 U,  
*NODE OUTPUT,NSET=AIR
 POR,
*ELEMENT OUTPUT,VAR=PRE
*OUTPUT,HISTORY,FREQ=1,VAR=PRE
*NODE OUTPUT,NSET=AIRCROWN
POR,
*NODE OUTPUT,NSET=AIRFOOT
POR,
*NODE OUTPUT,NSET=AIRSIDE
POR,
*NODE OUTPUT,NSET=SIDE
U,
*NODE OUTPUT,NSET=ROAD
U,
*END STEP
*******************************
*STEP,INC=50,NLGEOM
 4: RESPONSE TO VERTICAL EXCITATION
*STEADY STATE DYNAMICS, SUBSPACE PROJECTION=CONSTANT,INTER=EIGEN
 225.0, 235.0, 3
** 200.0, 260.0, 7 -> to produce manual figures
*CLOAD
 ROAD, 3, 1.
*NODE PRINT,FREQ=0
*EL PRINT,FREQ=0
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT
 U,  
*NODE OUTPUT,NSET=AIR
 POR,
*ELEMENT OUTPUT,VAR=PRE
*OUTPUT,HISTORY,FREQ=1,VAR=PRE
*NODE OUTPUT,NSET=AIRCROWN
POR,
*NODE OUTPUT,NSET=AIRFOOT
POR,
*NODE OUTPUT,NSET=AIRSIDE
POR,
*NODE OUTPUT,NSET=SIDE
U,
*NODE OUTPUT,NSET=ROAD
U,
*END STEP
*******************************
**STEP,INC=50,NLGEOM
** 5: RESPONSE TO HORIZONTAL EXCITATION
**STEADY STATE DYNAMICS, DIRECT,INTER=EIGEN
** 200.0, 260.0, 7 
***CLOAD
** ROAD, 1, 1.
***END STEP