This paper proposes a reduced-order model for the large scale circuits representing the multiconductor systems. It is based on the block Arnoldi algorithm for calculating congruence transformation matrix. By setting up the state equations in frequency domain for the multiconductor systems, the transfer functions are calculated and so the frequency response curves are obtained. A comparison is made between those curves obtained from the circuits with and without reduced-order treatment and a better agreement appears between them. By using inverse Laplace transform, the lightning transient responses in the multiconductor systems are given in time-domain. A better agreement is shown between calculated and experimental results, which confirm the validity of the proposed model.

1. Introduction

An external lightning protection system of a structure is usually constituted by numerous interconnected longitudinal and transverse conducting branches. It takes a form of three-dimensional multiconductor system and can be converted into an equivalent passive RLC network. For an actual tall structure, its external lightning protection system is large in physical dimension and complex in construction, which results in a large scale equivalent network. From the view point of lightning protection design, the need exists for obtaining lightning transient responses of the equivalent network. For the sake of saving storage capacity and computation time, a simplification should be made for the equivalent network; otherwise the circuit calculation is difficult to perform. From the existing research on lightning transient responses [16], we have almost not found the relevant procedure that can give a great simplification for the large scale equivalent network. Therefore, the order reduction model is necessary for the transient analysis of lightning protection systems. Recently, the investigation on passive network synthesis has been reported in literature [710], which is closely related to the order reduction model. The investigation work is of practical significance in analysis of large scale electrical networks. Based on the achievement of the existing work, a reduced-order algorithm (ROA) is introduced in this paper as a Krylov space method which provides a passive, stable, and accurate macromodel compared with the reduction techniques [11, 12]. In order to reduce the matrices directly in the system state equation, the transformation matrix can be calculated by using the block Arnoldi algorithm [13]. In the ROA, we have developed a robust framework by unveiling the connections between the passive reduction technique and other well known Krylov space based processes. In terms of the reduced-order model of the structure under study, amplitude-frequency and phase-frequency curves can be quantitatively calculated and so we can compare those curves with the frequency response curves of the original circuit system. Furthermore, the time function can be obtained with the help of the inverse Laplace transform and the lightning transient response is given by using the convolution theorem. The comparison is also made between calculated and measured results to confirm the validity of the proposed model.

2. State Equations of the Equivalent Circuit

Figure 1 shows the configuration of a cage-like multiconductor system formed by vertical and horizontal conducting branches. For setting up the state equations of the equivalent network, the electrical parameters of the conducting branches in the multiconductor system are represented by resistances, inductances, and capacitances. Since the distance between the two conducting branches is much greater than their respective radius, the electromagnetic couplings between them can be ignored. Therefore, the formulas for evaluating the electrical parameters are derived as follows.

Consider a horizontal conducting branch in the multiconductor system, as shown in Figure 2. Its capacitance can be expressed by [14] where .

Similarly, the capacitance of a vertical conducting branch, as shown in Figure 3, can be expressed by [14] where the coefficient is According to the electromagnetic analogy [15], the product of capacitance and inductance is equal to that of space permeability and permittivity . This gives the inductances for horizontal and vertical conducting branches Additionally, the conductance of an arbitrary conducting branch is estimated by where is the resistivity of conducting branch and the subscript takes as or .

In consideration of the propagation phenomenon of lightning current, each conducting branch of the multiconductor system is subdivided into a suitable number of segments. The length of each segment is assumed to be less than or equal to one-tenth of the wavelength corresponding to the maximum frequency likely to affect the system transient [16]. Each segment is represented by a π-type circuit, as shown in Figure 4, where and are the segment numbers of horizontal and vertical branches, respectively.

As a result, the multiconductor system is converted into an equivalent network constituting a large number of π-type circuits. On the basis of the electrical parameters calculated from , the capacitance, conductance, and inductance matrices can be formed for the equivalent network: where the subscripts , , and denote the total numbers of capacitances, conductances, and inductances of the equivalent network, respectively. The circuit components of the equivalent network are numbered in the sequence of capacitances, inductances, and lightning current source. In accordance with the component number, the incidence matrix of the equivalent network is written as The node voltages (i.e., the voltages across the capacitances) and the currents through the inductances are selected as the state variables: where the subscript denotes the total node number of the equivalent network. Thus, the state equation is set up for the equivalent network: where where and . Note that , , and are symmetric and nondefinite matrices.

Next, the so-called transfer function can be calculated, which describes the input-output behavior of a linear dynamical system in frequency domain [17]. Taking the Laplace transformation of , the transfer function is given by

The computation of is called simulation in frequency domain and of course not done by calculating . However, it is still expensive for large system. For this reason, the reduced-order model is necessarily proposed in this paper.

3. Reduced-Order Model

For simplifying the lightning transient calculation of the large scale RLC network, the ROA is employed here. The algorithm is based on the block Arnoldi method. Therefore, we introduce our notion of block Krylov subspaces for multiple starting vectors.

3.1. Model Order Reduction Techniques

A widely used class of model order reduction techniques is the class of Krylov subspace methods. Before block Krylov subspace methods are used, the two matrices and in have to be reduced, denoted by in the following. This can be done by rewriting as where Although and are sparse matrices, in general, the matrix is a dense matrix. However, block Krylov subspace methods involve only in the form of matrix-vector products and possibly [18].

