Modeling and Simulation of Monolithic Single-Supply Power Operational Ampliﬁers

: In this paper a simple PSpice (Personal Simulation Program with Integrated Circuit Emphasis) macro-model was developed, and veriﬁed for monolithic power operational ampliﬁers operated with a single-supply voltage. The proposed macro-model is developed using simpliﬁcation and build-up techniques for macro-modeling of operational ampliﬁers and simulates the basic static and dynamic characteristics, including input impedance, small-signal frequency responses at various voltage gains, output power versus supply voltage, slew-rate-limiting, voltage limiting, output offset voltage versus supply voltage ripples, and output resistance. Furthermore, the macro-model also takes into account the ground reference voltage in the ampliﬁer at a single power supply voltage. The model is implemented as a hierarchical structure suitable for the PSpice circuit simulation platform. The sub-circuit is built using standard PSpice components and analog behavioral modeling blocks. The accuracy of the model is veriﬁed by extracting the model parameters for single-supply power operational ampliﬁer TDA2005 from ST Microelectronics as example. The effectiveness of the model is validated by comparing the simulation results of the electrical parameters with the corresponding measured values obtained by experimental testing of sample circuits. The comparative analysis shows that the relative error of the modeled large-signal parameters is less than 15%. Moreover, an error of 15% is quite acceptable, considering the technological tolerances of the electrical parameters for this type of analog ICs.


Introduction
Despite the growing part of the digital devices in the modern electronic systems, there are still many applications that will always perform by analog circuits. One application of this type is when, for example, a small signal from a microphone or digital-to-analog converter needs to be amplified and to drive an outside load (e.g., actuator, loudspeaker, or wireline). In many cases, the amplifier has to transmit a significant electrical power (e.g., P L ≥ 1 W) to a certain load with high energy efficiency and without significant distortion. According to the value of the higher cutoff frequency f h , the power amplifiers are low-frequency (with f h ≤ 100 kHz) and high-frequency (with f h > 100 kHz) circuits.
The low-frequency monolithic power operational amplifiers are essential building blocks of the audio power amplifiers and the output drivers in many industrial controllers. To take full advantage of the monolithic powerful amplifiers and to provide the process of implementation of new electronic systems, a detailed analysis and evaluation is needed. In the process of circuit design, it is necessary to perform a preliminary evaluation of the various electronic circuit topologies in terms of static and dynamic modes of operation. For example, the achieving a certain operating frequency range at a certain value of the power gain is examined, as well as the ability to provide the output power on a load while maintaining the stability of the circuit. Moreover, the functionality of the designed electronic circuit must be simulated, the values of the passive components has to be adjusted and the obtained structure has to be experimentally tested to prove in advance the fulfillment of the requirements in the technical specification. Conducting accurate analog circuit simulations is conditioned by the simple and accurate models or macro-models of the passive and active components. The majority of the previously developed macromodels of operational amplifiers (op amps), suitable for the PSpice (Personal Simulation Program with Integrated Circuit Emphasis) circuit simulation program, have been intended for low power applications (ranging up to 100 mA) [1][2][3][4][5][6][7][8][9]. Moreover, the description of all models is adapted to the circuit simulators based on the SPICE (Simulation Program with Integrated Circuit Emphasis) program. In [2,6] a midpoint has been realized in the equivalent circuits, for relating all internal voltages to it, and to provide the proper operation of the macro-model at an asymmetrical power supply voltage. Also, in [7,9] through external resistor dividers, it has been possible to work with a single power supply voltage. However, these additional functional blocks do not provide a full possibility for the implementation of class-B and class-AB of the output stages in the power amplifiers, when only a single-supply voltage is used.
In this paper, based on the simulation models, represented in [1, 2,7], a simple electrical macro-model was developed, and verified for monolithic power amplifiers operated with a single-supply voltage. The developed macro-model included the basic static and dynamic effect of the real device and it was implemented as a hierarchical structure, suitable for the PSpice circuit simulation program. The paper is an extension of work [10] originally reported in the 28th International Conference on Mixed Design of Integrated Circuits and Systems-MIXDES 2021.
The paper is organized as follows. Section 2 presents the macro-model descriptionstructure of the commercially available monolithic power amplifiers, block diagram of the proposed macro-model and mathematical equations with the corresponding model parameters. Also, in Section 2 are represented the analytical expressions for determination of the modeling parameters of the proposed model. The verification of the developed macro-model is given in Section 3 and is performed by comparing the simulation results with experimental results at frequency and time domain for sample electronic circuits. For verification of the workability and efficiency of the macro-model, the model parameters are extracted for IC TDA2005 from STMicroelectronics as example. Finally, in Section 4 of the paper, the concluding remarks and guidelines for future work are given.

