# Multi-Field-Coupled Model and Solution of Active Electronically Scanned Array Antenna Based on Model Reconstruction

Authors: Dan Han; Jin Huang; Jinzhu Zhou; Mei Wang; Sen Xu; Haitao Li; Shen Li
Journal: International Journal of Antennas and Propagation
 (2018)

---

## Abstract

Active electronically scanned array antenna (AESA antenna) is capable of controlling the radiation pattern by controlling the feeding phase of the radiating elements. It has good performance and plays an important role in radar systems. With the development of AESA antenna towards high-frequency bands and high-density arrays, the structural-electromagnetic-thermal (SET) coupling becomes increasingly significant. It seriously restricts the realization of high performances of AESA antennas. However, the previously reported theoretical multi-field-coupled model for the coupling problem ignores the effect of the deformations of the feed system and array elements on the electrical performance. It only considers the positional deviations of the array elements in the coupling field. As a result, the accuracy of the numerical solution by the theoretical model is reduced. To overcome the above problems, this paper first establishes the field-circuit coupling model by introducing the deformation errors of the feed system into the existing theoretical model. Secondly, this paper proposes a new numerical solution for the multi-field-coupled problem of AESA antennas based on model reconstruction. And the model reconstruction includes the following: the NURBS (nonuniform rational B-spline) surface fitting algorithm that completes the mapping from finite element models to geometric models by the surface equations established by the node information and the local model reconstruction algorithm that determines the local geometric models by the positions and the directions. The NURBS surface fitting algorithm guarantees the accuracy of both the positions and shapes of array elements. The local model reconstruction algorithm ensures the accuracy of the amplitudes and phases of feed connectors. Finally, the numerical solution was applied to the 32-element AESA antenna and the simulations are close to the measurements.


---

## Body

## 1. Introduction

Active electronically scanned array antenna (AESA antenna) has high reliability, multiple functions, strong detection capability, and good stealth capability [1, 2]. In the past several decades, AESA antennas played an important role in civilian operations and the military arena [3]. To obtain better performance, the AESA antenna is developing towards higher frequency and higher density. The technological developments of meter wave and centimeter wave detection radars are relatively mature. Millimeter waves are widely used in communication systems. The frequency of high-speed and short-distance transmission equipment can reach up to dozens of gigahertz (GHz). The research of higher terahertz application has already begun. The frequency band of radio waves observed in radio astronomy can reach up to hundreds of gigahertz (GHz) [4]. Moreover, the assembly density of AESA antennas continues to increase with the transition from two-dimensional assembly to three-dimensional assembly. Accordingly, the volume of equipment must become increasingly small. For example, the size of radio frequency systems has been reduced from 0.03 m3 to 0.001 m3 [4]. Thus, the structural-electromagnetic-thermal (SET) coupling becomes more and more significant. It causes the gain loss, sidelobe level upgrade, and inaccurate beam pointing [5, 6], thus restricting the development of high-performance array antennas. Therefore, the multifield coupling is a crucial problem that needs to be solved [7].There are two major types of numerical solutions previously reported for the multifield coupling of AESA antennas. One is by the theoretical multi-field-coupled model [8, 9]. However, the model ignores the effect of the deformations of the feed system and array elements on the electrical performance. It only considers the positional deviations of the array elements in the coupling field [9], thus reducing the accuracy of the numerical solution. The other is by simulations. However, between different software in coupling fields, the data models are heterogeneous and the mesh of the finite element model does not match because of different analysis purposes and tools. To overcome the above problems in simulation, two approaches came into being. One approach is importing the heterogeneous mesh into the simulation software directly [10, 11]. However, the imported mesh must be simplified, refined, and homogenized to meet the analysis requirements and the process is complicated [12]. The other approach is fitting a new solid model from the heterogeneous mesh and then remeshing on the new model [13–15]. However, new errors are introduced in the fitting process and the accuracy of the result of the electromagnetic analysis is reduced [12].For the above problems, this paper first establishes the field-circuit coupling model by introducing the deformation errors of the feed system into the existing theoretical model. Secondly, this paper proposes a new simulation method based on model reconstruction that precisely reflects both the positional deviations and self-deformations of the array elements. One of the key benefits of the method is passing the solid model between different simulation tools in coupling fields and remesh on the new model, which avoids the tedious meshing process. Moreover, the NURBS (nonuniform rational B-spline) surface fitting algorithm guarantees the accuracy of both the positions and shapes of array elements. And the local model reconstruction algorithm ensures the accuracy of the amplitudes and phases of feed connectors. Thus, the accuracy of the electromagnetic analysis is improved. Figure1 shows the structure of the AESA antenna. And the antenna array and feed network that this paper concerns about are shown in the figure.Figure 1
Structural configuration of the AESA antenna.


## 2. Multi-Field-Coupled Model and Solution of AESA Antennas

The multifield coupling of the AESA antenna is primarily the sequential coupling among the temperature field, the structural displacement field, and the electromagnetic field. As shown in Figure2, structural deformations are caused by environmental loads and changes of structural parameters, also known as structural displacement field (path ①). Secondly, the structural displacement field affects the temperature field and the electromagnetic field through the above deformations (path ②). Finally, the influence of the temperature field has two forms: one is the impact on the performance of electromagnetic devices (such as T/R components) (path ③) and the other is the impact on the electromagnetic field by the changes of the structural displacement field (path ④) [16].Figure 2
Three-field coupling process of the AESA antenna.


### 2.1. Field-Circuit Coupling Model

