

  head 2 "Test of different wheel/rail couplings  wr_coupl_polach_2_varmu"
  head 3 "axl_111 pure long.creep; axl_112 pure lat.creep; axl_113 pure spin"
  head 4 "Cylindric wheel profile  mu= $mu_"

##
## Analyse specific input data
##
  idebug= 2                                       # Definition of different types of debug printing.
  idebug_file= out/wr_coupl_polach_2_varmu.out        # File to where the debug information will be written to

  fresp_param= Fourier_CG1  1.e-4  fstart .2   fstop 20.  fstep -1.05
  modal_param= Schur_fact1  1.e-4
  quasi_param= Damped_Tens  1.e-4 0.1   20    2.
  tsim_param = heun_c            tstart 0.  tstop  3.00  tstep 0.0001 tout .005

##
##  Track irregularities and designed track geometry
##
  func	  intpl_r   lat_trac	  -200    0.   0    0.
  func	  intpl_r   spv_trac	  -200 1435.   0 1435.
  func	  intpl_r   vert_trac	  -200    0.   0    0.
  func	  intpl_r   fi_trac	  -200    0.   0    0.

  func	  intpl_r   krok_R  -200 0   0 0
  func	  intpl_r   krok_h  -200 0   0 0
  func	  intpl_r   krok_z  -200 0   0 0

##
##  Wheel-rail geometry
##
# insert file $genkpf/kona_0.000_uf.kpfr
# in_substruct kpf_kona_0.000 [ 1 ]    # Same wheel-rail geometry in all wheels vehicle #1
#
  insert file $genkpf/kpf_rkona_cp1.kpfr
  in_substruct kpf_rkona_cp1 [ 1 0. ]    # Same wheel-rail geometry in all wheels vehicle #1

##
##  Define constants for the wheel-rail-coupling substructure wr_coupl
##
  func	const_block
   Vo = 10.   # Speed of the vehicles
   mu_= 0.5   # Friction coefficient, all wheels in the train
   ro_= 0.45  # Wheel radius, all wheels in the train
   bo_= 0.75  # Semi distance between the wheel profile origins, all wheels in the train
   ma_= 1887  # Mass of wheelset

   YMtrac= 1e-3         # From mm to m;   Scale factor 1.00; Lateral irregularities
   ZMtrac= 1e-3         # From mm to m;   Scale factor 1.00; Vertical irregularities
   GMtrac= 1e-3         # From mm to m;   Scale factor 1.00; Gauge irregularities
   CMtrac= 1e-3/(2*bo_) # From mm to rad; Scale factor 1.00; Cant irregularities

   gauge_average= 1435.    # Average gauge of the track
   gauge_dev_   = 0.000    # Deviation in gauge
  end_block

##
##  Create the wheelsets.
##  ---------------------
  lsys	 e_abs_bendrf  esys  Vo  0.  krok_R krok_h krok_z
                                     4 4    4 4    4 4
  substruct create_axles [
   lsys  l_local      lsa_$1  esys     0. 0. 0.
   mass  m_rigid_6    axl_$1  lsa_$1   0. 0. -ro_   ma_ ma_ ma_  720  100  720
   mass  m_rigid_6    trc_$1  lsa_$1   0. 0. 0.    1.e3 1.e3 1.e3 1.e3 1.e3 1.e3
   mass  fixpoint_6   grd_$1  lsa_$1   0. 0. 0.
   body  w_set_lsys   axl_$1  lsa_$1   ro_ bo_
   body  figure_lsys  trc_$1  lsa_$1   1 5
                                       .5 1 0  .5 -1 0   -.5 -1 0  -.5 1 0  .5 1 0
#
   constr fix_rigid_1  axl_$1  x   0.
   constr fix_rigid_1  axl_$1  y   0.
   constr fix_rigid_1  axl_$1  f   0.
   constr fix_free_1   axl_$1.vk   0.
   constr fix_rigid_1  axl_$1  p   0.
   constr fix_rigid_1  trc_$1  y   0.
#
   force rel_lsys1 vload_$1  axl_$1  0. 0. -ro_  0. 0. 200e3-ma_*9.81  0. 0. 0.
  ]
