[pt] Neste trabalho a análise do comportamento de sistemas rotativos do tipo eixo-rotormancal é estendida para incluir os efeitos da presença de mancais hidrodinâmicos na resposta dinâmica. Estes efeitos estão associados à não-linearidade da força de reação exercida pelos suportes sobre o eixo e dependem dos deslocamentos, velocidades transversais e da rotação própia do rotor. A modelagem estrutural do sistema é obtida empregando-se o método dos elementos finitos. O eixo é representado pelo modelo de viga de Timoshenko com dois nós, quatro graus-de-liberdade por nó, e a interpolação do campo de deslocamentos é obtida utilizando-se as funções de Hermite. Os rotores são modelados empregando-se elementos de inércia concentrada associada aos graus-de-liberdade de um ponto nodal do modelo. E, na representação dos mancais hidrodinâmicos utilizou-se a equação de Reynolds, com as hipóteses simplificadoras para mancais curtos, obtendo-se a solução para a distribuição de pressão do filme de óleo em forma fechada. Essa distribuição de pressão permite a obtenção dos coeficientes das matrizes e rigidez e de amortecimento associadas aos graus de liberdade do eixo no ponto nodal de representação do mancal. Para a integração temporal do sistema de equações diferencias utiliza-se o procedimento passo-a-passo, tendo-se implementado os métodos implícitos de Newmark e Wilson – teta, na forma incondicionalmente estável. Devido à não-linearidade das equações obtidas com a presença dos mancais hidrodinâmicos, em cada intervalo de tempo utiliza-se o procedimento de Newton-Raphson modificado para a correção da solução numérica obtida com outros resultados analíticos/numéricos disponíveis na literatura. Também, uma representação numérica para mancais hidrodinâmicos segmentados é apresentada, utilizando-se o desenvolvimento teórico para mancais simples. Neste caso a avaliação do procedimento numérico é fornecida comparando-se a solução numérica com resultados experimentais obtidos dos rotores de usina hidrogenada avaliada pelo CEPEL. Em ambos os procedimentos o rotor idealizado de jeffcott é empregado no estudo de casos. Verifica-se que os principais resultados associados aos efeitos da precessão auto-excitada (oil whirl), de chicoteamento (oil whip), e da estabilização dinâmica do sistema são reproduzidos pelos modelos numéricos utilizados. / [en] In this work a formulation for the analysis of shaft-rotor-bearing type rotating systems is extendend to accommodate the effects of hydrodynamic bearings in its dynamic response. These effects, which are associated to the nonlinear force on the shaft at the bearings, are dependent of the transverse displacements, transverse linear velocities an the angular veolicty of the shaft. The structure behavior is modeled by employing the finite element method. The shaft is represented by the two node timoshenko model for bearns, with four desgrees-of-freedom per node and Hermite interpolation functions to represent the displacement fields along the bearn axis. Rotors are modeled by using concentrated inertia elements associated to the shaft degrees-of-freedom at the bearing nodal point. In the numerical analysis considering the time integration of the system differential equation, a step-by-step procedure was employed with the newmark technique in this unconditionally stable form. Due to the nonlearities associated with the hydrodynamic bearings, the solution of the system of equations is obtained using a modified Newton-Raphson precedure at each time step for solution convergence. In the evaluation of the proposed computacional system, comparison with solutions obtained from analytical/numerical results available in the literature are used. Also, a numeric represemtation of tilting-pad bearings is presented using the theory for plain journal bearings, under the same simplified conditions. In this case an evaluation of the numerical procedure is given by comparing calculated solutions with experimental results obtained from the evaluation of a hydrogenaration plant provided by CEPEL-Brazilian Research Center For Eletrobras. In both plain an tilting-pad journal bearing numerical procedures, the idealized Jeffcott rotor is employed as a case study for different operating conditions. As a result, it is shown that the solutions associated to the main oil whirl and oil whip effects and afterwards dynamic stabilization are represented by the proposed numerical procedures employed.

[pt] Este trabalho descreve uma metodologia para modelagem de cargas elétricas, utilizando dados de tensão e corrente registrados durante distúrbios no sistema de potência. Estes modelos são utilizados na representação de subsistemas da rede elétrica em simulações computacionais que preveem o comportamento dinâmico do sistema de potência após perturbações em suas condições normais de operação.São apresentados resultados práticos da metodologia proposta, onde a carga é definida como um sistema cuja saída é a variação da potência consumida e a entrada é a variação da tensão, ambas medidas em barramentos de 69 kV da Companhia Hidroelétrica do São Francisco (CHESF), ponto de entrega de energia para concessionárias distribuidoras de energia do nordeste brasileiro. Estas distribuidoras são modeladas como cargas, supridas pelo sistema de transmissão da CHESF e todos os elementos consumidores de energia são agregados nestes modelos equivalentes, parametrizados para simular o maior número de contingências típicas medidas em cada um destes barramentos de carga.A técnica de estimação de parâmetros dos modelos de cargas é o Algoritmo Genético (AG) cujos resultados apresentaram precisão para a simulação de contingências de características bem distintas, caracterizando a abrangência alcançada no processo de identificação de sistemas.Ao final do trabalho são apresentadas curvas de desvios de potência ativa e reativa causadas por afundamentos de tensão, ambos registrados nos barramentos das subestações da CHESF. Estas curvas foram utilizadas para estimar os parâmetros dos modelos, obtidos individualmente para cada uma das subestações estudadas. / [en] This work describes a measurement-based load modeling methodology, using voltage and current data registered during power system disturbances. These load models are used on computational simulations for predicting power system stability after disturbances of system operational points. Practical results are presented of the proposed methodology, defining load as a system whose output is power deviation from its operational state and input is voltage sags, both measured at 69 kV bus bars of São Francisco Hydroelectric Company (CHESF), points of common coupling (PCC) between CHESF and local distribution utilities. Therefore, distribution utilities are seen as loads supplied by CHESF’s transmission system. All devices consuming power from the PCC are aggregated into an equivalent model, parameterized to simulate most typical contingencies measured by these 69kV load bars. Optimization technique used for load model parameter estimation is Genetic Algorithm (GA), showing his flexibility on implementation and good coverage and accuracy in the final results. At the end, it will be presented a set of active and reactive power curves during and after voltages sags, measured on CHESF’s substations. These curves were used as estimation data to parameterize load models for each substation chosen.