Macromodel Development
The basic technical requirements to the monolithic power amplifiers are: (1) to provide a certain electrical power to an external load, which in the general case has active and reactive components; (2) introduce minimal non-linear distortions into the wave-form of the amplified input signal; (3) to provide high energy efficiency (in this way, the greater economy is ensured, which is essential at high output power); (4) to have large input and small output resistance (this makes it easier to match the load with the output stage of the amplifier and achieves a high value of the power gain); (5) to have protection stages against short circuit or thermal overload of the output transistors. The implementation by the semiconductor manufacturers of the above requirements determines the efficiency and workability of the power amplifiers in a wide range of variations of the input voltage or power supply voltage.
Several semiconductor manufacturers, such as Maxim ® , Analog Devices ® , Philips ® , Sanyo ® , Samsung ® , SGS Thomson ® , ST Microelectronics ® , and Texas Instruments ® offer commercially available various types of monolithic power amplifiers. Some integrated circuits of the power amplifiers contain only an output stage, but often a preamplifier is included in the chip. Their basic electrical circuit is similar to those with discrete components [11][12][13][14][15]. The passive components that cannot be integrated into the chip (such as, capacitors with high capacitance, coils, and potentiometers) are connected externally to the circuit. One of the basic monolithic integrated power amplifiers is those of the ST Microelectronics from the TDA-series [14,15]. The basic DC and AC parameters and their typical values for the monolithic amplifiers from TDA-series are presented in Table 1.
In Figure 1 is a simplified circuit diagram of a TDA2030 power amplifier is given. As can be seen, the electrical circuit is divided into three sections (or stages): input stage, intermediate (or driver) stage and output stage. The input stage is a differential amplifier with a symmetrical input and single-ended output port. It provides the input impedance of the circuit and part of the open-loop voltage gain. The intermediate stage is an activeloaded common emitter amplifier, preceded by a voltage limiting circuit. From the output terminal of the intermediate stage, the signal is applied to the output stage. Its structure is similar to the basic output stage and includes two complementary Darlington transistors. The output stage provides the output voltage swing and output resistance. The transistors T 1 , and T 2 form the composite NPN-transistor according to the Darlington scheme, and the transistor T 3 , together with the Darlington transistors T 4 and T 5 , form a composite PNP-transistor. In addition to the main amplifier stages, the chip also includes thermal protection and current limiting circuits.
Based on the simplified circuit diagram in Figure 1, a block diagram of the proposed model ( Figure 2) is developed. In Figure 2, the block diagram includes basically an input stage, an intermediate (amplifying) stage, and an output stage connected in a cascade structure. Additionally, three blocks are connected to the structure, providing the midpoint of the model, the floating ground, and the supply voltage rejection ratio. According to the complexity of the amplitude-frequency characteristic of a concrete amplifier, the number of frequency-shaping stages can be increased, as well as the structure of the stage modeling the suppression of the pulsations in the supply voltage can be changed.
Based on the simplified circuit diagram in Figure 1, a block diagram of the proposed model ( Figure 2) is developed. In Figure 2, the block diagram includes basically an input stage, an intermediate (amplifying) stage, and an output stage connected in a cascade structure. Additionally, three blocks are connected to the structure, providing the midpoint of the model, the floating ground, and the supply voltage rejection ratio. According to the complexity of the amplitude-frequency characteristic of a concrete amplifier, the number of frequency-shaping stages can be increased, as well as the structure of the stage modeling the suppression of the pulsations in the supply voltage can be changed.      [14,15].
Based on the simplified circuit diagram in Figure 1, a block diagram of the proposed model ( Figure 2) is developed. In Figure 2, the block diagram includes basically an input stage, an intermediate (amplifying) stage, and an output stage connected in a cascade structure. Additionally, three blocks are connected to the structure, providing the midpoint of the model, the floating ground, and the supply voltage rejection ratio. According to the complexity of the amplitude-frequency characteristic of a concrete amplifier, the number of frequency-shaping stages can be increased, as well as the structure of the stage modeling the suppression of the pulsations in the supply voltage can be changed.

PSpice Implementation
To develop PSpice macro-model for monolithic power amplifiers, simplification and build-up techniques for macro-modeling of operational amplifiers was used [1]. Based on the block diagram, shown in Figure 2, the macro-model was implemented as a two-level hierarchical structure. The higher-level block diagram with the modeling parameters, defined without concrete values, is presented in Figure 3. Depending on the choice of a certain type of amplifier in dialog mode, concrete numerical values can be set for the modeling parameters. The equivalent electrical circuit (the lower level of hierarchy) of the macro-model is given in Figure 4. As can be seen, there is a correspondence between the names of the external terminals (or pins), and also in the description of the block diagram ( Figure 3) the name of the equivalent circuit in Figure 4 was set, so that it can be loaded and processed during of the simulation process.

