No contexto do método MAC e baseado em esquemas de diferenças finitas, este trabalho apresenta três estudos: i) uma análise de estabilidade, ii) o desenvolvimento de técnicas implícitas e, iii) a construção de métodos de projeção para escoamentos com superfície livre. Na análise de estabilidade, o principal resultado mostra que o método de Crank-Nicolson torna-se condicionalmente estável quando aplicado para uma malha deslocada com a discretiza ção explícita das condições de contorno do tipo Dirichlet. Entretanto, o mesmo método com condições de contorno implícitas é incondicionalmente estável. Para obter métodos mais estáveis, formulações implícitas são desenvolvidas para a equação da pressão na superfície livre, derivada da condição de tensão normal. Esta estratégia resulta no acoplamento dos campos de velocidade e pressão, o que exige a introdução de novos métodos de projeção. Os métodos de projeção assim desenvolvidos resultam em novas metodologias para escoamentos com superfície livre que são apropriados para o tratamento de problemas com baixo número de Reynolds. Além disso, mostra-se que os métodos propostos podem ser aplicados para fluidos viscoelásticos. Novas estratégias são derivadas para obter métodos de projeção de segunda ordem de precisão para escoamentos com superfícies livres. Além dos resultados teóricos sobre a estabilidade de esquemas numéricos, técnicas implícitas e métodos de projeção, testes computacionais são realizados e comparados para consolidação da teoria apresentada. Os resultados numéricos são obtidos no sistema FREEFLOW. A eficiência e robustez das técnicas desenvolvidas neste trabalho são demonstradas na solução de problemas tridimensionais complexos com superfície livre e baixo número de Reynolds, incluindo os problemas do jato oscilante e do inchamento do extrudado / In the context of the MAC method and based on finite difference schemes, this work presents three studies: i) a stability analysis, ii) the development of implicit techniques, and iii) the construction of projection methods for free surface flows. In the stability analysis, the main result shows a precise stability restriction on the Crank-Nicolson method when one uses a staggered grid with Dirichlet explicit boundary conditions. However, the same method with implicit boundary conditions becomes unconditionally stable. In order to obtain more stable methods, implicit formulations are applied for the pressure equation at the free surface, which is derived from the normal stress condition. This approach results in a coupling of the velocity and pressure fields; hence new projection methods for free surface flows need to be developed. The developed projection methods result in new methodologies for low Reynolds number free surface flows. It is also shown that the proposed methods can be applied for viscoelastic fluids. New strategies are derived for obtaining second-order accurate projection methods for free surface flows. In addition to the theoretical results on the stability of numerical schemes, implicit techniques and projection methods, computational tests are carried out and the results compared to consolidate the theory. The numerical results are obtained by the FREEFLOW system. The eficiency and robustness of the techniques in this work are demonstrated by solving complex tridimensional problems involving free surface and low Reynolds numbers, including the jet buckling and the extrudate swell problems

Made available in DSpace on 2015-03-04T18:51:21Z (GMT). No. of bitstreams: 1 Iury.pdf: 2129454 bytes, checksum: 6a7aff5ca085814119b9518b0aab9bef (MD5) Previous issue date: 2010-08-05 / Conselho Nacional de Desenvolvimento Cientifico e Tecnologico / This work aims at developing computational models capable of furnishing more realistic and less costly computationally for the problem of electrokinetic remediation of polluted clayey soils. Innovative results are obtained by improving the multiscale models previously developed by Lima and co-workers through the construction of perturbations of the local microscopic problems in conjuction with more realistic boundary conditions at the electrodes and with the development of precise estimates for the assymptotic behavior of the macroscopic solution. Considering the aliance of such techniques within the framework of the homogenization method of periodic structures we discretize the macroscopic model by the finite element method numerical simulations of an electroosmose experiment capable of predicting more realistic scenarios of electrokinetic remediation. / Este trabalho objetiva o desenvolvimento de modelos computacionais capazes de construir simulações numéricas mais realistas e menos custosas computacionalmente para o problema de descontaminação de solos argilosos por técnicas de eletrocinética. Resultados inovadores são obtidos aprimorando-se os modelos multi-escala desenvolvidos anteriormente por Lima e colaboradores via construção de soluções perturbativas dos problemas locais microscópicos aliada à condições de contorno mais realistas nos eletrodos e ao desenvolvimento de estimativas precisas para o comportamento assintótico da solução macroscópica. Por intermédio da conjunção destas técnicas imersas no contexto da teoria de homogeneização de estruturas periódicas discretizamos o modelo macroscópico pelo método dos elementos finitos e construimos simulações numéricas de um experimento de eletroosmose capazes de predizer cenários mais realistas em eletrorremediação de solos.