For the above three-field coupling process of AESA antennas, in the traditional three-field coupling model [8, 9], as shown in (1), the excitation source Imn′ by the feed system was ideally processed, often by applying a unit excitation source or an excitation source that satisfies the Taylor distribution. However, in the coupling field, the feed system also deforms. Moreover, the electrical performance is affected by the changed characteristic impedance and power ratio of the feed network, thereby affecting the excitation current. In addition, in the feed network, the λ/4 impedance converter has a narrow line width because of the high transmission impedance. And in order to reduce the coupling effect, the spacing of its lines is large. Therefore, in the coupling field, the sensitivity of deformation error of the λ/4 impedance converter is large. In this paper, the relation between the length variation Δl of the λ/4 impedance converter in the feed network and the excitation current Imn′ is deduced. And this relation is introduced into the existing coupling model to establish the field-circuit coupling model. Figure 3 shows the spatial relation of the AESA antenna.
(1)Eθ,ϕ=∑m=0M−1∑n=0N−1Emn′θ,ϕImn′×expjkmdx+∑i=0mΔxijδ,Tfxθ,ϕ+ndy+∑j=0nΔyijδ,Tfyθ,ϕ+Δzijδ,Tfzθ+jSmnT,where Eθ,ϕ is the pattern function, Emn′θ,ϕ is the pattern function of the radiating elements in an antenna array, Imn′ is the amplitude of excitation current, T is the structure temperature, δβ1,β2,…,βR is the structural displacement, βii=1,2,…,R are the structural design variables, fxθ,ϕ, fyθ,ϕ, and fzθ are functions of element position and direction determined by the array arrangement form, such as the hexagonal AESA antenna and rectangular AESA antenna, Δxijδ,T, Δyijδ,T, and Δzijδ,T are the deviation values of the element position determined by T and δ, SmnT is the array phase difference controlled by the phase shifter affected by T, and dx and dy are the intervals of the radiating elements.Figure 3
Spatial relation of the AESA antenna.Figure4(a) shows the circuit parameters of the power divider in the feed network. Next, the relation between Δl and Imn′ is deduced from two aspects. First, the deformation of the microstrip causes impedance mismatching of the power divider, which results in the nonzero reflection coefficient at the power input end Pin.
(2)Pin=P01−Γin2.Figure 4
Diagram of the power divider.
(a)
Diagram of circuit parameters
(b)
Three-port networkSecondly, the deformation of the microstrip causes the uneven distribution of the power on each branch, assuming that the input power on one of the branchesPinB=KPin. And the K was deduced according to the impedance converter theory as
(3)K=ZinAZinA+ZinB,and thus,
(4)PinA=1−KPin.At the same time, considering that there is also the reflection between theλ/4 impedance converter and the load, thus, the load power
(5)PA=PinA1−ΓinA2,PB=PinB1−ΓinB2.Finally, according to the relation between the current and the power, the element excitation current is obtained as(6)Imn′=P0×1−ZmnΔlmnZmn′Δlmn−Z0ZmnΔlmn+Zmn′ΔlmnZmnΔlmnZmn′Δlmn+Z0ZmnΔlmn+Zmn′Δlmn2×1−Zmn′Δlmn−ZlZmn′Δlmn+Zl2×Zmn′ΔlmnZmnΔlmn+Zmn′Δlmn1/2×Zl−1/2,where P0 is the total power at the input of the power divider, Z0 is the characteristic impedance of the microstrip at the input end, Zl is the impedance of the array element, Zmn and Zmn′ (also ZinA or ZinB in the figure) are the impedance of the λ/4 impedance converter looking into the load from the input end (the arrow pointing), and Zmn′ is the impedance of the λ/4 impedance converter looking into the source from the output end (the opposite direction of the arrow). Zmn is shown as
(7)Zmn=Z2Zl+jZ2tanβl+ΔlmnZ2+jZltanβl+Δlmn.Zmn′ and Zmn′ are the same, where Z2 is the characteristic impedance of the impedance converter. And the meaning of n′ in the formula is as follows:
(8)n′=n+1,nisanoddnumber,n−1,nisanevennumber.If the above power divider is equivalent to the three-port network, as shown in Figure4(b), the cascaded microwave full matrix is
(9)Thus, the scattering parameter matrix of the three-port network is(10)S=SΙΙ+SΙΙΙβ−SΙΙΙΙ−1SΙΙΙ,where SΙΙ, SΙΙΙ, SΙΙΙ, and SΙΙΙΙ are the blocked submatrix.


### 2.2. Numerical Solution for Multifield Coupling Problem

The multifield coupling of the AESA antenna is primarily sequential coupling. Thus, the solution is performed by sequential decoupling [16]. In the traditional method by theoretical model, the actual model was simplified during the numerical solution, which only considers the positional deviations of the array elements and ignores its self-deformations. This paper performs the numerical solution by using advanced simulation software, which is based on the model reconstruction that takes into account the deformation of the array element. This approach avoids the complicated meshing process and overcomes the problem of large fitting errors in the traditional numerical simulation using simulation software. Figure 5 shows the flow chart of the numerical solution for the multifield coupling problem, and the specific steps are as follows:
(1)
The finite element model of the AESA antenna is established. And the thermal load is applied to the model in the thermal analysis software. The temperature field and the deformed model of the antenna are obtained
(2)
The deformed model in (1) is introduced into the structural analysis software using the model reconstruction method in this paper (explained in Section3), and it is applied as the external load. The structural displacement analysis is performed. And the structural displacement field and the deformed model superimposed with (1) are obtained
(3)
The model superimposed in (2) is introduced into the electromagnetic analysis software still using the model reconstruction method in this paper, and the temperature field in (1) is applied to perform the final electromagnetic performanceFigure 5
The flow chart of the numerical solution.


## 2.1. Field-Circuit Coupling Model

