An automated procedure for determining asymptotic elastic stress fields at singular points

An analytical tool has been developed for determining the nature of the stress and displacement fields near a fairly general singular point in linear elasticity, using the asymptotic method due to Williams. The basic mathematics for this procedure is given in Section 11.2 of J.R.Barber, Elasticity, Kluwer, Dordrecht 2nd edn. (2002), 410pp.

A more detailed description of the procedure, including detailed instructions for using the analytical tool has been published in the Journal of Strain Analysis for Engineering Design and can be downloaded here.

The user is prompted to input the local geometry of the system, the material properties and the boundary conditions (and interface conditions in the case of composite bodies or problems involving contact between two or more bodies). The tool then computes the dominant eigenvalue and provides as output the equations defining the singular stress and displacement fields and contour plots of these fields. No knowledge of the asymptotic analysis procedure is required of the user.

The tool is written in the software code MATLAB v7.0 with the MATLAB GUI development environment (GUIDE) v2.5 and the MATLAB Symbolic Toolbox v3.1. It provides a graphic interface in which users can define their problem, determine the order of the corresponding singularity and generate the distribution of stress and displacement. Final results are provided in both text and graphic format.

To download the source code, click on this link and unzip the downloaded file. If and only if you have difficulty opening this link, try this one . After downloading and unzipping the file, open MATLAB and start the program with the command `ws'.

Two example programs are included: `williams.wat' and `bogy.wat', which solve the problems of the single wedge and the bi-material wedge respectively.

Please report any problems with the software or any suggestions for additional features or improvements to Donghee Lee and J.R.Barber

You are visitor number

Back to J.R.Barber's homepage.