Create and animate a wheel/rail-geometry function

gp	Animation files created by pre-processor KPF
kpff	Input data for program KPF
kpfr	Wheel/rail-geometry functions created by pre-processor KPF
r_prof	A directory containing rail profiles
runf	Model of a wheelset

Run the example kpff/uic60_i40_i20.kpff

Double-click file kpff/uic60_i40_i20.kpff and select open->run with the right mouse button. The execution of program KPF will now start, using file kpff/uic60_i40_i20.kpff as input.

As output three result files will be produced:

kpfr/uic60_i40_i20.kpfr	Wheel/rail-geometry function to be used as input in the main calculation programs in GENSYS
diags/uic60_i40_i20.pdf	A PDF file showing the wheel/rail-geometry functions
gp/uic60_i40_i20.gp	GPdat-file for animation of the wheel/rail-geometry functions

---

Examine the results in kpfr

File kpfr/uic60_i40_i20.kpfr is a plain text-file. Just mark the file by double-clicking the file and with the right mouse-button select open->op

---

View PDF results

File diags/uic60_i40_i20.pdf is a PDF file. View its contents by double-clicking the file and with the right mouse-button select open->run

For an explanation of the diagrams, see this link kpf.html#jOutput

---

Examine the results in postprocessor GPLOT

• Double-click directory runf/ in the subdirectory text area.
• Double-click file runf/kpfr_animation.tsimf in the file text area.
• Press the Gplot-pushbutton.
• Open file gp/uic60_i40_i20.gp by selecting "File->read_GPdat" in the pulldown menu bar
• Select "View->show 4 diags" from the pulldown menu bar
• The upper-left and upper-right drawing areas now shows a figure of the wheel/rail-geometries
• Please fill in wheel number 111l for the lower left drawing area and 111r for the lower right drawing area
• Set the radio buttons Drfn for the lower left and right drawing areas.
• Start the animation by pressing the "Fwd"- or "Back"- pushbutton.
• The animation can also be controlled by pulling the right scale in the bottom of the GPLOT-window.

Convert miniprof ban-files to rail-files

• Go to directory "r_prof/miniprof"
• Create an input data file for program kpf_conv_miniprof2kpf_rail:
  • In program genfile, click Start -> Program -> Kpf -> kpf_conv_miniprof2kpf_rail in the pull-down menu.
    Type "y" on the question "Do you want to create a Master input data file"
    
  • In a terminal write:

    kpf_conv_miniprof2kpf_rail -M
    

• Open the newly created file "Master.kpf_conv_miniprof2kpf_railf" for editing.
• Launch program kpf_conv_miniprof2kpf_rail, using the input data file "Master.kpf_conv_miniprof2kpf_railf"
• After execution all miniprof-files defined in command infil has been translated to rail-files.

After translating the rail profiles. Please, open and inspect the profiles with script gen_plot_rail.sh, to check that the profiles are correctly orientated. Eventually, it can be necessary to change the value of input data parameter TURN.
Files with extension .rail are associated with program gen_plot_rail.sh. To see the profile, mark the .rail-file and select genrun.

---

Convert miniprof whl-files to wheel-files

• Go to directory "w_prof/miniprof"
• Create an input data file for program kpf_conv_miniprof2kpf_wheel:
  • In program genfile, click Start -> Program -> Kpf -> kpf_conv_miniprof2kpf_wheel in the pull-down menu.
    Type "y" on the question "Do you want to create a Master input data file"
    
  • In a terminal write:

    kpf_conv_miniprof2kpf_wheel -M
    

• Open the newly created file "Master.kpf_conv_miniprof2kpf_wheelf" for editing.
• Launch program kpf_conv_miniprof2kpf_wheel, using the input data file "Master.kpf_conv_miniprof2kpf_wheelf"
• After execution all miniprof-files defined in command infil has been translated to wheel-files.

After translating the wheel profiles. Please, open and inspect the profiles with script gen_plot_wheel.sh. Files with extension .wheel are associated with program gen_plot_wheel.sh. To see the profile, mark the .rail-file and select genrun.

---

Launch multiple kpf-calculation by using program OPTI

• Open file "optif/manual_step.optif" for editing.
• Launch progam OPTI from the pulldown-menu "Gensys->Run_Current_File"
• Launch queue monitor sview and wait until all KPF-calculations are finished.

---

Make tables of the results

