Orientador: Laurence Duarte Colvara / Banca: Carlos Roberto Minussi / Banca: Wagner Peron Ferreira / Resumo: Este trabalho apresenta um estudo sobre a relação entre os principais tipos de bifurcações que ocorrem em sistemas elétricos de potência e em quais ocasiões elas podem ocorrer em máquinas síncronas com ou sem RAT (Regulador Automático de Tensão). Para explorar tais fenômenos, primeiramente o sistema é modelado, sendo utilizado para o caso MBI (Máquina - Barramento In nito) o modelo um eixo e, então, a matriz de estado é calculada para a análise dos autovalores. Para os sistemas multimáquinas estudados, são incluídos dois enrolamentos amortecedores nos eixos d ¡ q. São então apresentados os métodos de análise de estabilidade transitória convencionais, amplamente utilizados, conhecidos como método Tradicional e Método Direto. As condições para a ocorrência de bifurcações são analisadas utilizando os coe - cientes linearizados do modelo de He ron-Phillips para o caso MBI, onde é mostrado que se espera perder a estabilidade para o sistema com regulador automático de tensão através de uma bifurcação de Hopf e para o caso sem RAT através de uma bifurcação Sela-Nó. Por m, é analisado o ciclo-limite para o caso de uma máquina - barramento in nito e para sistemas multimáquinas através do modelo não-linear. A região de estabilidade é analisada no plano de fase, sendo mostrada a necessidade de incluir a variação de uxo no enrolamento de campo para uma análise correta da estabilidade. É também mostrado que este ciclo-limite pode reduzir a fronteira de estabilidade calculada pelo método convencional. / Abstract: The aim of this study is the relation among main types of bifurcations that occur in electrical power systems and the circumstances they can happen with the synchronous machines considered with or without AVR (Automatic Voltage Regulator). To explore such phenomena, the system is rst modeled with the synchronous machines described by the one axis model for the MIB (Machine - In nite Bus) case , and so the state matrix is computed for the analysis of its eigenvalues. For multimachine systems case two windings dampers are included in d-q axes. The conditions for the occurrence of bifurcations are analyzed using the coe cients of the He ron-Phillips model for MIB case, where it is shown that one expects the system with automatic voltage regulator lose synchronism through a Hopf bifurcation and for the case without RAT through a Saddle- Node Bifurcation. Finally, the nonlinear model is accounted for in order to consider the limit-cycle for the case of one machine - in nite bus case as well as for multimachine system. Since internal voltage a ects the boundary of the stability region it must be considered. Then the phase portrait does not su ce and the trajectories must to be observed in a sub space de ned with the internal voltage. It is also shown that this limit-cycle can reduce the boundary of stability calculated by means of the direct method. / Mestre

Orientador: Edevar Luvizotto Junior / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Civil, Arquitetura e Urbanismo / Made available in DSpace on 2018-08-11T21:41:10Z (GMT). No. of bitstreams: 1 Anjo_LuizFernandoResendedosSantos_D.pdf: 3588080 bytes, checksum: 122271523488dc8a40a0d509bd1b0cf8 (MD5) Previous issue date: 2008 / Resumo: Este texto descreve as etapas que objetivam a utilização da estrutura originalmenteproposta por Todini e Pilati no chamado método gradiente (MG), utilizado para análise em regime permanente em instalações a condutos forçados, na formulação de um modelo dinâmico inercial rígido (MDIR), para a análise de escoamentos transitórios lentos neste tipo de instalação. São apresentadas as bases teóricas para esta nova modelação, justificadas através do equacionamento geral do escoamento fluido em condutos forçados. Os resultados obtidos pelo MDIR são comparados com os resultados obtidos pelo programa EPANET que utiliza o método gradiente. Discussões a respeito da importância da incorporação do efeito de inércia são apresentadas através de um estudo de casos, no caso de modelação hidráulica, e para análises de qualidade decorrentes desta. / Abstract: This research describes the stages which aim the use of the original structure proposed by Todini and Pilati in the so called gradient method (GM), used in the analysis in steady state of pipe networks, in the formulation of an inertial rigid dynamic model (IRDM) to analyse slow transients in this type of installation. The theoretical bases are presented for this new method, justified by the general equation of fluids flow in pipes networks. The results obtained by the MDIR are compared to the results obtained in the EPANET program, which uses the gradient method. Discussion on the importance of the incorporation of the inertia effect are presented through a study case of hydraulic modelation and for quality analysis which derived from it. / Doutorado / Recursos Hidricos / Doutor em Engenharia Civil