Este trabalho visa a avaliar os efeitos dinâmicos em vigas esbeltas semi-infinitas com contato unilateral em apoio elástico sob flexão composta, e sua possível aplicação a dutos flexíveis de extração de petróleo offshore em configuração de catenária (Steel Catenary Risers) provocados por esforços de movimento da plataforma. O estudo se desenvolve com formulação de vigas semi-infinitas em duas dimensões (2D) com suporte unilateral elástico, considerando-se inicialmente apenas efeito de flexão, à qual se aplica o método das múltiplas escalas (MME) para obtenção dos modos de vibração, e posteriormente o método das variedades invariantes (MVI) para recuperação das relações modais entre as coordenadas generalizadas e as variáveis modais. Trata-se, a seguir, da formulação do problema com consideração da tração, tanto estática como dinâmica, e a projeção da equação de movimento completa segundo os modos obtidos no problema da flexão simples, para obtenção do modelo de ordem reduzida com efeito de tração, que permitiu estudar este último qualitativa e quantitativamente. A formulação inicial teve como ponto de partida a utilização de mudanças de variáveis que permitiu a transformação do problema de condições móveis em condições fixas de contorno, propiciando caracterizar o movimento vertical da viga na solução analítica do problema de vibrações livres, até o limite possível, para depois utilizar técnicas de integração numérica na solução do problema de vibrações forçadas. A motivação para estudo de tais efeitos está ligada à importância indiscutível que tem para a indústria petrolífera e para a economia mundial, além do crescente papel deste segmento em termos nacionais. A descoberta do Campo de Tupi, na Bacia de Santos, que representa a maior reserva de óleo e gás do País, alavancou pesquisas para extração de petróleo em águas ultra-profundas (profundidades da ordem de 2.500m) (fonte: www.petrobras.com.br), além da necessidade já estabelecida de maiores pesquisas em águas profundas (profundidades que podem chegar a 2.000m). Considere-se também a motivação da própria Petrobrás, cujo intuito é ter capacidade de extrair 4,5 milhões de barris por dia até 2020. Tal capacidade de produção deve exigir 45 sistemas de produção, sendo que cada sistema de produção necessita de 4 a 5 barcos de apoio, totalizando assim 200 barcos operando em 2020 (Fonte: O Estado de São Paulo, Ano 131, Nº 42786, Coluna Celso Ming). Apesar da crescente pesquisa em fontes alternativas de energia, a demanda não só por óleo e gás, mas também pela utilização de polímeros e componentes plásticos, ainda sustentará a extração de petróleo por no mínimo algumas décadas. Além disso, técnicas melhores na fixação de carbono (técnicas naturais e artificiais para seqüestro de carbono), permitem sonhar que a utilização de petróleo, desde a sua queima até a produção de material sintético possa se tornar cada vez menos impactante para o meio ambiente. / This work aims at analyzing the dynamical effects in slender reticulated beams on elastic support with unilateral contact under combined axial force and bending, and the possibility of its association with steel catenary risers connected to offshore structures, caused by movement of the platform. The study is developed based on 2D semi-infinite beams on unilateral elastic support, initially considering only bending effects, to which the multiple scales method (MSM) is applied to obtain vibration modes, and thereafter the invariant manifold method (IMM) to obtain the modal relationships between the generalized coordinates and the modal variables. It follows the formulation of the problem considering tension effects, both from static and dynamical sources, and the projection of the full equation of motion with respect to the vibration modes from the first formulation considering only bending effects, which allows us to obtain the reduced order model with tension effects, and its application to steel catenary risers (SCR). The starting point for the mathematical formulation is a variable transformation, which allows us to change a free boundary into a fixed boundary one. It is then possible to obtain the free-vibration analytical solution, as far as possible, and then use numerical integration schemes to obtain the forced vibration time responses. As motivation, we can firstly mention the unquestionable role of the oil and gas industry for the world economy, besides its increasing market share within the national economy. The discovery of the Tupi Field, in the Santos Basin, which represents the biggest oil and gas reserve of the Country, gave impulse to investments and research in all fields associated to the oil and gas area, in ultra deep water (depth around 2.500m), besides the well known need of research in deep water (depth around 2.000m). In spite of the growing research on alternative sources of energy, the world still demands not only oil and gas, but also polymers and plastic components, and this will still be the scenario for the next decades. Besides that, the increasing of knowledge in carbon sequestration may keep up the use of oil and gas for some more decades, hopefully with less environment impact.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / Como eventos de fissão induzida por nêutrons não ocorrem nas regiões nãomultiplicativas de reatores nucleares, e.g., moderador, refletor, e meios estruturais, essas regiões não geram potência e a eficiência computacional dos cálculos globais de reatores nucleares pode portanto ser aumentada eliminando os cálculos numéricos explícitos no interior das regiões não-multiplicativas em torno do núcleo ativo. É discutida nesta dissertação a eficiência computacional de condições de contorno aproximadas tipo albedo na formulação de ordenadas discretas (SN) para problemas de autovalor a dois grupos de energia em geometria bidimensional cartesiana. Albedo, palavra de origem latina para alvura, foi originalmente definido como a fração da luz incidente que é refletida difusamente por uma superfície. Esta palavra latina permaneceu como o termo científico usual em astronomia e nesta dissertação este conceito é estendido para reflexão de nêutrons. Este albedo SN nãoconvencional substitui aproximadamente a região refletora em torno do núcleo ativo do reator, pois os termos de fuga transversal são desprezados no interior do refletor. Se o problema, em particular, não possui termos de fuga transversal, i.e., trata-se de um problema unidimensional, então as condições de contorno albedo, como propostas nesta dissertação, são exatas. Por eficiência computacional entende-se analisar a precisão dos resultados numéricos em comparação com o tempo de execução computacional de cada simulação de um dado problema-modelo. Resultados numéricos para dois problemas-modelo com de simetria são considerados para ilustrar esta análise de eficiência. / As neutron fission events do not take place in the non-multiplying regions of nuclear reactors, e.g., moderator, reflector, and structural core, these regions do not generate power and the computational efficiency of nuclear reactor global calculations can hence be improved by eliminating the explicit numerical calculations within the non-multiplying regions around the active domain. Discussed here is the computational efficiency of approximate discrete ordinates (SN) albedo boundary conditions for two-energy group eigenvalue problems in X,Y geometry. Albedo, the Latin word for whiteness, was originally defined as the fraction of incident light reflected diffusely by a surface. This Latin word has remained the usual scientific term in astronomy and in this dissertation this concept is extended for the reflection of neutrons. The non-standard SN albedo substitutes approximately the reflector region around the active domain, as we neglect the transverse leakage terms within the nonmultiplying reflector. Should the problem have no transverse leakage terms, i.e., onedimensional slab geometry, then the offered albedo boundary conditions are exact. By computational efficiency we mean analyzing the accuracy of the numerical results versus the CPU execution time of each run for a given model problem. Numerical results to two symmetric test problems are shown to illustrate this efficiency analysis.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Um código computacional para escoamentos bifásicos incorporando metodologia híbrida entre oMétodo dos Elementos Finitos e a descrição Lagrangeana-Euleriana Arbitrária do movimento é usado para simular a dinâmica de um jato transversal de gotas na zona primária de quebra. Os corpos dispersos são descritos por meio de um método do tipo front-tracking que produz interfaces de espessura zero através de malhas formadas pela união de elementos adjacentes em ambas as fases e de técnicas de refinamento adaptativo. Condições de contorno periódicas são implementadas de modo variacionalmente consistente para todos os campos envolvidos nas simulações apresentadas e uma versão modificada do campo de pressão é adicionada à formulação do tipo um-fluido usada na equação da quantidade de movimento linear. Simulações numéricas diretas em três dimensões são executadas para diferentes configurações de líquidos imiscí veis compatíveis com resultados experimentais encontrados na literatura. Análises da hidrodinâmica do jato transversal de gotas nessas configurações considerando trajetórias, variação de formato de gota, espectro de pequenas perturbações, além de aspectos complementares relativos à qualidade de malha são apresentados e discutidos. / A two-phase flow computational code taking a hybrid Arbitrary Lagrangian-Eulerian desciption of movement along with the Finite Element Method is used to simulate the dynamics of an incompressible drop jet in crossflow in the primary breakup zone. Dispersed entities are described by means of a front-tracking method which produces zero-thickness interfaces through contiguous element meshing and adaptive refinement techniques. Periodic boundary conditions are implemented in a variationally consistent way for all the scalar fields involved in the presented simulations and amodified version of the pressure field is added to the one-fluid formulation employed in the momentum equation. Three-dimensional direct numerical simulations for different flow configurations of immiscible liquids pertinent to experimental results found in literature. Analyses of the hydrodynamics of the drop jet in crossflow in these configurations considering trajectories, drop shape variations, spectrum of small disturbances, besides additional aspects relating to mesh quality are presented and discussed.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / Como eventos de fissão induzida por nêutrons não ocorrem nas regiões nãomultiplicativas de reatores nucleares, e.g., moderador, refletor, e meios estruturais, essas regiões não geram potência e a eficiência computacional dos cálculos globais de reatores nucleares pode portanto ser aumentada eliminando os cálculos numéricos explícitos no interior das regiões não-multiplicativas em torno do núcleo ativo. É discutida nesta dissertação a eficiência computacional de condições de contorno aproximadas tipo albedo na formulação de ordenadas discretas (SN) para problemas de autovalor a dois grupos de energia em geometria bidimensional cartesiana. Albedo, palavra de origem latina para alvura, foi originalmente definido como a fração da luz incidente que é refletida difusamente por uma superfície. Esta palavra latina permaneceu como o termo científico usual em astronomia e nesta dissertação este conceito é estendido para reflexão de nêutrons. Este albedo SN nãoconvencional substitui aproximadamente a região refletora em torno do núcleo ativo do reator, pois os termos de fuga transversal são desprezados no interior do refletor. Se o problema, em particular, não possui termos de fuga transversal, i.e., trata-se de um problema unidimensional, então as condições de contorno albedo, como propostas nesta dissertação, são exatas. Por eficiência computacional entende-se analisar a precisão dos resultados numéricos em comparação com o tempo de execução computacional de cada simulação de um dado problema-modelo. Resultados numéricos para dois problemas-modelo com de simetria são considerados para ilustrar esta análise de eficiência. / As neutron fission events do not take place in the non-multiplying regions of nuclear reactors, e.g., moderator, reflector, and structural core, these regions do not generate power and the computational efficiency of nuclear reactor global calculations can hence be improved by eliminating the explicit numerical calculations within the non-multiplying regions around the active domain. Discussed here is the computational efficiency of approximate discrete ordinates (SN) albedo boundary conditions for two-energy group eigenvalue problems in X,Y geometry. Albedo, the Latin word for whiteness, was originally defined as the fraction of incident light reflected diffusely by a surface. This Latin word has remained the usual scientific term in astronomy and in this dissertation this concept is extended for the reflection of neutrons. The non-standard SN albedo substitutes approximately the reflector region around the active domain, as we neglect the transverse leakage terms within the nonmultiplying reflector. Should the problem have no transverse leakage terms, i.e., onedimensional slab geometry, then the offered albedo boundary conditions are exact. By computational efficiency we mean analyzing the accuracy of the numerical results versus the CPU execution time of each run for a given model problem. Numerical results to two symmetric test problems are shown to illustrate this efficiency analysis.

