Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2016 (2016), Article ID 7870462, 13 pages
Research Article

Dynamic Characterization of Typical Electrical Circuits via Structural Properties

Facultad de Ingeniería, UNAM, Edificio de Posgrado, Segundo Piso, Ciudad Universitaria, 04510 Ciudad de México, Mexico

Received 6 March 2016; Revised 22 July 2016; Accepted 31 July 2016

Academic Editor: Carlo Cosentino

Copyright © 2016 Sofía Avila-Becerril 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.


The characterization of a class of electrical circuits is carried out in terms of both stability properties and steady-state behavior. The main contribution is the interpretation of the electrical topology (how the elements that conform the circuits are interconnected) in terms of mathematical properties derived from the structure of their models. In this sense, at what extent the topology by itself defines the dynamic behavior of the systems is explained. The study is based on the graph theory allowing capturing, departing from the well-known Kirchhoff laws, the topology of the circuits into several matrices with specific structure. The algebraic analysis of these matrices permits identifying conditions that determine whether the system is stable in the sense of Lyapunov and the kind of steady-state behavior that it exhibits. The approach is mainly focused on typical topologies widely used in practice, namely, radial, ring, and mesh networks.

1. Introduction

Mathematical characterization of electrical circuits is a topic with a very long history. From the celebrated Kirchhoff laws, the maturity achieved in this field is widely recognized and reflected on classical references like [1, 2], among a vast number of them.

Fortunately, the richness of the available knowledge has become in an impelling force that continually opens new perspectives to take advantage of the known theory [35]. Indeed, modern applications, which range from microcircuits found in integrated circuit design to complex macronetworks like electrical power systems, are approached using the solid results formulated to characterize this kind of system.

One particular research activity that has received great attention is the control of electrical networks, whose interest is to modify the behavior of the circuit with the objective that some variables attain a prescribed value. In this context, several results have been reported about the reformulation of well-known features of circuits or even the establishment of new properties. Thus, it is possible to find results covering from characterization [6, 7] to control [8] and applications that include both small [5] and large [9] size systems.

Another quite interesting approach to control electrical circuits is based on the modification of their structure to achieve a given operation. Viewing the control problem from this perspective is not new; for example, in the field of electrical power systems [10], both voltage and frequency stability problems are usually confronted by adding to the network new elements, called compensators, adequately located. This compensation action means the addition of new capacitors, inductors, and resistors in specific locations and in such a way that the desired behavior of the power system is achieved. Under the latter approach, several recent results have been reported whose interest lies in exploiting the structure of a given circuit to reach a specific goal. For example, in [11], the problem of shaping a resistive circuit behavior through the interconnection of another resistive circuit is studied, leading to a methodology denoted by “partial synthesis by interconnection” which, in turn, belongs to a more general controller design approach called Control by Interconnection (CbI) design [12]. Roughly speaking, under this perspective, the objective is to look at the controller as one dynamical system that when interconnected with other ones (the plant) generates a new dynamical system with desired properties. Considering this scenario, if it is well understood how the dynamic behavior of a given circuit is determined by its structure, then this knowledge can be used to redesign or to control the system in order to satisfy some prescribed specifications.

The aim of this paper is to provide an answer to the question stated above for the class of networks most used in practice, namely, radial, ring, and mesh topologies, which in the sequel will be identified as typical electrical circuits. The main objective is to recast both the stability properties and the steady-state behavior exhibited by these networks in terms of structural features of their mathematical models, in particular the matrices involved in them which are characterized by the interconnection of the different elements that compose the circuit.

From a technical perspective, the identification of the aforementioned properties is carried out using arguments from graph and matrix theory [13, 14] while the stability properties are formulated in terms of Lyapunov theory. Concerning the steady-state behavior, it is considered as the value to which the trajectories of the system tend as time tends to infinity; that is, these can be time-varying or constant (equilibrium points) trajectories.

