*Article* **Antenna Current Calculation Based on Equivalent Transmission Line Model**

**Shusheng Wei <sup>1</sup> and Wusong Wen 2,\***


**Abstract:** This paper provides a new way for spatial current/field profiles for frequency-selective surface analytical approximation. It confirms that the per unit length radiation resistance of an equivalent transmission line model for line antenna has little influence on the normalized current distribution. The two-wire equivalent transmission line model (typically used for transmitting line antenna) is applied to the receiving line antenna. In this case, the corresponding incident field is decomposed into odd and even mode for asymmetric distribution. A one-wire equivalent transmission line model is then introduced for any antenna composed of relative narrow strips. The incident field does not need to be decomposed. According to the simulation, the transmission line loss has little influence on the current distribution.

**Keywords:** antenna current; transmission line model; frequency-selective surface analytical approximation

## **1. Introduction**

Numerical methods, such as method of moment (MoM), are usually applied to calculate the transmitting current of dipole antenna. However, the current distribution calculated is in numerical form, which cannot be directly used for frequency-selective surface (FSS), where the analytical current distribution is needed [1–7]. Two promising approximate analytical methods for FSS are the periodic method of moment (PMM) and multimodal network approach. However, it is necessary for both of them to obtain the current/field distribution in the scatterer. PMM solves this problem using MoM [8–10], while the multimodal network approach only considers simple rectangular scatters [11–13]. For a rectangular patch (or aperture) under oblique incidence in the principle scan plane of the structure, approximate closed-form expressions for aperture field and current profiles of patch surface were obtained. However, to the authors' knowledge, the approximate closed-form expressions for other structures have not been reported. Therefore, the application of the above approximate analytical methods depends on the development of current/field distribution in the scatterer of other types. In [14], the spatial current/field profile is calculated using a full-wave simulation, at a single- and low-frequency value. It is then assumed to be independent of the frequency in the considered range of interest. However, this method is restrictive, treating structures that are relatively simple and incidence angles that are relatively small. Besides, the numerical preparation is time consuming and means the analytical method is complex.

The calculation of spatial current/field profile for an FSS element is similar to the analysis of current distribution for line antenna [15,16]. The equivalent transmission line model was often adopted for the study of input impedance and mutual coupling of transmitting line antenna [17–19]. It was usually a two-wire transmission line model with lossy lumped feed. For single antennas, radiation lossy can be considered easily, as can the

**Citation:** Wei, S.; Wen, W. Antenna Current Calculation Based on Equivalent Transmission Line Model. *Micromachines* **2022**, *13*, 714. https://doi.org/10.3390/mi13050714

Academic Editor: Wensheng Zhao

Received: 8 April 2022 Accepted: 27 April 2022 Published: 30 April 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

mutual coupling lossy. A two-wire straight antenna refers to two collinear or equidistant antennas, which can be one or more segments. The segments may or may not be parallel. The two-wire straight antenna can only be equivalent to the two-wire transmission line model, and the requirement is that the antenna structures are symmetrical. Besides, the tangential component of the incident electric field on the two antennas should be even or odd symmetrical. For the asymmetric distribution of the incident field, the incident field can be decomposed into odd and even symmetrical. The total current in the antenna is then the superposition of currents in the above two modes. structures, the main problem is that the equations and parameters are restrictive in approximating the lossless case. Schelkunoff also studied the current distribution on the receiving antenna (actually a reflector antenna, which is a straight continuous fine wire antenna, which is exposed to the incident field and the corresponding tangential component is uniform) using a two-wire transmission line model. The model is restrictive for vertical incidence (the tangential field of incident field is uniform on the antenna). For oblique

Schelkunoff mainly used a one-wire transmission line model in [17] to study the input impedance, mutual coupling, antenna impedance, and the approximate solution of transmitting/receiving antenna current. Compared to the two-wire transmission line model, the one-wire transmission model is more flexible. The former makes it easier to study the dipole oscillator antenna. For the one-wire transmission line model of complex structures, the main problem is that the equations and parameters are restrictive in approximating the lossless case. Schelkunoff also studied the current distribution on the receiving antenna (actually a reflector antenna, which is a straight continuous fine wire antenna, which is exposed to the incident field and the corresponding tangential component is uniform) using a two-wire transmission line model. The model is restrictive for vertical incidence (the tangential field of incident field is uniform on the antenna). For oblique incidence, the tangential electric field of the incident field on the antenna varies with the spatial position. It results in the induced electromotive force distributed asymmetrically about the center of the antenna. The currents on the lines are not equivalent and reverse, so the two-wire model cannot be used. incidence, the tangential electric field of the incident field on the antenna varies with the spatial position. It results in the induced electromotive force distributed asymmetrically about the center of the antenna. The currents on the lines are not equivalent and reverse, so the two-wire model cannot be used. In this paper, the current distributions of transmitting and receiving antenna, based on the equivalent transmission line method, are analyzed. Both the two-wire and one-wire model are studied, together with the radiation lossy. An iterative scheme is then introduced to further improve the accuracy. The current distribution calculated is in analytic form, which can be directly used to approximate the analytical method for FSS. **2. Transmitting Antenna**