#
  in_substruct create_axles [ 111 ]
  in_substruct create_axles [ 112 ]
  in_substruct create_axles [ 113 ]
#
  mass	fixpoint_6   car_1    lsa_111	0. 0. 0.

##
##  Define pure creepage on the three different axles
##  -------------------------------------------------
  no_warning func operp axl_111.vk= -Vo/ro_ - .25 * time

  no_warning func operp axl_112.p =  .1 * time

  no_warning func operp axl_113.vp= -16 * time
  no_warning func operp axl_113.vx= -12 * time

##
##        Calculation of creep and creepforces between wheels and rails
##        ================================================================
   func const knwr.F0_= 100e3 # all wheels
#
   func  intpl_r knwr_      # Stiffness Tread-Rail all wheels
   -1.57     600e6
    0.       600e6
    1.57     600e6
#
   func  intpl_r cnwr_    # Damping Tread-Rail all wheels
    -1.57    60e3
     0.     250e3
     1.57    60e3

  func const  mu.O= 0.305         # Coefficient of friction at zero slip
  func const  mu.O= 0.5           # Coefficient of friction at zero slip
  func const  mu.A= 0.50          # Ratio mu.Inf/mu.O
  func const  mu.a= 1-mu.A
  func const  mu.B= -.16          # Coefficient of exponential friction decrease (s/m)
  func const  mu.B= -.16*10       # Coefficient of exponential friction decrease (s/m)

  substruct wr_coupl [
#
   func const  mu_$1l1= mu.O   func const  mu_$1r1= mu.O    # Current coefficient of friction
   func const  mu_$1l2= mu.O   func const  mu_$1r2= mu.O    # will be set later when current 
   func const  mu_$1l3= mu.O   func const  mu_$1r3= mu.O    # slip has been calculated
#
   mass m_rigid_6f  ral_$1r lsa_$1  0.  bo_  0.    0. 0. 0.  0. 0. 0.
   mass m_rigid_6f  ral_$1l lsa_$1  0. -bo_  0.    0. 0. 0.  0. 0. 0.
   constr fix_rigid_1  ral_$1r x 0.  constr fix_rigid_1  ral_$1r y 0.  constr fix_rigid_1  ral_$1r z 0.
   constr fix_rigid_1  ral_$1r f 0.  constr fix_rigid_1  ral_$1r k 0.  constr fix_rigid_1  ral_$1r p 0.
   constr fix_rigid_1  ral_$1l x 0.  constr fix_rigid_1  ral_$1l y 0.  constr fix_rigid_1  ral_$1l z 0.
   constr fix_rigid_1  ral_$1l f 0.  constr fix_rigid_1  ral_$1l k 0.  constr fix_rigid_1  ral_$1l p 0.
#
   constr fix_free_1  axl_$1.k  0.
   constr fix_rigid_1 trc_$1  x 0.  constr fix_rigid_1 trc_$1  z 0.
   constr fix_rigid_1 trc_$1  f 0.  constr fix_rigid_1 trc_$1  k 0.
   constr fix_rigid_1 trc_$1  p 0.
#
   func const axl_$1.vk=  -Vo/ro_ 
#
   func  operp whe_$1r.x= axl_$1.x - bo_ * axl_$1.p
   func  operp whe_$1l.x= axl_$1.x + bo_ * axl_$1.p
#
   func  operp whe_$1r.vx= axl_$1.vx - bo_ * axl_$1.vp
   func  operp whe_$1l.vx= axl_$1.vx + bo_ * axl_$1.vp
#
   coupl creep_polach_2 cp1_$1r
                axl_$1  0.0        bo_  0.0
                ral_$1r whe_$1r.x  bo_  0.0  esys m
                  0.    # tral$1r.y        # Lateral position track irregularity
                  0.    # tral$1r.z        # Vertical position track irregularity
                  0.    # tral$1r.vy       # Lateral velocity track irregularity
                  0.    # tral$1r.vz       # Vertical velocity track irregularity
                  1.    # mulfact_nux_$1r  # Longitudinal creep relaxation due to contaminated surfaces
                  1.    # mulfact_nuy_$1r  # Lateral      creep relaxation due to contaminated surfaces
                  1.    # mulfact_spin_$1r # Spin         creep relaxation due to contaminated surfaces
                  cp1_zfn_$1               # Wheel lift geometric function
                  cp1_drfn_$1              # Wheel radius geometric function
                  cp1_gamfn_$1             # Contact angle geometric function
                  cp1_rofn_$1              # Lateral curvature geometric function
                  cp1_poswfn_$1            # Lateral position of contact point on wheel 
                  cp1_posrfn_$1            # Lateral position of contact point on rail
                  knwr.F0_$1r              # Wheel/rail prestess force in normal direction of contact surface
                  knwr_$1r                 # Wheel/rail stiffness in normal direction of contact surface
                  cnwr_$1r                 # Wheel/rail damping   in normal direction of contact surface
                  2.05e11                  # The combined modulus of elasticity for body1 and body2
                  .25                      # The combined Poisson's ratio for body1 and body2
                  .95                      # kA
                  0.90                     # kS
                  mu_$1r1                  # The coefficient of friction in the contact surface
                  ro_$1r                   # Wheel rolling radius, will vary for OOR-wheels
                  1e99                     # Vertical radius of rail                                     
#
   coupl creep_polach_2 cp1_$1l
                axl_$1  0.0       -bo_  0.0
                ral_$1l whe_$1l.x -bo_  0.0  esys m
                  0.    # tral$1l.y        # Lateral position track irregularity
                  0.    # tral$1l.z        # Vertical position track irregularity
                  0.    # tral$1l.vy       # Lateral velocity track irregularity
                  0.    # tral$1l.vz       # Vertical velocity track irregularity
                  1.    # mulfact_nux_$1l  # Longitudinal creep relaxation due to contaminated surfaces
                  1.    # mulfact_nuy_$1l  # Lateral      creep relaxation due to contaminated surfaces
                  1.    # mulfact_spin_$1l # Spin         creep relaxation due to contaminated surfaces
                  cp1_zfn_$1               # Wheel lift geometric function
                  cp1_drfn_$1              # Wheel radius geometric function
                  cp1_gamfn_$1             # Contact angle geometric function
                  cp1_rofn_$1              # Lateral curvature geometric function
                  cp1_poswfn_$1            # Lateral position of contact point on wheel 
                  cp1_posrfn_$1            # Lateral position of contact point on rail
                  knwr.F0_$1l              # Wheel/rail prestess force in normal direction of contact surface
                  knwr_$1l                 # Wheel/rail stiffness in normal direction of contact surface
                  cnwr_$1l                 # Wheel/rail damping   in normal direction of contact surface
                  2.05e11                  # The combined modulus of elasticity for body1 and body2
                  .25                      # The combined Poisson's ratio for body1 and body2
                  .95                      # kA
                  0.90                     # kS
                  mu_$1l1                  # The coefficient of friction in the contact surface
                  ro_$1l                   # Wheel rolling radius, will vary for OOR-wheels
                  1e99                     # Vertical radius of rail                                     
#
#
   coupl creep_polach_2 cp2_$1r
                axl_$1  0.0        bo_  0.0
                ral_$1r whe_$1r.x  bo_  0.0  esys m
                  0.    # tral$1r.y        # Lateral position track irregularity
                  0.    # tral$1r.z        # Vertical position track irregularity
                  0.    # tral$1r.vy       # Lateral velocity track irregularity
                  0.    # tral$1r.vz       # Vertical velocity track irregularity
                  1.    # mulfact_nux_$1r  # Longitudinal creep relaxation due to contaminated surfaces
                  1.    # mulfact_nuy_$1r  # Lateral      creep relaxation due to contaminated surfaces
                  1.    # mulfact_spin_$1r # Spin         creep relaxation due to contaminated surfaces
                  cp2_zfn_$1               # Wheel lift geometric function
                  cp2_drfn_$1              # Wheel radius geometric function
                  cp2_gamfn_$1             # Contact angle geometric function
                  cp2_rofn_$1              # Lateral curvature geometric function
                  cp2_poswfn_$1            # Lateral position of contact point on wheel 
                  cp2_posrfn_$1            # Lateral position of contact point on rail
                  0.                       # Wheel/rail prestess force in normal direction of contact surface
                  knwr_$1r                 # Wheel/rail stiffness in normal direction of contact surface
                  cnwr_$1r                 # Wheel/rail damping   in normal direction of contact surface
                  2.05e11                  # The combined modulus of elasticity for body1 and body2
                  .25                      # The combined Poisson's ratio for body1 and body2
                  .95                      # kA
                  0.90                     # kS
                  mu_$1r2                  # The coefficient of friction in the contact surface
                  ro_$1r                   # Wheel rolling radius, will vary for OOR-wheels
                  1e99                     # Vertical radius of rail                                     
#
   coupl creep_polach_2 cp2_$1l
                axl_$1  0.0       -bo_  0.0
                ral_$1l whe_$1l.x -bo_  0.0  esys m
                  0.    # tral$1l.y        # Lateral position track irregularity
                  0.    # tral$1l.z        # Vertical position track irregularity
                  0.    # tral$1l.vy       # Lateral velocity track irregularity
                  0.    # tral$1l.vz       # Vertical velocity track irregularity
                  1.    # mulfact_nux_$1l  # Longitudinal creep relaxation due to contaminated surfaces
                  1.    # mulfact_nuy_$1l  # Lateral      creep relaxation due to contaminated surfaces
                  1.    # mulfact_spin_$1l # Spin         creep relaxation due to contaminated surfaces
                  cp2_zfn_$1               # Wheel lift geometric function
                  cp2_drfn_$1              # Wheel radius geometric function
                  cp2_gamfn_$1             # Contact angle geometric function
                  cp2_rofn_$1              # Lateral curvature geometric function
                  cp2_poswfn_$1            # Lateral position of contact point on wheel 
                  cp2_posrfn_$1            # Lateral position of contact point on rail
                  0.                       # Wheel/rail prestess force in normal direction of contact surface
                  knwr_$1l                 # Wheel/rail stiffness in normal direction of contact surface
                  cnwr_$1l                 # Wheel/rail damping   in normal direction of contact surface
                  2.05e11                  # The combined modulus of elasticity for body1 and body2
                  .25                      # The combined Poisson's ratio for body1 and body2
                  .95                      # kA
                  0.90                     # kS
                  mu_$1l2                  # The coefficient of friction in the contact surface
                  ro_$1l                   # Wheel rolling radius, will vary for OOR-wheels
                  1e99                     # Vertical radius of rail                                     
#
#
   coupl creep_polach_2 cp3_$1r
                axl_$1  0.0        bo_  0.0
                ral_$1r whe_$1r.x  bo_  0.0  esys m
                  0.    # tral$1r.y        # Lateral position track irregularity
                  0.    # tral$1r.z        # Vertical position track irregularity
                  0.    # tral$1r.vy       # Lateral velocity track irregularity
                  0.    # tral$1r.vz       # Vertical velocity track irregularity
                  1.    # mulfact_nux_$1r  # Longitudinal creep relaxation due to contaminated surfaces
                  1.    # mulfact_nuy_$1r  # Lateral      creep relaxation due to contaminated surfaces
                  1.    # mulfact_spin_$1r # Spin         creep relaxation due to contaminated surfaces
                  cp3_zfn_$1               # Wheel lift geometric function
                  cp3_drfn_$1              # Wheel radius geometric function
                  cp3_gamfn_$1             # Contact angle geometric function
                  cp3_rofn_$1              # Lateral curvature geometric function
                  cp3_poswfn_$1            # Lateral position of contact point on wheel 
                  cp3_posrfn_$1            # Lateral position of contact point on rail
                  0.                       # Wheel/rail prestess force in normal direction of contact surface
                  knwr_$1r                 # Wheel/rail stiffness in normal direction of contact surface
                  cnwr_$1r                 # Wheel/rail damping   in normal direction of contact surface
                  2.05e11                  # The combined modulus of elasticity for body1 and body2
                  .25                      # The combined Poisson's ratio for body1 and body2
                  .95                      # kA
                  0.90                     # kS
                  mu_$1r3                  # The coefficient of friction in the contact surface
                  ro_$1r                   # Wheel rolling radius, will vary for OOR-wheels
                  1e99                     # Vertical radius of rail                                     
#
   coupl creep_polach_2 cp3_$1l
                axl_$1  0.0       -bo_  0.0
                ral_$1l whe_$1l.x -bo_  0.0  esys m
                  0.    # tral$1l.y        # Lateral position track irregularity
                  0.    # tral$1l.z        # Vertical position track irregularity
                  0.    # tral$1l.vy       # Lateral velocity track irregularity
                  0.    # tral$1l.vz       # Vertical velocity track irregularity
                  1.    # mulfact_nux_$1l  # Longitudinal creep relaxation due to contaminated surfaces
                  1.    # mulfact_nuy_$1l  # Lateral      creep relaxation due to contaminated surfaces
                  1.    # mulfact_spin_$1l # Spin         creep relaxation due to contaminated surfaces
                  cp3_zfn_$1               # Wheel lift geometric function
                  cp3_drfn_$1              # Wheel radius geometric function
                  cp3_gamfn_$1             # Contact angle geometric function
                  cp3_rofn_$1              # Lateral curvature geometric function
                  cp3_poswfn_$1            # Lateral position of contact point on wheel 
                  cp3_posrfn_$1            # Lateral position of contact point on rail
                  0.                       # Wheel/rail prestess force in normal direction of contact surface
                  knwr_$1l                 # Wheel/rail stiffness in normal direction of contact surface
                  cnwr_$1l                 # Wheel/rail damping   in normal direction of contact surface
                  2.05e11                  # The combined modulus of elasticity for body1 and body2
                  .25                      # The combined Poisson's ratio for body1 and body2
                  .95                      # kA
                  0.90                     # kS
                  mu_$1l3                  # The coefficient of friction in the contact surface
                  ro_$1l                   # Wheel rolling radius, will vary for OOR-wheels
                  1e99                     # Vertical radius of rail                                     
#
  func mul    cp1_$1l.wx=  cp1_$1l.nux  cp1_$1l.Vref    # Longitudinal slip speed in contact point
  func mul    cp1_$1r.wx=  cp1_$1r.nux  cp1_$1r.Vref
  func mul    cp2_$1l.wx=  cp2_$1l.nux  cp2_$1l.Vref
  func mul    cp2_$1r.wx=  cp2_$1r.nux  cp2_$1r.Vref
  func mul    cp3_$1l.wx=  cp3_$1l.nux  cp3_$1l.Vref
  func mul    cp3_$1r.wx=  cp3_$1r.nux  cp3_$1r.Vref
#
  func mul    cp1_$1l.wy=  cp1_$1l.nuy  cp1_$1l.Vref    # Lateral slip speed in contact point
  func mul    cp1_$1r.wy=  cp1_$1r.nuy  cp1_$1r.Vref
  func mul    cp2_$1l.wy=  cp2_$1l.nuy  cp2_$1l.Vref
  func mul    cp2_$1r.wy=  cp2_$1r.nuy  cp2_$1r.Vref
  func mul    cp3_$1l.wy=  cp3_$1l.nuy  cp3_$1l.Vref
  func mul    cp3_$1r.wy=  cp3_$1r.nuy  cp3_$1r.Vref
#
  func operp  cp1_$1l.w= sqrt( pow2( cp1_$1l.wx ) + pow2( cp1_$1l.wy ) )
  func operp  cp1_$1r.w= sqrt( pow2( cp1_$1r.wx ) + pow2( cp1_$1r.wy ) )
  func operp  cp2_$1l.w= sqrt( pow2( cp2_$1l.wx ) + pow2( cp2_$1l.wy ) )
  func operp  cp2_$1r.w= sqrt( pow2( cp2_$1r.wx ) + pow2( cp2_$1r.wy ) )
  func operp  cp3_$1l.w= sqrt( pow2( cp3_$1l.wx ) + pow2( cp3_$1l.wy ) )
  func operp  cp3_$1r.w= sqrt( pow2( cp3_$1r.wx ) + pow2( cp3_$1r.wy ) )
#
  no_warning func operp  mu_$1l1= mu.O * ( mu.a * exp( mu.B * cp1_$1l.w ) + mu.A )    # Current coefficient of friction
  no_warning func operp  mu_$1r1= mu.O * ( mu.a * exp( mu.B * cp1_$1r.w ) + mu.A )    # mu= mu0*((1-A)*exp(-Bw) + A)
  no_warning func operp  mu_$1l2= mu.O * ( mu.a * exp( mu.B * cp2_$1l.w ) + mu.A )
  no_warning func operp  mu_$1r2= mu.O * ( mu.a * exp( mu.B * cp2_$1r.w ) + mu.A )
  no_warning func operp  mu_$1l3= mu.O * ( mu.a * exp( mu.B * cp3_$1l.w ) + mu.A )
  no_warning func operp  mu_$1r3= mu.O * ( mu.a * exp( mu.B * cp3_$1r.w ) + mu.A )
#
  s_var sngl cp1_$1r.Vref  s_var sngl mu_$1r1
  s_var sngl cp1_$1r.w     s_var sngl cp1_$1r.wx     s_var sngl cp1_$1r.wy
 ]
  in_substruct wr_coupl [ 111 ]
  in_substruct wr_coupl [ 112 ]
  in_substruct wr_coupl [ 113 ]

