Este trabalho apresenta um estudo experimental do efeito da razão de bloqueio em um canal aerodinâmico no Número de Strouhal e o fenômeno da biestabilidade. É estudado o escoamento turbulento sobre um cilindro e também sobre dois cilindros dispostos lado a lado, onde fez-se uso de cilindros de vários diâmetros a fim de variar a razão de bloqueio no canal. Dois tipos de espaçamento entre os tubos são estudados. A técnica experimental consiste na medição de flutuações de velocidades em um canal aerodinâmico utilizando a técnica de anemometria de fio quente. Os dados obtidos da medição no canal aerodinâmico são tratados com o uso de ferramentas estatísticas, espectrais e de ondaletas. Para um cilindro, os resultados confirmam o aumento do número de Strouhal com o aumento da razão de bloqueio, na faixa do número de Reynolds utilizada. Razões de bloqueio mais baixas podem acarretar diminuição do número de Strouhal, também em função do número de Reynolds. Os resultados confirmam a existência do fenômeno da biestabilidade no escoamento sobre dois cilindros dispostos lado a lado, assim como a influência da razão de bloqueio sobre o fenômeno, sendo mais predominante em um dos espaçamentos estudados. Para razões de bloqueio altas, o fenômeno da biestabilidade pode não ocorrer ou pouco ocorre. Para razões de bloqueio mais baixas podemos ter também a não ocorrência do fenômeno em função das dimensões dos tubos serem mais reduzidas, fazendo com que a captação do fenômeno também seja reduzida. / This paper presents an experimental study of the effect of blockage ratio in a aerodynamic channel on Strouhal number and the bistability phenomenon. It studied the turbulent flow around a single cylinder and also on two cylinders arranged side-by-side, where use has been made of cylinders of various diameters to vary the blockage ratio of the channel. Two types of pitch to diameter ratio are studied. The experimental technique consists of measuring velocity fluctuations in aerodynamic channel using the technique of hot-wire anemometry. The data obtained from measuring the aerodynamic channel are treated with the use of statistical tools, spectral and wavelet analysis. For a single cylinder, the results confirm the increase of the Strouhal number with increasing blockage ratio in the range of Reynolds number used. Lower blockage ratio can cause a decrease of the Strouhal number, also a function of Reynolds number. The results confirm the existence of the bistability phenomenon in the flow on two cylinders arranged side-by-side, as well as the influence of blockage ratio on the phenomenon, being more predominant in one of the row spacings. For higher blockage ratios, the bistability phenomenon may not occur or occurs shortly. For lower blockage ratio, due to the smaller dimensions of the tubes the phenomenon may not occur on the capture of the phenomenon is also reduced with the present experimental technique.

É conhecida ao longo da história da humanidade a grande quantidade de acidentes vinculados à ação de cargas cíclicas. Há mais de 100 anos, metodologias levam em conta o fenômeno de fadiga em projetos de novos componentes. As metodologias existentes para avaliar a vida em fadiga são empíricas e têm sido aplicadas com sucesso na determinação da vida do componente em fadiga quando o mesmo é submetido a carregamento uniaxial ou a carregamento multiaxial proporcional. Mas a evidencia experimental tem mostrado que quando a solicitação é multiaxial e não proporcional, as leis antes mencionadas deixam de ser adequadas. Dessa forma os critérios clássicos utilizados não preveem corretamente a vida do componente. Neste contexto o presente trabalho avalia algumas das metodologias de fadiga multiaxial não proporcionais disponíveis na bibliografia especializada. Os resultados obtidos com estes métodos são comparados entre sim e com os resultados obtidos utilizando a metodologia de fadiga multiaxial proporcional. Para comparar as metodologias citadas é apresentado um exemplo de solicitação simples, e também a análise do componente de uma máquina agrícola solicitada por um histórico de tensões típicas, obtida a partir do teste da máquina em campo. Finalmente são apresentadas as conclusões sobre as metodologias utilizadas e sobre a importância de levar em conta a não proporcionalidade de um estado de tensões multiaxial oscilantes. / It is known throughout human history the large amount of accidents linked to action of cyclic loading. For over 100 years, methodologies take into account the fatigue phenomenon of new components projects. Existing methodologies for assessing the fatigue life are empirical and have been successfully applied in the determination of fatigue life of the component when it is subjected to uniaxial loading or proportional multiaxial loading. But the experimental evidence has shown that when the request is multiaxial and not proportional, the laws mentioned above are no longer adequate. Thus the classical criteria used not correctly predict component life. In this context, this paper analyzes some of the methodologies of multiaxial fatigue non proportionate available in the relevant literature. The results obtained with these methods are compared with the results so obtained and using the proportional multiaxial fatigue methodology. To compare the methods mentioned is an example of simple request, and also the analysis component of an agricultural machine requested by a history of typical strains, obtained from the machine field testing on a test track. Finally conclusions on the methodologies used and the importance of taking into account the non-proportionality of a state of oscillating multiaxial stresses are presented.