In order to present the contributions of the paper, Kirchhoff’s laws of a generic circuit are first formulated in terms of basic cutsets and loopsets [15]. After this, the fundamental concept of admissible trajectories that refers to the dynamic behavior that can be exhibited by the system is introduced. With this knowledge, the stability analysis of these behaviors is developed going from simple (linear) to more complex (nonlinear) structures, placing special attention on the typical networks. The last part of the paper concerns the characterization of the model components that influences the steady-state behavior of the typical topologies. A couple of illustrative case studies are included to illustrate the usefulness of the presented results. These analyses are related to the operation often found in practice given by a Mesh network equipped with Direct Current (DC) and Alternate Current (AC) voltage sources. For the former, some consensus operation [16] is concluded, while for the latter it is evidenced that compensation procedures can be developed in a systematic way.

The rest of the paper is organized as follows. In Section 2, using well-known results from the graph theory, Kirchhoff’s laws are stated, in terms of basic cutsets and loopsets, for a model that includes a broad class of electrical circuits, while Section 3 is devoted to the presentation of its dynamic model. In Section 4, the stability analysis is carried out and the bases for the characterization of the steady-state behavior of typical networks are included in Section 5. The usefulness of the presented results is illustrated in Section 6. Some concluding remarks are discussed in Section 7.

2. Electrical Circuit Graphs

An electrical network can be defined as an oriented graph consisting of a finite set of nodes and a finite set of edges such that is a subset of , with pairs of where no self-loops are allowed. In this case, the set of nodes is the interconnection points of elements whereas the set of edges is associated, one to one, with lumped one-port (two-terminal) elements. Additionally, it is considered that the graph is connected in the sense that each node can be reached from any other node by tracing a path through the edges.

With each edge (each lumped one-port element), a voltage across its terminals and a current that flows through it are associated. The orientation of the edge coincides with the direction of a positive current and a decreasing voltage.

Once the lumped elements are interconnected, their port variables must satisfy the constraints stated by the Kirchhoff Current and Voltage Laws (KCL and KVL, resp.) that in this paper are stated in terms of basic cutsets and basic loopsets for a given spanning tree [14, 15]. A basic cutset (or fundamental cutset) is composed of one edge of the tree (branch) and some or all of the edges of the cotree (chords). A basic loopset (or fundamental cycle) is conformed by one chord and some or all of the branches in such a way that a closed loop is formed. In this way, there are two natural matrices associated with the graph, the basic cutset matrix and the basic loopset matrix , so that KCL and KVL are given, respectively, byIf the currents and voltages associated with the edges are ordered in such a way that with and being the currents and voltages of the tree and and being the currents and voltages of the cotree, respectively, it is possible to write down the network constraints [14] aswhere is called the fundamental loop matrix.

From this last expression, it can be noticed that the th row of indicates which of the cotree elements is incident to the th basic cutset and that, in correspondence, the th column shows which of the tree elements belongs to the th basic loopset.

At this point, it is important to state that the matrix determines the topological structure of the system and its dynamic behavior. The analysis of this matrix is the main interest of the paper. All the results presented throughout this contribution are related to its structure.

3. Electrical Circuit Dynamic

In this section, the dynamic equations for the electrical circuit presented in Section 2 are obtained following ideas of, for example, [11]. The analysis is restricted to networks where capacitor-only loops and inductor-only cutsets are not admitted [1]. In addition, the circuit elements are grouped such that voltage sources, all the capacitors, and some resistors appear at the tree, while inductors and the rest of the resistors are in the cotree (For the sake of simplicity, presentation of the current sources will be omitted from the analysis.), leading towith voltage tree and current cotree vectorswhere , , , such that , and , , with .

If the total stored energy of the circuit is defined as , it is well known [15] that the port variables of the capacitors, inductors, and resistors can be represented aswhere and are assumed to be bijective functions.

The partition introduced above results in the fact that matrix , in turn, can be partitioned aswhere the subscript stands for the interconnections between tree and cotree elements.

Remark 1. The partition of matrix presented in (7) is important from the point of view of the analysis developed in this paper. It allows, as will be clear below, identifying which part of the model influences the stability properties of the system and which part determines the steady-state behavior.

As already reported in the literature, substitution of (6a)–(6c) into (3) leads to the dynamical model given byunder the definitionswith matriceswhereand complemented by the algebraic constraint

