Mathematical Multibody Model of a Soft Mounted Induction Motor Regarding Forced Vibrations Due to Dynamic Rotor Eccentricities Considering Electromagnetic Field Damping ()

The paper presents a mathematical multibody model of a soft mounted induction motor with sleeve bearings regarding forced vibrations caused by dynamic rotor eccentricities considering electromagnetic field damping. The multibody model contains the mass of the stator, rotor, shaft journals and bearing housings, the electromagnetic forces with respect of electromagnetic field damping, stiffness and internal (rotating) damping of the rotor, different kinds of dynamic rotor eccentricity, stiffness and damping of the bearing housings and end shields, stiffness and damping of the oil film of the sleeve bearings and stiffness and damping of the foundation. With this multibody model, the bearing housing vibrations and the relative shaft vibrations in the sleeve bearings can be derived.

Cite this paper

Werner, U. (2017) Mathematical Multibody Model of a Soft Mounted Induction Motor Regarding Forced Vibrations Due to Dynamic Rotor Eccentricities Considering Electromagnetic Field Damping. *Journal of Applied Mathematics and Physics*, **5**, 346-364. doi: 10.4236/jamp.2017.52032.

1. Introduction

Fast running induction motors with high power ratings, $({P}_{N}>1\text{MW};$ ${n}_{N}\ge 2900\text{rpm})$ are often equipped with sleeve bearings, because of the high circumferential speed of the shaft journals, and are often mounted on soft foundations (Figure 1). A soft foundation may be realized by e.g. rubber elements under the motor feet to decouple the motor from the foundation. But also a steel frame foundation can be often characterized to be soft, because of the light weight construction. Vibrations of rotating machines are often an issue [1] - [10] .

To guarantee a safe operation, the vibrations at the sleeve bearings are often monitored [1] [2] [3] . Usually the relative shaft displacements between the shaft journals and the bearing shells are measured, using induction sensors. Additionally also the bearing housing vibrations may be measured by accelerations sensors (Figure 2).

Increasing requirements in standards and specifications of electrical machines regarding vibration limits [11] [12] [13] [14] ―e.g. IEC 60034-14, ANSI/API 541, ISO 10816, ISO 7919―require high sophisticated calculation methods. In addition to the mechanical excitation―e.g. mechanical unbalance [1] [2] [3] ― also magnetic forces occur which may lead to high vibrations [4] - [10] . In industry, these magnetic forces are nowadays still considered without the electromagnetic field damping effect, when analyzing the vibrations. The aim of the paper is now to present a multibody model for a soft mounted induction motor and to present a practical way how to consider electromagnetic field damping.

2. Dynamic Rotor Eccentricity

The three most important dynamic eccentricities for induction motors―eccen- tricity of rotor mass, bent rotor deflection and magnetic eccentricity―are here considered in the paper (Figure 3) [9] [10] .

・ Eccentricity of rotor mass ${\stackrel{^}{e}}_{u}$ which is e.g. caused by residual unbalance, which remains after the balancing process.

Figure 1. Induction motor with sleeve bearings on a soft foundation.

Figure 2. Sensors at the sleeve bearing housing.

Figure 3. Dynamic rotor eccentricities.

Figure 4. Magnetic forces at the rotor due to eccentricity.

・ Bent rotor deflection $\stackrel{^}{a}$ , which is e.g. caused by thermal bending of the rotor.

・ Magnetic eccentricity ${\stackrel{^}{e}}_{m}$ , which is e.g. caused by deviation of concentricity between the inner diameter of the rotor core and the outer diameter of the rotor core. The so caused mechanical unbalance is compensated by a placed unbalance, so that the centre of rotor mass U is not displaced from the rotation axis.

3. Electromagnetic Field Damping