Taludes escavados em filitos muitas vezes estão sujeitos à instabilidade, principalmente quando há outras famílias de descontinuidades presentes na estrutura. Este problema está presente na mina de Candiota, onde foi identificado, em um dos taludes, potencialidade para ruptura por tombamento. Este trabalho aborda uma discussão sobre métodos analíticos e numéricos aplicado a tombamento primário, com objetivo de avaliar o fator de segurança do talude em estudo por meio de método numérico por elementos finitos, utilizando o Phase2 da Rocscience, e métodos analíticos. Também são apresentadas doze simulações para avaliar se o movimento de tombamento em modelagens numéricas obedece a um efeito de escala. Com uso de métodos numéricos dois modelos foram construídos: um contínuo equivalente e outro contínuo (mas com as descontinuidades inseridas no modelo). No modelo com as descontinuidades a rocha foi tratada como elástica por meio do critério de Hoek-Brown e as descontinuidades foram tratadas como elasto-plásticas por meio do critério de Mohr-Coulomb. Os resultados das modelagens da mina de Candiota foram semelhantes e mostraram fatores de segurança em níveis de estabilidade, sendo que no modelo contínuo (mas com descontinuidades inseridas) se observou de forma mais clara as tensões cisalhantes induzidas nas descontinuidades quando ocorre o movimento de tombamento. O uso dos métodos analíticos mostrou-se ineficientes para o problema apresentado. Os resultados das doze modelagens sugerem que o aumento do espaçamento das descontinuidades influencia nos fatores de resistência dos taludes e na forma do movimento de tombamento. / Slopes in phyllites are often subject to instability, especially when there are other families of discontinuities in the structure. This problem is present in Candiota mine, which was identified in one of the slopes potential to toppling. This work show a discussion of analytical and numerical methods applied to toppling, to evaluate the slope safety factor being studied by numerical method finite element using the Phase2 of Rocscience, and analytical methods. Also shown are twelve simulations to evaluate if the toppling in numerical modeling follows a scaling effect. With use of numerical methods two models were built: an equivalent continuous and another solid (but with discontinuities inserted in the model). In the model with the discontinuities the rock was treated as elastic by the Hoek-Brown criterion and discontinuities were treated as elastic-plastic by the Mohr-Coulomb criterion. The results of Candiota mine were similar and showed safety factors in stability levels, and in the continuous model (with discontinuities) was observed more clearly the shear stress induced in discontinuities when the movement of overturning occurs. The use of analytical methods proved inefficient for the problem presented. The results of the twelve modeling suggests increasing the spacing of the discontinuities of slope influences the resistance factor and the form of the toppling movement.

