Os fundamentos teóricos desenvolvidos neste trabalho e que dão suporte aos métodos propostos garantem que as estimativas obtidas sejam sempre conservadoras (no sentido de que elas são sempre subconjuntos da região de estabilidade verdadeira) e, portanto, possuam elevado grau de confiança ao concluir sobre a estabilidade do sistema. Os métodos apresentados consistem em extensões dos métodos Closest UEP e CUEP, utilizados na análise de estabilidade transitória de sistemas elétricos de potência, para sistemas que admitem FEG. Embora os métodos Closest UEP e CUEP forneçam estimativas de forma rápida e precisa, sua aplicação está limitada à existência de uma Função Energia (FE) para o sistema, o que consiste em uma forte limitação. Muitos sistemas não admitem FE e, mesmo quando se pode provar a existência de uma FE, a impossibilidade de exibi-la impede a aplicação dos métodos citados. Outra contribuição deste trabalho consiste em um método computacional que permite a obtenção de uma FEG para sistemas polinomiais. O método apresentado também é aplicado a uma classe de problemas não polinomiais, provenientes da modelagem de sistemas elétricos de potência, mediante uma mudança não linear de variáveis que permite a construção de um sistema polinomial equivalente. Através dos métodos apresentados, visa-se disponibilizar métodos computacionais que permitam a obtenção de estimativas rápidas e precisas e que possam ser aplicados a uma ampla classe de sistemas: aqueles que admitem FEG. Com isso, almeja-se não somente contribuir para o desenvolvimento de métodos para análise de estabilidade de sistemas elétricos de potência mas, também, disponibilizá-los a outras áreas do conhecimento. / In this work, we develop computational methods to estimate stability regions and the relevant part of stability boundary of attracting sets of nonlinear dynamical systems. Such methods are based on Generalized Energy Function (GEF) theory and, therefore, can be applied to a larger class of problems than those based on Energy Functions (EF). The theoretical foundations developed in this work, which support the proposed methods, ensure that the estimates are always conservative (in the sense that they are subsets of the true stability region), providing high confidence level when asserting the stability of a system. The presented methods are extensions of the Closest UEP and the CUEP methods, used in the assessment of stability of electrical power systems, to the systems that admit GEF. Even though the Closest UEP and CUEP methods provide estimates in a fast and accurate way, they are only applicable to systems that admit EFs, which consists in a strong limitation for their usage. Many systems do not admit EF and, even if it is possible to prove the existence of an EF, the impossibility to exhibit it in the form of elementary mathematical functions prevents the application of such methods. Other contribution of this work is a computational method to obtain a GEF for polinomial systems. We also applied the presented method to a class of non polinomial systems arising from electrical power system models, after a nonlinear change of variables that provides an equivalent polinomial system. By means of the proposed methods, we aim to offer computational methods to allow fast and accurate stability region estimates which could be used in a broad class of dynamical systems: those that admit GEF. This way, we plan to contribute for the development of methods used in the assessment of stability of electrical power systems and make such tools available to systems from other areas of science.

Os fundamentos teóricos desenvolvidos neste trabalho e que dão suporte aos métodos propostos garantem que as estimativas obtidas sejam sempre conservadoras (no sentido de que elas são sempre subconjuntos da região de estabilidade verdadeira) e, portanto, possuam elevado grau de confiança ao concluir sobre a estabilidade do sistema. Os métodos apresentados consistem em extensões dos métodos Closest UEP e CUEP, utilizados na análise de estabilidade transitória de sistemas elétricos de potência, para sistemas que admitem FEG. Embora os métodos Closest UEP e CUEP forneçam estimativas de forma rápida e precisa, sua aplicação está limitada à existência de uma Função Energia (FE) para o sistema, o que consiste em uma forte limitação. Muitos sistemas não admitem FE e, mesmo quando se pode provar a existência de uma FE, a impossibilidade de exibi-la impede a aplicação dos métodos citados. Outra contribuição deste trabalho consiste em um método computacional que permite a obtenção de uma FEG para sistemas polinomiais. O método apresentado também é aplicado a uma classe de problemas não polinomiais, provenientes da modelagem de sistemas elétricos de potência, mediante uma mudança não linear de variáveis que permite a construção de um sistema polinomial equivalente. Através dos métodos apresentados, visa-se disponibilizar métodos computacionais que permitam a obtenção de estimativas rápidas e precisas e que possam ser aplicados a uma ampla classe de sistemas: aqueles que admitem FEG. Com isso, almeja-se não somente contribuir para o desenvolvimento de métodos para análise de estabilidade de sistemas elétricos de potência mas, também, disponibilizá-los a outras áreas do conhecimento. / In this work, we develop computational methods to estimate stability regions and the relevant part of stability boundary of attracting sets of nonlinear dynamical systems. Such methods are based on Generalized Energy Function (GEF) theory and, therefore, can be applied to a larger class of problems than those based on Energy Functions (EF). The theoretical foundations developed in this work, which support the proposed methods, ensure that the estimates are always conservative (in the sense that they are subsets of the true stability region), providing high confidence level when asserting the stability of a system. The presented methods are extensions of the Closest UEP and the CUEP methods, used in the assessment of stability of electrical power systems, to the systems that admit GEF. Even though the Closest UEP and CUEP methods provide estimates in a fast and accurate way, they are only applicable to systems that admit EFs, which consists in a strong limitation for their usage. Many systems do not admit EF and, even if it is possible to prove the existence of an EF, the impossibility to exhibit it in the form of elementary mathematical functions prevents the application of such methods. Other contribution of this work is a computational method to obtain a GEF for polinomial systems. We also applied the presented method to a class of non polinomial systems arising from electrical power system models, after a nonlinear change of variables that provides an equivalent polinomial system. By means of the proposed methods, we aim to offer computational methods to allow fast and accurate stability region estimates which could be used in a broad class of dynamical systems: those that admit GEF. This way, we plan to contribute for the development of methods used in the assessment of stability of electrical power systems and make such tools available to systems from other areas of science.

