Green’s functions and boundary element method formulation for 3D anisotropic media

The implementation of Wang’s theoretical solution is presented for elastostatic displacement Green’s function for three-dimensional solids of general anisotropy. Excerpts from the authors’ fortran code are included. A numerical algorithm for the calculation of the derivatives of the Green’s displacements and stresses is also introduced. These implementations have been incorporated into a boundary element method (BEM) code developed by the authors. The numerical results of Green’s displacements, stresses and stress derivatives are in perfect agreement with closed-form solutions for transversely isotropic solids. The BEM code results are also in very close agreement with both exact solutions and other BEM formulations, even if coarse mesh discretizations are used.