Esta dissertação apresenta um estudo da caracterização numérica do comportamento mecânico de um músculo esquelético estriado sujeito à combinação de diferentes contrações. A partir dos mecânismos fisiológicos da contração do tecido muscular esquelético, são apresentadas as principais características da produção de força do músculo em diferentes tipos de contrações (isométrica, concêntrica, excêntrica) e suas combinações. A representação numérica da resposta mecânica do músculo esquelético foi investigada por diferentes autores. Dos diversos modelos encontrados na literatura, alguns foram escolhidos para a continuidade deste estudo, segundo sua capacidade representativa e facilidade de implementação em métodos numéricos de solução para grandes deformações. Os modelos disponíveis na literatura foram implementados e avaliados para diferentes sequências de contrações, sendo também realizado um ajuste de parâmetros para dados experimentais de contrações isométricas e isométricas-excêntricas-isométricas. Nenhuma das propostas avaliadas apresentou resultados satisfatórios. A principal deficiência dos modelos foi a incapacidade de reproduzir a dependência do histórico de carregamentos, ou seja, da combinação de contrações. Desta maneira, uma nova proposta de modelo constitutivo foi desenvolvida e implementada para testes uniaxiais e em um código de elementos finitos para testes tridimensionais. A partir destas implementações, foram executados testes numéricos para diferentes sequências de contrações, avaliando a capacidade representativa do modelo proposto neste trabalho. O modelo proposto apresentou bons resultados para contrações isométricas e a combinação de contrações com diferentes níveis de alongamento, ou encurtamento, a uma mesma velocidade. A principal contribuição deste modelo é a capacidade de representar a resposta associada à fadiga muscular e o ganho, ou perda de força, observados experimentalmente. / This work present a striated skeletal muscle numerical characterization subjected to different contractions. From the physiological contraction mechanism of the skeletal muscle tissue the force produced in different kinds of contraction (isometric, concentric and eccentric), and their combination, are presented. The numerical representation of these response was investigated by different authors. From the diverse models found in literature, some were selected to be studied according to their representation capability and the implementation ease in numerical methods for large strains. The models available in literature were implemented and evaluated for different contractions combinations and a parameter identification for experimental results of two isometric contraction and a isometric-eccentric-isometric contraction. None of the proposed models presented satisfactory results. The main deficiency of these models were the incapability of reproducing the loading history dependence, in other words, the contraction combination. In this way, a new material model was proposed and implemented to uniaxial and tridimensional finite element method tests. From these, different contractions sequences evaluated the proposed model representation capabilities. The proposed model present good results to isometric contraction, as well as contraction combination with different stretch, or shortening, level in the same speed. The main contribution of this model is the capability of represent the the response associated to the muscle fatigue and the force gain or loss, experimentally observed.

Os processos de conversão de energia tendem a considerar cada vez mais restrições econômicas e ambientais, tornando-se necessário o entendimento da interação entre combustão e turbulência. Este trabalho tem como objetivo o desenvolvimento de soluções analíticas para a fração de mistura de uma chama difusiva, sob forma de um jato turbulento axissimétrico. Foi desenvolvida, também, uma metodologia analíticonumérica para a determinação das frações mássicas dos componentes, considerando uma reação de combustão de dois passos. Os resultados foram comparados com dados experimentais encontrados na literatura para uma chama de hidrogênio H2=N2 (50/50% em volume). De modo geral, os resultados obtidos foram satisfatórios frente aos dados experimentais, sendo a principal limitação o fato das expressões analíticas obtidas não serem capazes de representar o jato próximo à saída do bocal, sendo válidas a partir de, aproximadamente, x=d > 10, onde x é a coordenada ao longo do comprimento do jato. A principal vantagem do método empregado neste trabalho é a diminuição da complexidade do sistema de equações a ser resolvido numericamente. / The energy conversion processes tend to consider even more economical and environmental constraints, making it necessary to understand the interaction between combustion and turbulence. This study aims at the development of analytical solutions for the mixture fraction of a diffusive flame in the form of an axisymmetric turbulent jet. It was also considered an analytical-numerical approach for the determination of the mass fractions of the compounds, for a two-step reaction. The results were compared with data found in literature for a hydrogen flame H2=N2 (50/50 % by volume). Overall, the results were satisfactory when compared with the experimental data, however the principal limitation was the fact that the analytical expressions were not able to represent the jet near the nozzle exit, being the solution valid from, approximately, x=d > 10, where x is the coordinate along the jet length. The main advantage of the method employed in this work is the decrease in the complexity of the equations system to be solved numerically.

O estudo de parâmetros cinéticos para a queima de carvão se dá através de utilização de equipamentos laboratoriais como o forno de queda livre (DTF). O presente trabalho trata do desenvolvimento e análise de um modelo numérico unidimensional de um DTF. A modelagem é realizada a partir de dados da literatura, e as equações discretizadas que caracterizam os fenômenos do escoamento, da transferência de calor e da combustão do carvão foram programadas na linguagem FORTRAN 90. Os resultados são comparados com um modelo de referência e com dados experimentais, e apresentaram boa concordância. Um estudo de sensibilidade é realizado para entender o comportamento da queima de carvão frente a mudanças de alguns parâmetros de operação do DTF. A análise de sensibilidade mostrou a versatilidade do modelo e, consequentemente, do experimento. A partir da variação da concentração de O2, da temperatura de operação e das vazões de entrada pode-se obter uma gama de resultados. O uso do modelo aqui apresentado concomitante com resultados experimentais resulta em um estudo aprofundado dos parâmetros cinéticos de combustão de carvão. / The study of kinetic parameters for coal combustion occurs through use of laboratory equipment such as drop tube furnace (DTF). The present work shows the development and analysis of a one-dimensional numerical model of a DTF. The modeling is developed from literature works, and the discrete equations that characterize the flow, heat transfer and coal combustion are programmed in FORTRAN 90 language. The results are compared with a reference model and experimental data, and showed good agreement. A sensitivity study is performed to understand the behavior of coal combustion due to changes of some working parameters of the DTF. The sensitivity analysis showed the versatility of the model and thereby the experiment. From the variation of the O2 concentration, operating temperature and input flow rates a range of results can be obtained. Using the model presented here together with experimental results leads to a detailed study of the kinetic parameters of coal combustion.