For the above three-field coupling process of AESA antennas, in the traditional three-field coupling model [8, 9], as shown in (1), the excitation source Imn′ by the feed system was ideally processed, often by applying a unit excitation source or an excitation source that satisfies the Taylor distribution. However, in the coupling field, the feed system also deforms. Moreover, the electrical performance is affected by the changed characteristic impedance and power ratio of the feed network, thereby affecting the excitation current. In addition, in the feed network, the λ/4 impedance converter has a narrow line width because of the high transmission impedance. And in order to reduce the coupling effect, the spacing of its lines is large. Therefore, in the coupling field, the sensitivity of deformation error of the λ/4 impedance converter is large. In this paper, the relation between the length variation Δl of the λ/4 impedance converter in the feed network and the excitation current Imn′ is deduced. And this relation is introduced into the existing coupling model to establish the field-circuit coupling model. Figure 3 shows the spatial relation of the AESA antenna.
(1)Eθ,ϕ=∑m=0M−1∑n=0N−1Emn′θ,ϕImn′×expjkmdx+∑i=0mΔxijδ,Tfxθ,ϕ+ndy+∑j=0nΔyijδ,Tfyθ,ϕ+Δzijδ,Tfzθ+jSmnT,where Eθ,ϕ is the pattern function, Emn′θ,ϕ is the pattern function of the radiating elements in an antenna array, Imn′ is the amplitude of excitation current, T is the structure temperature, δβ1,β2,…,βR is the structural displacement, βii=1,2,…,R are the structural design variables, fxθ,ϕ, fyθ,ϕ, and fzθ are functions of element position and direction determined by the array arrangement form, such as the hexagonal AESA antenna and rectangular AESA antenna, Δxijδ,T, Δyijδ,T, and Δzijδ,T are the deviation values of the element position determined by T and δ, SmnT is the array phase difference controlled by the phase shifter affected by T, and dx and dy are the intervals of the radiating elements.Figure 3
Spatial relation of the AESA antenna.Figure4(a) shows the circuit parameters of the power divider in the feed network. Next, the relation between Δl and Imn′ is deduced from two aspects. First, the deformation of the microstrip causes impedance mismatching of the power divider, which results in the nonzero reflection coefficient at the power input end Pin.
(2)Pin=P01−Γin2.Figure 4
Diagram of the power divider.
(a)
Diagram of circuit parameters
(b)
Three-port networkSecondly, the deformation of the microstrip causes the uneven distribution of the power on each branch, assuming that the input power on one of the branchesPinB=KPin. And the K was deduced according to the impedance converter theory as
(3)K=ZinAZinA+ZinB,and thus,
(4)PinA=1−KPin.At the same time, considering that there is also the reflection between theλ/4 impedance converter and the load, thus, the load power
(5)PA=PinA1−ΓinA2,PB=PinB1−ΓinB2.Finally, according to the relation between the current and the power, the element excitation current is obtained as(6)Imn′=P0×1−ZmnΔlmnZmn′Δlmn−Z0ZmnΔlmn+Zmn′ΔlmnZmnΔlmnZmn′Δlmn+Z0ZmnΔlmn+Zmn′Δlmn2×1−Zmn′Δlmn−ZlZmn′Δlmn+Zl2×Zmn′ΔlmnZmnΔlmn+Zmn′Δlmn1/2×Zl−1/2,where P0 is the total power at the input of the power divider, Z0 is the characteristic impedance of the microstrip at the input end, Zl is the impedance of the array element, Zmn and Zmn′ (also ZinA or ZinB in the figure) are the impedance of the λ/4 impedance converter looking into the load from the input end (the arrow pointing), and Zmn′ is the impedance of the λ/4 impedance converter looking into the source from the output end (the opposite direction of the arrow). Zmn is shown as
(7)Zmn=Z2Zl+jZ2tanβl+ΔlmnZ2+jZltanβl+Δlmn.Zmn′ and Zmn′ are the same, where Z2 is the characteristic impedance of the impedance converter. And the meaning of n′ in the formula is as follows:
(8)n′=n+1,nisanoddnumber,n−1,nisanevennumber.If the above power divider is equivalent to the three-port network, as shown in Figure4(b), the cascaded microwave full matrix is
(9)Thus, the scattering parameter matrix of the three-port network is(10)S=SΙΙ+SΙΙΙβ−SΙΙΙΙ−1SΙΙΙ,where SΙΙ, SΙΙΙ, SΙΙΙ, and SΙΙΙΙ are the blocked submatrix.


## 2.2. Numerical Solution for Multifield Coupling Problem

The multifield coupling of the AESA antenna is primarily sequential coupling. Thus, the solution is performed by sequential decoupling [16]. In the traditional method by theoretical model, the actual model was simplified during the numerical solution, which only considers the positional deviations of the array elements and ignores its self-deformations. This paper performs the numerical solution by using advanced simulation software, which is based on the model reconstruction that takes into account the deformation of the array element. This approach avoids the complicated meshing process and overcomes the problem of large fitting errors in the traditional numerical simulation using simulation software. Figure 5 shows the flow chart of the numerical solution for the multifield coupling problem, and the specific steps are as follows:
(1)
The finite element model of the AESA antenna is established. And the thermal load is applied to the model in the thermal analysis software. The temperature field and the deformed model of the antenna are obtained
(2)
The deformed model in (1) is introduced into the structural analysis software using the model reconstruction method in this paper (explained in Section3), and it is applied as the external load. The structural displacement analysis is performed. And the structural displacement field and the deformed model superimposed with (1) are obtained
(3)
The model superimposed in (2) is introduced into the electromagnetic analysis software still using the model reconstruction method in this paper, and the temperature field in (1) is applied to perform the final electromagnetic performanceFigure 5
The flow chart of the numerical solution.


## 3. Model Reconstruction

The key to solve the multifield coupling problem of AESA antennas using the simulation software is to overcome the problems that the data models are heterogeneous and the mesh does not match caused by the different analysis tools. In this paper, the model reconstruction of antenna array includes the reconstruction of array elements and the reconstruction of SMA connectors. Figure6 shows the structure of the antenna array. And the array elements and feeds are shown in the figure. For the above, the NURBS (nonuniform Rational B-Spline) surface fitting algorithm that guarantees the accuracy of positions and shapes of array elements and the local model reconstruction algorithm that ensures the accuracy of amplitudes and phases of feed connectors are proposed. The flow chart of the reconstruction algorithm is shown in Figure 7.Figure 6
Structure of the antenna array.Figure 7
Reconstruction algorithm flow chart.


### 3.1. Array Element Reconstruction Based on NURBS Surface Fitting

