module CS_fatigue_variables #ifndef dp use, intrinsic :: iso_fortran_env, only: dp=>real64 #endif implicit none public real(dp) :: residual_sig_hoop !! residual hoop stress in strucutal material (Pa) real(dp) :: n_cycle !! Allowable number of cycles for CS stress model real(dp) :: n_cycle_min !! Minimum llowable number of cycles for CS stress model real(dp) :: t_crack_radial !! Initial depth of crack in thickness of conduit (m) real(dp) :: t_crack_vertical !! Inital vertical crack size (m) real(dp) :: t_structural_radial !! Thickness of CS conductor conduit (m) real(dp) :: t_structural_vertical !! Vertical thickness of CS conductor conduit (m) real(dp) :: bkt_life_csf !! Switch to pass bkt_life cycles to n_cycle_min real(dp) :: sf_vertical_crack !! Safety factor for vertical crack size (-) real(dp) :: sf_radial_crack !! Safety factor for radial crack size (-) real(dp) :: sf_fast_fracture !! safety factor for stress intensity factor (-) real(dp) :: paris_coefficient !! Paris equation material coefficient (-) real(dp) :: paris_power_law !! Paris equation material power law (-) real(dp) :: walker_coefficient !! walker coefficent (-) real(dp) :: fracture_toughness !! fracture toughness (MPa m^1/2) contains subroutine init_CS_fatigue_variables !! Initialise module variables implicit none residual_sig_hoop = 2.4D8 t_crack_radial = 6.0D-3 t_crack_vertical = 0.89D-3 n_cycle = 0.0D0 n_cycle_min = 2.0D4 t_structural_vertical = 0.022D0 t_structural_radial = 0.07D0 bkt_life_csf = 0 sf_vertical_crack = 2.0D0 sf_radial_crack = 2.0D0 sf_fast_fracture = 1.5D0 paris_coefficient = 65.0D-14 paris_power_law = 3.5D0 walker_coefficient = 0.436D0 fracture_toughness = 2.0e2 end subroutine init_CS_fatigue_variables end module CS_fatigue_variables