= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = STRESS INTENSITY FACTOR ANALYSIS PROGRAMS = = = = = = = = = = = = for = = = = = = = = = = = = STRESS FRACTURE MECHANICS = = = = = = = = = = = = = = = = Copyright (c) 2022 by J. Ervin Glover, Ph.D. = = = = All Rights Reserved = = = = = = = = = = = = = = = = ( ABSTRACTS ) = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = ======================================================================== ABSTRACTS OF STRESS INTENSITY FACTOR ANALYSIS PROGRAMS ======================================================================== Constraint variables in this document are defined as follows: K0 = stress intensity factor for a crack in an infinite sheet KI/K0 = complex stress functional value TIP A = left tip of a crack TIP B = right tip of a crack a = half-length of a crack (in inches) b = distance from the center of a crack to the boundary (in inches) a0 = initial value of the parameter 'a' af = final value of the parameter 'a' R = radius of a hole in a mechanical part (in inches) DELTA = increment of the crack-length (in inches) t1 = thickness of the cracked sheet (in inches) t2 = thickness of the uncracked sheet (in inches) E1 = Young's Modulus for a cracked sheet E2 = Young's Modulus for an uncracked sheet p = a fixed point in measured proximity to a crack x = horizontal distance from a fixed point (P) to the center of a crack (in inches) y = vertical distance from a fixed point (P) to the center of a crack (in inches) LF = Localized Force LM = Localized Moment ======================================================================== ======================================================================== ABSTRACTS OF STRESS INTENSITY FACTOR ANALYSIS PROGRAMS ======================================================================== -- SIF-1 -- SIF-1 determines, by piecewise linear interpolation, the various values of KI/K0 from a fixed value of 'b', an initial value of 'a', a final value of 'a', a fixed value of 'R' (the radius of a hole in a mechanical part) and a fixed increment (delta) of 'a', all of which are prompted for input at the computer terminal. (KI/K0 is for a tip of a crack near a circular hole in a sheet subjected to a uniform uniaxial stress). =================================================================== -- SIF-2 -- SIF-2 determines, by piecewise linear interpolation, the various values of KI/K0 from a fixed value of 'b', an initial value of 'a', a final value of 'a', a fixed increment 'delta' of 'a', the thick- ness 't1' of the cracked sheet, the thickness 't2' of the uncracked sheet, Young's Modulus 'E1' of the cracked sheet, and Young's Modulus 'E2' of the uncracked sheet, all of which are prompted for input at the computer terminal. (KI/K0 is for a tip of a crack near the junction of two uniformly stressed sheets). =================================================================== -- SIF-3 -- SIF-3 analyzes, by piecewise linear interpolation, the force perpendicular to a crack or the localized moment at a fixed point, and determines various values of KI/K0 from an initial value of 'a', a final value of 'a', a fixed increment 'delta' of 'a', the horizontal distance 'x' from a fixed point 'P' to the center of the crack, and the vertical distance 'y' from 'P' to the crack (in inches), all of which are prompted for input at the computer terminal. (KI/K0 is for a tip of a crack in a sheet subjected to a force perpendicular to the crack or subjected to a localized moment). =================================================================== -- SIF-4 -- SIF-4 determines, by piecewise linear interpolation, the various values of KI/K0 from a fixed value of 'b', an initial value of 'a', a final value of 'a', (and a fixed increment 'delta' of 'a', all of which are prompted for input at the computer terminal. (KI/K0 is for an internal crack perpendicular to the edge of a sheet subjected to a uniaxial tensile stress). =================================================================== -- SIF-4A -- SIF-4A determines, by piecewise linear interpolation, the various values of KI/K0 from a fixed value of 'b' and varying values of 'a', which are prompted for input at the computer terminal. (KI/K0 is for an internal crack perpendicular to the edge of a sheet subjected to a uniaxial tensile stress). =================================================================== -- SIF-4B -- SIF-4B determines, by piecewise linear interpolation, the various values of KI/K0 from varying values of'a' and 'b', which are prompted for input at the computer terminal. (KI/K0 is for an internal crack perpendicular to the edge of a sheet subjected to a uniaxial tensile stress). =================================================================== -- SIF-4C -- SIF-4C determines, by piecewise linear interpolation, the various values of KI/K0 from fixed values of 'a' and 'b' listed in an input file representing a given case and supplied by the engineer/user. (KI/K0 is for an internal crack perpendicular to the edge of a sheet subjected to a uniaxial tensile stress). =================================================================== -- SIF-5 -- SIF-5 determines, by piecewise linear interpolation, the various values of KI/K0 from a fixed value of 'b', an initial value of 'a', a final value of 'a', and a fixed increment 'delta' of 'a', all of which are prompted for input at the computer terminal. (KI/K0 is for a tip of a central crack in a finite width sheet subjected to a uniform bending moment). =================================================================== -- SIF-5A -- SIF-5A determines, by piecewise linear interpolation, the various values of KI/K0 from a fixed value of 'b' and varying values of 'a', which are prompted for input at the computer terminal. (KI/K0 is for a tip of a central crack in a finite width sheet subjected to a uniform bending moment). =================================================================== -- SIF-5B -- SIF-4B determines, by piecewise linear interpolation, the various values of KI/K0 from varying values of'a' and 'b', which are prompted for input at the computer terminal. (KI/K0 is for a central crack in a finite width sheet subjected to a uniform bending moment). =================================================================== -- SIF-5C -- SIF-5C determines, by piecewise linear interpolation, the various values of KI/K0 from fixed values of 'a' and 'b' listed in an input file representing a given case and supplied by the engineer/user. (KI/K0 is for a tip of a central crack in a finite width sheet subjected to a uniform bending moment). =================================================================== -- JIM6 -- JIM6 determines the various values of BETA6 = KI/K0 from an initial crack length 'a0' (in inches), a final crack length 'af' (in inches), an increment 'delta' of the crack length, and a radius 'R' of a hole in a mechanical part. (KI/K0 is for the BOWIE CORRECTION FOR A SINGLE FLAW AT A HOLE). =================================================================== -- JIM31 -- JIM31 determines the various values of BETA31 = KI/K0 from an initial half-crack depth 'depth' (in inches), an initial half-crack length 'rlen' (in inches), a final half-crack depth 'depthf', an increment of the crack depth 'delta', an assumed angle (80 degrees) and a half-plate thickness 'thick' (in inches). (KI/K0 is the NEWMAN/RAJU EMBEDDED ELLIPTICAL FLAW BETA SOLUTION). =================================================================== -- JIM32 -- JIM32 determines the various values of BETA32 = KI/K0 from an initial crack depth 'depth' (in inches), an initial half-crack length 'rlen' (in inches), a final crack depth 'depthf', an increment of the crack depth 'delta', an assumed angle (80 degrees) and a plate thickness 'thick' (in inches). (KI/K0 is the NEWMAN/RAJU SEMI-ELLIPTICAL SURFACE FLAW BETA SOLUTION). =================================================================== -- JIM33 -- JIM33 determines the various values of BETA33 = KI/K0 from an initial crack depth 'a' (in inches), an initial crack length 'c' (in inches), a final crack depth 'depthf', an increment of the crack depth 'delta', an assumed angle (80 degrees) and a plate thickness 't' (in inches). (KI/K0 is the NEWMAN/RAJU QUARTER-ELLIPTICAL CORNER CRACK BETA SOLUTION). =================================================================== -- JIM34 -- JIM34 determines the various values of BETA34 = KI/K0 from an initial half-crack depth 'a' (in inches), an initial half-crack length 'c' (in inches), a final half-crack depth 'depthf', an increment of the crack depth 'delta', an assumed angle (80 degrees), a half-plate thickness 't' (in inches), and a hole of radius 'R' in a mechanical part. (KI/K0 is the NEWMAN/RAJU SEMI-ELLIPTICAL SURFACE CRACK AT A HOLE BETA SOLUTION FOR ONE CRACK). =================================================================== -- JIM36 -- JIM36 determines the various values of BETA36 = KI/K0 from an initial half-crack depth 'a' (in inches), an initial half-crack length 'c' (in inches), a final half-crack depth 'depthf', an increment of the crack depth 'delta', an assumed angle (80 degrees), a half-plate thickness 't' (in inches), and a hole of radius 'R' in a mechanical part. (KI/K0 is the NEWMAN/RAJU QUARTER-ELLIPTICAL CORNER CRACK AT A HOLE BETA SOLUTION FOR ONE CRACK). =================================================================== STRESS INTENSITY FACTOR ADVISOR PROGRAM CASES Case # Program ------ ------- 1. ---------> SIF-4 (4A,4B,4C) FORTRAN 2. ---------> JIM32 FORTRAN 3. ---------> JIM31 FORTRAN 4. ---------> SIF-4 (4A,4B,4C) FORTRAN 5. ---------> JIM33 FORTRAN 6. ---------> JIM6 FORTRAN 7. ---------> JIM36 FORTRAN 8. ---------> JIM34 FORTRAN 9. ---------> SIF-3 FORTRAN 10. ---------> SIF-3 FORTRAN 11. ---------> SIF-5 (5A,5B,5C) FORTRAN 12. ---------> SIF-4 (4A,4B,4C) FORTRAN 13. ---------> SIF-1 FORTRAN 14. ---------> SIF-2 FORTRAN ===================================================================