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