disp2sif - Calculate mode I+II stress intensity factors (SIF) from surface displacements

The code implements Muskhelishvili's complex function approach to calculate 2D stress intensity factors (modes I+II) from crack tip displacement fields. Surface displacements can be used with a plane stress elastic model. The code uses power/Fourier series expansion and a conformal mapping of an ellipse onto a circle. Fortran complex intrinsics are used in the code.

Releases

18-JAN-2017: release 1.0. However, the code, documentation and tests date from 2003-2006, with minor modifications since then.

Publications using this code

• A. Steuwer, M. Rahman, A. Shterenlikht, M. E. Fitzpatrick, L. Edwards, P. J. Withers, The evolution of crack-tip stresses during a fatigue overload event, ACTA MATERIALIA, 58(2010):4039-4052. DOI: 10.1016/j.actamat.2010.03.013
• P. Lopez-Crespo, R. L. Burguete, E. A. Patterson, A. Shterenlikht, P. J. Withers, J. R. Yates, Study of a Crack at a Fastener Hole by Digital Image Correlation, EXPERIMENTAL MECHANICS, 49(2009):551-559. DOI: 10.1007/s11340-008-9161-1
• P. Lopez-Crespo, A. Shterenlikht, J. R. Yates, E. A. Patterson, P. J. Withers, Some experimental observations on crack closure and crack-tip plasticity, FATIGUE & FRACTURE OF ENGINEERING MATERIALS & STRUCTURES, 32(2009):418-429 DOI: 10.1111/j.1460-2695.2009.01345.x
• P. Lopez-Crespo, A. Shterenlikht, E. A. Patterson, J. R. Yates, P. J. Withers, The stress intensity of mixed mode cracks determined by digital image correlation, JOURNAL OF STRAIN ANALYSIS FOR ENGINEERING DESIGN, 43(2008):769-780. DOI: 10.1243/03093247JSA419
• P. Lopez-Crespo, A. Shterenlikht, E. A. Patterson, J. R. Yates, P. J. Withers, Fatigue crack monitoring using image correlation, Key Engineering Materials, 385-387(2008):341-344. DOI: 10.4028/www.scientific.net/KEM.385-387.341
• A. Shterenlikht, F. A. Diaz Garrido, P. Lopez-Crespo, P. J. Withers, E. A. Patterson, Mixed mode (I+II) stress intensity factor measurement by electronic speckle pattern interferometry and image correlation, Advances in Experimental Mechanics, 1-2(2004):107-112 DOI: 10.4028/www.scientific.net/AMM.1-2.107

Further information on the code

• A. Shterenlikht, P. J. Withers, P. Lopez-Crespo, J. R. Yates, E. A. Patterson, Mixed mode (I+II) SIF calculation using surface displacements, circa 2004-2005. An old presentation/report, slides1.pdf
• A. Shterenlikht, P. Lopez-Crespo, Calculation of mixed mode (I+II) SIFs using in-plane displacement fields near crack, circa 2004-2005. This is the most detailed report on the analytic model used in the code, report.pdf
• A. Shterenlikht, P. Lopez-Crespo, Calculation of mixed mode (I+II) SIFs using in-plane displacement fields near crack, circa 2004-2005. A shorter version of the report on the analytic model used in the code, repshort.pdf

Other Fortran related projects

• ABUMPACK - Abaqus UMAT and VUMAT for GTN, Rousselier, von Mises plasticity
• CASUP - Cellular Automata library for HPC and supercomputers
• Coarrays - Fortran 2008, 2018 coarrays course - native Fortran SPMD parallel programming
• CMPLX - Validation of Fortran 2008 complex intrinsics and minus zero on branch cuts
• OPTPACK - OPTimisation PACKage - multi-start global optimisation in Fortran 2003 driven with shell scripts
• SCNU - Side Cut Non-Uniform residual stress measurement method
• disp2sif - Calculation of mixed mode (I+II) stress intensity factors (SIF) from crack tip displacements
• erf - Error functions of complex arguments, erf(z), implemented in modern Fortran

validate this page