Updates from rel.9910 to rel.0003
Author: Ingemar Persson (DEsolver)
Store results directly in variable slask, in order to avoid writing
to file *.ca. This change saves approx. 20 % of simulations times
plus the time spent in program conv_ca.
Adjust step size in heun_u so the time will hit every tout exactly.
Read command INSERT directly in subroutine input, in order to
reduce CPU-time in the input data phase.
If idebug greater than 0 write information on names of the insert-files and
substructures as the input reading proceed.
New command CURRENT_DIRECTORY which makes it possible for the
queue system to find the files, if the execution takes place
on an another computer.
Use computed goto statement in subroutine de, func and coupl, in order to
reduce the CPU-time.
Read track irregularities stored by "func intpl_track_irr2".
Write the Jacobian to an external file named $ident.modjac.
A new air spring "coupl k_air3_exp" with a more stable calculation
of the smooth friction part has been introduced. The smooth friction
is similar to the behavior of a low-pass filter.
Special code written for the linearization process of coupl k_air3
and k_air3_exp in program modal and fresp.
A new substructure file wr_coupl_ne1 replacing t_irr_ne1.
The substructure file wr_coupl_ne1 has a better and more
accurate lookup table for calculating track forces.
New hard coded function wr_coupl_nl1, with linear creep-creep force
New hard coded function wr_coupl_npol1, calculating the creep forces
according to O.Polach, Winterthur.
A new substructure file wr_coupl_nr1 has also been introduced
which uses Kalkers program Rollen for calculating the
To the file calc.out, write vector var integer formatted,
if the value is an address.
If command Ideal_track is given instead of an input data file
intpl_track_irr and intpl_track_irr2 will create a track without
Write only "Writing timestep" every 30 time passing tout, instead of
every 10th time as earlier.
If no contact force between wheel and rail set a/b equal to 1
instead of 0 as earlier.
If tstep is greater than tout set tstep equal to tout.
In integrator heun_u, adjust tstep if it is too big in the input data file.
If tstep greater than 0.001 or set tstep=0.001.
Increase number of sub diagrams from 20(4x5) to 153(17x9).
New input data command PAPER_FORMAT.
Command PAPER_FORMAT can be given the following values: A0, A1, A2,
A3 and A4.
The orientation of the paper is controlled by command XAXIS_ALONG.
Increase number of input data lines in a page from 101 to 301.
The automatic grid starts relative to the location of the X-axis.
The size of the automatic grid has been changed, in rel.0003
the distance between the lines is smaller if the diagram is
If LYCM is manually given don't try to make it into a even
multiple of D_YAX_MARKINGS.
Introduced a new command limit_line with possibility to draw 4
limit lines in a diagram.
Change the order of plotting the wheel- and rails, in order to
make the diagram more clear.
The length of the coupling must be longer than 1e-5 m to be
regarded as a coupling with length.
Pulldown menu File+read_runf reads a complete new model into the memory
of program GPLOT.
Program Group KPF
Writing a table above the "Coefficient of contact angle difference"
similar to the table written over the effective conicity plot.
A new plot containing the "Contact angle difference parameter".
Corrected an error when printing a wheel/rail-geometry function with
Swedish associated names. The angle of contact for the left wheel was
Compensation of the measured wheel profiles due to the static load which
bends the axle of the wheelset.
Possibility to add roll rotation and lateral shift of the profile due
to static axle load on the designed axle.
Do not bother if not the X-axis is consequently increasing,
in measurements this type of errors can happen.
New input data command RAIL_FL_ORIGO, which defines the distance
between the origin of the rail and the flange contact point.
As default RAIL_FL_ORIGO is set to 32.5 mm which is the case for
normal gauge track.
A new program in GENSYS which calculates the surface of the wheel and
rail expressed in a coordinate system oriented in the surface of the
Program Group TRACK
In the high frequency region use always random phase angles, because
otherwise all phase angles tends to be zero in the high frequency
region, causing a Dirac pulse in the beginning of the track.
Program TRC_STRIX_TRACK is rewritten so it requires an input data
file, like other program in the GENSYS package.
A standard input data file will be copied into the users directory,
if program TRC_STRIX_TRACK is started without an input data file.
New program calculating the Q-number of a track written in
New program calculating the Q-number of a track written in
trac- or trax- format.
A menu of all programs will be shown if genrun is started without
any input data file.
If a program is started without input data file, the user will asked
if a standard input data file shall be copied and opened in the
The programs npick and predat are also reading options from file
.gen_conf or gen_conf.
Script kpfr2matlab.awk for translating a wheel-rail geometry file
*.kpfr into a *.m-file readable by program matlab.
In the ident argument for program MPLOT, the user can give the ident
text as usual or the name of the MPdat-file containing the ident.
A new part in the GENSYS documentation containing a
over substructures modeling the
wheel-rail contact mechanism.
A new part in the GENSYS documentation containing
input data examples
earlier located in directory $gensys/hints.
- Make the input reading more robust in all MISC-, KPF- and TRACK-programs.
Read all lines into a character string before decoding the line, in
order to avoid Fortran I/O-Errors.
Blank lines in the input data file are not considered.
Lines beginning with a #-sign are considered as comments.
In graphical programs set minimum line thickness to 3.
Output to postscript printer is now written in postscript level 2,
in order to make it possible to select A4 or A3 papers in the paper tray.
The GENSYS home page has opened at http://home.swipnet.se/gensys.