It is important to notice that if we concentrate on circuits with linear resistive elements, thenwhere , a diagonal matrix with entries in the resistances of the tree resistors, and also diagonal and composed of the resistances of the cotree. As usual, all circuit resistances are assumed to be strictly positive [11]. For this linear representation, it holds thatwhere and are identities matrices and, consequently, models (8) and (15) can be equivalently represented as a port-controlled Hamiltonian system with dissipation [17] of the formdefiningwhere , , and

The following remarks about the obtained model are in order.

Remark 2. The algebraic constraint (13) represents the current demanded to the voltage sources for a given operation point. Along the paper, it is considered that the sources are ideal; that is, they can supply any amount of current, and therefore the constraint is always satisfied.

Remark 3. Notice that the matrix , which includes the relationship between capacitors and inductances, appears in a skew-symmetric form in the model. It will be shown that this condition implies that it does not modify the energy balance equation of the circuit. However, it will be fundamental to determine the steady-state behavior of the circuit.

Remark 4. Also related to matrix , the time derivative of along the trajectories of (16), considering , is given byputting in evidence the fact that stability of the network depends on the matrices , , , and .

Remark 5. If the network is conformed by linear elements, the property above is stronger sinceThus, since , a solution of will be asymptotically stable if the system is zero state detectable considering as output the variable and is a minimum argument of the scalar function .

4. Structural Properties for Stability

The purpose of this section is to identify conditions for stability when is different from zero. In this context, it is necessary to identify the steady-state trajectories that are achievable by the system, denoted by admissible trajectories. These behaviors are solution ofwhere the existence of an input that generates the behavior has been implicitly assumed. A particular case is when is an equilibrium point. Under this condition, is constant and the steady-state operation (equilibrium point) is input dependent and is determined by

The admissible trajectories directly depend on the matrix presented in (7). Therefore, the results presented below concentrate the attention in finding structural properties of this matrix such that some stability properties are guaranteed.

The first result states the stronger result concerning the achievable admissible trajectories although the simplest structure for the network, namely, conditions for , is identified to guarantee tracking of a time-varying solution of (21) but this result applies only when the circuit is composed of elements that exhibit a linear constitutive relationship, that is, when the total stored energy takes the formwith positive and diagonal matrices and of inductances and capacitances, respectively. In this case, it is clear thatwhile the resistors satisfy (14).

Proposition 6. Consider a linear electrical network described by (16) and (23) with being a time-varying input such that its steady-state behavior, characterized byis well posed.
Under these conditions,with if

Proof. Under the linear characteristic assumed for the system, it is possible to describe the error dynamic bywith the energy-like functionand the identities , .
Since the purpose is to analyze the stability properties of the equilibrium point of (28), function (29) can be considered as Lyapunov function candidate. Hence, it is straightforward to show that its time derivative along the trajectories of (28) satisfieswhich can be written in an equivalent way aswith and , being symmetric positive definite matrices defined by whileThe proof is completed by noting that the maximal invariant set where is , which leads to and . Therefore, the necessity to guarantee that are the only solutions that satisfy these constraints brings out as the sufficient conditions that assure asymptotic stability of .

Remark 7. An interesting feature in the proof of the presented result is that the zero state detectability that could be invoked in (31) in order to conclude asymptotic stability has been recast, via (33), in terms of the structural properties of the system. In this scenario, to what extent the stability properties of the network are determined by its structure has been stated.

Remark 8. Even though Proposition 6 is related to linear systems, in the authors’ opinion, it is important since this kind of representation for electrical networks is still widely used under a great variety of applications.

In the proposition below, the result presented above is extended in the sense that the asymptotic stability properties are still guaranteed considering nonlinear capacitors and inductors. However, this is achieved at the expense that the admissible trajectories are no longer time-varying but equilibrium points.

Proposition 9. Consider an electrical network described by model (8) with being a constant input such that the steady-state behavior is well posed. In addition, assume the following: (A.1)The resistors involved in the circuit are characterized by linear constitutive relationships satisfying (15).(A.2)The equilibrium point that corresponds to locally satisfies .Under these conditions, the equilibrium point is locally asymptotically stable if

