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.

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

- 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

- 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

- Global optimisation in Fortran 77, 90 and 2003
- Solid mechanics, Abaqus UMAT and VUMAT in Fortran 77, 90 and 2003
- Micromechanics of materials for clusters and supercomputers in Fortran 2008 with coarrays
- Inverse method for residual stress calculation in Fortran 2003 with OpenMP
- Fortran 2008, 2015 coarrays course
- Calculation of mixed mode (I+II) stress intensity factors (SIF) from crack tip displacements
- Validation of Fortran 2008 complex intrinsics and minus zero on branch cuts
- Error functions of complex arguments, erf(z), implemented in modern Fortran

Help with usage, modification and further development of this code can provided via CMPLX Engineering Ltd, an engineering consultancy which specialises in Fortran, solid mechanics and supercomputing.