[pt] O tema abordado neste trabalho é a observação e análise, em regime permanente e dinâmico, da ocorrência de um fenômeno que já foi observado em condições reais de operação do sistema elétrico brasileiro, que é a relação oposta à usual entre a grandeza controlada e a grandeza controladora. Nestes caso, mesmo que haja margem de recursos para manter a tensão controlada, ela não é útil. Por exemplo, uma diminuição na relação de transformação num transformador de tapes variáveis, com intuito de aumentar a tensão controlada acaba por reduzí-la, até que os limites de troca de tapes sejam atingidos ou o sistema entre em colapso. Para demonstrar a existência do problema, foram executadas simulações, em regime permanente e dinâmico, e verificado o efeito do controle de tensão por um transformador com tapes variáveis e por compensadores estáticos de potência reativa, situações corriqueiras de um sistema de potência. Foram demonstradas situações em que foi possível verificar a mudança da região de operação. Para a análise em regime permanente foi utilizado um algoritmo de fluxo de carga, e para a análise dinâmica, uma simulação no domínio do tempo. Nas simulações envolvendo transformadores de tapes variáveis, foi possível verificar a existência o efeito reverso da ação de controle de tensão nas análises estática e dinâmica. Nas simulações utilizando o compensador estático de potência reativa houve divergência entre os resultados das duas análises. / [en] This work presents the observation and analysis, in steady state and dynamic performance, of the phenomenon already observed in real operation conditions of the Brazilian Electric System, which is the opposite relationship between the controlled value and the target value. In this case, even if the resources have margin to keep the voltage controlled, this is not useful. For example, the reduction in the turn ratio on load tap changer transformer, with the aim of increase in the controlled voltage, result in its reduction, until the tap changer limit is reached or the system is led to the collapse. To demonstrate the existence of this problem, steady state and dynamic performance simulations were done, and the voltage control effect by on load tap changer transformer and static var compensator , current situations in a power system. Several situations where is possible verify the operation region changing was demonstrated. For steady state analysis was used a load flow algorithm and, for the dynamic analysis, a time domain simulation. In the simulations with on load tap changer transformer, it was possible to verify the existence of the reverse effect of the voltage control action in the static and dynamic analysis. In the simulations using static var compensator, a divergence was found between the results in the two analyses.

[pt] Centros de operação ao redor do mundo têm verificado o aumento da complexidade ao longo dos anos, principalmente pela grande quantidade de novas informações que devem ser consideradas. No sistema elétrico brasileiro não é diferente, impactado pela participação cada vez mais ativa de fontes renováveis, como eólicas e fotovoltaicas, além de elos de corrente contínua. Seguindo uma tendência global, a consideração de valores de inércia sincronizada para a avaliação do desempenho da frequência e de fluxos sistêmicos se faz necessária, configurando-se como uma nova informação a ser exposta à operação em tempo real. Visando garantir mais segurança e confiabilidade ao ambiente de tempo real, o presente trabalho tem como objetivo traduzir os estudos realizados pela equipe de planejamento do ONS para a operação, através de um sistema capaz de monitorar parâmetros como a inércia sincronizada, reserva de potência operativa e outras informações associadas, ampliando a consciência situacional dos usuários. Além disso, o sistema deve promover simulações de contingências severas que possam provocar atuações de esquemas de corte de carga e indicar ações que possam ser realizadas para mitigar o problema. / [en] Operation centers around the world have verified an increase of the complexity throughout the years, mostly due to huge amounts of new information that must be considered. The Brazilian electric system is no different, impacted by the increasingly active participation of renewable sources, such as wind and photovoltaic, as well as direct current links. Following a global trend, the consideration of values of synchronous inertia for the evaluation of the performance of the frequency is necessary, configuring as a new information to be exposed to real-time operation. With the purpose of guaranteeing more safety and reliability to the real-time environment, this work has the goal of translating the studies carried out by the ONS planning team to the operation, through a system capable of monitoring parameters such as synchronous inertia, operating reserve and other associated informations, increasing users situational awareness. Besides that, the system must promote simulations of severe contingencies that may trigger load shedding schemes and indicate actions that can be taken to soften the problem.