*HEADING
 DESIGN SENSTIIVITY ANALYSIS FOR AXISYMMETRIC TIRE MODEL
 STEP 1: INFLATE TIRE TO 200 KPa
 STEP 2: FREQUENCY ANALYSIS ABOUT INFLATED STATE
 UNITS: KG, M 
*PARAMETER
tNominal=0.0
rubDen=1100.
beltE=172.2E+9
carcE=9.87E+9
*DESIGN PARAMETER
tNominal,rubDen,beltE,carcE
*DSA CONTROLS, SIZING FREQUENCY=1
*NODE,NSET=NTIRE,INP=dsatire_axi_half_node.inp
*INCLUDE,INP=dsatire_axi_half_psv.inp
*ELEMENT,TYPE=CAX4H,ELSET=TREAD
    1,   50,   55,   54,   49
    2,   45,   50,   49,   44
    3,   40,   45,   44,   39
    5,   35,   40,   39,   34
    7,   31,   35,   34,   30
*ELEMENT,TYPE=CAX3H,ELSET=TREAD
    4,   27,   31,   30
*ELEMENT,TYPE=CAX4H,ELSET=SIDE
   15,   27,   30,   28,   25
   16,   24,   27,   25,   22
   17,   21,   24,   22,   19
   18,   18,   21,   19,   16
   19,   15,   18,   16,   13
   20,   12,   15,   13,   10
   21,   30,   34,   32,   28
   29,    9,   12,   10,    7
   30,    6,    9,    7,    4
   31,    3,    6,    4,    1
*ELEMENT,TYPE=CAX4H,ELSET=BELT
   35,   49,   54,   51,   46
   36,   44,   49,   46,   41
   37,   39,   44,   41,   36
   38,   34,   39,   36,   32
*REBAR,ELEMENT=CONTINUUM,MATERIAL=BELT,
GEOMETRY=ISO,NAME=BELT1
 BELT, 0.2118683E-6, 1.16E-3,  70.0, 0.50,  3
*REBAR,ELEMENT=CONTINUUM,MATERIAL=BELT,
GEOMETRY=ISO,NAME=BELT2
 BELT, 0.2118683E-6, 1.16E-3, 110.0, 0.83,  3
*REBAR,ELEMENT=CONTINUUM,MATERIAL=CARCASS,
GEOMETRY=ISO,NAME=CARCASS
 BELT, 0.4208352E-6, 1.00E-3, 0.0, 0.0,  3
 SIDE, 0.4208352E-6, 1.00E-3, 0.0, 0.0,  3
*SOLID SECTION,ELSET=TREAD,MATERIAL=RUBBER
*SOLID SECTION,ELSET=SIDE,MATERIAL=RUBBER
*SOLID SECTION,ELSET=BELT,MATERIAL=RUBBER
*MATERIAL,NAME=RUBBER
*HYPERELASTIC,N=1,MODULI=LONG TERM
 1.0e6,
*VISCOELASTIC,TIME=PRONY
 0.3, 0.0, 0.1
*DENSITY
 <rubDen>
*MATERIAL,NAME=BELT
*ELASTIC,TYPE=ISO
 <beltE>, 0.3
*DENSITY
 5900.,
*MATERIAL,NAME=CARCASS
*ELASTIC,TYPE=ISO
 <carcE>,  0.3
*DENSITY
 1500.,
*NSET,NSET=RIM
 6, 3, 1
*NSET,NSET=SYM
 51, 54, 55
*ELSET,ELSET=SOLID,GENERATE
 1, 76,1
*ELSET,ELSET=EALL
 TREAD,SIDE,BELT
*SURFACE,NAME=INSIDE
 BELT, S3
 SIDE, S3
***********************************
*STEP,DSA=YES,NLGEOM=YES
 1:  INFLATION
*STATIC, LONG TERM
 0.25, 1.0
*CONTROLS,PARAMETERS=FIELD
1.e-5,
*BOUNDARY
 RIM, 1, 2
 RIM, 5, 
 SYM, 2,
 SYM, 5,
*DSLOAD
 INSIDE, P, 200.E3
*EL PRINT,FREQUENCY=50,REBAR,ELSET=SOLID
 S,
*NODE PRINT,FREQUENCY=100,TOTAL=YES
 U, RF
*OUTPUT,FIELD,FREQ=1
*ELEMENT OUTPUT
 S
*NODE OUTPUT
 U,RF
*DESIGN RESPONSE,FREQ=1
*ELEMENT RESPONSE
 S
*NODE RESPONSE
 U,RF
*END STEP
**************************************************
*STEP,DSA=YES
 4:  FREQUENCY
*FREQUENCY,EIGENSOL=LANCZOS
 5,
*OUTPUT,HISTORY
*ELEMENT OUTPUT,ELSET=EALL
 MASS
*DESIGN RESPONSE
*ELEMENT RESPONSE,ELSET=EALL
 MASS
*END STEP