The peridynamic theory is a generalization of classical continuum mechanics and takes into account the interaction between material points separated by a finite distance within a peridynamic horizon δ. The parameter δ corresponds to a length scale and is treated as a material property related to the microstructure of the body. Since the balance of linear momentum is written in terms of an integral equation that remains valid in the presence of discontinuities, the peridynamic theory is suitable for studying the material behavior in regions with singularities. The first part of this work concerns the evaluation of the properties of a linear elastic peridynamic material in the context of a three-dimensional state-based peridynamic theory, which uses the difference displacement quotient field in the neighborhood of a material point and considers both length and relative angle changes. This material model is based upon a free energy function that contains four material constants, being, therefore, different from other peridynamic models found in the literature, which contain only two material constants. Using convergence results of the peridynamic theory to the classical linear elasticity theory in the limit of small horizons and a correspondence argument between the free energy function and the strain energy density function from the classical theory, expressions were obtained previously relating three peridynamic constants to the classical elastic constants of an isotropic linear elastic material. To calculate the fourth peridynamic material constant, which couples both bond length and relative angle changes, the correspondence argument is used once again together with the strain field of a linearly elastic beam subjected to pure bending. The expression for the fourth constant is obtained in terms of the Poisson\'s ratio and the shear elastic modulus of the classical theory. The validity of this expression is confirmed through the consideration of other experiments in mechanics, such as bending of a beam by terminal loads and anti-plane shear of a circular cylinder. In particular, numerical results indicate that the expressions for the constants are independent of the experiment chosen. The second part of this work concerns an investigation of the behavior of a one-dimensional linearly elastic bar of length L in the context of the peridynamic theory; especially, near the ends of the bar, where it is expected that the behavior of the peridynamic bar may be very different from the behavior of a classical linear elastic bar. The bar is in equilibrium without body force, is fixed at one end, and is subjected to an imposed displacement at the other end. The bar has micromodulus C, which is related to the Young\'s modulus E in the classical theory through different expressions found in the literature. Depending on the expression for C, the displacement field may be singular near the ends, which is in contrast to the linear behavior of the displacement field observed in classical linear elasticity. In spite of the above, it is also shown that the peridynamic displacement field converges to its classical counterpart as the peridynamic horizon tends to zero. / A teoria peridinâmica é uma generalização da teoria clássica da mecânica do contínuo e considera a interação de pontos materiais devido a forças que agem a uma distância finita entre si, além da qual considera-se nula a força de interação. Por ter o balanço de momento linear formulado como uma equação integral que permanece válida na presença de descontinuidades, a teoria peridinâmica é adequada para o estudo do comportamento de materiais em regiões com singularidades. A primeira parte deste trabalho consiste no cálculo das propriedades de um material peridinâmico elástico linear no contexto de uma teoria peridinâmica de estado, linearmente elástica e tridimensional, que utiliza o campo quociente de deslocamento relativo na vizinhança de um ponto material e leva em conta mudanças relativas angulares e de comprimento. Esse modelo utiliza uma função energia livre que apresenta quatro constantes materiais, sendo, portanto, diferente de outros modelos peridinâmicos investigados na literatura, os quais contêm somente duas constantes materiais. Utilizando resultados de convergência da teoria peridinâmica para a teoria de elasticidade linear clássica no limite de pequenos horizontes e um argumento de correspondência entre as funções energia livre proposta e densidade de energia de deformação da teoria clássica, expressões para três constantes peridinâmicas foram obtidas em função das constantes de um material elástico e isotrópico da teoria clássica. O argumento de correspondêmcia, em conjunto com o campo de deformações de uma viga submetida à flexão pura, é utilizado para calcular a quarta constante peridinâmica do material, que relaciona mudanças angulares relativas e de comprimentos das ligações entre as partículas. Obtem-se uma expressão para a quarta constante em termos do coeficiente de Poisson e do módulo de elasticidade ao cisalhamento da teoria clássica. A validade dessa expressão é confirmada por meio da consideração de outros experimentos da mecânica, tais como flexão de um viga por cargas terminais e cisalhamento anti-plano de um eixo cilíndrico. Em particular, os resultados numéricos indicam que as expressões para as constantes são independentes do experimento escolhido. A segunda parte deste trabalho consiste em uma investigação do comportamento de uma barra unidimensional linearmente elástica de comprimento L no contexto da teoria peridinâmica; especialmente, próximo às extremidades da barra, onde espera-se que o comportamento da barra peridinâmica possa ser muito diferente do comportamento de uma barra elástica linear clássica. A barra está em equilíbrio e sem força de corpo, fixa em uma extremidade, e sujeita a deslocamento imposto na outra extremidade. A barra possui micromódulo C, o qual está relacionado ao módulo de Young E da teoria clássica por meio de diferentes expressões encontradas na literatura. Dependendo da expressão para C, o campo de deslocamento pode ser singular próximo às extremidades, o que contrasta com o comportamento linear do campo de deslocamento observado na elasticidade linear clássica. Apesar disso, é mostrado também que o campo de deslocamento peridinâmico converge para o campo de deslocamento da teoria clássica quando o horizonte peridinâmico tende a zero.

