Piezoelectric Transducers: Complete Electromechanical Model with Parameter Extraction

This paper presents a complete electromechanical (EM) model of piezoelectric transducers (PTs) independent of high or low coupling assumptions, vibration conditions, and geometry. The PT’s spring stiffness is modeled as part of the domain coupling transformer, and the piezoelectric EM coupling coefficient is modeled explicitly as a split inductor transformer. This separates the coupling coefficient from the coefficient used for conversion between mechanical and electrical domains, providing a more insightful understanding of the energy transfers occurring within a PT and allowing for analysis not previously possible. This also illustrates the role the PT’s spring plays in EM energy conversion. The model is analyzed and discussed from a circuits and energy harvesting perspective. Coupling between domains and how loading affects coupled energy are examined. Moreover, simple methods for experimentally extracting model parameters, including the coupling coefficient, are provided to empower engineers to quickly and easily integrate PTs in SPICE simulations for the rapid and improved development of PT interface circuits. The model and parameter extractions are validated by comparing them to the measured response of a physical cantilever-style PT excited by regular and irregular vibrations. In most cases, less than a 5–10% error between measured and simulated responses is observed.


Introduction
To increase ease of use, practicality, and desirability, many devices, especially in the Internet of Things, require minimal maintenance power methods.Numerous mechanical, biomedical, and civil systems produce vibrations, and piezoelectric transducers (PTs), with their high energy densities and off-the-shelf availability, can effectively and easily tap into this energy source [1][2][3][4][5].
PTs convert mechanical energy, in the form of motion or vibration, to electrical energy that can be harvested to charge a battery, v B , and power a microsystem, as illustrated in Figure 1.Ultimately, maximizing the harvester's output power, P O , is most important, but to do so, an understanding of how the PT delivers power, P PZ , is required [6][7][8].