In this paper, the current distributions of transmitting and receiving antenna, based on the equivalent transmission line method, are analyzed. Both the two-wire and one-wire model are studied, together with the radiation lossy. An iterative scheme is then introduced to further improve the accuracy. The current distribution calculated is in analytic form, which can be directly used to approximate the analytical method for FSS. *Equivalent Circuit of the Transformers* Symmetric diploes transmitting antenna (Figure 1) can be equivalent to the two-wire transmission line model (Figure 2a). For the lossless equivalent transmission line corre-

#### **2. Transmitting Antenna** sponding to equidistance symmetric diploes antenna (Figure 1a), the average characteris-

*Micromachines* **2022**, *13*, x FOR PEER REVIEW 2 of 10

then the superposition of currents in the above two modes.

lumped feed. For single antennas, radiation lossy can be considered easily, as can the mutual coupling lossy. A two-wire straight antenna refers to two collinear or equidistant antennas, which can be one or more segments. The segments may or may not be parallel. The two-wire straight antenna can only be equivalent to the two-wire transmission line model, and the requirement is that the antenna structures are symmetrical. Besides, the tangential component of the incident electric field on the two antennas should be even or odd symmetrical. For the asymmetric distribution of the incident field, the incident field can be decomposed into odd and even symmetrical. The total current in the antenna is

Schelkunoff mainly used a one-wire transmission line model in [17] to study the in-

put impedance, mutual coupling, antenna impedance, and the approximate solution of transmitting/receiving antenna current. Compared to the two-wire transmission line

#### *Equivalent Circuit of the Transformers* tic impedance is evaluated using [14,15]

Symmetric diploes transmitting antenna (Figure 1) can be equivalent to the two-wire transmission line model (Figure 2a). For the lossless equivalent transmission line corresponding to equidistance symmetric diploes antenna (Figure 1a), the average characteristic impedance is evaluated using [14,15] 0 <sup>120</sup> ln *r <sup>d</sup> <sup>Z</sup>* ε ρ <sup>=</sup> .

$$Z\_0 = \frac{120}{\sqrt{\varepsilon\_r}} \ln \frac{d}{\rho}.\tag{1}$$

(1)

**Figure 1.** Symmetric diploes transmitting **Figure 1.** Symmetric diploes transmitting line antenna ( line antenna (**a**) equidistance; ( **a**) equidistance; ( **b**) collinear. **b**) collinear.

**Figure 2.** Two-wire equivalent transmission line model of Figure 1. (**a**) Equivalent transmission line; (**b**) cascaded two-port network. **Figure 2.** Two-wire equivalent transmission line model of Figure 1. (**a**) Equivalent transmission line; (**b**) cascaded two-port network.

For the lossless equivalent transmission line corresponding to collinear symmetric diploes antenna (Figure 1b), the average characteristic impedance is calculated as [17,18] For the lossless equivalent transmission line corresponding to collinear symmetric diploes antenna (Figure 1b), the average characteristic impedance is calculated as [17,18]

$$Z\_0 = \frac{1}{l - l\_1} \int\_{l\_1}^{l} \frac{120}{\sqrt{\varepsilon\_r}} \ln \frac{2(z - l\_1)}{\rho} dz. \tag{2}$$
 
$$\text{dactance and capacity can then be obtained}$$

(2)

(3)

(4)

(5)

(6)

1 0 1 *l r l l* ε ρ The per-unit-length inductance and capacity can then be obtained

$$L = \sqrt{\mu\_0 \varepsilon\_0 \varepsilon\_r} Z\_0, \; \mathbb{C} = \frac{\sqrt{\mu\_0 \varepsilon\_0 \varepsilon\_r}}{Z\_0}. \tag{3}$$

*L Z* = µεε0 0 *<sup>r</sup> C Z* µεε = The propagation constant is *r*<sup>1</sup> = *jβ* = *jω* √ *LC* = *jω* √*µ*0*ε*0*ε<sup>r</sup>* .

 00 0 *<sup>r</sup>* , 0 . The propagation constant is <sup>1</sup> 0 0 *<sup>r</sup> r j j LC j* = = = β ω ω µεε . According to the theory of the transmission line [20,21], the equivalent transmission According to the theory of the transmission line [20,21], the equivalent transmission line model in Figure 2a can be equivalent to the cascaded two-port network in Figure 2b. *V*<sup>i</sup> and *I*<sup>i</sup> denote the corresponding port voltage and current, respectively. A voltage source *V*<sup>1</sup> = *V*e*<sup>s</sup>* is connected at the left-hand-side port and the right-hand-side is an open circuit (*I*<sup>2</sup> = 0). The transmission matrix *A*<sup>1</sup> can be expressed as

