*HEADING
 JOINTC WITH TFV DEPENDENCIES; MODELING AN AUTOMOBILE SUSPENSION  
*RESTART, WRITE, FREQUENCY=5
**
** Auto suspension sub-assembly modeled in sprung condition.
** Coil spring and tire spring have initial preload as
** given by their respective spring load-deflection curves.
** Weight of frame (by gravity load) rests on the coil spring top
** and the two wishbone bearing (pivot) points.
**
** FRAME
**
** Model frame as lumped mass at auto centerline
*NODE, NSET=FRAME
 1001,       0.00,       0.00,       0.00
*ELEMENT, TYPE=MASS, ELSET=FRAME
 1001, 1001
** Mass used is one quarter of total auto mass
*MASS, ELSET=FRAME
 0.423, 
*BOUNDARY
 1001, 1,2, 0.0
 1001, 4,6, 0.0
**   
** WHEEL
**
** 3001 is lower rim, 3002 is wheel center
*NODE, NSET=WHEEL
 3001,       0.00,    -779.50,    -315.00
 3002,       0.00,    -779.50,       0.00
*ELEMENT, TYPE=B31, ELSET=WHEEL
 3001, 3001, 3002
*BEAM SECTION, SECTION=CIRC, MATERIAL=STEEL, ELSET=WHEEL
 24.0, 
  0.0,      1.0,      0.0
*BOUNDARY
 WHEEL, 6,6, 0.0
**
** TIRE
**
*NODE, NSET=GROUND
 2000,       0.00,       0.00,    -335.00
*NODE, NSET=TIRE
 2001,       0.00,    -779.50,    -335.00
 2091,     167.50,    -779.50,    -290.11
 2092,    -167.50,    -779.50,    -290.11
*ORIENTATION, NAME=XYZ
 1.0, 0.,0.,0.,1.,0.
 3, 0.0
*ELEMENT, TYPE=JOINTC, ELSET=TIRE
 2001, 3001, 2001
*JOINT, ELSET=TIRE, ORIENTATION=XYZ
*SPRING,DEPENDENCIES=1
 1, 
 175.,,,0.
 175.,,,1.
*SPRING,DEPENDENCIES=1
 2, 
 175.,,,0.
 175.,,,1.
*SPRING, NONLINEAR,DEPENDENCIES=1
 3, 
 -10000., -24.21,,0.
      0., -19.21,,0.
   4150.,   0.00,,0.
  10630.,  30.00,,0.
 -10000., -24.21,,1.
      0., -19.21,,1.
   4150.,   0.00,,1.
  10630.,  30.00,,1.
** for visualization only
*ELEMENT, TYPE=T3D3, ELSET=PLOTA
 2091, 2091, 2001, 2092
*SOLID SECTION, MATERIAL=STEEL, ELSET=PLOTA
 0.0005, 
*MPC
 BEAM, 2091, 2001
 BEAM, 2092, 2001
**
** Tie bottom of tire to ground (makes global rotation possible)
*MPC
 BEAM, 2001, 2000
**
** AXLE
**
*NODE, NSET=AXLE
 4001,       0.00,    -705.30,       0.00
 4002,       0.00,    -736.30,       0.00
 4003,       2.56,    -750.00,    -116.00
 4004,     152.00,    -695.61,     -13.81
 4005,     152.00,    -695.61,     -13.81
*ELEMENT, TYPE=B31, ELSET=AXLE1
 4020, 3002, 4001
 4021, 4001, 4002
*BEAM SECTION, SECTION=CIRC, MATERIAL=STEEL, ELSET=AXLE1
 13.0, 
*ELEMENT, TYPE=B31, ELSET=AXLE2
 4022, 4002, 4003
*BEAM SECTION, SECTION=CIRC, MATERIAL=STEEL, ELSET=AXLE2
 15.0, 
  1.0,      0.0,      0.0
*ELEMENT, TYPE=B31, ELSET=AXLE3
 4023, 4002, 4004
*BEAM GENERAL SECTION, SECTION=GENERAL, ELSET=AXLE3, DENSITY=7.8E-9
 391.0,  17000.0,      0.0,     9400.0,    22600.0  
   0.0,      0.0,     -1.0
   2.1E+05,  8.077E+04