[pt] Este trabalho discute técnicas para simulação de redes elétricas com ênfase em operação em tempo real. As redes possuem linhas de transmissão como elementos a parâmetros distribuídos, enquanto redes RLC correspondem aos elementos concentrados, se bem que muitos dos resultados não se restringem só a estes tipos de elementos. A modelagem de linhas de transmissão é considerada primeiro, sob o ponto de vista da comparação entre dois modelos da literatura, caracterizando sua adequação para simulação em tempo real. Desta comparação resulta uma combinação de procedimentos como adequada à simulação. Na seqüência, é apresentado um novo algoritmo de integração numérica, baseado na descrição de estado de rede e em filtragem digital multitaxa. O novo algoritmo mostrou que pode reduzir o esforço computacional da simulação da rede, ao mesmo tempo que melhorou a qualidade das mesmas. O próximo objetivo relaciona-se com as oscilações numéricas, causadas pelo uso da Integração Trapezoidal em algumas configurações de rede. Um método de análise, baseado em teoria de rede e processamento digital de sinais é desenvolvido. A Partir dele são identificadas as configurações de redes, passíveis de gerações numéricas na sua contraparte digital. Métodos convencionais de eliminação das oscilações numéricas são analisados, incluindo no método o sistema multitaxa. Em todos os casos o sistema novo apresentou melhores resultados, sendo que seu uso apropriado, pode, por se, eliminar as oscilações numéricas. Alguns casos ilustrativos de redes de sistemas de potência estão incluídos, sendo desenvolvidos em linguagem MATLAB. Eles indicam que a combinação do modelo sugerido para linha de transmissão com o novo método de integração numérica multitaxa, permite a simulação da rede com reduzido esforço computacional e portanto aumentar a possibilidade de simulação de redes complexas em tempo real. / [en] This work discusses techniques for electrical network digital simulation, emphasizing real time operation. The networks contain transmission lines as distributed elements and RLC elements as lumped parameters, although many of the results are not restricted to these elements. Transmission line modeling is considered first, comparing two of the methods found in literature from the real time simulation point of view. The convenience of their use is subject of analysis, leading, later on, to suggest a blend of approaches for real time simulation. In the sequence of the work, it is introduced a new numerical integration algorithm, based on network state space description and multirate digital filtering. The new integration method reduces the simulation computational complexity while simultaneously improving its accuracy. Next, the focus is on numerical oscillations caused by the Trapezoidal Integration approach in some network configurations. It is presented a method to identify electrical network configurations that might generate numerical oscillations in their digital counterparts. The method is a result of both, network and digital signal theory. Conventional and new techniques are applied in order to eliminate the numerical oscillations. The new multirate approach improves convetional techniques results, while properly used eliminates by itself the numerical oscillations. Some power systems network illustrative cases are included, using the new integration method and the proposed digital transmission line model. It is shown that this combination produces improved performance, final accuracy and reduced computational effort, resulting appropriate for real time simulations.

[pt] O trabalho desenvolvido estabelece uma metodologia de estudo para a determinação dos requisitos de tensão de restabelecimento transitória de disjuntores através de simulações em computador digital com a utilização do programa de cálculo de transitórios elétricos Electromagnetic Transients Program (EMTP) de Bonneville Power Administration (BPA). / [en] It is established in this paper a study method to determine the transient recovery voltage requirements for circuit breakers, by means of simulations in digital computer, using the Electromagnetic Transients Program (EMTP) de Bonneville Power Administration (BPA).

