## ## ## Default vehicle data ## ===================================== ## ## Basic vehicle geometry ## ---------------------- func const aca= 6.0 # Axle longitudinal semi-distance ## ## Vehicle bodies ## -------------- ## Car-body (c) func const hfloor= 1.1 # floor level func const hccg_tara= 1.0 # centre of gravity position, vertical func const mc_tara= 6000 # weight of empty body func const Jfc_tara= mc_tara*(1.5^2+0.3^2)/3 # roll inertia, m*(1.5^2+h^2)/3 func const Jkc_tara= mc_tara*(0.3^2+(aca+2)^2)/3 # pitch inertia, m*(h^2+a^2)/3 func const Jpc_tara= mc_tara*((aca+2)^2+1.5^2)/3 # yaw inertia, m*(a^2+1.5^2)/3 # func const hccg_load= 2.2 # centre of gravity position, vertical func const mc_fact= mc_load/37000 func const Jfc_load= mc_load*(1.5^2+0.3^2)/3 # roll inertia, m*(1.5^2+h^2)/3 func const Jkc_load= mc_load*(0.3^2+(aca+2)^2)/3 # pitch inertia, m*(h^2+a^2)/3 func const Jpc_load= mc_load*((aca+2)^2+1.5^2)/3 # yaw inertia, m*(a^2+1.5^2)/3 # func const mc= mc_tara+mc_load func const hccg= (mc_tara*hccg_tara+mc_load*hccg_load)/mc func const Jfc= Jfc_tara+mc_tara*(hccg-hccg_tara)^2+Jfc_load+mc_load*(hccg-hccg_load)^2 func const Jkc= Jkc_tara+mc_tara*(hccg-hccg_tara)^2+Jkc_load+mc_load*(hccg-hccg_load)^2 func const Jpc= Jpc_tara+Jpc_load ## ## Axles (a) func const ro_= 0.50 # centre of gravity pos., vert. (wheel radius) func const ma= 1000 # mass func const Jfa= 800 # moment of inertia, roll func const Jka= 140 # moment of inertia, pitch func const Jpa= 800 # moment of inertia, yaw ## ## Track pieces (t) func const mt 1000 # mass ## ## Add coupling stiffnesses ## ============================================================== ## ## UIC Linkage ## ----------- 9876543210987654321098765432109876543210987654321012345678901234567890123456789012345678901234567890123456789012345678 ! ! car ! _______ ______________________ _______________________ ____________________ ________ ! \ | \ | | / \ / ! \ | \ | | / | / ! \ | \ | | / | / ! \ | __ \ | | / ___ | / ! \ | /__/___________________________________________________________________________________________\__\ | / ! \ | / / |_______________________________________________________________________________| \ \ | / ! \ | / / |_____________________________________________________________| \ \ | / ! \ | / / \ |_____________________________________________| / \ \ | / ! \ | / / \ |___________________________________| / \ \ | / ! \ \/ / \ |_________________________| / \ \/ / ! \ / / \ | _____________ | / \ \ / ! \ / /\ \ | / \ | / / \ \ / ! \ |__/ | \ | / \ | / | \__| / ! \ | \ | | | | / | / ! \______/ \ | | | | / \______/ ! \ | | axl | | / ! \ | | | | / ! \ | \ / | / ! \ | \_____________/ | / ! \ | | / ! \__| |__/ ! ! ! ## ## Lateral semi-distance ## ---------------------------------------------------------------- func const bkzca= 1.0 ## ## ## Longitudinal Stiffness data coupling type k+kf ## ---------------------------------------------------------------- func const_block kxca = 78e3+422e3*mc_fact FxcafrA_= 2*(2000+9000*mc_fact) # Friction force, longitudinal end_block coupl p_lin kxca_ = 0. kxca coupl p_lin kfxca_= 0. 7.5*kxca coupl p_lin cxca_ = 0. 2*0.05*sqrt(mc/4*kxca) # Material viscous damping ## ## Lateral Stiffness data coupling type kf_exp2 ## ---------------------------------------------------------------- func const_block kyca = 123e3+377e3*mc_fact FycafrA_= 2*(2450+7750*mc_fact) # Friction force, lateral end_block coupl p_lin kyca_ = 0. kyca coupl p_lin kfyca_= 0. 25*kyca coupl p_lin cyca_ = 0. 2*0.05*sqrt(mc/4*kyca) # Material viscous damping ## ## Vertical Stiffness data coupling type kf_exp2 ## ---------------------------------------------------------------- func const_block kzca = 1600e3 FzcafrA_= 2*(3000+7600*mc_fact) # Asymptotic vertical friction force FfA FzcafrB_= -FzcafrA_ # Asymptotic vertical friction force FfB kzcafrA_= 20*(3200e3+1600e3*mc_fact) # Initial friction stiffness kfA kzcafrB_= 20*(3200e3+1600e3*mc_fact) # Initial friction stiffness kfB end_block coupl p_lin kzca_= -mc_tara/4*9.81 kzca coupl p_lin czca_= 0. 2*0.05*sqrt(mc/4*kzca) # Material viscous damping ## ## Longitudinal Stops +- 20 mm play ## ---------------------------------------------------------------- func const bkxcast= 1. coupl p_nlin kxcast_ 0. -4.00000E-02 -4.00000E+06 -2.20000E-02 -4.00000E+04 -2.00000E-02 0.00000E+00 0.00000E+00 0.00000E+00 2.00000E-02 0.00000E+00 2.20000E-02 4.00000E+04 4.00000E-02 4.00000E+06 ## ## Lateral Stops +- 20 mm play ## ---------------------------------------------------------------- func const bkycast= 1. coupl p_nlin kycastl_ 0. -1.00000E-01 0.00000E+00 0.00000E+00 0.00000E+00 2.00000E-02 0.00000E+00 2.20000E-02 1.68000E+04 4.00000E-02 5.50000E+04 # coupl p_nlin kycastr_ 0. -4.00000E-02 -5.50000E+04 -2.20000E-02 -1.68000E+04 -2.00000E-02 0.00000E+00 1.00000E-01 0.00000E+00 ## ## Contact point: The Normal Problem Linear stiffness ## --------------------------------------------------- ! ! ! | | ! | whe | ! | | ! | | ! \___________________________________ | ! | \ | ! | / \ | ! / | \ | ! knwr_ \ _| \ / ! / | \_______/ ! \ _| ! | / knfr_ ! ____|____/ ! / \ ! | ral | ! | | ! |___ __| ! \ / ! | | ! | | ! | | ! | | ! ______/ \______ ! |_________________| ! ! func const pknwr= 1 # Indicator for linear stiffness func const knwr.F0_= (mc+2*ma)*9.81/4 # func intpl_r knwr_= -1.57 2400e6 # Stiffness Tread-Rail all wheels 0. 2400e6 # The stiffness can vary depending on contact angle 1.57 2400e6 # func intpl_r knfr_= -1.57 2400e6 # Stiffness Flange-Rail all wheels 0. 2400e6 # The stiffness can vary depending on contact angle 1.57 2400e6 ## ## Rail fastners: Concrete sleepers ## -------------------------------- ! kyrt_ _______ ! _______ |----/\/\/\/----| | | ! / \ | | | | ! [ ral ]--------| |-------| | ! [__ __] | |------ | | | ! | | |----| |------| | | ! | | |------ | | ! | | cyrt_ | | ! | | | | ! _____/ \_____ | | ! |_______________| | | ! | | | ! ____|____ | | ! | | | | ! |_|_| \ | | ! czrt_ | | / kzrt_ | | ! |___| \ | | ! | / | | ! |_________| | | ! | | | ! ________________________|____________________________________| |_______________ ! ! ! ! trc ! ! ! _____________________________________________________________________________________ ! ! func const kzrt.F0_= (mc+2*ma)*9.81/4 # func const kyrt_= 42e6 func const cyrt_= 400e3 func const kzrt_= 75e6 func const czrt_=1600e3 ## ## Ballast stiffness: Typical ## ------------------------------------------------------- ! ! __ __ ! ( ) ( ) |/ ! __________||___________________________________||__________ kytg_ |/ ! | | |----/\/\/\/----| |/ ! | | | | |/ ! | trc_ |--------| |-------|/ grd_ ! | | | |------ | |/ ! |___________________________________________________________| |----| |------| |/ ! |------ |/ ! cytg_ |/ coupl p_lin kytg_= 0 30e6 coupl p_lin cytg_= 0 300e3