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