##
##  Modal parameters of the wheelsets.
##  ----------------------------------
# mass  m_flex_1 axl_111  70   0.10               ## Torsionsform
# coupl m_flex_1 cp1_111r end_2  0.0636396 0. 0.  # 0. 0.14142 0.
# coupl m_flex_1 cp1_111l end_2  -.0636396 0. 0.  # 0. 0.14142 0.
# s_var sngl axl_111.c1
# s_var sngl axl_111.v1

##
##  Store variables for postprocessing
##  ----------------------------------
#
  s_var scalar_0 mu_
#
  s_var sngl ro_
  s_var sngl bo_
  s_var sngl Vo
#
  s_var var_0 esys.x
  s_var var_0 esys.y
  s_var var_0 esys.z
  s_var var_0 esys.f
  s_var var_0 esys.k
  s_var var_0 esys.p
#
  s_var var_0 esys.vx
  s_var var_0 esys.vy
  s_var var_0 esys.vz
  s_var var_0 esys.vf
  s_var var_0 esys.vk
  s_var var_0 esys.vp
#
  s_var var_0 esys.ax
  s_var var_0 esys.ay
  s_var var_0 esys.az
  s_var var_0 esys.af
  s_var var_0 esys.ak
  s_var var_0 esys.ap
#
  substruct store_axle [
   s_var var_0 axl_$1.x
   s_var var_0 axl_$1.y
   s_var var_0 axl_$1.z
   s_var var_0 axl_$1.f
   s_var var_0 axl_$1.k
   s_var var_0 axl_$1.p
#
   s_var var_0 axl_$1.vx
   s_var var_0 axl_$1.vy
   s_var var_0 axl_$1.vz
   s_var var_0 axl_$1.vf
   s_var var_0 axl_$1.vk
   s_var var_0 axl_$1.vp
#
   s_var var_0 axl_$1.Fx
   s_var var_0 axl_$1.Fy
   s_var var_0 axl_$1.Fz
   s_var var_0 axl_$1.Mf
   s_var var_0 axl_$1.Mk
   s_var var_0 axl_$1.Mp
#
   s_var var_0 axl_$1.Ax
   s_var var_0 axl_$1.Ay
   s_var var_0 axl_$1.Az
   s_var var_0 axl_$1.Af
   s_var var_0 axl_$1.Ak
   s_var var_0 axl_$1.Ap
#
  s_var var_0 trc_$1.x     s_var var_0 trc_$1.y     s_var var_0 trc_$1.z
  s_var var_0 trc_$1.f     s_var var_0 trc_$1.k     s_var var_0 trc_$1.p
  s_var var_0 trc_$1.vx    s_var var_0 trc_$1.vy    s_var var_0 trc_$1.vz
  s_var var_0 trc_$1.vf    s_var var_0 trc_$1.vk    s_var var_0 trc_$1.vp
  ]
