#
#
head 2 "Test of different wheel/rail couplings cwr_coupl= $cwr_coupl"
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= 1
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 tstart 0. tstop 500 tstep 1.00 tout 1.00
# tsim_param= heun_u .5 tstart 0. tstop 2.50 tstep 0.001 tout .005
# tsim_param= heun tstart 0. tstop 2.50 tstep 0.005 tout .005
# tsim_param= heun tstart 0. tstop 5.00 tstep 0.005 tout .005
tsim_param= heun_u .05 .2e-3 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 10 0 20 .005 100 .005
func intpl_r krok_R -200 0 0 0 #.001 .005 100 .005
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.kpfr
# in_substruct kpf_rkona [ 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
# ro_= 0.10 # Wheel radius, all wheels in the train
bo = 0.75 # Semi distance between the wheel profile origins, all wheels in the train
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_ 1887 1887 1887 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.010
# constr fix_rigid_1 axl_$1 z 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-1887*9.81 0. 0. 0.
# coupl p_lin czbap= 0. 500e3
# coupl c czba_$1r axl_$1 0 bo -ro_
# grd_$1 0 bo -ro_
# czbap esys z
# coupl c czba_$1l axl_$1 0 -bo -ro_
# grd_$1 0 -bo -ro_
# czbap esys z
]
#
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.
#
# coupl k_lin ip_coupl axl_111 x 0. 0. -ro_
# trc_111 x 0. 0. -ro_ 17
# coupl p_kfrkc p_name +-`F0 +-`ke +-`Ffm +-`x2 +-`kv +-`c
# coupl p_kfrkc p_rubb 0 100 10 .01 0 1
# coupl kfrkc c_rubb axl_111 0 0 -ro_
# trc_111 0 0 -ro_
# p_rubb esys y
#
#coupl p_lin stifx 0. 500e3
#coupl p_lin stify 0. 500e3
#coupl p_lin stifz -100e3 1e6
#coupl k stiff axl_111 0 0 oper_init( -ro_ - .1 )
# trc_111 0 0 oper_init( -ro_ + .1 )
# stifx esys y
#force rel_lsys1 ipforc axl_111 0. 0. 0. 100 0 0 0 0 0
# coupl k_air3_exp p_airip axl_111 0 0 -ro_
# trc_111 0 0 -ro_
# esys m
# stifx 0.
# stify 0.
# stifz
# 10 .01 10 .01
# 0 1
# 0 1 1 10
##
## Define pure creepage on the three different axles
## -------------------------------------------------
# func operp axl_111.vx= -.05
# func operp axl_111.vx= time * -.25
func operp axl_111.vk= time * -.25 - Vo / ro_
# func operp axl_112.vy= -.25
# func operp axl_112.vy= time * -.25
func operp axl_112.p = time * .1
# func operp axl_113.vp= time * -1.e+2
# func operp axl_113.vx= time * -.75e+2
# func operp axl_113.vp= time * -50
# func operp axl_113.vx= time * -37.5
func operp axl_113.vp= time * -3.4
func operp axl_113.vx= time * -2.55
# func operp axl_113.vp= 0.
# func operp axl_113.vx= 0.
constr fix_rigid_1 axl_111 f 0.
constr fix_rigid_1 axl_112 f 0.
constr fix_rigid_1 axl_113 f 0.
# func operp axl_111.vx= 500 * .5e-3
# func operp axl_112.vy= 500 * .5e-3
# func operp axl_113.vp= 500 * 1.e-1
# func operp axl_113.vx= 500 * .75e-1
##
## Calculate creepage and creep-forces
## -----------------------------------
# func char cwr_coupl= wr_coupl_nl1
# func char cwr_coupl= wr_coupl_ne1
# func char cwr_coupl= wr_coupl_nr1
# func char cwr_coupl= wr_coupl_nra1
# func char cwr_coupl= wr_coupl_pe1
# func char cwr_coupl= wr_coupl_pr1
# func char cwr_coupl= wr_coupl_lt1
func char cwr_coupl= wr_coupl_ts1
if_then_char_init cwr_coupl .eq. wr_coupl_ne1 .or.
cwr_coupl .eq. wr_coupl_nr1 .or.
cwr_coupl .eq. wr_coupl_nra1 # {
#
func const knwr.F0_= 100e3 # all wheels
#
func intpl_r knwr_ # Stiffness Tread-Rail all wheels
-1.57 42e6
0. 180e6
1.57 42e6
func intpl_r cnwr_ # Damping Tread-Rail all wheels
-1.57 60e3
0. 250e3
1.57 60e3
func intpl_r knfr_ # Stiffness Flange-Rail all wheels
-1.57 42e6
0. 180e6
1.57 42e6
func intpl_r cnfr_ # Damping Flange-Rail all wheels
-1.57 60e3
0. 250e3
1.57 60e3
endif # }
#
if_then_char_init cwr_coupl .eq. wr_coupl_nl1
insert file $gensys/calc/insert_files/wr_coupl_nl1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_npol1
insert file $gensys/calc/insert_files/wr_coupl_npol1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_ne1
insert file $gensys/calc/insert_files/wr_coupl_ne1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_nr1
insert file $gensys/calc/insert_files/wr_coupl_nr1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_nra1
insert file $gensys/calc/insert_files/wr_coupl_nra1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
##
##
##
if_then_char_init cwr_coupl .eq. wr_coupl_pe1 .or.
cwr_coupl .eq. wr_coupl_pr1 .or.
cwr_coupl .eq. wr_coupl_pra1 .or.
cwr_coupl .eq. wr_coupl_lt1 # {
#
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 600e3
# 0. 600e3
# 1.57 600e3
func intpl_r knfr_ # Stiffness Flange-Rail all wheels
-1.57 600e6
0. 600e6
1.57 600e6
# func intpl_r cnfr_ # Damping Flange-Rail all wheels
# -1.57 600e3
# 0. 600e3
# 1.57 600e3
func const_block
pknwr = 1 # Linjär styvhet i kontaktpunkten
kyrt = 42e6
kzrt = 180e6
kzrt.F0_= knwr.F0_
cyrt = 205e3
czrt = 250e3
end_block
endif # }
if_then_char_init cwr_coupl .eq. wr_coupl_lt1 # {
func intpl_r cnwr_ # Damping Tread-Rail all wheels
-1.57 60e3
0. 250e3
1.57 60e3
func intpl_r cnfr_ # Damping Flange-Rail all wheels
-1.57 60e3
0. 250e3
1.57 60e3
endif # }
##
if_then_char_init cwr_coupl .eq. wr_coupl_pe1
insert file $gensys/calc/insert_files/wr_coupl_pe1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_pr1
insert file $gensys/calc/insert_files/wr_coupl_pr1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
substruct wr_coupl_lt1 [
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_lookuptable_1 cpt_$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
cpt_1.zfn # Wheel lift geometric function
cpt_1.drfn # Wheel radius geometric function
cpt_1.gamfn # Contact angle geometric function
cpt_1.rofn # Lateral curvature geometric function
cpt_1.poswfn # Lateral position of contact point on wheel
cpt_1.posrfn # 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 in body1 and body2
.3 # The combined Poisson's ratio in body1 and body2
mu_$1r # The coefficient of friction in the contact surface
ro_$1r # Wheel rolling radius, will vary for OOR-wheels
#
coupl creep_lookuptable_1 cpt_$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
cpt_1.zfn # Wheel lift geometric function
cpt_1.drfn # Wheel radius geometric function
cpt_1.gamfn # Contact angle geometric function
cpt_1.rofn # Lateral curvature geometric function
cpt_1.poswfn # Lateral position of contact point on wheel
cpt_1.posrfn # 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 in body1 and body2
.3 # The combined Poisson's ratio in body1 and body2
mu_$1l # The coefficient of friction in the contact surface
ro_$1l # Wheel rolling radius, will vary for OOR-wheels
]
if_then_char_init cwr_coupl .eq. wr_coupl_lt1
in_substruct wr_coupl_lt1 [ 111 ]
in_substruct wr_coupl_lt1 [ 112 ]
in_substruct wr_coupl_lt1 [ 113 ]
endif
#
# Stiffnesses in Tanel Springs
# ------------------------------------------
func const kz_fastsim= 2.1379e+07
func const kz_tanel_1= kz_fastsim*0.2
func const kx_fastsim= 34.45e6
func const kx_tanel_1= kx_fastsim*0.2
func const kx_tanel_2= kx_fastsim*0.15
func const ky_fastsim= 30.48e6
func const ky_tanel_1= ky_fastsim*0.2
func const ky_tanel_2= ky_fastsim*0.15
func const cz_fastsim= 20.
##
## Styvheter för rälens borstar
## ----------------------------
# func const kz_fastsim= 2.1379e+07
# func const kz_inter_1= kz_fastsim*0.2
## func const kx_fastsim= 7.3460e+06
# func const kx_fastsim= 3.1413e+07
# func const kx_fastsim= 2.65e+07
# func const kx_inter_1= kx_fastsim*0.2
# func const kx_inter_2= kx_fastsim*0.15
## func const ky_fastsim= 7.4972e+06
## func const ky_fastsim= 3.1018e+07
# func const ky_fastsim= 2.65e+07
# func const ky_inter_1= ky_fastsim*0.2
# func const ky_inter_2= ky_fastsim*0.15
#
# Hjul och räl profiler
# ----------------------------
func intpl_r rail_prof_right= -.030 .004606079858
-.027 .003713967783
-.024 .002922711204
-.021 .002229861410
-.018 .001633338981
-.015 .001131400334
-.012 .000722610832
-.009 .000405823463
-.006 .000180162292
-.003 .000045010130
.000 .000000000000
.003 .000045010130
.006 .000180162292
.009 .000405823463
.012 .000722610832
.015 .001131400334
.018 .001633338981
.021 .002229861410
.024 .002922711204
.027 .003713967783
.030 .004606079858
func intpl_r wheel_prof_right= -1. -.0#3
1. 0.0#3
func intpl_r rail_prof_left = -.030 .004606079858
-.027 .003713967783
-.024 .002922711204
-.021 .002229861410
-.018 .001633338981
-.015 .001131400334
-.012 .000722610832
-.009 .000405823463
-.006 .000180162292
-.003 .000045010130
.000 .000000000000
.003 .000045010130
.006 .000180162292
.009 .000405823463
.012 .000722610832
.015 .001131400334
.018 .001633338981
.021 .002229861410
.024 .002922711204
.027 .003713967783
.030 .004606079858
func intpl_r wheel_prof_left = -1. -.0#3
1. 0.0#3
substruct wr_coupl_ts1 [
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_tanel_springs_1 cpt_$1r
# axl_$1 0. bo 0.
# ral_$1r whe_$1r.x bo 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
# 0.7 # mulfact_spin_$1r # Spin creep relaxation due to contaminated surfaces
# kz_fastsim # kz_fastsim
# kz_inter_1 # kz_inter_1
# kx_fastsim # cpb_kx
# kx_inter_1 # cpb_kx1
# kx_inter_2 # cpb_kx1
# ky_fastsim # cpb_ky
# ky_inter_1 # cpb_ky1
# ky_inter_2 # cpb_ky1
# 0. # Wheel/rail prestess force in normal direction of contact surface
# 0. # Wheel/rail stiffness in normal direction of contact surface
# 0. # Wheel/rail damping in normal direction of contact surface
# mu_$1r # The coefficient of friction in the contact surface
# ro_$1r # Wheel rolling radius, will vary for OOR-wheels
#
# 21 # Number of brushes on the surface of the rail in longitudinal direction
# .003 # Distance between brushes in longitudinal direction
#
# -1 # Orientation
# rail_prof_right
# wheel_prof_right
#
##
# coupl creep_tanel_springs_1 cpt_$1l
# axl_$1 0. -bo 0.
# ral_$1l whe_$1l.x -bo 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
# 0.7 # mulfact_spin_$1l # Spin creep relaxation due to contaminated surfaces
# kz_fastsim # kz_fastsim
# kz_inter_1 # kz_inter_1
# kx_fastsim # cpb_kx
# kx_inter_1 # cpb_kx1
# kx_inter_2 # cpb_kx1
# ky_fastsim # cpb_ky
# ky_inter_1 # cpb_ky1
# ky_inter_2 # cpb_ky1
# 0. # Wheel/rail prestess force in normal direction of contact surface
# 0. # Wheel/rail stiffness in normal direction of contact surface
# 0. # Wheel/rail damping in normal direction of contact surface
# mu_$1l # The coefficient of friction in the contact surface
# ro_$1l # Wheel rolling radius, will vary for OOR-wheels
#
# 21 # Number of brushes on the surface of the rail in longitudinal direction
# .003 # Distance between brushes in longitudinal direction
#
# 1 # Orientation
# rail_prof_left
# wheel_prof_left
#
#
func tral_interp_spline $1 # Interpolate track irregularities
#
coupl creep_tanel_springs_1 cpt_$1r
axl_$1 0. bo 0.
ral_$1r whe_$1r.x bo 0. esys m
tral$1r.y # Lateral position track irregularity
tral$1r.z # Vertical position track irregularity
tral$1r.vy # Lateral velocity track irregularity
tral$1r.vz # Vertical velocity track irregularity
ro_$1r # Wheel rolling radius, will vary for OOR-wheels
-1 # Orientation
mu_$1r # The coefficient of friction in the contact surface
1. # mulfact_nux_$1r # Longitudinal creep relaxation due to contaminated surfaces
1. # mulfact_nuy_$1r # Lateral creep relaxation due to contaminated surfaces
.78 # mulfact_spin_$1r # Spin creep relaxation due to contaminated surfaces
kz_fastsim # kz_fastsim
kz_tanel_1 # kz_tanel_1
kx_fastsim # cpb_kx
kx_tanel_1 # cpb_kx1
kx_tanel_2 # cpb_kx1
ky_fastsim # cpb_ky
ky_tanel_1 # cpb_ky1
ky_tanel_2 # cpb_ky1
cz_fastsim # Wheel/rail damping in normal direction of contact surface
wheel_prof_right
21 # Number of brushes on the surface of the rail in longitudinal direction
.003 # Distance between brushes in longitudinal direction
rail_prof_right
#
coupl creep_tanel_springs_1 cpt_$1l
axl_$1 0. -bo 0.
ral_$1l whe_$1l.x -bo 0. esys m
tral$1l.y # Lateral position track irregularity
tral$1l.z # Vertical position track irregularity
tral$1l.vy # Lateral velocity track irregularity
tral$1l.vz # Vertical velocity track irregularity
ro_$1l # Wheel rolling radius, will vary for OOR-wheels
1 # Orientation
mu_$1l # The coefficient of friction in the contact surface
1. # mulfact_nux_$1l # Longitudinal creep relaxation due to contaminated surfaces
1. # mulfact_nuy_$1l # Lateral creep relaxation due to contaminated surfaces
.78 # mulfact_spin_$1l # Spin creep relaxation due to contaminated surfaces
kz_fastsim # kz_fastsim
kz_tanel_1 # kz_tanel_1
kx_fastsim # cpb_kx
kx_tanel_1 # cpb_kx1
kx_tanel_2 # cpb_kx1
ky_fastsim # cpb_ky
ky_tanel_1 # cpb_ky1
ky_tanel_2 # cpb_ky1
cz_fastsim # Wheel/rail damping in normal direction of contact surface
wheel_prof_left
21 # Number of brushes on the surface of the rail in longitudinal direction
.003 # Distance between brushes in longitudinal direction
rail_prof_left
]
if_then_char_init cwr_coupl .eq. wr_coupl_ts1
in_substruct wr_coupl_ts1 [ 111 ]
in_substruct wr_coupl_ts1 [ 112 ]
in_substruct wr_coupl_ts1 [ 113 ]
endif
##
## Modal parameters of the wheelsets.
## ----------------------------------
# mass m_flex_1 axl_111 70 0.10 ## Torsionsform
# coupl m_flex_1 cpt_111r end_2 0.0636396 0. 0. # 0. 0.14142 0.
# coupl m_flex_1 cpt_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
## ----------------------------------
#
func print06_init cpt_111r.a_b # Print the initial value of a variable to standard output
#
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 cpt_$1.z
s_var var_0 cpt_$1.dr
s_var var_0 cpt_$1.gam
s_var var_0 cpt_$1.c
s_var var_0 cpt_$1.a_b
s_var var_0 cpt_$1.posw
s_var var_0 cpt_$1.bo
#
s_var var_0 cpt_$1.nux
s_var var_0 cpt_$1.nuy
s_var var_0 cpt_$1.spin
#
s_var var_0 cpt_$1.vx
s_var var_0 cpt_$1.vy
s_var var_0 cpt_$1.vz
s_var var_0 cpt_$1.vf
s_var var_0 cpt_$1.vk
s_var var_0 cpt_$1.vp
s_var var_0 cpt_$1.vux
s_var var_0 cpt_$1.vuy
s_var var_0 cpt_$1.vuz
s_var var_0 cpt_$1.Vref
s_var var_0 cpt_$1.trl_h
#
s_var var_0 cpt_$1.Fn
s_var var_0 cpt_$1.Fny
s_var var_0 cpt_$1.Fx
s_var var_0 cpt_$1.Fy
s_var var_0 cpt_$1.Fz
s_var var_0 cpt_$1.F1x
s_var var_0 cpt_$1.F1y
s_var var_0 cpt_$1.F1z
s_var var_0 cpt_$1.F2x
s_var var_0 cpt_$1.F2y
s_var var_0 cpt_$1.F2z
s_var var_0 cpt_$1.M1f
s_var var_0 cpt_$1.M1k
s_var var_0 cpt_$1.M1p
s_var var_0 cpt_$1.M2f
s_var var_0 cpt_$1.M2k
s_var var_0 cpt_$1.M2p
#
s_var var_0 cpt_$1.fwk
s_var var_0 cpt_$1.wArchH
s_var var_0 cpt_$1.wArchH_1
s_var var_0 cpt_$1.wArchH_2
s_var var_0 cpt_$1.wArchH_3
s_var var_0 cpt_$1.wArchH_4
s_var var_0 cpt_$1.wArchH_5
s_var var_0 cpt_$1.wArchH_6
s_var var_0 cpt_$1.wArchH_7
s_var var_0 cpt_$1.wArchH_8
s_var var_0 cpt_$1.wArchH_9
s_var var_0 cpt_$1.wArchH_10
s_var var_0 cpt_$1.wArchH_11
s_var var_0 cpt_$1.wArchH_12
s_var var_0 cpt_$1.wArchH_13
s_var var_0 cpt_$1.wArchH_14
s_var var_0 cpt_$1.wArchH_15
s_var var_0 cpt_$1.wArchH_16
s_var var_0 cpt_$1.wArchS
s_var var_0 cpt_$1.wArchS_1
s_var var_0 cpt_$1.wArchS_2
s_var var_0 cpt_$1.wArchS_3
s_var var_0 cpt_$1.wArchS_4
s_var var_0 cpt_$1.wArchS_5
s_var var_0 cpt_$1.wArchS_6
s_var var_0 cpt_$1.wArchS_7
s_var var_0 cpt_$1.wArchS_8
s_var var_0 cpt_$1.wArchS_9
s_var var_0 cpt_$1.wArchS_10
s_var var_0 cpt_$1.wArchS_11
s_var var_0 cpt_$1.wArchS_12
s_var var_0 cpt_$1.wArchS_13
s_var var_0 cpt_$1.wArchS_14
s_var var_0 cpt_$1.wArchS_15
s_var var_0 cpt_$1.wArchS_16
#
s_var var_0 cpt_$1.slip_area
]
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 ]
#
initval read_gpdat gp/quasi_ts1.gp 1
#
eof
Calculates initial values for runf/pure_creep_ts1.tsimf
#
#
head 2 "Test of different wheel/rail couplings cwr_coupl= $cwr_coupl"
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= 1
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 tstart 0. tstop 500 tstep 1.00 tout 1.00
# tsim_param= heun_u .5 tstart 0. tstop 2.50 tstep 0.001 tout .005
# tsim_param= heun tstart 0. tstop 2.50 tstep 0.005 tout .005
# tsim_param= heun tstart 0. tstop 5.00 tstep 0.005 tout .005
tsim_param= heun_u .05 .2e-3 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 10 0 20 .005 100 .005
func intpl_r krok_R -200 0 0 0 #.001 .005 100 .005
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.kpfr
# in_substruct kpf_rkona [ 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
# ro_= 0.10 # Wheel radius, all wheels in the train
bo = 0.75 # Semi distance between the wheel profile origins, all wheels in the train
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_ 1887 1887 1887 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.010
# constr fix_rigid_1 axl_$1 z 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-1887*9.81 0. 0. 0.
# coupl p_lin czbap= 0. 500e3
# coupl c czba_$1r axl_$1 0 bo -ro_
# grd_$1 0 bo -ro_
# czbap esys z
# coupl c czba_$1l axl_$1 0 -bo -ro_
# grd_$1 0 -bo -ro_
# czbap esys z
]
#
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.
#
# coupl k_lin ip_coupl axl_111 x 0. 0. -ro_
# trc_111 x 0. 0. -ro_ 17
# coupl p_kfrkc p_name +-`F0 +-`ke +-`Ffm +-`x2 +-`kv +-`c
# coupl p_kfrkc p_rubb 0 100 10 .01 0 1
# coupl kfrkc c_rubb axl_111 0 0 -ro_
# trc_111 0 0 -ro_
# p_rubb esys y
#
#coupl p_lin stifx 0. 500e3
#coupl p_lin stify 0. 500e3
#coupl p_lin stifz -100e3 1e6
#coupl k stiff axl_111 0 0 oper_init( -ro_ - .1 )
# trc_111 0 0 oper_init( -ro_ + .1 )
# stifx esys y
#force rel_lsys1 ipforc axl_111 0. 0. 0. 100 0 0 0 0 0
# coupl k_air3_exp p_airip axl_111 0 0 -ro_
# trc_111 0 0 -ro_
# esys m
# stifx 0.
# stify 0.
# stifz
# 10 .01 10 .01
# 0 1
# 0 1 1 10
##
## Define pure creepage on the three different axles
## -------------------------------------------------
# func operp axl_111.vx= -.05
# func operp axl_111.vx= time * -.25
func operp axl_111.vk= time * -.25 - Vo / ro_
# func operp axl_112.vy= -.25
# func operp axl_112.vy= time * -.25
func operp axl_112.p = time * .1
# func operp axl_113.vp= time * -1.e+2
# func operp axl_113.vx= time * -.75e+2
# func operp axl_113.vp= time * -50
# func operp axl_113.vx= time * -37.5
func operp axl_113.vp= time * -3.4
func operp axl_113.vx= time * -2.55
# func operp axl_113.vp= 0.
# func operp axl_113.vx= 0.
constr fix_rigid_1 axl_111 f 0.
constr fix_rigid_1 axl_112 f 0.
constr fix_rigid_1 axl_113 f 0.
# func operp axl_111.vx= 500 * .5e-3
# func operp axl_112.vy= 500 * .5e-3
# func operp axl_113.vp= 500 * 1.e-1
# func operp axl_113.vx= 500 * .75e-1
##
## Calculate creepage and creep-forces
## -----------------------------------
# func char cwr_coupl= wr_coupl_nl1
# func char cwr_coupl= wr_coupl_ne1
# func char cwr_coupl= wr_coupl_nr1
# func char cwr_coupl= wr_coupl_nra1
# func char cwr_coupl= wr_coupl_pe1
# func char cwr_coupl= wr_coupl_pr1
# func char cwr_coupl= wr_coupl_lt1
func char cwr_coupl= wr_coupl_ts1
if_then_char_init cwr_coupl .eq. wr_coupl_ne1 .or.
cwr_coupl .eq. wr_coupl_nr1 .or.
cwr_coupl .eq. wr_coupl_nra1 # {
#
func const knwr.F0_= 100e3 # all wheels
#
func intpl_r knwr_ # Stiffness Tread-Rail all wheels
-1.57 42e6
0. 180e6
1.57 42e6
func intpl_r cnwr_ # Damping Tread-Rail all wheels
-1.57 60e3
0. 250e3
1.57 60e3
func intpl_r knfr_ # Stiffness Flange-Rail all wheels
-1.57 42e6
0. 180e6
1.57 42e6
func intpl_r cnfr_ # Damping Flange-Rail all wheels
-1.57 60e3
0. 250e3
1.57 60e3
endif # }
#
if_then_char_init cwr_coupl .eq. wr_coupl_nl1
insert file $gensys/calc/insert_files/wr_coupl_nl1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_npol1
insert file $gensys/calc/insert_files/wr_coupl_npol1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_ne1
insert file $gensys/calc/insert_files/wr_coupl_ne1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_nr1
insert file $gensys/calc/insert_files/wr_coupl_nr1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_nra1
insert file $gensys/calc/insert_files/wr_coupl_nra1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
##
##
##
if_then_char_init cwr_coupl .eq. wr_coupl_pe1 .or.
cwr_coupl .eq. wr_coupl_pr1 .or.
cwr_coupl .eq. wr_coupl_pra1 .or.
cwr_coupl .eq. wr_coupl_lt1 # {
#
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 600e3
# 0. 600e3
# 1.57 600e3
func intpl_r knfr_ # Stiffness Flange-Rail all wheels
-1.57 600e6
0. 600e6
1.57 600e6
# func intpl_r cnfr_ # Damping Flange-Rail all wheels
# -1.57 600e3
# 0. 600e3
# 1.57 600e3
func const_block
pknwr = 1 # Linjär styvhet i kontaktpunkten
kyrt = 42e6
kzrt = 180e6
kzrt.F0_= knwr.F0_
cyrt = 205e3
czrt = 250e3
end_block
endif # }
if_then_char_init cwr_coupl .eq. wr_coupl_lt1 # {
func intpl_r cnwr_ # Damping Tread-Rail all wheels
-1.57 60e3
0. 250e3
1.57 60e3
func intpl_r cnfr_ # Damping Flange-Rail all wheels
-1.57 60e3
0. 250e3
1.57 60e3
endif # }
##
if_then_char_init cwr_coupl .eq. wr_coupl_pe1
insert file $gensys/calc/insert_files/wr_coupl_pe1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_pr1
insert file $gensys/calc/insert_files/wr_coupl_pr1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
substruct wr_coupl_lt1 [
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_lookuptable_1 cpt_$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
cpt_1.zfn # Wheel lift geometric function
cpt_1.drfn # Wheel radius geometric function
cpt_1.gamfn # Contact angle geometric function
cpt_1.rofn # Lateral curvature geometric function
cpt_1.poswfn # Lateral position of contact point on wheel
cpt_1.posrfn # 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 in body1 and body2
.3 # The combined Poisson's ratio in body1 and body2
mu_$1r # The coefficient of friction in the contact surface
ro_$1r # Wheel rolling radius, will vary for OOR-wheels
#
coupl creep_lookuptable_1 cpt_$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
cpt_1.zfn # Wheel lift geometric function
cpt_1.drfn # Wheel radius geometric function
cpt_1.gamfn # Contact angle geometric function
cpt_1.rofn # Lateral curvature geometric function
cpt_1.poswfn # Lateral position of contact point on wheel
cpt_1.posrfn # 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 in body1 and body2
.3 # The combined Poisson's ratio in body1 and body2
mu_$1l # The coefficient of friction in the contact surface
ro_$1l # Wheel rolling radius, will vary for OOR-wheels
]
if_then_char_init cwr_coupl .eq. wr_coupl_lt1
in_substruct wr_coupl_lt1 [ 111 ]
in_substruct wr_coupl_lt1 [ 112 ]
in_substruct wr_coupl_lt1 [ 113 ]
endif
#
# Stiffnesses in Tanel Springs
# ------------------------------------------
func const kz_fastsim= 2.1379e+07
func const kz_tanel_1= kz_fastsim*0.2
func const kx_fastsim= 34.45e6
func const kx_tanel_1= kx_fastsim*0.2
func const kx_tanel_2= kx_fastsim*0.15
func const ky_fastsim= 30.48e6
func const ky_tanel_1= ky_fastsim*0.2
func const ky_tanel_2= ky_fastsim*0.15
func const cz_fastsim= 20.
##
## Styvheter för rälens borstar
## ----------------------------
# func const kz_fastsim= 2.1379e+07
# func const kz_inter_1= kz_fastsim*0.2
## func const kx_fastsim= 7.3460e+06
# func const kx_fastsim= 3.1413e+07
# func const kx_fastsim= 2.65e+07
# func const kx_inter_1= kx_fastsim*0.2
# func const kx_inter_2= kx_fastsim*0.15
## func const ky_fastsim= 7.4972e+06
## func const ky_fastsim= 3.1018e+07
# func const ky_fastsim= 2.65e+07
# func const ky_inter_1= ky_fastsim*0.2
# func const ky_inter_2= ky_fastsim*0.15
#
# Hjul och räl profiler
# ----------------------------
func intpl_r rail_prof_right= -.030 .004606079858
-.027 .003713967783
-.024 .002922711204
-.021 .002229861410
-.018 .001633338981
-.015 .001131400334
-.012 .000722610832
-.009 .000405823463
-.006 .000180162292
-.003 .000045010130
.000 .000000000000
.003 .000045010130
.006 .000180162292
.009 .000405823463
.012 .000722610832
.015 .001131400334
.018 .001633338981
.021 .002229861410
.024 .002922711204
.027 .003713967783
.030 .004606079858
func intpl_r wheel_prof_right= -1. -.0#3
1. 0.0#3
func intpl_r rail_prof_left = -.030 .004606079858
-.027 .003713967783
-.024 .002922711204
-.021 .002229861410
-.018 .001633338981
-.015 .001131400334
-.012 .000722610832
-.009 .000405823463
-.006 .000180162292
-.003 .000045010130
.000 .000000000000
.003 .000045010130
.006 .000180162292
.009 .000405823463
.012 .000722610832
.015 .001131400334
.018 .001633338981
.021 .002229861410
.024 .002922711204
.027 .003713967783
.030 .004606079858
func intpl_r wheel_prof_left = -1. -.0#3
1. 0.0#3
substruct wr_coupl_ts1 [
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_tanel_springs_1 cpt_$1r
# axl_$1 0. bo 0.
# ral_$1r whe_$1r.x bo 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
# 0.7 # mulfact_spin_$1r # Spin creep relaxation due to contaminated surfaces
# kz_fastsim # kz_fastsim
# kz_inter_1 # kz_inter_1
# kx_fastsim # cpb_kx
# kx_inter_1 # cpb_kx1
# kx_inter_2 # cpb_kx1
# ky_fastsim # cpb_ky
# ky_inter_1 # cpb_ky1
# ky_inter_2 # cpb_ky1
# 0. # Wheel/rail prestess force in normal direction of contact surface
# 0. # Wheel/rail stiffness in normal direction of contact surface
# 0. # Wheel/rail damping in normal direction of contact surface
# mu_$1r # The coefficient of friction in the contact surface
# ro_$1r # Wheel rolling radius, will vary for OOR-wheels
#
# 21 # Number of brushes on the surface of the rail in longitudinal direction
# .003 # Distance between brushes in longitudinal direction
#
# -1 # Orientation
# rail_prof_right
# wheel_prof_right
#
##
# coupl creep_tanel_springs_1 cpt_$1l
# axl_$1 0. -bo 0.
# ral_$1l whe_$1l.x -bo 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
# 0.7 # mulfact_spin_$1l # Spin creep relaxation due to contaminated surfaces
# kz_fastsim # kz_fastsim
# kz_inter_1 # kz_inter_1
# kx_fastsim # cpb_kx
# kx_inter_1 # cpb_kx1
# kx_inter_2 # cpb_kx1
# ky_fastsim # cpb_ky
# ky_inter_1 # cpb_ky1
# ky_inter_2 # cpb_ky1
# 0. # Wheel/rail prestess force in normal direction of contact surface
# 0. # Wheel/rail stiffness in normal direction of contact surface
# 0. # Wheel/rail damping in normal direction of contact surface
# mu_$1l # The coefficient of friction in the contact surface
# ro_$1l # Wheel rolling radius, will vary for OOR-wheels
#
# 21 # Number of brushes on the surface of the rail in longitudinal direction
# .003 # Distance between brushes in longitudinal direction
#
# 1 # Orientation
# rail_prof_left
# wheel_prof_left
#
#
func tral_interp_spline $1 # Interpolate track irregularities
#
coupl creep_tanel_springs_1 cpt_$1r
axl_$1 0. bo 0.
ral_$1r whe_$1r.x bo 0. esys m
tral$1r.y # Lateral position track irregularity
tral$1r.z # Vertical position track irregularity
tral$1r.vy # Lateral velocity track irregularity
tral$1r.vz # Vertical velocity track irregularity
ro_$1r # Wheel rolling radius, will vary for OOR-wheels
-1 # Orientation
mu_$1r # The coefficient of friction in the contact surface
1. # mulfact_nux_$1r # Longitudinal creep relaxation due to contaminated surfaces
1. # mulfact_nuy_$1r # Lateral creep relaxation due to contaminated surfaces
.78 # mulfact_spin_$1r # Spin creep relaxation due to contaminated surfaces
kz_fastsim # kz_fastsim
kz_tanel_1 # kz_tanel_1
kx_fastsim # cpb_kx
kx_tanel_1 # cpb_kx1
kx_tanel_2 # cpb_kx1
ky_fastsim # cpb_ky
ky_tanel_1 # cpb_ky1
ky_tanel_2 # cpb_ky1
cz_fastsim # Wheel/rail damping in normal direction of contact surface
wheel_prof_right
21 # Number of brushes on the surface of the rail in longitudinal direction
.003 # Distance between brushes in longitudinal direction
rail_prof_right
#
coupl creep_tanel_springs_1 cpt_$1l
axl_$1 0. -bo 0.
ral_$1l whe_$1l.x -bo 0. esys m
tral$1l.y # Lateral position track irregularity
tral$1l.z # Vertical position track irregularity
tral$1l.vy # Lateral velocity track irregularity
tral$1l.vz # Vertical velocity track irregularity
ro_$1l # Wheel rolling radius, will vary for OOR-wheels
1 # Orientation
mu_$1l # The coefficient of friction in the contact surface
1. # mulfact_nux_$1l # Longitudinal creep relaxation due to contaminated surfaces
1. # mulfact_nuy_$1l # Lateral creep relaxation due to contaminated surfaces
.78 # mulfact_spin_$1l # Spin creep relaxation due to contaminated surfaces
kz_fastsim # kz_fastsim
kz_tanel_1 # kz_tanel_1
kx_fastsim # cpb_kx
kx_tanel_1 # cpb_kx1
kx_tanel_2 # cpb_kx1
ky_fastsim # cpb_ky
ky_tanel_1 # cpb_ky1
ky_tanel_2 # cpb_ky1
cz_fastsim # Wheel/rail damping in normal direction of contact surface
wheel_prof_left
21 # Number of brushes on the surface of the rail in longitudinal direction
.003 # Distance between brushes in longitudinal direction
rail_prof_left
]
if_then_char_init cwr_coupl .eq. wr_coupl_ts1
in_substruct wr_coupl_ts1 [ 111 ]
in_substruct wr_coupl_ts1 [ 112 ]
in_substruct wr_coupl_ts1 [ 113 ]
endif
##
## Modal parameters of the wheelsets.
## ----------------------------------
# mass m_flex_1 axl_111 70 0.10 ## Torsionsform
# coupl m_flex_1 cpt_111r end_2 0.0636396 0. 0. # 0. 0.14142 0.
# coupl m_flex_1 cpt_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
## ----------------------------------
#
func print06_init cpt_111r.a_b # Print the initial value of a variable to standard output
#
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 cpt_$1.z
s_var var_0 cpt_$1.dr
s_var var_0 cpt_$1.gam
s_var var_0 cpt_$1.c
s_var var_0 cpt_$1.a_b
s_var var_0 cpt_$1.posw
s_var var_0 cpt_$1.bo
#
s_var var_0 cpt_$1.nux
s_var var_0 cpt_$1.nuy
s_var var_0 cpt_$1.spin
#
s_var var_0 cpt_$1.vx
s_var var_0 cpt_$1.vy
s_var var_0 cpt_$1.vz
s_var var_0 cpt_$1.vf
s_var var_0 cpt_$1.vk
s_var var_0 cpt_$1.vp
s_var var_0 cpt_$1.vux
s_var var_0 cpt_$1.vuy
s_var var_0 cpt_$1.vuz
s_var var_0 cpt_$1.Vref
s_var var_0 cpt_$1.trl_h
#
s_var var_0 cpt_$1.Fn
s_var var_0 cpt_$1.Fny
s_var var_0 cpt_$1.Fx
s_var var_0 cpt_$1.Fy
s_var var_0 cpt_$1.Fz
s_var var_0 cpt_$1.F1x
s_var var_0 cpt_$1.F1y
s_var var_0 cpt_$1.F1z
s_var var_0 cpt_$1.F2x
s_var var_0 cpt_$1.F2y
s_var var_0 cpt_$1.F2z
s_var var_0 cpt_$1.M1f
s_var var_0 cpt_$1.M1k
s_var var_0 cpt_$1.M1p
s_var var_0 cpt_$1.M2f
s_var var_0 cpt_$1.M2k
s_var var_0 cpt_$1.M2p
#
s_var var_0 cpt_$1.fwk
s_var var_0 cpt_$1.wArchH
s_var var_0 cpt_$1.wArchH_1
s_var var_0 cpt_$1.wArchH_2
s_var var_0 cpt_$1.wArchH_3
s_var var_0 cpt_$1.wArchH_4
s_var var_0 cpt_$1.wArchH_5
s_var var_0 cpt_$1.wArchH_6
s_var var_0 cpt_$1.wArchH_7
s_var var_0 cpt_$1.wArchH_8
s_var var_0 cpt_$1.wArchH_9
s_var var_0 cpt_$1.wArchH_10
s_var var_0 cpt_$1.wArchH_11
s_var var_0 cpt_$1.wArchH_12
s_var var_0 cpt_$1.wArchH_13
s_var var_0 cpt_$1.wArchH_14
s_var var_0 cpt_$1.wArchH_15
s_var var_0 cpt_$1.wArchH_16
s_var var_0 cpt_$1.wArchS
s_var var_0 cpt_$1.wArchS_1
s_var var_0 cpt_$1.wArchS_2
s_var var_0 cpt_$1.wArchS_3
s_var var_0 cpt_$1.wArchS_4
s_var var_0 cpt_$1.wArchS_5
s_var var_0 cpt_$1.wArchS_6
s_var var_0 cpt_$1.wArchS_7
s_var var_0 cpt_$1.wArchS_8
s_var var_0 cpt_$1.wArchS_9
s_var var_0 cpt_$1.wArchS_10
s_var var_0 cpt_$1.wArchS_11
s_var var_0 cpt_$1.wArchS_12
s_var var_0 cpt_$1.wArchS_13
s_var var_0 cpt_$1.wArchS_14
s_var var_0 cpt_$1.wArchS_15
s_var var_0 cpt_$1.wArchS_16
#
s_var var_0 cpt_$1.slip_area
]
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 ]
#
# initval read_gpdat gp/quasi_ts1.gp 1
#
eof
Creates a diagram of the creepage creep force relationship
## ## Input data for program MPLOT ## Page ------------------------- page_head 1 Longitudinal creepage ygridint= 1 xgridint= 2 overwrite= global y_mid= auto yint/cm= auto xax_yval= bot x_left= 0 xint/cm= auto Diagram 11 x_left= 0 xint/cm= 1e-3 y_bot = 0 yint/cm= 4e3 Curve var_id= pure_creep_ts1 yvar= cpt_111r.Fx xvar= cpt_111r.nux # Curve var_id= pure_creep_lt1 yvar= cpt_111r.Fx xvar= cpt_111r.nux # Curve var_id= pure_creep_pe1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_pr1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_pra1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_ne1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_nr1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_nra1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_npol1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_nl1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux EndPage Page ------------------------- page_head 1 Lateral creepage ygridint= 1 xgridint= 2 overwrite= global y_mid= auto yint/cm= auto xax_yval= bot x_left= 0 xint/cm= auto Diagram 11 # xint/cm= .0025 # y_bot= 0. yint/cm= 10e3 x_left= 0 xint/cm= 1e-3 y_bot = 0 yint/cm= 4e3 Curve var_id= pure_creep_ts1 yvar= cpt_112r.Fy xvar= cpt_112r.nuy # Curve var_id= pure_creep_lt1 yvar= cpt_112r.Fy xvar= cpt_112r.nuy # Curve var_id= pure_creep_pe1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_pr1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_pra1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_ne1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_nr1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_nra1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_npol1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_nl1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy EndPage Page ------------------------- page_head 1 Spin creepage ygridint= 1 xgridint= 2 overwrite= global y_mid= auto yint/cm= auto xax_yval= bot x_left= 0 xint/cm= auto Diagram 11 # xint/cm= .0025 # y_bot= 0. yint/cm= 10e3 y_bot = 0 yint/cm= 2e3 x_left= 0 xint/cm= .05 Curve var_id= pure_creep_ts1 yvar= cpt_113r.Fy xvar= cpt_113r.spin # Curve var_id= pure_creep_lt1 yvar= cpt_113r.Fy xvar= cpt_113r.spin # Curve var_id= pure_creep_pe1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_pr1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_pra1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_ne1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_nr1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_nra1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_npol1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_nl1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin EndPage