# An Approximation Mathematical Formula of Pattern Analysis for Distorted Reflector Antennas considering Surface Normal Vector Variation

Authors: Shuxin Zhang; Guigeng Yang; Yiqun Zhang
Journal: International Journal of Antennas and Propagation
 (2016)

---

## Abstract

An approximation mathematical formula to compute the distorted radiation pattern for reflector antennas is presented. In this approximation formula, besides the phase error caused by the structural deformation being added in the far field integral, the surface normal vector variation is also taken into consideration. The formula is derived by expanding the surface normal vector into a first-order Taylor series and the phase error into a second-order Taylor series. By assembling the integrals including the contributions of both surface normal vector variation and phase error, the far field electrical vector expressed as a function of structural nodal displacements is further obtained in a matrix form. Simulation results of a distorted reflector show the application of this formula.


---

## Body

## 1. Introduction

Reflector antennas are widely used in radar, telecommunication, radio astronomy, and other applications. Reflector surface distortions are inevitable in outer working conditions owing to thermal, gravitational, and dynamic effects and manufacturing tolerance [1]. With the increasing frequency, the stringent requirements on reflector surface accuracy become demanding [2]. Towards the evaluation of surface accuracy on electromagnetic performance, several tolerance analysis techniques have been proposed in scientific literature for random error [3–6] and systematic error [7–11] in centered reflectors [7–10], offset reflectors [11], and membrane reflectors [11]. By performing tolerance analysis, the reflector antenna can be designed and manufactured with considering the surface tolerance.The process of antenna design belongs to the area of multidisciplinary design, which contains two main steps of electromagnetic design and structural design. In the previous antenna design, the structural design procedure was implemented by simply reducing the surface root-mean-square (RMS) error to fulfill the requirements of antenna electromagnetic performance, which usually could not obtain the final satisfactory electromagnetic performance as observed in [10, 12]. Thus, the antenna design has been transformed from just reducing surface RMS distortion to an integrated structural-electromagnetic optimization design with a multidisciplinary optimization model [12–15]. The integrated structural-electromagnetic design concept is employed in various applications, which directly combines the electromagnetic performance gain or sidelobe levels with structural inputs other than simply reducing the surface RMS error. Some studies have been performed in the antenna design based on the integrated structural-electromagnetic concept, such as the shape control of antennas [12, 13], multidisciplinary optimization [14, 15], and electromechanical coupling analysis [9, 10]. In the integrated optimization design procedure, the antenna far field pattern will be calculated in the optimization iteration repeatedly, which requires an accurate and fast pattern analysis method for a given reflector with various deformations. An approximation radiation integral method for distorted reflector antennas was proposed in [16] using the decomposition orthogonal basis functions. The proposed method benefits the repeated calculations with allowable computation accuracy. Furthermore, referencing to the three-term Taylor series of phase error in [16], a matrix form formula between the surface nodal displacements obtained by structural finite element method (FEM) and far field electrical vector was derived in [17]. The far field pattern can be directly and easily calculated by substituting the surface nodal displacements into the matrix form, which benefits the iterative pattern analysis procedure of antenna design.As for the distorted reflector antennas, phase error caused by the surface nodal displacement is usually added in the integral to evaluate its effect on far field radiation pattern [9–11]. The phase error is introduced based on the assumption that the reflector is located at the far field zone of the primary feed, and the amplitude variation of incident magnetic field is neglected. Strictly speaking, the surface distortion introduces not only phase error but also surface normal vector variation. In the previous study, the surface normal variation is neglected, which is based on the fact that, for small-amplitude, smoothly varying errors, the surface normal vector will not deviate substantially from the surface of an undistorted reflector [16]. The questions appear as what percentage of the surface normal variation contributes in the far field integral for distorted reflectors so as to be neglected, and what is the final expression for the matrix form in [17] when considering both the surface normal variation and phase error. The aim of this study, which is a continuation of the work in [17], is to further develop the matrix form formula considering both the phase error and the surface normal vector variation caused by surface nodal displacements.In this study, the reflector surface is divided into many small elements as in structural FEM analysis, and the far field pattern is obtained by superposition of each element’s far field contribution. The surface normal vector variation caused by surface nodal displacement is expanded into a first-order Taylor series during the implementation, and the shape functions in [17] are also employed to express the integration point displacements as functions of nodal displacements. Both the phase error and the surface normal variation are included in these integrals, and the matrix form in [17] by assembling the integrals on elements is thus extended.The main difference between this formula and the one in [17] lies in the normal vector variation. The approximated formula in [17] just takes the phase error caused by surface nodal displacements into account, while it fails to consider the surface normal vector variation. The derived formula in this study will take both the phase error and normal vector variation into account theoretically, which can be considered as an extended formula for the one in [17].The study is organized as follows: Section2 outlines the derivation of the approximation mathematical formula considering both the surface normal variation and phase error, a distorted reflector is calculated to show the proportion of the surface normal variation contribution in Section 3, and the conclusion is drawn in Section 4.


