Home Reference Manuals Return Track Group Documentation

Users Manual for Program TRC_IPSD

   Introduction
   Input data commands
   Example

Introduction

Program TRC_IPSD creates track irregularities which fulfills a given PSD-spectra. In command RNDFIG the user can choose a random number which will be used when selecting the phase angles between different waves in the track irregularity spectra. As an alternative the user can in command INFIL feed an old existent track into program TRC_IPSD and the old track will work as a base when generating the new track which fulfills required PSD-spectra.

Input data commands

Input data is read in free format, valid separators between the input values are <space>, <comma>, <tab>, <equal sign> or <carriage return>. The commands can be written both in lower and upper case letters. The operation of the program is controlled by the commands described below; some of the commands also need arguments.

DRAW_ERRI_LEVELS: Draw high and low ERRI-levels in the PSD-diagrams.
The lines are drawn in a turquoise color. DRAW_ERRI_LEVELS can be set equal to "y" for yes or "n" for no.
Declared= Character*1 Default= "y"
GAUGE: Average gauge of the generated track.
Declared Real(8) Default= 1.435 [m]
INFIL: Existing track data file written in trac, trax or trax_wdesign format. Phase angles between different waves from an existing track will be used when creating the new track UTFIL. In case you do not have measured track irregularities, phase angles can be chosen randomly if you set INFIL equal to "no".
Declared= Character*132 Default= "no"
PSD_LAT: Vector describing the PSD-spectra in lateral direction. The PSD-spectra is given as node coordinates along the spectra: Y1, X1, Y2, X2, Y3,,, etc. The Y values are given in the unit [m*m/1/m]. The X axle is spatial frequency and the values are given in the unit [1/m]. The whole range from infinite long waves down to the shortest possible waves must be given in the spectra, i.e. X1 must be equal to 0. and Xn must be less or equal to 2/DX
Declared= Real(8) (100) Default= 0.
PSD_VERT: Vectors describing the PSD-spectra in vertical direction. The PSD-spectra are given as node coordinates along the spectra: Y1, X1, Y2, X2, Y3,,, etc. The Y values are given in the unit [m*m/1/m]. The X axle is spatial frequency and the values are given in the unit [1/m]. The whole range from infinite long waves down to 1[m] long waves must be given in spectra above, i.e. X1 must be equal to 0. and Xn must be equal to 1.
Declared= Real(8) (100) Default= 0.
PSD_GAUGE: Vectors describing the PSD-spectra in gauge direction. The PSD-spectra are given as node coordinates along the spectra: Y1, X1, Y2, X2, Y3,,, etc. The Y values are given in the unit [m*m/1/m]. The X axle is spatial frequency and the values are given in the unit [1/m]. The whole range from infinite long waves down to 1[m] long waves must be given in spectra above, i.e. X1 must be equal to 0. and Xn must be equal to 1.
Declared= Real(8) (100) Default= 0.
PSD_FI: Vectors describing the PSD-spectra in cant direction. The PSD-spectra are given as node coordinates along the spectra: Y1, X1, Y2, X2, Y3,,, etc. The Y values are given in the unit [rad*rad/1/m]. The X axle is spatial frequency and the values are given in the unit [1/m]. The whole range from infinite long waves down to 1[m] long waves must be given in spectra above, i.e. X1 must be equal to 0. and Xn must be equal to 1.
Declared= Real(8) (100) Default= 0.
RNDFIG: Initial value when calculating the random series numbers.
The random numbers are used as phase angles between different waves in the track irregularity spectra. RNDFIG must be greater than 0.001 and less than 0.999. Please avoid setting RNDFIG equal to 0.5.
Declared Real(8) Default= 0.1
UTFIL: File containing the output results.
The file contains a new track fulfilling the spectra given in PSD_LAT, PSD_VERT, PSD_GAUGE and PSD_FI. If the file extension of UTFIL equals "trax_wdesign", the output file will be written in trax_wdesign format. Otherwise the output will be written in trax-format.
Declared= Character*132 Default= 'trc_ipsd.trax'
XSTART: Start coordinate when reading track irregularity data from INFIL.
The track read from INFIL should start at least 250[m] before XSTART to avoid disturbances from the end. If INFILE starts too close to XSTART a warning message will be written.
Declared= Real(8)
Default= -1.E99 (Automatically set XSTART 250[m] after beginning of INFIL)
XSTOP: End coordinate when reading track irregularity data from INFIL.
The track read from INFIL should end at least 250[m] after XSTOP to avoid disturbances from the end. If INFILE starts too close to XSTOP a warning message will be written.
Declared= Real(8)
Default= 1.E99 (Automatically set XSTOP 250[m] before end of INFIL)
DX: Equidistant step between the points in the X-axis.
Declared= Real(8) Default= 0.5 [m]

Example #1

Create a PSD spectra with constant amplitude for all frequencies


##
##      Input data for program TRC_IPSD
##

 RNDFIG= 0.1
 INFIL = $gentrc/K0_4a21-n-k.trax  # Read phase angles from an existed trac-file.
 UTFIL = trc_ipsd.trax

 DX    = 0.5
 XSTART=-1.e99  XSTOP= 1.e99  GAUGE= 1.435

