Procedures for modeling the performance of a transmitting array require either an explicit or implicit definition of the velocity of each radiating element as a function of the electrical input signal and the acoustic pressure on the radiating surface. To date, most array modeling work has been restricted to linear systems where the velocities are computed from the electrical input signal and the linear equivalent circuit for each element in the array. When an array contains nonlinear amplifiers or transducers, this equivalent circuit, in general, will be a function both of the input signal and the pressure on the radiating surface, and the linear analysis techniques cannot be used. Two methods, repeated matrix inversion and iterative lookup, have been developed and applied to the hydroacoustic amplifier/transducer module. In both methods, the nonlinear operating profile characteristics of the coupled amplifier/transducer module are known, either by prediction using a simulation program for nonlinear devices or by measurement. The operating profile is defined by several equivalent Thevenin circuits in the former method and by tabular data in the latter. Several examples show that for larger arrays the latter method requires much less computer time and storage than the matrix method.