This work aims at developing computational models capable of furnishing more realistic and less costly computationally for the problem of electrokinetic remediation of polluted clayey soils. Innovative results are obtained by improving the multiscale models previously developed by Lima and co-workers through the construction of perturbations of the local microscopic problems in conjuction with more realistic boundary conditions at the electrodes and with the development of precise estimates for the assymptotic behavior of the macroscopic solution. Considering the aliance of such techniques within the framework of the homogenization method of periodic structures we discretize the macroscopic model by the finite element method numerical simulations of an electroosmose experiment capable of predicting more realistic scenarios of electrokinetic remediation. / Este trabalho objetiva o desenvolvimento de modelos computacionais capazes de construir simulações numéricas mais realistas e menos custosas computacionalmente para o problema de descontaminação de solos argilosos por técnicas de eletrocinética. Resultados inovadores são obtidos aprimorando-se os modelos multi-escala desenvolvidos anteriormente por Lima e colaboradores via construção de soluções perturbativas dos problemas locais microscópicos aliada à condições de contorno mais realistas nos eletrodos e ao desenvolvimento de estimativas precisas para o comportamento assintótico da solução macroscópica. Por intermédio da conjunção destas técnicas imersas no contexto da teoria de homogeneização de estruturas periódicas discretizamos o modelo macroscópico pelo método dos elementos finitos e construimos simulações numéricas de um experimento de eletroosmose capazes de predizer cenários mais realistas em eletrorremediação de solos.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Um código computacional para escoamentos bifásicos incorporando metodologia híbrida entre oMétodo dos Elementos Finitos e a descrição Lagrangeana-Euleriana Arbitrária do movimento é usado para simular a dinâmica de um jato transversal de gotas na zona primária de quebra. Os corpos dispersos são descritos por meio de um método do tipo front-tracking que produz interfaces de espessura zero através de malhas formadas pela união de elementos adjacentes em ambas as fases e de técnicas de refinamento adaptativo. Condições de contorno periódicas são implementadas de modo variacionalmente consistente para todos os campos envolvidos nas simulações apresentadas e uma versão modificada do campo de pressão é adicionada à formulação do tipo um-fluido usada na equação da quantidade de movimento linear. Simulações numéricas diretas em três dimensões são executadas para diferentes configurações de líquidos imiscí veis compatíveis com resultados experimentais encontrados na literatura. Análises da hidrodinâmica do jato transversal de gotas nessas configurações considerando trajetórias, variação de formato de gota, espectro de pequenas perturbações, além de aspectos complementares relativos à qualidade de malha são apresentados e discutidos. / A two-phase flow computational code taking a hybrid Arbitrary Lagrangian-Eulerian desciption of movement along with the Finite Element Method is used to simulate the dynamics of an incompressible drop jet in crossflow in the primary breakup zone. Dispersed entities are described by means of a front-tracking method which produces zero-thickness interfaces through contiguous element meshing and adaptive refinement techniques. Periodic boundary conditions are implemented in a variationally consistent way for all the scalar fields involved in the presented simulations and amodified version of the pressure field is added to the one-fluid formulation employed in the momentum equation. Three-dimensional direct numerical simulations for different flow configurations of immiscible liquids pertinent to experimental results found in literature. Analyses of the hydrodynamics of the drop jet in crossflow in these configurations considering trajectories, drop shape variations, spectrum of small disturbances, besides additional aspects relating to mesh quality are presented and discussed.

