Abstract

A novel noniterative method for pattern synthesis of conformal antenna array is put forward. First, a new pattern function formula for general conformal antenna array is derived. Then, according to the new pattern function formula, a new block matrix equation (BME) and a simple matrix equation (SME) are obtained. SME has the same form as the equation of linear array pattern synthesis. The new method can be applied to pattern synthesis of any conformal array. Moreover, due to SME having the same form as the matrix equation of linear array pattern synthesis, the new simple matrix equation can be dealt with by the approaches for linear array pattern synthesis, which significantly expands the existing approaches for conformal antenna array pattern synthesis. Two different conformal array configurations are taken as the examples to demonstrate the advantages of the new method. Results of the simulations show that the new method can flexibly and effectively be applied to synthesize patterns for various conformal array architectures.

1. Introduction

Design and analysis of conformal array have always been a difficult work since conformal array was presented several decades ago for its complicated architecture. Some special and frequently used conformal array configurations have been mainly explored by researchers. The cylindrical conformal antenna arrays were studied in references [1, 2]. In references [3, 4], the conical conformal antenna arrays were investigated. A method of radiation pattern computation of pyramidal conformal antenna array was developed in reference [5]. In references [69], the spherical conformal antenna arrays were researched.

The antenna elements of conformal array are mounted on the surface of the host and the shape of the surface was not allowed to be changed, which can make conformal array meet the requirements of a lot of special individual applications and also make every antenna element axis have different orientation. Hence, the conformal array’s radiated beam pattern cannot be calculated from array beam pattern multiplication theorem—array element directivity function multiplying array factor.

Pattern synthesis of conformal antenna array is much more complex than linear antenna array, and the methods of conformal array pattern synthesis are much less than those for linear arrays. In references [10, 11], Euler’s rotation approach was utilized to generate conformal beam pattern. Based on Euler’s rotation approach, rectangular microstrip patch antenna on spherical structure was analyzed [6], frequency-invariant pattern of conformal array antenna with low cross-polarisation was synthesized [12], an adaptive wideband beam forming algorithm was presented to form beam for conformal array [13], and optimum pattern synthesis of nonuniform spherical array was discussed in reference [14]. Three smart approaches were presented in the article [1517] to synthesize pattern for conformal array. In reference [18], the constructive analytical phasing (CAP) method based on a new coordinate rotation formula was presented to calculate the pattern of conformal array. In references [1921], the geometric algebra, finite-difference time-domain (FDTD), and the moment method were developed for conformal array pattern synthesis. In reference [22], a multidimension to one-dimension transformation approach was developed for the similar antenna array pattern synthesis.

Most of these methods mentioned above for conformal array pattern synthesis were based on Euler’s rotation approach, in which the antenna element pattern was derived through its coordinate rotation operation and the process was complicated. All the solutions of these methods for conformal array were acquired by iterative computation formulas. These computation formulas for conformal array were different from those for linear array, and their computation processes were much more complicated. Moreover, lots of methods for linear array cannot be applied in the process for conformal array pattern synthesis. This situation significantly impacts the effectiveness of these methods for conformal array. Moreover, there is convergence problem in smart approaches. The reported general success rate of smart methods is not more than 60%, and the success case of smart methods for the element number of array more than 50 has not been reported due to the huge complexity and a significant drop in the success rate. The success rate drops very fast while the element number of array increases.

To overcome these problems, this article presents a new simple analytical pattern formula for conformal array pattern synthesis. A novel pattern synthesis block matrix formula for conformal array is put forward first. Then, a new method is presented to transform the matrix formula into a form similar to the matrix equation of the linear array, so all the abundant methods for linear array can be utilized to synthesize the beam pattern for conformal array, and in particular, many noniterative computation formulas, such as the least square method and the Fourier transform (FT) method, can be applied for conformal array pattern synthesis, which can significantly improve the diversity of the method of pattern synthesis for conformal array and reduce the amount of the computation. In addition, the new method is a general approach that can be applied to any conformal array architecture.

The rest of the contents of this paper are organized as follows: in Section 2, new mathematical model and analytical formula for pattern synthesis of conformal array are presented; in Section 3, a new block matrix equation for pattern synthesis of conformal array is developed; in Section 4, a new simple matrix equation for pattern synthesis of conformal arrays is presented. It can be easily solved by using methods as similar as those for linear array pattern synthesis; examples are taken to demonstrate the effect of the new method in Section 5; and Section 6 draws conclusions.

2. Conformal Array Pattern Function Formula

The configuration of an arbitrary conformal antenna array is considered, as shown in Figure 1. Assuming the element number of the array to be M, the array’s antenna element location is , and the current excitation of the element is . Then, the beam pattern function of the conformal array can be expressed as follows:wherein is the antenna element’s pattern function, φ is the azimuth angle of a far field point P(r, θ, φ), θ is the elevation angle of point P, r is the distance of point P, and is the beam pattern function of the conformal array.