$$A\_1 = \begin{bmatrix} \cos[\beta(l-l\_1)] & -jZ\_0 \sin[\beta(l-l\_1)] \\ \frac{\sin[\beta(l-l\_1)]}{jZ\_0} & \cos[\beta(l-l\_1)] \end{bmatrix}. \tag{4}$$

<sup>2</sup> *I* = 0 ). The transmission matrix *A*<sup>1</sup> can be expressed as According to the microwave network theory, we obtain

*V V* <sup>1</sup> = *<sup>s</sup>*

when *l* = 0.05

λ

viation for MoM when 1.0 *<sup>r</sup>*

$$
\begin{bmatrix} V\_2 \\ I\_2 \end{bmatrix} = A\_1 \begin{bmatrix} V\_1 \\ I\_1 \end{bmatrix}. \tag{5}
$$

$$
\text{ry condition } V\_1 = \tilde{V}\_s \text{ and } I\_2 = 0 \text{ into Equation (5), we can get voltages and currents}
$$

1 0 cos *l l jZ* β Substituting the boundary condition *V*<sup>1</sup> = *V*e*<sup>s</sup>* and *I*<sup>2</sup> = 0 into Equation (5), we can obtain the coefficients of the voltages and currents

$$V\_2 = \frac{\tilde{V}\_s}{\cos[\beta(l-l\_1)]}, \ I\_1 = j \tan[\beta(l-l\_1)] \frac{\tilde{V}\_s}{Z\_0}.\tag{6}$$

2 1

1 2 1 *A I I* <sup>=</sup> . Based on the above port voltages and currents, the currents distribution of the equivalent transmission lines (Figure 1) can be written as

$$I(z) = I\_1 \cos[\beta(z - l\_1)] + \frac{V\_1}{jZ\_0} \sin[\beta(z - l\_1)].\tag{7}$$

1 1 1

The simulated normalized current distribution of the proposed method and MoM are shown in Figure 3. It can be observed that the results comply well with each other

0

 β

≠ . It should be pointed out that this method can be extended

. (7)

λ

. In fact, there is de-

*jZ* = −+ −

( ) <sup>2</sup> <sup>1</sup> cos *Vs <sup>V</sup>* β *l l* <sup>=</sup> <sup>−</sup> , 1 1 ( ) 0 tan *Vs I j ll <sup>Z</sup>* = − β . Based on the above port voltages and currents, the currents distribution of the equiv-The simulated normalized current distribution of the proposed method and MoM are shown in Figure 3. It can be observed that the results comply well with each other when *l* = 0.05*λ*. However, they differ from each other when *l* = 0.25*λ*. In fact, there is deviation for MoM when *ε<sup>r</sup>* 6= 1.0. It should be pointed out that this method can be extended to multi-section line antenna. In this case, the equivalent transmission line model

β

. However, they differ from each other when *l* = 0.25

ε

alent transmission lines (Figure 1) can be written as

0.6

Normalized Current

0.8

1.0

ysis of the complicated line antenna.

0.25.

to multi-section line antenna. In this case, the equivalent transmission line model and two-

1.0

to multi-section line antenna. In this case, the equivalent transmission line model and twoport network are cascaded by the corresponding multi sections. It would affiliate the anal-

*Micromachines* **2022**, *13*, x FOR PEER REVIEW 4 of 10

Proposed MoM

*Micromachines* **2022**, *13*, x FOR PEER REVIEW 4 of 10

and two-port network are cascaded by the corresponding multi sections. It would affiliate the analysis of the complicated line antenna. port network are cascaded by the corresponding multi sections. It would affiliate the analysis of the complicated line antenna. 0.8 

**Figure 3.** Normalized current distribution of the proposed method and MoM. (**a**) = 0.05; (**b**) = **Figure 3.** Normalized current distribution of the proposed method and MoM. (**a**) *l* = 0.05*λ*; (**b**) *l* = 0.25*λ*. As discussed in the introduction, the necessary condition for a single straight antenna to be equivalent to a two-wire transmission line model is that the incident field of the

#### **3. Two-Wire Transmission Line Model of Receiving Antenna** antenna is symmetrical. Mode decomposition can be applied to asymmetric situations. In

**3. Two-Wire Transmission Line Model of Receiving Antenna** As discussed in the introduction, the necessary condition for a single straight antenna to be equivalent to a two-wire transmission line model is that the incident field of the antenna is symmetrical. Mode decomposition can be applied to asymmetric situations. In As discussed in the introduction, the necessary condition for a single straight antenna to be equivalent to a two-wire transmission line model is that the incident field of the antenna is symmetrical. Mode decomposition can be applied to asymmetric situations. In this section, the equivalent transmission model is applied to the analysis of receiving antenna, shown in Figure 4. this section, the equivalent transmission model is applied to the analysis of receiving antenna, shown in Figure 4.

**Figure 4.** Two-wire straight antenna exposed to incident electric field.