#
  substruct store_wheel [
   s_var var_0 ral_$1.y
   s_var var_0 ral_$1.z
   s_var var_0 ral_$1.vy
   s_var var_0 ral_$1.vz
#
   s_var var_0 cp_$1.eta
#
   s_var var_0 tral$1.y
   s_var var_0 tral$1.z
   s_var var_0 tral$1.vy
   s_var var_0 tral$1.vz

   s_var var_0 cp1_$1.z
   s_var var_0 cp1_$1.dr
   s_var var_0 cp1_$1.gam
   s_var var_0 cp1_$1.c
   s_var var_0 cp1_$1.a_b
   s_var var_0 cp1_$1.posw
   s_var var_0 cp1_$1.bo
#
   s_var var_0 cp1_$1.nux
   s_var var_0 cp1_$1.nuy
   s_var var_0 cp1_$1.spin
#
   s_var var_0 cp1_$1.vux
   s_var var_0 cp1_$1.vuy
   s_var var_0 cp1_$1.vuz
   s_var var_0 cp1_$1.Vref
#
   s_var var_0 cp1_$1.Fn
   s_var var_0 cp1_$1.Fny
   s_var var_0 cp1_$1.Fx
   s_var var_0 cp1_$1.Fy
   s_var var_0 cp1_$1.Fz
   s_var var_0 cp1_$1.F1x
   s_var var_0 cp1_$1.F1y
   s_var var_0 cp1_$1.F1z
   s_var var_0 cp1_$1.F2x
   s_var var_0 cp1_$1.F2y
   s_var var_0 cp1_$1.F2z
   s_var var_0 cp1_$1.M1f
   s_var var_0 cp1_$1.M1k
   s_var var_0 cp1_$1.M1p
   s_var var_0 cp1_$1.M2f
   s_var var_0 cp1_$1.M2k
   s_var var_0 cp1_$1.M2p
#
   s_var var_0 cp1_$1.fwk
   s_var var_0 cp1_$1.wArchH
   s_var var_0 cp1_$1.wArchS
  ]
  in_substruct store_axle  [ 111  ]
  in_substruct store_axle  [ 112  ]
  in_substruct store_axle  [ 113  ]
  in_substruct store_wheel [ 111r ]
  in_substruct store_wheel [ 112r ]
  in_substruct store_wheel [ 113r ]
  in_substruct store_wheel [ 111l ]
  in_substruct store_wheel [ 112l ]
  in_substruct store_wheel [ 113l ]
#
# in_substruct store [ 111l ]
# in_substruct store [ 112l ]
# in_substruct store [ 113l ]
#
 eof