PSpice Implementation
To develop PSpice macro-model for monolithic power amplifiers, simplification and build-up techniques for macro-modeling of operational amplifiers was used [1]. Based on the block diagram, shown in Figure 2, the macro-model was implemented as a two-level hierarchical structure. The higher-level block diagram with the modeling parameters, defined without concrete values, is presented in Figure 3. Depending on the choice of a certain type of amplifier in dialog mode, concrete numerical values can be set for the modeling parameters. The equivalent electrical circuit (the lower level of hierarchy) of the macro-model is given in Figure 4. As can be seen, there is a correspondence between the names of the external terminals (or pins), and also in the description of the block diagram ( Figure 3) the name of the equivalent circuit in Figure 4 was set, so that it can be loaded and processed during of the simulation process.  In particular, the structure of the equivalent circuit ( Figure 4) includes both standard PSpice components and analog behavioral modeling (ABM) blocks. The usage of ABM blocks, such as various controlled sources, provides the implementation of linear and nonlinear mathematical functions in the PSpice program. As can be seen, the equivalent electrical circuit is represented as properly connected controlled sources, independently controlled sources, resistors, capacitors, and a small group of ideal diodes. The external interface terminals of the circuit are non-inverting input (input+), inverting input (input-), an output terminal (out), supply voltage ripple rejection terminal (SVRR), positive power supply (+Vs), and analog ground (GND). As can be seen, there is no ground reference in any of the signal-processing blocks. The only exception is the functional block, which simulates the attenuation of the power supply voltage ripples. All internally generated node voltages are referred to as the voltage at node nVgnd, which is the floating ground of the equivalent circuit. The node voltage nVgnd is produced by the VCVS EH controlled by the midpoint and power supply voltage. The midpoint, called nVref in the model, was

