Home
Reference Manuals
Return
##
## Input data for program MPLOT
##
##
## In nedit, mark the file
## $gensys/homepage/doc_html/mplot.html
## and select view_html with mouse button 3
## -----------------------------------------
##
ilaser= 6 # Create plots in postscript color
Create_scalar new Xeval_start= 150
Create_scalar new Xeval_stop = 500
Create_scalar new Xplot_start= 0
##
## Ride Comfort according to Ride Index Wz.
## ----------------------------------------
Ftwz car1b1.ay FTname=car1b1.ay.ft ityp=lpv tstart=Xeval_start tstop=Xeval_stop Tsname=lsb_11.pn
Ftwz car1.m.ay FTname=car1.m.ay.ft ityp=lpv tstart=Xeval_start tstop=Xeval_stop Tsname=lsc_1.pn
Ftwz car1b2.ay FTname=car1b2.ay.ft ityp=lpv tstart=Xeval_start tstop=Xeval_stop Tsname=lsb_12.pn
#
Ftwz car1b1.az FTname=car1b1.az.ft ityp=vpv tstart=Xeval_start tstop=Xeval_stop Tsname=lsb_11.pn
Ftwz car1.m.az FTname=car1.m.az.ft ityp=vpv tstart=Xeval_start tstop=Xeval_stop Tsname=lsc_1.pn
Ftwz car1b2.az FTname=car1b2.az.ft ityp=vpv tstart=Xeval_start tstop=Xeval_stop Tsname=lsb_12.pn
##
## Filter acceleration according to ERRI Question B 153 Report no.18.
## ------------------------------------------------------------------
Trans car1b1.ax Fname= car1b1.ax.ER type= ERRI153_WD
Trans car1.m.ax Fname= car1.m.ax.ER type= ERRI153_WD
Trans car1b2.ax Fname= car1b2.ax.ER type= ERRI153_WD
Stat car1b1.ax.ER tstart= Xeval_start tstop=Xeval_stop Tname= lsb_11.pn
Stat car1.m.ax.ER tstart= Xeval_start tstop=Xeval_stop Tname= lsc_1.pn
Stat car1b2.ax.ER tstart= Xeval_start tstop=Xeval_stop Tname= lsb_12.pn
#
Trans car1b1.ay Fname= car1b1.ay.ER type= ERRI153_WD
Trans car1.m.ay Fname= car1.m.ay.ER type= ERRI153_WD
Trans car1b2.ay Fname= car1b2.ay.ER type= ERRI153_WD
Stat car1b1.ay.ER tstart= Xeval_start tstop=Xeval_stop Tname= lsb_11.pn
Stat car1.m.ay.ER tstart= Xeval_start tstop=Xeval_stop Tname= lsc_1.pn
Stat car1b2.ay.ER tstart= Xeval_start tstop=Xeval_stop Tname= lsb_12.pn
#
Trans car1b1.az Fname= car1b1.az.ER type= ERRI153_WB
Trans car1.m.az Fname= car1.m.az.ER type= ERRI153_WB
Trans car1b2.az Fname= car1b2.az.ER type= ERRI153_WB
Stat car1b1.az.ER tstart= Xeval_start tstop=Xeval_stop Tname= lsb_11.pn
Stat car1.m.az.ER tstart= Xeval_start tstop=Xeval_stop Tname= lsc_1.pn
Stat car1b2.az.ER tstart= Xeval_start tstop=Xeval_stop Tname= lsb_12.pn
##
## Calculate max values for pantograph sway
## -----------------------------------------------------------
Stat car1b1.y tstart= Xeval_start tstop=Xeval_stop Tname= lsb_11.pn Ptile= .15 25 50 75 99.85
Stat car1b2.y tstart= Xeval_start tstop=Xeval_stop Tname= lsb_12.pn Ptile= .15 25 50 75 99.85
##
## Vehicle turnover according to the intercept method.
## ---------------------------------------------------
Func operp Qdiff= Q111l + Q112l - Q111r - Q112r
Func abs Qdiffa= Qdiff
Func operp Qsum= Q111l + Q112l + Q111r + Q112r
Func operp bt_unfilt= Qdiffa / Qsum
Filt lpass2_0 1.5 0.7071 bt_unfilt bt
Stat bt tstart= Xeval_start tstop=Xeval_stop Tname= lsb_11.pn Ptile= .15 25 50 75 99.85
##
## Calculate quasistatical Y-forces
## -----------------------------------------------------------------
Stat Y111l tstart= Xeval_start tstop=Xeval_stop Tname= lsa_111.pn Ptile= .15 25 50 75 99.85
Stat Y111r tstart= Xeval_start tstop=Xeval_stop Tname= lsa_111.pn Ptile= .15 25 50 75 99.85
Stat Y112l tstart= Xeval_start tstop=Xeval_stop Tname= lsa_112.pn Ptile= .15 25 50 75 99.85
Stat Y112r tstart= Xeval_start tstop=Xeval_stop Tname= lsa_112.pn Ptile= .15 25 50 75 99.85
##
## Calculate quasistatical Q-forces
## -----------------------------------------------------------------
Stat Q111l tstart= Xeval_start tstop=Xeval_stop Tname= lsa_111.pn Ptile= .15 25 50 75 99.85
Stat Q111r tstart= Xeval_start tstop=Xeval_stop Tname= lsa_111.pn Ptile= .15 25 50 75 99.85
Stat Q112l tstart= Xeval_start tstop=Xeval_stop Tname= lsa_112.pn Ptile= .15 25 50 75 99.85
Stat Q112r tstart= Xeval_start tstop=Xeval_stop Tname= lsa_112.pn Ptile= .15 25 50 75 99.85
##
## Filter the vertical track forces with a 4:order Butterworthfilter
## -----------------------------------------------------------------
substruct Buttw_4 [
filt lpass2_0 1.5 0.3827 $1 $1_f2
filt lpass2_0 1.5 0.9239 $1_f2 $1_f4
stat $1_f4 tstart=Xeval_start tstop=Xeval_stop Tname=$2 ]
#
in_substruct Buttw_4 [ Q111l lsa_111.pn ]
in_substruct Buttw_4 [ Q111r lsa_111.pn ]
in_substruct Buttw_4 [ Q112l lsa_112.pn ]
in_substruct Buttw_4 [ Q112r lsa_112.pn ]
##
## Lowpass the track-shift force in a sliding 2m-min-window
## ------------------------------------------------------------
Filt min 2. S111 Fname=S111.2m Tname=lsa_111.pn
Filt min 2. S112 Fname=S112.2m Tname=lsa_112.pn
Filt min 2. S121 Fname=S121.2m Tname=lsa_121.pn
Filt min 2. S122 Fname=S122.2m Tname=lsa_122.pn
Stat S111.2m tstart=Xeval_start tstop=Xeval_stop Tname=lsa_111.pn Ptile= .15 25 50 75 99.85
Stat S112.2m tstart=Xeval_start tstop=Xeval_stop Tname=lsa_112.pn Ptile= .15 25 50 75 99.85
Stat S121.2m tstart=Xeval_start tstop=Xeval_stop Tname=lsa_121.pn Ptile= .15 25 50 75 99.85
Stat S122.2m tstart=Xeval_start tstop=Xeval_stop Tname=lsa_122.pn Ptile= .15 25 50 75 99.85
Stat S111 tstart=Xeval_start tstop=Xeval_stop Tname=lsa_111.pn Ptile= .15 25 50 75 99.85
Stat S112 tstart=Xeval_start tstop=Xeval_stop Tname=lsa_112.pn Ptile= .15 25 50 75 99.85
Stat S121 tstart=Xeval_start tstop=Xeval_stop Tname=lsa_121.pn Ptile= .15 25 50 75 99.85
Stat S122 tstart=Xeval_start tstop=Xeval_stop Tname=lsa_122.pn Ptile= .15 25 50 75 99.85
#
Func operp S111.2m.max= S111.2mMED + 3 * S111.2mSTD
Func operp S112.2m.max= S112.2mMED + 3 * S112.2mSTD
Func operp S121.2m.max= S121.2mMED + 3 * S121.2mSTD
Func operp S122.2m.max= S122.2mMED + 3 * S122.2mSTD
#
Func max S_med3std= S111.2m.max S112.2m.max S121.2m.max S122.2m.max
##
## Lowpass the flange climbing ratio force by a sliding 50ms-average-window
## ------------------------------------------------------------------------
Filt mean .050 Y/Q111l Fname=Y/Q111l.50
Filt mean .050 Y/Q112l Fname=Y/Q112l.50
Filt mean .050 Y/Q111r Fname=Y/Q111r.50
Filt mean .050 Y/Q112r Fname=Y/Q112r.50
Stat Y/Q111l.50 tstart=Xeval_start tstop=Xeval_stop Tname=lsa_111.pn Ptile= .15 25 50 75 99.85
Stat Y/Q111r.50 tstart=Xeval_start tstop=Xeval_stop Tname=lsa_111.pn Ptile= .15 25 50 75 99.85
Stat Y/Q112l.50 tstart=Xeval_start tstop=Xeval_stop Tname=lsa_112.pn Ptile= .15 25 50 75 99.85
Stat Y/Q112r.50 tstart=Xeval_start tstop=Xeval_stop Tname=lsa_112.pn Ptile= .15 25 50 75 99.85
#
Func max Y/Qmax= Y/Q111l.50MAX Y/Q111r.50MAX Y/Q112l.50MAX Y/Q112r.50MAX
##
## Print the results to file *.print
## -----------------------------------------------------------
print scalar S_med3std Y/Qmax
car1b2.yMAX car1b2.yMIN btMIN
car1b1.ay.ERRMS car1.m.ay.ERRMS car1b2.ay.ERRMS
car1b1.az.ERRMS car1.m.az.ERRMS car1b2.az.ERRMS
car1b1.ayWZ car1.m.ayWZ car1b2.ayWZ
car1b1.azWZ car1.m.azWZ car1b2.azWZ
# NMV_b1 NMV_.m NMV_b2
##
## Plot the accelerations in the car-body
## -----------------------------------------------------------
Page
xvar= lsa_111.pn
x_left= Xplot_start
Diagram 11
head1= 'Curve_1= car1.m.ay Wz=$car1.m.ayWZ RMS=$car1.m.ay.ERRMS'
head2= 'Curve_2= car1b1.ay Wz=$car1b1.ayWZ RMS=$car1b1.ay.ERRMS'
head3= 'Curve_3= car1b2.ay Wz=$car1b2.ayWZ RMS=$car1b2.ay.ERRMS'
Curve yvar= car1.m.ay
Curve yvar= car1b1.ay
Curve yvar= car1b2.ay
EndDiagram
Diagram 21
head1= 'Curve_1= car1.m.az Wz=$car1.m.azWZ RMS=$car1.m.az.ERRMS'
head2= 'Curve_2= car1b1.az Wz=$car1b1.azWZ RMS=$car1b1.az.ERRMS'
head3= 'Curve_3= car1b2.az Wz=$car1b2.azWZ RMS=$car1b2.az.ERRMS'
Curve yvar= car1.m.az
Curve yvar= car1b1.az
Curve yvar= car1b2.az
EndDiagram
EndPage
##
## Plot the displacements of the car-body
## -----------------------------------------------------------
Page
xvar= lsa_111.pn
x_left= Xplot_start
Diagram 11
Curve yvar= car_1.x
Diagram 12
Curve yvar= car_1.y
Diagram 13
Curve yvar= car_1.z
Diagram 21
Curve yvar= car_1.f
Diagram 22
Curve yvar= car_1.k
Diagram 23
Curve yvar= car_1.p
EndPage
##
## Plot the displacements of the bogies
## -----------------------------------------------------------
Page
xvar= lsa_111.pn
x_left= Xplot_start
Diagram 11
Curve yvar= bog_11.x
Curve yvar= bog_12.x
Diagram 12
Curve yvar= bog_11.y
Curve yvar= bog_12.y
Diagram 13
Curve yvar= bog_11.z
Curve yvar= bog_12.z
Diagram 21
Curve yvar= bog_11.f
Curve yvar= bog_12.f
Diagram 22
Curve yvar= bog_11.k
Curve yvar= bog_12.k
Diagram 23
Curve yvar= bog_11.p
Curve yvar= bog_12.p
EndPage
##
## Plot the displacements of the wheelsets in leading bogie
## -----------------------------------------------------------
Page
xvar= lsa_111.pn
x_left= Xplot_start
Diagram 11
Curve yvar= axl_111.x
Curve yvar= axl_112.x
Diagram 12
Curve yvar= axl_111.y
Curve yvar= axl_112.y
Diagram 13
Curve yvar= axl_111.z
Curve yvar= axl_112.z
Diagram 21
Curve yvar= axl_111.f
Curve yvar= axl_112.f
Diagram 22
Curve yvar= axl_111.k
Curve yvar= axl_112.k
Diagram 23
Curve yvar= axl_111.p
Curve yvar= axl_112.p
EndPage
##
## Plot Y- and Q- forces
## -----------------------------------------------------------
Page
xvar= lsa_111.pn
x_left= Xplot_start
Diagram 11
head1= 'Wheel guiding forces, first wheelset'
head2= 'Curve_1= Y111l Max=$Y111lMAX Mean=$Y111lMED P99.85=$FDY111lP5'
head3= 'Curve_2= Y111r Max=$Y111rMAX Mean=$Y111rMED P99.85=$FDY111rP5'
yint/cm= 20e3
Curve yvar= Y111l
Curve yvar= Y111r
EndDiagram
#
Diagram 21
head1= 'Wheel guiding forces, second wheelset'
head2= 'Curve_1= Y112l Max=$Y112lMAX Mean=$Y112lMED P99.85=$FDY112lP5'
head3= 'Curve_2= Y112r Max=$Y112rMAX Mean=$Y112rMED P99.85=$FDY112rP5'
yint/cm= 20e3
Curve yvar= Y112l
Curve yvar= Y112r
EndDiagram
#
Diagram 12
head1= 'Filtered vertical wheel forces, first wheelset'
head2= 'Curve_1= Q111l Max=$Q111lMAX Mean=$Q111lMED P99.85=$FDQ111lP5'
head3= 'Curve_2= Q111r Max=$Q111rMAX Mean=$Q111rMED P99.85=$FDQ111rP5'
yint/cm= 20e3 y_mid= 100e3
Curve yvar= Q111l_f4
Curve yvar= Q111r_f4
EndDiagram
#
Diagram 22
head1= 'Filtered vertical wheel forces, second wheelset'
head2= 'Curve_1= Q112l Max=$Q112lMAX Mean=$Q112lMED P99.85=$FDQ112lP5'
head3= 'Curve_2= Q112r Max=$Q112rMAX Mean=$Q112rMED P99.85=$FDQ112rP5'
yint/cm= 20e3 y_mid= 100e3
Curve yvar= Q112l_f4
Curve yvar= Q112r_f4
EndDiagram
#
EndPage
##
## Plot the lateral track forces
## -----------------------------------------------------------
Page
xvar= lsa_111.pn
x_left= Xplot_start
Diagram 11
head1= 'Track-shift forces, filtered in a sliding 2m-min-window, first bogie'
head2= 'Curve_1= S111.2m Max=$S111.2mMAX Mean=$S111.2mMED P99.85=$FDS111.2mP5'
head3= 'Curve_2= S112.2m Max=$S112.2mMAX Mean=$S112.2mMED P99.85=$FDS112.2mP5'
draw_limit= 75e3 # limit_expl= 'Max allowed track-shift force 75e3'
y_bot= -20e3 yint/cm= 20e3
Curve yvar= S111.2m
Curve yvar= S112.2m
EndDiagram
#
Diagram 21
head1= 'Track-shift forces, filtered in a sliding 2m-min-window, second bogie'
head2= 'Curve_1= S121.2m Max=$S121.2mMAX Mean=$S121.2mMED P99.85=$FDS121.2mP5'
head3= 'Curve_2= S122.2m Max=$S122.2mMAX Mean=$S122.2mMED P99.85=$FDS122.2mP5'
draw_limit= 75e3 # limit_expl= 'Max allowed track-shift force 75e3'
y_bot= -20e3 yint/cm= 20e3
Curve yvar= S121.2m
Curve yvar= S122.2m
EndDiagram
#
Diagram 12
head1= 'Flange climb ratio, filtered in a sliding 2m-mean-window, first wheelset'
head2= 'Curve_1= Y/Q111l.50 Max=$Y/Q111l.50MAX Mean=$Y/Q111l.50MED P99.85=$FDY/Q111l.50P5'
head3= 'Curve_2= Y/Q111r.50 Max=$Y/Q111r.50MAX Mean=$Y/Q111r.50MED P99.85=$FDY/Q111r.50P5'
draw_limit= 0.8 # limit_expl= 'Max allowed flange-climbing ratio = 0.8'
y_bot= -.2 yint/cm= .2
Curve yvar= Y/Q111l.50
Curve yvar= Y/Q111r.50
EndDiagram
#
Diagram 22
head1= 'Flange climb ratio, filtered in a sliding 2m-mean-window, second wheelset'
head2= 'Curve_1= Y/Q112l.50 Max=$Y/Q112l.50MAX Mean=$Y/Q112l.50MED P99.85=$FDY/Q112l.50P5'
head3= 'Curve_2= Y/Q112r.50 Max=$Y/Q112r.50MAX Mean=$Y/Q112r.50MED P99.85=$FDY/Q112r.50P5'
draw_limit= 0.8 # limit_expl= 'Max allowed flange-climbing ratio = 0.8'
y_bot= -.2 yint/cm= .2
Curve yvar= Y/Q112l.50
Curve yvar= Y/Q112r.50
EndDiagram
#
EndPage
##
## Plot pantograph sway
## -----------------------------------------------------------
Page
xvar= lsa_111.pn
x_left= Xplot_start
Diagram 11
head1= 'Pantograph sway'
head2= 'Curve_1= car1b1.y Min=$car1b1.yMIN Mean=$car1b1.yMED P0.15=$FDcar1b1.yP1'
head3= 'Curve_1= car1b2.y Min=$car1b2.yMIN Mean=$car1b2.yMED P0.15=$FDcar1b2.yP1'
Curve yvar= car1b1.y
Curve yvar= car1b2.y
EndDiagram
EndPage
##
## Vehicle overturning
## -----------------------------------------------------------
Page
xvar= lsa_111.pn
x_left= Xplot_start
Diagram 11
head1= 'Vehicle overturning'
head2= 'Curve_1= bt Max=$btMAX Mean=$btMED P99.85=$FDbtP5'
Curve yvar= bt
EndDiagram
EndPage
##
## Plot motion in lateral stops
## -----------------------------------------------------------
Stat kycbs11.d tstart=Xeval_start tstop=Xeval_stop Tname=lsb_11.pn Ptile= .15 25 50 75 99.85
Stat kycbs12.d tstart=Xeval_start tstop=Xeval_stop Tname=lsb_12.pn Ptile= .15 25 50 75 99.85
Page
xvar= lsa_111.pn
x_left= Xplot_start
Diagram 11
head1= 'Motion in lateral stops'
head2= 'Curve_1= kycbs11.d Max=$kycbs11.dMAX Mean=$kycbs11.dMED P99.85=$FDkycbs11.dP5'
head3= 'Curve_1= kycbs12.d Max=$kycbs12.dMAX Mean=$kycbs12.dMED P99.85=$FDkycbs12.dP5'
Curve yvar= kycbs11.d
Curve yvar= kycbs12.d
EndDiagram
EndPage
stop
##
## Plot the 2m-filtered track-shift forces
## -----------------------------------------------------------
Page
xvar= lsa_111.pn
x_left= Xplot_start
limit_expl = 'Max allowed track-shift force according to prudhomme'
draw_limit = 75e3
Diagram 11
Curve yvar= S111.2m
Diagram 12
Curve yvar= S112.2m
Diagram 21
Curve yvar= S121.2m
Diagram 22
Curve yvar= S122.2m
EndPage
##
## Plot the 50ms-filtered flange-climbing ratio
## -----------------------------------------------------------
Page
xvar= lsa_111.pn
x_left= Xplot_start
limit_expl = 'Max allowed flange-climbing ratio'
draw_limit = 0.8
Diagram 11
Curve yvar= Y/Q111l.50
Diagram 12
Curve yvar= Y/Q112l.50
Diagram 21
Curve yvar= Y/Q111r.50
Diagram 22
Curve yvar= Y/Q112r.50
EndPage