International Journal of Antennas and Propagation
Volume 2016 (2016), Article ID 4824873, 6 pages
http://dx.doi.org/10.1155/2016/4824873

Research Article

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

Shuxin Zhang, Guigeng Yang, and Yiqun Zhang

Key Laboratory of Electronic Equipment Structure Design, Ministry of Education, Xidian University, Xi’an 710071, China

Received 22 March 2016; Revised 10 June 2016; Accepted 19 June 2016

Academic Editor: Toni Björninen

Copyright © 2016 Shuxin Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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.

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: Section 2 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]

where is the surface normal vector, is the incident magnetic field, the vector locates the integration point, is the imaginary unit, is the wavenumber, is the free-space wavelength at the working frequency, the unit vector 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 in the -direction shown in Figure 1, 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

where is the normal vector variation, is the undistorted spherical component in the feed coordinate system as shown in Figure 1, and defines the observation direction. Both and 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 is defined as the reflector unprimed coordinate system and is for the feed coordinate system. Shown in Figure 1, denotes the observation point and is the surface point vector in global coordinate system.

Figure 1: A distorted reflector.

The surface nodal displacements shown in Figure 1 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  RMS amplitude [16]:

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

where is the surface nodal -displacement of the th node and is the gradient of normal vector with respect to the th nodal displacement.

Because of the similarity of the equilateral triangles, two triangles and are employed in Figure 2 to show the gradient of surface normal vector with respect to nodal displacement, where the symbols , , , and denote the triangular nodes. The gradients of triangular normal vector with respect to each nodal displacement are listed in Table 1, where 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

Figure 2: Two equilateral triangles in aperture plane.

The same as in [17], the triangular element in aperture plane shown in Figure 3 is chosen in the analysis, where the symbols , , and represent the triangle nodes and is an arbitrary point inside the element. Supposing that the nodal displacements in -direction are labeled as , , and for nodes , , and , respectively, the displacement of an integration point in (3) within the element can be expressed as a function of nodal displacements multiplied by the element shape functions, and

where , , and are the element shape functions with respect to the nodes , , and , respectively.

Figure 3: A triangular element in FEM analysis [17].

Substituting (3), (4), and (5) into (2) and neglecting the high-order terms, the integral on the triangular element (labeled with as superscript) can be expanded as

where the superscript designates the element , is the undistorted contribution in the integral, and and are the linear and square terms in [17]; the first three terms in the right side are not repeated here for brevity, and

where , , and are the nodal displacements of nodes , , and , where , , and , , , and , respectively, and

where is the integral region of element , is the gradient of surface normal vector with respect to th nodal displacement, are the element shape function with respect to th node in element . From its expression, it can be seen that 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:

where is the nodal displacements vector in the element , the operator converts a matrix into a vector, and are the same coefficients as in [17], and

Assuming that there are elements in the FEM analysis and observation points in the far field requirements, the extended final matrix form can be expressed by assembling the whole elements integrals as

where is the -direction displacements of all nodes, is the product term between each two nodal displacements in one element including square terms of each nodal displacement, and

where , is the undistorted electrical far field vector, and 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 , , , and and the undistorted electrical far field 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  RMS amplitude, a simulation is performed in this study. The contributions of the matrices and in (12) are computed with the former matrices and . For the distorted reflector in [16, 17] of which the diameter , the focal length is adopted in this simulation. The feed parameters are also the same as the ones in [16, 17] with 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

where the variables and are the same as in [17] which describes the polar coordinate components in the projected aperture plane of the reflector.

Figures 4–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 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  RMS amplitude. It can be understood from this that, in the equivalent induced current density () 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.214

Figure 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 about  RMS amplitude.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by the China National 973 Program under Grant no. 2015CB857100 and the National Natural Science Foundation of China under Grant nos. 51490660, 51405364, 51305323, 51475348, and 51475349.

References

Y. Rahmat-Samii and R. L. Haupt, “Reflector antenna developments: a perspective on the past, present and future,” IEEE Antennas and Propagation Magazine, vol. 57, no. 2, pp. 85–95, 2015.

