This paper reports theoretical analysis and experimental results on the dynamics of piezoelectric MEMS mechanical logic relays. The multiple degree of freedom analytical model, based on modal decomposition, utilizes modal parameters obtained from finite element analysis and an analytical model of piezoelectric actuation. The model accounts for exact device geometry, damping, drive waveform variables, and high electric field piezoelectric nonlinearity. The piezoelectrically excited modal force is calculated directly and provides insight into design optimization for switching speed. The model accurately predicts the propagation delay dependence on actuation voltage of mechanically distinct relay designs. The model explains the observed discrepancies in switching speed of these devices relative to single degree of freedom switching speed models and suggests the strong potential for improved switching speed performance in relays designed for mechanical logic and RF circuits through the exploitation of higher order vibrational modes.