Proof. In this case, the total stored energy is a nonlinear function. Therefore, under assumption (A.1), the dynamic behavior of the system is described by (16) while its equilibria are characterized by the solutions ofIf assumption (A.2) holds, then, following ideas reported in [18], it is possible to consider the Lyapunov function candidate given bywhose time derivative along the trajectories of (16), under the condition , yieldswith , defined in (17a), where we exploited the fact that, in spite of the nonlinear nature of , the identities introduced in (6a) and (6b) are valid.
The proof is completed following a similar procedure to Proposition 6 since it is possible to get thatwith

Remark 10. Notice that the function qualifies as Lyapunov function candidate only for equilibrium points. Current research is carried out trying to find an alternative analysis to deal with the tracking problem but, to the authors’ knowledge, this imposes an open problem in the related literature.

Remark 11. If is a global minimum of , with radially unbounded, then the equilibrium point is globally asymptotically stable.

Motivated by the identified role played by matrices and in the presented results, the next two propositions are focused on a practically important class of circuits which is characterized by a particular interconnection, first among resistors and inductors and second among resistors and capacitors.

Property 1. If the number of tree resistors is equal to the number of inductors and they are one-to-one series connected, thenwith as already defined and being a zero matrix.

Proof. If there exists the same number of tree resistors as inductors, then . On the other hand, from (3), the constraintmust be satisfied. Hence, the fact that tree resistor currents must be, one-to-one, equal to the inductor currents leads to the results stated in (40).

Property 2. If the number of cotree resistors is equal to the number of capacitors and they are one-to-one parallel connected, thenwith as already defined and being a zero matrix.

Proof. If there exists the same number of cotree resistors as capacitors, then . In this case, also obtained from (3), the constraint that must be satisfied iswhere, due to Property 1, , leading to the expressionThe proof concludes by noting that each entry of must be equal to its corresponding entry in , a condition that holds only if (42) is true.

Remark 12. An immediate implication of the presented properties is that the networks that enjoy them automatically satisfy the conditions stated in Propositions 6 and 9 to guarantee asymptotic stability of the admissible trajectories.

Remark 13. In the context of this paper, the stated properties are important since all the considered typical networks enjoy them, as will be shown in Section 5.

In order to present the last result of this section, writing model (8) under the effect of Properties 1 and 2 is convenient. Hence, the dynamical behavior of the electrical network is described bywhere

Due to the achieved drastic simplification in the model structure, which illustrates how a proper interconnection of the elements determines to a large extent the dynamic behavior of the system, previous results can be extended by relaxing the assumptions about and , which under the new conditions only depend on the system state.

Proposition 14. Consider the electrical network described by (45) with being a constant input such that the steady-state behavior is well posed. Assume (A.2) holds and in addition assume the following: (A.3)The maps and define incremental output strictly passive operators in the sense thathold for and considering as input.Under these conditions, the equilibrium point is locally asymptotically stable.

Proof. The equilibria of the system are characterized byThus, if (A.2) holds, it is possible to consider the Lyapunov function candidate , reported in [18] and defined in (36), which leads, under the condition , toUsing identities (6a) and (6b), this expression can be equivalently written as Since (A.3) holds, it is clear that with maximal invariance set defined aswhich concludes the proof.

Remark 15. It must be clear that if the dissipation functions and define linear mappings, introduction of Properties 1 and 2 does not destroy the results presented either in Proposition 6 or in Proposition 9.

5. Topological Structure of Typical Networks

In this section, the dynamic characterization started in Section 4 is complemented for the class of typical electrical networks. The main objective of the analysis is to explicitly state the structure of the matrix for the aforementioned typical networks. In this sense, the usefulness of this study is twofold:(i)With this structure at hand, it is possible to determine in a systematic way the stability properties and steady-state behavior of the networks.(ii)An interpretation is provided, given via structural characterization, to the widely recognized advantages exhibited by these networks, for example, reliability properties.

In order to present the results mentioned above, first a generic network that captures in a unified way the characteristics of the three typical networks is considered. Once the topology of this network is identified, each particular case is presented.

The first condition that the model of the generic network satisfies is that Properties 1 and 2 stated in Section 4 hold. Hence, its dynamic behavior is represented by model (45) leaving and as the two matrices that remain to be characterized; that is, it is necessary to consider the interconnections between capacitors and inductors, for the former matrix, and sources and inductors, for the latter.

