Search

The interior stress problem is solved numerically for a single-edge notched specimen under uniaxial load. The algorithm is based on a modification of a Fredholm second-kind integral equation with compact operators due to Muskhelishvili. Several singular basis functions for each of the seven corners in the geometry enable high uniform resolution of the stress field with a modest number of discretization points. As a consequence, notch stress intensity factors can be computed directly from the solution. This is an improvement over other procedures where the stress field is not resolved in the corners and where notch stress intensity factors are computed in a roundabout way via a path-independent integral. Numerical examples illustrate the superior stability and economy of the new scheme.

We present a fast algorithm for the calculation of elastostatic fields in locally isotropic composites. The method uses an integral equation approach due to Sherman, combined with the fast multipole method and an adaptive quadrature technique. Accurate solutions can be obtained with inclusions of arbitrary shape at a cost roughly proportional to the number of points needed to resolve the interface. Large-scale problems, with hundreds of thousands of interface points can be solved using modest computational resources.