PSpice Implementation
To develop PSpice macro-model for monolithic power amplifiers, simplification and build-up techniques for macro-modeling of operational amplifiers was used [1]. Based on the block diagram, shown in Figure 2, the macro-model was implemented as a two-level hierarchical structure. The higher-level block diagram with the modeling parameters, defined without concrete values, is presented in Figure 3. Depending on the choice of a certain type of amplifier in dialog mode, concrete numerical values can be set for the modeling parameters. The equivalent electrical circuit (the lower level of hierarchy) of the macro-model is given in Figure 4. As can be seen, there is a correspondence between the names of the external terminals (or pins), and also in the description of the block diagram ( Figure 3) the name of the equivalent circuit in Figure 4 was set, so that it can be loaded and processed during of the simulation process.  In particular, the structure of the equivalent circuit ( Figure 4) includes both standard PSpice components and analog behavioral modeling (ABM) blocks. The usage of ABM blocks, such as various controlled sources, provides the implementation of linear and nonlinear mathematical functions in the PSpice program. As can be seen, the equivalent electrical circuit is represented as properly connected controlled sources, independently controlled sources, resistors, capacitors, and a small group of ideal diodes. The external interface terminals of the circuit are non-inverting input (input+), inverting input (input-), an output terminal (out), supply voltage ripple rejection terminal (SVRR), positive power supply (+Vs), and analog ground (GND). As can be seen, there is no ground reference in any of the signal-processing blocks. The only exception is the functional block, which simulates the attenuation of the power supply voltage ripples. All internally generated node voltages are referred to as the voltage at node nVgnd, which is the floating ground of the equivalent circuit. The node voltage nVgnd is produced by the VCVS EH controlled by the midpoint and power supply voltage. The midpoint, called nVref in the model, was In particular, the structure of the equivalent circuit ( Figure 4) includes both standard PSpice components and analog behavioral modeling (ABM) blocks. The usage of ABM blocks, such as various controlled sources, provides the implementation of linear and non-linear mathematical functions in the PSpice program. As can be seen, the equivalent electrical circuit is represented as properly connected controlled sources, independently controlled sources, resistors, capacitors, and a small group of ideal diodes. The external interface terminals of the circuit are non-inverting input (input+), inverting input (input-), an output terminal (out), supply voltage ripple rejection terminal (SVRR), positive power supply (+Vs), and analog ground (GND). As can be seen, there is no ground reference in any of the signal-processing blocks. The only exception is the functional block, which simulates the attenuation of the power supply voltage ripples. All internally generated node voltages are referred to as the voltage at node nVgnd, which is the floating ground of the equivalent circuit. The node voltage nVgnd is produced by the VCVS E H controlled by the midpoint and power supply voltage. The midpoint, called nVref in the model, was generated by two equal resistors (R r1 and R r2 ) connected between the supply rail and the external ground reference.
The input stage consists of r id and C id , modeling the input resistance and capacitance between the two input terminals (1 and 7) at a small signal. r + iCM and r − iCM are the input resistances, reflecting the behavior at a common-mode input signal. The input resistance r id is defined as the resistance measured directly between the inverting and non-inverting input of the amplifier: . In the case of an amplifier circuit with asymmetric input, a circuit of a non-inverting amplifier with a grounded inverting input and an active non-inverting input are possible, and vice versa-a circuit of an inverting amplifier with a grounded non-inverting input.
The input resistance in these cases is represented by The input common-mode resistance is represented, according to For the majority of the amplifiers r + iCM ≈ r − iCM = 2r iCM . The input bias current was modeled by two ideal current sources (I + iB and I − iB ), connected between the two input terminals (1 and 7) and midpoint nVref. The input offset voltage and offset drift as a result of the power supply ripples is represented with a voltagecontrolled voltage source (VCVS) E os , connected between nodes 2 and 7. The offset voltage modelling power supply ripples (PSRR-power-supply rejection ratio), which are accurately referred to the VCVS E os , come from two separate frequency shaping stages in the model. The first additionally defined stage consisted of a capacitor C 0SVR and a grounded resistor R 0SVR . The C 0SVR is connected between the power supply voltage (node Vsplus) and the input of the second shaping stage. The voltage drop over the R 0SVR presents the voltage at a node 6 that depends on the ripples of the supply voltage (V S ± ∆V S ). The transfer function of the stage is expressed by The second stage consisted of a voltage-controlled current source (VCCS) G SVR controlled by the power supply voltage, and the RLC-stage. The current source was chosen to be linear one-port generator, following the equation The current I G SVR will flow through the elements R 1SVR , L 1SVR , and C SVR = C 1SVR + C 3 (where C 3 is an external grounded capacitor, connected to node SVR), towards the internal ground. In such a way, the voltage V(SVR,0) will depend on the amplitude and the frequency of the ripples of the supply line, and the transfer function of the stage can be expressed by where k is DC proportionality constant, T 1 is the time-constant of the polynomic in the numerator, ξ is the damping factor, and T 2 is the time-constant of the polynomic in the denominator.
From the comparison of the left and right sides of the above equation, the following formulas were obtained for the main parameters: k = k 1,GSVR R 1SVR , T 1 = √ L 1SVR C SVRtime constant determining the minimum in the frequency response, ξ = (2R 1SVR ) −1 √ L 1SVR /C SVR and T 2 = R 1SVR C SVR , where L 1SVR is equivalent series inductance, C SVR is equivalent series capacitance, and R 1SVR is equivalent parallel resistance. The proposed stage from Figure 4 provides a serial resonance, determined by the inductance L 1SVR , and the capacitance C SVR , with a minimum value of the resonance resistance. The voltage generated at node SVRR is used for the equation of VCVS E os as follows V Eos = k 0,Eos + k 1,Eos V(SVRR, 0) The component k 0,Eos in Equation (6) determines the input offset voltage constant (V os ), the coefficient k 1,Eoos presents the change in V os as a result of changing the supply voltage. Then, PSRR will be determined by the following expression where GAI N = V OUT /V(1, 7) is the differential voltage gain of the amplifier and GAI N PSR where K 1 [V −1 ] is the multiplication factor of the stage. The transfer function (9) is represented in graphical form in Figure 5. Approximately linear region of this function is −0.5K −1 Energies 2021, 14, x FOR PEER REVIEW 8 of 18 , elements can be taken only up to the third order term in the Taylor series, wherein the formula (9) yields . The maximum output current 1 I charged the ca-  The simplification of expression (9) can be performed if the function is decomposed in the Taylor order, i.e., at |x| < π/2. If V(1, 7) << K 1 −1 , elements can be taken only up to the third order term in the Taylor series, wherein the formula (9) yields If V(1, 7) > K 1 −1 , the V(3) ≈ I 1 R p1 . The maximum output current I 1 charged the capacitor C p1 , which created a change of the output voltage per unit time of the amplifier. This rate of change of the voltage per unit time defined the maximum dV/dt (voltage change rate) or slew rate (or SR) of the amplifier: For the model's small-signal open-loop DC voltage gain, is obtained To obtain minimal nonlinear distortions in the waveform of the output voltage , it is necessary that the input voltage V(1, 7) is lower than 0.5K −1 1 . Since the mathematical function tanh is symmetric with respect to the origin of the coordinate system, and with respect to the characteristics of the differential pair of transistors in the input stage (transistors have to be as close as possible to technological parameters, which is achieved by applying a single technological process of obtaining the electrical circuit), for the developed macro-model it is possible to obtain smaller values of the nonlinear distortions.
The voltage-limiting network consisted of two ideal diodes (D 1 and D 2 ) and two voltage sources (V 1 and V 2 ), connected in series. This network limited the voltage at node 3 and the other internal voltages of the macro-model up to values below the power supply voltage during an input-overdrive condition: where ϕ T is the thermal voltage and I S is the saturation current for the two diodes. The first (dominant) and the second pole frequency are represented in the model by two frequency shaping stages. The dominant pole was modeled by the VCCS G 1 , the resistor R p1 and the capacitor C p1 . The second (non-dominant) pole was modeled by the VCCS G 2 , the resistor R p2 and the capacitor C p2 . The DC voltage gain of the second frequency-shaping block was equal to unity, because the g m2 of the VCCS G 2 is equal to the reciprocal of the resistance R p2 , connected from each node of the VCCS to the floating node of the circuit. In a wide frequency range, the open-loop amplitude-frequency response of the model can be represented by where f p1 = 1/2πR p1 C p1 is the first pole frequency and f p2 = 1/2πR p2 C p2 is the secondary pole frequency.
The output stage of the model consisted of controlled sources and passive components provided the output signal. It simulated the output resistance and the DC and dynamic of the current consumed by the external power supply device. The VCCS G o in the output block drives the resistor R o , connected between the node 5 and the floating ground. The G o acted as an active current generator and provided the desired voltage drop over its parallel resistor. The DC voltage gain of the output block was equal to unity, because the g mo of the VCCS G o was equal to the reciprocal of the resistance R o .
The limitation of the output current was presented in the macro-model through the controlled voltage source E lim , the resistor R sense , and the pair of diodes D L1 and D L2 [16]. When the voltage drop across the resistor R sense becomes approximately equal to 0.6 V, one of the two diodes (D L1 or D L2 ) turns on and the current is limited to I omax = 0.6V/R sense .
When there is no load at the output, the model draws only a quiescent current from the power-supply rail by the independent current source I CC in the model, thus behaves similarly to an ideal class-B output stage with a gain equal to unity. When an AC voltage is applied to the input port, as a result of the amplification, the output voltage will Energies 2021, 14, 4611 9 of 17 have greater amplitude, which is obtained by increasing the power consumption. In the proposed model, the dynamic behavior of the supply current is produced through a pair of current-controlled current sources (CCCSs) (F 1 and F 2 ) with unity gain and a pair of ideal diodes (D 3 and D 4 ), similar to the op amp models, presented in [2,3,6]. The F 1 is controlled by the magnitude of the output current through the load, which in turn generates a current, driving the two diodes. Depending on the magnitude and the direction of the current through the F 2 , an AC current flows from the power supply to the ground, which average value I CC,av is proportional to the real consumed current. Another technique for modeling the supply current variations is by using relatively complex high-order (n > 2) nonlinear current source F 2 [4].