If the magnetic centre M of the rotor is displaced from the centre of the stator bore (Figure 4), additionally electromagnetic fields―eccentricity fields―occur [4] - [9] . These additional fields produce a radial magnetic force ${\stackrel{^}{F}}_{mr}$ in direction of the smallest air gap. If the rotor angular frequency differs to the angular frequencies of these eccentricity fields, these fields induce a voltage into the rotor cage. The so produced harmonic rotor currents create electromagnetic fields, which lower the magnitude of the origin eccentricity fields. Therefore, the radial magnetic force ${\stackrel{^}{F}}_{mr}$ is reduced and an additional magnetic force ${\stackrel{^}{F}}_{mt}$ is generated, in tangential direction [7] [8] [9] .

These electromagnetic forces act on the rotor but in opposite direction also at the stator. For forced vibration caused by dynamic rotor eccentricity the whirling angular frequency ${\omega}_{F}$ is equal to the rotational angular frequency $\Omega $ :

${\omega}_{F}=\Omega $ (1)

Referring to [7] [8] [9] , the radial electromagnetic force can be described by an electromagnetic spring constant ${c}_{md}$ and the tangential electromagnetic force by an electromagnetic damper constant ${d}_{m}$ (with ${\omega}_{F}\ne 0$ ), depending on the pole-pair number $p$ :