#
#        White PSD spectra
#
  PSD_LAT  = 1.E-6,     0,    1.E-6,     1.      # 1[mm] in lateral direction
  PSD_VERT = 1.E-6,     0,    1.E-6,     1.      # 1[mm] in vertical direction
  PSD_GAUGE= 1.E-6,     0,    1.E-6,     1.      # 1[mm] in gauge
  PSD_FI   = .444444e-6 0,    .444444e-6 1.      # 1[mm] in cant

Example #2
Create a PSD spectra according to ERRI QUESTION B176

In the ERRI report B176 the following PSD-spectras was presented:


  Sa=               Aa*0.8246^2              
      (omega^2+0.0206^2) * (omega^2+0.8246^2)

  Sv=               Av*0.8246^2              
      (omega^2+0.0206^2) * (omega^2+0.8246^2)

  Sc=               Av/bo^2*0.8246^2*omega^2                      
      (omega^2+0.0206^2) * (omega^2+0.8246^2) * (omega^2+0.4380^2)

Where:

Sa	=	PSD-spectra in lateral direction (alignment)
Sv	=	PSD-spectra in vertical direction
Sc	=	PSD-spectra of cant (expressed in radians)

The coefficients Aa and Av are as follows:

Low level of irregularity:

Aa	=	2.119e-7 [ m rad ]
Av	=	4.032e-7 [ m rad ]

High level of irregularity:

Aa	=	6.125e-7 [ m rad ]
Av	=	1.080e-6 [ m rad ]


##
##      Input data for program TRC_IPSD
##

 RNDFIG= 0.1
 INFIL = $gentrc/K0_4a21-n-k.trax  # Read phase angles from an existed trac-file.
 UTFIL = trc_ipsd.trax

 DX    = 0.5
 XSTART=-1.e99  XSTOP= 1.e99  GAUGE= 1.435

#
#        ERRI B176 level high
#
# Sa= 2*pi*6.125e-7*0.8246**2 / ((2*pi*freq)**2+0.0206**2) / ((2*pi*freq)**2+0.8246**2)   Lateral
# Sv= 2*pi*1.080e-6*0.8246**2 / ((2*pi*freq)**2+0.0206**2) / ((2*pi*freq)**2+0.8246**2)   Vertical
# Sc= 2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*freq)**2 / ((2*pi*freq)**2+0.0206**2) / ((2*pi*freq)**2+0.8246**2) / ((2*pi*freq)**2+0.4380**2)
#
  PSD_LAT  =
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.      )**2+0.0206**2) / ((2*pi*  0.      )**2+0.8246**2)`   0.             # Infinite wave length;  Unit [m^2/(1/m)]
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.003183)**2+0.0206**2) / ((2*pi*  0.003183)**2+0.8246**2)`   0.003183       # 314.16m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.010000)**2+0.0206**2) / ((2*pi*  0.010000)**2+0.8246**2)`   0.010000       # 100.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2)`   0.020000       #  50.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.050000)**2+0.0206**2) / ((2*pi*  0.050000)**2+0.8246**2)`   0.050000       #  20.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.100000)**2+0.0206**2) / ((2*pi*  0.100000)**2+0.8246**2)`   0.100000       #  10.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.200000)**2+0.0206**2) / ((2*pi*  0.200000)**2+0.8246**2)`   0.200000       #   5.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.500000)**2+0.0206**2) / ((2*pi*  0.500000)**2+0.8246**2)`   0.500000       #   2.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  1.000000)**2+0.0206**2) / ((2*pi*  1.000000)**2+0.8246**2)`   1.000000       #   1.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  2.000000)**2+0.0206**2) / ((2*pi*  2.000000)**2+0.8246**2)`   2.000000       #   0.50m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi* 10.000000)**2+0.0206**2) / ((2*pi* 10.000000)**2+0.8246**2)`  10.000000       #   0.10m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*100.000000)**2+0.0206**2) / ((2*pi*100.000000)**2+0.8246**2)` 100.000000       #   0.01m waves
#
  PSD_VERT =
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.      )**2+0.0206**2) / ((2*pi*  0.      )**2+0.8246**2)`   0.             # Infinite wave length;  Unit [m^2/(1/m)]
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.003183)**2+0.0206**2) / ((2*pi*  0.003183)**2+0.8246**2)`   0.003183       # 314.16m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.010000)**2+0.0206**2) / ((2*pi*  0.010000)**2+0.8246**2)`   0.010000       # 100.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2)`   0.020000       #  50.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.050000)**2+0.0206**2) / ((2*pi*  0.050000)**2+0.8246**2)`   0.050000       #  20.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.100000)**2+0.0206**2) / ((2*pi*  0.100000)**2+0.8246**2)`   0.100000       #  10.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.200000)**2+0.0206**2) / ((2*pi*  0.200000)**2+0.8246**2)`   0.200000       #   5.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.500000)**2+0.0206**2) / ((2*pi*  0.500000)**2+0.8246**2)`   0.500000       #   2.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  1.000000)**2+0.0206**2) / ((2*pi*  1.000000)**2+0.8246**2)`   1.000000       #   1.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  2.000000)**2+0.0206**2) / ((2*pi*  2.000000)**2+0.8246**2)`   2.000000       #   0.50m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi* 10.000000)**2+0.0206**2) / ((2*pi* 10.000000)**2+0.8246**2)`  10.000000       #   0.10m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*100.000000)**2+0.0206**2) / ((2*pi*100.000000)**2+0.8246**2)` 100.000000       #   0.01m waves
#
  PSD_GAUGE=
     `0.05*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2)`   0.             # Infinite wave length;  Unit [m^2/(1/m)]
     `0.05*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2)`   0.003183       # 314.16m waves
     `0.05*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2)`   0.010000       # 100.00m waves
     `0.05*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2)`   0.020000       #  50.00m waves
     `0.25*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.050000)**2+0.0206**2) / ((2*pi*  0.050000)**2+0.8246**2)`   0.050000       #  20.00m waves
     `0.50*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.100000)**2+0.0206**2) / ((2*pi*  0.100000)**2+0.8246**2)`   0.100000       #  10.00m waves
     `1.00*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.200000)**2+0.0206**2) / ((2*pi*  0.200000)**2+0.8246**2)`   0.200000       #   5.00m waves
     `1.00*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.500000)**2+0.0206**2) / ((2*pi*  0.500000)**2+0.8246**2)`   0.500000       #   2.00m waves
     `1.00*2*pi*6.125e-7*0.8246**2 / ((2*pi*  1.000000)**2+0.0206**2) / ((2*pi*  1.000000)**2+0.8246**2)`   1.000000       #   1.00m waves
     `1.00*2*pi*6.125e-7*0.8246**2 / ((2*pi*  2.000000)**2+0.0206**2) / ((2*pi*  2.000000)**2+0.8246**2)`   2.000000       #   0.50m waves
     `1.00*2*pi*6.125e-7*0.8246**2 / ((2*pi* 10.000000)**2+0.0206**2) / ((2*pi* 10.000000)**2+0.8246**2)`  10.000000       #   0.10m waves
     `1.00*2*pi*6.125e-7*0.8246**2 / ((2*pi*100.000000)**2+0.0206**2) / ((2*pi*100.000000)**2+0.8246**2)` 100.000000       #   0.01m waves
