Home    Search    Reference Manuals    Return   

##
##
##      Parameters for different types of analysis
##      ==========================================================
  fresp_param= Fourier_CG1 1.e-4  fstart= .05  fstop= 20.   fstep= -1.04
  modal_param= Schur_fact1 1.e-4
  quasi_param= Damped_Tens 1.e-4  .2 10  16.
  tsim_param = heun_u              tstart= .0   tstep= 0.001  tout= 0.020
#
# write_ca_file         # During simulation, write continuously to file $ident.ca

##  
##  Define substructure files, saving variables for postprocessing
##  ==============================================================
  insert file $gensys/calc/insert_files/save_car1.ins   # Save variables for carbody
  insert file $gensys/calc/insert_files/save_bol1.ins   # Save variables for bolster beam
  insert file $gensys/calc/insert_files/save_bog1.ins   # Save variables for bogie
  insert file $gensys/calc/insert_files/save_axl1.ins   # Save variables for axle
  insert file $gensys/calc/insert_files/save_creep1.ins # Save variables for the contact points
  insert file $gensys/calc/insert_files/save_mass_1.ins # Save variables for masses

##
##  Calculate wear and RCF indexes
##  ==============================================================
 substruct wear_RCF [
  func fl_wear_w $1  mu_cp1$1 mu_cp2$1 mu_cp3$1
  s_var sngl cpa_$1.Fnu    s_var sngl cp1_$1.Fnu    s_var var_0 cp2_$1.Fnu    s_var var_0 cp3_$1.Fnu
  s_var sngl cpa_$1.FMnu   s_var sngl cp1_$1.FMnu   s_var var_0 cp2_$1.FMnu   s_var var_0 cp3_$1.FMnu
#
  func lpass2q_0  cp1_$1.FMnuq=  cp1_$1.FMnu  3 .7
  func max2       cp1_$1.FMnul=  cp1_$1.FMnu  cp1_$1.FMnuq
  s_var sngl cp1_$1.FMnul
#
  if_then_init .exist. cp2_$1.FMnu
   func lpass2q_0  cp2_$1.FMnuq=  cp2_$1.FMnu  3 .7
   func max2       cp2_$1.FMnul=  cp2_$1.FMnu  cp2_$1.FMnuq
  else
   func const cp2_$1.FMnul= 0.
  endif
  s_var sngl cp2_$1.FMnul
#
  if_then_init .exist. cp3_$1.FMnu
   func lpass2q_0  cp3_$1.FMnuq=  cp3_$1.FMnu  3 .7
   func max2       cp3_$1.FMnul=  cp3_$1.FMnu  cp3_$1.FMnuq
  else
   func const cp3_$1.FMnul= 0.
  endif
  s_var sngl cp3_$1.FMnul

##
##                                    Yield    Dang  Residual Wear    Contact     Wear
##                       Name   whe_  Stress   Van   Stress   Limit   ForceLlim   Cp1       Cp2           Cp3
  func rolling_fatigue_3 RCF_$1  $1   300e6    0.32    0.    125 175   4000  cp1_$1.FMnul  cp2_$1.FMnul  cp3_$1.FMnul   
##
## YieldStress= 300e6 after some hardening
##
## Surface fatigue index    Subsurface fatigue index  Deep subsurface fatigue index 
  s_var sngl RCF_$1.FPs   s_var sngl RCF_$1.FPb    s_var sngl RCF_$1.FPd
  s_var sngl RCF_$1.FPs1  s_var sngl RCF_$1.FIs1
  s_var sngl RCF_$1.FPs2  s_var sngl RCF_$1.FIs2
  s_var sngl RCF_$1.FPs3  s_var sngl RCF_$1.FIs3

  func operp RCF_$1.FPbm= RCF_$1.FPb / 450e6 - 1
  s_var sngl RCF_$1.FPbm
 ]