With any questions, bug reports, feedback and any other comments please submit a ticket.

24-JAN-2016: SCNURS 1.0, initial public release.

- 2016: H.K. Kim, M.J. Pavier, A. Shterenlikht, Measurement of Highly Non-uniform Residual Stress Fields in Thin Plate Using a New Side Cut Destructive Method, Proc. Residual Stresses 2016 – ICRS-10, Eds. Thomas M. Holden, Ondrej Muransky, Lyndon Edwards, Materials Research Proceedings 2(2016):115-120 DOI: 10.21741/9781945291173-20
- 2016: A. Shterenlikht, H. K. Kim, Forward and inverse solutions for 2D semi-infinite strip with self-equilibrated end loading, README, the most detailed explanation of the method and the code, work in progress.
- 2016: A. Shterenlikht, Finding roots of sinx+x=0 and sinx-x=0, short note, PDF
- 2016: H.K. Kim, A. Shterenlikht, M.J. Pavier, A. Velichko, N.A. Alexander, On stability of a new side cut destructive method for measuring non-uniform residual stress in thin plates, Int. J Solids Struct, 100-101(2016):223-233, DOI: 10.1016/j.ijsolstr.2016.08.019
- 2015: H. K. Kim, BSSM Young stress analysis competition talk.
- 2015: H. K. Kim, M. J. Pavier, A. Shterenlikht, Highly non-uniform residual stress measurement with reduced plastic error, Mech Eng Dept, The University of Bristol, UK, poster.
- 2015: H. K. Kim, H. E. Coules, M. J. Pavier, A. Shterenlikht, Measurement of Highly Non-Uniform Residual Stress Fields with Reduced Plastic Error, Exp. Mech. 55(7):1211-1224, DOI:10.1007/s11340-015-0025-1, PDF.
- 2014: H. K. Kim, M. J. Pavier, A. Shterenlikht, Plasticity and stress heterogeneity influence on mechanical stress relaxation residual stress measurements, Proc. of 9th European Conference on Residual Stresses (ECRS), Troyes, France, 7-10 July 2014, Advanced Materials Research, 996:249-255, 2014, DOI:10.4028/www.scientific.net/AMR.996.249, paper and talk.
- 2013: H. K. Kim, M. J. Pavier, A. Shterenlikht, Measuring locally non-uniform in-plane residual stress with straight cuts and DIC, Proc. 9th Int. Conf. Advances in Exp. Mechanics, 3-5 September 2013, Cardiff, Wales, UK, The British Society for Strain Measurements (BSSM), paper and talk.
- 2000: I. A. Razumovskii, A. L. Shterenlikht, Determining the locally-nonuniform residual-stress fields in plane parts by the sectioning method, J. Machinery Manufacture 4:40-45, translated from Russian, Allerton Press, USA, ISSN: 1052-6188.

A 2D problem is assumed (x_{1}, x_{2}).
Residual stress is assumed to be constant through thickness,
x_{3}.
The plate length is assumed to be much greater then width.
The method is as follows:

- Cut on x1=0 plane. Note that cut can be progressed along x2 or along x3, with different effects on plastic flow.
- Measure relaxation disp. fields, u1 and u2, e.g. with DIC.
- Use Mathieu series 2D elastic solution to the measured displacements.

See above our papers on the analysis of plastic flow on the measured RS fields.

Enforcing the BC gives these 2 equations:

$\mathrm{sin}2x+2x=0$

$\mathrm{sin}2x-2x=0$

Each equation has an infinite number of complex roots, which can be plotted on the complex plane:

Self-equlibrated loading is applied on the end. The forward problem is solved either analytically, or with FE. The result is u1 and u2 disp. fields. Then the inverse problem is solved and the series coefficients are calculated from the displacements using LLS or rank-deficient minimum norm LLS. We use Lapack.

In all plots the applied stress fields are shown with symbols. The reconstructed fields are shown with lines.

The inverse problem is ill-posed i.e. ill-conditioned. These are various manifestations.

5151 data points, 2 displacement values for each point, so 5151*2 = 10302 equations in total.

The matrix is 10302 x 40.

The matrix is 10302 x 200.

The matrix is 10302 x 320.

The matrix is 10302 x 400.

The matrix is 10302 x 40.

The matrix is 10302 x 160.

The matrix is 10302 x 196.

The matrix is 10302 x 200.

The matrix is 10302 x 400.

The matrix is 10302 x 40.

The matrix is 10302 x 160.

The matrix is 10302 x 196.

The matrix is 10302 x 200.

The matrix is 10302 x 320.

The matrix is 10302 x 40.

The matrix is 10302 x 160.

The matrix is 10302 x 196.

The matrix is 10302 x 200.

The matrix is 10302 x 320.

The matrix is 10302 x 400.

The matrix is 10302 x 20.

The matrix is 10302 x 320.

The matrix is 10302 x 400.

The matrix is 10302 x 40.

The matrix is 10302 x 160.

The matrix is 10302 x 196.

The matrix is 10302 x 200.

- 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