Este trabalho visa a avaliar os efeitos dinâmicos em vigas esbeltas semi-infinitas com contato unilateral em apoio elástico sob flexão composta, e sua possível aplicação a dutos flexíveis de extração de petróleo offshore em configuração de catenária (Steel Catenary Risers) provocados por esforços de movimento da plataforma. O estudo se desenvolve com formulação de vigas semi-infinitas em duas dimensões (2D) com suporte unilateral elástico, considerando-se inicialmente apenas efeito de flexão, à qual se aplica o método das múltiplas escalas (MME) para obtenção dos modos de vibração, e posteriormente o método das variedades invariantes (MVI) para recuperação das relações modais entre as coordenadas generalizadas e as variáveis modais. Trata-se, a seguir, da formulação do problema com consideração da tração, tanto estática como dinâmica, e a projeção da equação de movimento completa segundo os modos obtidos no problema da flexão simples, para obtenção do modelo de ordem reduzida com efeito de tração, que permitiu estudar este último qualitativa e quantitativamente. A formulação inicial teve como ponto de partida a utilização de mudanças de variáveis que permitiu a transformação do problema de condições móveis em condições fixas de contorno, propiciando caracterizar o movimento vertical da viga na solução analítica do problema de vibrações livres, até o limite possível, para depois utilizar técnicas de integração numérica na solução do problema de vibrações forçadas. A motivação para estudo de tais efeitos está ligada à importância indiscutível que tem para a indústria petrolífera e para a economia mundial, além do crescente papel deste segmento em termos nacionais. A descoberta do Campo de Tupi, na Bacia de Santos, que representa a maior reserva de óleo e gás do País, alavancou pesquisas para extração de petróleo em águas ultra-profundas (profundidades da ordem de 2.500m) (fonte: www.petrobras.com.br), além da necessidade já estabelecida de maiores pesquisas em águas profundas (profundidades que podem chegar a 2.000m). Considere-se também a motivação da própria Petrobrás, cujo intuito é ter capacidade de extrair 4,5 milhões de barris por dia até 2020. Tal capacidade de produção deve exigir 45 sistemas de produção, sendo que cada sistema de produção necessita de 4 a 5 barcos de apoio, totalizando assim 200 barcos operando em 2020 (Fonte: O Estado de São Paulo, Ano 131, Nº 42786, Coluna Celso Ming). Apesar da crescente pesquisa em fontes alternativas de energia, a demanda não só por óleo e gás, mas também pela utilização de polímeros e componentes plásticos, ainda sustentará a extração de petróleo por no mínimo algumas décadas. Além disso, técnicas melhores na fixação de carbono (técnicas naturais e artificiais para seqüestro de carbono), permitem sonhar que a utilização de petróleo, desde a sua queima até a produção de material sintético possa se tornar cada vez menos impactante para o meio ambiente. / This work aims at analyzing the dynamical effects in slender reticulated beams on elastic support with unilateral contact under combined axial force and bending, and the possibility of its association with steel catenary risers connected to offshore structures, caused by movement of the platform. The study is developed based on 2D semi-infinite beams on unilateral elastic support, initially considering only bending effects, to which the multiple scales method (MSM) is applied to obtain vibration modes, and thereafter the invariant manifold method (IMM) to obtain the modal relationships between the generalized coordinates and the modal variables. It follows the formulation of the problem considering tension effects, both from static and dynamical sources, and the projection of the full equation of motion with respect to the vibration modes from the first formulation considering only bending effects, which allows us to obtain the reduced order model with tension effects, and its application to steel catenary risers (SCR). The starting point for the mathematical formulation is a variable transformation, which allows us to change a free boundary into a fixed boundary one. It is then possible to obtain the free-vibration analytical solution, as far as possible, and then use numerical integration schemes to obtain the forced vibration time responses. As motivation, we can firstly mention the unquestionable role of the oil and gas industry for the world economy, besides its increasing market share within the national economy. The discovery of the Tupi Field, in the Santos Basin, which represents the biggest oil and gas reserve of the Country, gave impulse to investments and research in all fields associated to the oil and gas area, in ultra deep water (depth around 2.500m), besides the well known need of research in deep water (depth around 2.000m). In spite of the growing research on alternative sources of energy, the world still demands not only oil and gas, but also polymers and plastic components, and this will still be the scenario for the next decades. Besides that, the increasing of knowledge in carbon sequestration may keep up the use of oil and gas for some more decades, hopefully with less environment impact.