*<sup>z</sup>* <sup>=</sup> <sup>0</sup> <sup>−</sup>*<sup>l</sup> <sup>l</sup>* <sup>1</sup> *l* <sup>1</sup> −*l* **Figure 4.** Two-wire straight antenna exposed to incident electric field. The two-wire lossless transmission line model of Figure 4 is shown in Figure 5, where The two-wire lossless transmission line model of Figure 4 is shown in Figure 5, where *V*(*z*) is the voltage between the two lines. The voltage and current on the line are governed by the transmission line equations

*z*

*z*

<sup>∂</sup> (8)

<sup>∂</sup> (8)

$$\frac{\partial V(z)}{\partial z} = -jw\mathbf{L} \cdot I(z) + 2E\_z^i(z) \tag{8}$$

$$\frac{\partial I(z)}{\partial z} = -jw\mathbf{C} \cdot V(z) \tag{9}$$

=− ⋅ +

*V z*

*V z*

∂

∂

( ) ( ) <sup>2</sup> ( ) *<sup>i</sup>*

*jwL I z E z <sup>z</sup>*

=− ⋅ +

*jwL I z E z <sup>z</sup>*

( ) ( ) *I z*

( ) ( ) *I z jwC V z <sup>z</sup>*

erned by the transmission line equations

erned by the transmission line equations

− +

*j Cdz* ω

**Figure 5.** Two-wire transmission line model of Figure 4. **Figure 5.** Two-wire transmission line model of Figure 4. *Il Il*(−= += ) ( ) <sup>0</sup> . (10)

*Micromachines* **2022**, *13*, x FOR PEER REVIEW 5 of 10

dent field on the antenna is

 ω

ϕ

*V z*( )

where <sup>0</sup> 0 0 β ω µε<sup>=</sup> *<sup>j</sup>* ,

( ) *<sup>v</sup>* ϕ

*z* and

The corresponding transmission line parameters can be calculated following Section 2. The boundary condition is The corresponding transmission line parameters can be calculated following Section 2. The boundary condition is Substituting (10) into (8) and (9) yields the general expression of potential and current

Substituting (10) into (8) and (9) yields the general expression of potential and current

( ) ( )

*I z Ae Be z*

−

0

*Z*

*i i*

1 1

*rz rz*

obtained according to the boundary condition of the transmission line. Suppose the inci-

$$I(-l) = I(+l) = 0.\tag{10}$$

*Il Il* (−= += ) ( ) <sup>0</sup> . (10) Substituting (10) into (8) and (9) yields the general expression of potential and current ( ) ( ) 1 1 1 *rz rz i I z Ae Be z* ϕ − = −+ (11)

*I z*( ) ( ) *I z dI z* ( ) + ( ) *<sup>i</sup> E z dz <sup>z</sup> j Ldz*

ω

*i*

ϕ

$$\begin{array}{l} V(z) = Ae^{-r\_1 z} + Be^{r\_1 z} + \varphi\_v(z) \\ I(z) = \frac{1}{Z\_0} [Ae^{-r\_1 z} - Be^{r\_1 z} + \varphi\_i(z)] \end{array} \tag{11}$$

0 *Z* = −+ ( ) *<sup>v</sup>* ϕ *z* and ϕ*<sup>i</sup>* (*z*) are the terms corresponding to ( ) *<sup>i</sup> E z <sup>z</sup>* . The coefficients *A* and *B* can be obtained according to the boundary condition of the transmission line. Suppose the inci*ϕv*(*z*) and *ϕi*(*z*) are the terms corresponding to *E i z* (*z*). The coefficients *A* and *B* can be obtained according to the boundary condition of the transmission line. Suppose the incident field on the antenna is dent field on the antenna is ( ) <sup>0</sup> sin <sup>0</sup> cos *<sup>i</sup> j z Ez E e* β θ θ<sup>−</sup> = ⋅ . (12)

$$E^{\dot{i}}(z) = E\_0 \cos \theta \cdot e^{-j\beta\_0 \sin \theta z}. \tag{12}$$

( ) <sup>0</sup> sin <sup>0</sup> cos *<sup>i</sup> j z Ez E e* β θ θ <sup>−</sup> = ⋅ . (12) where *β*<sup>0</sup> = *jω* √*µ*0*ε*0, *<sup>θ</sup>* is the incident angle. Since the incident field is asymmetrically distributed, mode decomposition is required. where <sup>0</sup> 0 0 βµε <sup>=</sup> *<sup>j</sup>* , θ is the incident angle. Since the incident field is asymmetrically distributed, mode decomposition is required.

$$E\_{\varepsilon}^{i}(z) = \frac{E^{i}(z) + E^{i}(-z)}{2} = E\_{0}\cos\theta \cdot \cos(\beta\_{0}\sin\theta z) \tag{13}$$