Parameter Extraction
Based on the theoretical analysis of the proposed equivalent circuit in this subsection are proposed analytical expressions for determining the modeling parameters of the proposed model (Figure 4) for simulations of monolithic power amplifiers. For better clarity, the expressions are systematized in a table (see Table 2) and they are given separately for each of the individual building blocks. When determining the values of the model parameters, some of them are obtained directly from the datasheet typical values, while others are typical estimated values from device design. Table 2. Analytical expressions for the modeling parameters.

Input Stage
(1) The parameters r id , C id , r + iCM ≈ r − iCM = 2r iCM , I iB , k 0,Eos = V io and I omax of the input stage is obtained from the datasheet and are substituted directly into the macro-model.
(2) Power-supply rejection stage: R 0SVR = 1.10 6 Ω, where f min is the lower frequency with a minimum value of the PSRR. The damping factor is ξ = 1/2Q (where Q is the quality factor), T 1 = 1/2π f max ( f max is the frequency, at which the PSRR has maximum value), C 1SVR = 10pF-parasitic terminal capacitance, R 1SVR = (2ξT 1 )/C 3 (where C 3 is the value from the datasheet), L 1SVR = T 2 1 /C 1SVR -series inductor, T 2 = R 1SVR C 1SVR and k 1,SVR = 1/PSRR max R 1SVR , where PSRR max is the maximum value of the PSRR at f max .

Intermediate (Second) Stage
(3) The current I 1 according to the datasheet is set to be equal to 1 mA, 0.1 mA, or 0.01 mA. The value of the I 1 has to be less than the quiescent current. (4) Voltage gain stage and dominant pole frequency shaping stage where f p1 ≈ f 1 /A d0 is the dominant pole frequency, f 1 is a unity-gain bandwidth and K 1 = A d0 /(I 1 R p1 ); (5) Non-dominant pole shaping stage R p2 = 1.10 6 Ω, C p2 = 1/2π f p2 R p2 and g m2 = 1/R p2 is the transconductance of the VCCS G 2 ; (6) Voltage limiting stage The values of the voltage sources in a voltage-limiting network are: V 1 = +V s − V OUT,max + ϕ T ln(I 1 /I S ) and V 2 = −V OUT,min + ϕ T ln(I 1 /I S ), where ϕ T = 0.02585 V at T = 27 • C and I S = 1.10 −12 A for both diodes.