The AESA antenna undergoes various forms of deformation via the different external loads during service. The NURBS surface provides the unified mathematical form for the surfaces representing arbitrary shapes, thereby facilitating the transfer for models through the surface equations and avoiding the geometric modeling errors because of the heterogeneous data models between different coupled fields. At the same time, NURBS technology has been adopted by many excellent CAD/CAM software and has extremely high accuracy for surface fitting [17], thereby improving the fitting accuracy of array element surfaces in the model reconstruction of AESA antennas. The expression of the p×q-ary NURBS surface is [18]
(11)Su,v=∑i=0m∑j=0nwi,jdi,jNi,puNj,qv∑i=0m∑j=0nwi,jNi,puNj,qv,where di,j,i=0,1,…,m;j=0,1,…,n are the control vertices that correspond to a topological rectangular array that forms a control mesh, wi,j is the weight associated with the vertices, Ni,pu and Nj,qv are the canonical B-spline basis functions determined by the vector U=u0,u1,…,um−p+1 of the direction u and the vector V=v0,v1,…,vn−q+1 of the direction v, respectively, according to the de Boor-Cox recursion formula. The recursion formula of Ni,pu is defined as
(12)Ni,0u=1,ui≤u≤ui+1;0,Ni,pu=u−uiui+p−uiNi,p−1u+ui+p+1−uui+p+1−ui+1Ni+1,p−1u.Postulating0/0=0, the recursion formula of Nj,qv is similar.The specific steps of NURBS surface fitting are as follows:(1)
Assuming that the array of the NURBS surface isp×q-ary and assuming that wi,j=1
(2)
The node information of the finite element model was extracted as the data points fitted to the NURBS surface. Parameter values were assigned to each data pointpi,j that determines the node vector U and V [19]
(3)
The cross-section curves were constructed from the cross-section data pointspi,j on the nodal vector U using the NURBS surface inverse algorithm [18], and the control vertices d¯i,ji=0,1,…,m;j=0,1,…,n−q+1 were obtained
(4)
The cross-section curves were constructed from the cross-section data pointsd¯i,j on the nodal vector V using the NURBS surface inverse algorithm, and the control vertices di,ji=0,1,…,m;j=0,1,…,n were obtained
(5)
The surface was corrected locally and the new weightwi,j∗ was calculated by the weight changing method [19]
(6)
The NURBS surface was obtained based on the calculated control pointsdi,j, the weights wi,j∗, the node vectors U and V, and the surface arrays p and qAfter obtaining the fitting surface, the surface equation is input into the next physics field in the coupling analysis through the programmed interface program and a new surface is created in the simulation software using the surface equation. In addition, the node information of the model boundary extracted from the previous physical field is used as the constraint condition of the surface equation, which is used to reconstruct the boundary of the array element models. The reconstructed array element model is described in Section3.3.


### 3.2. Local Model Reconstruction for the SMA Connector

In order to ensure the accuracy of amplitudes and phases of the reconstructed feed model (SMA connector), this paper discusses the model reconstruction algorithm from two aspects: the location and the pointing direction of SMA connectors. Those two are the important factors affecting the impedance matching between the radiating patch and the feed lines and thus affecting the feeding amplitudes and phases.Determining the location of the SMA connector, the output port of the connector is approximated as a plane circle. The maximum values and minimum values on the plane circle by the coordinatesx, y, and z of the nodes are searched, as shown in Figure 8. According to the similar triangle principle, the median value is taken as the point position of the connector. The approximate conditions of the feed port are as follows:
(1)
The SMA connector has large rigidity and the actual deformation in service is minimal
(2)
Not considering the ductility of the materialFigure 8
Determination of the location of the SMA connector.Determining the pointing direction of the SMA connector: the surface equation has been obtained in Section3.1 and written in the general form as Fx,y,z=0. Assuming that the continuity condition is satisfied at P0x0,y0,z0 and substituting the coordinates of the position of the SMA connector (also the approximate coordinates of the center of the circle), the normal direction equation of the surface is
(13)x−x0FxP0=y−y0FyP0=z−z0FzP0.The partial derivativesFx, Fy, and Fz are constructed by substituting the coefficients of the surface equation into the above formula. Moreover, the normal direction at P0x0,y0,z0 is obtained as the new pointing direction of the SMA connector. The reconstructed SMA connector model is described in Section 3.3.


### 3.3. Results and Analysis

In the traditional numerical solutions, some scholars analyzed the multifield coupling of AESA antennas by fitting the solid model in the analysis software and then remeshing on the new model. For example, the Finite Element Modeler (FEM) software is used to complete this work. Figure9 shows the antenna element model fitted by the FEM software. In the fitting process by the software, the finite element model must first be converted into the form of patch in .x_t format through the surface information of the model. Secondly, the entity model is sewn with the patch, which is a necessary step for FEM to handle the model. In this process, some mesh information and connection information between components are lost and the fitting errors are introduced. Due to the above reasons, the accuracy of the final model is reduced as follows:
(1)
As shown in the figure, the cylindrical surface of the SMA connector consists of multiple patches that affect the integrity of the model. This is because of the sparse meshing in the previous structural analysis. However, if the mesh is encrypted, the amount of calculations will be increased. When the model is complex or there are many components, the calculation efficiency will be seriously affected. At the same time, in order to take into account the fitting accuracy and the amount of calculation, the process of mesh adjustment is extremely complicated
(2)
As shown in the figure, the upper-end face of the SMA connector is recognized as nonplanar by the electromagnetic analysis software, which makes it impossible to apply the excitation current in the electromagnetic analysis software. This is because of the deformation of the antenna array
(3)
As shown in the figure, the lower-end face of the SMA connector and the substrate cannot be completely attached to each other, which affects the feeding effectFigure 9
Antenna element model fitted by the FEM software.In the model reconstruction of this paper, as long as the location and the pointing direction of the SMA connector are accurately determined according to Section3.2, the standard cylinder in the electromagnetic analysis software can be used to complete the reconstruction of the SMA connector model. Therefore, the model is complete and regular and its lower-end surface is plane. In addition, the substrate and the SMA connector need to be subtracted by the Boolean subtraction operation in the analysis software. After the subtraction, the curvature of each point on the upper-end surface of the SMA connector coincides with the curvature of each point on the substrate and a perfect fit between the two is guaranteed. At the same time, the NURBS surface fitting algorithm ensures the accuracy of the array element model. Figure 10(a) shows the antenna element model generated by the model reconstruction method in this paper, and Figure 10(b) shows the side view.Figure 10
Antenna element model generated by the model reconstruction method.
(a)
Element model generated by the model reconstruction method
(b)
Side viewThe results of the electrical performance obtained by the FEM fitting method and the model reconstruction method were compared with the experimental test results, respectively. Figure11(a) shows the radiation pattern and Figure 11(b) shows the return loss. The results show that the model reconstruction method has higher analysis accuracy.Figure 11
Electrical performance results.
(a)
Radiation pattern
(b)
Return loss


