Modeling procedures for the stress analysis of riveted lap joints in aircraft structure are proposed in this paper using both analytical and numerical methods. In the analytical method, a complex variational approach is employed to determine the stresses in joined plates containing single or multiple loaded holes. The effects of finite geometry are taken into account by the variational formulations. An iterative scheme is carried out to handle the deformation compatibility between all joined members, hi the numerical method, finite element analyses are conducted using the commercial packages MSC/Patran and MSC/Nastran. Gap elements are used to simulate the rivet-hole interactions. Both linear and nonlinear deformations are considered. The two modeling procedures are complementary to each other. While the analytical method is efficient for parametric studies, the numerical method is capable of dealing with relatively more complicated geometry and loading conditions. Results of several example problems are presented and reasonably good agreement between the two modeling procedures is demonstrated.