Consider that the inductors are divided into three types, namely, -inductors that belong to a trajectory that connects a source with a capacitor, -inductors that belong to a trajectory that connects a source with another source, and -inductors that belong to a trajectory that connects a capacitor with another capacitor, such that . Moreover, assume that the capacitors belong at least to one of the following classes:(C.1)The th capacitor, , shares cutset with -inductors.(C.2)The th capacitor, , shares cutset with -inductors.However, the voltages source satisfies the following: (C.3)The th voltage source, , shares cutset with inductors.

Under the above classification, the variables associated with capacitors and inductors can be organized such that the rows of matrix can be divided into two parts: the first corresponding to capacitors that hold simultaneously conditions (C.1) and (C.2) and the second corresponding to those capacitors that hold only with condition (C.2). As will be clear below, capacitors that satisfy only condition (C.1) are included as a particular case of the former.

On the other hand, the columns of are divided into three blocks each one corresponding to -, -, and -inductors, respectively. Moreover, the inductors connected to the th capacitor as in (C.1) are grouped for each of the capacitors connected to them.

Following the stated organization, matrix takes the formwhere one has the following: (i)The minus sign appears since, from (3), the constraint that must be satisfied isconsidering that .(ii), , are vectors filled with ones denoting the condition stated in (C.1). In this case, it has been assumed that there exist capacitors of this kind. In addition, it holds that .(iii)The zero columns reflect the fact that any capacitor can be connected to type inductors.(iv)Row vectors include the possibility that a given capacitor can be simultaneously connected to and type inductors. If the th capacitor is connected to one of the -inductors, a appears in the corresponding entry; otherwise, a appears. If the capacitor only holds with (C.1), the vector is the zero vector.(v)In the rows that go from to , only vectors appear since they correspond to (C.2) class capacitors.

Two consequences of structure (52) are presented below.

Property 3. Each column of the matrix is composed of one and one and the rest of the entries equal to zero.

Proof. The columns of matrix denote how the capacitors are connected to other capacitors. Since for each pair of capacitors this connection is carried out by one -inductor, the current leaves one cutset and enters another, which is represented by the in one row and the in the other row.

The second property is actually a corollary of Property 3.

Property 4. The vector , that is, the vector filled with ones of dimension , is a left eigenvector of matrix satisfying .

Proof. The proof can be carried out by direct computation.

Concerning matrix , its structure comes from (13) under the condition ; that is, the constraintmust be satisfied. Thus, the columns of are divided into three blocks corresponding to , , and type inductors, respectively, where the third one is zero due to the fact that sources do not belong to cutsets where -inductors are involved. Hence, this matrix takes the formwhere one has the following: (i)The entries different from zero of each row of stand for the connection of sources with -inductors.(ii)In concordance with the partition of matrix , which in turn comes from (C.1), the th row of , , must be divided into sections leading towhere each , , has only one entry equal to if the th source is connected to the th capacitor. Otherwise, the vector is zero.(iii)Since two sources cannot be connected to the same -inductor, each column of also has only one entry different from zero.(iv)The sum of the entries different from zero of the th row of equals .(v)The entries different from zero of each column of stand for the connection of sources with other sources.(vi)The sum of the entries different from zero of the th row of equals . Therefore, of (C.3).(vii)Matrix is a zero matrix that exhibits the fact that sources cannot be related to -inductors.

Due to structure (55), the following properties hold.

Property 5. The vector , that is, the vector filled with ones of dimension , is a right eigenvector of matrix satisfying .

Proof. The columns of matrix denote how a source is connected to another source. Since for each pair of sources this connection is carried out by one -inductor, the current leaves one cutset and enters another, which is represented by the fact that, in each column of the matrix, only a in one row and a in another row appear. The proof is ended by direct computation.

Property 6. The following identity holds:where , with as previously defined, and .

Proof. The proof comes as a direct consequence of the previous propositions.

The last part of this section shows that the general network presented above can be specialized to three typical circuits.

(1 ) Radial Network. A radial network is the simplest topology that can be found in the literature and can appear under several scenarios [10, 19, 20]. It is equipped with only one source connected to all the loads, as illustrated in Figure 1.

Figure 1: Example of radial topology.

