Home    Search    Reference Manuals    Return   


##
##
## Vehicle properties Master file of a four axle bogie vehicle
## ===========================================================

##
##  Basic vehicle geometry
##  ----------------------
  func const acb_=    6.5    # bogie pivot semi-distance (longitudinal)
  func const aba_=    1.5    # wheelset semi-distance within a bogie (long.)

##
##  Vehicle bodies
##  --------------
## Car-body (c)
  func const accg_=  0.0     # centre of gravity position, longitudinal
  func const hccg_=  2.0     # centre of gravity position, vertical
  func const mc_=   52000                             # mass
  func const Jfc_=  mc_*(1.50^2+1.538^2)/3            # moment of inertia, roll
  func const Jkc_=  mc_*(1.50^2+(acb_+aba_+.58)^2)/3  # moment of inertia, pitch
  func const Jpc_=  Jkc_                              # moment of inertia, yaw
  func const hfloor= 1.3     # floor level
##
## Bogie frames (b)
  func const hbcg_=   0.7    # centre of gravity position, vertical
  func const mb_ =  10000    # mass
  func const Jfb_=   3000    # moment of inertia, roll
  func const Jkb_=  10000    # moment of inertia, pitch
  func const Jpb_=  15000    # moment of inertia, yaw
##
## Axles (a)
  func const ro_ =   0.50    # centre of gravity pos., vert. (wheel radius)
  func const ma_ =   2000    # mass
  func const Jfa_=   1200    # moment of inertia, roll
  func const Jka_=    250    # moment of inertia, pitch
  func const Jpa_=   1200    # moment of inertia, yaw
##
## Track pieces (t)
  func const myt_= 1000
  func const mzt_= 1000
  func const Jft_= 1000

##
##  Add coupling stiffnesses
##  ==============================================================


##
## Secondary suspension: Coil springs between car and bogie
## ------------------------------------------------------------
  func  const  kzcb.B=  1.0             # lateral semi distance [m]
  func  const  kzcb.H=  1.0             # top of spring [m]
  func  const  kzcb.hs= 0.32            # height of spring [m]
  coupl p_lin  kxcb_=   0.0     600e3   # stiffness, longitudinal shear [N/m]
  coupl p_lin  kycb_=   0.0     600e3   # stiffness, lateral shear [N/m]
  coupl p_lin  kzcb_11= -(mc_*(acb_1+accg_1)/(acb_1+acb_1))*9.81/2  900e3  # stiffness, vertical compression [N/m]
  coupl p_lin  kzcb_12= -(mc_*(acb_1-accg_1)/(acb_1+acb_1))*9.81/2  900e3  # stiffness, vertical compression [N/m]

##
## Secondary suspension: Anti-roll bars
## -----------------------------------------------------
  func const  kfcb.H= 1.0       # Height above top of rail [m]
  coupl p_lin kfcb_ = 0. 2.5e6  # Roll stiffness [N/m]

##
## Secondary suspension: Traction rods
## -----------------------------------------------------
  func  const ktr.Ac= 7.5       # rod position in car-body, longitudinal
  func  const ktr.Ab= -.2       # rod position in bogie,   longitudinal
  func  const ktr.Bc= 0.0       # rod position in car-body, lateral
  func  const ktr.Bb= 0.0       # rod position in bogie,   lateral
  func  const ktr.Hc= 0.4       # rod position in car-body, vertical
  func  const ktr.Hb= 0.3       # rod position in bogie,   vertical
  coupl p_lin ktr_=   0.  25e6  # stiffness
  coupl p_lin ctr_=   0. 100e3  # parallel viscous damping

##
## Secondary suspension: Lateral bumpstops
## -----------------------------------------------------
  func  const    kycbs.H=   1.0         # Height above top of rail [m]
  coupl p_nlin_s kycbs_= 0. 0.030,0.0   # Breakpoint #1 [m],[N]
                            0.040,15e3  # Breakpoint #2 [m],[N]
                            0.090,2e6   # Breakpoint #3 [m],[N]

##
## Secondary suspension: Vertical bumpstops
## -----------------------------------------------------
  func  const    kzcbs.A=   0.0         # Longitudinal distance    [m]
  func  const    kzcbs.B=   1.0         # Lateral semi-distance    [m]
  func  const    kzcbs.H=   1.0         # Height above top of rail [m]
  coupl p_nlin_s kzcbs_= 0. 0.025,0.0   # Breakpoint #1 [m],[N]
                            0.030,15e3  # Breakpoint #2 [m],[N]
                            0.080,2e6   # Breakpoint #3 [m],[N]