$$E\_o^i(z) = \frac{E^i(z) - E^i(-z)}{2} = -jE\_0 \cos \theta \cdot \sin(\beta\_0 \sin \theta \cdot z) \tag{14}$$

( ) ( ) ( ) 0 0 cos sin sin ( ) <sup>2</sup> *i i i o Ez E z E z jE* θ βθ *<sup>z</sup>* − − = =− ⋅ ⋅ (14) The total current on the antenna is then obtained by the superposition of the above two modes. The corresponding simulated results are shown in Figure 6. The total current on the antenna is then obtained by the superposition of the above two modes. The corresponding simulated results are shown in Figure 6.

**Figure 6.** Normalized current distribution of two-wire transmission line model. *l*<sup>1</sup> = 0, *<sup>ρ</sup>* <sup>=</sup> <sup>1</sup> <sup>×</sup> <sup>10</sup>−<sup>7</sup> m, *<sup>f</sup>* <sup>=</sup> 10 GHz, *<sup>ε</sup><sup>r</sup>* <sup>=</sup> 1.6, *<sup>θ</sup>* <sup>=</sup> <sup>0</sup> ◦ . (**a**) *l* = 0.25*λ*; (**b**) *l* = 0.5*λ*.

#### **4. One-Wire Transmission Line Model of Receiving Antenna**

In this section, the one-wire transmission line model for a single straight antenna is analyzed. In this case, the mode decomposition of the incident field is eliminated. A single straight antenna is a continuous straight-line antenna, which could be one section or line equations.

*Micromachines* **2022**, *13*, x FOR PEER REVIEW 6 of 10

*Micromachines* **2022**, *13*, x FOR PEER REVIEW 6 of 10

**4. One-Wire Transmission Line Model of Receiving Antenna**

**4. One-Wire Transmission Line Model of Receiving Antenna**

1 × 10−7 m, = 10 GHz, = 1.6, = 0°

1 × 10−7 m, = 10 GHz, = 1.6, = 0°

multiple non-parallel sections. It can be equivalent to the one-wire equivalent transmission line model. Each section corresponds to a potential and a current equation. A single straight antenna with one section is illustrated in Figure 7, and *E i z* (*z*) is the tangential component in the incident electric field. The corresponding one-wire equivalent transmission line model is shown in Figure 8. For the transmission line excited by the distribution voltage source, the electric potential and current on the line are governed by the transmission line equations. ( ) ( ) ( ) *<sup>i</sup> z V z jwL I z E z <sup>z</sup>* ∂ =− ⋅ + <sup>∂</sup> (15) ( ) ( ) *I z jwC V z <sup>z</sup>* <sup>∂</sup> =− ⋅ <sup>∂</sup> (16) ( ) ( ) ( ) *<sup>i</sup> z V z jwL I z E z <sup>z</sup>* ∂ =− ⋅ + <sup>∂</sup> (15) ( ) ( ) *I z jwC V z <sup>z</sup>* <sup>∂</sup> =− ⋅ <sup>∂</sup> (16)

$$\frac{\partial V(z)}{\partial z} = -jwL \cdot I(z) + E\_z^i(z) \tag{15}$$

$$\frac{\partial I(z)}{\partial z} = -jw\mathbb{C} \cdot V(z) \tag{16}$$

*V*(*z*) and *I*(*z*) are potential and current, respectively. The boundary condition is *Il Il* (−= += ) ( ) 0 (17)

$$I(-l) = I(+l) = 0\tag{17}$$

**Figure 6.** Normalized current distribution of two-wire transmission line model. 1 = 0, =

**Figure 6.** Normalized current distribution of two-wire transmission line model. 1 = 0, =

. (**a**) = 0.25; (**b**) = 0.5.

. (**a**) = 0.25; (**b**) = 0.5.

In this section, the one-wire transmission line model for a single straight antenna is

In this section, the one-wire transmission line model for a single straight antenna is

analyzed. In this case, the mode decomposition of the incident field is eliminated. A single straight antenna is a continuous straight-line antenna, which could be one section or multiple non-parallel sections. It can be equivalent to the one-wire equivalent transmission line model. Each section corresponds to a potential and a current equation. A single straight antenna with one section is illustrated in Figure 7, and ( ) *<sup>i</sup> E z <sup>z</sup>* is the tangential component in the incident electric field. The corresponding one-wire equivalent transmission

analyzed. In this case, the mode decomposition of the incident field is eliminated. A single straight antenna is a continuous straight-line antenna, which could be one section or multiple non-parallel sections. It can be equivalent to the one-wire equivalent transmission line model. Each section corresponds to a potential and a current equation. A single straight antenna with one section is illustrated in Figure 7, and ( ) *<sup>i</sup> E z <sup>z</sup>* is the tangential component in the incident electric field. The corresponding one-wire equivalent transmission line model is shown in Figure 8. For the transmission line excited by the distribution volt-

age source, the electric potential and current on the line are governed by the transmission