Na área de simulação de processos, existe uma visível tendência da migração das ferramentas seqüenciais modulares, que hoje são as mais amplamente utilizadas, para as baseadas em equações. Uma das principais vantagens do paradigma baseado em equações ou simultâneo é que este se mostra eficiente na solução de problemas de simulação, otimização, estimação de parâmetros e reconciliação de dados, todos baseados em um mesmo conjunto de modelos, evitando retrabalho de modelagem. Porém, a tecnologia simultânea também apresenta algumas deficiências, onde destacam-se os problemas de robustez tanto na modelagem quanto na obtenção de resultados numéricos. Este trabalho tem como objetivo reunir e desenvolver técnicas que permitam reduzir estas deficiências. Para tanto, as técnicas conhecidas para depuração de sistemas de equações que representam problemas estacionários e dinâmicos foram estudadas em detalhe. Pôde-se observar que para o caso estático os métodos disponíveis para depuração de modelos, já se apresentam em um nível bem desenvolvido. Já para o caso dinâmico, onde há uma maior complexidade, as técnicas conhecidas encontram-se em um nível de desenvolvimento muito menor. Neste ponto encontram-se as principais contribuições deste trabalho. / In the field of process simulation the movement from the sequential modular tools, which are currently the most widely used, to the equation based approach is clear. One of the key advantages of the equation based or simulatneous approach is that using a single model one can solve simulation, optimization, parameter estimation, and optimization problems. This fact avoids modeling rework for each application. However, the simultaneous technology has problems regarding modeling and solving robustness. This work aims to group and develop methods capable of minimize these deficiencies. In order to achieve this goal, available debugging approaches for both steady-state and dynamic system of equations were studied in detail. For the steady-state case well stablished debugging techniques are known. For dynamic models, where the complexity is higher, the analysis and debugging methods are much less mature. This was the source for the major contributions of this work.

No presente trabalho são apresentadas soluções analíticas das equações de difusão de nêutrons bidimensionais com dois grupos de energia, a saber, nêutrons rápidos e térmicos em uma placa com propriedades homogêneas. Alem disso, são resolvidos detalhadamente os problemas onde a placa homogênea é substituída por duas e quatro regiões, tornando-os não-homogêneos. A partir da aplicação da transformada de Laplace e da Técnica da Transformada Integral Generalizada (GITT), respectivamente, é resolvida em uma forma analítica o problema de autovalor resultante para o fluxo de nêutrons. No problema heterogêneo são usados filtros para homogenizar as condições de contorno não-homogêneas. Esta é a condição para a aplicação da GITT. Os três problemas mencionados acima são resolvidos aplicando primeiramente a GITT, o qual reduz a dimensão da equação de difusão, seguida da aplicação da transformada de Laplace, o qual reduz a ordem da equação. Deste procedimento, resulta um sistema de equações algébricas dependente das constantes de integração. 0 sistema é resolvido usando a técnica da eliminação de Gauss. Os fluxos transformados pela GITT são recuperados invertendo-se analiticamente a transformada de Laplace usando a expansão de Heaviside, os quais ainda dependem das constantes de integração. A partir da aplicação das condições de contorno e de interface (para os problemas não-homogêneos) obtém-se um sistema de equações algébricas homogêneas, de onde é determinado o fator de multiplicação efetivo Keff pelo método da bissecção. As constantes de integração são determinadas fazendo use da potencia prescrita da placa. Assim, os fluxos de nêutrons transformados pela GITT ficam determinados e os fluxos de nêutrons rápidos e térmicos são recuperados através da formula da inversa da GITT, usando a expansão do potencial. Resultados são comparados com a solução do método de diferenças finitas. / In the present work we present analytical solutions of the bi-dimensional neutron diffusion equation with two energy groups, i.e. fast and thermal neutrons in a sheet with homogeneous properties. Further we solve the detailed problem where the homogeneous sheet is substituted by two and four regions, rendering the problem a non-homogeneous one. Upon application of the Laplace transform and Generalized Integral Transform Tecnique (GITT), respectively, we solve in an analytical fashion the resulting eigenvalue problem for the neutron flux. In the heterogeneous problem, we use filter functions in order to homogenize the non-homogeneous boundary conditions. This is a condition for the application of GITT. We solve the three problems mentioned above applying first GITT, which reduces the dimension of the diffusion equation followed by the Laplace transform, which reduces the order of the equation. This procedure yields a non-homogeneous algebraic system depending on integration constants. The system is solved using the elimination technique by Gauss. The transformed fluxes by GITT are recovered upon inverting analytically the Laplace transform using Heaviside's expansion which depend still on the integration constants. Upon application of the boundary and interface conditions (for the non-homogeneous problem) one obtains a system of homogeneous algebraic equations, where we determine the effective multiplication factor keff by the bisection method. The integration constants are determined making use of the predefined power of the sheet. Thus the neutron fluxes transformed by GITT are determined and the fast and thermal neutron flux are recovered by the inverse formula of GITT, using the potential expansion. Results are compared to the solution by the finite difference method.