## 2. The Approximation Mathematical Formula

The radiation integral for Physical Optics (PO) formula of reflector antennas is expressed as [16] (1)
T
→
θ
,
ϕ
=
∫
σ
2
N
→
×
H
→
r
→
′
e
j
k
r
→
′
·
r
^
d
σ
,where N
→ is the surface normal vector, H
→ is the incident magnetic field, the vector r
→
′ locates the integration point, j
=
-
1 is the imaginary unit, k
=
2
π
/
λ is the wavenumber, λ is the free-space wavelength at the working frequency, the unit vector r
^ is in the observation direction, and σ is the projected aperture region of reflector antennas.Provided that a nodal displacement from an ideal undistorted reflector is defined asε
→
z
′
=
Δ
z
z
^ in the z-direction shown in Figure 1, Δ
z is a scalar function that defines the nodal deformation. Strictly speaking, in the distorted antenna analysis the surface deformation introduces not only the phase error but also the surface normal vector variation. Thus, the far field integral becomes (2)
T
→
θ
,
ϕ
=
∫
σ
2
N
→
+
Δ
N
→
×
H
→
r
→
′
e
j
k
r
→
′
·
r
^
e
j
k
Δ
z
cos
⁡
θ
s
+
cos
⁡
θ
d
σ
,where Δ
N
→ is the normal vector variation, θ
s is the undistorted spherical component in the feed coordinate system as shown in Figure 1, and θ defines the observation direction. Both N
→ and H
→ are in the ideal undistorted nominal state. As described in Figure 1, the solid line represents the ideal surface and the dashed line is for the distorted surface. A surface nodal displacement is shown in Figure 1. The global coordinate system (
x
,
y
,
z
) is defined as the reflector unprimed coordinate system and (
x
s
,
y
s
,
z
s
) is for the feed coordinate system. Shown in Figure 1, r
→ denotes the observation point and r
→
′ is the surface point vector in global coordinate system.Figure 1
A distorted reflector.The surface nodal displacements shown in Figure1 will introduce phase error and surface normal vector variation in far field radiation integral as illustrated in (2). In [16], a three-term Taylor series approximation about phase error is adopted to express the phase error as a second-order polynomial expansion of nodal displacements, and its accuracy can be sufficient enough for errors profiles up to about 0.1
λ RMS amplitude [16]: (3)
e
j
k
Δ
z
cos
⁡
θ
s
+
cos
⁡
θ
≈
∑
i
=
0
2
1
i
!
j
k
Δ
z
cos
⁡
θ
s
+
cos
⁡
θ
i
.In this error analysis, the reflector surface is divided into a series of small triangular elements, and the far field pattern is thus calculated by superposition of each element’s far field vector contribution. In this derivation, the small elements are required to be equilateral triangles projected in the aperture plane, such that the element normal vector variation can be written as a linear combination of triangular nodal displacements in a first-order Taylor series, and the coefficients of these series depend only on the triangular side length and the nodal position with respect to the element center point. The normal vector variation in (2) can be expressed as (4)
Δ
N
→
=
∑
i
=
1
3
∂
N
→
∂
z
i
Δ
z
i
,where Δ
z
i is the surface nodal z-displacement of the ith node and ∂
N
→
/
∂
z
i is the gradient of normal vector with respect to the ith nodal displacement.Because of the similarity of the equilateral triangles, two trianglesA
B
C and A
′
B
C are employed in Figure 2 to show the gradient of surface normal vector with respect to nodal displacement, where the symbols A, B, C, and A
′ denote the triangular nodes. The gradients of triangular normal vector with respect to each nodal displacement are listed in Table 1, where l is the triangular side length.Table 1
Gradients of normal vector with respect to triangular nodal displacement.
Element
Node
Gradient with respect to nodal displacement
Value
A
B
C
A
∂
N
→
∂
z
A
(
0
,
-
2
/
3
,
0
)
l
B
∂
N
→
∂
z
B
(
1,1
/
3
,
0
)
l
C
∂
N
→
∂
z
C
(
-
1,1
/
3
,
0
)
l
A
′
B
C
A
′
∂
N
→
∂
z
A
′
(
0,2
/
3
,
0
)
l
B
∂
N
→
∂
z
B
(
1
,
-
1
/
3
,
0
)
l
C
∂
N
→
∂
z
C
(
-
1
,
-
1
/
3
,
0
)
lFigure 2
Two equilateral triangles in aperture plane.The same as in [17], the triangular element A
B
C in aperture plane shown in Figure 3 is chosen in the analysis, where the symbols A, B, and C represent the triangle nodes and P is an arbitrary point inside the element. Supposing that the nodal displacements in z-direction are labeled as Δ
z
a, Δ
z
b, and Δ
z
c for nodes A, B, and C, respectively, the displacement of an integration point P in (3) within the element can be expressed as a function of nodal displacements multiplied by the element shape functions, and(5)
Δ
z
=
ξ
a
,
ξ
b
,
ξ
c
Δ
z
a
Δ
z
b
Δ
z
c
,where ξ
a, ξ
b, and ξ
c are the element shape functions with respect to the nodes A, B, and C, respectively.Figure 3
A triangular elementA
B
C in FEM analysis [17].Substituting (3), (4), and (5) into (2) and neglecting the high-order terms, the integral on the triangular element A
B
C (labeled with e as superscript) can be expanded as (6)
T
→
e
=
T
→
0
e
+
T
→
1
e
+
T
→
2
e
+
T
→
3
e
+
T
→
4
e
,where the superscript e designates the element A
B
C, T
→
0
e is the undistorted contribution in the integral, and T
→
1
e and T
→
2
e are the linear and square terms in [17]; the first three terms in the right side are not repeated here for brevity, and (7)
T
→
3
e
=
∑
i
Δ
z
i
T
→
3
,
i
e
,
T
→
4
e
=
∑
u
,
v
Δ
z
u
Δ
z
v
T
→
4
,
u
v
e
,where Δ
z
i, Δ
z
u, and Δ
z
v are the nodal displacements of nodes i, u, and v, where i
=
a, b, c and u, v
=
a, b, and c, respectively, and (8)
T
→
3
,
i
e
=
∫
σ
A
B
C
2
∂
N
→
∂
z
i
×
H
→
r
→
′
e
j
k
r
→
′
·
r
^
d
σ
A
B
C
,
(9)
T
→
4
,
u
v
e
=
∫
σ
A
B
C
2
∂
N
→
∂
z
u
×
H
→
r
→
′
e
j
k
r
→
′
·
r
^
j
k
ξ
v
cos
⁡
θ
s
+
cos
⁡
θ
d
σ
A
B
C
,where σ
A
B
C is the integral region of element A
B
C, ∂
N
→
/
∂
z
u is the gradient of surface normal vector with respect to uth nodal displacement, ξ
v are the element shape function with respect to vth node in element A
B
C. From its expression, it can be seen that T
→
4
,
u
v
e
≠
T
→
4
,
v
u
e because of the different sequences of normal vector gradient and shape function in (9).The same as in [17], the integral in (6) can be rewritten as a matrix form: (10)
Δ
T
→
e
=
h
1
e
+
h
3
e
·
Δ
z
e
+
h
2
e
+
h
4
e
·
v
e
c
Δ
z
e
·
Δ
z
e
T
,where Δ
z
e
=
{
Δ
z
a
,
Δ
z
b
,
Δ
z
c
}
T is the nodal displacements vector in the element A
B
C, the operator v
e
c
(
·
) converts a matrix into a vector, h
1
e and h
2
e are the same coefficients as in [17], and(11)
h
3
e
=
T
→
3
,
a
e
,
T
→
3
,
b
e
,
T
→
3
,
c
e
,
h
4
e
=
T
→
4
,
a
a
e
,
T
→
4
,
a
b
e
,
T
→
4
,
a
c
e
,
T
→
4
,
b
a
e
,
T
→
4
,
b
b
e
,
T
→
4
,
b
c
e
,
T
→
4
,
c
a
e
,
T
→
4
,
c
b
e
,
T
→
4
,
c
c
e
.Assuming that there arem elements in the FEM analysis and n observation points in the far field requirements, the extended final matrix form can be expressed by assembling the whole elements integrals as (12)
Δ
E
→
=
-
j
k
η
e
-
j
k
r
4
π
r
I
-
-
-
r
^
r
^
·
H
1
+
H
3
Δ
z
+
H
2
+
H
4
Δ
z
2
,where Δ
z is the z-direction displacements of all nodes, Δ
z
2 is the product term between each two nodal displacements in one element including square terms of each nodal displacement, and(13)
Δ
E
→
=
Δ
E
→
θ
1
,
ϕ
1
,
Δ
E
→
θ
2
,
ϕ
2
,
…
,
Δ
E
→
θ
n
,
ϕ
n
T
,
H
1
=
A
e
=
1
m
h
1
e
,
H
2
=
A
e
=
1
m
h
2
e
,
H
3
=
A
e
=
1
m
h
3
e
,
H
4
=
A
e
=
1
m
h
4
e
,
Δ
z
=
A
e
=
1
m
Δ
z
e
,
Δ
z
2
=
A
e
=
1
m
v
e
c
Δ
z
e
·
Δ
z
e
T
,where Δ
E
→
=
E
→
-
E
→
0, E
→
0 is the undistorted electrical far field vector, and A
e
=
1
m is the finite element assembly operator [17, 18].The matrix form in [17] is thus extended in (12) considering both the phase error and the surface normal vector variation. The same as the properties of the matrix form in [17], the extended formula can also benefit the iterative optimization procedure, and once the matrices H
1, H
2, H
3, and H
4 and the undistorted electrical far field E
→
0 are determined, the distorted radiation pattern will be easily obtained by a simple matrix operation with the stored data.The extended formula can also be employed in the application of pattern analysis for different surface-error profiles as addressed in [17]. For example, in the integrated structural-electromagnetic optimization design [14, 15], the antenna surface is determined through different various loads and structural design variables, which makes the repeated pattern analysis required for a given reflector with varied error profiles. It is also needed to analyze the electromagnetic performance of a reflector with different surface-error profiles as in the contoured beam pattern synthesis. Besides, this formula can also be applied into the uncertainty bounds analysis of reflector antennas with surface random errors. The main application of this approximated formula lies in the requirement of repeated pattern analysis, namely, pattern analysis for a given reflector with various surface-error distributions. Once the prestored data is computed, the antenna pattern can be easily calculated by recalling the stored matrices and performing a simple matrix operation for a given reflector.