Este trabalho está dividido em duas partes distintas que constituem contribuições inéditas ao estudo da estabilidade em sistemas elétricos de potência. A primeira parte do trabalho é a mais importante e trata do problema da identificação das soluções de baixa tensão críticas do fluxo de carga. Esta parte do trabalho se presta a análise de estabilidade de tensão a pequenas perturbações. Os últimos capítulos deste trabalho apresentam também uma proposta de função energia estendida que modela as cargas dependentes da tensão segundo o modelo ZIP de carga, considerando a estrutura da rede preservada. Assim, a função energia proposta pode ser utilizada para analisar tanto a estabilidade de tensão como a estabilidade de ângulo em sistemas de potência. Esta proposta também é inédita na literatura. Embora a função energia proposta tenha sido aplicada apenas a sistemas de dimensão reduzidas, os resultados apresentados neste trabalho nos levam a acreditar que essa mesma função energia pode ser utilizada na análise de estabilidade de sistemas de potência de grandes dimensões. Já o método proposto para identificação das soluções de baixa tensão das equações de fluxo de carga se utiliza de um sistema dinâmico auxiliar das equações de fluxo de carga. O sistema dinâmico auxiliar utilizado não tem significado físico, mas pode ser escolhido de tal forma que a solução usual das equações de fluxo de carga seja um ponto de equilíbrio estável do sistema dinâmico auxiliar, eque as soluções de baixa tensão do fluxo de carga sejam pontos de equilíbrio instáveis do sistema dinâmico auxiliar. Dessa forma, é possível calcular as soluções de baixa tensão do fluxo de carga, calculando-se os pontos de equilíbrio instáveis do sistema dinâmico auxiliar. Assim, é possível utilizar partes da teoria de sistemas dinâmicos para estudar as soluções das equações de fluxo de carga. Baseado nestes princípios, foi desenvolvido um programa para calcular trajetórias do sistema dinâmico auxiliar, que se iniciam e se mantêm nas vizinhanças da fronteira da área de atração do ponto de equilíbrio estável do SEP. Dessa forma é possível afirmar que a trajetória calculada tende a convergir para a solução crítica das equações de fluxo de carga. O programa foi inicialmente concebido para calcular as soluções de baixa tensão de sistemas elétricos sem perdas. Em seguida o programa desenvolvido foi adaptado para calcular as soluções de baixa tensão de sistemas de potência completos, incluindo também as resistências das linhas de transmissão. Esta última versão do programa foi testada para os sistemas IEEE 39 e IEEE 118 barras, e os resultados obtidos se mostraram bastante satisfatórios. Assim, o método proposto é uma ferramenta original e eficaz para a solução do problema de calcular a solução crítica das equações de fluxo de carga de sistemas elétricos de potência. / This work may be divided into two distinct parts. Both of them are new contributions to stability analysis of power systems. In the first part it is proposed a new method to calculate the critical load flow low voltage solutions, and it is the main part of this work. Meanwhile, the last two chapters of this work presents a proposed extended energy function that consider the common load ZIP models. It allows the analysis of angle and voltage stability for power systems subjected to large disturbances. This work proposes a method to calculate the low voltage solutions (LVS) of the load flow equations of an electrical power system. The proposed method identifies the LVS involved in the saddle-node bifurcation leading the power system to a voltage collapse. This solution is known as the critical low voltage solution. In order to perform the proposed calculation, an auxiliary dynamical gradient system is used. It is shown that the equilibrium points of that associated auxiliary dynamical gradient system are the solutions of the load flow equations. In such manner, the paper proposes identifying the critical LVS calculating the equilibrium points of an auxiliary dynamical gradient system. The proposed method was tested on the Stagg 5-bus, on the IEEE 39-bus and on IEEE 118-bus test systems, and the results are presented at the end of the text.