Este trabalho apresenta metodologias para a construção e a sintonia de estimadores de estados não-lineares visando aplicações práticas. O funcionamento de um estimador de estados não-linear está calcado em quatro etapas básicas: (a) sintonia; (b) predição; (c) atualização da matriz de covariância de estados; (d) filtragem e suavização dos estados. As principais contribuições deste trabalho para cada uma destas etapas podem ser resumidas como segue: (a) Sintonia. A sintonia adequada da matriz de covariância do ruído de processos é fundamental na aplicação dos estimadores de estado com modelos sujeitos a incertezas paramétricas e estruturais. Sendo assim, foi proposto um novo algoritmo para a sintonia desta matriz que considera dois novos métodos para a determinação da matriz de covariância dos parâmetros. Este algoritmo melhorou significativamente a precisão da estimação dos estados na presença dessas incertezas, com potencialidade para ser usado na atualização de modelos em linha em práticas industriais. (b) Predição. Uma das etapas mais importantes para a aplicação do estimador de estados é a formulação dos modelos usados. Desta forma, foi mostrado como a formulação do modelo a ser usada em um estimador de estados pode impactar na observabilidade do sistema e na sintonia das matrizes de covariância. Também são apresentadas as principais recomendações para formular um bom modelo. (c) Atualização da matriz de covariância dos estados. A robustez numérica das matrizes de covariância dos estados usadas em estimadores de estados sem e com restrições é ilustrada através de dois exemplos da engenharia química que apresentam multiplicidade de soluções. Mostrou-se que a melhor forma de atualizar os estados consiste na resolução de um problema de otimização sujeito a restrições onde as estimativas fisicamente inviáveis dos estados são evitadas. Este também preserva a gaussianidade dos ruídos evitando que estes sejam mal distribuídos. (d) Filtragem e suavização dos estados. Entre as formulações estudadas, observou-se também que a melhor relação entre a acuracidade das estimativas e a viabilidade de aplicação prática é obtida com a formulação do filtro de Kalman estendido sujeita a restrições (denominada Constrained Extended Kalman Filter - CEKF), uma vez que esta demanda menor esforço computacional que a estimação de horizonte móvel, apresentando um desempenho comparável exceto no caso de estimativas ruins da condição inicial dos estados. Como uma solução alternativa eficiente para a estimação de horizonte móvel neste último caso, foi proposto um novo estimador baseado na inclusão de uma estratégia de suavização na formulação do CEKF, referenciado como CEKF & Smoother (CEKF&S). / This work presents approaches to building and tuning nonlinear state estimators aiming practical applications. The implementation of a nonlinear state estimator is supported by four basic steps: (a) tuning; (b) forecast; (c) state covariance matrix update; (d) states filtering and smoothing. The main contributions of this work for each one of these stages can be summarized as follows: (a) Tuning. An appropriate choice of the process-noise covariance matrix is crucial in applying state estimators with models subjected to parametric and structural uncertainties. Thus, a new process-noise covariance matrix tuning algorithm is presented in this work which incorporates two new methods for the parameter covariance matrix computation. The algorithm has improved significantly the state estimation accuracy when the presence of such uncertainties, with potential to be applied in on-line model update in industrial practice. (b) Forecast. One of the most important stages in applying state estimators is the used model formulation. In this way, it has been shown that the model formulation to be used in state estimator can impact on the system observability and noisecovariance matrices tuning. In this work it is also presented the main recommendations to formulate an appropriated model. (c) State covariance matrix update. The numerical robustness of the state covariance matrices used in unconstrained and constrained state estimators is illustrated by two chemical engineering examples tending to multiple solutions. It has been shown that the best technique to update the states consists in solving an optimization problem subjected to constraints, since it prevents from physically unfeasible states. It also preserves the noise gaussianity preventing from bad noise distribution. (d) States filtering and smoothing. Among the studied formulations, it was also noticed that the better relationship between performance and practical application is obtained with an extended Kalman filter formulation subjected to constraints (called Constrained Extended Kalman Filter - CEKF) because it requires small computational effort than MHE with comparable performance, except in case of poor guesses of the initial state. As an efficient solution for moving horizon estimation in the last case, it was proposed a new estimator based on the addition of a smoother strategy into the CEKF formulation, referred as CEKF & Smoother (CEKF&S).