## 3.1. Array Element Reconstruction Based on NURBS Surface Fitting

The AESA antenna undergoes various forms of deformation via the different external loads during service. The NURBS surface provides the unified mathematical form for the surfaces representing arbitrary shapes, thereby facilitating the transfer for models through the surface equations and avoiding the geometric modeling errors because of the heterogeneous data models between different coupled fields. At the same time, NURBS technology has been adopted by many excellent CAD/CAM software and has extremely high accuracy for surface fitting [17], thereby improving the fitting accuracy of array element surfaces in the model reconstruction of AESA antennas. The expression of the p×q-ary NURBS surface is [18]
(11)Su,v=∑i=0m∑j=0nwi,jdi,jNi,puNj,qv∑i=0m∑j=0nwi,jNi,puNj,qv,where di,j,i=0,1,…,m;j=0,1,…,n are the control vertices that correspond to a topological rectangular array that forms a control mesh, wi,j is the weight associated with the vertices, Ni,pu and Nj,qv are the canonical B-spline basis functions determined by the vector U=u0,u1,…,um−p+1 of the direction u and the vector V=v0,v1,…,vn−q+1 of the direction v, respectively, according to the de Boor-Cox recursion formula. The recursion formula of Ni,pu is defined as
(12)Ni,0u=1,ui≤u≤ui+1;0,Ni,pu=u−uiui+p−uiNi,p−1u+ui+p+1−uui+p+1−ui+1Ni+1,p−1u.Postulating0/0=0, the recursion formula of Nj,qv is similar.The specific steps of NURBS surface fitting are as follows:(1)
Assuming that the array of the NURBS surface isp×q-ary and assuming that wi,j=1
(2)
The node information of the finite element model was extracted as the data points fitted to the NURBS surface. Parameter values were assigned to each data pointpi,j that determines the node vector U and V [19]
(3)
The cross-section curves were constructed from the cross-section data pointspi,j on the nodal vector U using the NURBS surface inverse algorithm [18], and the control vertices d¯i,ji=0,1,…,m;j=0,1,…,n−q+1 were obtained
(4)
The cross-section curves were constructed from the cross-section data pointsd¯i,j on the nodal vector V using the NURBS surface inverse algorithm, and the control vertices di,ji=0,1,…,m;j=0,1,…,n were obtained
(5)
The surface was corrected locally and the new weightwi,j∗ was calculated by the weight changing method [19]
(6)
The NURBS surface was obtained based on the calculated control pointsdi,j, the weights wi,j∗, the node vectors U and V, and the surface arrays p and qAfter obtaining the fitting surface, the surface equation is input into the next physics field in the coupling analysis through the programmed interface program and a new surface is created in the simulation software using the surface equation. In addition, the node information of the model boundary extracted from the previous physical field is used as the constraint condition of the surface equation, which is used to reconstruct the boundary of the array element models. The reconstructed array element model is described in Section3.3.


## 3.2. Local Model Reconstruction for the SMA Connector

In order to ensure the accuracy of amplitudes and phases of the reconstructed feed model (SMA connector), this paper discusses the model reconstruction algorithm from two aspects: the location and the pointing direction of SMA connectors. Those two are the important factors affecting the impedance matching between the radiating patch and the feed lines and thus affecting the feeding amplitudes and phases.Determining the location of the SMA connector, the output port of the connector is approximated as a plane circle. The maximum values and minimum values on the plane circle by the coordinatesx, y, and z of the nodes are searched, as shown in Figure 8. According to the similar triangle principle, the median value is taken as the point position of the connector. The approximate conditions of the feed port are as follows:
(1)
The SMA connector has large rigidity and the actual deformation in service is minimal
(2)
Not considering the ductility of the materialFigure 8
Determination of the location of the SMA connector.Determining the pointing direction of the SMA connector: the surface equation has been obtained in Section3.1 and written in the general form as Fx,y,z=0. Assuming that the continuity condition is satisfied at P0x0,y0,z0 and substituting the coordinates of the position of the SMA connector (also the approximate coordinates of the center of the circle), the normal direction equation of the surface is
(13)x−x0FxP0=y−y0FyP0=z−z0FzP0.The partial derivativesFx, Fy, and Fz are constructed by substituting the coefficients of the surface equation into the above formula. Moreover, the normal direction at P0x0,y0,z0 is obtained as the new pointing direction of the SMA connector. The reconstructed SMA connector model is described in Section 3.3.


## 3.3. Results and Analysis

In the traditional numerical solutions, some scholars analyzed the multifield coupling of AESA antennas by fitting the solid model in the analysis software and then remeshing on the new model. For example, the Finite Element Modeler (FEM) software is used to complete this work. Figure9 shows the antenna element model fitted by the FEM software. In the fitting process by the software, the finite element model must first be converted into the form of patch in .x_t format through the surface information of the model. Secondly, the entity model is sewn with the patch, which is a necessary step for FEM to handle the model. In this process, some mesh information and connection information between components are lost and the fitting errors are introduced. Due to the above reasons, the accuracy of the final model is reduced as follows:
(1)
As shown in the figure, the cylindrical surface of the SMA connector consists of multiple patches that affect the integrity of the model. This is because of the sparse meshing in the previous structural analysis. However, if the mesh is encrypted, the amount of calculations will be increased. When the model is complex or there are many components, the calculation efficiency will be seriously affected. At the same time, in order to take into account the fitting accuracy and the amount of calculation, the process of mesh adjustment is extremely complicated
(2)
As shown in the figure, the upper-end face of the SMA connector is recognized as nonplanar by the electromagnetic analysis software, which makes it impossible to apply the excitation current in the electromagnetic analysis software. This is because of the deformation of the antenna array
(3)
As shown in the figure, the lower-end face of the SMA connector and the substrate cannot be completely attached to each other, which affects the feeding effectFigure 9
Antenna element model fitted by the FEM software.In the model reconstruction of this paper, as long as the location and the pointing direction of the SMA connector are accurately determined according to Section3.2, the standard cylinder in the electromagnetic analysis software can be used to complete the reconstruction of the SMA connector model. Therefore, the model is complete and regular and its lower-end surface is plane. In addition, the substrate and the SMA connector need to be subtracted by the Boolean subtraction operation in the analysis software. After the subtraction, the curvature of each point on the upper-end surface of the SMA connector coincides with the curvature of each point on the substrate and a perfect fit between the two is guaranteed. At the same time, the NURBS surface fitting algorithm ensures the accuracy of the array element model. Figure 10(a) shows the antenna element model generated by the model reconstruction method in this paper, and Figure 10(b) shows the side view.Figure 10
Antenna element model generated by the model reconstruction method.
(a)
Element model generated by the model reconstruction method
(b)
Side viewThe results of the electrical performance obtained by the FEM fitting method and the model reconstruction method were compared with the experimental test results, respectively. Figure11(a) shows the radiation pattern and Figure 11(b) shows the return loss. The results show that the model reconstruction method has higher analysis accuracy.Figure 11
Electrical performance results.
(a)
Radiation pattern
(b)
Return loss


