*HEADING
 DESIGN SENSITIVITY ANALYSIS FOR TIRE MODEL
 STEP 1: INFLATE TIRE TO 200 KPa
 STEP 2: FOOTPRINT ANALYSIS (DISPLACEMENT CONTROL)
 STEP 3: FOOTPRINT ANALYSIS (LOAD CONTROL)
 STEP 4: FREQUENCY ANALYSIS AROUND FOOTPRINT STATE
 UNITS: KG, M
*PARAMETER
tNominal = 0.
denRub   = 1100.
Ebelt    = 172.2E+09
Ecarc    = 9.87E+9
*DESIGN PARAMETER
tNominal,denRub,Ebelt,Ecarc
*INCLUDE,INP=dsatire_model.inp
*INCLUDE,INP=dsatire_psv.inp
*NODE,NSET=ROAD
 9999,   0.0,   0.0, -0.02
*ELSET,ELSET=FOOT,GEN
 1001, 4801, 200
 1002, 4802, 200
 1003, 4803, 200
 1004, 4804, 200
 1005, 4805, 200 
 1007, 4807, 200
11001,14801, 200
11002,14802, 200
11003,14803, 200
11004,14804, 200
11005,14805, 200
11006,14806, 200
11007,14807, 200
*SURFACE,TYPE=CYLINDER,NAME=SROAD
 0., 0.,-0.31657,   1., 0.,-0.31657
 0., 1.,-0.31657
 START, -0.3, 0.
 LINE,   0.3, 0.
*RIGID BODY,REF NODE=9999,ANALYTICAL SURFACE=SROAD
*SURFACE,NAME=STREAD
 FOOT, S3
*CONTACT PAIR,INTERACTION=SRIGID
 STREAD, SROAD
*SURFACE INTERACTION,NAME=SRIGID
*FRICTION
 0.0
*ELSET,ELSET=SECT,GENERATE
 2800, 3200, 1
*ELSET,ELSET=EALL
 TREAD,SIDE,BELT,SOLID
*NSET,NSET=SECT,GENERATE
 2800, 3400, 1
*NSET,NSET=FOOT,ELSET=FOOT
*NSET,NSET=NOUTP,GENERATE
 1055, 5055, 200
**************************************************
*STEP,INC=100,NLGEOM=YES,DSA=YES,UNSYMM=YES
 1:  INFLATION
*STATIC, LONG TERM
 0.25, 1.0
*CONTROLS,PARAMETER=FIELD
 1.e-5
*DSA CONTROLS
**averages over increments for similar axisymmetric problem
tNominal, CD, 7.6E-05
denRub,   CD, 1.7E-04
Ebelt,    CD, 7.2E+06
Ecarc,    CD, 4.5E+05
*BOUNDARY
 RIM, 1, 3
*DSLOAD
 INSIDE, P, 200.E3
*OUTPUT,FIELD,FREQ=0
*ELEMENT OUTPUT,ELSET=FOOT
 S
*NODE OUTPUT,NSET=FOOT
 U,RF
*CONTACT OUTPUT
 CSTRESS
*END STEP
**************************************************
*STEP,INC=100,NLGEOM=YES,UNSYMM=YES
 2:  FOOTPRINT (Displacement controlled)
*STATIC, LONG TERM
 0.2, 1.0
*BOUNDARY,OP=NEW
 RIM,  1, 3
 ROAD, 1, 2
 ROAD, 4, 6
 ROAD, 3, , 0.02
*OUTPUT,FIELD,FREQ=1
*ELEMENT OUTPUT,ELSET=FOOT
 S
*NODE OUTPUT,NSET=FOOT
 U,RF
*CONTACT OUTPUT
 CSTRESS
*DESIGN RESPONSE,FREQ=1
*ELEMENT RESPONSE,ELSET=FOOT
 S
*NODE RESPONSE,NSET=FOOT
 U,RF
*CONTACT RESPONSE
 CSTRESS
*END STEP
**************************************************
*STEP,INC=100,NLGEOM=YES,UNSYMM=NO
 3:  FOOTPRINT (Load controlled)
*STATIC, LONG TERM
 1.0, 1.0
*BOUNDARY,OP=NEW
 RIM,  1, 3
 ROAD, 1, 2
 ROAD, 4, 6
*CLOAD,OP=NEW
 ROAD, 3, 3300.
*END STEP
**************************************************
*STEP,DSA=YES
 4:  FREQUENCY
*FREQUENCY,EIGENSOL=LANCZOS
 5,
*DSA CONTROLS
**based on similar axisymmetric problem
tNominal, CD, 1.5E-06
denRub,   CD, 1.7E-04
Ebelt,    CD, 2.6E+05
Ecarc,    CD, 1.5E+05
*OUTPUT,HISTORY
*DESIGN RESPONSE
*END STEP