Assuming that all antenna elements of the conformal antenna array are homogeneous but with different axis orientations and locations and that the antenna element’s axis orientation is (, ), we set up a new local coordinate system for the antenna element with z’ coordinate along its axis orientation (, ). In the antenna element’s new local coordinate system, the space direction coordinate can be written as (, ). The antenna element’s pattern function is denoted as f(, ). According to Euler’s right-handed rotation method in the x-y-z rotation order with the rotation angles , , and , in turn [10], and can be derived as follows:

In far field,wherein λ is the wavelength of the radiated signal. Therefore, equation (1) can be re-expressed as follows:

Letwhere in equation (6), the superscript T denotes transpose operation. Hence, equation (1) can be expressed as follows:

Equation (8) is the pattern function of the conformal array. It can be transformed into a simple matrix equation to solve.

3. Block Matrix Equation of Conformal Array Pattern Synthesis

We set θ in the range of [0°, 180°] and φ in [0°, 180°]. Their discrete values are θ1, θ2, ..., θk1, ..., θK1 and φ1, φ2, ..., φk2, ..., φK2, respectively.

Letwherein F is the discrete matrix of the pattern function in far field.

Letwherein Sm is the discrete matrix of the function .

Let

So, equation (8) can be re-expressed as follows:

This is the matrix equation of conformal array pattern synthesis. But this equation is a block matrix equation (BME).

4. Simple Matrix Equation of Conformal Array Pattern Synthesis

Directly solving equation (12) is difficult; we transform this equation into an easily solved new equation and then solve it. Lettherefore

Similarly, lettherefore

According to equation (12), we can get

Let

Similarly, let

Therefore, equation (17) can be re-expressed as follows:

It is apparent that equation (20) is a general simple matrix equation, not a block matrix equation, and b is a column vector. Equation (20) has the similar form as the matrix equation of the linear array pattern synthesis. So, equation (20) can be solved by using the existing methods for linear array pattern synthesis. On the one hand, equation (20) can be used to calculate the beam pattern of conformal array from A and I. On the other hand, when the expected pattern b is prescribed, the antenna element excitation vector I can be solved from equation (20). Equation (20) is named as the simple matrix equation (SME) method for conformal array pattern synthesis.

5. Examples of the New Method

Two examples are taken to demonstrate the performance of the new method and its effect in this part.

Example 1. A cylinder conformal array configuration is considered in Figure 2. There are N = 30 elements equidistantly distributed as the distance λ/2 along a circle centered at the origin point O on the x-o-y plane with radius a= 1 m. This circle is the first layer of the array. The first layer is parallelly moved up along  + z coordinate axis direction to get other 9 layers of the cylinder array by (k − 1)λ/2 distance with k referring to the kth layer. It is worth mentioning that every other layer is a copy of the first layer in every aspect aside from having different z coordinates. λ is the central wavelength of the signal transmitted by the array. φ is the azimuth angle of a far field point P(r, θ, φ), θ is the elevation angle of point P, and r is the distance of point P from the origin point O. In the first layer, it is easy to get the spherical coordinate location of the mth antenna element, that is, (a, π/2, 2π(m − 1)/N). Assuming that all antenna elements are homogeneous half wavelength patch dipole antennae but with different axis orientations and locations, let the mth antenna element of the first layer axis orientation be (βm, γm). We set up a new local coordinate system x’-y’-z’ for the mth antenna element of the first layer with z’ coordinate along its axis orientation (βm, γm). In the first layer mth antenna element’s new local coordinate system, the space direction coordinate can be written as (, ). According to Euler’s right-handed rotation method in the x-y-z rotation order with the rotation angles αm = 0, βm = π/6, and γm = 2π(m − 1)/N in turn [10], and can be derived from equations (2) and (3). The first layer mth antenna element pattern function f(, ) can be expressed as follows:Because every other layer is a copy of the first layer in every aspect aside from having different z coordinates, for mth antenna element of every layer in every layer’s local coordinate system, the pattern function of the element is the same as equation (21).
The new method is applied to synthesize the beam pattern for the cylinder array. We set θ in the range of [0°, 180°] and φ in [0°, 180°]. Their discrete values θ1, θ2, ..., θk1, ..., θK1 and φ1, φ2, ..., φk2, ..., φK2 all are equidistantly distributed with the adjacent distance 1°. When the target pattern is given, the antenna excitation vector can be obtained from equation (20) [22] as follows:where the superscript H denotes conjugation transpose operation, andwhere diag denotes diagonalization operation, and abs denotes getting the absolute value operation.
Here, the target pattern b is generated by using MATLAB; the matrix A is obtained according to the new method mentioned above. After I has been obtained from equation (22), the pattern formed by I can be obtained from equation (20).
MATLAB software is applied to simulate this example. The three-dimension (3D) figures of the target pattern and the synthesized pattern using the new method are shown in Figures 3(a) and 3(b), respectively.
From Figure 3, it can be seen that the new method well accomplishes the beam pattern synthesis mission of the cylinder conformal array. The location and the attenuation of the synthesized beam well accord with the target beam.
The two-dimension (2D) figures of simulation results are shown in Figures 4(a) and 4(b), where SME denotes the result of the simple matrix equation method, TARG refers to the target beam pattern, and CAP refers to the result of constructive analytical phasing (CAP) method in reference [18] as a comparison object. In order to further verify the new method, ANSYS HFSS software is applied to simulate the example. HFSS simulation results are also shown in Figures 4(a) and 4(b) with the data marked by HFSS. Figure 4(a) is the results of φ = 60°. Figure 4(b) is the results of θ = 60°.
ANSYS HFSS software is used to simulate this example with the signal wavelength λ = 4sin(π/N) = 7.3 mm. The rogers5880 material substrate with thickness d1 = 0.508 mm is used to bear the half wavelength dipole microwave patch antenna with the copper thickness d2 = 35 μm. The dielectric permittivity εr of rogers5880 material is 2.2, and its loss tangent tanδ is 0.001.
From Figure 4, it can be seen that the result of the new method marked by SME has the higher attenuation than CAP and HFSS simulation results. The difference between SME and HFSS should be due to the near field mutual coupling among the antenna elements of the array. The near field mutual coupling will not be discussed here in order to focus on the new method.

