Home
Reference Manuals
Return
#
#
# Substructure file t_irr_p1.ins rev.971113
#
# ----------------------------------------------------------
# Stiff axle.
#
# Track irregularities will be interpolated by cubic splines
# read from fields generated by command "func intpl_r",
# "func fifo_mem" or "fifo_mem2".
#
# Individual wheel/rail geometry functions for each wheel, will be
# read if the geometry functions contains the wheel number.
#
# All operations are made in one big function, in order to speed up
# the calculations.
#
# ----------------------------------------------------------
#
# This substructure do the following:
# 1) Set constraints to the axle in pitch motion, and in all
# directions for the track except in lateral direction.
# 2) Interpolate the track irregularities by cubic splines
# 3) Calculate cp_$2r.eta and cp_$2l.eta as the lateral
# displacement between wheel and rail
# 4) Interpolate the wheel-rail geometry functions linearly
# 5) Calculate the deformations and normal forces in the
# contact points.
# 6) Calculate creep and creep forces in the contact points.
#
# ----------------------------------------------------------
#
#
# Variables which must be defined before this substructure file can
# be used:
#
# Vo = Nominal speed of the vehicle in m/s
#
# Ro$2h = Nominal rolling radius, axle number $2, right side.
# ro$2r = If Ro$2h not can be found
# Ro$2v = Nominal rolling radius, axle number $2, left side.
# ro$2l = If Ro$2v not can be found
# Ro$1 = If not the variables above can be found
# ro$1 = If Ro$1 not can be found
# Ro = If ro$1 not can be found
# ro = If Ro not can be found
#
# Boh = Lateral distance between track center line and nominal
# rolling circle of right wheel
# Bo = If Boh not can be found
# bo = If Bo not can be found
# 0.75 = If bo not can be found
#
# Bov = Lateral distance between track center line and nominal
# rolling circle of left wheel
# -Bo = If Bov not can be found
# -bo = If Bo not can be found
# -.75 = If bo not can be found
#
# kph$2.my = Friction coefficient between wheel and rail, tread right side.
# cpt_$2r.mu= If kph$2.my not can be found
# my = If cpt_$2r.mu not can be found
# mu = If my not can be found
#
# kpv$2.my = Friction coefficient between wheel and rail, tread left side.
# cpt_$2l.mu= If kpv$2.my not can be found
# my = If cpt_$2l.mu not can be found
# mu = If my not can be found
#
# kfh$2.my = Friction coefficient between wheel and rail, flange left side.
# cpf_$2r.mu= If kfh$2.my not can be found
# myf = If cpf_$2r.mu not can be found
# muf = If myf not can be found
# my = If muf not can be found
# mu = If my not can be found
#
# kfv$2.my = Friction coefficient between wheel and rail, flange left side.
# cpf_$2l.mu= If kfv$2.my not can be found
# myf = If cpf_$2l.mu not can be found
# muf = If myf not can be found
# my = If muf not can be found
# mu = If my not can be found
#
#
# mulfact_nux_tread = Longitudinal creepage reduction factor on tread
# mulfact_nux = If mulfact_nux_tread not can be found
# mulfact_nuy_tread = Lateral creepage reduction factor on tread
# mulfact_nux = If mulfact_nuy_tread not can be found
# mulfact_spin_tread = Spin creepage reduction factor on tread
# mulfact_nux = If mulfact_spin_tread not can be found
# mulfact_nux_flange = Longitudinal creepage reduction factor on flange
# mulfact_nux = If mulfact_nux_flange not can be found
# mulfact_nuy_flange = Lateral creepage reduction factor on flange
# mulfact_nux = If mulfact_nuy_flange not can be found
# mulfact_spin_flange= Spin creepage reduction factor on flange
# mulfact_nux = If mulfact_spin_flange not can be found
# If one or more of the above multiplication factors are missing,
# they will be replaced with the value 1.
#
#
# Surface stiffness between rail and wheel
# ========================================
# knwr$2h_r= Stiffness perpendicular to the contact surface
# on tread right wheel, defined by "func intpl_r".
# Ex: func intpl_r knwr$2h_r -1.57 Horizontal stiffness
# 0. Vertical stiffness
# 1.57 Horizontal stiffness
# knwr$2r_r= If knwr$2h_r not can be found
# knwr$1_r = If knwr$2r_r not can be found
# knwr_r = If knwr$1_r not can be found
#
# knwr$2h_r.F0= Pre-stress force in spring knwr$2h_r.
# knwr$2r_r.F0= If knwr$2h_r.F0 not can be found
# knwr$1_r.F0 = If knwr$2r_r.F0 not can be found
# knwr_r.F0 = If knwr$1_r.F0 not can be found
#
#
# knwr$2v_r= Stiffness perpendicular to the contact surface
# on tread left wheel, defined by "func intpl_r".
# knwr$2l_r= If knwr$2v_r not can be found
# knwr$1_r = If knwr$2l_r not can be found
# knwr_r = If knwr$1_r not can be found
#
# knwr$2v_r.F0= Pre-stress force in spring knwr$2v_r.
# knwr$2l_r.F0= If knwr$2v_r.F0 not can be found
# knwr$1_r.F0 = If knwr$2l_r.F0 not can be found
# knwr_r.F0 = If knwr$1_r.F0 not can be found
#
# pknwr = The value to which the rail-wheel deformation shall be
# raised to when calculating the rail-wheel contact force.
# In Hertz theory this value should be set to 1.5.
# If not pknwr is defined pknwr will be set to 1. i.e.
# a linear stiffness between rail-wheels are assumed.
# Normally the rail-wheel stiffness is so high compared to
# other stiffnesses in the model that it only work as a
# method of calculating how much vertical load is resting
# on the tread and how much vertical load is resting
# on the flange.
#
#
# Surface stiffness between rail and wheel
# ========================================
# kyrt$2h= Lateral stiffness between right rail and track
# kyrt$2r= If kyrt$2h not can be found
# kyrt$1 = If kyrt$2r not can be found
# kyrt = If kyrt$1 not can be found
#
# kzrt$2h= Vertical stiffness between right rail and track
# kzrt$2r= If kzrt$2h not can be found
# kzrt$1 = If kzrt$2r not can be found
# kzrt = If kzrt$1 not can be found
#
# kzrt$2h.F0= Pre-stress force in spring kzrt$2h
# kzrt$2r.F0= If kzrt$2h.F0 not can be found
# kzrt$1.F0 = If kzrt$2r.F0 not can be found
# kzrt.F0 = If kzrt$1.F0 not can be found
# 0. = If kzrt.F0 not can be found
#
# cyrt$2h= Lateral damping between right rail and track
# cyrt$2r= If cyrt$2h not can be found
# cyrt$1 = If cyrt$2r not can be found
# cyrt = If cyrt$1 not can be found
#
# czrt$2h= Vertical damping between right rail and track
# czrt$2r= If czrt$2h not can be found
# czrt$1 = If czrt$2r not can be found
# czrt = If czrt$1 not can be found
#
# kyrt$2v= Lateral stiffness between left rail and track
# kyrt$2l= If kyrt$2v not can be found
# kyrt$1 = If kyrt$2l not can be found
# kyrt = If kyrt$1 not can be found
#
# kzrt$2v= Vertical stiffness between left rail and track
# kzrt$2l= If kzrt$2v not can be found
# kzrt$1 = If kzrt$2l not can be found
# kzrt = If kzrt$1 not can be found
#
# kzrt$2v.F0= Pre-stress force in spring kzrt$2v
# kzrt$2l.F0= If kzrt$2v.F0 not can be found
# kzrt$1.F0 = If kzrt$2l.F0 not can be found
# kzrt.F0 = If kzrt$1.F0 not can be found
# 0. = If kzrt.F0 not can be found
#
# cyrt$2v= Lateral damping between left rail and track
# cyrt$2l= If cyrt$2v not can be found
# cyrt$1 = If cyrt$2l not can be found
# cyrt = If cyrt$1 not can be found
#
# czrt$2v= Vertical damping between left rail and track
# czrt$2l= If czrt$2v not can be found
# czrt$1 = If czrt$2l not can be found
# czrt = If czrt$1 not can be found
#
# sv_trac = Average gauge of the track defined in the data field spv_trac
# gauge_average = If sv_trac not can be found
# sv_v$1 = Deviation in gauge for vehicle number $1
# A positive value in sv_v$1 entails a wide gauge
# A negative value in sv_v$1 entails a narrow gauge
# gauge_dev_v$1 = If sv_v$1 not can be found
# sv_v = If gauge_dev_v$1 not can be found
# gauge_dev = If sv_v not can be found
#
# YMtrac = Multiplying factor for lateral track irregularities
# ZMtrac = Multiplying factor for vertical track irregularities
# FItrac = Multiplying factor for cant irregularities
# CMtrac = If FItrac not can be found
# SVtrac = Multiplying factor for gauge irregularities
# GMtrac = If SVtrac not can be found
#
# lsa_$2.pn= Position along the track of the linear local
# coordinate system lsa_$2
# axl_$2.? = Variables containing the motions of the axle (wheelset)
# trc_$2.? = Variables containing the motions of the track
#
# Track irregularities will be interpolated from data fields with the
# following names:
#
# lat_trac = Lateral irregularities (positive direction = right)
# vert_trac = Vertical irregularities (positive direction = down)
# fi_trac = Cant irregularities (positive direction = positive
# rotation round the x-axle)
# spv_trac = Gauge irregularities (positive direction = wide gauge)
# gauge_trac= If spv_trac not can be found
#
#
# Wheel-rail geometry functions describing the tread of the right
# wheel and rail shall be defined by command `func intpl_r` and must
# have one of the following names:
#
# drkp_r$2h = Change in rolling radius, tread, right wheel
# drkp_r$2r = If drkp_r$2h not can be found
# cpt_$2r.drfn = If drkp_r$2r not can be found
# drkp_r$1 = If cpt_$2r.drfn not can be found
# cpt_$1.drfn = If drkp_r$1 not can be found
# drkp_r = If cpt_$1.drfn not can be found
# cpt_.drfn = If drkp_r not can be found
#
# gamkph_r$2h = Angle of contact surface, tread, right wheel
# gamkph_r$2r = If gamkph_r$2h not can be found
# cpt_$2r.gamfn= If gamkph_r$2r not can be found
# gamkph_r$1 = If cpt_$2r.gamfn not can be found
# cpt_$1.gamfn = If gamkph_r$1 not can be found
# gamkph_r = If cpt_$1.gamfn not can be found
# cpt_.gamfn = If gamkph_r not can be found
#
# zkp_r$2h = Wheel lift, tread, right wheel
# zkp_r$2r = If zkp_r$2h not can be found
# cpt_$2r.zfn = If zkp_r$2r not can be found
# zkp_r$1 = If cpt_$2r.zfn not can be found
# cpt_$1.zfn = If zkp_r$1 not can be found
# zkp_r = If cpt_$1.zfn not can be found
# cpt_.zfn = If zkp_r not can be found
#
# rohr_r$2h = Wheel-rail curvature difference, tread, right wheel
# rohr_r$2r = If rohr_r$2h not can be found
# cpt_$2r.rofn = If rohr_r$2r not can be found
# rohr_r$1 = If cpt_$2r.rofn not can be found
# cpt_$1.rofn = If rohr_r$1 not can be found
# rohr_r = If cpt_$1.rofn not can be found
# cpt_.rofn = If rohr_r not can be found
#
# poskph_r$2h = Position of contact surface of wheel, tread, right wheel
# poskph_r$2r = If poskph_r$2h not can be found
# cpt_$2r.poswfn= If poskph_r$2r not can be found
# poskph_r$1 = If cpt_$2r.poswfn not can be found
# cpt_$1.poswfn = If poskph_r$1 not can be found
# poskph_r = If cpt_$1.poswfn not can be found
# cpt_.poswfn = If poskph_r not can be found
# 0. = If cpt_.poswfn not can be found
#
#
# Wheel-rail geometry functions describing the tread of the left
# wheel and rail shall be defined by command `func intpl_r` and must
# have one of the following names:
#
# drkp_r$2v = Change in rolling radius, tread, left wheel
# drkp_r$2l = If drkp_r$2v not can be found
# cpt_$2l.drfn = If drkp_r$2l not can be found
# drkp_r$1 = If cpt_$2l.drfn not can be found
# cpt_$1.drfn = If drkp_r$1 not can be found
# drkp_r = If cpt_$1.drfn not can be found
# cpt_.drfn = If drkp_r not can be found
#
# gamkpv_r$2v = Angle of contact surface, tread, left wheel
# gamkpv_r$2l = If gamkph_r$2v not can be found
# cpt_$2l.gamfn= If gamkph_r$2l not can be found
# gamkpv_r$1 = If cpt_$2l.gamfn not can be found
# cpt_$1.gamfn = If gamkph_r$1 not can be found
# gamkpv_r = If cpt_$1.gamfn not can be found
# cpt_.gamfn = If gamkph_r not can be found
#
# zkp_r$2v = Wheel lift, tread, left wheel
# zkp_r$2l = If zkp_r$2v not can be found
# cpt_$2l.zfn = If zkp_r$2l not can be found
# zkp_r$1 = If cpt_$2l.zfn not can be found
# cpt_$1.zfn = If zkp_r$1 not can be found
# zkp_r = If cpt_$1.zfn not can be found
# cpt_.zfn = If zkp_r not can be found
#
# rohr_r$2v = Wheel-rail curvature difference, tread, left wheel
# rohr_r$2l = If rohr_r$2v not can be found
# cpt_$2l.rofn = If rohr_r$2l not can be found
# rohr_r$1 = If cpt_$2l.rofn not can be found
# cpt_$1.rofn = If rohr_r$1 not can be found
# rohr_r = If cpt_$1.rofn not can be found
# cpt_.rofn = If rohr_r not can be found
#
# poskph_r$2v = Position of contact surface of wheel, tread, left wheel
# poskph_r$2l = If poskph_r$2v not can be found
# cpt_$2l.poswfn= If poskph_r$2l not can be found
# poskph_r$1 = If cpt_$2l.poswfn not can be found
# cpt_$1.poswfn = If poskph_r$1 not can be found
# poskph_r = If cpt_$1.poswfn not can be found
# cpt_.poswfn = If poskph_r not can be found
# 0. = If cpt_.poswfn not can be found
#
# Wheel-rail geometry functions describing the flange of the right
# wheel and rail shall be defined by command `func intpl_r` and must
# have one of the following names:
#
# drkp_f$2h = Change in rolling radius, flange right wheel
# drkp_f$2r = If drkp_f$2h not can be found
# cpf_$2r.drfn = If drkp_f$2r not can be found
# drkp_f$1 = If cpf_$2r.drfn not can be found
# cpf_$1.drfn = If drkp_f$1 not can be found
# drkp_f = If cpf_$1.drfn not can be found
# cpf_.drfn = If drkp_f not can be found
#
# gamkph_f$2h = Angle of contact surface, flange right wheel
# gamkph_f$2r = If gamkph_f$2h not can be found
# cpf_$2r.gamfn= If gamkph_f$2r not can be found
# gamkph_f$1 = If cpf_$2r.gamfn not can be found
# cpf_$1.gamfn = If gamkph_f$1 not can be found
# gamkph_f = If cpf_$1.gamfn not can be found
# cpf_.gamfn = If gamkph_f not can be found
#
# zkp_f$2h = Wheel lift, flange right wheel
# zkp_f$2r = If zkp_f$2h not can be found
# cpf_$2r.zfn = If zkp_f$2r not can be found
# zkp_f$1 = If cpf_$2r.zfn not can be found
# cpf_$1.zfn = If zkp_f$1 not can be found
# zkp_f = If cpf_$1.zfn not can be found
# cpf_.zfn = If zkp_f not can be found
#
# rohr_f$2h = Wheel-rail curvature difference, flange right wheel
# rohr_f$2r = If rohr_f$2h not can be found
# cpf_$2r.rofn = If rohr_f$2r not can be found
# rohr_f$1 = If cpf_$2r.rofn not can be found
# cpf_$1.rofn = If rohr_f$1 not can be found
# rohr_f = If cpf_$1.rofn not can be found
# cpf_.rofn = If rohr_f not can be found
#
# poskph_f$2h = Position of contact surface of wheel, flange right wheel
# poskph_f$2r = If poskph_f$2h not can be found
# cpf_$2r.poswfn= If poskph_f$2r not can be found
# poskph_f$1 = If cpf_$2r.poswfn not can be found
# cpf_$1.poswfn = If poskph_f$1 not can be found
# poskph_f = If cpf_$1.poswfn not can be found
# cpf_.poswfn = If poskph_f not can be found
# 0. = If cpf_.poswfn not can be found
#
# Wheel-rail geometry functions describing the flange of the left
# wheel and rail shall be defined by command `func intpl_r` and must
# have one of the following names:
#
# drkp_f$2h = Change in rolling radius, flange left wheel
# drkp_f$2r = If drkp_f$2h not can be found
# cpf_$2r.drfn = If drkp_f$2r not can be found
# drkp_f$1 = If cpf_$2r.drfn not can be found
# cpf_$1.drfn = If drkp_f$1 not can be found
# drkp_f = If cpf_$1.drfn not can be found
# cpf_.drfn = If drkp_f not can be found
#
# gamkph_f$2h = Angle of contact surface, flange left wheel
# gamkph_f$2r = If gamkph_f$2h not can be found
# cpf_$2r.gamfn= If gamkph_f$2r not can be found
# gamkph_f$1 = If cpf_$2r.gamfn not can be found
# cpf_$1.gamfn = If gamkph_f$1 not can be found
# gamkph_f = If cpf_$1.gamfn not can be found
# cpf_.gamfn = If gamkph_f not can be found
#
# zkp_f$2h = Wheel lift, flange left wheel
# zkp_f$2r = If zkp_f$2h not can be found
# cpf_$2r.zfn = If zkp_f$2r not can be found
# zkp_f$1 = If cpf_$2r.zfn not can be found
# cpf_$1.zfn = If zkp_f$1 not can be found
# zkp_f = If cpf_$1.zfn not can be found
# cpf_.zfn = If zkp_f not can be found
#
# rohr_f$2h = Wheel-rail curvature difference, flange left wheel
# rohr_f$2r = If rohr_f$2h not can be found
# cpf_$2r.rofn = If rohr_f$2r not can be found
# rohr_f$1 = If cpf_$2r.rofn not can be found
# cpf_$1.rofn = If rohr_f$1 not can be found
# rohr_f = If cpf_$1.rofn not can be found
# cpf_.rofn = If rohr_f not can be found
#
# poskph_f$2v = Position of contact surface of wheel, flange left wheel
# poskph_f$2l = If poskph_f$2v not can be found
# cpf_$2l.poswfn= If poskph_f$2l not can be found
# poskph_f$1 = If cpf_$2l.poswfn not can be found
# cpf_$1.poswfn = If poskph_f$1 not can be found
# poskph_f = If cpf_$1.poswfn not can be found
# cpf_.poswfn = If poskph_f not can be found
# 0. = If cpf_.poswfn not can be found
#
#
# Output:
# =======
# This substructure file will generate the following variables:
#
# tral$2.y = Lateral track irregularities, track center line
# tral$2r.y = Lateral track irregularities, right rail
# tral$2l.y = Lateral track irregularities, left rail
# tral$2r.z = Vertical track irregularities, right rail
# tral$2l.z = Vertical track irregularities, left rail
# tral$2r.vy = First order derivative of lateral track irregularities, right rail
# tral$2l.vy = First order derivative of lateral track irregularities, left rail
# tral$2r.vz = First order derivative of vertical track irregularities, right rail
# tral$2l.vz = First order derivative of vertical track irregularities, left rail
# tral$2.f = Track irregularities in roll direction
# tral$2r.k = Pitch irregularities right rail
# tral$2l.k = Pitch irregularities left rail
# tral$2r.p = Yaw irregularities right rail
# tral$2l.p = Yaw irregularities left rail
#
# cp_$2r.eta = Lateral shift between wheel and rail right wheel
# cpt_$2r.ksi= Longitudinal position of the contact point, tread right wheel
# cpt_$2r.dr = Change in wheel radius due to wheel-rail geometry, tread right wheel
# cpt_$2r.gam= Angle of contact point, tread right wheel
# cpt_$2r.z = Vertical motion of wheel due to wheel-rail geometry, tread right wheel
# cpt_$2r.irx= Longitudinal wheel-rail curvature difference, tread right wheel
# cpt_$2r.iry= Lateral wheel-rail curvature difference, tread right wheel
# cpt_$2r.bo = Lateral distance track center line to contact point, tread right wheel
#
# cpt_$2r.nux = Longitudinal creepage, the division has been made with regard to the
# longitudinal velocity of the car-body, tread right wheel
# cpt_$2r.nuy = Lateral creepage, the division has been made with regard to the
# longitudinal velocity of the car-body, tread right wheel
# cpt_$2r.spin= Spin creepage, the division has been made with regard to the
# longitudinal velocity of the car-body, tread right wheel
# cpt_$2r.nuxm= Modified long. creepage, due to greasy contact surface, tread right wheel
# cpt_$2r.nuym= Modified lat. creepage, due to greasy contact surface, tread right wheel
# cpt_$2r.spim= Modified spin creepage, due to greasy contact surface, tread right wheel
#
# cpt_$2r.a/b = The a/b-ratio of the contact ellipse, tread right wheel
# cpt_$2r.c = Geom. average radius of the contact ellipse c=sqrt(a*b), tread right wheel
#
# cpt_$2r.Fn = Contact force perpendicular to the contact surface, tread right wheel
# cpt_$2r.Fny = Lateral creep force in plane with the contact surface, tread right wheel
# cpt_$2r.Fx = Longitudinal force between wheel and track, tread right wheel
# cpt_$2r.Fy = Lateral force between wheel and track, tread right wheel
# cpt_$2r.Fz = Vertical force between wheel and track, tread right wheel
#
# ral_$2r.y = Lateral position of the massless rail-head, tread right wheel
# ral_$2r.vy = Lateral velocity of the massless rail-head, tread right wheel
# ral_$2r.z = Vertical position of the massless rail-head, tread right wheel
# ral_$2r.vz = Vertical velocity of the massless rail-head, tread right wheel
#
# The same naming convention has also been used on the wheels on the left hand side
# of the vehicle. The names of the variables on the left side can be created by
# changing the name ral_$2r.* into ral_$2l.*.
#
# If the wheel-rail geometry functions contain a separate description of the flange
# (two-point contact), an equal amount of variables as for the tread will be generated.
# The names of the variables will be the same as for the tread, just the introduction
# cpt_* will be changed into cpf_*.
# ---------------------------------------------------------------------------------------
#
#
#
# $1=car #, $2=axle #
substruct trac_irr [
note
note
note Start of substructure file t_irr_p1.ins Rev.980212
note -------------------------------------------------------
#
#
# Set constraints to axle and track
# ---------------------------------
constr fix_free_1 axl_$2.k 0.
#
constr fix_rigid_1 trc_$2 x 0. constr fix_rigid_1 trc_$2 z 0.
constr fix_rigid_1 trc_$2 f 0. constr fix_rigid_1 trc_$2 k 0.
constr fix_rigid_1 trc_$2 p 0.
#
func t_irr_p1 $2 $1
#
#
cpt_$2r trc_$2 cpt_$2r.ksi cpt_$2r.bo 0.
axl_$2 cpt_$2r.ksi cpt_$2r.bo 0.
# cpt_$2r trc_$2 cpt_$2r.ksi { Boh | Bo | bo } 0.
# axl_$2 cpt_$2r.ksi { Boh | Bo | bo } 0.
0. fzwrnr$2 cpt_$2r.gam cpt_$2r.nuxm cpt_$2r.nuym cpt_$2r.spim
# 1/Rh+1/Rr 1/Ryh+1/Ryr my
cpt_$2r.irx cpt_$2r.iry { kph$2.my | cpt_$2r.mu | myh | mutr | my | mu }
# E nu jvkt.rikt
2.05e11 .25 z
#
cpt_$2l trc_$2 cpt_$2l.ksi cpt_$2l.bo 0.
axl_$2 cpt_$2l.ksi cpt_$2l.bo 0.
# cpt_$2l trc_$2 cpt_$2l.ksi { Bov | -Bo | -bo } 0.
# axl_$2 cpt_$2l.ksi { Bov | -Bo | -bo } 0.
0. fzwrnl$2 cpt_$2l.gam cpt_$2l.nuxm cpt_$2l.nuym cpt_$2l.spim
# 1/Rh+1/Rr 1/Ryh+1/Ryr my
cpt_$2l.irx cpt_$2l.iry { kpv$2.my | cpt_$2l.mu | myv | mutl | my | mu }
# E nu jvkt.rikt
2.05e11 .25 z
#
#
cpf_$2r trc_$2 cpf_$2r.ksi cpf_$2r.bo 0.
axl_$2 cpf_$2r.ksi cpf_$2r.bo 0.
# cpf_$2r trc_$2 cpf_$2r.ksi { Boh | Bo | bo } 0.
# axl_$2 cpf_$2r.ksi { Boh | Bo | bo } 0.
0. fzfrnr$2 cpf_$2r.gam cpf_$2r.nuxm cpf_$2r.nuym cpf_$2r.spim
# 1/Rh+1/Rr 1/Ryh+1/Ryr
cpf_$2r.irx cpf_$2r.iry
{ kfh$2.my | cpf_$2r.mu | myfh | mufr | myf | muf | my | mu }
# E nu jvkt.rikt
2.05e11 .25 z
#
cpf_$2l trc_$2 cpf_$2l.ksi cpf_$2l.bo 0.
axl_$2 cpf_$2l.ksi cpf_$2l.bo 0.
# cpf_$2l trc_$2 cpf_$2l.ksi { Bov | -Bo | -bo } 0.
# axl_$2 cpf_$2l.ksi { Bov | -Bo | -bo } 0.
0. fzfrnl$2 cpf_$2l.gam cpf_$2l.nuxm cpf_$2l.nuym cpf_$2l.spim
# 1/Rh+1/Rr 1/Ryh+1/Ryr
cpf_$2l.irx cpf_$2l.iry
{ kfv$2.my | cpf_$2l.mu | myfv | mufl | myf | muf | my | mu }
# E nu jvkt.rikt
2.05e11 .25 z
#
# func print06_all cpt_$2r.nux
# func print06_all cpt_$2r.nuy
# func print06_all cpt_$2r.spi
# func print06_all cpt_$2r.Fnx
# func print06_all cpt_$2l.Fn
# func print06_all cpt_$2l.Fnx
# func print06_all cpt_$2l.Fny
# func print06_all cpt_$2r.Fn
# func print06_all cpt_$2r.Fny
#
note
note -------------------------------------------
note End of substructure file t_irr_p1.ins
note -------------------------------------------
]