Abstract
A study regarding to a development of the versatile numerical analysis program applicable to general 3-dimensional elastic problems based on the Body Force Method was examined. In the present analysis, a planner triangular surface element is chosen as the boundary element. In order to overcome the difficulty arising in a boundary integral of the fundamental solution with singular functions, a polar transform scheme is adopted. Through numerical examples, it was found that the present method is practical and convenient to analyze a stress concentration problem gregarding to a spherical cavity in an infinite solid.