A higher order, sub-parametric, laminated, 3D solid finite element was used for the analysis of very thick laminated composite plates. The geometry of this element is defined by four nodes in the X-Y plane which define a prism of material through the thickness of the laminate. There are twenty-four degrees of freedom at each node; translations at the upper and lower surfaces of the laminate in each of the three coordinate directions, and the derivatives of these translations with respect to each coordinate. This choice of degrees of freedom leads to displacement and strain compatibility at the corners. Stacking sequence effects are accounted for by explicitly integrating the strain energy density through the thickness of the element. The laminated solid element was combined with a gap-contact element to analyze thick laminated composite lugs loaded through flexible pins. The resulting model accounts for pin bending effects that produce non-uniform bearing stresses through the thickness of the lug. A thick composite lug experimental test program was performed, and provided data that was used to validate the analytical model. Two lug geometries and three stacking sequences were tested.