This paper presents a methodology for determining eigenfunctions and frequencies of the Euler-Bernoulli model for elastic beams that can include damping and devices located at intermediate or end points of the beam. The eigenfunctions or vibration modes of the beam are obtained by using solution basis generated by the dynamic solution of a fourth-order differential equation, through a block matrix formulation of the boundary and compatibility conditions. The use of the dynamic basis has been often used to reduce the calculations in obtaining the modes and frequencies. Forced responses are obtained with the Galerkin method by modifying the classical modal analysis with the inclusion of new conditions of orthogonality between modes that are suitable for problems with viscous damping or non-classical boundary conditions. / Este trabalho apresenta uma metodologia para determinar as autofunções e as frequências de um modelo Euler-Bernoulli para vigas elásticas que podem incluir amortecimento e dispositivos localizados num ponto intermediário ou nos extremos da viga. As autofunções ou modos de vibração da viga são obtidos usando uma base de solução gerada pela solução dinâmica de uma equação diferencial de quarta ordem, através de uma formulação matricial em blocos para as condições de contorno e de compatibilidade. O uso da base dinâmica tem sido frequentemente utilizada para reduzir os cálculos na obtenção dos modos e das frequências. Respostas forçadas são obtidas usando o método de Galerkin, modificando a análise modal clássica com a inclusão de novas condições de ortogonalidade entre os modos que são adequadas para problemas com amortecimento viscoso ou com condições de contorno não-clássicas