Linear, quadratic, and cubic isoparametric hexahedral solid elements have been added to the element library of NASTRAN. These elements are available for static, dynamic, buckling, and heat-transfer analyses. Because the isoparametric element matrices are generated by direct numerical integration over the volume of the element, variations in material properties, temperatures, and stresses within the elements are represented in the computations. In order to compare the accuracy of the new elements, three similar models of a slender cantilever were developed, one for each element. All elements performed well. As expected, however, the linear element model yielded excellent results only when shear behavior predominated. In contrast, the results obtained from the quadratic and cubic element models were excellent in both shear and bending