Research Article  Open Access
Mathematical Modeling of Transient Responses in a Large Scale Multiconductor System
Abstract
This paper proposes a reducedorder 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 reducedorder treatment and a better agreement appears between them. By using inverse Laplace transform, the lightning transient responses in the multiconductor systems are given in timedomain. 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 threedimensional 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 [1–6], 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 [7–10], 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 reducedorder 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 reducedorder model of the structure under study, amplitudefrequency and phasefrequency 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 cagelike 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 onetenth 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 socalled transfer function can be calculated, which describes the inputoutput 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 reducedorder model is necessarily proposed in this paper.
3. ReducedOrder 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 matrixvector 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 reducedorder 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 Arnolditype 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 Arnolditype 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 ReducedOrder 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 reducedorder matrices are
These types of transformations are known as congruence transformations. From and , the reducedorder 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 reducedorder model in the following form: where , , , and are the dominant pole and corresponding residue, respectively, is the number of reduceorder, and represents any direct coupling between ports.
4. TimeDomain Simulation of the ReducedOrder Model
For a complete simulation of circuit response, the linear elements should be simulated along with the reducedorder 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 reducedorder 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 timedomain analysis. Hence, the circuit response in timedomain 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 timedomain 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 reducedscale 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 reducedorder model proposed above, the amplitudefrequency curves and phasefrequency curves of original and reducedorder 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 reducedorder model.
For the sake of computing the circuit response of the reducedscale 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 reducedorder 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 reducedorder macromodeling algorithm has been proposed for calculating the lightning transient responses in the lightning protection systems of structures. The reducedorder model given here uses the block Arnoldi algorithm to reduce the matrix dimension. On the basis of matrix operation and discretization, the amplitudefrequency curves and phasefrequency curves can be obtained by means of the MATLAB program. These curves under including and excluding the reducedorder 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 timedomain 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.
Acknowledgment
This work was financially supported by the National Natural Science Foundation of China under Contract no. 50977002.
References
 K. J. Kerns, I. L. Wemple, and A. T. Yang, “Stable and efficient reduction of substrate model networks using congruence transforms,” in Proceedings of the 1995 IEEE/ACM International Conference on ComputerAided Design (ICCAD '95), pp. 207–214, November 1995. View at: Google Scholar
 L. T. Pillage and R. A. Rohrer, “Asymptotic waveform evaluation for timing analysis,” IEEE Transactions on ComputerAided Design of Integrated Circuits and Systems, vol. 9, no. 4, pp. 352–366, 1990. View at: Publisher Site  Google Scholar
 P. Feldmann and R. W. Freund, “Efficient linear circuit analysis by Pade approximation via the Lanczos process,” IEEE Transactions on ComputerAided Design of Integrated Circuits and Systems, vol. 14, no. 5, pp. 639–649, 1995. View at: Publisher Site  Google Scholar
 L. M. Silveira, M. Kamon, and J. White, “Efficient reducedorder modeling of frequencydependent coupling inductances associated with 3D interconnect structures,” in Proceedings of the 32nd Design Automation Conference (DAC '95), pp. 376–380, June 1995. View at: Google Scholar
 E. J. Grimme, Krylov projection methods for model reduction [Ph.D. dissertation], University of Illinois at UrbanaChampaign, Champaign, Ill, USA, 1997.
 D. L. Boley, “Krylov space methods on statespace control models,” Circuits, Systems, and Signal Processing, vol. 13, no. 6, pp. 733–758, 1994. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Z. Q. Chen, K. Wang, M. Yin, C. Li et al., “Synthesis of nport resistive networks containing 2n terminals,” Int. J. Circ. Theor. Appl, 2013. View at: Google Scholar
 K. Wang, M. Z. Q. Chen, and Y. Hu, “Synthesis of, biquadratic impedances with at most four passive elements,” Journal of the Franklin Institute, vol. 351, pp. 1251–1267, 2014. View at: Google Scholar
 M. Z. Q. Chen, K. Wang, Z. Shu, and C. Li, “Realizations of a special class of admittances with strictly lower complexity than canonical forms,” IEEE Transactions on Circuits and Systems. I. Regular Papers, vol. 60, no. 9, pp. 2465–2473, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 M. Z. Q. Chen, K. Wang, Y. Zou, and J. Lam, “Realization of a special class of admittances with one damper and one inerter for mechanical control,” IEEE Transactions on Automatic Control, vol. 58, no. 7, pp. 1841–1846, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 A. Odabasioglu, M. Celik, and L. T. Pileggi, “PRIMA: passive reducedorder interconnect macromodeling algorithm,” IEEE Transactions on ComputerAided Design of Integrated Circuits and Systems, vol. 17, no. 8, pp. 645–654, 1998. View at: Publisher Site  Google Scholar
 A. Odabasioglu, M. Celik, and L. T. Pileggi, “Practical considerations for passive reduction of RLC circuits,” in Proceedings of the 1999 IEEE/ACM International Conference on ComputerAided Design (ICCAD '99), pp. 214–219, November 1999. View at: Google Scholar
 R. W. Freund, “Reducedorder modeling techniques based on Krylov subspaces and their use in circuit simulation,” in Applied and Computational Control, Signals, and Circuits, Vol. 1, Appl. Comput. Control Signals Circuits, pp. 435–498, Birkhäuser, Boston, Mass, USA, 1999. View at: Google Scholar  MathSciNet
 U. Y. Yosseli, E. S. Kochanov, and M. G. Stryncki, Calculation of Capactances, Energy Press, Moscow, Russian, 1981.
 C. Z. Feng, Electromagnetic Field Theory, Higher Education Press, Beijing, China, 2002.
 S. Cristina and A. Orlandi, “Calculation of the induced effects due to a lightning stroke,” IEE Proceedings B, vol. 139, no. 4, pp. 374–380, 1992. View at: Publisher Site  Google Scholar
 M. Günther, M. Hoschek, and P. Rentrop, “Differentialalgebraic equations in electric circuit simulation,” International Journal of Electronics and Communications, vol. 57, pp. 101–107, 2000. View at: Google Scholar
 R. W. Freund, “Krylovsubspace methods for reducedorder modeling in circuit simulation,” Journal of Computational and Applied Mathematics, vol. 123, no. 12, pp. 395–421, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 K. J. Kerns and A. T. Yang, “Stable and efficient reduction of large, multiport RC networks by pole analysis via congruence transformations,” IEEE Transactions on ComputerAided Design of Integrated Circuits and Systems, vol. 16, no. 7, pp. 734–744, 1997. View at: Publisher Site  Google Scholar
 W. E. Arnoldi, “The principle of minimized iteration in the solution of the matrix eigenvalue problem,” Quarterly of Applied Mathematics, vol. 9, pp. 17–29, 1951. View at: Google Scholar  MathSciNet
 S. Y. Kim, N. Gopal, and L. T. Pillage, “Timedomain macromodels for VLSI interconnect analysis,” IEEE Transactions on ComputerAided Design of Integrated Circuits and Systems, vol. 13, no. 10, pp. 1257–1270, 1994. View at: Publisher Site  Google Scholar
 V. Raghavan, J. E. Bracken, and R. A. Rohrer, “AWESpice: a general tool for the accurate and efficient simulation of interconnect problems,” in Proceedings of the 29th ACM/IEEE Design Automation Conference, pp. 87–92, June 1992. View at: Google Scholar
 I. A. Metwally, F. H. Heidler, and W. J. Zischank, “Magnetic fields and loop voltages inside reduced and fullscale structures produced by direct lightning strikes,” IEEE Transactions on Electromagnetic Compatibility, vol. 48, no. 2, pp. 414–426, 2006. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Caiwang Sheng and Xiaoqing Zhang. 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.