The main advantage of this class of circuits lies in its simplicity, which in terms of the characterization carried out in this paper is reflected in the fact that it is conformed only by -inductors. Thus, and since there are as many inductors as capacitors, .

Under these conditions, a radial network is topologically characterized by

(2) Ring Network. This configuration looks for improving reliability by connecting each load to a second source, aiming that in case of a failure in one source the loads continue being connected to an energy supplier. One example of this kind of network is shown in Figure 2.

Figure 2: Example of ring topology.

The topological properties of this kind of circuit are as follows:(i)There are no -inductors.(ii)Only two sources are included leading to the fact that and that the number of capacitors that simultaneously satisfy (C.1) and (C.2) is also 2.(iii)Due to the last item, , hence, the first two row vectors of matrix are zero.

Taking into account the listed features, the matrices of a ring network take the formwhere in this case while .

(3) Mesh Network. This is the most complex but at the same time the most reliable network. Its main advantage comes from the fact that every load is connected to all sources of the circuit. Actually, it can be viewed as a radial circuit with links between sources, allowing, in addition, the sources to be connected between them. A typical configuration of this class is presented in Figure 3.

Figure 3: Example of mesh topology.

Regarding the general structure introduced before, only few particularizations can be done:(i)The number of sources equals the number of capacitors; then, .(ii)Since all the capacitors are connected to all sources, for all .(iii)There are no -inductors or (C.2) class capacitors.

As a consequence of the properties listed above, the matrices that topologically characterize a mesh network arewhere is the number of sources and is the number of sources connected to other sources, whilewith the particular feature that all the partitions of the th row of include an element different from zero since all the sources are connected to all capacitors.

Equations (58)–(61) fully characterize the topological structure of the three typical networks. It is the authors’ belief that this unified model can be used to further analyze current topologies or to develop new designs.

6. Steady-State Characterization: Two Case Studies

The aim of this section is to exploit the structure of the fundamental loop matrix under two different scenarios to exhibit the systematization that can be achieved to determine the dynamic behavior of a given network. Although the main objective is to point out the advantages of the approach, it is important to mention that the two case studies are closely related to situations often found in practice.

6.1. Lossless DC Network

The first case study considers a network operating under constant voltage sources, with possibly nonlinear capacitors and inductors, with linear cotree resistances and assuming that the tree resistances are equal to zero; that is, . This scenario is frequently considered in practice, for example, in electrical power systems, when it is considered that the transmission lines are very long [10]. In addition, it is considered that Properties 1 and 2 hold.

Under the aforementioned conditions, model (16) reduces to

Regarding the stability properties of the network, Proposition 9 can be directly applied although in this case, due to the lossless assumption, only stability of the equilibrium point can be concluded, since

However, it is easy to verify, by direct substitution of in (62), that the maximal invariant set of the system corresponds to , a condition that proves that asymptotic stability is attained.

The characterization of the steady-state behavior of the circuit is defined by the equilibria of the system, which are the solutions of

At this point, it is clear that, for a given , the characterization can be carried out for the capacitor voltages or the inductor currents. Motivated by the stability analysis, if the capacitor voltages are chosen, it is possible to write thatsince, from Property 3, is row full rank.

Departing from this last expression, in the following proposition, we illustrated that, exploiting the structures for and , it is straightforward to conclude the steady-state operation achieved by a given network. This result is illustrated for the case of a mesh circuit.

Proposition 16. Consider a mesh electrical circuit characterized by (60) and (61). Assume the following: (i)The network is lossless; that is, .(ii)Propositions 6 and 9 hold.(iii)The vector of voltage sources is composed of constant values.Under these conditions, capacitor voltage achieves average consensus [16] in the sense thatwithbeing the steady-state average value of .

Proof. From (60), it is directly obtained thatleading to , where is the identity matrix.
On the other hand, from (61), all the partitions of the columns of matrix contain an element different from zero; thus, the productleads to the result that , that is, an matrix filled with ones.
Putting all the obtained results together in (65), it holds thatwhich is equivalent to (66).

Remark 17. A direct corollary of the last proposition, which can be alternatively proved from Property 6, refers to the case when , with . Under this condition, voltage capacitor consensus [16] is achieved in the sense that .