(7) Output resistance and voltage gain of the output stage
The R o = r OUT is the output resistance from the datasheet; g mo = 1/R o is the transconductance of the output VCCS G o ; (8) Current limiting stage and quiescent current For the maximum output current the sense resistor is R sense = 0.6 V/I omax ; The total quiescent current that flows from the power supply rail to the ground in the model is

Results and Discussion
To verify the efficiency of the proposed macro-model the simulation results of the electrical parameters were compared with the corresponding results obtained by the experimental testing of sample analog circuit, using commercially available monolithic power amplifier TDA2005 [17]. The modeling parameters were extracted for TDA2005 by analyzing semiconductor data books and through some characterization measurements. All the modeling parameters for the proposed model that are in good agreement with the experimental results related to the TDA2005 model at V + S = 14.4 V, R L = 4 Ω and T A = 25 • C, are: k 0,Eos = 2.5 µV, r id = 70 kΩ, r + iCM = r − iCM = 140 kΩ, C id = 5 pF, I + iB = I − iB = 100 nA, PSRR max = 562.34 (or 55 dB at f max = 400 Hz), C 0SVR = 16 nF (at f min = 10 Hz), C 1SVR = 10 pF, L 1SVR = 15.84 mH, R 1SVR = 79.6 Ω, I 1 = 1 mA, To obtain the amplitude-frequency responses for the various voltage gains, the test circuit, shown in Figure 6 was used. The circuit for the simulation with OrCAD PSpice ® 9.0 (previously MSIM ® 8) was based on the test condition, given in the datasheet. For the experimental testing, the same circuit was implemented on an FR4 printed circuit board laminate with through-hole device passive components. The recommended area of aluminum plate heatsink was S HS ≥ 400 cm 2 [15], where the junction-to-case thermal resistance is R th j−c = 3 • C/W [17].
The total quiescent current that flows from the power supply rail to the ground in the model is

Results and Discussion
To verify the efficiency of the proposed macro-model the simulation results of the electrical parameters were compared with the corresponding results obtained by the experimental testing of sample analog circuit, using commercially available monolithic power amplifier TDA2005 [17]. The modeling parameters were extracted for TDA2005 by analyzing semiconductor data books and through some characterization measurements. All the modeling parameters for the proposed model that are in good agreement with the experimental results related to the TDA2005 model at , and mA I CC 75  . To obtain the amplitude-frequency responses for the various voltage gains, the test circuit, shown in Figure 6 was used. The circuit for the simulation with OrCAD PSpice ® 9.0 (previously MSIM ® 8) was based on the test condition, given in the datasheet. For the experimental testing, the same circuit was implemented on an FR4 printed circuit board laminate with through-hole device passive components. The recommended area of aluminum plate heatsink was 2 400cm S HS  [15], where the junction-to-case thermal resistance [17]. Figure 7  . The maximum relative error between the simulation results of the cutoff frequencies and the measured values was not greater than 10%. Figure 6. Test circuit [17] for simulation and experimental testing at various voltage gains A V , through the feedback resistor R 2 .   Figure 8 illustrates the frequency response for the simulated PSRR at VS = 14.4 V and RL = 4 Ω. The comparative analysis showed that the simulated response fitted quite closely the actual response [17] (the relative error was below 8%) with correct amount of maximal value of PSRR and the slope of characteristic in the transition area.
The workability of the proposed macro-model was validated by examining the transient characteristics at a large-signal using the electronic circuits with various external load resistance. For a large-signal, the simulation and experimental results of the energy parameters at RL = 2 Ω ± 5%, and 4 Ω ± 5%, are plotted in Figures 9 and 10. To obtain the electrical parameters, the test circuit with R2 = 1.2 kΩ, given in Figure 6, was used. The sinusoidal input signal with a frequency of 1 kHz was obtained by SFG-2010. The amplitude of the input signal had an amplitude, at which the output signal had no visible nonlinear distortions (below 5%). The supply voltage was swept from 8 V to 18 V by using HY3005D-3. The achieved relative error between the measured values for the output power Additionally, at a large signal in Figure 11a,b, the results of a simulation study of the nonlinear distortions at Vs  14.4 V and 63 ) are shown. For this purpose, a transient analysis in combination with a Fourier analysis was performed in the PSpice program environment. An input signal with a sinusoidal form was applied to the input terminals of the circuit to provide approximately the output power for which a value was given for the nonlinear distortions in the datasheet of the real device. For the frequency range from 50 Hz to 15 kHz, a final simulation time equal to 100 ms was selected. In this case, the print step and the step celling had the same values equal to 1 μs. In this way, a waveform of the output signal can be provided with a print step of no more  than PERIOD OF THE INPUT SIGNAL/100. As can be seen, the relative error, for some values of the input frequency, was greater than 20% (as the maximum value of the total harmonic distortion (THD) was below 1%), therefore the results showed that the proposed macro-model can simulate the trend of change of nonlinear distortions in the process of preliminary studies of sample electronic circuits. The average value of the absolute deviation (taking into account the algebraic sign for each value of the coefficient) at a load equal to 2 Ω is 13.069%, and at a load equal to 4 Ω is 20.821%. In many practical cases, the modeled direction of change of nonlinear distortions is sufficient, taking into account the technological tolerances of the parameters of the integral and discrete components, as well as the instability of the supply voltage. For the maximum output power with a typical value of 20 W, when bridging two power operational amplifiers, according to the test conditions in the datasheet of the real amplifier, as a result of the simulation study for each of the output voltages a THD was equal to 3.9%. According to the datasheet, the THD was up to 10%. The workability of the proposed macro-model was validated by examining the transient characteristics at a large-signal using the electronic circuits with various external load resistance. For a large-signal, the simulation and experimental results of the energy parameters at R L = 2 Ω ± 5%, and 4 Ω ± 5%, are plotted in Figures 9 and 10. To obtain the electrical parameters, the test circuit with R 2 = 1.2 kΩ, given in Figure 6, was used. The sinusoidal input signal with a frequency of 1 kHz was obtained by SFG-2010. The amplitude of the input signal had an amplitude, at which the output signal had no visible non-linear distortions (below 5%). The supply voltage was swept from 8 V to 18 V by using HY3005D-3. The achieved relative error between the measured values for the output power P L,max = V 2 o,max /2R L and the simulations for the proposed model, was up to 10%. The error for the energy efficiency η = P L,max /P CC,max (P CC,max = I CC,av × V + S , where I CC,av is the average value of the supply current) was below 10% (see Figure 10a,b) at V s > 12 V and R L = 2 Ω and 4 Ω. For V s ≤ 12 V, and R L = 4 Ω, Figure 10b the supply current had lower value in comparison with the I CC,av at R L = 2 Ω (Figure 10a), and the error was up to 15% due to the relatively greater influence of the nonlinearities of the volt-ampere characteristics of the transistors.   Error, % P L , W

