#
#
head 2 "Test of different wheel/rail couplings cwr_coupl= $cwr_coupl"
head 3 "axl_111 pure long.creep; axl_112 pure lat.creep; axl_113 pure spin"
head 4 "Cylindric wheel profile mu= $mu_"
##
## Analyse specific input data
##
idebug= 1
fresp_param= Fourier_CG1 1.e-4 fstart .2 fstop 20. fstep -1.05
modal_param= Schur_fact1 1.e-4
quasi_param= Damped_Tens 1.e-4 0.1 20 2.
#
# tsim_param = heun tstart 0. tstop 500 tstep 1.00 tout 1.00
# tsim_param = heun_u .5 tstart 0. tstop 2.50 tstep 0.001 tout .005
# tsim_param = heun tstart 0. tstop 2.50 tstep 0.005 tout .005
# tsim_param = heun tstart 0. tstop 5.00 tstep 0.005 tout .005
tsim_param = heun_u .05 .2e-3 tstart 0. tstop 3.00 tstep 0.0001 tout .005
##
## Track irregularities and designed track geometry
##
func intpl_r lat_trac -200 0. 0 0.
func intpl_r spv_trac -200 1435. 0 1435.
func intpl_r vert_trac -200 0. 0 0.
func intpl_r fi_trac -200 0. 0 0.
# func intpl_r krok_R -200 0 10 0 20 .005 100 .005
func intpl_r krok_R -200 0 0 0 #.001 .005 100 .005
func intpl_r krok_h -200 0 0 0
func intpl_r krok_z -200 0 0 0
##
## Wheel-rail geometry
##
# insert file $genkpf/kona_0.000_uf.kpfr
# in_substruct kpf_kona_0.000 [ 1 ] # Same wheel-rail geometry in all wheels vehicle #1
#
insert file $genkpf/kpf_rkona.kpfr
in_substruct kpf_rkona [ 1 0. ] # Same wheel-rail geometry in all wheels vehicle #1
##
## Define constants for the wheel-rail-coupling substructure wr_coupl
##
func const_block
Vo = 10. # Speed of the vehicles
mu_= 0.5 # Friction coefficient, all wheels in the train
ro_= 0.45 # Wheel radius, all wheels in the train
bo = 0.75 # Semi distance between the wheel profile origins, all wheels in the train
YMtrac= 1e-3 # From mm to m; Scale factor 1.00; Lateral irregularities
ZMtrac= 1e-3 # From mm to m; Scale factor 1.00; Vertical irregularities
GMtrac= 1e-3 # From mm to m; Scale factor 1.00; Gauge irregularities
CMtrac= 1e-3/(2*bo) # From mm to rad; Scale factor 1.00; Cant irregularities
gauge_average= 1435. # Average gauge of the track
gauge_dev_ = 0.000 # Deviation in gauge
end_block
##
## Create the wheelsets.
## ---------------------
lsys e_abs_bendrf esys Vo 0. krok_R krok_h krok_z
4 4 4 4 4 4
substruct create_axles [
lsys l_local lsa_$1 esys 0. 0. 0.
mass m_rigid_6 axl_$1 lsa_$1 0. 0. -ro_ 1887 1887 1887 720 100 720
mass m_rigid_6 trc_$1 lsa_$1 0. 0. 0. 1.e3 1.e3 1.e3 1.e3 1.e3 1.e3
mass fixpoint_6 grd_$1 lsa_$1 0. 0. 0.
body w_set_lsys axl_$1 lsa_$1 ro_ bo
body figure_lsys trc_$1 lsa_$1 1 5
.5 1 0 .5 -1 0 -.5 -1 0 -.5 1 0 .5 1 0
#
constr fix_rigid_1 axl_$1 x 0.
constr fix_rigid_1 axl_$1 y 0.010
# constr fix_rigid_1 axl_$1 z 0.
# constr fix_rigid_1 axl_$1 f 0.
constr fix_free_1 axl_$1.vk 0.
constr fix_rigid_1 axl_$1 p 0.
constr fix_rigid_1 trc_$1 y 0.
#
force rel_lsys1 vload_$1 axl_$1 0. 0. -ro_ 0. 0. 200e3-1887*9.81 0. 0. 0.
# coupl p_lin czbap= 0. 50e3
# coupl c czba_$1r axl_$1 0 bo -ro_
# grd_$1 0 bo -ro_
# czbap esys z
# coupl c czba_$1l axl_$1 0 -bo -ro_
# grd_$1 0 -bo -ro_
# czbap esys z
]
#
in_substruct create_axles [ 111 ]
in_substruct create_axles [ 112 ]
in_substruct create_axles [ 113 ]
#
mass fixpoint_6 car_1 lsa_111 0. 0. 0.
#
# coupl k_lin ip_coupl axl_111 x 0. 0. -ro_
# trc_111 x 0. 0. -ro_ 17
# coupl p_kfrkc p_name +-`F0 +-`ke +-`Ffm +-`x2 +-`kv +-`c
# coupl p_kfrkc p_rubb 0 100 10 .01 0 1
# coupl kfrkc c_rubb axl_111 0 0 -ro_
# trc_111 0 0 -ro_
# p_rubb esys y
#
#coupl p_lin stifx 0. 500e3
#coupl p_lin stify 0. 500e3
#coupl p_lin stifz -100e3 1e6
#coupl k stiff axl_111 0 0 oper_init( -ro_ - .1 )
# trc_111 0 0 oper_init( -ro_ + .1 )
# stifx esys y
#force rel_lsys1 ipforc axl_111 0. 0. 0. 100 0 0 0 0 0
# coupl k_air3_exp p_airip axl_111 0 0 -ro_
# trc_111 0 0 -ro_
# esys m
# stifx 0.
# stify 0.
# stifz
# 10 .01 10 .01
# 0 1
# 0 1 1 10
##
## Define pure creepage on the three different axles
## -------------------------------------------------
# func operp axl_111.vx= -.05
# func operp axl_111.vx= time * -.25
func operp axl_111.vk= time * -.25
# func operp axl_112.vy= -.25
# func operp axl_112.vy= time * -.25
func operp axl_112.p = time * .1
# func operp axl_113.vp= time * -1.e+2
# func operp axl_113.vx= time * -.75e+2
# func operp axl_113.vp= time * -50
# func operp axl_113.vx= time * -37.5
func operp axl_113.vp= time * -3.4
func operp axl_113.vx= time * -2.55
# func operp axl_113.vp= 0.
# func operp axl_113.vx= 0.
constr fix_rigid_1 axl_111 f 0.
constr fix_rigid_1 axl_112 f 0.
constr fix_rigid_1 axl_113 f 0.
# func operp axl_111.vx= 500 * .5e-3
# func operp axl_112.vy= 500 * .5e-3
# func operp axl_113.vp= 500 * 1.e-1
# func operp axl_113.vx= 500 * .75e-1
##
## Calculate creepage and creep-forces
## -----------------------------------
# func char cwr_coupl= wr_coupl_nl1
# func char cwr_coupl= wr_coupl_npol1
# func char cwr_coupl= wr_coupl_ne1
# func char cwr_coupl= wr_coupl_nr1
# func char cwr_coupl= wr_coupl_nra1
# func char cwr_coupl= wr_coupl_pe1
func char cwr_coupl= wr_coupl_pe2
# func char cwr_coupl= wr_coupl_pr1
# func char cwr_coupl= wr_coupl_pra1
# func char cwr_coupl= wr_coupl_lt1
if_then_char_init cwr_coupl .eq. wr_coupl_nl1 .or.
cwr_coupl .eq. wr_coupl_npol1 .or.
cwr_coupl .eq. wr_coupl_ne1 .or.
cwr_coupl .eq. wr_coupl_nr1 .or.
cwr_coupl .eq. wr_coupl_nra1 # {
#
func const knwr.F0_= 100e3 # all wheels
#
func intpl_r knwr_ # Stiffness Tread-Rail all wheels
-1.57 42e6
0. 180e6
1.57 42e6
func intpl_r cnwr_ # Damping Tread-Rail all wheels
-1.57 60e3
0. 250e3
1.57 60e3
func intpl_r knfr_ # Stiffness Flange-Rail all wheels
-1.57 42e6
0. 180e6
1.57 42e6
func intpl_r cnfr_ # Damping Flange-Rail all wheels
-1.57 60e3
0. 250e3
1.57 60e3
endif # }
#
if_then_char_init cwr_coupl .eq. wr_coupl_nl1
insert file $gensys/calc/insert_files/wr_coupl_nl1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_npol1
insert file $gensys/calc/insert_files/wr_coupl_npol1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_ne1
insert file $gensys/calc/insert_files/wr_coupl_ne1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_nr1
insert file $gensys/calc/insert_files/wr_coupl_nr1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_nra1
insert file $gensys/calc/insert_files/wr_coupl_nra1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
##
##
##
if_then_char_init cwr_coupl .eq. wr_coupl_pe1 .or.
cwr_coupl .eq. wr_coupl_pe2 .or.
cwr_coupl .eq. wr_coupl_pr1 .or.
cwr_coupl .eq. wr_coupl_pra1 .or.
cwr_coupl .eq. wr_coupl_lt1 # {
#
func const knwr.F0_= 100e3 # all wheels
#
func intpl_r knwr_ # Stiffness Tread-Rail all wheels
-1.57 600e6
0. 600e6
1.57 600e6
# func intpl_r cnwr_ # Damping Tread-Rail all wheels
# -1.57 600e3
# 0. 600e3
# 1.57 600e3
func intpl_r knfr_ # Stiffness Flange-Rail all wheels
-1.57 600e6
0. 600e6
1.57 600e6
# func intpl_r cnfr_ # Damping Flange-Rail all wheels
# -1.57 600e3
# 0. 600e3
# 1.57 600e3
func const_block
pknwr = 1 # Linjär styvhet i kontaktpunkten
kyrt_ = 42e6*100.
kzrt_ = 180e6*100.
kzrt.F0_= knwr.F0_
cyrt_ = 205e3*10#31.6
czrt_ = 250e3*10#31.6
end_block
endif # }
if_then_char_init cwr_coupl .eq. wr_coupl_lt1 # {
func intpl_r cnwr_ # Damping Tread-Rail all wheels
-1.57 60e3
0. 250e3
1.57 60e3
func intpl_r cnfr_ # Damping Flange-Rail all wheels
-1.57 60e3
0. 250e3
1.57 60e3
endif # }
##
if_then_char_init cwr_coupl .eq. wr_coupl_pe1
insert file $gensys/calc/insert_files/wr_coupl_pe1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_pe2
insert file $gensys/calc/insert_files/wr_coupl_pe2.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_pr1
insert file $gensys/calc/insert_files/wr_coupl_pr1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
if_then_char_init cwr_coupl .eq. wr_coupl_pra1
insert file $gensys/calc/insert_files/wr_coupl_pra1.ins
in_substruct wr_coupl [ 1 111 ]
in_substruct wr_coupl [ 1 112 ]
in_substruct wr_coupl [ 1 113 ]
endif
##
## Modal parameters of the wheelsets.
## ----------------------------------
# mass m_flex_1 axl_111 70 0.10 ## Torsionsform
# coupl m_flex_1 cpt_111r end_2 0.0636396 0. 0. # 0. 0.14142 0.
# coupl m_flex_1 cpt_111l end_2 -.0636396 0. 0. # 0. 0.14142 0.
# s_var sngl axl_111.c1
# s_var sngl axl_111.v1
##
## Store variables for postprocessing
## ----------------------------------
#
func print06_init cpt_111r.a_b # Print the initial value of a variable to standard output
#
s_var scalar_0 mu_
#
s_var sngl ro_
s_var sngl bo
s_var sngl Vo
#
s_var var_0 esys.x
s_var var_0 esys.y
s_var var_0 esys.z
s_var var_0 esys.f
s_var var_0 esys.k
s_var var_0 esys.p
#
s_var var_0 esys.vx
s_var var_0 esys.vy
s_var var_0 esys.vz
s_var var_0 esys.vf
s_var var_0 esys.vk
s_var var_0 esys.vp
#
s_var var_0 esys.ax
s_var var_0 esys.ay
s_var var_0 esys.az
s_var var_0 esys.af
s_var var_0 esys.ak
s_var var_0 esys.ap
#
substruct store_axle [
s_var var_0 axl_$1.x
s_var var_0 axl_$1.y
s_var var_0 axl_$1.z
s_var var_0 axl_$1.f
s_var var_0 axl_$1.k
s_var var_0 axl_$1.p
#
s_var var_0 axl_$1.vx
s_var var_0 axl_$1.vy
s_var var_0 axl_$1.vz
s_var var_0 axl_$1.vf
s_var var_0 axl_$1.vk
s_var var_0 axl_$1.vp
#
s_var var_0 axl_$1.Fx
s_var var_0 axl_$1.Fy
s_var var_0 axl_$1.Fz
s_var var_0 axl_$1.Mf
s_var var_0 axl_$1.Mk
s_var var_0 axl_$1.Mp
#
s_var var_0 axl_$1.Ax
s_var var_0 axl_$1.Ay
s_var var_0 axl_$1.Az
s_var var_0 axl_$1.Af
s_var var_0 axl_$1.Ak
s_var var_0 axl_$1.Ap
#
s_var var_0 trc_$1.x s_var var_0 trc_$1.y s_var var_0 trc_$1.z
s_var var_0 trc_$1.f s_var var_0 trc_$1.k s_var var_0 trc_$1.p
s_var var_0 trc_$1.vx s_var var_0 trc_$1.vy s_var var_0 trc_$1.vz
s_var var_0 trc_$1.vf s_var var_0 trc_$1.vk s_var var_0 trc_$1.vp
]
#
substruct store_wheel [
s_var var_0 ral_$1.y
s_var var_0 ral_$1.z
s_var var_0 ral_$1.vy
s_var var_0 ral_$1.vz
#
s_var var_0 cp_$1.eta
#
s_var var_0 tral$1.y
s_var var_0 tral$1.z
s_var var_0 tral$1.vy
s_var var_0 tral$1.vz
s_var var_0 cpt_$1.z
s_var var_0 cpt_$1.dr
s_var var_0 cpt_$1.gam
s_var var_0 cpt_$1.c
s_var var_0 cpt_$1.a_b
s_var var_0 cpt_$1.posw
s_var var_0 cpt_$1.bo
#
s_var var_0 cpt_$1.nux
s_var var_0 cpt_$1.nuy
s_var var_0 cpt_$1.spin
#
s_var var_0 cpt_$1.vux
s_var var_0 cpt_$1.vuy
s_var var_0 cpt_$1.vuz
s_var var_0 cpt_$1.Vref
#
s_var var_0 cpt_$1.Fn
s_var var_0 cpt_$1.Fny
s_var var_0 cpt_$1.Fx
s_var var_0 cpt_$1.Fy
s_var var_0 cpt_$1.Fz
s_var var_0 cpt_$1.F1x
s_var var_0 cpt_$1.F1y
s_var var_0 cpt_$1.F1z
s_var var_0 cpt_$1.F2x
s_var var_0 cpt_$1.F2y
s_var var_0 cpt_$1.F2z
s_var var_0 cpt_$1.M1f
s_var var_0 cpt_$1.M1k
s_var var_0 cpt_$1.M1p
s_var var_0 cpt_$1.M2f
s_var var_0 cpt_$1.M2k
s_var var_0 cpt_$1.M2p
#
s_var var_0 cpt_$1.fwk
s_var var_0 cpt_$1.wArchH
s_var var_0 cpt_$1.wArchH_1
s_var var_0 cpt_$1.wArchH_2
s_var var_0 cpt_$1.wArchH_3
s_var var_0 cpt_$1.wArchH_4
s_var var_0 cpt_$1.wArchH_5
s_var var_0 cpt_$1.wArchH_6
s_var var_0 cpt_$1.wArchH_7
s_var var_0 cpt_$1.wArchH_8
s_var var_0 cpt_$1.wArchH_9
s_var var_0 cpt_$1.wArchH_10
s_var var_0 cpt_$1.wArchH_11
s_var var_0 cpt_$1.wArchH_12
s_var var_0 cpt_$1.wArchH_13
s_var var_0 cpt_$1.wArchH_14
s_var var_0 cpt_$1.wArchH_15
s_var var_0 cpt_$1.wArchH_16
s_var var_0 cpt_$1.wArchS
s_var var_0 cpt_$1.wArchS_1
s_var var_0 cpt_$1.wArchS_2
s_var var_0 cpt_$1.wArchS_3
s_var var_0 cpt_$1.wArchS_4
s_var var_0 cpt_$1.wArchS_5
s_var var_0 cpt_$1.wArchS_6
s_var var_0 cpt_$1.wArchS_7
s_var var_0 cpt_$1.wArchS_8
s_var var_0 cpt_$1.wArchS_9
s_var var_0 cpt_$1.wArchS_10
s_var var_0 cpt_$1.wArchS_11
s_var var_0 cpt_$1.wArchS_12
s_var var_0 cpt_$1.wArchS_13
s_var var_0 cpt_$1.wArchS_14
s_var var_0 cpt_$1.wArchS_15
s_var var_0 cpt_$1.wArchS_16
]
in_substruct store_axle [ 111 ]
in_substruct store_axle [ 112 ]
in_substruct store_axle [ 113 ]
in_substruct store_wheel [ 111r ]
in_substruct store_wheel [ 112r ]
in_substruct store_wheel [ 113r ]
in_substruct store_wheel [ 111l ]
in_substruct store_wheel [ 112l ]
in_substruct store_wheel [ 113l ]
#
# in_substruct store [ 111l ]
# in_substruct store [ 112l ]
# in_substruct store [ 113l ]
#
eof
Creates a diagram of the creepage creep force relationship
## ## Input data for program MPLOT ## Page ------------------------- page_head 1 "Longitudinal creepage" ygridint= 1 xgridint= 2 overwrite= global y_mid= auto yint/cm= auto xax_yval= bot x_left= 0 xint/cm= auto Diagram 11 x_left= 0 xint/cm= 1e-3 y_bot = 0 yint/cm= 4e3 # Curve var_id= pure_creep_ts1 yvar= cpt_111r.Fx xvar= cpt_111r.nux # Curve var_id= pure_creep_fasim_1 yvar= cpt_111r.Fx xvar= cpt_111r.nux # Curve var_id= pure_creep_lt1 yvar= cpt_111r.Fx xvar= cpt_111r.nux # Curve var_id= pure_creep_pe0 yvar= cp1_111r.Fx xvar= cp1_111r.nux # Curve var_id= pure_creep_pe1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux Curve var_id= pure_creep_pe2 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_pe3 yvar= cp1_111r.Fx xvar= cp1_111r.nux # Curve var_id= pure_creep_pr1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_pra1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_ne1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_nr1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_nra1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_npol1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux # Curve var_id= pure_creep_nl1 yvar=-cpt_111r.Fx xvar=-cpt_111r.nux EndPage Page ------------------------- page_head 1 "Lateral creepage" ygridint= 1 xgridint= 2 overwrite= global y_mid= auto yint/cm= auto xax_yval= bot x_left= 0 xint/cm= auto Diagram 11 # xint/cm= .0025 # y_bot= 0. yint/cm= 10e3 x_left= 0 xint/cm= 1e-3 y_bot = 0 yint/cm= 4e3 # Curve var_id= pure_creep_ts1 yvar= cpt_112r.Fy xvar= cpt_112r.nuy # Curve var_id= pure_creep_fasim_1 yvar= cpt_112r.Fy xvar= cpt_112r.nuy # Curve var_id= pure_creep_lt1 yvar= cpt_112r.Fy xvar= cpt_112r.nuy # Curve var_id= pure_creep_pe0 yvar= cp1_112r.Fy xvar= cp1_112r.nuy # Curve var_id= pure_creep_pe1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy Curve var_id= pure_creep_pe2 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_pe3 yvar= cp1_112r.Fy xvar= cp1_112r.nuy # Curve var_id= pure_creep_pr1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_pra1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_ne1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_nr1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_nra1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_npol1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy # Curve var_id= pure_creep_nl1 yvar=-cpt_112r.Fy xvar=-cpt_112r.nuy EndPage Page ------------------------- page_head 1 "Spin creepage" ygridint= 1 xgridint= 2 overwrite= global y_mid= auto yint/cm= auto xax_yval= bot x_left= 0 xint/cm= auto Diagram 11 # xint/cm= .0025 # y_bot= 0. yint/cm= 10e3 y_bot = 0 yint/cm= 2e3 x_left= 0 xint/cm= .05 # Curve var_id= pure_creep_ts1 yvar= cpt_113r.Fy xvar= cpt_113r.spin # Curve var_id= pure_creep_fasim_1 yvar= cpt_113r.Fy xvar= cpt_113r.spin # Curve var_id= pure_creep_lt1 yvar= cpt_113r.Fy xvar= cpt_113r.spin # Curve var_id= pure_creep_pe0 yvar= cp1_113r.Fy xvar= cp1_113r.spin # Curve var_id= pure_creep_pe1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin Curve var_id= pure_creep_pe2 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_pe3 yvar= cp1_113r.Fy xvar= cp1_113r.spin # Curve var_id= pure_creep_pr1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_pra1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_ne1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_nr1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_nra1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_npol1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin # Curve var_id= pure_creep_nl1 yvar=-cpt_113r.Fy xvar=-cpt_113r.spin EndPage