##
## Secondary suspension: Lateral viscous dampers
## --------------------------------------------------------
  func const    cycb.A = 0.34           # damper position, longitudinal
  func const    cycb.Bc= 0.9            # position in car-body, lateral
  func const    cycb.Bb= 0.4            # position in bogie,   lateral
  func const    cycb.Hc= 0.7            # position in car-body, vertical
  func const    cycb.Hb= 0.6            # position in bogie,   vertical
  coupl p_lin   cycb_=   0.0  40e3      # damping coefficient

##
## Secondary suspension: Vertical viscous dampers
## --------------------------------------------------------
  func const    czcb.A = 0.34           # damper position, longitudinal
  func const    czcb.B = 1.25           # damper position, lateral
  func const    czcb.Hc= 1.0            # position in car-body, vertical
  func const    czcb.Hb= 0.5            # position in bogie,   vertical
  coupl p_lin   czcb_=   0.0  40e3      # damping coefficient

##
## Secondary suspension: Yaw viscous dampers
## --------------------------------------------------------
  func const    cccb.Ac=  7.3       # position in car-body, longitudinal
  func const    cccb.Ab= -0.4       # position in bogie,   longitudinal
  func const    cccb.Bc=  1.3       # position in car-body, lateral
  func const    cccb.Bb=  1.3       # position in bogie,   lateral
  func const    cccb.Hc=  0.7       # position in car-body, vertical
  func const    cccb.Hb=  0.6       # position in bogie,   vertical
  coupl p_lin   kccb_= 0. 25e6      # series stiffness
  coupl p_nlin  cccb_= 0. -1.032 -26e3  # Blow-off compression
                          -0.032 -16e3  # Damping coeff. compression
                           0.      0.
                           0.032  16e3  # Damping coeff. expansion
                           1.032  26e3  # Blow-off expansion

##
##  Primary suspension: springs
##  ----------------------------------------------------------------------------
  func const  kmba.B=   1.0   # spring position, lateral
#
  func const  kmba.F0_11= (2*kzcb_11.F0-mb_11*9.81)/4
  coupl p_lin36 kmba_11=  0.    0. kmba.F0_11 0.  0.    0.
                         20e6   0.    0.    0.    0.    0.
                          0.   20e6   0.    0.    0.    0.
                          0.    0. 1200e3   0.    0.    0.
                          0.    0.    0.    0.    0.    0.
                          0.    0.    0.    0.    0.    0.
                          0.    0.    0.    0.    0.    0.
#
  func const  kmba.F0_12= (2*kzcb_12.F0-mb_12*9.81)/4
  coupl p_lin36 kmba_12=  0.    0. kmba.F0_12 0.  0.    0.
                         20e6   0.    0.    0.    0.    0.
                          0.   20e6   0.    0.    0.    0.
                          0.    0. 1200e3   0.    0.    0.
                          0.    0.    0.    0.    0.    0.
                          0.    0.    0.    0.    0.    0.
                          0.    0.    0.    0.    0.    0.
#
  coupl p_lin36 cmba_=    0.    0.    0.    0.    0.    0.
                         20e3   0.    0.    0.    0.    0.
                          0.   20e3   0.    0.    0.    0.
                          0.    0.    5e3   0.    0.    0.
                          0.    0.    0.    0.    0.    0.
                          0.    0.    0.    0.    0.    0.
                          0.    0.    0.    0.    0.    0.

##
##  Primary suspension: Lateral bumpstops
##  ----------------------------------------------------------------------------
  func  const    kybas.H= 0.4           # Height above top of rail [m]
  coupl p_nlin_s kybas_=  0.            # symmetric non-linear stiffness
                          0.025 0.0     # breakpoint #1
                          0.050 1e6     # breakpoint #2

##
##  Primary suspension: Vertical bumpstops
##  ----------------------------------------------------------------------------
  func  const    kzbas.A= 0.0           # Longitudinal distance    [m]
  func  const    kzbas.B= 1.1           # Lateral semi-distance    [m]
  func  const    kzbas.H= 0.4           # Height above top of rail [m]
  coupl p_nlin_s kzbas_=  0.            # symmetric non-linear stiffness
                          0.025 0.0     # breakpoint #1
                          0.050 1e6     # breakpoint #2

##
##  Primary suspension: Vertical viscous damper
##  ----------------------------------------------------------------------------
  func const  czba.A = 0.0      # damper position, longitudinal
  func const  czba.Bb= 1.1      # damper position, lateral in bogie
  func const  czba.Ba= 1.1      # damper position, lateral in axle
  func const  czba.Hb= 1.0      # damper position, vertical in bogie
  func const  czba.Ha= ro_      # damper position, vertical in axle
  coupl p_lin czba_=   0. 60e3  # viscous damping