Simulated
Error at R L = 4 Error at R L = 2 Measured  Error, % Error at R L = 2 Measured Error, % Error at R L = 4 Measured Additionally, at a large signal in Figure 11a,b, the results of a simulation study of the nonlinear distortions at V s ≤ 14.4 V and A V = 364.63 (≈ 50 dB) are shown. For this purpose, a transient analysis in combination with a Fourier analysis was performed in the PSpice program environment. An input signal with a sinusoidal form was applied to the input terminals of the circuit to provide approximately the output power for which a value was given for the nonlinear distortions in the datasheet of the real device. For the frequency range from 50 Hz to 15 kHz, a final simulation time equal to 100 ms was selected. In this case, the print step and the step celling had the same values equal to 1 µs.
In this way, a waveform of the output signal can be provided with a print step of no more than PERIOD OF THE INPUT SIGNAL/100. As can be seen, the relative error, for some values of the input frequency, was greater than 20% (as the maximum value of the total harmonic distortion (THD) was below 1%), therefore the results showed that the proposed macro-model can simulate the trend of change of nonlinear distortions in the process of preliminary studies of sample electronic circuits. The average value of the absolute deviation (taking into account the algebraic sign for each value of the coefficient) at a load equal to 2 Ω is 13.069%, and at a load equal to 4 Ω is 20.821%. In many practical cases, the modeled direction of change of nonlinear distortions is sufficient, taking into account the technological tolerances of the parameters of the integral and discrete components, as well as the instability of the supply voltage. For the maximum output power with a typical value of 20 W, when bridging two power operational amplifiers, according to the test conditions in the datasheet of the real amplifier, as a result of the simulation study for each of the output voltages a THD was equal to 3.9%. According to the datasheet, the THD was up to 10%. To obtain large-signal transient responses for the various load resistances the test circuit, given in Figure 12 is used. For this purpose, a signal with a pulse form of the voltage and amplitude at which a maximum swing of the output voltage is obtained was applied to the input terminals of the circuit. For the purpose of the simulation study, a transient analysis with a final time of 100 ms (frequency of the input signal 1 kHz and duty cycle 0.5) was set, and the print step and the step celling had the same values equal to 0.1 μs. For the transient response, the simulation and the experimental results of the output signal at RL = 4 Ω ±  and 8 Ω ± , are plotted in Figures 13 and 14. The output signal, obtained from the simulation study, was reported from the beginning of the time moment equal to 50 ms, because due to the influence of the coupling capacitors in the circuit there was a transient process of obtaining the steady-state value of the output voltage's amplitude. The measured values were Figure 13a) and (Figure 14a), respectively and the corresponding simulated value for the slew rates were Figure 13b) and (Figure 14b), respectively. The simulated results for the positive slew rate at both external load resistances are given in Figures 13c and 14c, respectively. The relative error between the simulation results of the slew rates and the measured values was below 5%. The voltage drop ( ) in the flat region of the output pulse was due to the influence of the coupling capacitors in the input and output network of the circuit, and the value of the load. As can be seen, for a load resistance of   8 , the simulated value was 2.72% ( Figure 14d) and the measured value was 2.25%. The Error, % Error at R L = 2 Measured Error, % Error at R L = 4 Measured To obtain large-signal transient responses for the various load resistances the test circuit, given in Figure 12 is used. For this purpose, a signal with a pulse form of the voltage and amplitude at which a maximum swing of the output voltage is obtained was applied to the input terminals of the circuit. For the purpose of the simulation study, a transient analysis with a final time of 100 ms (frequency of the input signal 1 kHz and duty cycle 0.5) was set, and the print step and the step celling had the same values equal to 0.1 µs. For the transient response, the simulation and the experimental results of the output signal at R L = 4 Ω ± 5%, and 8 Ω ± 5%, are plotted in Figures 13 and 14. The output signal, obtained from the simulation study, was reported from the beginning of the time moment equal to 50 ms, because due to the influence of the coupling capacitors in the circuit there was a transient process of obtaining the steady-state value of the output voltage's amplitude. The measured values were SR = 10.2 V/µs-R L = 4 Ω (Figure 13a) and SR = 10.6 V/µs-R L = 8 Ω (Figure 14a), respectively and the corresponding simulated value for the slew rates were SR = 9.98 V/µs at R L = 4 Ω (Figure 13b) and SR = 10.1 V/µs at R L = 8 Ω (Figure 14b), respectively. The simulated results for the positive slew rate at both external load resistances are given in Figures 13c and 14c, respectively. The relative error between the simulation results of the slew rates and the measured values was below 5%. The voltage drop (δ = ∆U 0 /U 0max ) in the flat region of the output pulse was due to the influence of the coupling capacitors in the input and output network of the circuit, and the value of the load. As can be seen, for a load resistance of R L = 8 Ω, the voltage drop had a smaller value. At R L = 4 Ω the simulated value of the voltage drop was equal to 5.41% (Figure 13d) and the measured value was 4.95%. For R L = 8 Ω, the simulated value was 2.72% ( Figure 14d) and the measured value was 2.25%. The maximum value of the relative error was up to 10%, that is acceptable comparing with the manufacturing tolerances, especially for the commercially available coupling capacitors.      Table 3 presents a comparison of the created model with known models of various types of operational amplifiers. As can be seen, all models and macro-models simulate the basic DC and AC parameters, for small signal operation and in most of them at dual power supply voltage from symmetrical external voltage source. When working with high  Table 3 presents a comparison of the created model with known models of various types of operational amplifiers. As can be seen, all models and macro-models simulate the basic DC and AC parameters, for small signal operation and in most of them at dual power supply voltage from symmetrical external voltage source. When working with high output current and voltage, and single-supply voltage, for some of the macro-models [7,9] are added external stages and function blocks allowing modeling of nonlinear distortions and the floating ground at the single-supply voltage without simulating simultaneous effects at large-signal on the operational amplifiers.

