Multimass modeling of phonation should be based on valid aerodynamic characteristics. Measured glottal wall pressures over a wide range of glottal geometries and transglottal pressures were used to obtain average driving forces on the surfaces of two-mass models of phonation. A Plexiglas model of the larynx was used for this purpose. Predictions from the classic equations of Ishizaka and Flanagan were compared to these driving forces and movies of the calculated motions of the two-mass models were produced. The data show that the magnitudes of the negative pressures on the lower glottal region tend to be considerably less than those predicted by the Ishizaka–Flanagan equations, which leads to an increase in the number of cycles (longer time) for the vocal folds to reach steady state, lower peak volume flows, and larger open quotients. Subsequently, analytic expressions were used to represent the intraglottal pressures, with special attention to the entrance losses and the form used to describe the Poisseuille frictional effects, and their accuracy was assessed by a comparison with the measured flow rates and pressures. Results from the analytic representations were compared with those generated from the Ishizaka and Flanagan formulation. [Work supported by NIH Grant No. 2 R01 DC03577.]