Example 2. A semispherical conformal array configuration is considered and shown in Figure 5. The semisphere is centered at the origin point O with radius a= 1 m. Let φ be the azimuth angle and θ be the elevation angle of a space direction. As shown in Figure 5, the benchmark line is the arc on the semisphere when θ from 0 to 90° and φ = 0°. There are N = 11 elements equidistantly distributed as the distance λ/2 along the benchmark line shown as 1th to Nth in Figure 5. The equator circle is the circle on the semisphere and x-o-y plane. It is easy to see that the Nth antenna element is on the equator circle. The mth latitude circle on the semisphere parallels to the equator circle and goes through the mth antenna element on the benchmark line with m from 1 to N − 1. On the equator circle and all the latitude circles, there are antenna elements equidistantly distributed as the distance λ/2 along each circle. λ is the central wavelength of the signal transmitted by the array. Assuming all antenna elements are homogeneous half wavelength patch dipole antennae but with different axis orientations and locations, let the kth antenna element’s spherical coordinates be (a, , and ). Assuming its axis orientation being (βk, γk) and βk = , γk = . According to Euler’s right-handed rotation method in the x-y-z rotation order with the rotation angles αk = 0, βk, and γk in turn, the array kth antenna element pattern function can be obtained from equation (22).
Similar to Example 1, the target pattern b is generated by using MATLAB; the matrix A is obtained according to the new method mentioned above. After I has been obtained from equation (22), the pattern formed by I can be obtained from equation (20).
The 3D figures of the target pattern and the synthesized pattern using the new method are shown in Figures 6(a) and 6(b), respectively. From Figure 6, it can be learned that the new method well accomplishes the beam pattern synthesis mission of the semisphere conformal array. The location and the attenuation of the synthesized beam well accord with the target beam.
The 2D figures of simulation results are shown in Figures 7(a) and 7(b), where SME denotes the result of the simple matrix equation method, TARG refers to the target beam pattern, and CAP refers to the result of constructive analytical phasing (CAP) method in reference [18] as a comparison object. ANSYS HFSS software simulation results are also shown in Figures 7(a) and 7(b) with the data marked by HFSS. Figure 7(a) is the results of φ = 60°. Figure 7(b) is the results of θ = 60°.
ANSYS HFSS software is still used to simulate this example with the signal wavelength λ = 4sin(π/4/N) = 5.4 mm. The rogers5880 material substrate with thickness d1 = 0.508 mm is used to bear the half wavelength dipole microwave patch antenna with the copper thickness d2 = 35 μm.
From Figure 7, it can be seen that the results of the new method marked by SME have higher attenuation than CAP and HFSS simulation results. The difference between SME and HFSS should be due to the near-field mutual coupling among the antenna elements of the array. The near field mutual coupling will not be discussed here in order to focus on the new method.
The SME method is a noniterative method. The results of SME can be obtained through equation (22) straightly. Therefore, the new method has much less computation complexity than published conformal array pattern synthesis approaches that all are iterative methods. Both of the abovementioned examples are simulated on MATLAB software platform on a laptop with Intel i7-7700HQ and 16G memory. It takes no more than one minute to complete the simulation process in each example. The computing process is swift.
The results of the examples show that the new method has excellent effect for conformal array pattern synthesis with little amount of computation. In addition, the adjacent spacing between the discrete values of θ and φ can be arbitrarily set wider to further reduce the amount of computation, and satisfactory results can also be achieved in the light of our experiments.

6. Conclusions

An analytical function formula is developed for conformal antenna array pattern synthesis by using Euler’s rotation approach. By a new way, the function is discretized to a block matrix equation (BME). Based on the BME, a simple matrix equation (SME) method is developed for conformal array pattern synthesis. Due to the similar form with the matrix equation of the linear array pattern synthesis, SME can be solved by many approaches for linear array pattern synthesis. The result of SME can be straightly obtained through a noniterative matrix equation. Therefore, the new method has much less computation complexity than published conformal array pattern synthesis approaches that all are iterative methods. Two examples are taken to demonstrate the effectiveness and flexibility of the new method. The results of the simulations show that the new method is highly efficient with less amount of computation.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.