The stress model of the hybrid-Tre!tz "nite element formulation is applied to the linear elastostatic analysis of solids. The stresses are approximated in the domain of the element and displacements on its boundary. Complete, linearly independent, hierarchical polynomial approximation functions are used in both domain and boundary approximations. The displacement basis is de"ned independently on each inter-element surface. Continuity at the edges and on the corners of the elements is not enforced a priori. The stress basis is constrained to solve locally the Beltrami governing di!erential equation. It is derived from the associated Papkovitch}Neuber elastic displacement solution. Generalized variables are used to ensure that the approximations are independent of the geometric description of the elements. The solving system is derived directly from the fundamental relations of elastostatics. The solving system is symmetric, when the same property applies to the local elasticity condition, sparse, described by boundary integral arrays and well suited to p-re"nement and parallel processing. The numerical implementation of these equations is discussed and numerical tests are presented to illustrate the performance of the "nite element formulation. thus obtained can be recognized in the embracing works of Pian and Tong [3] and Brezzi and Fortin [4].The formulation used here was identi"ed later as the stress model of the hybrid-Tre!tz "nite element formulation and applied to di!erent two-dimensional structural analysis problems [5}8]. It is, therefore, directly related with the displacement frame Tre!tz element developed by Jirousek and Stein and their co-workers [9, 10]. The theoretical basis of the formulation is established in [11], where a comparative analysis with the related T-elements is also presented.In the terminology adopted here, a "nite element derived from the direct approximation of the stress "eld is termed a stress "nite element model. The stress "eld is directly related with two fundamental conditions, namely the equilibrium condition in the domain of the element and the di!usivity (Neumann) condition on the associated boundary tractions (Cauchy stresses). Three classes of formulations can be established depending on the equilibrium constraints enforced a priori on the stress basis, namely the hybrid-mixed, the hybrid and the hybrid-Tre!tz formulations of the "nite element stress model [12,13].The hybrid-mixed formulation is the most general. No constraints (besides linear independence and completeness) are placed on the stress approximation basis. Besides the stresses, the displacements are also approximated, and independently, in the domain and on the boundary of the element. These bases are used to enforce on average the equilibrium and the di!usivity conditions, respectively. The hybrid formulation is obtained by constraining a priori the stress approximation to satisfy locally the equilibrium condition. Consequently, the displacement approximation in the domain of the element becomes redundant an...