**Figure 7.** Single straight antenna exposed to incident electric field. **Figure 7.** Single straight antenna exposed to incident electric field. **Figure 7.** Single straight antenna exposed to incident electric field.

**Figure 8.** One-wire equivalent transmission line model of Figure 7.

**Figure 8.** One-wire equivalent transmission line model of Figure 7**.** Compared with the two-wire equivalent transmission line model, the one-wire **Figure 8.** One-wire equivalent transmission line model of Figure 7**.** Compared with the two-wire equivalent transmission line model, the one-wire model is a more general method. It can be used for any antenna composed of relatively Compared with the two-wire equivalent transmission line model, the one-wire model is a more general method. It can be used for any antenna composed of relatively narrow strips or slots. This will benefit the analysis of FSS, where the element currents are usually unknown. The normalized current distribution of Figure 7 is shown in Figure 9, together with the results of the two-wire model.

model is a more general method. It can be used for any antenna composed of relatively narrow strips or slots. This will benefit the analysis of FSS, where the element currents are

narrow strips or slots. This will benefit the analysis of FSS, where the element currents are

usually unknown. The normalized current distribution of Figure 7 is shown in Figure 9,

**Figure 9.** Normalized current distribution of one-wire transmission line model. 1 = 0, = 1 × 10−7 m, = 10 GHz, = 1.6, = 0° . (**a**) = 0.25; (**b**) = 0.5. **Figure 9.** Normalized current distribution of one-wire transmission line model. *l*<sup>1</sup> = 0, *<sup>ρ</sup>* <sup>=</sup> <sup>1</sup> <sup>×</sup> <sup>10</sup>−<sup>7</sup> m, *<sup>f</sup>* <sup>=</sup> 10 GHz, *<sup>ε</sup><sup>r</sup>* <sup>=</sup> 1.6, *<sup>θ</sup>* <sup>=</sup> <sup>0</sup> ◦ . (**a**) *l* = 0.25*λ*; (**b**) *l* = 0.5*λ*.

#### **5. Loss of Equivalent Transmission Line 5. Loss of Equivalent Transmission Line**

antenna

where 120 / *Zw r* =

π ε

The equivalent transmission line models discussed above are lossless. In this section, the antenna lossy is analyzed. In this case, the current in the receiving antenna needs to be obtained first and the per unit length radiation resistance can then be calculated. However, the current is practically an unknown quantity. One way to obtain it is to approximate the current in the form of sine function. However, it is restricted to some special cases (e.g., electrically small antenna and vertical incident). An alternative way is to use the lossless receiving current *I z* <sup>0</sup> ( ) as the approximation, which has been discussed in Sections 3 and 5. *I z* <sup>0</sup> ( ) produces the scattered fields, namely the re-radiated field. The radiation energy is the radiation lossy of the antenna. Consequently, the radiation resistance, potential and current of the equivalent transmission line model can then be solved. Through the iterative algorithm of the above process, a more accurate current distribution The equivalent transmission line models discussed above are lossless. In this section, the antenna lossy is analyzed. In this case, the current in the receiving antenna needs to be obtained first and the per unit length radiation resistance can then be calculated. However, the current is practically an unknown quantity. One way to obtain it is to approximate the current in the form of sine function. However, it is restricted to some special cases (e.g., electrically small antenna and vertical incident). An alternative way is to use the lossless receiving current *I*0(*z*) as the approximation, which has been discussed in Sections 3 and 5. *I*0(*z*) produces the scattered fields, namely the re-radiated field. The radiation energy is the radiation lossy of the antenna. Consequently, the radiation resistance, potential and current of the equivalent transmission line model can then be solved. Through the iterative algorithm of the above process, a more accurate current distribution of the receiving antenna can be obtained.

of the receiving antenna can be obtained. Suppose *r* is parallel to <sup>1</sup> *r* , θ θ = <sup>1</sup> . When calculating the magnitude of the field, assume <sup>1</sup> *r r* ≈ . When calculating the phase of a field, assume <sup>1</sup> *r rz* ≈ − cosθ . The far-field Suppose *r* is parallel to *r*1, *θ* = *θ*1. When calculating the magnitude of the field, assume *r*<sup>1</sup> ≈ *r*. When calculating the phase of a field, assume *r*<sup>1</sup> ≈ *r* − *z* cos *θ*. The far-field radiation produced by the current element is

$$d\mathbb{E}\_{\theta} \approx j \frac{\mathfrak{A}0\mathfrak{f}}{\sqrt{\varepsilon\_r}} \frac{\mathbf{I}(z)dz}{r} \sin\theta e^{-j\mathfrak{f}(r-z\cos\theta)}\tag{18}$$

*r* The far-field radiation electric field of the antenna is the integral of (18) on the whole The far-field radiation electric field of the antenna is the integral of (18) on the whole antenna

ε