• Program KPF creates a directory named "m". In this directory program KPF writes results formatted as matlab/octave-files.
• With the UNIX-command grep and sort the user can create summary tables of the results e.g.:
  
  Information about the wheel profiles:
  Flange thickness measured at an height of 10mm from nominal running circle:

   grep '^ *Sd[rl]= '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r 

  Flange thickness measured at an height of 3mm from nominal running circle:

   grep '^ *Sd3[rl]= '  m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r 

  Flange height:

   grep '^ *Sh[rl]= '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r 

  Max flange angle:

   grep '^ *beta[rl]= ' m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r 

  Flange flank qR:

   grep '^ *qR[rl]= '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r 

  Flange contact position:

   grep '^ *S[rl] = '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 4 -r 

  Flange gauge:

   grep '^ *SR = '      m/*.m | column -t --table-columns-limit 2 | sort -g -k 4 -r 

  False flange:

   grep '^ *Th[rl]= '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r 

  Gradient Index Profile wheel left:

   grep '^ *GIPw_L= '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r 

  Gradient Index Profile wheel right:

   grep '^ *GIPw_R= '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r 

  Information about the rail profiles:
  Gauge measuring position GAUGE_MEAS_INTERVAL under Top_of_Rail:

   grep '^ *G[rl] = '     m/*.m | column -t --table-columns-limit 2 | sort -g -k 4 -r 

  Track gauge measured at GAUGE_MEAS_INTERVAL under Top_of_Rail:

   grep '^ *Ga = '        m/*.m | column -t --table-columns-limit 2 | sort -g -k 4 -r 

  Gauge measuring position 3mm under Top_of_Rail:

   grep '^ *G[rl]3= '     m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r 

  Gradient Index Profile rail left:

   grep -A 10 '^ *GIPr_L= ' m/*.m | grep '\(1435\)' | column -t --table-columns-limit 2 | sort -g -k 2 -r 

  Gradient Index Profile rail right:

   grep -A 10 '^ *GIPr_R= ' m/*.m | grep '\(1435\)' | column -t --table-columns-limit 2 | sort -g -k 2 -r 

  Information about the wheel/rail-combination:
  Minimum curve left without flange contact:

   grep -A 10 '^ *RSI_Re_L= '   m/*.m | grep '\(1435\)' | column -t --table-columns-limit 2 | sort -g -k 2 

  Minimum curve right without flange contact:

   grep -A 10 '^ *RSI_Re_R= '   m/*.m | grep '\(1435\)' | column -t --table-columns-limit 2 | sort -g -k 2 -r 

  Max contact stress:

   grep '^ *SigmaMax_[rl]= '    m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r 

  Wheelset lateral offset when both left and right wheel have the same effective rolling radius:

   grep -A 10 '^ *Conicity_Yoff     = ' m/*.m | grep '\(1435\)' | column -t --table-columns-limit 2 | sort -g -k 2 -r 

  Conicity according to UIC519 at gauge 1435mm:

   grep -A 10 '^ *Conicity_UIC519   = ' m/*.m | grep '\(1435\)' | awk '{print $1 " " $5}' | column -t --table-columns-limit 2 | sort -g -k 2 -r 

---

Results created by the above grep-commands

Flange thickness measured at an height of 10mm from nominal running circle:

grep '^ *Sd[rl]= '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r
m/uic60_i40_i20.m:  Sdr=   32.4997;  % Flange thickness at 10mm, right wheel
m/uic60_i40_i20.m:  Sdl=   32.4997;  % Flange thickness at 10mm, left wheel
m/hig3_S1002.m:     Sdr=   32.4997;  % Flange thickness at 10mm, right wheel
m/hig3_S1002.m:     Sdl=   32.4997;  % Flange thickness at 10mm, left wheel
m/hig2_S1002.m:     Sdr=   32.4997;  % Flange thickness at 10mm, right wheel
m/hig2_S1002.m:     Sdl=   32.4997;  % Flange thickness at 10mm, left wheel
m/hig1_S1002.m:     Sdr=   32.4997;  % Flange thickness at 10mm, right wheel
m/hig1_S1002.m:     Sdl=   32.4997;  % Flange thickness at 10mm, left wheel
m/hig3_0611.m:      Sdr=   32.0844;  % Flange thickness at 10mm, right wheel
m/hig3_0611.m:      Sdl=   32.0844;  % Flange thickness at 10mm, left wheel
m/hig2_0611.m:      Sdr=   32.0844;  % Flange thickness at 10mm, right wheel
m/hig2_0611.m:      Sdl=   32.0844;  % Flange thickness at 10mm, left wheel
m/hig1_0611.m:      Sdr=   32.0844;  % Flange thickness at 10mm, right wheel
m/hig1_0611.m:      Sdl=   32.0844;  % Flange thickness at 10mm, left wheel
m/hig3_0011.m:      Sdr=   31.4205;  % Flange thickness at 10mm, right wheel
m/hig3_0011.m:      Sdl=   31.4205;  % Flange thickness at 10mm, left wheel
m/hig2_0011.m:      Sdr=   31.4205;  % Flange thickness at 10mm, right wheel
m/hig2_0011.m:      Sdl=   31.4205;  % Flange thickness at 10mm, left wheel
m/hig1_0011.m:      Sdr=   31.4205;  % Flange thickness at 10mm, right wheel
m/hig1_0011.m:      Sdl=   31.4205;  % Flange thickness at 10mm, left wheel

Flange thickness measured at an height of 3mm from nominal running circle:

grep '^ *Sd3[rl]= '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r
m/hig3_0611.m:      Sd3r=   43.3294;  % Flange thickness at 3mm, right wheel
m/hig3_0611.m:      Sd3l=   43.3294;  % Flange thickness at 3mm, left wheel
m/hig2_0611.m:      Sd3r=   43.3294;  % Flange thickness at 3mm, right wheel
m/hig2_0611.m:      Sd3l=   43.3294;  % Flange thickness at 3mm, left wheel
m/hig1_0611.m:      Sd3r=   43.3294;  % Flange thickness at 3mm, right wheel
m/hig1_0611.m:      Sd3l=   43.3294;  % Flange thickness at 3mm, left wheel
m/uic60_i40_i20.m:  Sd3r=   42.9773;  % Flange thickness at 3mm, right wheel
m/uic60_i40_i20.m:  Sd3l=   42.9773;  % Flange thickness at 3mm, left wheel
m/hig3_S1002.m:     Sd3r=   42.9773;  % Flange thickness at 3mm, right wheel
m/hig3_S1002.m:     Sd3l=   42.9773;  % Flange thickness at 3mm, left wheel
m/hig2_S1002.m:     Sd3r=   42.9773;  % Flange thickness at 3mm, right wheel
m/hig2_S1002.m:     Sd3l=   42.9773;  % Flange thickness at 3mm, left wheel
m/hig1_S1002.m:     Sd3r=   42.9773;  % Flange thickness at 3mm, right wheel
m/hig1_S1002.m:     Sd3l=   42.9773;  % Flange thickness at 3mm, left wheel
m/hig3_0011.m:      Sd3r=   42.5711;  % Flange thickness at 3mm, right wheel
m/hig3_0011.m:      Sd3l=   42.5711;  % Flange thickness at 3mm, left wheel
m/hig2_0011.m:      Sd3r=   42.5711;  % Flange thickness at 3mm, right wheel
m/hig2_0011.m:      Sd3l=   42.5711;  % Flange thickness at 3mm, left wheel
m/hig1_0011.m:      Sd3r=   42.5711;  % Flange thickness at 3mm, right wheel
m/hig1_0011.m:      Sd3l=   42.5711;  % Flange thickness at 3mm, left wheel

Flange height:

grep '^ *Sh[rl]= '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r
m/hig3_0611.m:      Shr=   29.0499;  % Flange height rel to nom running circle, right wheel
m/hig3_0611.m:      Shl=   29.0499;  % Flange height rel to nom running circle, left wheel
m/hig2_0611.m:      Shr=   29.0499;  % Flange height rel to nom running circle, right wheel
m/hig2_0611.m:      Shl=   29.0499;  % Flange height rel to nom running circle, left wheel
m/hig1_0611.m:      Shr=   29.0499;  % Flange height rel to nom running circle, right wheel
m/hig1_0611.m:      Shl=   29.0499;  % Flange height rel to nom running circle, left wheel
m/hig3_0011.m:      Shr=   28.7425;  % Flange height rel to nom running circle, right wheel
m/hig3_0011.m:      Shl=   28.7425;  % Flange height rel to nom running circle, left wheel
m/hig2_0011.m:      Shr=   28.7425;  % Flange height rel to nom running circle, right wheel
m/hig2_0011.m:      Shl=   28.7425;  % Flange height rel to nom running circle, left wheel
m/hig1_0011.m:      Shr=   28.7425;  % Flange height rel to nom running circle, right wheel
m/hig1_0011.m:      Shl=   28.7425;  % Flange height rel to nom running circle, left wheel
m/uic60_i40_i20.m:  Shr=   28.0000;  % Flange height rel to nom running circle, right wheel
m/uic60_i40_i20.m:  Shl=   28.0000;  % Flange height rel to nom running circle, left wheel
m/hig3_S1002.m:     Shr=   28.0000;  % Flange height rel to nom running circle, right wheel
m/hig3_S1002.m:     Shl=   28.0000;  % Flange height rel to nom running circle, left wheel
m/hig2_S1002.m:     Shr=   28.0000;  % Flange height rel to nom running circle, right wheel
m/hig2_S1002.m:     Shl=   28.0000;  % Flange height rel to nom running circle, left wheel
m/hig1_S1002.m:     Shr=   28.0000;  % Flange height rel to nom running circle, right wheel
m/hig1_S1002.m:     Shl=   28.0000;  % Flange height rel to nom running circle, left wheel

Max flange angle:

grep '^ *beta[rl]= ' m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r
m/hig3_0611.m:      betar=   73.8875;  % Max flange angle, right wheel. Found at Y= 39.82
m/hig3_0611.m:      betal=   73.8875;  % Max flange angle, left  wheel. Found at Y= 39.82
m/hig2_0611.m:      betar=   73.8875;  % Max flange angle, right wheel. Found at Y= 39.82
m/hig2_0611.m:      betal=   73.8875;  % Max flange angle, left  wheel. Found at Y= 39.82
m/hig1_0611.m:      betar=   73.8875;  % Max flange angle, right wheel. Found at Y= 39.82
m/hig1_0611.m:      betal=   73.8875;  % Max flange angle, left  wheel. Found at Y= 39.82
m/hig3_0011.m:      betar=   73.7554;  % Max flange angle, right wheel. Found at Y= 40.42
m/hig3_0011.m:      betal=   73.7554;  % Max flange angle, left  wheel. Found at Y= 40.42
m/hig2_0011.m:      betar=   73.7554;  % Max flange angle, right wheel. Found at Y= 40.42
m/hig2_0011.m:      betal=   73.7554;  % Max flange angle, left  wheel. Found at Y= 40.42
m/hig1_0011.m:      betar=   73.7554;  % Max flange angle, right wheel. Found at Y= 40.42
m/hig1_0011.m:      betal=   73.7554;  % Max flange angle, left  wheel. Found at Y= 40.42
m/uic60_i40_i20.m:  betar=   69.9996;  % Max flange angle, right wheel. Found at Y= 39.32
m/uic60_i40_i20.m:  betal=   69.9996;  % Max flange angle, left  wheel. Found at Y= 39.32
m/hig3_S1002.m:     betar=   69.9996;  % Max flange angle, right wheel. Found at Y= 39.32
m/hig3_S1002.m:     betal=   69.9996;  % Max flange angle, left  wheel. Found at Y= 39.32
m/hig2_S1002.m:     betar=   69.9996;  % Max flange angle, right wheel. Found at Y= 39.32
m/hig2_S1002.m:     betal=   69.9996;  % Max flange angle, left  wheel. Found at Y= 39.32
m/hig1_S1002.m:     betar=   69.9996;  % Max flange angle, right wheel. Found at Y= 39.32
m/hig1_S1002.m:     betal=   69.9996;  % Max flange angle, left  wheel. Found at Y= 39.32

Flange flank qR:

grep '^ *qR[rl]= '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r
m/uic60_i40_i20.m:  qRr=   10.7930;  % Flange flank from flange_thick_height to a node
m/uic60_i40_i20.m:  qRl=   10.7930;  % 2mm under flange top, right and left wheel respectively
m/hig3_S1002.m:     qRr=   10.7930;  % Flange flank from flange_thick_height to a node
m/hig3_S1002.m:     qRl=   10.7930;  % 2mm under flange top, right and left wheel respectively
m/hig2_S1002.m:     qRr=   10.7930;  % Flange flank from flange_thick_height to a node
m/hig2_S1002.m:     qRl=   10.7930;  % 2mm under flange top, right and left wheel respectively
m/hig1_S1002.m:     qRr=   10.7930;  % Flange flank from flange_thick_height to a node
m/hig1_S1002.m:     qRl=   10.7930;  % 2mm under flange top, right and left wheel respectively
m/hig3_0011.m:      qRr=   10.4079;  % Flange flank from flange_thick_height to a node
m/hig3_0011.m:      qRl=   10.4079;  % 2mm under flange top, right and left wheel respectively
m/hig2_0011.m:      qRr=   10.4079;  % Flange flank from flange_thick_height to a node
m/hig2_0011.m:      qRl=   10.4079;  % 2mm under flange top, right and left wheel respectively
m/hig1_0011.m:      qRr=   10.4079;  % Flange flank from flange_thick_height to a node
m/hig1_0011.m:      qRl=   10.4079;  % 2mm under flange top, right and left wheel respectively
m/hig3_0611.m:      qRr=   10.2785;  % Flange flank from flange_thick_height to a node
m/hig3_0611.m:      qRl=   10.2785;  % 2mm under flange top, right and left wheel respectively
m/hig2_0611.m:      qRr=   10.2785;  % Flange flank from flange_thick_height to a node
m/hig2_0611.m:      qRl=   10.2785;  % 2mm under flange top, right and left wheel respectively
m/hig1_0611.m:      qRr=   10.2785;  % Flange flank from flange_thick_height to a node
m/hig1_0611.m:      qRl=   10.2785;  % 2mm under flange top, right and left wheel respectively

Flange contact position:

grep '^ *S[rl] = '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 4 -r
m/uic60_i40_i20.m:  Sr =  712.4996;  % Lateral distances to flange contact patches
m/uic60_i40_i20.m:  Sl =  712.4996;  % taken at height 10.00 above nominal running circle
m/hig3_S1002.m:     Sr =  712.4996;  % Lateral distances to flange contact patches
m/hig3_S1002.m:     Sl =  712.4996;  % taken at height 10.00 above nominal running circle
m/hig2_S1002.m:     Sr =  712.4996;  % Lateral distances to flange contact patches
m/hig2_S1002.m:     Sl =  712.4996;  % taken at height 10.00 above nominal running circle
m/hig1_S1002.m:     Sr =  712.4996;  % Lateral distances to flange contact patches
m/hig1_S1002.m:     Sl =  712.4996;  % taken at height 10.00 above nominal running circle
m/hig3_0611.m:      Sr =  712.0844;  % Lateral distances to flange contact patches
m/hig3_0611.m:      Sl =  712.0844;  % taken at height 10.00 above nominal running circle
m/hig2_0611.m:      Sr =  712.0844;  % Lateral distances to flange contact patches
m/hig2_0611.m:      Sl =  712.0844;  % taken at height 10.00 above nominal running circle
m/hig1_0611.m:      Sr =  712.0844;  % Lateral distances to flange contact patches
m/hig1_0611.m:      Sl =  712.0844;  % taken at height 10.00 above nominal running circle
m/hig3_0011.m:      Sr =  711.4205;  % Lateral distances to flange contact patches
m/hig3_0011.m:      Sl =  711.4205;  % taken at height 10.00 above nominal running circle
m/hig2_0011.m:      Sr =  711.4205;  % Lateral distances to flange contact patches
m/hig2_0011.m:      Sl =  711.4205;  % taken at height 10.00 above nominal running circle
m/hig1_0011.m:      Sr =  711.4205;  % Lateral distances to flange contact patches
m/hig1_0011.m:      Sl =  711.4205;  % taken at height 10.00 above nominal running circle

Flange gauge:

grep '^ *SR = '      m/*.m | column -t --table-columns-limit 2 | sort -g -k 4 -r
m/uic60_i40_i20.m:  SR = 1424.9993;  % Flange gauge. I.e. total distance Sr+Sl
m/hig3_S1002.m:     SR = 1424.9993;  % Flange gauge. I.e. total distance Sr+Sl
m/hig2_S1002.m:     SR = 1424.9993;  % Flange gauge. I.e. total distance Sr+Sl
m/hig1_S1002.m:     SR = 1424.9993;  % Flange gauge. I.e. total distance Sr+Sl
m/hig3_0611.m:      SR = 1424.1687;  % Flange gauge. I.e. total distance Sr+Sl
m/hig2_0611.m:      SR = 1424.1687;  % Flange gauge. I.e. total distance Sr+Sl
m/hig1_0611.m:      SR = 1424.1687;  % Flange gauge. I.e. total distance Sr+Sl
m/hig3_0011.m:      SR = 1422.8411;  % Flange gauge. I.e. total distance Sr+Sl
m/hig2_0011.m:      SR = 1422.8411;  % Flange gauge. I.e. total distance Sr+Sl
m/hig1_0011.m:      SR = 1422.8411;  % Flange gauge. I.e. total distance Sr+Sl

False flange:

grep '^ *Th[rl]= '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r
m/uic60_i40_i20.m:  Thr=    0.0000;  % Height of false flange, right wheel
m/uic60_i40_i20.m:  Thl=    0.0000;  % Height of false flange, left wheel
m/hig3_S1002.m:     Thr=    0.0000;  % Height of false flange, right wheel
m/hig3_S1002.m:     Thl=    0.0000;  % Height of false flange, left wheel
m/hig3_0611.m:      Thr=    0.0000;  % Height of false flange, right wheel
m/hig3_0611.m:      Thl=    0.0000;  % Height of false flange, left wheel
m/hig3_0011.m:      Thr=    0.0000;  % Height of false flange, right wheel
m/hig3_0011.m:      Thl=    0.0000;  % Height of false flange, left wheel
m/hig2_S1002.m:     Thr=    0.0000;  % Height of false flange, right wheel
m/hig2_S1002.m:     Thl=    0.0000;  % Height of false flange, left wheel
m/hig2_0611.m:      Thr=    0.0000;  % Height of false flange, right wheel
m/hig2_0611.m:      Thl=    0.0000;  % Height of false flange, left wheel
m/hig2_0011.m:      Thr=    0.0000;  % Height of false flange, right wheel
m/hig2_0011.m:      Thl=    0.0000;  % Height of false flange, left wheel
m/hig1_S1002.m:     Thr=    0.0000;  % Height of false flange, right wheel
m/hig1_S1002.m:     Thl=    0.0000;  % Height of false flange, left wheel
m/hig1_0611.m:      Thr=    0.0000;  % Height of false flange, right wheel
m/hig1_0611.m:      Thl=    0.0000;  % Height of false flange, left wheel
m/hig1_0011.m:      Thr=    0.0000;  % Height of false flange, right wheel
m/hig1_0011.m:      Thl=    0.0000;  % Height of false flange, left wheel

Gradient Index Profile wheel left:

grep '^ *GIPw_L= '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r
m/hig3_0611.m:      GIPw_L=    7.0243;  % GIPw, left  wheel
m/hig2_0611.m:      GIPw_L=    7.0243;  % GIPw, left  wheel
m/hig1_0611.m:      GIPw_L=    7.0243;  % GIPw, left  wheel
m/hig3_0011.m:      GIPw_L=    6.6089;  % GIPw, left  wheel
m/hig2_0011.m:      GIPw_L=    6.6089;  % GIPw, left  wheel
m/hig1_0011.m:      GIPw_L=    6.6089;  % GIPw, left  wheel
m/uic60_i40_i20.m:  GIPw_L=    6.4467;  % GIPw, left  wheel
m/hig3_S1002.m:     GIPw_L=    6.4467;  % GIPw, left  wheel
m/hig2_S1002.m:     GIPw_L=    6.4467;  % GIPw, left  wheel
m/hig1_S1002.m:     GIPw_L=    6.4467;  % GIPw, left  wheel

Gradient Index Profile wheel right:

grep '^ *GIPw_R= '   m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r
m/hig3_0611.m:      GIPw_R=    7.0243;  % GIPw, right wheel
m/hig2_0611.m:      GIPw_R=    7.0243;  % GIPw, right wheel
m/hig1_0611.m:      GIPw_R=    7.0243;  % GIPw, right wheel
m/hig3_0011.m:      GIPw_R=    6.6089;  % GIPw, right wheel
m/hig2_0011.m:      GIPw_R=    6.6089;  % GIPw, right wheel
m/hig1_0011.m:      GIPw_R=    6.6089;  % GIPw, right wheel
m/uic60_i40_i20.m:  GIPw_R=    6.4467;  % GIPw, right wheel
m/hig3_S1002.m:     GIPw_R=    6.4467;  % GIPw, right wheel
m/hig2_S1002.m:     GIPw_R=    6.4467;  % GIPw, right wheel
m/hig1_S1002.m:     GIPw_R=    6.4467;  % GIPw, right wheel

Gauge measuring position GAUGE_MEAS_INTERVAL under Top_of_Rail:

grep '^ *G[rl] = '     m/*.m | column -t --table-columns-limit 2 | sort -g -k 4 -r
m/hig3_S1002.m:     Gr =  717.5004;  % Gauge measuring positions, right and left rail
m/hig3_S1002.m:     Gl =  717.5004;  % taken at 14.00 under top of rail.
m/hig3_0611.m:      Gr =  717.5004;  % Gauge measuring positions, right and left rail
m/hig3_0611.m:      Gl =  717.5004;  % taken at 14.00 under top of rail.
m/hig3_0011.m:      Gr =  717.5004;  % Gauge measuring positions, right and left rail
m/hig3_0011.m:      Gl =  717.5004;  % taken at 14.00 under top of rail.
m/uic60_i40_i20.m:  Gl =  717.5002;  % taken at 14.00 under top of rail.
m/uic60_i40_i20.m:  Gr =  717.5001;  % Gauge measuring positions, right and left rail
m/hig2_S1002.m:     Gr =  717.5001;  % Gauge measuring positions, right and left rail
m/hig2_S1002.m:     Gl =  717.5001;  % taken at 14.00 under top of rail.
m/hig2_0611.m:      Gr =  717.5001;  % Gauge measuring positions, right and left rail
m/hig2_0611.m:      Gl =  717.5001;  % taken at 14.00 under top of rail.
m/hig2_0011.m:      Gr =  717.5001;  % Gauge measuring positions, right and left rail
m/hig2_0011.m:      Gl =  717.5001;  % taken at 14.00 under top of rail.
m/hig1_S1002.m:     Gr =  717.4984;  % Gauge measuring positions, right and left rail
m/hig1_S1002.m:     Gl =  717.4984;  % taken at 14.00 under top of rail.
m/hig1_0611.m:      Gr =  717.4984;  % Gauge measuring positions, right and left rail
m/hig1_0611.m:      Gl =  717.4984;  % taken at 14.00 under top of rail.
m/hig1_0011.m:      Gr =  717.4984;  % Gauge measuring positions, right and left rail
m/hig1_0011.m:      Gl =  717.4984;  % taken at 14.00 under top of rail.

Track gauge measured at GAUGE_MEAS_INTERVAL under Top_of_Rail:

grep '^ *Ga = '        m/*.m | column -t --table-columns-limit 2 | sort -g -k 4 -r
m/hig3_S1002.m:     Ga = 1435.0007;  % Track gauge. I.e. total distance Gr+Gl
m/hig3_0611.m:      Ga = 1435.0007;  % Track gauge. I.e. total distance Gr+Gl
m/hig3_0011.m:      Ga = 1435.0007;  % Track gauge. I.e. total distance Gr+Gl
m/uic60_i40_i20.m:  Ga = 1435.0004;  % Track gauge. I.e. total distance Gr+Gl
m/hig2_S1002.m:     Ga = 1435.0000;  % Track gauge. I.e. total distance Gr+Gl
m/hig2_0611.m:      Ga = 1435.0000;  % Track gauge. I.e. total distance Gr+Gl
m/hig2_0011.m:      Ga = 1435.0000;  % Track gauge. I.e. total distance Gr+Gl
m/hig1_S1002.m:     Ga = 1434.9968;  % Track gauge. I.e. total distance Gr+Gl
m/hig1_0611.m:      Ga = 1434.9968;  % Track gauge. I.e. total distance Gr+Gl
m/hig1_0011.m:      Ga = 1434.9968;  % Track gauge. I.e. total distance Gr+Gl

Gauge measuring position 3mm under Top_of_Rail:

grep '^ *G[rl]3= '     m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r
m/uic60_i40_i20.m:  Gr3=   12.6790;  % Gauge measuring position to the 3mm-point, right rail
m/hig1_S1002.m:     Gr3=   10.9410;  % Gauge measuring position to the 3mm-point, right rail
m/hig1_S1002.m:     Gl3=   10.9410;  % Gauge measuring position to the 3mm-point, left rail
m/hig1_0611.m:      Gr3=   10.9410;  % Gauge measuring position to the 3mm-point, right rail
m/hig1_0611.m:      Gl3=   10.9410;  % Gauge measuring position to the 3mm-point, left rail
m/hig1_0011.m:      Gr3=   10.9410;  % Gauge measuring position to the 3mm-point, right rail
m/hig1_0011.m:      Gl3=   10.9410;  % Gauge measuring position to the 3mm-point, left rail
m/uic60_i40_i20.m:  Gl3=    9.7481;  % Gauge measuring position to the 3mm-point, left rail
m/hig3_S1002.m:     Gr3=    9.7050;  % Gauge measuring position to the 3mm-point, right rail
m/hig3_S1002.m:     Gl3=    9.7050;  % Gauge measuring position to the 3mm-point, left rail
m/hig3_0611.m:      Gr3=    9.7050;  % Gauge measuring position to the 3mm-point, right rail
m/hig3_0611.m:      Gl3=    9.7050;  % Gauge measuring position to the 3mm-point, left rail
m/hig3_0011.m:      Gr3=    9.7050;  % Gauge measuring position to the 3mm-point, right rail
m/hig3_0011.m:      Gl3=    9.7050;  % Gauge measuring position to the 3mm-point, left rail
m/hig2_S1002.m:     Gr3=    9.4086;  % Gauge measuring position to the 3mm-point, right rail
m/hig2_S1002.m:     Gl3=    9.4086;  % Gauge measuring position to the 3mm-point, left rail
m/hig2_0611.m:      Gr3=    9.4086;  % Gauge measuring position to the 3mm-point, right rail
m/hig2_0611.m:      Gl3=    9.4086;  % Gauge measuring position to the 3mm-point, left rail
m/hig2_0011.m:      Gr3=    9.4086;  % Gauge measuring position to the 3mm-point, right rail
m/hig2_0011.m:      Gl3=    9.4086;  % Gauge measuring position to the 3mm-point, left rail

Gradient Index Profile rail left:

grep -A 10 '^ *GIPr_L= ' m/*.m | grep '\(1435\)' | column -t --table-columns-limit 2 | sort -g -k 2 -r
m/hig1_S1002.m-     7.87     % (1435)
m/hig1_0611.m-      7.87     % (1435)
m/hig1_0011.m-      7.87     % (1435)
m/hig3_S1002.m-     7.71     % (1435)
m/hig3_0611.m-      7.71     % (1435)
m/hig3_0011.m-      7.71     % (1435)
m/uic60_i40_i20.m-  7.63     % (1435)
m/hig2_S1002.m-     7.30     % (1435)
m/hig2_0611.m-      7.30     % (1435)
m/hig2_0011.m-      7.30     % (1435)

Gradient Index Profile rail right:

grep -A 10 '^ *GIPr_R= ' m/*.m | grep '\(1435\)' | column -t --table-columns-limit 2 | sort -g -k 2 -r
m/uic60_i40_i20.m-  9.76     % (1435)
m/hig1_S1002.m-     7.87     % (1435)
m/hig1_0611.m-      7.87     % (1435)
m/hig1_0011.m-      7.87     % (1435)
m/hig3_S1002.m-     7.71     % (1435)
m/hig3_0611.m-      7.71     % (1435)
m/hig3_0011.m-      7.71     % (1435)
m/hig2_S1002.m-     7.30     % (1435)
m/hig2_0611.m-      7.30     % (1435)
m/hig2_0011.m-      7.30     % (1435)

Minimum curve left without flange contact:

grep -A 10 '^ *RSI_Re_L= '   m/*.m | grep '\(1435\)' | column -t --table-columns-limit 2 | sort -g -k 2
m/uic60_i40_i20.m-  -7.041E+03                % (1435)
m/hig1_0611.m-      -7.576E+02                % (1435)
m/hig1_0011.m-      -7.077E+02                % (1435)
m/hig2_S1002.m-     -6.697E+02                % (1435)
m/hig1_S1002.m-     -6.329E+02                % (1435)
m/hig3_S1002.m-     -2.325E+02                % (1435)
m/hig2_0011.m-      -2.125E+02                % (1435)
m/hig3_0011.m-      -2.071E+02                % (1435)
m/hig3_0611.m-      -2.061E+02                % (1435)
m/hig2_0611.m-      -1.853E+02                % (1435)

Minimum curve right without flange contact:

grep -A 10 '^ *RSI_Re_R= '   m/*.m | grep '\(1435\)' | column -t --table-columns-limit 2 | sort -g -k 2 -r
m/hig1_0611.m-      7.576E+02                % (1435)
m/hig1_0011.m-      7.077E+02                % (1435)
m/hig2_S1002.m-     6.697E+02                % (1435)
m/hig1_S1002.m-     6.329E+02                % (1435)
m/hig3_S1002.m-     2.325E+02                % (1435)
m/hig2_0011.m-      2.125E+02                % (1435)
m/hig3_0011.m-      2.071E+02                % (1435)
m/hig3_0611.m-      2.061E+02                % (1435)
m/uic60_i40_i20.m-  2.019E+02                % (1435)
m/hig2_0611.m-      1.853E+02                % (1435)

Max contact stress:

grep '^ *SigmaMax_[rl]= '    m/*.m | column -t --table-columns-limit 2 | sort -g -k 3 -r
m/hig2_0611.m:      SigmaMax_r= 2.5606E+09;               % Max contact pressure, right wheel
m/hig2_0611.m:      SigmaMax_l= 2.5606E+09;               % Max contact pressure, left wheel
m/uic60_i40_i20.m:  SigmaMax_l= 2.3408E+09;               % Max contact pressure, left wheel
m/hig2_0011.m:      SigmaMax_r= 2.0774E+09;               % Max contact pressure, right wheel
m/hig2_0011.m:      SigmaMax_l= 2.0774E+09;               % Max contact pressure, left wheel
m/hig3_0611.m:      SigmaMax_r= 1.9591E+09;               % Max contact pressure, right wheel
m/hig3_0611.m:      SigmaMax_l= 1.9591E+09;               % Max contact pressure, left wheel
m/hig3_S1002.m:     SigmaMax_r= 1.8078E+09;               % Max contact pressure, right wheel
m/hig3_S1002.m:     SigmaMax_l= 1.8078E+09;               % Max contact pressure, left wheel
m/hig1_0611.m:      SigmaMax_r= 1.7267E+09;               % Max contact pressure, right wheel
m/hig1_0611.m:      SigmaMax_l= 1.7267E+09;               % Max contact pressure, left wheel
m/hig3_0011.m:      SigmaMax_r= 1.6831E+09;               % Max contact pressure, right wheel
m/hig3_0011.m:      SigmaMax_l= 1.6831E+09;               % Max contact pressure, left wheel
m/hig1_0011.m:      SigmaMax_r= 1.6791E+09;               % Max contact pressure, right wheel
m/hig1_0011.m:      SigmaMax_l= 1.6791E+09;               % Max contact pressure, left wheel
m/uic60_i40_i20.m:  SigmaMax_r= 1.6415E+09;               % Max contact pressure, right wheel
m/hig1_S1002.m:     SigmaMax_r= 1.2766E+09;               % Max contact pressure, right wheel
m/hig1_S1002.m:     SigmaMax_l= 1.2766E+09;               % Max contact pressure, left wheel
m/hig2_S1002.m:     SigmaMax_r= 1.2016E+09;               % Max contact pressure, right wheel
m/hig2_S1002.m:     SigmaMax_l= 1.2016E+09;               % Max contact pressure, left wheel

Wheelset lateral offset when both left and right wheel have the same effective rolling radius:

grep -A 10 '^ *Conicity_Yoff     = ' m/*.m | grep '\(1435\)' | column -t --table-columns-limit 2 | sort -g -k 2 -r
m/uic60_i40_i20.m-  2.7961         % (1435)
m/hig3_S1002.m-     0.0000         % (1435)
m/hig3_0611.m-      0.0000         % (1435)
m/hig3_0011.m-      0.0000         % (1435)
m/hig2_S1002.m-     0.0000         % (1435)
m/hig2_0611.m-      0.0000         % (1435)
m/hig2_0011.m-      0.0000         % (1435)
m/hig1_S1002.m-     0.0000         % (1435)
m/hig1_0611.m-      0.0000         % (1435)
m/hig1_0011.m-      0.0000         % (1435)

Conicity according to UIC519 at gauge 1435mm:

grep -A 10 '^ *Conicity_UIC519   = ' m/*.m | grep '\(1435\)' | awk '{print $1 " " $5}' | column -t --table-columns-limit 2 | sort -g -k 2 -r
m/hig2_S1002.m-     0.119
m/hig3_S1002.m-     0.093
m/hig1_S1002.m-     0.093
m/hig2_0611.m-      0.072
m/hig2_0011.m-      0.072
m/hig1_0611.m-      0.066
m/hig3_0611.m-      0.065
m/hig3_0011.m-      0.065
m/hig1_0011.m-      0.065
m/uic60_i40_i20.m-  0.034