Edifícios altos têm pouca rigidez lateral devido a suas dimensões. Os principais carregamentos laterais nestas estruturas são de sismos e vento. Portanto, suas propriedades dinâmicas e a correta estimativa das mesmas são fundamentais. Todavia, mesmo havendo muita pesquisa nesta área, a predição de frequências naturais de vibração livre de edifícios altos com precisão não é trivial. Métodos computacionais e fórmulas de normas ao redor do mundo divergem, ao passo que cidades brasileiras estão crescendo com a construção de maiores estruturas com frequências naturais estimadas abaixo de valores esperados. Neste trabalho, as incertezas acerca da estimativa de frequências naturais são investigadas através da comparação de resultados de diferentes metodologias. Primeiramente, frequências naturais estimadas ainda na fase de projeto de 8 prédios altos construídos ou em construção no Brasil são comparadas com formulações típicas baseadas em dados experimentais. Em um segundo momento, 3 edifícios altos têm suas primeiras frequências naturais estimadas através de: um método simplificado proposto neste trabalho; 2 programas de MEF diferentes; um modelo em MEF incluindo diagonais de compressão a fim de replicar o comportamento da alvenaria de fechamento e, finalmente, através da aplicação de um fator para calibrar a frequência natural obtida via programação em MEF. Tal fator é oriundo de uma formulação elaborada neste estudo para corrigir a frequência natural estimada sem a consideração do acréscimo de rigidez causado pela alvenaria de fechamento. Os resultados evidenciam que as baixas frequências naturais são semelhantes a medições de resposta sísmica. Também é relatada uma boa convergência para diferentes métodos computacionais enquanto os resultados das análises com alvenaria não estrutural revelam uma ótima convergência. Este estudo apresenta um melhor entendimento das razões para as baixas frequências naturais dos arranha-céus brasileiros enquanto expõe avanços na estimativa de frequências naturais com alta precisão. / Tall buildings have low lateral stiffness due to their dimensions and the main lateral forces on these structures are wind and seismic loads. Therefore, an accurate estimation of their dynamic properties is essential. Nevertheless, even with plenty of studies in this area, predicting natural frequencies of tall buildings is not trivial. Numeric computational methods and simplified procedures around the world diverge, whereas Brazilian cities are growing with the construction of higher buildings with lower estimated values for natural frequencies than the expected. In this work, the uncertainties are investigated through the comparison of different methods of estimating natural frequencies. Firstly, the natural frequencies estimated during project phase of 8 buildings, recently constructed or under construction in Brazil, are compared with formulas based on experimental measurements. Secondly, the natural frequencies of 3 buildings are estimated with different methods: two FEM models made with two computer programs commercially available, a simplified method proposed in this work, a FEM model with diagonal elements representing the masonry infill and a method based on applying a correction factor on the natural frequency obtained with another method. This last method is proposed here in order to consider the effect of infill masonry on the stiffness. The results evince that the low frequencies are similar to the ones measured on seismic response. The different computational methods have a good convergence while the analysis including the infill masonry has an outstanding convergence. This study reveals a better understanding of the reasons for low natural frequencies of the Brazilian tall buildings while presents improvements on estimating natural frequencies with high precision.