$p>1:\{\begin{array}{l}{c}_{md}=\frac{{c}_{m}}{2}\cdot \left({\alpha}_{p+1}+{\alpha}_{p-1}\right)\hfill \\ {d}_{m}=-\frac{1}{{\omega}_{F}}\cdot \frac{{c}_{m}}{2}\cdot \left({\delta}_{p+1}-{\delta}_{p-1}\right)\hfill \\ \text{with}:{c}_{m}=\frac{\text{\pi}\cdot R\cdot l}{2\cdot {\mu}_{0}\cdot {\delta}^{\u2033}}\cdot {\stackrel{^}{B}}_{p}^{2}\hfill \end{array}$ (2)

$p=1:\{\begin{array}{l}{c}_{md}={c}_{m}\cdot {\alpha}_{p+1}\hfill \\ {d}_{m}=-\frac{1}{{\omega}_{F}}\cdot {c}_{m}\cdot {\delta}_{p+1}\hfill \\ \text{with}:{c}_{m}=\frac{1}{2}\cdot \frac{\text{\pi}\cdot R\cdot l}{2\cdot {\mu}_{0}\cdot {\delta}^{\u2033}}\cdot {\stackrel{^}{B}}_{p}^{2}\hfill \end{array}$ (3)

The constant ${c}_{m}$ describes the magnetic spring constant, without electromagnetic field damping, $l$ the length of the core, $R$ the radius of the stator bore, ${\mu}_{0}$ the permeability of air, ${\delta}^{\u2033}$ the equivalent magnetic air gap width, ${\stackrel{^}{B}}_{p}$ the amplitude of fundamental air gap field, ${\alpha}_{p+1}$ and ${\alpha}_{p-1}$ the real parts and ${\delta}_{p+1}$ and ${\delta}_{p-1}$ the imaginary parts of the complex field damping value. For 2-pole motors $\left(p=1\right)$ the components ${\alpha}_{p-1}$ and ${\delta}_{p-1}$ do not exist, neglecting the homopolar flux. Without electromagnetic field damping, the field damping coefficients become [7] [8] [9] :

${\alpha}_{p+1}={\alpha}_{p-1}=1;{\delta}_{p+1}={\delta}_{p-1}=0$ (4)

With the ordinal number $\nu =p\pm 1$ for an eccentricity field wave, the electromagnetic field damping coefficients can be calculated as follows [7] [8] [9] :

${\alpha}_{\nu}=1-{K}_{\nu}\cdot {s}_{\nu}^{2};\text{}{\delta}_{\nu}=-{K}_{\nu}\cdot {\beta}_{\nu}\cdot {s}_{\nu}$ with:

${\beta}_{\nu}=\frac{{R}_{2,\nu}}{{\omega}_{1}\left({L}_{2h,\nu}+{L}_{2\sigma ,\nu}\right)};\text{}{K}_{\nu}=\frac{1}{{\beta}_{\nu}^{2}+{s}_{\nu}^{2}}\cdot \frac{{\xi}_{\text{Schr},\nu}^{2}\cdot {\zeta}_{K,\nu}^{2}}{1+\frac{{L}_{2\sigma ,\nu}}{{L}_{2h,\nu}}}$ (5)

${R}_{2,\nu}$ presents the resistance of a rotor bar and ring segment, ${\omega}_{1}$ the electrical stator angular frequency, ${L}_{2h,\nu}$ the main field inductance of a rotor mesh, ${L}_{2\sigma ,\nu}$ the leakage inductance of a bar and ring segment, ${\xi}_{\text{Schr},\nu}$ the screwing factor and ${\zeta}_{K,\nu}$ the coupling factor. A very important parameter is here the harmonic slip ${s}_{\nu}$ , which can be described by [7] [8] [9] :

${s}_{\nu}=\frac{\frac{{\omega}_{\nu}}{\nu}-\Omega}{\frac{{\omega}_{1}}{\nu}}\text{with}:\text{}\Omega =\frac{{\omega}_{1}}{p}\left(1-s\right)$ (6)

Here, s presents the fundamental slip of the induction motor, ${\omega}_{1}$ the electrical stator angular frequency and ${\omega}_{\nu}/\nu $ the angular frequencies of the eccentricity fields, depending on the kind of eccentricity:

・ Static eccentricity : ${\omega}_{\nu}={\omega}_{1}$ ,

・ Dynamic eccentricity as a circular forward whirl: ${\omega}_{\nu}={\omega}_{1}\pm {\omega}_{F}$ ,

・ Dynamic eccentricity as a circular backward whirl: ${\omega}_{\nu}={\omega}_{1}\mp {\omega}_{F}$ .

In order to consider electromagnetic field damping by a simple magnetic spring element ${c}_{md}$ and a simple magnetic damper element ${d}_{m}$ , the determination has to be made, that the calculation of ${c}_{md}$ and ${d}_{m}$ is here only based on circular forward orbits [9] . This definition presents the highest electromagnetic influence, when considering electromagnetic field damping. Because of the fact, that for forced vibrations due to dynamic eccentricity the whirling frequency is equal to the rotary angular frequency $\left({\omega}_{F}=\Omega \right)$ , and that only circular forward orbits are considered for calculating the magnetic spring and damper value, the harmonic slip ${s}_{\nu}$ becomes equal to the fundamental slip $s$ [7] [8] [9] :

${s}_{\nu}=s$ (7)

4. Multibody Model

The vibration model is on the one side an enhancement of the model in [9] , where only the rotor dynamic for rigid foundation is analyzed and on the other side an enhancement of the model in [10] , where no electromagnetic field damping, no rotating damping of the rotor, no damping of the bearing housing and no mass of the bearing housings and shaft journals is considered. The innovation of the presented model is now that all these influences are now united in one single multibody model. The model is a plane multibody model, which consists of two main masses, the rotor mass ${m}_{w}$ , and the stator mass ${m}_{s}$ , which has the inertia ${\theta}_{sx}$ and is concentrated in the centre of gravity S (Figure 5).

Additional masses are the mass of the shaft journal
${m}_{v}$ and the mass of the bearing housing
${m}_{b}$ . The rotor, rotating with the rotary angular frequency
$\Omega $ , presents a concentrated mass and has no inertia moments (no gyroscopic effect is considered).The movement of the shaft journal in the sleeve bearing is described by the shaft journal centre point V. The point B, which is positioned in the axial middle of the sleeve bearing shell, describes the movement of the bearing housing. The rotor mass is linked to the stator mass by the stiffness
$c$ and internal (rotating) damping
${d}_{i}$ of the rotor, the oil film stiffness matrix C_{v} and oil film damping matrix D_{v} of the sleeve bearings, which suppose to be equal for both sides, as well as the bearing house and end shield stiffness and damping matrix C_{b} and D_{b}. The stator structure can be defined to be rigid, compared to the soft foundation. The foundation stiffness matrix C_{f} and the foundation damping matrix D_{f} connect the stator feet, F_{L} (left side) and F_{R} (right side), to the ground. The foundation stiffness and damping on the right side and on the left side is identical and the foundation stiffness values c_{fy} and c_{fz} and the foundation damping values d_{fy} and d_{fz} are the values for each motor side.

Figure 5. Multibody model.

The electromagnetism is considered by the electromagnetic spring and damper matrix ${C}_{m}$ and ${D}_{m}$ , where also electromagnetic field damping is included. Excitations are all three kinds of dynamic rotor eccentricity―eccentricity of rotor mass, bent rotor deflection and magnetic eccentricity―but are not pictured in Figure 5, because of the complexity. All used coordinate systems are fixed.

5. Stiffness and Damping Coefficients

The oil film stiffness and damping coefficients c_{ij} and
${d}_{ij}\left(i,j=z,y\right)$ of the sleeve bearing can be calculated by solving the Reynolds differential equation [15] [16] , and are depending on the rotary angular frequency
$\Omega $ :

${c}_{ij}={c}_{ij}\left(\Omega \right);\text{}{d}_{ij}={d}_{ij}\left(\Omega \right)$ (8)

The stiffness of the rotor $c$ is constant. According to [3] , the internal material damping of the rotor ${d}_{i}$ is described here by the mechanical loss factor $\mathrm{tan}{\delta}_{i}$ of the rotor, depending on the whirling angular frequency ${\omega}_{F}$ , which is here identically to the rotary angular frequency $\Omega $ :

${d}_{i}\left(\Omega \right)=\frac{c\cdot \mathrm{tan}{\delta}_{i}}{\Omega}$ (9)

The same approach is used for the bearing housing with end shield and the foundation. The stiffness of the bearing housing with end shield $\left({c}_{bz};{c}_{by}\right)$ and of the foundation $\left({c}_{fz};{c}_{fy}\right)$ is constant. The damping of the bearing housing with end shield $\left({d}_{bz};{d}_{by}\right)$ and of the foundation $\left({d}_{fz};{d}_{fy}\right)$ can be again described by the mechanical loss factor of the bearing housing with end shield $\mathrm{tan}{\delta}_{b}$ and of the foundation $\mathrm{tan}{\delta}_{f}$ :

${d}_{bz}\left(\Omega \right)=\frac{{c}_{bz}\cdot \mathrm{tan}{\delta}_{b}}{\Omega};\text{}{d}_{by}\left(\Omega \right)=\frac{{c}_{by}\cdot \mathrm{tan}{\delta}_{b}}{\Omega}$ (10)

${d}_{fz}\left(\Omega \right)=\frac{{c}_{fz}\cdot \mathrm{tan}{\delta}_{f}}{\Omega};\text{}{d}_{fy}\left(\Omega \right)=\frac{{c}_{fy}\cdot \mathrm{tan}{\delta}_{f}}{\Omega}$ (11)

The electromagnetic stiffness coefficient ${c}_{md}$ and damping coefficient ${d}_{m}$ are depending on the harmonic slip ${s}_{\nu}$ , which is here equal to the fundament slip $s$ , and on the whirling angular frequency ${\omega}_{F}$ , which is here equal to the rotary angular frequency $\Omega $ . If the motor is converter driven, the angular rotor frequency $\Omega $ as well as the fundament slip $s$ may variate arbitrarily. Therefore ${c}_{md}$ and ${d}_{m}$ become:

${c}_{md}={c}_{md}\left(\Omega ,s\right);\text{}{d}_{m}={d}_{m}\left(\Omega ,s\right)$ (12)

6. Mathematical Description

6.1. Derivation of the Differential Equation System

The forces at the rotor mass, at the shaft journals, at the bearing housings and at the stator mass can be derived in the fixed coordinate systems

$\left({y}_{W},{z}_{W};{y}_{V},{z}_{V};{y}_{B},{z}_{B};{y}_{S},{z}_{S}\right)$ (Figure 6). The rotating coordinate system

$\left({y}_{rw},{z}_{rw}\right)$ in Figure 6(a) is used for transferring the rotating damping of the rotor shaft from the rotating coordinate system into the fixed coordinate system $\left({y}_{W},{z}_{W}\right)$ [9] . The fixed coordinate systems in Figures 6(a)-(d) are used for

Figure 6. Vibration system split into subsystems.

deriving the equilibrium of forces and moments, for each single system.

Because of the small displacements of the stator mass $\left({z}_{s},{y}_{s},{\phi}_{s}\right)$ related to the dimensions of the machine $\left(h,b,\Psi \right)$ , linearization is possible [10] :

${z}_{fL}={z}_{s}-{\phi}_{s}\cdot b;{z}_{fR}={z}_{s}+{\phi}_{s}\cdot b$ (13)

${y}_{fL}={y}_{fR}={y}_{s}-{\phi}_{s}\cdot h$ (14)

To derive the inhomogeneous differential equation system, each single system ―Figures 6(a)-(d)―has to be analyzed. In Figures 6(a)-(c) the equilibrium of forces in vertical direction (z-direction) and in horizontal direction(y-direction) has to be determined for each single system. In Figure 6(d) additionally to the equilibrium of forces, the equilibrium of moments at the point S has to be determined. Based on these 9 differential equations, following inhomogeneous differential equation system can be derived:

$M\cdot \stackrel{\xa8}{q}+D\cdot \stackrel{\dot{}}{q}+C\cdot q={f}_{u}+{f}_{a}+{f}_{m}$ (15)

Coordinate vector $q$ :

$q={\left[{z}_{s};{z}_{w};{y}_{s};{y}_{w};{\phi}_{s};{z}_{v};{z}_{b};{y}_{v};{y}_{b}\right]}^{\text{T}}$ (16)

Mass matrix M:

$M=\left[\begin{array}{ccccccccc}{m}_{s}& 0& 0& 0& 0& 0& 0& 0& 0\\ 0& {m}_{w}& 0& 0& 0& 0& 0& 0& 0\\ 0& 0& {m}_{s}& 0& 0& 0& 0& 0& 0\\ 0& 0& 0& {m}_{w}& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& {\Theta}_{sx}& 0& 0& 0& 0\\ 0& 0& 0& 0& 0& 2{m}_{v}& 0& 0& 0\\ 0& 0& 0& 0& 0& 0& 2{m}_{b}& 0& 0\\ 0& 0& 0& 0& 0& 0& 0& 2{m}_{v}& 0\\ 0& 0& 0& 0& 0& 0& 0& 0& 2{m}_{b}\end{array}\right]$ (17)

Damping matrix $D$ :

$\begin{array}{l}D=[\begin{array}{cccc}2\left({d}_{fz}+{d}_{bz}\right)+{d}_{m}& -{d}_{m}& 0& 0\\ -{d}_{m}& {d}_{m}+{d}_{i}& 0& 0\\ 0& 0& 2\left({d}_{fy}+{d}_{by}\right)+{d}_{m}& -{d}_{m}\\ 0& 0& -{d}_{m}& {d}_{m}+{d}_{i}\\ 0& 0& -2{d}_{fy}\cdot h& 0\\ 0& -{d}_{i}& 0& 0\\ -2{d}_{bz}& 0& 0& 0\\ 0& 0& 0& -{d}_{i}\\ 0& 0& -2{d}_{by}& 0\end{array}\\ \text{}\begin{array}{c}0\\ 0\\ -2{d}_{fy}\cdot h\\ 0\\ 2\left({d}_{fy}{h}^{2}+{d}_{fz}{b}^{2}\right)\\ 0\\ 0\\ 0\\ 0\end{array}\begin{array}{cccc}0& -2{d}_{bz}& 0& 0\\ -{d}_{i}& 0& 0& 0\\ 0& 0& 0& -2{d}_{by}\\ 0& 0& -{d}_{i}& 0\\ 0& 0& 0& 0\\ 2{d}_{zz}+{d}_{i}& -2{d}_{zz}& 2{d}_{zy}& -2{d}_{zy}\\ -2{d}_{zz}& 2\left({d}_{zz}+{d}_{bz}\right)& -2{d}_{zy}& 2{d}_{zy}\\ 2{d}_{yz}& -2{d}_{yz}& 2{d}_{yy}+{d}_{i}& -2{d}_{yy}\\ -2{d}_{yz}& 2{d}_{yz}& -2{d}_{yy}& 2\left({d}_{yy}+{d}_{by}\right)\end{array}]\end{array}$ (18)

Stiffness matrix $C$ :

$\begin{array}{l}C=[\begin{array}{cccc}2\left({c}_{fz}+{c}_{bz}\right)-{c}_{md}& {c}_{md}& 0& 0\\ {c}_{md}& c-{c}_{md}& 0& \Omega {d}_{i}\\ 0& 0& 2\left({c}_{fy}+{c}_{by}\right)-{c}_{md}& {c}_{md}\\ 0& -\Omega {d}_{i}& {c}_{md}& c-{c}_{md}\\ 0& 0& -2{c}_{fy}h& 0\\ 0& -c& 0& -\Omega {d}_{i}\\ -2{c}_{bz}& 0& 0& 0\\ 0& \Omega {d}_{i}& 0& -c\\ 0& 0& -2{c}_{by}& 0\end{array}\\ \text{}\begin{array}{ccccc}0& 0& -2{c}_{bz}& 0& 0\\ 0& -c& 0& -\Omega {d}_{i}& 0\\ -2{c}_{fy}h& 0& 0& 0& -2{c}_{by}\\ 0& \Omega {d}_{i}& 0& -c& 0\\ 2\left({c}_{fy}{h}^{2}+{c}_{fz}{b}^{2}\right)& 0& 0& 0& 0\\ 0& 2{c}_{zz}+c& -2{c}_{zz}& 2{c}_{zy}+\Omega {d}_{i}& -2{c}_{zy}\\ 0& -2{c}_{zz}& 2\left({c}_{zz}+{c}_{bz}\right)& -2{c}_{zy}& 2{c}_{zy}\\ 0& 2{c}_{yz}-\Omega {d}_{i}& -2{c}_{yz}& 2{c}_{yy}+c& -2{c}_{yy}\\ 0& -2{c}_{yz}& 2{c}_{yz}& -2{c}_{yy}& 2\left({c}_{yy}+{c}_{by}\right)\end{array}]\end{array}$ (19)

For the calculation of the forced vibrations, the complex form is used. Therefore the excitation vectors can be described as follows:

・ Mass eccentricity:

${f}_{u}={\stackrel{^}{f}}_{u}\cdot {\text{e}}^{j\cdot \left(\Omega \cdot t+{\phi}_{u}\right)}$ (20)

・ Bent rotor deflection:

${f}_{a}={\stackrel{^}{f}}_{a}\cdot {\text{e}}^{j\cdot \left(\Omega \cdot t+{\phi}_{a}\right)}$ (21)

・ Magnetic eccentricity

${f}_{m}={\stackrel{^}{f}}_{m}\cdot {\text{e}}^{j\cdot \left(\Omega \cdot t+{\phi}_{m}\right)}$ (22)

with the amplitude vectors:

${\stackrel{^}{f}}_{u}=\left[\begin{array}{c}0\\ {\stackrel{^}{e}}_{u}\cdot {m}_{w}\cdot {\Omega}^{2}\\ 0\\ -j\cdot {\stackrel{^}{e}}_{u}\cdot {m}_{w}\cdot {\Omega}^{2}\\ 0\\ 0\\ 0\\ 0\\ 0\end{array}\right];\text{}{\stackrel{^}{f}}_{a}=\left[\begin{array}{c}0\\ \stackrel{^}{a}\cdot c\\ 0\\ -j\cdot \stackrel{^}{a}\cdot c\\ 0\\ -\stackrel{^}{a}\cdot c\\ 0\\ j\cdot \stackrel{^}{a}\cdot c\\ 0\end{array}\right];\text{}{\stackrel{^}{f}}_{m}=\left[\begin{array}{c}-{\stackrel{^}{e}}_{m}\cdot {c}_{md}\\ {\stackrel{^}{e}}_{m}\cdot {c}_{md}\\ j\cdot {\stackrel{^}{e}}_{m}\cdot {c}_{md}\\ -j\cdot {\stackrel{^}{e}}_{m}\cdot {c}_{md}\\ 0\\ 0\\ 0\\ 0\\ 0\end{array}\right]$ (23)

6.2. Solution of the Differential Equation System

With the complex form for each particular excitation:

${q}_{\kappa}={\stackrel{^}{q}}_{\kappa}\cdot {\text{e}}^{j\cdot \left(\Omega t+{\phi}_{\kappa}\right)};\kappa =u,a,m$ (24)

the complex amplitude vector for each single excitation can be calculated by:

$\begin{array}{l}{\stackrel{^}{q}}_{\kappa}={\left[-M\cdot {\Omega}^{2}+D\cdot j\cdot \Omega +C\right]}^{-1}\cdot {\stackrel{^}{f}}_{\kappa};\text{\hspace{0.17em}}\text{with:}\\ {\stackrel{^}{q}}_{\kappa}=\left[\begin{array}{c}\begin{array}{c}{\stackrel{^}{z}}_{s,\kappa}\\ {\stackrel{^}{z}}_{w,\kappa}\\ {\stackrel{^}{y}}_{s,\kappa}\end{array}\\ {\stackrel{^}{y}}_{w,\kappa}\\ {\stackrel{^}{\phi}}_{s,\kappa}\\ {\stackrel{^}{z}}_{v,\kappa}\\ {\stackrel{^}{z}}_{b,\kappa}\\ {\stackrel{^}{y}}_{v,\kappa}\\ {\stackrel{^}{y}}_{b,\kappa}\end{array}\right]=[\begin{array}{c}\begin{array}{c}\left|{\stackrel{^}{z}}_{s,\kappa}\right|\cdot {\text{e}}^{j\cdot {\alpha}_{{z}_{s,\kappa}}}\\ \left|{\stackrel{^}{z}}_{w,\kappa}\right|\cdot {\text{e}}^{j\cdot {\alpha}_{{z}_{w,\kappa}}}\\ \left|{\stackrel{^}{y}}_{s,\kappa}\right|\cdot {\text{e}}^{j\cdot {\alpha}_{{y}_{s,\kappa}}}\end{array}\\ \left|{\stackrel{^}{y}}_{w,\kappa}\right|\cdot {\text{e}}^{j\cdot {\alpha}_{{y}_{w,\kappa}}}\\ \left|{\stackrel{^}{\phi}}_{s,\kappa}\right|\cdot {\text{e}}^{j\cdot {\alpha}_{{\phi}_{s,\kappa}}}\\ \left|{\stackrel{^}{z}}_{v,\kappa}\right|\cdot {\text{e}}^{j\cdot {\alpha}_{{z}_{v,\kappa}}}\\ \left|{\stackrel{^}{z}}_{b,\kappa}\right|\cdot {\text{e}}^{j\cdot {\alpha}_{{z}_{b,\kappa}}}\\ \left|{\stackrel{^}{y}}_{v,\kappa}\right|\cdot {\text{e}}^{j\cdot {\alpha}_{{y}_{v,\kappa}}}\end{array}\end{array}$