A boundary integral Trefftz formulation with symmetric collocation

The Galerkin and collocation methods are combined in the implementation of a boundary integral formulation based on the Trefftz method for linear elastostatics. A finite element approach is used in the derivation of the formulation. The domain is subdivided in regions or elements, which need not be bounded, simply connected or convex. The stress field is directly approximated in each element using a complete solution set of the governing Beltrami condition. This stress basis is used to enforce on average, in the Galerkin sense, the compatibility and elasticity conditions. The boundary of each element is, in turn, subdivided into boundary elements whereon the displacements are independently approximated using Dirac functions. This basis is used to enforce by collocation the static admissibility conditions, which reduce to the Neumann conditions as the stress approximation satisfies locally the domain equilibrium condition. The resulting solving system is symmetric and sparse. The coefficients of the structural matrices and vectors are defined either by regular boundary integral expressions or determined by direct collocation of the trial functions.