An efficient synthesis approach is developed which permits the calculation of the steady-state frequency response (FRF) of an assembly that is comprised of linear components and nonlinear joints. Receptance or accelerance matrices are used to characterize the linear components, which permits condensation of the system to just the joint degrees-of-freedom which is an immediate computational savings. Integrated into the technique is a frequency continuation scheme which permits studies across a wide range of frequencies — a key improvement to the continuation technique is the use of spline fits for the component FRF curves which provides significant cost savings over using modal expansion series at each iteration and also permits experimental data to be used if available. Finally, the calculated nonlinear joint forces are then used to study the power flow, structural and airborne, in the assembly, as well as to provide insight into the dynamic behavior within the components.