#
  PSD_FI   =
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.      )**2 / ((2*pi*  0.      )**2+0.0206**2) / ((2*pi*  0.      )**2+0.8246**2) / ((2*pi*  0.      )**2+0.4380**2)`   0.             # Infinite wave length;  Unit [rad^2/(1/m)]
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.010000)**2 / ((2*pi*  0.010000)**2+0.0206**2) / ((2*pi*  0.010000)**2+0.8246**2) / ((2*pi*  0.010000)**2+0.4380**2)`   0.010000       # 100.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.020000)**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2) / ((2*pi*  0.020000)**2+0.4380**2)`   0.020000       #  50.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.050000)**2 / ((2*pi*  0.050000)**2+0.0206**2) / ((2*pi*  0.050000)**2+0.8246**2) / ((2*pi*  0.050000)**2+0.4380**2)`   0.050000       #  20.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.100000)**2 / ((2*pi*  0.100000)**2+0.0206**2) / ((2*pi*  0.100000)**2+0.8246**2) / ((2*pi*  0.100000)**2+0.4380**2)`   0.100000       #  10.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.200000)**2 / ((2*pi*  0.200000)**2+0.0206**2) / ((2*pi*  0.200000)**2+0.8246**2) / ((2*pi*  0.200000)**2+0.4380**2)`   0.200000       #   5.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.500000)**2 / ((2*pi*  0.500000)**2+0.0206**2) / ((2*pi*  0.500000)**2+0.8246**2) / ((2*pi*  0.500000)**2+0.4380**2)`   0.500000       #   2.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  1.000000)**2 / ((2*pi*  1.000000)**2+0.0206**2) / ((2*pi*  1.000000)**2+0.8246**2) / ((2*pi*  1.000000)**2+0.4380**2)`   1.000000       #   1.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  2.000000)**2 / ((2*pi*  2.000000)**2+0.0206**2) / ((2*pi*  2.000000)**2+0.8246**2) / ((2*pi*  2.000000)**2+0.4380**2)`   2.000000       #   0.50m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi* 10.000000)**2 / ((2*pi* 10.000000)**2+0.0206**2) / ((2*pi* 10.000000)**2+0.8246**2) / ((2*pi* 10.000000)**2+0.4380**2)`  10.000000       #   0.10m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*100.000000)**2 / ((2*pi*100.000000)**2+0.0206**2) / ((2*pi*100.000000)**2+0.8246**2) / ((2*pi*100.000000)**2+0.4380**2)` 100.000000       #   0.01m waves

Users Manual for Program TRC_IPSD

Table of Contents

Introduction

Input data commands

Example #1

Example #2 Create a PSD spectra according to ERRI QUESTION B176

Example #2
Create a PSD spectra according to ERRI QUESTION B176