[ES] La simulación de transitorios forma parte del proceso de licenciamiento de una central nuclear. Esto implica que los códigos, así como los modelos utilizados deben estar verificados y validados. Normalmente, esta simulación se realiza con códigos termohidráulicos de planta que tienen una definición de la cinética del reactor muy simplificada con cinética puntual o unidimensional. Una mejora importante en la simulación de transitorios base de diseño se basa en la utilización de códigos acoplados termohidráulico-neutrónicos, que permiten obtener resultados sobre la evolución de la potencia del reactor en tres dimensiones. Los códigos neutrónicos 3D necesitan parámetros de la cinética y secciones eficaces también en 3D ajustados al punto del ciclo que se quiere simular y que abarquen las condiciones que se alcancen durante el transitorio. Por otro lado, para poder verificar tanto los códigos como los modelos es necesario llevar a cabo una serie de simulaciones de diferentes transitorios. De esta manera, se comprueba cómo funciona el código acoplado en diferentes condiciones de operación y simulación. Esta tesis contribuye al conocimiento del uso de códigos termohidráulico-neutrónicos acoplados en la simulación de transitorios base de diseño (Design Basis Accidents -DBAs). Los códigos mejorados y verificados son los códigos termohidráulicos RELAP5, TRAC-BF1 y TRACE y el código neutrónico PARCS. Los parámetros neutrónicos necesarios en PARCS se han obtenido aplicando una metodología que simplifica el modelo del núcleo. Esta metodología, ya desarrollada e implementada, denominada SIMTAB, se ha mejorado, tanto en las posibilidades de aplicación de la misma como en la optimización y actualización de la programación del código fuente. Los transitorios analizados con los códigos RELAP5/PARCS acoplados son: transitorio por expulsión de barra de control y transitorio de inyección de boro en un reactor PWR. Con los códigos TRAC-BF1/PARCS acoplados se ha analizado el transitorio por disparo de turbina en la C. N. Peach Bottom. Para llevar a cabo las simulaciones con TRAC-BF1/PARCS se ha implementado el acoplamiento de ambos códigos, puesto que originalmente el código TRAC-BF1 no estaba preparado para ello. El análisis de inestabilidades en reactores BWR se ha realizado con RELAP5/PARCS en dos reactores BWR: C. N. Peach Bottom y C. N. Ringhals 1. Para ello se ha desarrollado una metodología de análisis que abarca desde la definición del modelo termohidráulico y del modelo neutrónico hasta el análisis de las señales simuladas obtenidas con PARCS. La metodología también incluye la aplicación de diferentes perturbaciones basadas en los modos Lambda y en el análisis de las señales reales de planta. Se ha llevado a cabo un estudio del modelo para el cálculo de la concentración de Boro en los códigos termohidráulicos y se ha mejorado este modelo en el código TRAC-BF1, incorporando un nuevo método de resolución en el código fuente. El modelo para el cálculo del calor de desintegración también se ha revisado y mejorado en los códigos TRAC-BF1 y PARCS. En ambos casos se ha implementado el modelo ANS 2005. El análisis de sensibilidad e incertidumbre está ligado a los resultados de los códigos de mejor estimación como los mejorados en esta tesis. Este análisis se ha realizado sobre los transitorios de expulsión de barra en un reactor PWR y el transitorio de caída de barra en un reactor BWR con RELAP5/PARCS. Los resultados de estos trabajos aportan una metodología de aplicación para la simulación correcta de transitorios con códigos acoplados. Además, ha servido para detectar y subsanar deficiencias en los códigos, y de esta manera disponer de unos códigos de mejor estimación preparados para el análisis de transitorios base de diseño. / [CA] La simulació de transitoris forma part del procés de llicenciament d'una central nuclear. Això implica que els codis, així com els models utilitzats han d'estar verificats i validats. Normalment, aquesta simulació es realitza amb codis termohidràulics de planta que tenen una definició de la cinètica del reactor molt simplificada amb cinètica puntual o unidimensional. Una millora important en la simulació de transitoris base de disseny es basa en la utilització de codis acoblats termohidràulic-neutrònics, que permeten obtindre resultats sobre l'evolució de la potència del reactor en tres dimensions. Els codis neutrònics 3D necessiten paràmetres de la cinètica i seccions eficaces també en 3D ajustats al punt del cicle que es vol simular i que abasten les condicions que s'aconseguisquen durant el transitori. D'altra banda, per a poder verificar tant els codis com els models és necessari dur a terme una sèrie de simulacions de diferents transitoris. D'aquesta manera, es comprova com funciona el codi acoblat en diferents condicions d'operació i simulació. Aquesta tesi contribueix al coneixement de l'ús de codis termohidràulic-neutrònics acoblats en la simulació de transitoris base de disseny. Els codis millorats i verificats són els codis termohidràulics RELAP5, TRAC-BF1 i TRACE i el codi neutrònic PARCS. Els paràmetres neutrònics necessaris en PARCS s'han obtingut aplicant una metodologia que simplifica el model del nucli. Aquesta metodologia, ja desenvolupada i implementada, denominada SIMTAB, s'ha millorat, tant en les possibilitats d'aplicació de la mateixa com en l'optimització i actualització de la programació del codi font. Els transitoris analitzats amb els codis RELAP5/PARCS acoblats són: transitori per expulsió de barra de control i transitori d'injecció de bor en un reactor PWR. Amb els codis TRAC-BF1/PARCS acoblats s'ha analitzat el transitori per disparament de turbina en la C. N. Peach Bottom. Per a dur a terme les simulacions amb TRAC-BF1/PARCS s'ha implementat l'acoblament de tots dos codis, ja que originalment el codi TRAC-BF1 no estava preparat per a això. L'anàlisi d'inestabilitats en reactors BWR s'ha realitzat amb RELAP5/PARCS en dos reactors BWR: C. N. Peach Bottom i C. N. Ringhals 1. Per a això s'ha desenvolupat una metodologia d'anàlisi que abasta des de la definició del model termohidràulic i del model neutrònic fins a l'anàlisi dels senyals simulats. La metodologia també inclou l'aplicació de diferents pertorbacions basades en els modes Lambda i en l'anàlisi dels senyals reals de planta. S'ha dut a terme un estudi del model per al càlcul de la concentració de Bor en els codis termohidràulics i s'ha millorat aquest model en el codi TRAC-BF1, incorporant un nou mètode de resolució en el codi font. El model per al càlcul de la calor de desintegració també s'ha revisat i millorat en els codis TRAC-BF1 i PARCS. En tots dos casos s'ha implementat el model ANS 2005. L'anàlisi de sensibilitat i incertesa està lligat als resultats dels codis de millor estimació com els millorats en aquesta tesi. Aquesta anàlisi s'ha realitzat sobre els transitoris d'expulsió de barra en un reactor PWR i el transitori de caiguda de barra en un reactor BWR amb RELAP5/PARCS. Els resultats d'aquests treballs aporten una metodologia d'aplicació per a la simulació correcta de transitoris amb codis acoblats. A més, ha servit per a detectar i esmenar deficiències en els codis, i d'aquesta manera disposar d'uns codis de millor estimació preparats per a l'anàlisi de transitoris base de disseny. / [EN] The simulation of transients is part of the licensing process of a nuclear power plant. This implies that the codes as well as the models used must be verified and validated. Normally, this simulation is performed with thermalhydraulic plant codes that have a very simplified definition of reactor kinetics with point or one-dimensional kinetics. An important improvement in the simulation of design-basis transients rely on the use of thermohydraulic-neutronic coupled codes, which allow to obtain results of the evolution of the reactor power in three dimensions. The 3D neutron codes need parameters of the kinetics and cross-sections also in 3D adjusted to the point of the cycle to be simulated that must cover the conditions reached during the transient. On the other hand, to be able to verify both the codes and the models it is necessary to carry out a series of simulations of different transients. In this way, it is checked how the coupled code works in different operating and simulation conditions. This thesis contributes to increase the knowledge of the use of thermalhydraulic-neutronic coupled codes in the simulation of design basis accidents (DBAs). The improved and verified codes are the thermalhydraulic codes RELAP5, TRAC-BF1 and TRACE and the neutronic code PARCS. The necessary neutronic parameters in PARCS have been obtained by applying a methodology that simplifies the core model. This methodology, already developed and implemented, called SIMTAB, has been improved in this thesis in its application possibilities and also in the optimization and updating of the source code. The transients analyzed with RELAP5/PARCS coupled code are: control rod ejection transient and boron injection transient in a PWR reactor. With TRAC-BF1/PARCS coupled code, the transient analyzed is the turbine trip transient in Peach Bottom NPP. To carry out the simulations with TRAC-BF1/PARCS, the coupling of both codes has been implemented before, since originally the TRAC-BF1 code was not prepared for it. The analysis of instabilities in BWR reactors has been carried out with RELAP5/PARCS in two BWR reactors: Peach Bottom NPP and Ringhals 1 NPP. A methodology has been developed which cover from the definition of the thermalhydraulic model and the neutron model to the simulated signal analysis. The methodology also includes the application of different disturbances based on Lambda modes and the analysis of real plant signals. A study of the model for the calculation of the Boron concentration in thermalhydraulic codes has been carried out. This model has been improved in the TRAC-BF1 code, incorporating a new resolution method in the source code. The model for the calculation of the decay heat has also been revised and improved in TRAC-BF1 and PARCS codes. In both cases, the ANS 2005 model has been implemented. The sensitivity and uncertainty analysis is linked to the results of the best estimate codes such as those improved in this thesis. This analysis has been carried out on the control rod ejection transients in a PWR reactor and the control rod drop transient in a BWR reactor with RELAP5/PARCS. The results of these works provide an application methodology for the correct simulation of transients with coupled codes. In addition, it has been used to detect and correct deficiencies in the codes, and therefore, to have better estimate codes prepared for the analysis of design-basis transients. / Barrachina Celda, TM. (2020). Aportaciones y Mejoras en los Códigos Termohidráulicos y Neutrónicos de Estimación Óptima RELAP5, TRAC-BF1, TRACE Y PARCS [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/158745 / TESIS

Orientador: Edevar Luvizotto Junior / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Civil, Arquitetura e Urbanismo / Made available in DSpace on 2018-08-07T09:53:33Z (GMT). No. of bitstreams: 1 ArandaSaldana_EdwinAntonio_M.pdf: 1687985 bytes, checksum: 0c5e968a32c41fdf3cc83eca866e6ea6 (MD5) Previous issue date: 2006 / Resumo: O presente trabalho apresenta as bases para a elaboração de um modelo hidráulico computacional baseado na filosofia de programação orientada por objetos (POO) para o tratamento de problemas de escoamento de fluidos em condutos forçados em regime permanente e transitório. Na modelação hidráulica foi empregado um modelo dinâmico inercial elástico com solução através do Método das Características (MOC). Os estudos de casos apresentados mostram o potencial da ferramenta proposta para a análise de escoamentos de instalações hidráulicas em condutos forçados / Abstract: The present work shows the basis for the elaboration of a computational hydraulic model based on the object oriented programming (POO) philosophy for the management of flow problems such as flow rate in pipe systems in steady and unsteady states. The hydraulic modeling is based on the inertial dynamic model with solved through the Method of Characteristics (MOC). The different cases of studies presented in this dissertation show the potential of this tool proposal for the analysis of pipe flow installations in hydraulic systems / Mestrado / Recursos Hidricos / Mestre em Engenharia Civil

Orientador: Paulo Vatavuk / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Civil, Arquitetura e Urbanismo / Made available in DSpace on 2018-08-11T08:47:49Z (GMT). No. of bitstreams: 1 Souza_EdsonVictorde_M.pdf: 1243495 bytes, checksum: 65a1b651d2f304cfca69e91ab51e9514 (MD5) Previous issue date: 2008 / Resumo: O estudo do transitório hidráulico em condutos forçados tem grande relevância em projetos de sistemas de tubulações. A quantificação das pressões máximas e mínimas é de fundamental interesse para o projetista, a fim de que este possa dimensionar a tubulação e/ou prever dispositivos de proteção, cuja finalidade é amortecer as variações de carga, prejudiciais à vida útil da instalação. Neste trabalho optou-se pelo estudo de um destes dispositivos, o reservatório hidropneumático, que protege as instalações contra sobrepressões e subpressões, impedindo assim a descontinuidade do escoamento quando o sistema de bombeamento é desligado ou quando ocorre um fechamento rápido de válvula de bloqueio. O objetivo principal da presente pesquisa é o de formular tabelas e gráficos de aplicação para o pré-dimensionamento de reservatórios hidropneumáticos para a atenuação de transitórios hidráulicos em sistemas de condutos forçados, facilitando, desta maneira, o trabalho dos projetistas na escolha deste tipo de dispositivo de proteção de sistemas hidráulicos / Abstract: The evaluation of hydraulic transients is very important in the design of pipe networks. The determination of the maximum and minimum pressures is of great interest to the designer, so that he can size the pipes and/or design protection devices, with the objective of decrease head variations, that can decrease the instalation life. In this work, one of those devices was chosen for study, the air vessel, that protects the instalation against overpressures and underpressures, avoiding the interruption of the flow when the pumping system is turned off or when a valve closes too fast. The present work has the main objective of making charts and graphs of application for pre-sizing of air vessels to the mitigation of Waterhammer in pressure pipes by facilitating, in this way, the work of designers in choosing this type of protection of hydraulic systems / Mestrado / Recursos Hidricos / Mestre em Engenharia Civil

Orientadores: Sandra Lucia da Cruz, João A. F. R. Pereira / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-07-31T22:23:21Z (GMT). No. of bitstreams: 1 Souza_ArlanLucasde_D.pdf: 3739716 bytes, checksum: 40314e4f8359d3c06b6df822abda9f09 (MD5) Previous issue date: 2002 / Resumo: Uma nova técnica de detecção de vazamentos em tubulação é apresenta­da. A técnica é baseada na análise espectral de sinais de pressão medidos em trechos da tubulação que favorecem a formação de ondas estacionárias. Sendo baseada no espectro do sinal, e não na sua representação temporal, a análise possibilita a identificação do vazamento mesmo quando este ocorre durante a realização de operações que naturalmente encobririam o transiente de pressão gerado pelo vazamento. Testes experimentais foram realizados numa tubulação de 1250 m de comprimento em várias condições de fluxo e configurações de vazamento durante a partida e a parada de uma bomba centrífuga. Transientes de pressão foram obtidos por quatro transdutores conectados a um computador padrão PC comandado por um software de aquisição de dados e controle desenvolvido. Os testes mostraram que a análise espectral dos transientes de pressão, juntamente com o conhecimento dos pontos de reflexão na tubulação, consiste em um método simples e eficiente de se detectar vazamentos / Abstract: The development and test of a technique for leak detection in pipelines is presented. The technique is based on the spectral analysis of pressure signals measured in pipeline sections where the formation of stationary waves is favoured, allowing leakage detection during the start stop of pumps. Experimental tests were performed in a 1250 m long pipeline for various operational conditions of the pipeline (liquid flow rate and leakage configuration). Pressure transients were obtained by four transducers connected to a PC compute. The obtained results show that the spectral analysis of pressure transients, together with the knowledge of reflection points provide a simple and efficient way of identifying leaks during the start stop of pumps in pipelines / Doutorado / Sistemas de Processos Quimicos e Informatica / Doutor em Engenharia Química

El presente trabajo estudia como la participación de los grupos de gasto cambia ante variaciones de los ingresos, teniendo en cuenta que los ingresos son clasificados en transitorios y permanentes. La literatura supone que los hogares deben de considerar la naturaleza de esta clasificación de ingresos y, en consecuencia, los ingresos transitorios deberían de ser asignados a bienes de inversión, mientras que los ingresos permanentes deben de ser asignados a los bienes de consumo. Los resultados muestran que, ante variaciones de estos dos tipos de ingresos, la composición de los grupos de gasto varía, sin embargo, la categoría de alimentos sigue siendo la más representativa. Asimismo, se concluye que, los ingresos transitorios son asignados a los bienes de inversión, pero también son asignados a los bienes de consumo, mostrando que los hogares peruanos todavía no disciernen entre estos dos tipos de ingresos. / This paper studies how the participation of expenditure groups changes in response to changes in income, taking into account that income is classified into transitory and permanent. The literature assumes that households should consider the nature of this income classification and, consequently, transitory income should be allocated to investment goods, while permanent income should be allocated to consumption goods. The results show that, in the face of variations in these two types of income, the composition of the expenditure groups varies; however, the food category continues to be the most representative. It is also concluded that transitory income is allocated to investment goods, but it is also allocated to consumption goods, showing that Peruvian households still do not distinguish between these two types of income. / Trabajo de investigación

[pt] Para simular, prever e controlar os sistemas de energia elétrica, engenheiros precisam de ferramentas computacionais para modelar os componentes dessa rede interconectada altamente complexa. Muitos esforços ao longo do século passado foram dedicados a desenvolver modelos matemáticos para geradores, linhas de transmissão, compensadores de potência reativa, transformadores e assim por diante. Os principais componentes dos sistemas de potência são representados precisamente através de modelos matemáticos, mas as cargas ainda são uma fonte de incerteza nas simulações, devido à sua característica de aleatoriedade. Modelos de carga conservadores superestimam a resposta de potência a desvios de tensão, enquanto modelos de carga excessivamente otimistas podem subestimar as margens de estabilidade, deixando o sistema muito próximo do seu limite operacional. É preciso estabelecer representações de cargas tão próximas da realidade quanto possível, a fim de explorar os recursos de rede de modo mais eficiente. Este trabalho fornece uma metodologia para modelagem de carga, investigando e resumindo as etapas do processo, que podem ser implementadas de diversas maneiras. O tratamento de dados, a escolha de uma representação matemática do modelo de carga e sua estimação de parâmetros são apresentados através de estudos de caso reais, tanto para uma aplicação focada na dinâmica do sistema elétrica, quanto para uma aplicação em regime permanente. Discute-se como otimização e conceitos de inferência estatística podem ser ferramentas efetivas para alcançar melhores aproximações sobre como a carga responderá a perturbações causadas por variações de tensão, sejam estas variações espontâneas, ou devido a ações de controle, ou causadas por curtos-circuitos. / [en] To simulate, predict and control Electric Power Systems (EPS), engineers need tools to model the components of this highly complex interconnected network. Many efforts over the past century were dedicated to develop mathematical models for generators, transmission lines, reactive power compensators, transformers and so on. The main components of the power systems are precisely represented by mathematical models, but the loads are still a source of uncertainty in the simulations, due to their random characteristics. It is well known that conservative load models super estimate power response to voltage deviations, and, on the other hand, over-optimistic load models may underestimate stability margins, leading a system to operate too close to its limit. It is necessary to stablish load representations as close to reality as possible, in order to fully exploit grid resources. This work provides a methodology for load modeling, investigating and summarizing the steps of the process, whose can be implemented in a wide variety of ways. Data treatment, the choice of a load model representation and their parameters estimation are presented through real case studies, both for dynamic simulation and a steady state application. It is discussed how optimization and statistical inference concepts can be effective tools to reach better approximations on how load will respond to disturbances caused by voltage variations, whether these were spontaneous, due to control actions, or caused by short-circuits.