## 4. Simulations and Experiment Results

The following is the case of multifield coupling analysis of the 32-element AESA antenna whose operating frequency is 5.8 GHz in C-band. The measurements and the comparison of simulations by different methods were shown.


### 4.1. Simulation Process Based on Model Reconstruction

#### 4.1.1. Analysis of the Structural Displacement Field and Temperature Field

As described in Section2.2, the analysis of the structural displacement field and the temperature field was performed in ANSYS. In this case, the left-end face of the antenna was fixed and the displacement of 50 mm (close to a wavelength of the 5.8 GHz antenna) was applied to the right-end face. At the same time, the temperature load was applied in the form of body load. Superimposing the force load and temperature load, the structure displacement field distribution and temperature field distribution of the antenna were obtained. And the deformation of the antenna array was extracted for the subsequent model reconstruction and electrical performance analysis, as shown in Figure 12. Figure 13 shows the material parameter testing experiment. The element types and material properties are shown in Table 1.Figure 12
Antenna array.Figure 13
Testing experiment of material parameters.Table 1
Element types and material attribute parameters.
Component
Material attribute
Element type
Elastic modulusMpa
Poisson’s ratio
Densitykg/m3
Framework
Rogers5880
SOLID185
33,500
0.38
1130
Substrate
Rogers4350
SOLID185
1200
0.36
565
Patch
Copper
SOLID185
1.08e11
0.33
8900
Feed core
Copper
SOLID185
1.08e11
0.33
8900
Feed housing
Rogers5880
SOLID185
33,500
0.38
1130


#### 4.1.2. Model Reconstruction and Electrical Performance Analysis

According to the method described in Section3, the model was reconstructed in the electromagnetic analysis software by the node information of the finite element model in the temperature and structural displacement fields. Meshing and electrical performance analysis were performed on the reconstruction model. Figure 14 shows the antenna entity model reconstructed in the electromagnetic field. The model is complete and accurate and can be directly used for simulation analysis. Figure 14(a) shows the front view, where the reconstructed model of the substrate and the patch is regular and accurate and the two are completely attached to each other. Figure 14(b) shows the rear view, where the position of the SMA connector is accurate and the direction of the coaxial line is the same as the normal of the surface. The analysis results of electrical performance are described in Section 4.2.Figure 14
Reconstructed antenna array model.
(a)
Front view of the reconstructed antenna array
(b)
Rear view of the reconstructed antenna array


### 4.2. Comparisons of Simulations and Experiment Results

The electrical performance of the antenna array was measured, as shown in Figure15. Figure 16 shows the antenna array model [20].Figure 15
Electrical performance tester.Figure 16
Antenna array model.


#### 4.2.1. Analysis of Electrical Performance without Considering Field Coupling

Without considering the field coupling, the radiation pattern obtained by the method in this paper was compared with experimental test results, as shown in Figure17. By comparison, the relative error of the gain between the two is 1.92%, the relative error of the beam pointing is 0%, and the relative error of the sidelobe level is less than 9.90%, as shown in Table 2. The experimental results illustrate the accuracy of this method.Figure 17
Radiation pattern without considering field coupling.Table 2
Comparison of measured and simulation results.
Gain (dB)
Beam direction (°)
Sidelobe (dB)
Simulation
20.44
0
6.56
6.63
Experiment
20.84
0
7.28
6.91
Relative error (%)
1.92
0
9.90
4.05


#### 4.2.2. Analysis of Electrical Performance in the Coupling Field

In the coupling field, the results of the model reconstruction method and other traditional methods were compared with the measured results as shown in Figure18. These methods were compared from the five aspects: analysis accuracy, calculation time, applicability for AESA antenna, automation capability, and coupling analysis capability. First, Table 3 shows that the model reconstruction method has the highest accuracy compared with the other traditional methods. As follows, the relative error of gain is 2.24%, the relative error of beam pointing (ratio relative to the beam width) is 2.17%, and the relative error of the sidelobe is less than 8.71%. Secondly, Table 4 shows the comparisons from other aspects. It shows that the model reconstruction method is superior to other methods: (1) The model reconstruction method avoids the tedious meshing process and saves much time in computing. (2) The model reconstruction method has good applicability for the structure of AESA antenna, which is capable to describe both connection of components and self-deformation of elements accurately. (3) The model reconstruction method has great automation capability, thereby improving the efficiency of the solution. (4) The model reconstruction method could to be used to analyze both multifield coupling and field-circuit coupling.Figure 18
Radiation pattern in the coupling field.Table 3
Comparison between solution results using different methods.
Gain (dB)
Beam direction (°)
Left sidelobe (dB)
Right sidelobe (dB)
Measured result
19.68
7.46
6.43
5.74
Model reconstruction
20.12
7.73
6.29
6.24
Relative error compared with experiment (%)
2.24
2.17
1.97
8.71
Software fitting
17.45
8.70
7.26
6.89
Relative error compared with experiment (%)
11.33
11.48
12.91
20.03
SET model
20.80
7.80
6.75
6.75
Relative error compared with experiment (%)
5.69
3.06
4.98
17.60Table 4
Comparison of different methods in multifield coupling analysis.
Model reconstruction
Software fitting
SET model
Calculation time (s)
3020
43200
3650
Applicability for AESA antenna
Describing both connection of components and self-deformation of elements
Unable to describe connection of components
Unable to describe self-deformation of elements
Automation capability
Fully automated
Unable to be automated
Unable to be automated in the process of model simplification and equivalence
Coupling analysis capability
Both multifield coupling and field-circuit coupling
Difficult to analyze multifield coupling and unable on analyzing field-circuit coupling
Unable to analyze field-circuit coupling