Remark 18. It is easy to verify that for radial and ring topologies consensus on the voltages is attained.

Remark 19. If inductive loses are included, , consensus is no longer preserved. Instead, the entries of the matrix depend on the values of the tree resistors. However, it seems that this situation allows for designing compensation techniques, that is, adding new lumped elements to the circuit, such that a prescribed behavior is accomplished.

6.2. AC Steady-State Behavior

The second approached case study consists of an electric network operating under sinusoidal voltage sources. In this case, for facility of analysis, all the passive elements are considered to be linear although, in contrast to the first case study, inductive loses are included. Hence, the circuit dynamic is described by (16) together with (23) and, for convenience, is represented aswhere with , , and as previously defined.

Since the steady-state behavior is now time-varying, the admissible trajectories are given as solution of

As usual [2], to carry out the analysis, it is assumed that voltage and currents are of the form with the phasor . Therefore, admissible trajectories are defined byleading toConsidering that the circuit satisfies Properties 1 and 2 (i.e., identities (40) and (42) hold), the model reduces towhere , is the capacitance matrix, and is the inductance matrix. Hence, it is obtained thatwithand diagonal admittance matrices

In the next proposition, for illustrative purposes, the steady-state behavior of capacitor voltages is characterized assuming that the network is of the mesh type, giving continuity to the first case study. It is convenient to introduce some partitions of several matrices that facilitate the computations involved in the result.

In accordance with the results presented in Section 5, the inductance matrix takes the formwith , . In addition, the former must be divided into matrices aswhere since for a mesh network .

From Properties 1 and 2, the dimension of is , due to the fact that , while ; that is, it equals the dimension of .

Under the aforementioned partitions, matrix introduced in (80) is given bywhere with , , and are submatrices of of dimension corresponding to .

Concerning matrix , whose structure is presented in (61), the submatrix is represented aswith , , matrices composed of the vectors introduced in (56) exhibiting only one entry different from zero on each of their columns.

Proposition 20. Consider a mesh electrical circuit characterized by (60) and (61). Assume the following: (i)Propositions 6 and 9 hold.(ii)The vector of voltage sources is composed of sinusoidal functions.Under these conditions, the input/output relationship between the capacitor voltage phasors and the voltage source phasors is given bywhere , .

Proof. Defining and factorizing , the first equation of (77) can be equivalently written aswhere, using (60), it is straightforward to computewhile, recalling (61), it is easy to show thatwhere has been partitioned in a congruent way.
The proof concludes by substituting (87) and (88) into (86).

The usefulness of the presented result lies in the very well defined structure of (85) which allows for systematic manipulation of the capacitance, inductance, and resistance values in order to force a prescribed steady-state behavior for the capacitor voltages. For example, consider the mesh circuit of Figure 3 which has associated the matrices

In this case, the expressions for the capacitor voltages are given by with , , , and . It is clear that for a given value of the voltages sources a proper choice of the circuit parameters can lead to consensus in the capacitor voltages.

Remark 21. Evidently, the idea to manipulate the capacitance, inductance, and resistance values in order to attain a given steady-state behavior is not new. Actually, in many applications, this procedure is recognized as compensation and is related to the addition of new elements parallel or series connected with the original ones. The advantage offered by the approach presented in this paper is that the compensation analysis can be carried out in a systematic way.

7. Concluding Remarks

In this paper a dynamic characterization of a class of electrical circuits has been presented. The approached circuit topologies are the most used in practice, namely, radial, ring, and mesh typical networks. The characterization contemplates both stability properties and steady behavior and its main feature is that it is based on the structural properties of the networks which have been obtained by using arguments from the graph theory. It has been shown that the structure of these circuits, initially conceived for reliability issues, strongly defines their dynamic behavior. The usefulness of the results reported in this paper lies in the possibility of carrying out the characterization in a very systematic way. In addition, it has been shown that these reported results offer an alternative to deal with problems like compensator location to attain a prescribed behavior.

Competing Interests

The authors declare that they have no competing interests.