*BOUNDARY
 AXLE, 5,5, 0.0
**
** WISHBONE  
**
*NODE, NSET=WBONE
 5001,     -73.00,    -351.00,     -90.00
 5002,     325.00,    -351.00,     -90.00
 5011,     -17.00,    -540.00,    -128.00
 5012,     -73.00,    -398.60,    -110.70
 5013,     325.00,    -398.60,    -110.70
 5080,       2.56,    -750.00,    -116.00
*ELEMENT, TYPE=B31, ELSET=WB1
 5001, 5080, 5011
 5002, 5011, 5012
 5003, 5011, 5013
*BEAM GENERAL SECTION, SECTION=GENERAL, ELSET=WB1, DENSITY=7.8E-9
 542.0, 449591.0,       0.0,   196281.0,   639175.0
   0.0,      0.0,     -1.0
   2.1E+05,  8.077E+04
*ELEMENT, TYPE=B31, ELSET=WB2
 5033, 5012, 5001
 5034, 5013, 5002
*BEAM GENERAL SECTION, SECTION=GENERAL, ELSET=WB2, DENSITY=7.8E-9
 342.0,  81077.0,      0.0,   81077.0,  162155.0  
   0.0,      0.0,     -1.0
   2.1E+05,  8.077E+04
*ELEMENT, TYPE=B31, ELSET=WB3
 5035, 5012, 5013
*BEAM SECTION, SECTION=CIRC, MATERIAL=STEEL, ELSET=WB3
 13.0, 
**
*ELSET, ELSET=WBONE
 WB1, WB2, WB3
** Pin axle and wishbone together
*MPC
 PIN, 4003, 5080
**
** COIL SPRING
**
*NODE, NSET=CTOP
 6001,      -8.50,    -514.00,     214.00
*ELEMENT, TYPE=SPRINGA, ELSET=CSPRG
 6001, 5011, 6001
*SPRING, ELSET=CSPRG, NONLINEAR, DEPENDEN=2

 -13600.0, -100.00,0.,0.,0.
  -8900.0,    0.00,0.,0.,0.
  -4600.0,  100.00,0.,0.,0.
 -13600.0, -100.00,1.,0.,0.
  -8900.0,    0.00,1.,0.,0.
  -4600.0,  100.00,1.,0.,0.
 -13600.0, -100.00,0.,1.,0.
  -8900.0,    0.00,0.,1.,0.
  -4600.0,  100.00,0.,1.,0.
 -13600.0, -100.00,1.,1.,0.
  -8900.0,    0.00,1.,1.,0.
  -4600.0,  100.00,1.,1.,0.
 -13600.0, -100.00,0.,0.,1.
  -8900.0,    0.00,0.,0.,1.
  -4600.0,  100.00,0.,0.,1.
 -13600.0, -100.00,1.,0.,1.
  -8900.0,    0.00,1.,0.,1.
  -4600.0,  100.00,1.,0.,1.
 -13600.0, -100.00,0.,1.,1.
  -8900.0,    0.00,0.,1.,1.
  -4600.0,  100.00,0.,1.,1.
 -13600.0, -100.00,1.,1.,1.
  -8900.0,    0.00,1.,1.,1.
  -4600.0,  100.00,1.,1.,1.
** Tie top of spring to frame
*MPC
 BEAM, 6001, 1001
** for visualization only
*ELEMENT, TYPE=T3D2, ELSET=PLOTB 
 6091, 5011, 6001
*SOLID SECTION, MATERIAL=STEEL, ELSET=PLOTB 
 0.0005, 
**
** BEARING POINT A  
**
*NODE, NSET=A-PVT
 7001,     -73.00,    -351.00,     -90.00
*ELEMENT, TYPE=JOINTC, ELSET=APVT
 7001, 7001, 5001
** Tie fixed side of jointc to frame
*MPC 
 BEAM, 7001, 1001
