We investigate a classical notion denoted as 'Motional impedance (Zmot)' first developed by Kennelly in 1912, then explored by many other researchers early in the 20th century. Due to long history to establish impedance properties, Zmot has not been clearly understood with respect to classical impedance properties which makes retroactive study of the Zmot worthy. In this study, we project Zmot's unique characteristics onto anti-reciprocity where every electro-mechanical system as modeled by a gyrator. This work also clarifies that Zmot is a type of transfer impedance, thus is not necessarily to be a minimum-phase nor a Positive-Real function. Based on simulation, this study shows a shunt eddy current loss on the electrical side of the system is one possible source of the negative real part of Zmot which leads negative real part in Zmot plots. By taking a Balance Armature Receiver (BAR) as a specific example of the electro-mechanical systems, we also explore the (non)linear approximation principal of the BAR based on the magnetic force and hysteresis characteristics. We believe that this study puts anti-reciprocal physics on both the empirical and theoretical basis.
Motional impedance, Anti-reciprocal systems
The concept of motional impedance (Zmot) was first suggested by Arthur E. Kennelly in 1912a. It is a unique characteristic of every electro-mechanical system such as earphones, loudspeakers, hybrid cars, etc. The Zmot is defined as the component of the electrical impedance proportional to the output velocity. From measurements of the electrical impedance at the electrical input, one may accurately estimate the mechanical or acoustic load impedance. This property alone makes it worth of study [1]. aKennelly was born in 1861 in India and worked with Edison starting in 1887 at the age of 26 years. He then became a professor of electrical engineering at both Harvard University and the Massachusetts Institute of Technology (MIT). In 1827, Georg Ohm generalized measurement relationship between applied voltage and current varying resistance (as known as Ohm's law suggested only for DC circuits initially). Kennelly was the first person to coin the term as impedance for alternating current (AC) circuits [2]b and the first to utilize complex numbers in the circuit theory. Campbell's work follows the work of Heaviside [3] who first suggested loaded lines which leads to the generalized wave filter invention based on quasistatic (QS) theory [4]. bThe properties of the impedance have been established over a long-time span, Ohm (1827, DC impedance) → Kennelly (1893, AC impedance) → Brune (1931, establishment of impedance properties). Note that Heaviside [3] defined the term 'impedance' in 1982, a few months later in April 1893, Kennelly extended the definition of impedance to AC circuits [2]. In 1931, Brune first characterized a driving-point impedance as a positive-real (PR) quantity (positive definite operator in matrix form), requiring that every physically realizable impedance cannot have a negative resistance [5]. Brune's impedance is consistently studied with classic circuit theories such as Kirchhoff's circuit laws (KCL and KVL) under the quasistatic (QS) condition. The QS assumption requires no system delay (only bandlimited delay). Every QS impedance network is represented using lumped circuit elements such as resistors, inductors, and capacitors [6]. All Brune PR impedances are minimum phase (MP), because every PR function must be MP. Thus, every Brune impedance is PR, MP, and QS [7]. Kennelly defined Zmot as the difference between the two driving-point impedances, Zmot = Zin|free -Zin|blocked. (1) Kennelly and Affel in 1915 published a second paper where Zmot is characterized as a circle passing through the origin, with its diameter depressed by a certain angle (depressed compared to the circle in undamped impedance) in the impedance plane (Z(ω) = ℜ(ω)+jℑ(ω)). Kennelly and Affel described the circle in terms of the electrical and mechanical properties of the system. However, no justification for why the Zmot (ω) subscribed a circle was given. In 1954 Hunt further discusses the circle feature, again with no rationale for this functional form. Kennelly's third paper about Zmot [8] was focused on power concept of Zmot, and introduced motional power diagram, to provide a better physical understanding of the Zmot circle (in their view, power seemed a better concept to understand the system, compared to impedance). The motional power diagram is drawn based on magnetomotive force generated by the vibration of the diaphragm in the permanent magnetic field. They explained the motional power circle by means of 'active mechanical power', which is defined as a difference between electrical power and hysteresis power. In 1921, he went on to describe some techniques to measure various acoustic constants introduced in his previous three papers [9]. Historically, this work can be viewed as the first mechanical/acoustic impedance calibration purely from the electrical measurements. This idea was extended to compute acoustic impedance based on purely electrical data [10]. Given the practical importance of this technique, it needed to be fully characterized. This method is discussed in Ramo, et al. [1]. Besides Kennelly, Wegel [11] also considered Zmot. This paper is credited by Hunt as the inspiration of Hunt's 1954 two-port matrix representation. Wegel takes account of the general theory of receiver structures using a simple schematic having four coils. As applications, he takes four different specific cases of a receiver: A simple receiver, a receiver with eddy currents in the core, a simple induction-type receiver, and an electrodynamic receiver. He describes the receiver eddy current, which decreases proportional to square root of the frequency due to the diffusion and magnetic flux into the conducting iron core. However, the author did not derive any specific formula for this phenomenon, as it was simply an experimental observation. Vanderkooy [12] has defined the eddy-current formula from Maxwell's equation. He demonstrated that the eddy-current is a secondary current flow generated from the primary magnetic field. Frederick Vinton Hunt, an acoustic engineer, was influenced by Kennelly's discoveries. Hunt published Electroacoustics in 1954, in which, he analyzed and synthesized the electroacoustic (or electromechanical) system, modeling it as a 2 × 2 matrix using scalar forms of Ampere's law and Faraday's lawc. Hunt's two-port matrix parameters can be modeled using transmission and impedance matrix methods [7]. cIt was done by distinguishing two constants for a 90˚ time phase shift and for a 90˚ spatial phase shift [17]. F = BlkI and ɸ = Blku, where F, I, ɸ, u, B, l are force, current, voltage, velocity, magnetic intensity, and length of wire, respectively. Understanding that definition of an "impedance" is a ratio between potential and flow of a field (i.e., electrical, acoustic, or mechanical fields), one can expand the definition of impedance in multi-port network systems such as electromagnetic, electroacoustic and etc. Therefore, other than the classical driving point impedance, transfer impedance (i.e., voltage over velocity, current over velocity, or force over current and etc) was naturally developed, defined, and used by researchers to critically analyze and understand multidisciplinary systems. The Zmot is a kind of transfer impedance. Retroactively applying Brune's work, therefore Zmot does not obey the PR property. However, to our knowledge, characteristics of the Zmot, such as its negative real parts and circular shape, have yet to be clarified with regard to its physical properties [13-15]. Therefore, we aim to access the Zmot given its historical, conceptual, and physical analysis, and simulation. To engage in this goal, we will reduce the complexity of the proposed BAR model [16] to its essential elements for simplifying simulation purpose. The remainder of this study is structured as follows: Section 2 introduces theoretical concepts with definitions of Zmot as defined by Kennelly and Hunt. Section 3 provides a simple electro-magnetic model. To make use of the model, the result of the simulation follows in Section 4. Finally, we explain a possible physical reasons to make the negative real parts and the circle shape of the Zmot.
In this section, we review important theories needed to understand Zmot and anti-reciprocal networks. These include the definition and essential properties of Hunt's two-port equations which are required to define Zmot. Definition of motional impedance (Zmot) Following Wegel [11], Hunt [17] modeled an electro-mechanic system into a simple 2 × 2 impedance matrix relationship. There are four two-port network parameters, Ze(s), Zm(s), Tem(s), and Tme(s): [Φ(ω)F(ω)] = [Ze(s)Tme(s)Tem(s)Zm(s)][I(ω)U(ω)], (2) Ze = ΦI when U = 0, (3) Tem = ΦU when I =0, (4) Tme = FI when U= 0, (5) Zm = FU when I = 0, (6) Where Φ, I, F, and U are signals, the voltage, current, force, and velocity in frequency domain respectively, and 's = σ + jω' is the Laplace frequency [16]. Along with Eq. 2, the two-port 'electro-mechanic' transducer equation can alternatively be represented in ABCD matrix form, as given by [Φ(ω)I(ω)] = [A(s)C(s)B(s)D(s)][F(ω)−U(ω)] (7) Conversion between Eq. 7 and Eq. 2 has following relationship [A(s)C(s)B(s)D(s)] = 1Tme(s)[Ze(s)1ΔzZm(s)] (8) Where ∆Z = ZeZm-TemTme. Note that ∆Z = ±1 for a reciprocal and an anti-reciprocal systems respectively. For the impedance matrix to exist C = 1Tme ≠ 0, which means that Tme cannot be ∞. When the force 'F' is zero in Eq. 2, ΦI = Ze+TemUI, (9a) and UI = -TmeZm. (9b) Thus the 'unconstrained' electrical input impedance is Zin|F = 0 = ΦI|F = 0 = Ze−TemTmeZm = Ze+Zmot. (10) The Zmot is defined as the difference between two mechanical boundary conditions on the electrical impedance (Zin)d: dThe electrical conditions, 'open' and 'short', are synonymous with the two mechanical terms, 'blocked' and 'free (unconstrained)', respectively. 1) Zin with freely oscillating (vibrating) mechanical side (F = 0: Unconstrained condition), 2) Zin with the mechanical system is not allowed to move (Ze, U = 0: Blocked condition, Eq. 3), Zmot = −TemTmeZm = −TemTmeYm = Zin|F = 0−Ze = Zin|F = 0−Zin|U = 0. (11) Manipulating Eq. 8, Eq. 11 becomes Zmot = ΔzZm−Ze = −Ym−Ze, (12) Since ∆Z = -1 for the anti-reciprocal system. Equation 11 and 12 reinforce that Zmot does not have to be a Minimum-Phase (MP) nor PR function as the off-diagonal elements of Hunt's matrix, Tem and Tme, can augment delay lines. Also a difference between two PR functions does not have to be a PR function. More about these postulates follow next section. Core losses in EM transducer: Eddy current loss and Hysteresis loss Starting from Maxwell's equations we derive the core losses of motors by the use of Poynting's theorem (neglecting the source charge density and displacement current terms). By manipulating one of vector identities one may calculate the energy flow (W) into the volume (V) between two time points (t1, t2), W = -∫t2t1∫V[E(t)×H(t)]⋅dSdt = ∫t2t1∫v⎡⎣⎢H(t)⋅dB(t)dt+E(t)⋅J(t)∇×H⎤⎦⎥dVdt. (13) Here the displacement current term may be ignored for a conducting medium (i.e., conductance of a material, σ ≠ 0), leaving only J(t)e. For example, one can calculate a propagation cutoff frequency of two waves (diffusion and normal) in a medium, set the two terms in equal to each other. If e ∇×H(t) = σE(t)+∂D(t)∂t0 = J(t) μσ∂H∂tdiffusion = μ∈∂2H∂t2normal, (14) μσ(jω)H = μ∈(jω)2H, (15) σ= ∈(jω), (16) the cutoff frequency (fc) is fc = σ2π∈. (17) The fc of copper, for example, is about 4300 [GHz] (σ = 5.96 × 107, εr = 250,000, ε0 = 8.854 × 10-12), meaning that the wave below this frequency is diffusive. The corresponding wave length (λc) can be calculated as λc = ccopperfc = 3×1084.3×1012250,000−−−−−−−√≈0.14μ [m], (18) Where ccopper = c0∈r√. As the material is homogeneous, Eq. 13 reduces to W = ∫V⎡⎣⎢⎢⎢⎢∫B2B1HdBHysteresis loss+1σ∫t2t1J2(t)dtEddy current loss⎤⎦⎥⎥⎥⎥ dV, (19) Where B1 and B2 are the flux densities associated with t1 and t2, respectively. When the current (i.e., current density J) is generated from the system source, it induces the magnetic field (H) which results in two losses (Eq. 13). The air and the armature share the same induced H; however, due to the permeability differences between two mediums (µ0 = 4π × 10-7 and µi ≈ 2.5 × 10-1 for the vacuum and iron, respectively), the effective magnetic densities (B = µH) of each medium are different by a factor of 2 × 105. To increase the coupling magnetic force in Eq. 20, one can choose a softer material (that means choose the material having a smaller µ) to minimize the hysteresis loss in the system. fmag = Ψ22μiA, (20) Where ψ is the magnetic flux defined as B times the effective area (A) [17]. Equation 20 comes from the Maxwell stress tensor. Notice that Eq. 20 goes as ψ2 which leads a non-linear behavior in electromagnetic systems [17,18]. In case of a BAR, ψ in Eq. 20 can be expressed as a sum of ψDC and ψAC. fmag = 12μiA⎛⎝⎜⎜Ψ2DCi)+2ΨDCΨACii)+Ψ2ACiii)⎞⎠⎟⎟, (21) Where ψDC and ψAC are fluxes due to permanent magnets and a coil respectively. The i) Term in Eq. 21 is the steady pull on the armature which results in balancing the armature at the pivot. Assuming the ψAC to be sinusoidal, the iii) Term gives rise to a second harmonic force on the system with a negligible effect. Therefore, the remaining ii) Term takes the principal role of fmag, so that the force can be linearly approximated. To investigate the hysteresis characteristic (a.k.a. B and H curve) of the BAR, its operational path may be drawn as a small lens shape on the full hysteresis loop of the ferromagnetic material (i.e., armature) [17]. By approximating this elliptic shape as a linear line, one assumes a constant permeability which may bring computational benefits. The area of the lens corresponds to a hysteresis loss term, due to the movement of the magnetically polarized domains, in response to the magnetic field (i.e., Force ∝ flux2). To reduce the eddy current loss, one can consider the skin depth (δ) of the material which is inversely exponentially proportional to the eddy current (meaning that the strength of the eddy current is stronger on the surface and exponentially weaker in the core of the material), δ = 2ωσμ−−−−√, (22) Where ω and σ are the angular frequency and conductivity of the material [12]. The strength of eddy current becomes smaller with distance from the coil and its density attenuates exponentially with the skin depth. Therefore, choosing a hard material (large in its µ) reduces the eddy current loss. This is the opposite of the case of the hysteresis loss. Based on Eq. 20, if we choose a hard material, we expect weaker coupling magnetic force across the system.
For the further application, we introduce a simple electro-mechanic network model including a semi-inductor. The goal here is to demonstrate some condition that ℜZmot < 0 based on the simplified electro-mechanic model. To suggest a simplest possible electro-magnetic model with the physical explanation, we take an example of the BAR transducer. Figure 1 explains DC status of the BAR considering only the armature and two magnets. Assuming that the armature consists of many magnetic dipoles, the initial net polarity of the armature without magnets is neutral. However, the dipoles line up from north to south once it is influenced by permanent magnets, and the armature behaves like a magnet with DC magnetic density B0. More details about BAR operation principle are illustrated in Figure 2; the BAR's behavior when I = 0 and I ≠ 0. Not shown in this figure, however, when there is a current (I ≠ 0) in the coil which surrounds the armature, there is an eddy current which does not couple to the armature's movement. When there is no current in the coil (left of Figure 2), the armature is balancedf, so all fluxes and forces are also balanced. When current flows through the coil, the armature vibrates based on the direction of the AC current. Let us take the case when the current is greater than zero as depicted in the right part of Figure 2. In this case, the flux on the upper part of the armature, ψup, becomes a sum of ψDC-upper and ψAC and the flux on the lower part of the armature, ψlow, becomes a difference of those two fluxes so that ψup is larger than ψlow in this case. Therefore, one can draw relatively greater flux loop on the upper part of the BAR than the lower part so that the magnetic force (Eq. 20) of the upper part is larger to make the armature go up. fThe name of the BAR, balanced armature receiver, is originated from this DC status of the device. Finally, Figure 3 shows the simplifying process of modeling the BAR to help conceptual understanding, starting from the picture (the right bottom picture, Knowles ED7045). Illustrative representation of the picture is shown on the left part of the Figure 3: A coil is surrounding the armature in the middle of the device. When there is a current flow (I), it generates the magnetic field (H). There is an eddy current loss when the induced H couples to the armature. A part of H penetrates into the armature based to its skin depth generating eddy currents which direction is opposite to the the conducting current. This phenomenon is independent of the permanent magnets meaning that the eddy current cannot be seen from the mechanical side of the device. The armature behaves like a magnet due to the permanent magnets sandwiching the armature as shown in Figure 1. It alters its polarity to vibrate the armature based on the direction of the AC current. During the polarity changing process, this system experiences hysteresis loss, the energy required to rotate the domains of magnetic dipoles in the armature. Note that the magnetic density of armature is defined based on the permeability of armature (B = µaH) which is much greater than the permeability of air (µ0). As the diaphragm is attached to the tip of the armature, when the armature moves, the diaphragm also moves. The simple electro-mechanic model has been reduced from the Kim and Allen's original work (Figure 4: The electro-acoustic network model, Kim and Allen [16]. Related theories are discussed in section 2.1. Left sided figure in Figure 4 shows a oversimplified two-port network containing only essential components for better and easier understanding of the physical electro-mechanic system. In this simple model, any acoustic or resistive components are eliminated. In this figure we have four components: A semi-inductor, an inductor in the electrical port, a mass in the mechanical port, and a gyrator that links two ports. The two circuits in Figure 4 represent equivalent circuits via the mobility (dual) analogy. In both, very low and high frequencies the capacitor 'm' is opened. The parallel relation of semi-inductor and inductor enables the semi-inductor's high frequency dominance Vanderkooy [12]. The mid frequency is governed by the inductor L and the capacitor m. If we ignore the semi-inductor in Figure 4, the system looks like a Helmhorltz resonator with neck mass L and barrel compliance m. Therefore, these two components act like a resonator in the system. To realize this system into a matrix form, we can use ABCD matrix cascading method which results in Eq. 23. [Φ(ω)I(ω)] = [11ks√01][10sL1][01GG0][10sm1][F(ω)−U(ω)], (23) Where L, K, G, and m are the inductance, the semi-inductance, the gyration coefficient, and the mass of the system respectively. Let's isolating the ABCD matrix part in Eq. 23 and setting L, K, G, and m to be '1' for a simple to make the algebra simple calculation, the equation is reduced to [11s√01][10s1][0110][10s1] = [11s√sss√+1][011s]. (24) Finally, the ABCD matrix of the system in Figure 4 is [Φ(ω)I(ω)] = [T1][F(ω)−U(ω)] = [A(s)C(s)B(s)D(s)][F(ω)−U(ω)] = [sss√+11+s21s√+s2s√+s][F(ω)−U(ω)], (25) Where ∆T1 = -1. Converting Eq. 25 into an impedance matrix, Z1 = [Z11Z21Z12Z22], (26) Where Z11 = sss√+1 = ss√s+s√ (≡s||s√), (27) Z12 = −1ss√+1 = -s√s+s√, (28) Z21 = 1ss√+1 = s√s+s√, (29) Z22 = 1s√+s2s√+sss√+1 = 1+s2+ss√s+s√. (30) By substituting 's' with 'jω' one can easily find that all impedances of this system (Eq. 27, 28, 29, and 30) are complex quantities, meaning that all have both real and imaginary parts in each frequency point. The results shown in Eq. 27 - Eq. 30 are a counter example that does not follow the traditional approach of a lossless LC network. In the other words, a lossy network has been realized without having a resistor in a system. We will show in the next section that this is due to existence of the semi-inductor in a system by comparing a case where the semi-inductor does not exist. Using Eq. 11, Zmot of this system can be calculated as Zmot1 = 1(ss√+1)(1s√+s2s√+s) = ss√+s+s2+2s2s√+s3. (31) For computational benefits, we can convert Eq. 31 to an admittance (Ymot) to investigate the real part of Zmot, Ymot1 = 1+(s√)−1+s+2ss√+s2 = 1+(jω−−√)−1+jω+2jωjω−−√+(jω)2 = (1−ω2−2ωω−−√2–√+ω−−√2–√ω)+j(2ωω−−√2–√−ω−−√2–√ω+ω). (32) Since ω is always greater than 0, the real part of Eq. 32 can have negative real parts if the equation satisfies (1−ω2−2ωω−−√2–√+ω−−√2–√ω) < 0. (33) For example, if we have an angular frequency ω = 1[rad/sec], Eq. 33 is satisfied (1−1−2–√−2√2 = -12√ < 0).
We can generalize if Ymot is none positive then Zmot is also not positive. In this specific example, any angular frequency (ω) which satisfies Eq. 33 can have negative resistance in Zmot. This Zmot is not a positive definite quantity, which means it does not conserve energy of the network. This mathematically supports that it is legal to have negative real parts in Zmot which is caused by the square root components in the system such as eddy-current.
Figure 5 shows the motional impedance and input impedances with both open and short circuit conditions. To help understand better, one can think the open circuit impedance when a system is demagnetized, and the short circuit condition is the systems (i.e., a transducer) free oscillation in vacuum.
In this section, we discuss interesting features of the Zmot, its circular shape (in the complex plane) and investigate the reason to make it a circle. Based on the claim that we made, the negative real part of the Zmot is due to the diffusion of the material, we briefly discuss the two kinds transformer's core losses via the basic electro-magnetic theories including the eddy current losses.
The negative real parts of Zmot Figure 6 demonstrates a case where a difference of two input impedances (Zin with different boundary conditions) goes negative. For example, taking Z1 = Z2 = 100Ω. Based on the definition of Zmot (Eq. 11), sub-tracting the open circuit impedance from the short circuit impedance results in -50Ω (Zin|Φ2 = 0 - Zin|I2 = 0 = Z1||Z2-Z1 = 50Ω-100Ω). This simplest example tells us a lot about the nature of Zmot, as well as modeling the electro-mechanic system.
Next consider an electro-mechanic system. If there is no shunt resistance (i.e., Z1) in the electrical side of the system, Zmot cannot have negative real part as we see in Figure 6. Current is split at the parallel junction. Once the divided partial current flows through the shunted component(s) to ground, the other components (across the parallel junction), cannot see the split current. The other components are not affected by this split current which does not flow inside them. Therefore, the current is loss as the other components point of view.
To apply this theory to our system explaining the negative real parts of Zmot, the one possible physical place for the shunt component loss is the eddy current, the diffusing current into magnetic core such as armature in our specific case. It has been shown experimentally since Kennelly and Pierce [19], that Zmot has negative real parts. This fact supports the view that a shunt loss in electrical side of the system must contribute to this loss (semi-inductor) when modeling the system [16].
The circular shape of Zmot
In the polar impedance plane (ℜ vs. ℑ), Zmot is a circle passing through the origin [20]. The unusual shape may be explained by the physical nature of the anti-reciprocal electromechanic systems. The left side circuit (1) in Figure 7 describes a (typical) mechanical electromechanic network. The series of a damper, a mass, and a stiffness of the system are represented as circuit components R, L, and C, respectively. The Zmot is defined as a mechanical characteristic observed on the electrical side; therefore, simulation of these three main mechanical elements on the electrical side is our main concern.
Two circuits shown in Figure 7 are functionally equivalent. Part (1) is physically intuitive due to using a gyrator, and (2) is a dual version of (1) via the mobility analogy [21]. Figure 8 simulates the two circuit cases in Figure 7; the blue line (1) Without the gyrator (purely mechanical case) and the red line (2) Decoding the gyrator using mobility method to see the mechanical behavior on the electrical input side. The upper and lower plots in the left plane represent the magnitude and phase of the input impedance, and the right polar plot shows the real and imaginary parts of the impedance.
In Figure 8, the red circle on the polar plot (Zdual) shows Zmot, which is the dual of ZM namely,
ZM = R+1sC+sL|R,L,C = 1 = 1+1jω+jω = ⎧⎩⎨⎪⎪∞1−∞ω→∞ω→1ω→−∞, (34) Zdual = 1R||sC||1sL|R,L,C = 1 = 11+jω+1jω = ⎧⎩⎨⎪⎪010ω→0ω→1ω→∞. (35) The reason we have a circle shape of Zmot is because we are observing mechanical behavior across the gyrator. Note that Fc stands for the transition frequency between C (low-frequency) and L (high-frequency) for both original and dual of magnitude and phase plots. In polar plots, when ℑZ → +∞, Z is dominated by L and in case of ℑZ → -∞, Z depends on C.
One may suggest a refined model of Zmot based on Figure 7. The only difference between the real experimental data of Zmot and the simulation in Figure 8 is the angular rotation (clockwise direction) of the circle pivoted at the origin of the circle, which will introduce the negative real part in Zmot. One way to realize this model is to add a phase delay in the system (e-jΦ(ω)) along with mechanical circuits.
Rotating the circle toward the negative real part is related to any shunt loss in the electrical part of the system.
It is interesting to review history of impedance concept in retroactive manner. Study on impedance was built over a long-time span from Ohm to Brune. Since the AC impedance concept was developed by Kennelly [2], it took 28 years to rigorously define the properties of the impedance by Brune [22].
In this study, anti-reciprocity, a unique two port networks postulate, has been investigated via a specific real-world example, the BAR model. Zmot has been simulated and its result has demonstrated that it is not a physically realizable PR impedance with respect to Brunes impedance postulates [22]. In this sense, Zmot is a transfer impedance not a driving point impedance. Moreover, we have studied the physics behind of Zmots non-PR characteristic which lies in one of the core losses embedded in an EM system; the Eddy current loss (a theoretically appreciated shunt loss in electrical side of the system) which is independent of the coupling between the current and the force in the system. The Zmot reflects a purely mechanical behavior, a resonance, when it is looked at electrical side, thus it has a circular shape. And we have shown that the nature of anti-reciprocity makes a loop in the impedance of the electro-mechanic system by applying the mobility analogy with a gyrator which represents an anti-reciprocal characteristic. In summary, one of the most unique anti-reciprocal system characteristic, Zmot has been fully and uniquely analyzed based on both the empirical and theoretical foundations for the first time since its invention by A.E. Kennelly in 1912.