The proper definition of block Krylov subspaces is necessary to be involved and the use of these subspaces results in much more efficient reduced-order modeling techniques. The block Krylov space generated by matrices and is defined as

3.2. The Algorithm of Calculating Congruence Transformation Matrix

The congruence transformation matrix is calculated based on the above definition of . Here, is an orthonormal matrix and subjected to ( is an identity matrix) [19]. The classical Arnoldi process [20] generates orthonormal basis vectors for the sequence of Krylov subspaces , , induced by and . In order to saving storage capacity and computation time, a smaller value of the order is better. In fact, exact deflation at step of the Arnoldi-type process occurs if and only if . Similarly, inexact deflation occurs if and only if , but . Therefore, the order will be determined according to , where is a suitably chosen deflation tolerance. A simple Arnoldi-type algorithm is presented in the following [18].(0)Set ().Set .For do the following:(1)compute and check if the deflation criterion () is fulfilled.If yes, is deflated by doing the following.Set . If , set and stop.Set for .Return to step .(2)Set and .(3)Compute .(4)For do the following:set and .(5)For do the following:set and .

3.3. A Reduced-Order Algorithm

After computing the transformation matrix , we can easily find the ROA matrices. Applying the change of variable in and multiplying both sides by give For the macromodel, the reduced-order matrices are

These types of transformations are known as congruence transformations. From and , the reduced-order transfer function is given by Define The reduced system can be formulated as

By applying any eigendecomposition routine on , we can obtain where is a diagonal matrix. Since is a real matrix, any complex eigenvalue or eigenvector has its conjugate. Therefore, we can use real and matrices in the following transformations: where and .   is the transformation matrix that can be defined as follows:if is real, ;if , . will contain 2 × 2 blocks on the diagonal for complex eigenvalues, and 1 × 1 blocks for the real eigenvalues. By substituting into , we can get The final transfer function is given by As a result, we can obtain the final reduced-order model in the following form: where , , , and are the dominant pole and corresponding residue, respectively, is the number of reduce-order, and represents any direct coupling between ports.

4. Time-Domain Simulation of the Reduced-Order Model

For a complete simulation of circuit response, the linear elements should be simulated along with the reduced-order macromodels. Here, we describe a method to include the ROA macromodel. The method is based on a -parameter description of the macromodel [21]. In order to calculate the -parameters of the reduced-order system, the eigendecomposition steps [ and ] are used. After obtaining the final form of the transfer function , the time function can be calculated by using inverse Laplace transform of partial fraction expansion.

According to , we need to transform the rational fraction into proper fraction. So, we can get the following: where the remainder term is a proper fraction.

The denominator polynomials can be factored via partial fraction expansion, and the roots of the equation are solved. The roots of the equations have two kinds of situation, such as simple root and conjugate complex.

(i) If the equation has simple roots, we can set the simple roots as . Therefore, the expansion form of is where , .

After determining the coefficients of , the time function is

(ii) If the equation has two conjugate complex roots (, ), we thus have We set and ; thus

Subsequent to the calculation of the time function , convolution is needed for time-domain analysis. Hence, the circuit response in time-domain at the ports can be given by which requires complexity, where is the number of time points during simulation. For this reason, recursive convolution [22] and time-domain -parameter macromodels [21] were developed, where the complexity is linear with the number of time points. The algorithms of these approaches are used by MATLAB in detail.

5. Numerical Example

A reduced-scale structure of lightning protection system is shown in Figure 5. The diameter of each steel bar is 8 mm. According to the above model and algorithm, the programs are written by MATLAB. Table 1 summarizes the reduction results of matrix dimension and the computational time by using ROA.

In order to verify the validity of the reduced-order model proposed above, the amplitude-frequency curves and phase-frequency curves of original and reduced-order model are computed, as shown in Figures 6 and 7.

From the two figures above, we can see that a better agreement appears between original and reduced-order model.

For the sake of computing the circuit response of the reduced-scale model, the injection lightning current needs to be known. The waveform of the current source is shown in Figure 8. Therefore, the circuit response of reduced-order model is computed by the program we have written. The waveforms of the lightning transient responses and can be obtained and they are compared with that measured experimentally [23] as shown in Figures 9 and 10. It can be seen from Figures 8 and 9 that very good agreement is also found between computed and measured results.

6. Conclusions

A novel reduced-order macromodeling algorithm has been proposed for calculating the lightning transient responses in the lightning protection systems of structures. The reduced-order model given here uses the block Arnoldi algorithm to reduce the matrix dimension. On the basis of matrix operation and discretization, the amplitude-frequency curves and phase-frequency curves can be obtained by means of the MATLAB program. These curves under including and excluding the reduced-order treatment are compared with each other and a better agreement is shown between them. The inverse Laplace transform of partial fraction expansion has also been performed for time-domain analysis. The computed lightning transient responses in a lightning protection system conform reasonably with those from experimental measurement, which confirms the validity of the proposed model.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This work was financially supported by the National Natural Science Foundation of China under Contract no. 50977002.