*JOINT, ELSET=APVT, ORIENTATION=XYZ
*SPRING, NONLINEAR,DEPENDENCIES=1
 1, 
 -16000.0, -4.8,,0.
 -10000.0, -3.8,,0.
  -6000.0, -2.6,,0.
  -3000.0, -1.1,,0.
      0.0,  0.0,,0.
   3000.0,  1.1,,0.
   6000.0,  2.6,,0.
  10000.0,  3.8,,0.
  16000.0,  4.8,,0.
 -16000.0, -4.8,,1.
 -10000.0, -3.8,,1.
  -6000.0, -2.6,,1.
  -3000.0, -1.1,,1.
      0.0,  0.0,,1.
   3000.0,  1.1,,1.
   6000.0,  2.6,,1.
  10000.0,  3.8,,1.
  16000.0,  4.8,,1.
*SPRING, NONLINEAR,DEPENDENCIES=1
 2, 
 -10000.0, -1.80,,0.
  -6000.0, -1.25,,0.
  -2000.0, -0.53,,0.
      0.0,  0.00,,0.
   2000.0,  0.53,,0.
   6000.0,  1.25,,0.
  10000.0,  1.80,,0.
 -10000.0, -1.80,,1.
  -6000.0, -1.25,,1.
  -2000.0, -0.53,,1.
      0.0,  0.00,,1.
   2000.0,  0.53,,1.
   6000.0,  1.25,,1.
  10000.0,  1.80,,1.
*SPRING, NONLINEAR,DEPENDENCIES=1
 3, 
 -10000.0, -2.40,,0.
  -5000.0, -1.45,,0.
  -2000.0, -0.68,,0.
      0.0,  0.00,,0.
   2000.0,  0.68,,0.
   5000.0,  1.45,,0.
  10000.0,  2.40,,0.
 -10000.0, -2.40,,1.
  -5000.0, -1.45,,1.
  -2000.0, -0.68,,1.
      0.0,  0.00,,1.
   2000.0,  0.68,,1.
   5000.0,  1.45,,1.
  10000.0,  2.40,,1.
*SPRING,DEPENDENCIES=1
 4, 
 1.142E5,,,0.
 1.142E5,,,1.
*SPRING,DEPENDENCIES=1
 5, 
 2.8568E6,,,0.
 2.8568E6,,,1.
*SPRING,DEPENDENCIES=1
 6, 
 2.8568E6,,,0.
 2.8568E6,,,1.
**
** BEARING POINT B
**
*NODE, NSET=B-PVT
 7002,     325.00,    -351.00,     -90.00
*ELEMENT, TYPE=JOINTC, ELSET=BPVT
 7002, 7002, 5002
** Tie fixed side of jointc to frame
*MPC
 BEAM, 7002, 1001
*JOINT, ELSET=BPVT, ORIENTATION=XYZ
*SPRING, NONLINEAR,DEPENDENCIES=1
 1, 
 -1600.0, -6.00,,0.
 -1000.0, -2.50,,0.
  -700.0, -1.50,,0.
  -500.0, -1.00,,0.
  -300.0, -0.50,,0.
     0.0,  0.00,,0.
   300.0,  0.50,,0.
   500.0,  1.00,,0.
   700.0,  1.50,,0.
  1000.0,  2.50,,0.
  1600.0,  6.00,,0.
 -1600.0, -6.00,,1.
 -1000.0, -2.50,,1.
  -700.0, -1.50,,1.
  -500.0, -1.00,,1.
  -300.0, -0.50,,1.
     0.0,  0.00,,1.
   300.0,  0.50,,1.
   500.0,  1.00,,1.
   700.0,  1.50,,1.
  1000.0,  2.50,,1.
  1600.0,  6.00,,1.