Introduction
To increase ease of use, practicality, and desirability, many devices, especially in the Internet of Things, require minimal maintenance power methods.Numerous mechanical, biomedical, and civil systems produce vibrations, and piezoelectric transducers (PTs), with their high energy densities and off-the-shelf availability, can effectively and easily tap into this energy source [1][2][3][4][5].
PTs convert mechanical energy, in the form of motion or vibration, to electrical energy that can be harvested to charge a battery, v , and power a microsystem, as illustrated in Figure 1.Ultimately, maximizing the harvester's output power, P , is most important, but to do so, an understanding of how the PT delivers power, P , is required [6][7][8].systems well.In a PT, an understanding of how mechanical motion produces electrical energy, electrical loading affects output power, and energy couples between domains is necessary when designing harvesters effectively for specific motion inputs.This requires an insightful PT model with parameters that can effectively approximate the physical device.The model shown in [9][10][11] provides a circuit representation of both the mechanical and electrical domains of a PT, along with a method of coupling between the two to capture loading effects.However, the transformer's coupling coefficient and turns ratio are combined into a single coefficient, which removes insight and analysis potential and can lead to inaccuracies.It is also known that PTs exhibit a dual resonance (parallel and series) behavior that is not captured [12][13][14][15].In other words, that model does not adequately model PT output impedance.Moreover, mechanical engineers often use mechanical methods, material properties, and/or finite element analysis to determine model parameters, which may be inaccessible to many electrical engineers [9][10][11].
Other models exclude the mechanical domain and approximate i PZ to be proportional to the PT's tip velocity (for cantilevered PTs).In these cases, the PT is assumed to be very weakly coupled, so the electrical load is said to have negligible effect on the PT's motion [1,4,7,[16][17][18][19][20].It is not clear though, what constitutes weak coupling or how a PT is determined to be weakly coupled.DC-DC converters using PTs, such as those described in [12][13][14][15], require good impedance models for accurate converter and control design.Such applications often use the Van Dyke model, but this model is lumped, leading to potential error in converter design and PT impedance understanding, and its derivation could be expanded on [21].This model is typically used for "strongly coupled" PTs, which poses similar ambiguity as before.
This paper presents a complete electromechanical (EM) model for PTs independent of strong or weak coupling assumptions, vibration conditions, and transducer geometries, making it generally applicable to many PT types, sizes, and applications.Because the model details enable insight into the power and energy transfers and losses occurring within PTs, it is particularly useful in energy harvesting contexts.PT impedances and resonances are also analyzed and discussed, exhibiting the model's use for PT resonators.The model is also reduced to its Norton and Thevenin equivalents, which are applicable in all coupling contexts.Simple methods for experimentally extracting model parameters are presented and discussed to empower designers to quickly and easily integrate PTs in SPICE simulations for the rapid and improved development of PT interface circuits.

Complete Electromechanical Model
The complete EM model of a PT presented in this section strives to capture and distinguish between the aspects of the PT responsible for motion, energy capture, domain conversion, and electrical behavior.Equivalents of the complete model are also presented for ease of use.The model is developed using a cantilever style PT but is generally valid.

Complete Model
In Figure 2, the PT mass, M T , has kinetic energy when in motion with velocity, v T .Similarly, a capacitor has energy when a voltage develops across it.So, M T can be modeled as a capacitor, where C T ≡ M T [kg], and the mass's velocity, in [m/s], is analogous to the capacitor voltage, v T [22,23]: Sensors 2024, 24, x FOR PEER REVIEW 3 of 21 A mass moves if force is applied to it and a capacitor develops a voltage if current is applied to it.Since  represents both voltage and velocity, its time derivative can represent the transducer's acceleration, a , and since  ≡ M , the current, i , into  is analogous to an applied mechanical force, f , on M in units of kg m/s or N .Therefore, dv T A mass moves if force is applied to it and a capacitor develops a voltage if current is applied to it.Since v T represents both voltage and velocity, its time derivative can represent the transducer's acceleration, a T , and since C T ≡ M T , the current, i S , into C T is analogous to an applied mechanical force, f S , on M T in units of kg m/s 2 or [N].Therefore, Springs store potential energy when compressed (or bent for a cantilever) a distance, d X , from rest.Inductors store energy when magnetized by ϕ L , thus an inductor can model a spring with spring constant, K T : where L T is inversely equal to K T [N/m] and i LT represents force applied by the spring.
Frictional sources, such as air resistance, dissipate energy and dampen mechanical motion in the same way resistors dissipate energy in electrical systems: Mechanical dampers, D T [kg/s], apply a force, f D , in response to motion, v T .So, in keeping with the convention established above, mechanical damping, D T , is analogous to the inverse of resistance, R T , (i.e., conductance, G T ) where i R is the current through the resistor and v T is the voltage across it, as seen in ( 4) [24].
The LC circuit analogous to the spring mass system described above will have the same resonant frequency, f R : The energy transfer between the mechanical domain (MD) and electrical domain (ED) of a PT is typically modeled with a transformer [9][10][11].Since PTs and piezoelectric (PZ) materials are not capable of capturing all mechanical energy present for conversion to the ED, a coupling coefficient, k CL , exists: k CL is intrinsic to the PT and is what makes a PT a PT.It represents the inherent strength of the connection between domains-the larger k CL , the better the connection.k CL can only be between 0 and 1, inclusive.A k CL of 0 corresponds to a material that is not piezoelectric as this implies no coupling exists between the ED and MD [1,25].Figure 3 is the complete EM PT model with the MD modeled on the left.In the state of the art (SoA), k CL is often lumped in with the transformer [9][10][11].In Figure 3, however, it is represented as the ratio between L TC and L T , as (6) shows, where L TC is the portion of L T that perfectly couples with the ED.k CL can also be thought of as the percentage of v T that couples to current in the ED, i EC , when the output is shorted.
Sensors 2024, 24, x FOR PEER REVIEW 4 of 21 In a conventional transformer, voltage on the secondary is a multiple of primary voltage by a factor of the turns ratio,  [26].Of course, there are no "turns" in a PT, so  is referred to here as the translation coefficient, representing the conversion between the MD (force/velocity) and the ED (current/voltage): In a conventional transformer, voltage on the secondary is a multiple of primary voltage by a factor of the turns ratio, k T [26].Of course, there are no "turns" in a PT, so k T is referred to here as the translation coefficient, representing the conversion between the MD (force/velocity) and the ED (current/voltage): In a PT, it is generally accepted that velocity begets current, i EC , so the two dependent current sources are used in conjunction with the transformer to couple v T to i EC [1,4,7,9,[16][17][18][19][20].The dependent current sources perform a voltage-to-current conversion which reciprocates impedances (i.e., resistance to conductance, inductance to capacitance, etc.) and dividers (i.e., voltage dividers to current dividers and vice versa).Lastly, C PZ models the PT's ability to store electrical energy, and R PZ models leakage.C PZ is often dominated by the capacitance that exists between the PT's electrodes [10,27].The parameters in Figure 3 are the most fundamental to the model, so they have been italicized in the text to assist the reader in tracking them throughout the paper.
Presenting the model as it is described here (as opposed to representing mass as inductance and stiffness as capacitance as is often performed in the SoA, such as in [9][10][11]) is vital for a few reasons.First, it allows the inductor to be incorporated into the transformer, illustrating the role the PT stiffness plays in domain coupling/energy transfer by representing the spring as the medium for EM conversion [9].It also enables explicit visualization of the transformer coupling coefficient by using split inductors to clearly show that only a fraction of the velocity, v TC , (i.e., mechanical energy) couples to the ED.Lastly, by using split inductors to represent transformer coupling, the transformer turns ratio can also be explicitly defined.Separating k CL from k T is important because if the coupling coefficient were lumped in with k T , it would be incorrectly applied to all impedance translations, as will be shown in Section 2.2.It also introduces the dual resonance behavior into the model, as will be shown in Section 2.3.
This model makes the same assumptions about PTs as those used in and accepted by the SoA, such as in [9][10][11][12][13][14][15][16][17][18][19][20].These assumptions are: (1) A PT which behaves as a spring-massdamper system in the mechanical domain is used; (2) Mechanical velocity is in phase with PT short circuit current; (3) The PT's electrical domain is predominately capacitive.The transformer and split inductors are methods for fundamentally representing conversion and coupling between domains-regardless of its geometry and material, a PT converts energy between domains and only some fraction of energy couples.Therefore, this model and the following analysis can be applied to any PT that satisfies the three assumptions, which encompasses the majority of PTs used in energy harvesting and many PTs used as resonators [1,[9][10][11][12][13][14][15][16][17][18][19][20].

Electrical Model
To simplify analysis, it is helpful to use a circuit equivalent to Figure 3 which does not include the transformer or dependent sources.The MD impedances, Z T , can be referred to the ED, Z ′ Z , using typical transformer techniques as seen in [26]: As mentioned above, the dependent current sources reciprocate impedances and dividers.The results of the above translations are shown in Figure 4.By representing the model this way, coupling between the MD and ED is modeled with a coupling capacitor, C LC ′ , and a leakage capacitor, C L0 ′ , which are related to k CL by (6).Note that if k CL and k T were lumped, all translated impedances, instead of just L T0 and L TC , would be affected by k CL .When k CL is greater than 50%, C LC ′ is larger than C L0 ′ , so more than half the mechanical energy present can couple to the ED.When less than 50%, the opposite is true.
The results of the above translations are shown in Figure 4.By representing the model this way, coupling between the MD and ED is modeled with a coupling capacitor, C ′, and a leakage capacitor, C ′, which are related to  by (6).Note that if  and  were lumped, all translated impedances, instead of just  and  , would be affected by  .When  is greater than 50%, C ′ is larger than C ′, so more than half the mechanical energy present can couple to the ED.When less than 50%, the opposite is true.

Lumped Electrical Model
Reducing the model to its Norton and/or Thevenin equivalent, shown in Figure 5, is helpful for simple and easy application.To create these equivalent circuits, an understanding of the PT's output impedance, Z , is required.Z can be determined by shorting v in Figure 4 and evaluating impedance when looking left from  [28]: At very low frequencies, L ′ is effectively a short relative to C ′ and C ′.As frequency increases though, its impedance will no longer be negligible, and it will interact with C ′ and C ′, resulting in a mechanical, series resonance, f , given in (5).The quality factor, Q , at f is:

Lumped Electrical Model
Reducing the model to its Norton and/or Thevenin equivalent, shown in Figure 5, is helpful for simple and easy application.To create these equivalent circuits, an understanding of the PT's output impedance, Z PZ , is required.
The results of the above translations are shown in Figure 4.By representing the model this way, coupling between the MD and ED is modeled with a coupling capacitor, C ′, and a leakage capacitor, C ′, which are related to  by (6).Note that if  and  were lumped, all translated impedances, instead of just  and  , would be affected by  .When  is greater than 50%, C ′ is larger than C ′, so more than half the mechanical energy present can couple to the ED.When less than 50%, the opposite is true.

Lumped Electrical Model
Reducing the model to its Norton and/or Thevenin equivalent, shown in Figure 5, is helpful for simple and easy application.To create these equivalent circuits, an understanding of the PT's output impedance, Z , is required.Z can be determined by shorting v in Figure 4 and evaluating impedance when looking left from  [28]: At very low frequencies, L ′ is effectively a short relative to C ′ and C ′.As frequency increases though, its impedance will no longer be negligible, and it will interact with C ′ and C ′, resulting in a mechanical, series resonance, f , given in (5).The quality factor, Q , at f is: Z PZ can be determined by shorting v S in Figure 4 and evaluating impedance when looking left from v PZ [28]:

At very low frequencies, L C
′ is effectively a short relative to C LC ′ and C L0 ′ .As frequency increases though, its impedance will no longer be negligible, and it will interact with C LC ′ and C L0 ′ , resulting in a mechanical, series resonance, f R , given in (5).The quality factor, Q R , at f R is: Above, f R , L C ′ shunts the surrounding capacitors until frequency increases enough for L C ′ to effectively open relative to them, resulting in an EM, parallel resonance, f R ′ : with a quality factor, Q R ′ , of In this paper, "⊕" indicates the electrical series combination of capacitors (which are mathematically parallel).

Below f R , L C
′ shorts, so C PZ and C LC ′ dominate: until frequency drops low enough for them to be open relative to R PZ , which occurs at p PZ , shown in Figure 6.Past f R ′ , capacitors are shorted and L C ′ is open, so Z PZ is dominated by the series combination of C LC ′ and C L0 ′ in parallel with C PZ : Note that f R ′ and Q R ′ depend on C PZ .Specifically, C PZ affects C LC ′ (or L TC ), indicating that it directly alters the spring's mechanical stiffness [9].f R and Q R , however, are independent of the ED.Z PZ 's response across frequency is shown in Figure 6, but note that each attribute of the response is variable.f R and f R ′ will be more separated at large k CL s and come closer together as k CL reduces.As C PZ gets larger relative to C L0 ′ and C LC ′ , it will dominate the response, and the prominence (i.e., quality factors) of and spacing (i.e., k F , shown in ( 11)) between the f R and f R ′ peaks will reduce.When f R and f R ′ get too close together, Z PZ 's peaks may no longer correspond to f R and f R ′ .
with a quality factor, Q ′, of In this paper, "⊕" indicates the electrical series combination of capacitors (which are mathematically parallel).
Below f , L ′ shorts, so  and C ′ dominate: until frequency drops low enough for them to be open relative to  , which occurs at p , shown in Figure 6.Past f ′, capacitors are shorted and L ′ is open, so Z is dominated by the series combination of C ′ and C ′ in parallel with  : Note that f ′ and Q ′ depend on  .Specifically,  affects C ′ (or  ), indicating that it directly alters the spring's mechanical stiffness [9].f and Q , however, are independent of the ED.
Z 's response across frequency is shown in Figure 6, but note that each attribute of the response is variable.f and f ′ will be more separated at large  s and come closer together as  reduces.As  gets larger relative to C ′ and C ′, it will dominate the response, and the prominence (i.e., quality factors) of and spacing (i.e., k , shown in ( 11)) between the f and f ′ peaks will reduce.When f and f ′ get too close together, Z 's peaks may no longer correspond to f and f ′.
Using (15), k CL can be approximated as: which is sometimes used for PZ resonators [25,29].This is useful but may not be applicable in many other contexts as the condition may not be satisfied.The Norton and Thevenin impedances of Figures 3 and 4 are equal to Z PZ .The Norton equivalent current source, i NO , is equal to the current through a short circuit (SC) across the output, v PZ , which in this case, is equal to the current through C LC ′ [28].With the output shorted, L C ′ , R R ′ , and C L ′ are in series.Using the translations in Section 2.2, the current through that series combination is proportional to v T by k T : The current through C LC ′ is a current divided fraction, k CL , of C L ′ , as seen in ( 6).The Thevenin equivalent voltage source is related to i NO and Z PZ : Sensors 2024, 24, 4367 7 of 20

Additional Resonant Tones
A PT can resonate at multiple frequencies depending on its geometry.The shape it deforms into at each resonant frequency is known as a mode of vibration, but in the context of electrical engineering, it adds a tone to the frequency spectrum [30].
Each tone adds resonance in the same way as the fundamental, so adding LCR loops coupled to the ED via transformers accounts for them, as shown in Figure 7.Note that the current/force source in the additional tones is equal to the fundamental's source, and that the mechanical parameters, k T , and k CL are likely different for each tone [10].
The current through C ′ is a current divided fraction,  , of C ′, as seen in ( 6).The Thevenin equivalent voltage source is related to i and Z :

Additional Resonant Tones
A PT can resonate at multiple frequencies depending on its geometry.The shape it deforms into at each resonant frequency is known as a mode of vibration, but in the context of electrical engineering, it adds a tone to the frequency spectrum [30].
Each tone adds resonance in the same way as the fundamental, so adding LCR loops coupled to the ED via transformers accounts for them, as shown in Figure 7.Note that the current/force source in the additional tones is equal to the fundamental's source, and that the mechanical parameters,  , and  are likely different for each tone [10].While the  components of each tone experience the same current, the different   mean different current is induced to each tone.The tones' unique impedances also mean the ac part of  affects each tone's voltage differently.When translated, each tone adds an LCR branch which couples to the electrical load via C ′, as shown in Figure 8.In this case, each branch has the same voltage, but the voltage coupled to each tone is unique because the C ′s are different.While the L EC components of each tone experience the same current, the different k T s mean different current is induced to each tone.The tones' unique impedances also mean the ac part of i LE affects each tone's voltage differently.When translated, each tone adds an LCR branch which couples to the electrical load via C LC ′ , as shown in Figure 8.In this case, each branch has the same voltage, but the voltage coupled to each tone is unique because the C LC ′ s are different.
Sensors 2024, 24, x FOR PEER REVIEW 8 of 21 Interestingly, this means that tones add load capacitance relative to other tones.In other words, if looking at a particular tone, i, any tone greater or smaller than i will add to i's load capacitance, as shown in Figure 9, (20) and (21).For tones lower than i, L ′ is open circuited (OC), so the impedance of those tones is dominated by the series combination of C ′ and C ′ of each tone: The impedance of each tone appears parallel to each other at the output, so the combined load impedance they contribute is: For tones higher than i, L ′ is a short, so those tones contribute a load impedance, Z , dominated by C ′ of each tone: Interestingly, this means that tones add load capacitance relative to other tones.In other words, if looking at a particular tone, i, any tone greater or smaller than i will add to i's load capacitance, as shown in Figure 9, (20) and (21).

Coupling Coefficients
Sensors 2024, 24, x FOR PEER REVIEW 8 of 21 Interestingly, this means that tones add load capacitance relative to other tones.In other words, if looking at a particular tone, i, any tone greater or smaller than i will add to i's load capacitance, as shown in Figure 9, (20) and (21).For tones lower than i, L ′ is open circuited (OC), so the impedance of those tones is dominated by the series combination of C ′ and C ′ of each tone: The impedance of each tone appears parallel to each other at the output, so the combined load impedance they contribute is: For tones higher than i, L ′ is a short, so those tones contribute a load impedance, Z , dominated by C ′ of each tone: For tones lower than i, L Cn ′ is open circuited (OC), so the impedance of those tones is dominated by the series combination of C LCn ′ and C L0n ′ of each tone:

Coupling Coefficients
The impedance of each tone appears parallel to each other at the output, so the combined load impedance they contribute is: For tones higher than i, L Cn ′ is a short, so those tones contribute a load impedance, Z LD(H) , dominated by C LCn ′ of each tone: Sensors 2024, 24, 4367

Coupling Coefficients
The goal of a PT coupling coefficient is to express the fraction of energy that reaches the ED from the MD (or vice versa).The energy that does not reach the ED is not necessarily lost, as it could remain stored in the MD [31].To determine the fraction of energy in the MD that reaches the ED, two components of cross-domain energy transfer need to be discussed.First is the fraction of mechanical energy that is available to the ED.The second arises from the fact that only a fraction of that available energy can be captured and used by the ED.
The fraction of mechanical energy available to the ED is directly related to k CL since it represents the PT's inherent ability to convert energy between domains.To capture energy in the ED, though, an electrical load, such as C PZ , needs to be present, but this reduces the fraction of energy available to the ED since impedance is added to the coupling network [31].In other words, k CL is effectively reduced in the presence of an electrical load.In Figure 4, the series combination of C PZ and C LC ′ forms a current divider with C L0 ′ , so only a fraction of the power entering the capacitor network, P LT , reaches the EM branch, P EM , thus less mechanical energy is available to the ED.This effect is captured by the mechanical coupling coefficient: Mechanically, C PZ effectively increases the spring's stiffness (i.e., reduces L TC ), reducing the energy it can store for a given force packet, meaning less is available for coupling [9,31].k CM is less than or equal to k CL .
When OC, the energy stored in C PZ can be considered as having reached the ED since this is energy that the ED can use.Because L TC is finite, the transformer in Figure 3 requires some energy to transform voltage when loaded.This energy does not reach the ED, but is also not lost, it simply remains in the MD [26].In Figure 4, C PZ forms a voltage divider with C LC ′ , so only a fraction of P EM reaches C PZ .In other words, some energy remains in C LC ′ and does not reach C PZ .This effect is captured with the electrical coupling coefficient: which ideally equals 1. Mechanically, this occurs because the coupled portion of the spring's stiffness is shared between the MD and ED, so an applied force will result in energy that is partially stored in the coupled mechanical spring, L TC , and partially in the "electrical spring", C PZ [9].k CL represents the maximum fraction of mechanical energy available to the ED.k CM pertains to what fraction of that maximum mechanical energy is actually available to the ED when a load is present.k CE pertains to what fraction of the available energy is captured by the ED.Note that k CM and k CE are load dependent.The product of these two is the fraction of mechanical energy that is useable/harvestable in the ED and is referred to here as the PT coupling coefficient: Theoretically, k C cannot be larger than k CL , regardless of loading conditions, as k CL is a material limit.

Coupling Extremes
As k CL decreases, it is harder for actions in the ED to affect motion (i.e., it takes a higher v PZ to affect v LT ).From this perspective, a small k CL implies a "weakly coupled" PT.This is not to say, though, that the ED cannot have significant effects on motion when k CL is low.If C PZ is small enough, even a small i EC can produce a large v PZ .It is important Sensors 2024, 24, 4367 9 of 20 to note that coupling effects are ultimately determined by the combined effects of L T , k CL , and k T 2 relative to the electrical load.Therefore, a small k C , which accounts for all those parameters, always means the PT is weakly coupled.If C LC ′ is very small, it becomes more difficult for energy to reach the ED because C PZ would need to be significantly smaller for k CE to approach 1 (assuming C PZ is the only electrical load).When k CE approaches 0, only a small fraction of energy is shared between domains, so the MD and ED are essentially disconnected.A very small k CL , though, does not necessarily mean that C LC ′ is small, as it is also a function of L T and k T 2 .
k CE approaching 0, however, always means C LC ′ is much smaller than C PZ , so: In other words, C PZ effectively shorts the MD.Since capacitance typically corresponds to large impedances at frequencies of interest for energy harvesting, a Norton equivalent is preferred to model this extreme.This is often performed in the SoA by saying the PT is "weakly coupled" without much further justification, but now ( 23) and ( 25) can be used definitively [1,4,7,[16][17][18][19][20].Note that PTs traditionally considered weakly coupled may not actually fulfill this approximation.Note also, that changing the load could change the applicability of the approximation.
If k CL approaches 1, the MD becomes fully connected to the ED, so Z PZ includes f R and f R ′ .k CM approaching 1 means the same thing but also implies C PZ is much bigger than C LC ′ and C L0 ′ .At f R , L C ′ and C LC ′ effectively short and C L0 ′ opens, so: since R PZ is typically big.R R ′ is usually small, so for these cases, a Thevenin equivalent would be preferred.At f R ′ , the resonant network exhibits a huge impedance parallel to R PZ , so: Since R PZ is normally large, this case is also best modeled with a Norton equivalent.

Coupling with a Capacitive Load
When a PT is loaded with a capacitor, C LD , such that the total output capacitance, C O , is: and assuming R PZ is very large, an interesting tradeoff exists between k CM and k CE .To maximize k CM , C O should be much larger than C LC ′ and C L0 ′ , but to maximize k CE , C O should be much smaller than C LC ′ and C ′ L0 .So, if k CM is maximized, k CE , and therefore k C , go to 0, and vice versa.This implies an optimal C O exists to maximize k C , C O(MKC) .
Using ( 21)-( 24) and approximating R PZ as large, an expression for k C with respect to C O can be made: k C is evaluated at f R ′ as this is the frequency at which the maximum voltage occurs and the frequency the beam freely vibrates at.To find C O(MKC) , the derivative of (29) with respect to C O is equated to 0: Solving for C O in (30) gives the optimal C O to maximize k C for a particular k CL : Plugging (31) back into (29) gives an expression for the maximum PT coupling coefficient that can be achieved for a PT, k C(MAX) , with a particular k CL and capacitive load: This highlights the difference between a PT's intrinsic ability to couple between domains, k CL , and the fraction of energy that can be practically captured in the ED, k C .It also highlights the effect electrical loading can have on available electrical energy, since in this case, k C(MAX) is always less than k CL , unless k CL equals 1.The k C(MAX) line for a given k CL can be seen in Figure 10.If C O is larger than C O(MKC) , then the PT's motion is increased, so more energy is available, but the ability to capture that energy in the ED is decreased, underdamping the system.If C O is smaller than C O(MKC) , the ability to capture energy in the ED is increased, but mechanical motion is decreased, overdamping the system.Optimizing k C will always improve electrical energy capture, but as k CL becomes very small (much less than 1% for example), the conditions required to achieve optimization could become impractical.
Sensors 2024, 24, x FOR PEER REVIEW 11 of 21 also highlights the effect electrical loading can have on available electrical energy, since in this case, k is always less than  , unless  equals 1.The k line for a given  can be seen in Figure 10.If C is larger than C , then the PT's motion is increased, so more energy is available, but the ability to capture that energy in the ED is decreased, underdamping the system.If C is smaller than C , the ability to capture energy in the ED is increased, but mechanical motion is decreased, overdamping the system.Optimizing k will always improve electrical energy capture, but as  becomes very small (much less than 1% for example), the conditions required to achieve optimization could become impractical.

Parameter Extraction
The model and analysis presented in the previous sections are applicable to many PT types and dimensions, as discussed in Section 2.1.This section describes a method for extracting each parameter of the model specifically from a real cantilever-style PT with a tip mass to illustrate a way to use the model practically.The methods can be expanded upon in future work to develop extraction methods for other types of PTs and contexts.The method begins by assuming a single-tone PT and then provides techniques for accounting for higher tones.

Mechanical Parameters
The mechanical parameters are arguably the simplest to extract, so they are determined first.If the PT's tip mass, M , is much larger than the PT's mass, then [22,23]: Therefore,  can be determined by weighing the tip mass.Since the PT and mass form a spring-mass resonator,  is related to  by f .To measure f , SC the PT's output so the ED is not loading the MD and measure the period of d and/or  of an initially displaced beam using a displacement sensor, as shown in Figure 11.Measuring the response of an initially displaced beam allows the PT to resonate freely at its f .

Parameter Extraction
The model and analysis presented in the previous sections are applicable to many PT types and dimensions, as discussed in Section 2.1.This section describes a method for extracting each parameter of the model specifically from a real cantilever-style PT with a tip mass to illustrate a way to use the model practically.The methods can be expanded upon in future work to develop extraction methods for other types of PTs and contexts.The method begins by assuming a single-tone PT and then provides techniques for accounting for higher tones.

Mechanical Parameters
The mechanical parameters are arguably the simplest to extract, so they are determined first.If the PT's tip mass, M T , is much larger than the PT's mass, then [22,23]: Therefore, C T can be determined by weighing the tip mass.Since the PT and mass form a spring-mass resonator, L T is related to C T by f R .To measure f R , SC the PT's output so the ED is not loading the MD and measure the period of d X and/or v T of an initially displaced beam using a displacement sensor, as shown in Figure 11.Measuring the response of an initially displaced beam allows the PT to resonate freely at its f R . and f can then be used to determine  : The PT-mass system exhibits damped resonance, so its behavior is of the form: where  represents the initial force in the spring resulting from the initial displacement, d .By choosing two different points along the d curve, such as (d , t ) and (d , t ) as shown in Figure 11, two different instances of (35) are determined.Taking the ratio of these instances and rearranging the expression to isolate  results in: . (36)

Electrical Parameters
Extracting the electrical parameters is a little more involved.Recall from Section 2.3 that past f ′, L ′ is open, so the measured impedance, Z , at these frequencies would be: (37) Note that  cannot be measured directly, but C can be.
There are two main challenges with using an impedance analyzer (IA) to determine  : output impedance and minimum frequency. is typically very large (on the order of MΩ), so the IA's output impedance needs to be larger than that to give an accurate measure of  .Plus, since  is so large, p typically occurs around or below 1Hz, so the IA needs to operate at those frequencies to provide a measurement of  .If an IA is on hand which satisfies these requirements, that is the best method for determining  .If such an IA is not accessible, the PT's discharge rate, shown in Figure 12, can be used since a PT's output is largely capacitive.C T and f R can then be used to determine L T : The PT-mass system exhibits damped resonance, so its behavior is of the form: where i LT(I) represents the initial force in the spring resulting from the initial displacement, d X(I) .By choosing two different points along the d X curve, such as (d X1 , t 1 ) and (d X2 , t 2 ) as shown in Figure 11, two different instances of (35) are determined.Taking the ratio of these instances and rearranging the expression to isolate R T results in: . (36)

Electrical Parameters
Extracting the electrical parameters is a little more involved.Recall from Section 2.3 that past f R ′ , L C ′ is open, so the measured impedance, Z O , at these frequencies would be: Note that C PZ cannot be measured directly, but C HF can be.
There are two main challenges with using an impedance analyzer (IA) to determine R PZ : output impedance and minimum frequency.R PZ is typically very large (on the order of MΩ), so the IA's output impedance needs to be larger than that to give an accurate measure of R PZ .Plus, since R PZ is so large, p PZ typically occurs around or below 1 Hz, so the IA needs to operate at those frequencies to provide a measurement of R PZ .If an IA is on hand which satisfies these requirements, that is the best method for determining R PZ .If such an IA is not accessible, the PT's discharge rate, shown in Figure 12, can be used since a PT's output is largely capacitive.Begin by immobilizing the beam-this ensures Z is dominated by C across all frequencies greater than p by forcing  to 0. Then connect the PT to a power supply set to a particular voltage and allow it to charge.Disconnect the PT from the supply and measure the voltage across the PT's terminals,  , with an oscilloscope.The time it takes for the voltage to drop from the initial voltage,  , to 37% of  is one time con- Begin by immobilizing the beam-this ensures Z O is dominated by C HF across all frequencies greater than p PZ by forcing v T to 0. Then connect the PT to a power supply set to a particular voltage and allow it to charge.Disconnect the PT from the supply and measure the voltage across the PT's terminals, v PZ , with an oscilloscope.The time it takes for the voltage to drop from the initial voltage, v PZ(I) , to 37% of v PZ(I) is one time constant, τ RC .R PZ can be determined using: Rearranging (37) and substituting in (17) for k T gives an expression for C PZ : where i PZ,PK(SC) is the measured, peak, SC current corresponding to a particular measured, peak, SC velocity, v T,PK(SC) .v T,PK(SC) can be measured using a d X sensor.There are two unknowns in (39): C PZ and k CL , but it is known that k CL is between 0 and 1.To determine C PZ and k CL , software aided curve fitting is performed on the OC v PZ , shown in Figure 13.
Curve fitting computes a selection of variables to yield the least error within given boundary conditions [32,33].( 39) is used as a boundary condition to ensure the results apply to the PT. ( 39) and v PZ(OC) curve fitting essentially act as two equations to solve for the two unknowns: C PZ and k CL .Only k CL is adjusted during fitting since all other parameters have been determined, and because the bounds are well defined, there should only be one unique solution.After k CL is determined from the fit, it is plugged into (39) to find C PZ .k T can then be found with (17) using k CL , i PZ,PK(SC) and v T,PK(SC) .Begin by immobilizing the beam-this ensures Z is dominated by C across all frequencies greater than p by forcing  to 0. Then connect the PT to a power supply set to a particular voltage and allow it to charge.Disconnect the PT from the supply and measure the voltage across the PT's terminals,  , with an oscilloscope.The time it takes for the voltage to drop from the initial voltage,  , to 37% of  is one time constant, τ . can be determined using: Rearranging (37) and substituting in (17) for k gives an expression for  : where  , is the measured, peak, SC current corresponding to a particular measured, peak, SC velocity,  , . , can be measured using a d sensor.There are two unknowns in (39):  and  , but it is known that  is between 0 and 1.To determine  and  , software aided curve fitting is performed on the OC  , shown in Figure 13.Curve fitting computes a selection of variables to yield the least error within given boundary conditions [32,33].( 39) is used as a boundary condition to ensure the results apply to the PT.(39) and  curve fitting essentially act as two equations to solve for the two unknowns:  and  .Only  is adjusted during fitting since all other parameters have been determined, and because the bounds are well defined, there should only be one unique solution.After  is determined from the fit, it is plugged into (39) to find  . can then be found with ( 17) using  ,  , and  , .The relationship between  and  is described in Section 2.1.To restate,  can be used to determine  with: and  with: The relationship between k CL and L T is described in Section 2.1.To restate, k CL can be used to determine L TC with: and L T0 with:

Approximating 2nd Tone
As discussed in Section 2.4, higher tones can affect PT behavior, so it can be important to account for them.Approximating the second tone is discussed here, but the methods could be extended to higher tones.Unfortunately, it is difficult to extract these parameters since C T2 is not equal to M T .
If operating the PT with a constant applied force of frequency f ′ R , no energy is supplied by the higher tones, so they simply act as capacitive loads.In this case, use C HF1 , shown in Figure 14, in (39).This way, C PZ will include the capacitive effects contributed by the higher tones, providing the correct loading with respect to the MD but overestimating the ED's energy storage capacity.could be extended to higher tones.Unfortunately, it is difficult to extract these parameters since C is not equal to M .
If operating the PT with a constant applied force of frequency f ′, no energy is supplied by the higher tones, so they simply act as capacitive loads.In this case, use C , shown in Figure 14, in (39).This way,  will include the capacitive effects contributed by the higher tones, providing the correct loading with respect to the MD but overestimating the ED's energy storage capacity.If working with a step function input, the higher tones supply energy to the ED, so it is ideal to account for their full models.While f ′ can be measured, C is unknown, so L cannot be determined.C and  , can also be measured, but  , can be difficult to measure since a high-speed, high-sensitivity displacement sensor is required.So, to approximate R , C , L , k , and k , parameter estimation/curve fitting can be used [32,33]. has a similar limiting condition as before, but this time, the first tone capacitance is accounted for, and more is unknown: Note that while curve fitting here should provide a good approximation of the effects of the second tone, since so many variables are being fitted and the limiting conditions are not as rigorous or as well-defined as for the first tone, the resulting second tone parameters may not necessarily correspond to the actual device parameters, and it is possible that there could be more than one unique solution.However, if there is a relationship between C and M and a good d sensor is on hand, extraction could be performed like what is described in A and B.

Test Setup
To validate the model described in Section 2, the methods described in Section 4 are performed on the MIDE S129 PT to extract its model parameters using the test setup shown in Figure 15.SPICE simulations are then performed on the model with extracted parameters and the results are compared with measured data.If working with a step function input, the higher tones supply energy to the ED, so it is ideal to account for their full models.While f ′ R2 can be measured, C T2 is unknown, so L T2 cannot be determined.C HF2 and i PZ,PK(SC) can also be measured, but v T,PK(SC) can be difficult to measure since a high-speed, high-sensitivity displacement sensor is required.So, to approximate R T2 , C T2 , L T2 , k T2 , and k CL2 , parameter estimation/curve fitting can be used [32,33].C PZ has a similar limiting condition as before, but this time, the first tone capacitance is accounted for, and more is unknown: Note that while curve fitting here should provide a good approximation of the effects of the second tone, since so many variables are being fitted and the limiting conditions are not as rigorous or as well-defined as for the first tone, the resulting second tone parameters may not necessarily correspond to the actual device parameters, and it is possible that there could be more than one unique solution.However, if there is a relationship between C T2 and M T and a good d X sensor is on hand, extraction could be performed like what is described in A and B.

Test Setup
To validate the model described in Section 2, the methods described in Section 4 are performed on the MIDE S129 PT to extract its model parameters using the test setup shown in Figure 15.SPICE simulations are then performed on the model with extracted parameters and the results are compared with measured data.The PT is mounted on a Bruel & Kjaer shaker with a 36 g weight attached to its tip.The frequency response analyzer (FRA) is used as the IA.The picoammeter is used to measure i and to SC the PT.The oscilloscope measures  .The PT is 2.795 × 0.407 × 0.029 in.

Parameter Extraction
The MIDE S129 PT only weighs about 1.4 g, so with a 36 g tip mass, (33) holds true.So,  is 36 mF ("milli" is used because the SI unit of mass is kg).d of the PT tip is measured using the Keyence d sensor shown in Figure 15.While shorted, the PT is bent to an initial displacement, released, and allowed to vibrate freely.Figure 11 shows the measured d and  , for the PT which is used to extract  and  .The PT is mounted on a Bruel & Kjaer shaker with a 36 g weight attached to its tip.The frequency response analyzer (FRA) is used as the IA.The picoammeter is used to measure i EC and to SC the PT.The oscilloscope measures v PZ .The PT is 2.795 × 0.407 × 0.029 in.

Parameter Extraction
The MIDE S129 PT only weighs about 1.4 g, so with a 36 g tip mass, (33) holds true.So, C T is 36 mF ("milli" is used because the SI unit of mass is kg).d X of the PT tip is measured using the Keyence d X sensor shown in Figure 15.While shorted, the PT is bent to an initial displacement, released, and allowed to vibrate freely.Figure 11 shows the measured d X and v T , for the PT which is used to extract L T and R T .
R PZ is determined using τ PZ measured in Figure 12 and the immobilized capacitance of the PT.The PT is immobilized by clamping it to a table and is charged using the power supply shown in Figure 15.The discharge is measured using the oscilloscope in Figure 15.Since the MIDE S129 with a 36 g tip mass has multiple tones, as seen in Figure 16, C HF1 is used in (39) and C HF2 is used when approximating the second tone parameters.Recall from Section 4.3, though, that using C HF1 in (39) means the extracted C PZ includes capacitive loading effects from the higher tones in Figure 16.The PT is mounted on a Bruel & Kjaer shaker with a 36 g weight attached to its tip.The frequency response analyzer (FRA) is used as the IA.The picoammeter is used to measure i and to SC the PT.The oscilloscope measures  .The PT is 2.795 × 0.407 × 0.029 in.

Parameter Extraction
The MIDE S129 PT only weighs about 1.4 g, so with a 36 g tip mass, (33) holds true.So,  is 36 mF ("milli" is used because the SI unit of mass is kg).d of the PT tip is measured using the Keyence d sensor shown in Figure 15.While shorted, the PT is bent to an initial displacement, released, and allowed to vibrate freely.Figure 11 shows the measured d and  , for the PT which is used to extract  and  .
is determined using τ measured in Figure 12 and the immobilized capacitance of the PT.The PT is immobilized by clamping it to a table and is charged using the power supply shown in Figure 15.The discharge is measured using the oscilloscope in Figure 15.Since the MIDE S129 with a 36 g tip mass has multiple tones, as seen in Figure 16, C is used in (39) and C is used when approximating the second tone parameters.Recall from Section 4.3, though, that using C in (39) means the extracted  includes capacitive loading effects from the higher tones in Figure 16.The measured  , and  , are shown in Figure 17.All variables required for (39) have now been measured and/or determined, so curve fitting can commence.The results of the curve fit are shown in Figure 13.All extracted model parameters are shown in Table 1.Plugging these into (23), k for this PT loaded by its  is 10%, which is not small enough to assume it is weakly coupled, as (25) shows (k would need to be much less than 1%), contrary to what would be performed in the SoA.The measured i PZ,PK(SC) and v T,PK(SC) are shown in Figure 17.All variables required for (39) have now been measured and/or determined, so curve fitting can commence.The results of the curve fit are shown in Figure 13.All extracted model parameters are shown in Table 1.Plugging these into (23), k CE for this PT loaded by its C PZ is 10%, which is not small enough to assume it is weakly coupled, as (25) shows (k CE would need to be much less than 1%), contrary to what would be performed in the SoA.

Resonant Vibrations
In this section, the PT's response to a constant applied vibration at f ′ is considered and compared to simulated model results.In this case, the force source, i , is not known, but  is, so the Norton model shown in Figure 5 is used for simulations.

Resonant Vibrations
In this section, the PT's response to a constant applied vibration at f R ' is considered and compared to simulated model results.In this case, the force source, i S , is not known, but v T(SC) is, so the Norton model shown in Figure 5 is used for simulations.
Consider first the SC case, shown in Figure 17, where measured and simulated i PZ are compared.The error between simulated and measured data is found with: and its average is plotted at the bottom of Figure 17.The areas highlighted in gray indicate where percent error exceeds ±15%.Note that error spikes every time i PZ crosses 0 because if measured data equals 0 and simulated data do not (which happens with the smallest difference in period), error spikes to infinity.
Figure 18 shows simulated and measured data for the OC PT.Note the difference in period, t R and t R ′ , between the SC and OC cases, respectively, which results from C PZ 's inclusion when OC, as discussed in Section 2.3 [9].t R ′ is the period of an f R ′ vibration.Also note that i PZ is in phase with v T while v PZ is in phase with d X [1,4,7,9,[16][17][18][19][20].Error V is smaller than Error I because v PZ(OC) is used for curve fitting, so any parameter errors will be more apparent in the i PZ(SC) plots.In Figure 19, resistors are applied across the PT to test the model's ability to predict accurate responses under various loading conditions.Energy burned by a load resistor, E , is measured across a one-second time interval when the PT experiences the same  as in Figure 17.Error between simulated and measured response likely comes largely from errors in measurements used for parameter extraction.Measuring C , for example, can be difficult since it is flanked by resonances, so the impedance line is slightly variable.Small rounding differences when measuring f can also lead to large error differences since it is squared when calculating  .Certain assumptions/approximations made during parameter extraction, such as neglecting the PT's mass and approximating impedance as purely capacitive at C can also contribute error.Errors introduced by measurement equipment, such as quantization noise and electronic noise, and the test setup, such as parasitics, can also add to total error.In Figure 19, resistors are applied across the PT to test the model's ability to predict accurate responses under various loading conditions.Energy burned by a load resistor, E R , is measured across a one-second time interval when the PT experiences the same v T(SC) as in Figure 17.In Figure 19, resistors are applied across the PT to test the model's ability to predi accurate responses under various loading conditions.Energy burned by a load resisto E , is measured across a one-second time interval when the PT experiences the sam  as in Figure 17.Error between simulated and measured response likely comes largely from errors measurements used for parameter extraction.Measuring C , for example, can be diffi cult since it is flanked by resonances, so the impedance line is slightly variable.Sma rounding differences when measuring f can also lead to large error differences since is squared when calculating  .Certain assumptions/approximations made during p rameter extraction, such as neglecting the PT's mass and approximating impedance Error between simulated and measured response likely comes largely from errors in measurements used for parameter extraction.Measuring C HF1 , for example, can be difficult since it is flanked by resonances, so the impedance line is slightly variable.Small rounding differences when measuring f R can also lead to large error differences since it is squared when calculating L T .Certain assumptions/approximations made during parameter extraction, such as neglecting the PT's mass and approximating impedance as purely capacitive at C HF1 can also contribute error.Errors introduced by measurement equipment, such as quantization noise and electronic noise, and the test setup, such as parasitics, can also add to total error.

Step Response
This section discusses the PT's response to a step input.The PT is bent and clamped to have an initial displacement, then released to resonate freely.This can be modeled in Figure 3 as i S having a constant, non-zero value long enough for the system to reach steady state.At steady state, all the current (force) goes through L T (the spring), which models the PT's spring storing energy while displaced.To model releasing the PT, drop i S to 0 very quickly.Note this is not a perfect step response as that would require an instantaneous beam release.Instead, there is some "ramp" as i S drops to 0 in a finite amount of time.
Since the input is approximately a step function, it is ideal to model higher tones.Table 2 shows the approximated 2nd tone parameters found using methods described in Section 4.3.Remember from Section 4.3 that these 2nd tone parameters may not correspond to actual device parameters but provide good approximations for device behavior and are included here to demonstrate the model's ability to account for higher tones.C PZ also reduces to 19.6 nF since it no longer includes the second tone's load, per (42).A perfect step function contains many frequencies, but as the step becomes less ideal (i.e., it takes more time to drop to 0), power in the higher frequencies reduces [34].In this experiment, an ideal beam release was not achieved, so tones higher than the second do not receive enough energy to have significant effects and are thus neglected.Figure 20 shows the response for the SC PT released from an initial displacement, d X(I) , which corresponds to an initial force, i S(I) .Note that i S(I) can be determined with (35) because all the current (force) goes through L T .Since the PT is SC, it oscillates at f R .The effects of the second tone can also be seen, especially in Figure 20b, and its frequency, f R2 , can be measured.The errors with and without the second tone approximation are shown.After the second tone's effects subside, the errors are the same and always less than ±10%.However, including the second tone greatly reduces error while its effect is present, indicating the model's ability to accurately account for higher tones.
Figure 21 shows the response for the OC PT released from approximately the same d X(I) , but this time the PT oscillates at f R ′ since it is OC.Again, the difference in error with and without the second-order approximation is only present for the first few cycles, but it is significant during those cycles.The reason the error seems to decrease with time is most likely a result of slight frequency error.Interestingly, the simulated voltage is higher (more negative) if the second tone is not accounted for (as is corroborated with the Error V plot).This implies that the second tone drains C PZ , effectively opposing the first tone [9].
Figure 20 shows the response for the SC PT released from an initial displacement, d , which corresponds to an initial force,  .Note that  can be determined with (35) because all the current (force) goes through  .Since the PT is SC, it oscillates at f .The effects of the second tone can also be seen, especially in Figure 20b, and its frequency, f , can be measured.The errors with and without the second tone approximation are shown.After the second tone's effects subside, the errors are the same and always less than 10%.However, including the second tone greatly reduces error while its effect is present, indicating the model's ability to accurately account for higher tones.Figure 21 shows the response for the OC PT released from approximately the same d , but this time the PT oscillates at f ′ since it is OC.Again, the difference in error with and without the second-order approximation is only present for the first few cycles, but it is significant during those cycles.The reason the error seems to decrease with time is most likely a result of slight frequency error.Interestingly, the simulated voltage is higher (more negative) if the second tone is not accounted for (as is corroborated with the Error plot).This implies that the second tone drains  , effectively opposing the first tone [9]. Figure 22 tests the two-tone model's ability to predict accurate responses under various loading conditions when the PT is initially displaced.Energy burned by the resistor is measured after the beam has resonated for 9 s with a d of about 1 mm.The measured and simulated data show good agreement with only about an 8% error and the same MPP R .Excluding the second tone's effects yields about a 9% error.This difference is likely small because the 9-s window is much larger than the fraction of a second that the second tone is present.If the input were closer to an ideal step function, though, the difference is expected to be larger.

Conclusions
The complete EM model presented here provides insight into PT energy transfers which was not previously possible by explicitly defining PT parameters not previously represented individually nor completely.The model is generally applicable to many PT types and dimensions and does not rely on coupling assumptions, removing guesswork from designing PT interface circuits.  is split into its coupled and uncoupled parts, enabling explicit and accurate definition, visualization, and grasp of   and   .Incorporating   into the transformer illustrates its role and related losses in EM energy transfer.Total coupling between domains, limitations to the amount of energy harvestable, and  Figure 22 tests the two-tone model's ability to predict accurate responses under various loading conditions when the PT is initially displaced.Energy burned by the resistor is measured after the beam has resonated for 9 s with a d X(I) of about 1 mm.The measured and simulated data show good agreement with only about an 8% error and the same MPP R LOAD .Excluding the second tone's effects yields about a 9% error.This difference is likely small because the 9-s window is much larger than the fraction of a second that the second tone is present.If the input were closer to an ideal step function, though, the difference is expected to be larger.(more negative) if the second tone is not accounted for (as is corroborated with the Error plot).This implies that the second tone drains  , effectively opposing the first tone [9].  .Excluding the second tone's effects yields about a 9% error.This difference is likely small because the 9-s window is much larger than the fraction of a second that the second tone is present.If the input were closer to an ideal step function, though, the difference is expected to be larger.

Conclusions
The complete EM model presented here provides insight into PT energy transfers which was not previously possible by explicitly defining PT parameters not previously represented individually nor completely.The model is generally applicable to many PT types and dimensions and does not rely on coupling assumptions, removing guesswork

Conclusions
The complete EM model presented here provides insight into PT energy transfers which was not previously possible by explicitly defining PT parameters not previously represented individually nor completely.The model is generally applicable to many PT types and dimensions and does not rely on coupling assumptions, removing guesswork from designing PT interface circuits.L T is split into its coupled and uncoupled parts, enabling explicit and accurate definition, visualization, and grasp of k CL and k T .Incorporating L T into the transformer illustrates its role and related losses in EM energy transfer.Total coupling between domains, limitations to the amount of energy harvestable, and loading effects are shown to be dependent on k CL , k T 2 , the load, and the PT's stiffness.It is shown that when k CE approaches 0, the model can be reduced to a current source parallel to C PZ (typical for weakly coupled PTs).When k C approaches 1, a Thevenin equivalent is preferred.
When constant sinusoidal vibrations are applied to a PT, it is shown that the other tones simply load the fundamental, whereas when an irregular vibration is applied, the other tones can oppose or support energy generation.This model also adequately captures a PT's impedance behavior making it useful in many contexts, removing the need to switch to different models for different applications, and, in a sense, unifying concepts of various existing PT models into one complete model.
Parameter extraction methods for cantilever style PTs with a tip mass are presented and performed to illustrate a way to use the model practically.Using the extraction methods, designers can effectively and easily model off-the-shelf PTs, enabling easy integration of PTs into SPICE simulations to assist in harvester or interface circuit development.In most cases, less than 5-10% error is observed between measurements and the model with extracted parameters.The model's ability to capture responses of PTs used as both resonators and energy sources, which are generally accepted in the SoA, further supports it.
The model described here is generalized, but work can be performed to apply it to specific PTs by deriving expressions for each component in Figure 3 in terms of material, dielectric, and piezoelectric constants.Concepts from this paper and [9] and an understanding of transformer turns ratios and mutual inductance can help accomplish this.The resulting expressions would be material and dimension dependent and would provide more insight into the physical parameters contributing to energy generation and loss.

Figure 6 .
Figure 6.Simulated impedance vs. frequency for a single tone PT on a logarithmic scale.

Figure 6 .
Figure 6.Simulated impedance vs. frequency for a single tone PT on a logarithmic scale.If C PZ is very small relative to C L0 ′ and C LC ′ , Q and k F will increase until limited by other factors such as mechanical impedances, k T , and k CL , and f R ′ 's expression reduces to:

Figure 7 .
Figure 7. Complete model with additional tones.

Figure 7 .
Figure 7. Complete model with additional tones.

Figure 8 .
Figure 8. Electrical model with additional tones.

Figure 9 .
Figure 9. Electrical model of the ith tone with simplified tonal loading.

Figure 8 .
Figure 8. Electrical model with additional tones.

Figure 8 .
Figure 8. Electrical model with additional tones.

Figure 9 .
Figure 9. Electrical model of the ith tone with simplified tonal loading.

Figure 9 .
Figure 9. Electrical model of the ith tone with simplified tonal loading.

Figure 10 .
Figure 10.Simulated k at various k and C combinations.

Figure 10 .
Figure 10.Simulated k C at various k CL and C O combinations.

Figure 11 .
Figure 11.Measured/simulated d and v for initially displaced SC PT.

Figure 11 .
Figure 11.Measured/simulated d X and v T for initially displaced SC PT.

Figure 13 .
Figure 13.Measured v and modeled v after curve fitting with resulting percent error.

Figure 13 .
Figure 13.Measured v PZ(OC) and modeled v PZ(OC) after curve fitting with resulting percent error.

Figure 14 .
Figure 14.Simulated Z for multi-tone PT to show C measurements.

Figure 14 .
Figure 14.Simulated Z O for multi-tone PT to show C HF measurements.

Figure 16 .
Figure 16.Measured Z vs. frequency of MIDE S129 with 36 g tip mass using FRA.

Figure 16 .
Figure 16.Measured Z O vs. frequency of MIDE S129 with 36 g tip mass using FRA.

Sensors 2024 , 21 Figure 17 .
Figure 17.Measured v , d , and i and simulated i for constantly vibrating PT and average percent error.

Figure 17 .
Figure 17.Measured v T(SC) , d X(SC) , and i PZ(SC) and simulated i PZ(SC) for constantly vibrating PT and average percent error.

21 Figure 18 .
Figure 18.Measured v , d and v and simulated v for constantly vibrating PT and average percent error.
Figure 19 estimates the maximum power point (MPP) resistance to be about the same as what was measured and shows approximately a 9% error between the MPP energies, E .Deviations between simulated and measured data at very low R are attributed to the noise floor limits of the oscilloscope.Similarly, deviation at high R is most likely due to the limitation of the ammeter.

Figure 18 .
Figure 18.Measured v T(OC) , d X(OC) and v PZ(OC) and simulated v PZ(OC) for constantly vibrating PT and average percent error.

Figure 18 .
Figure 18.Measured v , d and v and simulated v for constantly vibrating P and average percent error.
Figure 19 estimates the maximum power point (MPP) resistance be about the same as what was measured and shows approximately a 9% error betwee the MPP energies, E .Deviations between simulated and measured data at very lo R are attributed to the noise floor limits of the oscilloscope.Similarly, deviation high R is most likely due to the limitation of the ammeter.

Figure 19 .
Figure 19.Measured/simulated E R vs. load for constantly vibrating PT.

Figure 20 .
Figure 20.(a) Measured d , v , and i for initially displaced SC PT and average percent i error compared with simulations.(b) Zoomed in on the gray region of (a).

Figure 20 .
Figure 20.(a) Measured d X , v T , and i PZ for initially displaced SC PT and average percent i PZ error compared with simulations.(b) Zoomed in on the gray region of (a).

Figure 21 .
Figure 21.(a) Measured d , v , and v for initially displaced OC PT and average percent v error compared with simulations.(b) Zoomed in on the gray region of (a).

Figure 22 .
Figure 22.Measured/simulated E vs. load for initially displaced PT.

Figure 21 .
Figure 21.(a) Measured d X , v T , and v PZ for initially displaced OC PT and average percent v PZ error compared with simulations.(b) Zoomed in on the gray region of (a).

Figure 21 .
Figure 21.(a) Measured d , v , and v for initially displaced OC PT and average percent v error compared with simulations.(b) Zoomed in on the gray region of (a).

Figure 22
Figure22tests the two-tone model's ability to predict accurate responses under various loading conditions when the PT is initially displaced.Energy burned by the resistor is measured after the beam has resonated for 9 s with a d of about 1 mm.The measured and simulated data show good agreement with only about an 8% error and the same MPP R .Excluding the second tone's effects yields about a 9% error.This difference is likely small because the 9-s window is much larger than the fraction of a second that the second tone is present.If the input were closer to an ideal step function, though, the difference is expected to be larger.

Figure 22 .
Figure 22.Measured/simulated E vs. load for initially displaced PT.

Figure 22 .
Figure 22.Measured/simulated E OUT vs. load for initially displaced PT.