Conclusions
An electrical macro-model for monolithic power amplifiers was developed, implemented, and verified. The proposed model is intended for the PSpice-based simulation platform. The proposed model showed a good agreement between the simulation and the experimental testing, and the relative error for the basic static and dynamic small-signal parameters was below 10%. The maximum value of the relative error for the energy efficiency at a large signal was up to 15% (for supply voltages between 8 and 12 V) and below 10% at V s > 12 V, which guarantees a sufficient degree of accuracy, considering the technological tolerances of the electrical parameters of the components and some variations of the power supply voltage. Further, the also takes into account the internal floating ground with midpoint, offset drift as a result of the power supply ripples, dynamic behavior of the supply current, nonlinear distortion, and internal voltage-and current-limiting networks. This makes the developed model especially useful for analog single-supply processing circuit design at a block level, as thereby to use for predicting and optimize the values of the passive components for the performance of the functional blocks; for instance, to define and analyze the energy parameters, according to the power dissipation.
One possible direction for future work of the author is to develop a behavioral model by using VHDL-AMS or Verilog-A languages. As a result, the created models will allow describing the structure of the model with a smaller number of passive elements and con-trolled sources, as well as a smaller number of nodes. For example, one possible advantage of using VHDL-AMS is that the complete set of equations, that represent the behavior of the real device, can be implemented by combining behavioral elements and structural description as a net-list. By using this technique, the obtained behavioral models can have better convergence during the computation process, especially to simulate the transient responses at a large signal. Data Availability Statement: Some or all data, model generated, or used during the study are available in a repository or online.