$$E\_{\theta} = \int\_{-l}^{l} dE\_{\theta} dz = j \frac{\Im 0 \not\!\!/ }{\sqrt{\varepsilon\_{r}}} \frac{e^{-j\beta r}}{r} \sin \theta \int\_{-l}^{l} I(z) e^{j\beta z \cos \theta} dz. \tag{19}$$
 
$$\dots \quad . \qquad . \qquad . \qquad . \qquad . \qquad . \qquad . \qquad . \qquad . \qquad . \qquad . \qquad . \qquad . \qquad . \qquad . \qquad . \qquad . \qquad .$$

*r* ε According to the approximation property of the plane wave in the far field, the active According to the approximation property of the plane wave in the far field, the active power of the antenna is

*l l*

$$P\_r = \frac{1}{2} \oint\_S \left[ \stackrel{\rightarrow}{E} \times \stackrel{\rightarrow}{H}^\* \right] \cdot d\stackrel{\rightarrow}{s} = \frac{1}{2} \int\_0^{2\pi} \int\_0^{\pi} \frac{1}{Z\_\theta} |E\_\theta|^2 r^2 \sin\theta d\theta d\varphi \tag{20}$$

$$\text{where}$$

$$R = \text{erf} \qquad \text{erf} \qquad \text{erf} \qquad \text{or} \qquad \text{erf} \qquad \text{or} \qquad \text{(20)} : \text{\textbullet (20)} : \text{\textbullet 11}$$

where *Z<sup>w</sup>* = 120*π*/ √ *εr* . Substituting (19) into (20) yields

$$P\_r = \frac{5\pi\beta^2}{\sqrt{\varepsilon\_r}} \int\_0^\pi \left| \int\_{-l}^l I(z)e^{j\beta z \cos\theta} dz \right|^2 \sin^3\theta d\theta \tag{21}$$

*r* The total lossy on the transmission line (Figure 10, the average actual power) is

$$P\_r' = \frac{1}{2} \int\_{-l}^{l} |I(z)|^2 R\_1 dz \tag{22}$$

ε

*R*1

The total lossy on the transmission line (Figure 10, the average actual power) is

<sup>1</sup> ( ) <sup>2</sup> *l <sup>r</sup> <sup>l</sup> P I z R dz* <sup>−</sup>

( )

*l l*

−

∫

<sup>15</sup> sin *<sup>l</sup> <sup>r</sup> j z l*

0

π β

∫ ∫

*r*

ε−

2 2

( )

β θ

*I z dz*

cos 3

*I z e dz d*

Suppose that all radiation lossy of an antenna are equivalent to the total power dissipated on the per unit length radiation resistance *R*<sup>1</sup> of the transmission line, we can solve

2 1

′ <sup>=</sup> ∫ (22)

θ θ

(23)

**Figure 10.** The radiation field of the receiving antenna current. **Figure 10.** The radiation field of the receiving antenna current.

1

=

*R*

Consequently, the equivalent transmission line model of Figures 2, 5 and 8 should be revised. Taking Figure 7 as an example, the lossy is shown in Figure 11. The corresponding Suppose that all radiation lossy of an antenna are equivalent to the total power dissipated on the per unit length radiation resistance *R*<sup>1</sup> of the transmission line, we can solve *R*<sup>1</sup>

$$R\_1 = \frac{\frac{15\beta\_r^2}{\sqrt{\epsilon\_r}} \int\_0^\pi \left| \int\_{-l}^l I(z)e^{j\beta z \cos\theta} dz \right|^2 \sin^3\theta d\theta}{\int\_{-l}^l |I(z)| dz} \tag{23}$$

( ) ( ) *I z jwC V z <sup>z</sup>* <sup>∂</sup> = ⋅ <sup>∂</sup> (25) Consequently, the equivalent transmission line model of Figures 2, 5 and 8 should be revised. Taking Figure 7 as an example, the lossy is shown in Figure 11. The corresponding transmission line equations and parameters should be revised as well

ω

$$\frac{\partial V(z)}{\partial z} = -(jwL + \mathcal{R}\_1) \cdot I(z) + E\_z^i(z) \tag{24}$$

 ω

$$\frac{\partial I(z)}{\partial z} = jw\mathbf{C} \cdot V(z) \tag{25}$$

$$Z\_t = \sqrt{(j\text{u}\mathbb{L} + \mathbb{R}\_1)/j\text{u}\mathbb{C}}\tag{28}$$

0 5 10 15 20

Transmitting antenna

Two-wire even mode Two-wire odd mode

One-wire

Frequency (GHz)

$$r = j\beta\_r = \sqrt{(jwL + R\_1)j\omega\mathcal{C}} = j\omega\sqrt{\mu\_0\varepsilon\_0\varepsilon\_r} \tag{27}$$

β

**Figure 11.** One-wire equivalent lossy transmission line model of Figure 7.

The simulation results show that the lossy has little influence on the normalized cur-

**Figure 11.** One-wire equivalent lossy transmission line model of Figure 7. The simulation results show that the lossy has little influence on the normalized current distribution. The per unit length radiation resistance *R*<sup>1</sup> is shown in Figure 12.

