Home Search Reference Manuals Return
## ## Input data for program MPLOT ## ## ----------------------------------------- ## 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