## 4.1. Simulation Process Based on Model Reconstruction

### 4.1.1. Analysis of the Structural Displacement Field and Temperature Field

As described in Section2.2, the analysis of the structural displacement field and the temperature field was performed in ANSYS. In this case, the left-end face of the antenna was fixed and the displacement of 50 mm (close to a wavelength of the 5.8 GHz antenna) was applied to the right-end face. At the same time, the temperature load was applied in the form of body load. Superimposing the force load and temperature load, the structure displacement field distribution and temperature field distribution of the antenna were obtained. And the deformation of the antenna array was extracted for the subsequent model reconstruction and electrical performance analysis, as shown in Figure 12. Figure 13 shows the material parameter testing experiment. The element types and material properties are shown in Table 1.Figure 12
Antenna array.Figure 13
Testing experiment of material parameters.Table 1
Element types and material attribute parameters.
Component
Material attribute
Element type
Elastic modulusMpa
Poisson’s ratio
Densitykg/m3
Framework
Rogers5880
SOLID185
33,500
0.38
1130
Substrate
Rogers4350
SOLID185
1200
0.36
565
Patch
Copper
SOLID185
1.08e11
0.33
8900
Feed core
Copper
SOLID185
1.08e11
0.33
8900
Feed housing
Rogers5880
SOLID185
33,500
0.38
1130


### 4.1.2. Model Reconstruction and Electrical Performance Analysis

According to the method described in Section3, the model was reconstructed in the electromagnetic analysis software by the node information of the finite element model in the temperature and structural displacement fields. Meshing and electrical performance analysis were performed on the reconstruction model. Figure 14 shows the antenna entity model reconstructed in the electromagnetic field. The model is complete and accurate and can be directly used for simulation analysis. Figure 14(a) shows the front view, where the reconstructed model of the substrate and the patch is regular and accurate and the two are completely attached to each other. Figure 14(b) shows the rear view, where the position of the SMA connector is accurate and the direction of the coaxial line is the same as the normal of the surface. The analysis results of electrical performance are described in Section 4.2.Figure 14
Reconstructed antenna array model.
(a)
Front view of the reconstructed antenna array
(b)
Rear view of the reconstructed antenna array


## 4.1.1. Analysis of the Structural Displacement Field and Temperature Field

As described in Section2.2, the analysis of the structural displacement field and the temperature field was performed in ANSYS. In this case, the left-end face of the antenna was fixed and the displacement of 50 mm (close to a wavelength of the 5.8 GHz antenna) was applied to the right-end face. At the same time, the temperature load was applied in the form of body load. Superimposing the force load and temperature load, the structure displacement field distribution and temperature field distribution of the antenna were obtained. And the deformation of the antenna array was extracted for the subsequent model reconstruction and electrical performance analysis, as shown in Figure 12. Figure 13 shows the material parameter testing experiment. The element types and material properties are shown in Table 1.Figure 12
Antenna array.Figure 13
Testing experiment of material parameters.Table 1
Element types and material attribute parameters.
Component
Material attribute
Element type
Elastic modulusMpa
Poisson’s ratio
Densitykg/m3
Framework
Rogers5880
SOLID185
33,500
0.38
1130
Substrate
Rogers4350
SOLID185
1200
0.36
565
Patch
Copper
SOLID185
1.08e11
0.33
8900
Feed core
Copper
SOLID185
1.08e11
0.33
8900
Feed housing
Rogers5880
SOLID185
33,500
0.38
1130


## 4.1.2. Model Reconstruction and Electrical Performance Analysis

According to the method described in Section3, the model was reconstructed in the electromagnetic analysis software by the node information of the finite element model in the temperature and structural displacement fields. Meshing and electrical performance analysis were performed on the reconstruction model. Figure 14 shows the antenna entity model reconstructed in the electromagnetic field. The model is complete and accurate and can be directly used for simulation analysis. Figure 14(a) shows the front view, where the reconstructed model of the substrate and the patch is regular and accurate and the two are completely attached to each other. Figure 14(b) shows the rear view, where the position of the SMA connector is accurate and the direction of the coaxial line is the same as the normal of the surface. The analysis results of electrical performance are described in Section 4.2.Figure 14
Reconstructed antenna array model.
(a)
Front view of the reconstructed antenna array
(b)
Rear view of the reconstructed antenna array


## 4.2. Comparisons of Simulations and Experiment Results

The electrical performance of the antenna array was measured, as shown in Figure15. Figure 16 shows the antenna array model [20].Figure 15
Electrical performance tester.Figure 16
Antenna array model.


### 4.2.1. Analysis of Electrical Performance without Considering Field Coupling

Without considering the field coupling, the radiation pattern obtained by the method in this paper was compared with experimental test results, as shown in Figure17. By comparison, the relative error of the gain between the two is 1.92%, the relative error of the beam pointing is 0%, and the relative error of the sidelobe level is less than 9.90%, as shown in Table 2. The experimental results illustrate the accuracy of this method.Figure 17
Radiation pattern without considering field coupling.Table 2
Comparison of measured and simulation results.
Gain (dB)
Beam direction (°)
Sidelobe (dB)
Simulation
20.44
0
6.56
6.63
Experiment
20.84
0
7.28
6.91
Relative error (%)
1.92
0
9.90
4.05


### 4.2.2. Analysis of Electrical Performance in the Coupling Field