0

10000

20000

Resistance

The equivalent transmission line model is established for the current analysis of line

antenna. The two-wire model is discussed for transmitting antenna and receiving antenna (the incident field should be decomposed for asymmetric case). The one-wire model is introduced for any receiving antenna composed of narrow strips. The corresponding mode decomposition is eliminated. The transmission line radiation resistance hardly ef-

**Author Contributions:** Conceptualization, S.W. and W.W.; methodology, S.W.; software, S.W.; validation, W.W.; formal analysis, W.W.; investigation, W.W.; resources, S.W.; data curation, S.W.; writing—original draft preparation, S.W.; writing—review and editing, S.W.; visualization, S.W.; supervision, W.W.; project administration, W.W.; funding acquisition, W.W. All authors have read

**Funding:** This work was supported by the key project of Science and Technology Innovation Pro-

gram of Army Engineering University of PLA under Grant KYCQJQZL2119.

30000

40000

rent distribution. The per unit length radiation resistance *R*<sup>1</sup> is shown in Figure 12.

(**a**) (**b**)

0 5 10 15 20

Transmitting antenna

Two-wire even mode Two-wire odd mode

One-wire

Frequency (GHz)

fects the normalized current distribution.

and agreed to the published version of the manuscript..

basis function method and spectral rotation approaches. *IEEE Antennas Wirel. Propag. Lett.* **2013**, *12*, 1404–1407.

**Conflicts of Interest:** The authors declare no conflict of interest.

1. Orta, R.; Tascone, R.; Zich, R. A unified formulation for the analysis of general frequency selective surfaces. *Electromagnetics*

2. Pelletti, C.; Bianconi, G.; Mittra, R.; Monorchio, A. Analysis of finite conformal frequency selective surfaces via the characteristic

**Figure 12.** Per unit length radiation resistance. (**a**) = 0.25; (**b**) = 0.5.

0

**6. Conclusions**

**References**

**2007**, *5*, 307–329.

10000

20000

Resistance

30000

40000

50000

The simulation results show that the lossy has little influence on the normalized cur-

rent distribution. The per unit length radiation resistance *R*<sup>1</sup> is shown in Figure 12.

*I z*( ) *I z dI z* ( ) + ( ) ( ) *<sup>i</sup> E z dz <sup>z</sup> Rdz j Ldz*

*V z*( ) *V z dV z* ( ) <sup>+</sup> ( )

ω

*z l* =

**Figure 12.** Per unit length radiation resistance. (**a**) = 0.25; (**b**) = 0.5. **Figure 12.** Per unit length radiation resistance. (**a**) *l* = 0.25*λ*; (**b**) *l* = 0.5*λ*.

*j Cdz* ω

**Figure 11.** One-wire equivalent lossy transmission line model of Figure 7.

*z l* = − − +

#### **6. Conclusions 6. Conclusions**

The equivalent transmission line model is established for the current analysis of line antenna. The two-wire model is discussed for transmitting antenna and receiving antenna (the incident field should be decomposed for asymmetric case). The one-wire model is introduced for any receiving antenna composed of narrow strips. The corresponding mode decomposition is eliminated. The transmission line radiation resistance hardly effects the normalized current distribution. The equivalent transmission line model is established for the current analysis of line antenna. The two-wire model is discussed for transmitting antenna and receiving antenna (the incident field should be decomposed for asymmetric case). The one-wire model is introduced for any receiving antenna composed of narrow strips. The corresponding mode decomposition is eliminated. The transmission line radiation resistance hardly effects the normalized current distribution.

**Author Contributions:** Conceptualization, S.W. and W.W.; methodology, S.W.; software, S.W.; validation, W.W.; formal analysis, W.W.; investigation, W.W.; resources, S.W.; data curation, S.W.; writing—original draft preparation, S.W.; writing—review and editing, S.W.; visualization, S.W.; supervision, W.W.; project administration, W.W.; funding acquisition, W.W. All authors have read and agreed to the published version of the manuscript.. **Author Contributions:** Conceptualization, S.W. and W.W.; methodology, S.W.; software, S.W.; validation, W.W.; formal analysis, W.W.; investigation, W.W.; resources, S.W.; data curation, S.W.; writing—original draft preparation, S.W.; writing—review and editing, S.W.; visualization, S.W.; supervision, W.W.; project administration, W.W.; funding acquisition, W.W. All authors have read and agreed to the published version of the manuscript..

**Funding:** This work was supported by the key project of Science and Technology Innovation Program of Army Engineering University of PLA under Grant KYCQJQZL2119. **Funding:** This work was supported by the key project of Science and Technology Innovation Program of Army Engineering University of PLA under Grant KYCQJQZL2119.

**Conflicts of Interest:** The authors declare no conflict of interest. **Conflicts of Interest:** The authors declare no conflict of interest.

#### **References References**