Part of this work was supported by DGAPA-UNAM under Grant IN116516.


  1. R. K. Brayton and J. K. Moser, “A theory of nonlinear networks. I,” Quarterly of Applied Mathematics, vol. 22, pp. 1–33, 1964. View at Google Scholar · View at MathSciNet
  2. C. Desoer and E. Kuh, Basic Circuit Theory, McGraw-Hill, New York, NY, USA, 1969.
  3. L. Weiss and W. Mathis, “A hamiltonian formulation for complete nonlinear RLC-networks,” IEEE Transactions on Circuits and Systems. I. Fundamental Theory and Applications, vol. 44, no. 9, pp. 843–846, 1997. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  4. B. M. Maschke, A. J. Van der Schaft, and P. C. Breedveld, “An intrinsic Hamiltonian formulation of the dynamics of LC-circuits,” IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications, vol. 42, no. 2, pp. 73–82, 1995. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  5. D. Jeltsema and J. M. Scherpen, “A dual relation between port-Hamiltonian systems and the Brayton-Moser equations for nonlinear switched RLC circuits,” Automatica, vol. 39, no. 6, pp. 969–979, 2003. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  6. A. J. Van der Schaft and B. M. Maschke, “Port-hamiltonian systems on graphs,” SIAM Journal on Control and Optimization, vol. 51, no. 2, pp. 906–937, 2013. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  7. T. Stykel, “Balancing-related model reduction of circuit equations using topological structure,” in Model Reduction for Circuit Simulation, vol. 74 of Lecture Notes in Electrical Engineering, pp. 53–83, Springer, Dordrecht, Netherlands, 2011. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  8. R. Ortega, D. Jeltsema, and J. M. Scherpen, “Power shaping: a new paradigm for stabilization of nonlinear RLC circuits,” IEEE Transactions on Automatic Control, vol. 48, no. 10, pp. 1762–1767, 2003. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  9. S. Fiaz, D. Zonetti, R. Ortega, J. M. Scherpen, and A. J. van der Schaft, “A port-Hamiltonian approach to power network modeling and analysis,” European Journal of Control, vol. 19, no. 6, pp. 477–485, 2013. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  10. P. Kundur, N. Balu, and M. Lauby, Power System Stability and Control, vol. 7, McGraw-Hill, New York, NY, USA, 1994.
  11. A. Van der Schaft, “Characterization and partial synthesis of the behavior of resistive circuits at their terminals,” Systems & Control Letters, vol. 59, no. 7, pp. 423–428, 2010. View at Publisher · View at Google Scholar · View at Scopus
  12. R. Ortega, A. van der Schaft, F. Castanos, and A. Astolfi, “Control by interconnection and standard passivity-based control of Port-Hamiltonian systems,” IEEE Transactions on Automatic Control, vol. 53, no. 11, pp. 2527–2542, 2008. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  13. J. Bondy and S. Murty, Graph Theory with Applications, North-Holland Publishing, 1976.
  14. B. Bollobás, Modern Graph Theory, vol. 184 of Graduate Texts in Mathematics, Springer Science & Business Media, 1998. View at Publisher · View at Google Scholar · View at MathSciNet
  15. P. Wellstead, Introduction to Physical System Modelling, Academic Press, London, UK, 1979. View at MathSciNet
  16. H. Bai, M. Arcak, and J. Wen, Cooperative Control Design: A Systematic, Passivity-Based Approach, Springer, 2011.
  17. A. Van der Schaft, L2-Gain and Passivity in Nonlinear Control, Springer, New York, NY, USA, 1999.
  18. B. Jayawardhana, R. Ortega, E. García-Canseco, and F. Castaños, “Passivity of nonlinear incremental systems: application to PI stabilization of nonlinear RLC circuits,” Systems & Control Letters, vol. 56, no. 9-10, pp. 618–622, 2007. View at Publisher · View at Google Scholar · View at Scopus
  19. Westinghouse Electric Corporation, Electrical Transmission and Distribution Reference Book, Westinghouse Electric Corporation, East Pittsburgh, Pa, USA, 1965.
  20. U. Eminoglu and M. H. Hocaoglu, “A network topology-based voltage stability index for radial distribution networks,” International Journal of Power and Energy Systems, vol. 29, no. 2, pp. 131–143, 2009. View at Publisher · View at Google Scholar · View at Scopus