*SPRING, NONLINEAR,DEPENDENCIES=1
 2, 
 -6800.0,  -4.00,,0.
 -4800.0,  -3.75,,0.
 -3600.0,  -3.50,,0.
 -2800.0,  -3.25,,0.
 -1700.0,  -2.75,,0.
 -1000.0,  -2.25,,0.
  -700.0,  -1.75,,0.
     0.0,   0.00,,0.
   700.0,   1.75,,0.
  1000.0,   2.25,,0.
  1700.0,   2.75,,0.
  2800.0,   3.25,,0.
  3600.0,   3.50,,0.
  4800.0,   3.75,,0.
  6800.0,   4.00,,0.
 -6800.0,  -4.00,,1.
 -4800.0,  -3.75,,1.
 -3600.0,  -3.50,,1.
 -2800.0,  -3.25,,1.
 -1700.0,  -2.75,,1.
 -1000.0,  -2.25,,1.
  -700.0,  -1.75,,1.
     0.0,   0.00,,1.
   700.0,   1.75,,1.
  1000.0,   2.25,,1.
  1700.0,   2.75,,1.
  2800.0,   3.25,,1.
  3600.0,   3.50,,1.
  4800.0,   3.75,,1.
  6800.0,   4.00,,1.
*SPRING,NONLINEAR,DEPENDENCIES=1
 3, 
 -6400.0, -2.00,,0.
 -5200.0, -1.75,,0.
 -4160.0, -1.50,,0.
 -2400.0, -1.00,,0.
  1000.0,  0.00,,0.
  4400.0,  1.00,,0.
  6160.0,  1.50,,0.
  7200.0,  1.75,,0.
  8400.0,  2.00,,0.
 -6400.0, -2.00,,1.
 -5200.0, -1.75,,1.
 -4160.0, -1.50,,1.
 -2400.0, -1.00,,1.
  1000.0,  0.00,,1.
  4400.0,  1.00,,1.
  6160.0,  1.50,,1.
  7200.0,  1.75,,1.
  8400.0,  2.00,,1.
*SPRING,DEPENDENCIES=1
 4, 
 4.3082E4,,,0.
 4.3082E4,,,1.
*SPRING,DEPENDENCIES=1
 5, 
 1.1467E6,,,0.
 1.1467E6,,,1.
*SPRING,DEPENDENCIES=1
 6, 
 3.1311E5,,,0.
 3.1311E5,,,1.
**
*MATERIAL, NAME=STEEL
*ELASTIC
 2.1E+05,0.3
*DENSITY
 7.8E-9, 
*ELSET, ELSET=ALL
 1001, 
**
** Idealized triangular speed bump 400mm long and 100mm high
** Time data based on 5km/hr auto speed
*AMPLITUDE, NAME=BUMP, TIME=STEP TIME
 0.0, 0.0,  0.05, 0.0,  0.20, 100.0,  0.35, 0.0
 0.40, 0.0
*AMPLITUDE, NAME=STEPFCN
0.0,1.0, 1.0,1.0
**
**
*NSET,NSET=NALL
FRAME,WHEEL,GROUND,TIRE,AXLE,WBONE,CTOP,A-PVT,B-PVT,
*STEP, NLGEOM, INC=100
 Auto Weight/Tire Compression/Coil spring Equilibrium
*STATIC
 0.2, 0.2
*BOUNDARY
 2000, 1,6, 0.0
*DLOAD
 ALL, GRAV, 9815.0, 0.0, 0.0, -1.0
*EL PRINT,FREQUENCY=0
*NODE PRINT,FREQUENCY=0
*NSET, NSET=PLOT
 2000, 3002, 1001
*NODE FILE, NSET=PLOT
 U, 
*FIELD,VARIABLE=1,AMP=STEPFCN
NALL,0.5
*FIELD,VARIABLE=2,AMP=STEPFCN
NALL,0.5
*END STEP
**
**
*STEP, NLGEOM, INC=200
 Apply "bump" to tire bottom with boundary condition
*DYNAMIC, HAFTOL=5000
 5.0E-4, 0.40, , 0.01
*BOUNDARY, AMPLITUDE=BUMP
 2000, 3,3, 1.0 
*EL FILE, ELSET=WBONE
 XC, 
 UC, 
 VC, 
 HC, 
 HO, 
 RI, 
 MASS, VOL
*FIELD,VARIABLE=1,AMP=STEPFCN
NALL,0.5
*FIELD,VARIABLE=2,AMP=STEPFCN
NALL,0.5
*OUTPUT,FIELD,VAR=PRESELECT,FREQ=1
*OUTPUT,HISTORY,VAR=PRESELECT,FREQ=1
*END STEP