Y. Rahmat-Samii and A. C. Densmore, “Technology trends and challenges of antennas for satellite communication systems,” IEEE Transactions on Antennas and Propagation, vol. 63, no. 4, pp. 1191–1204, 2015.

J. Ruze, “Antenna tolerance theory—a review,” Proceedings of the IEEE, vol. 54, no. 4, pp. 633–640, 1966.

A. Greve and B. G. Hooghoudt, “Quality evaluation of radio reflector surfaces,” Astronomy & Astrophysics, vol. 93, no. 1, pp. 76–78, 1981.

N. Li, B. Duan, and F. Zheng, “Effect of the random error on the radiation characteristic of the reflector antenna based on two-dimensional fractal,” International Journal of Antennas and Propagation, vol. 2012, Article ID 543462, 6 pages, 2012.

P. Rocca, N. Anselmi, and A. Massa, “Interval arithmetic for pattern tolerance analysis of parabolic reflectors,” IEEE Transactions on Antennas and Propagation, vol. 62, no. 10, pp. 4952–4960, 2014.

P. Rocca, L. Manica, and A. Massa, “Interval-based analysis of pattern distortions in reflector antennas with bump-like surface deformations,” IET Microwaves, Antennas & Propagation, vol. 8, no. 15, pp. 1277–1285, 2014.

B. Gonzalez-Valdes, J. A. Martinez-Lorenzo, C. M. Rappaport, and A. G. Pino, “A new physical optics based approach to subreflector shaping for reflector antenna distortion compensation,” IEEE Transactions on Antennas and Propagation, vol. 61, no. 1, pp. 467–472, 2013.

P. Li, B. Y. Duan, W. Wang, and F. Zheng, “Electromechanical coupling analysis of ground reflector antennas under solar radiation,” IEEE Antennas and Propagation Magazine, vol. 54, no. 5, pp. 40–57, 2012.

B. Y. Duan and C. S. Wang, “Reflector antenna distortion analysis using MEFCM,” IEEE Transactions on Antennas and Propagation, vol. 57, no. 10, pp. 3409–3413, 2009.

K. Bahadori and Y. Rahmat-Samii, “Characterization of effects of periodic and aperiodic surface distortions on membrane reflector antennas,” IEEE Transactions on Antennas and Propagation, vol. 53, no. 9, pp. 2782–2791, 2005.

S. L. Padula, H. M. Adelman, M. C. Bailey, and R. T. Haftka, “Integrated structural electromagnetic shape control of large space antenna reflectors,” AIAA Journal, vol. 27, no. 6, pp. 814–819, 1989.

S. Zhang, J. Du, B. Duan, G. Yang, and Y. Ma, “Integrated structural-electromagnetic shape control of cable mesh reflector antennas,” AIAA Journal, vol. 53, no. 5, pp. 1395–1398, 2015.

J. S. Liu and L. Hollaway, “Integrated structure-electromagnetic optimization of large reflector antenna systems,” Structural Optimization, vol. 16, no. 1, pp. 29–36, 1998.

C. S. Wang, B. Y. Duan, and Y. Y. Qiu, “On distorted surface analysis and multidisciplinary structural optimization of large reflector antennas,” Structural and Multidisciplinary Optimization, vol. 33, no. 6, pp. 519–528, 2007.

W. T. Smith and R. J. Bastian, “An approximation of the radiation integral for distorted reflector antennas using surface-error decomposition,” IEEE Transactions on Antennas and Propagation, vol. 45, no. 1, pp. 5–10, 1997.

S. Zhang, B. Duan, G. Yang, Y. Zong, and Y. Zhang, “An approximation of pattern analysis for distorted reflector antennas using structural-electromagnetic coupling model,” IEEE Transactions on Antennas and Propagation, vol. 61, no. 9, pp. 4844–4847, 2013.

B. Hassani and E. Hinton, “A review of homogenization and topology opimization II—analytical and numerical solution of homogenization equations,” Computers and Structures, vol. 69, no. 6, pp. 719–738, 1998.