*HEADING
 CROSS SECTION FOR CIRCULAR DISK
 REFERENCE FILE FOR:     pstrebar
*PREPRINT,MOD=Y
*RESTART,WRITE
*NODE,NSET=AXLE
 99, 0.0, 0.5
*NODE,NSET=NTIRE
 101,   1.0, 0.0
 106,   2.0, 0.0
 201,   1.0, 1.0
 206,   2.0, 1.0
*NGEN,NSET=NSIDE1
  101, 106
*NGEN,NSET=NSIDE2
  201, 206
*NSET,NSET=NTIRE
 NSIDE1, NSIDE2
*NSET,NSET=NRIM
 101, 201
*NSET,NSET=NTREAD
 106, 206
*ELEMENT,TYPE=CAX4RH,ELSET=ETIRE
 101,  101,102,202,201
*ELEMENT,TYPE=SFMAX1,ELSET=ESHELL
 115,  105, 205
 123,  103, 104
 124,  104, 105
 133,  203, 204
 134,  204, 205
*ELGEN,ELSET=ETIRE
 101, 5
*ELSET,ELSET=ERIM
 101, 
*ELSET,ELSET=ETREAD
 105, 
*SURFACE,NAME=RIM,TYPE=ELEM
 101, S4
*ELSET,ELSET=ES
 124
*COUPLING,CONSTRAINT=RIM,REFNODE=AXLE,SURFACE=RIM
*KINEMATIC
 1,6
*SECTION CONTROL,HOUR=ENH,NAM=XX
*SOLID SECTION,MATERIAL=HYPER1,ELSET=ETIRE,CON=XX
*SURFACE SECTION,ELSET=ESHELL
 0.2,
*REBAR LAYER
 RB1,0.1,1.0,0,HYPER2,90,1 
 RB2,0.1,1.0,0,MARLOW,90,1 
*MATERIAL,NAME=HYPER1
*HYPERELASTIC
 80.0, 20.0
*DENSITY
 0.036, 
*MATERIAL,NAME=HYPER2
*HYPERELASTIC
 80.0, 20.0, 0.001
*VISCOELASTIC,TIME=PRONY
 0.4, 0.2, 0.1
*DENSITY
 0.036, 
*MATERIAL,NAME=MARLOW
*HYPERELASTIC,POIS=0.5,MARLOW
*UNIAXIAL TEST DATA
 -1813.50,  -.600          
 -1277.02,  -.545          
 -925.863,  -.490          
 -683.478,  -.435          
 -508.787,  -.380          
 -378.258,  -.325          
 -277.679,  -.270          
 -198.080,  -.215          
 -133.588,  -.160          
 -80.2357,  -.105          
 -75.8354,  -.100         
 -67.2343,  -.090         
 -58.8886,  -.080         
 -50.7862,  -.070         
 -42.9157,  -.060         
 -35.2664,  -.050         
 -27.8280,  -.040         
 -20.5909,  -.030         
 -13.5462,  -.020         
 -6.68520,  -.010         
   0     ,    0  
   6.5172,   .010         
  12.8737,   .020         
  19.0762,   .030         
  25.1312,   .040         
  31.0450,   .050         
  36.8234,   .060         
  42.4719,   .070         
  47.9959,   .080         
  53.4002,   .090         
  58.6897,   .100         
   66.416,   .115          
   92.853,   .170          
  116.737,   .225          
  138.532,   .280          
  158.606,   .335          
  177.246,   .390          
  194.686,   .445          
  211.111,   .500          
*VISCOELASTIC,TIME=PRONY
 0.4, 0.0, 0.1
*****
*STEP,INC=100,NLGEOM
 DO NOTHING
*STATIC
*BOUNDARY
 NSIDE1, 2
*NODE PRINT,FREQ=0
*EL PRINT,FREQ=0
*NODE FILE,NSET=NRIM
 U, RF
*ELFILE,ELSET=ETREAD
 S,
*END STEP