## 3. Numerical Assessments

The approximation method in [16] expanded the phase error caused by the surface displacements into three-term Taylor series, and the second-order approximation showed very good agreement with the exact computation. To find out what percentage of surface normal vector variation is taken in the far field integral for error profiles up to about 0.1
λ RMS amplitude, a simulation is performed in this study. The contributions of the matrices H
3 and H
4 in (12) are computed with the former matrices H
1 and H
2. For the distorted reflector in [16, 17] of which the diameter D
=
100
λ, the focal length F
=
100
λ is adopted in this simulation. The feed parameters are also the same as the ones in [16, 17] with Q
x
=
Q
y
=
8 to radiate a linearly polarized wave. The equilateral triangular side length projected in the aperture is λ. The three deformations in [16, 17] are also employed in this simulation, and their functions are (14)
Δ
z
ρ
′
,
ϕ
′
=
0.05
λ
sin
⁡
2
π
ρ
′
2
50
λ
2
,
(15)
Δ
z
ρ
′
,
ϕ
′
=
0.05
λ
cos
⁡
10
π
ρ
′
2
50
λ
2
,
(16)
Δ
z
ρ
′
,
ϕ
′
=
0.05
λ
cos
⁡
5
ϕ
′
,where the variables ρ
′ and ϕ
′ are the same as in [17] which describes the polar coordinate components in the projected aperture plane of the reflector.Figures4–6 show the far field patterns under three different deformations for the distorted reflector, respectively. Each figure has the patterns for the first-order approximation in [17] (dashed line), the extended first-order approximation (solid line), the second-order approximation in [17] (dot line), and the extended second-order approximation (the marked line) in x
z plane. Table 2 lists the maximum directivities of the patterns. The simulation results of the extended first-order and second-order approximations show very close agreement with the pervious study in [17], and the maximum difference is in the order of 0.001 dB, which indicates that the surface normal variation takes little percentage in the distorted reflector pattern analysis and it can be neglected in the approximation for errors profiles up to about 0.1
λ RMS amplitude. It can be understood from this that, in the equivalent induced current density (J
→
e
f
f
=
2
N
→
×
H
→) of the cross product of surface normal vector and the incident magnetic field, the phase of induced current density is same as the incident magnetic field while its amplitude variation can be neglected when distortion happens.Table 2
Maximum directivity of each pattern.
Unit: dB
Deformation (14)
Deformation (15)
Deformation (16)
First-order
49.072
49.023
49.023
Extended first-order
49.071
49.026
49.023
Second-order
48.273
48.269
48.214
Extended second-order
48.272
48.271
48.214Figure 4
Far field pattern comparisons of the extended formula with the previous one under deformation (14).Figure 5
Far field pattern comparisons of the extended formula with the previous one under deformation (15).Figure 6
Far field pattern comparisons of the extended formula with the previous one under deformation (16).


## 4. Conclusion

An extended approximation mathematical formula to compute the radiation pattern for distorted reflector antennas was proposed considering both the phase error and surface normal vector variation. The surface normal vector variation was added in the integral by a first-order Taylor series. The matrix form combining the surface nodal displacements and far field electrical vector was thus extended by assembling the integrals on structural elements. The extended formula also has the property of the pervious one that once the matrices are precomputed, the distorted pattern can be easily obtained by recalling the stored matrices and performing a simple matrix operation for a given reflector. Simulation results showed that the surface normal vector variation takes little percentage in the distorted reflector pattern analysis and its contribution is in the order of 0.001 dB which can be neglected in the approximation for errors profiles up to about0.1
λ RMS amplitude.



---

*Source: 4824873-2016-07-25.xml*