In the coupling field, the results of the model reconstruction method and other traditional methods were compared with the measured results as shown in Figure18. These methods were compared from the five aspects: analysis accuracy, calculation time, applicability for AESA antenna, automation capability, and coupling analysis capability. First, Table 3 shows that the model reconstruction method has the highest accuracy compared with the other traditional methods. As follows, the relative error of gain is 2.24%, the relative error of beam pointing (ratio relative to the beam width) is 2.17%, and the relative error of the sidelobe is less than 8.71%. Secondly, Table 4 shows the comparisons from other aspects. It shows that the model reconstruction method is superior to other methods: (1) The model reconstruction method avoids the tedious meshing process and saves much time in computing. (2) The model reconstruction method has good applicability for the structure of AESA antenna, which is capable to describe both connection of components and self-deformation of elements accurately. (3) The model reconstruction method has great automation capability, thereby improving the efficiency of the solution. (4) The model reconstruction method could to be used to analyze both multifield coupling and field-circuit coupling.Figure 18
Radiation pattern in the coupling field.Table 3
Comparison between solution results using different methods.
Gain (dB)
Beam direction (°)
Left sidelobe (dB)
Right sidelobe (dB)
Measured result
19.68
7.46
6.43
5.74
Model reconstruction
20.12
7.73
6.29
6.24
Relative error compared with experiment (%)
2.24
2.17
1.97
8.71
Software fitting
17.45
8.70
7.26
6.89
Relative error compared with experiment (%)
11.33
11.48
12.91
20.03
SET model
20.80
7.80
6.75
6.75
Relative error compared with experiment (%)
5.69
3.06
4.98
17.60Table 4
Comparison of different methods in multifield coupling analysis.
Model reconstruction
Software fitting
SET model
Calculation time (s)
3020
43200
3650
Applicability for AESA antenna
Describing both connection of components and self-deformation of elements
Unable to describe connection of components
Unable to describe self-deformation of elements
Automation capability
Fully automated
Unable to be automated
Unable to be automated in the process of model simplification and equivalence
Coupling analysis capability
Both multifield coupling and field-circuit coupling
Difficult to analyze multifield coupling and unable on analyzing field-circuit coupling
Unable to analyze field-circuit coupling


## 4.2.1. Analysis of Electrical Performance without Considering Field Coupling

Without considering the field coupling, the radiation pattern obtained by the method in this paper was compared with experimental test results, as shown in Figure17. By comparison, the relative error of the gain between the two is 1.92%, the relative error of the beam pointing is 0%, and the relative error of the sidelobe level is less than 9.90%, as shown in Table 2. The experimental results illustrate the accuracy of this method.Figure 17
Radiation pattern without considering field coupling.Table 2
Comparison of measured and simulation results.
Gain (dB)
Beam direction (°)
Sidelobe (dB)
Simulation
20.44
0
6.56
6.63
Experiment
20.84
0
7.28
6.91
Relative error (%)
1.92
0
9.90
4.05


## 4.2.2. Analysis of Electrical Performance in the Coupling Field

In the coupling field, the results of the model reconstruction method and other traditional methods were compared with the measured results as shown in Figure18. These methods were compared from the five aspects: analysis accuracy, calculation time, applicability for AESA antenna, automation capability, and coupling analysis capability. First, Table 3 shows that the model reconstruction method has the highest accuracy compared with the other traditional methods. As follows, the relative error of gain is 2.24%, the relative error of beam pointing (ratio relative to the beam width) is 2.17%, and the relative error of the sidelobe is less than 8.71%. Secondly, Table 4 shows the comparisons from other aspects. It shows that the model reconstruction method is superior to other methods: (1) The model reconstruction method avoids the tedious meshing process and saves much time in computing. (2) The model reconstruction method has good applicability for the structure of AESA antenna, which is capable to describe both connection of components and self-deformation of elements accurately. (3) The model reconstruction method has great automation capability, thereby improving the efficiency of the solution. (4) The model reconstruction method could to be used to analyze both multifield coupling and field-circuit coupling.Figure 18
Radiation pattern in the coupling field.Table 3
Comparison between solution results using different methods.
Gain (dB)
Beam direction (°)
Left sidelobe (dB)
Right sidelobe (dB)
Measured result
19.68
7.46
6.43
5.74
Model reconstruction
20.12
7.73
6.29
6.24
Relative error compared with experiment (%)
2.24
2.17
1.97
8.71
Software fitting
17.45
8.70
7.26
6.89
Relative error compared with experiment (%)
11.33
11.48
12.91
20.03
SET model
20.80
7.80
6.75
6.75
Relative error compared with experiment (%)
5.69
3.06
4.98
17.60Table 4
Comparison of different methods in multifield coupling analysis.
Model reconstruction
Software fitting
SET model
Calculation time (s)
3020
43200
3650
Applicability for AESA antenna
Describing both connection of components and self-deformation of elements
Unable to describe connection of components
Unable to describe self-deformation of elements
Automation capability
Fully automated
Unable to be automated
Unable to be automated in the process of model simplification and equivalence
Coupling analysis capability
Both multifield coupling and field-circuit coupling
Difficult to analyze multifield coupling and unable on analyzing field-circuit coupling
Unable to analyze field-circuit coupling


## 5. Conclusion

In this paper, the field-circuit coupling model was established by introducing the deformation errors of the feed system into the existing theoretical model. In addition, the new numerical solution method for the multi-field-coupled problem of AESA antennas based on model reconstruction was proposed. In the model reconstruction, two algorithms were presented: the NURBS (nonuniform rational B-spline) surface fitting algorithm that completes the mapping from finite element models to geometric models by the surface equations established by the node information and the local model reconstruction algorithm that determines the local geometric models by the positions and the directions. The numerical solution was applied to the 32-element AESA antenna. And the comparisons between the simulations and measured results shows the following: (1) Compared with the traditional fitting simulation software, the proposed method avoids the fussy meshing process and saves much time in computing. (2) Compared with the traditional coupled SET model, the proposed method has high accuracy. This is because the NURBS surface fitting algorithm guarantees the accuracy of both the positions and shapes of array elements and the local model reconstruction algorithm ensures the accuracy of the amplitudes and phases of feed connectors. Furthermore, the method proposed in this paper provides a new idea for the numerical solution of the multi-field-coupled problem of AESA antennas.



---

*Source: 3161928-2018-11-27.xml*