A piezoelectric rotary ultrasonic motor is modeled for the purpose of predicting, a priori, motor performance as a function of design parameters. The Rayleigh-Ritz assumed mode energy method is used to model the distributed piezoceramics and the traveling-wave dynamics of the stator. Natural frequencies and dynamic modeshapes are obtained for a generally configured motor. Nonlinear normal and tangential interface forces between the rotor and stator are incorporated into the forcing function along with the linear piezoelectric forcing. Given the applied torque, applied axial loading, and piezo drive voltages as inputs to the model, general motor performance measures are obtained--namely speed, input power, output power, and efficiency. The approach presented here provides a general framework for modeling these motors as well as a design tool for optimizing prototypes with the added flexibility of allowing for a wide variety of geometries and materials.