## ## ## Create wheel-rail connection ## Creep forces are calculated in a lookup table, according to ## Kalker's simplified theory. ## --------------------------------------------------------------- ## ## substruct wr_coupl [ # $1= axle # note note note Start of substructure file wr_coupl_pe3.ins Rev.091001 note ---------------------------------------------------------- # # # Set constraints and initial values to axle and track # ---------------------------------------------------- constr fix_free_1 whe_$1r.k = 0. constr fix_free_1 whe_$1l.k = 0. initval set_var whe_$1r.vk= -Vo/ro_$1r initval set_var whe_$1l.vk= -Vo/ro_$1l # 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 wr_coupl_pe3 # $1 # Name/number of the wheelset lsa_$1 # Name of the linear local coordinate system # lat_trac vert_trac spv_trac fi_trac # Track irregularities memory fields YMtrac ZMtrac GMtrac CMtrac # Multiplication factors for track irregularities gauge_average # Average gauge of spv_trac gauge_dev_$1 # Modify average gauge for a different conicity 1 # Type of contact normal to the contact surface # whe_$1r # Body wheel right side whe_$1l # Body wheel left side trc_$1 # Body track ro_$1r # Nominal wheel radius right wheel kyrt_$1r # Lateral stiffness rail - track, right side kzrt_$1r # Vertical stiffness rail - track, right side kzrt.F0_$1r # Vertical prestress force rail - track, right side cyrt_$1r # Lateral damping rail - track, right side czrt_$1r # Vertical damping rail - track, right side bo_$1r # Lateral semi-distance to nominal running circle, right side ro_$1l # Nominal wheel radius left wheel kyrt_$1l # Lateral stiffness rail - track, left side kzrt_$1l # Vertical stiffness rail - track, left side kzrt.F0_$1l # Vertical prestress force rail - track, left side cyrt_$1l # Lateral damping rail - track, left side czrt_$1l # Vertical damping rail - track, left side -bo_$1l # Lateral semi-distance to nominal running circle, left side # cp1_$1r # Name of contact point #1 right side trc_$1 cp1_$1r.ksi cp1_$1r.bo 0. # Contact on track whe_$1r 0. cp1_$1r.bo 0. # Contact on wheel/wheelset mu_cp1$1r # Coefficient of friction 2.05e11 .3 # Modulus of elasticity and Poisson's ratio mulfact_nux_cp1$1r mulfact_nuy_cp1$1r # Creepage reduction due to contaminated rail surface mulfact_spin_cp1$1r knwr.F0_cp1$1r knwr_cp1$1r # Prestress force and stiffness normal to the surface # cp1_$1l trc_$1 cp1_$1l.ksi cp1_$1l.bo 0. # Contact on track whe_$1l 0. cp1_$1l.bo 0. # Contact on wheel/wheelset mu_cp1$1l # Coefficient of friction 2.05e11 .3 # Modulus of elasticity and Poisson's ratio mulfact_nux_cp1$1l mulfact_nuy_cp1$1l # Creepage reduction due to contaminated rail surface mulfact_spin_cp1$1l knwr.F0_cp1$1l knwr_cp1$1l # Prestress force and stiffness normal to the surface # # cp2_$1r trc_$1 cp2_$1r.ksi cp2_$1r.bo 0. # Contact on track whe_$1r 0. cp2_$1r.bo 0. # Contact on wheel/wheelset mu_cp2$1r # Coefficient of friction 2.05e11 .3 # Modulus of elasticity and Poisson's ratio mulfact_nux_cp2$1r mulfact_nuy_cp2$1r # Creepage reduction due to contaminated rail surface mulfact_spin_cp2$1r knwr.F0_cp2$1r knwr_cp2$1r # Prestress force and stiffness normal to the surface # cp2_$1l trc_$1 cp2_$1l.ksi cp2_$1l.bo 0. # Contact on track whe_$1l 0. cp2_$1l.bo 0. # Contact on wheel/wheelset mu_cp2$1l # Coefficient of friction 2.05e11 .3 # Modulus of elasticity and Poisson's ratio mulfact_nux_cp2$1l mulfact_nuy_cp2$1l # Creepage reduction due to contaminated rail surface mulfact_spin_cp2$1l knwr.F0_cp2$1l knwr_cp2$1l # Prestress force and stiffness normal to the surface # # cp3_$1r trc_$1 cp3_$1r.ksi cp3_$1r.bo 0. # Contact on track whe_$1r 0. cp3_$1r.bo 0. # Contact on wheel/wheelset mu_cp3$1r # Coefficient of friction 2.05e11 .3 # Modulus of elasticity and Poisson's ratio mulfact_nux_cp3$1r mulfact_nuy_cp3$1r # Creepage reduction due to contaminated rail surface mulfact_spin_cp3$1r knwr.F0_cp3$1r knwr_cp3$1r # Prestress force and stiffness normal to the surface # cp3_$1l trc_$1 cp3_$1l.ksi cp3_$1l.bo 0. # Contact on track whe_$1l 0. cp3_$1l.bo 0. # Contact on wheel/wheelset mu_cp3$1l # Coefficient of friction 2.05e11 .3 # Modulus of elasticity and Poisson's ratio mulfact_nux_cp3$1l mulfact_nuy_cp3$1l # Creepage reduction due to contaminated rail surface mulfact_spin_cp3$1l knwr.F0_cp3$1l knwr_cp3$1l # Prestress force and stiffness normal to the surface ## ## ## Animate track irregularities in GPLOT ## ====================================== if_then_init .exist. rail_$1r.z no_warning func copy rail_$1r.y= tral$1r.y no_warning func copy rail_$1r.z= tral$1r.z no_warning func copy rail_$1r.k= tral$1r.k no_warning func copy rail_$1r.p= tral$1r.p # no_warning func copy rail_$1l.y= tral$1l.y no_warning func copy rail_$1l.z= tral$1l.z no_warning func copy rail_$1l.k= tral$1l.k no_warning func copy rail_$1l.p= tral$1l.p endif # func print06_all cpt_$1r.nux # func print06_all cpt_$1r.nuy # func print06_all cpt_$1r.spi # func print06_all cpt_$1r.Fnx # func print06_all cpt_$1l.Fn # func print06_all cpt_$1l.Fnx # func print06_all cpt_$1l.Fny # func print06_all cpt_$1r.Fn # func print06_all cpt_$1r.Fny # note note ---------------------------------------------------- note End of substructure file wr_coupl_pe3.ins note ---------------------------------------------------- ]