Home    Search    Reference Manual    Return   



Updates from rel.9910 to rel.0003


Author: Ingemar Persson (DEsolver)




Program CALC
  1. 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.
  2. Adjust step size in heun_u so the time will hit every tout exactly.
  3. 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.
  4. Use computed goto statement in subroutine de, func and coupl, in order to reduce the CPU-time.
  5. Read track irregularities stored by "func intpl_track_irr2".
  6. Write the Jacobian to an external file named $ident.modjac.
  7. 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.
  8. Special code written for the linearization process of coupl k_air3 and k_air3_exp in program modal and fresp.
  9. 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.
  10. New hard coded function wr_coupl_nl1, with linear creep-creep force relationship.
  11. New hard coded function wr_coupl_npol1, calculating the creep forces according to O.Polach, Winterthur.
  12. A new substructure file wr_coupl_nr1 has also been introduced which uses Kalkers program Rollen for calculating the creep forces.
  13. To the file calc.out, write vector var integer formatted, if the value is an address.
  14. If command Ideal_track is given instead of an input data file intpl_track_irr and intpl_track_irr2 will create a track without track irregularities.
  15. Write only "Writing timestep" every 30 time passing tout, instead of every 10th time as earlier.
  16. If no contact force between wheel and rail set a/b equal to 1 instead of 0 as earlier.
  17. If tstep is greater than tout set tstep equal to tout.
  18. 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.



Program MPLOT
  1. Increase number of sub diagrams from 20(4x5) to 153(17x9).
  2. 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.
  3. Increase number of input data lines in a page from 101 to 301.
  4. The automatic grid starts relative to the location of the X-axis.
  5. The size of the automatic grid has been changed, in rel.0003 the distance between the lines is smaller if the diagram is smaller.
  6. If LYCM is manually given don't try to make it into a even multiple of D_YAX_MARKINGS.
  7. Introduced a new command limit_line with possibility to draw 4 limit lines in a diagram.



Program GPLOT
  1. Change the order of plotting the wheel- and rails, in order to make the diagram more clear.
  2. The length of the coupling must be longer than 1e-5 m to be regarded as a coupling with length.
  3. Pulldown menu File+read_runf reads a complete new model into the memory of program GPLOT.



Program Group KPF

Program KPF
  1. Writing a table above the "Coefficient of contact angle difference" similar to the table written over the effective conicity plot.
  2. A new plot containing the "Contact angle difference parameter".
  3. Corrected an error when printing a wheel/rail-geometry function with Swedish associated names. The angle of contact for the left wheel was missing.
  4. Compensation of the measured wheel profiles due to the static load which bends the axle of the wheelset.
  5. Possibility to add roll rotation and lateral shift of the profile due to static axle load on the designed axle.


Program KPF_ROT
  1. Do not bother if not the X-axis is consequently increasing, in measurements this type of errors can happen.
  2. 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.


Program PRE_CONTACT
  1. 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 contact point.



Program Group TRACK

Program TRC_IPSD
  1. 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
  1. 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.


Program TRC_QNUMB_MAUZ
  1. New program calculating the Q-number of a track written in mauz-format.


Program TRC_QNUMB_TRACK
  1. New program calculating the Q-number of a track written in trac- or trax- format.



Scripts
  1. A menu of all programs will be shown if genrun is started without any input data file.
  2. 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 users directory.
  3. The programs npick and predat are also reading options from file .gen_conf or gen_conf.
  4. Script kpfr2matlab.awk for translating a wheel-rail geometry file *.kpfr into a *.m-file readable by program matlab.
  5. 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.



Documentation
  1. A new part in the GENSYS documentation containing a table over substructures modeling the wheel-rail contact mechanism.
  2. A new part in the GENSYS documentation containing input data examples earlier located in directory $gensys/hints.



Misc
  1. 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.
  2. In graphical programs set minimum line thickness to 3.
  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.
  4. The GENSYS home page has opened at http://home.swipnet.se/gensys.