Spelling suggestions: "subject:"[een] FINITE ELEMENTS"" "subject:"[enn] FINITE ELEMENTS""
81 |
Estudo e modelagem de sistemas de detecção de danos em estruturas mecânicas baseados na impedância eletromecânica /Antunes, Rothschild Alencastro. January 2019 (has links)
Orientador: Jozué Vieira Filho / Resumo: As técnicas de detecção de danos baseadas na Impedância Eletromecânica (EMI) baseiam-se na capacidade dos materiais piezoeléctricos em atuar como sensores e atuadores e contribuem para o desenvolvimento de sistemas de Structural Health Monitoring (SHM). As técnicas clássicas baseadas na EMI utilizam um transdutor Pb-Lead Zirconate Titanite (PZT ) ligado à estrutura monitorada e medem a assinatura de impedância do PZT. No entanto, as técnicas baseadas na EMI dependem de diferentes fatores, como faixa de frequência, número de PZT, temperatura ambiente, tipo de estrutura, entre outros. Assim, para demonstrar a eficácia dos métodos baseados na EMI, faz-se necessário realizar experimentos práticos, o que não é uma tarefa trivial, considerando tais fatores. Portanto, neste trabalho são estudados e propostos procedimentos para criar modelos numéricos, usando elementos finitos (FE), de técnicas baseadas na EMI usando o software PZFlex®. Além disso, os modelos desenvolvidos são usados para propor uma técnica inovadora de compensação de temperatura em sistemas baseados na EMI. Foram modeladas, simuladas e analisadas algumas estruturas clássicas como placa-de-alumínio/PZT e tubo-de-aço/PZT. Os resultados das simulações foram comparados com os equivalentes obtidos com modelos experimentais reais e mostraram-se fortemente correlacionados, indicando que o modelo proposto pode ser uma ferramenta poderosa para o desenvolvimento de técnicas de SHM baseadas na EMI. Foram realizadas simulações ... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: The damage detection techniques based on the Electromechanical Impedance (EMI) rely on the ability of piezoelectric materials in acting as sensors and actuators and contribute to the development of Structural Health Monitoring (SHM) systems. The classical EMI-based techniques use a Pb-Lead Zirconate Titanite (PZT) transducer bonded to the monitored structure and measure the impedance signature of the PZT. However, the techniques based on EMI depend on different factors, such as frequency range, number of PZTs, environmental temperature, type of structure, among others. Thus, to demonstrate the effectiveness of EMI-based methods, it is necessary to carry out practical experiments, which is not a trivial task considering such factors. Therefore, in this work, it is studied and proposed procedures to create numerical models of techniques based on the EMI, using finite elements (FE) and the PZFlex® software. In addition, the developed models are used to propose an innovative temperature compensation technique for EMI-based systems. Some classical structures were modelled, simulated and analyzed, such as aluminum plate/PZT and steel pipe/PZT. The results of the simulations were compared with the equivalents obtained by experimental models and showed to be strongly correlated, indicating that the proposed model can be a powerful tool for the development of EMI-based SHM techniques. Simulations were performed to analyze the behavior of the signatures under the effect of temperature,... (Complete abstract click electronic access below) / Doutor
|
82 |
[en] STRUCTURAL INTEGRITY OF PLAIN DENTS IN PIPELINES / [pt] ANÁLISE DE INTEGRIDADE ESTRUTURAL DE DUTOS COM MOSSAS SIMPLESMAIRA AMANDA VARGAS AVILA 26 December 2007 (has links)
[pt] O presente trabalho analisa o impacto de uma mossa simples
transversal na
integridade estrutural de espécimes tubulares de paredes
finas com dimensões e
características semelhantes aos dutos usados para
transporte de óleo e gás. O
estudo da mossa simples envolveu as etapas de criação,
recuperação elástica e
avaliação da sua resistência estática e cíclica. A
abordagem do problema
fundamentou-se na análise de dados experimentais e na
modelagem numérica por
elementos finitos. Na análise por elementos finitos foram
consideradas: não
linearidade por contato, plasticidade, grandes
deslocamentos e grandes
deformações. Adicionalmente, foram feitas comparações
entre os resultados e os
procedimentos usados na indústria, de modo a avaliar
criticamente eventuais
diferenças e semelhanças entre estes. Como resultado, a
avaliação da ruptura
estática permitiu verificar a pouca influência de mossas,
com a geometria
estudada, na resistência à ruptura. Entretanto, os códigos
de projeto estudados
não aceitam mossas com profundidades e deformações menores
que as estudadas.
A avaliação de resistência cíclica permitiu verificar a
influência deste tipo de dano
na vida à fadiga. No final foi elaborada uma rotina
simples, baseada em
fundamentos teóricos, que permite quantificar a vida à
fadiga de mossas, com
características similares, a partir do seu fator de
concentração de tensão. / [en] This work analyzes the impact of transverse plain dents in
the structural
integrity of thin-thickness pipe specimens with similar
dimensions and properties
than typical pipelines. The dent analyzes included the
following topics: creation,
burst and fatigue evaluation and elastic recuperation. The
problem was based in
the study of experimental data and numerical modeling by
finite elements. The
numerical modeling considered: non linearity by contact,
plasticity, large
displacement and large strain. Finally, comparisons were
made between the
results and the typical industry procedures, in order to
identify eventual
differences and similarities. The burst evaluation
indicated the low influence of
the studied plain dents in the burst resistance.
Otherwise, the project standards do
not accept dents with the depth and strain levels like
those studied. The cyclic
analyze verified the influence of this type of damage in
the fatigue life. At the
end of this study, a simple procedure was developed based
in theoric
fundamentals to quantify the fatigue life of plain dents,
with similar characteristics
studied above, using their intensity stress factor.
|
83 |
[en] APPLICATION OF FINITE ELEMENT METHOD IN THE ANALYSIS OF FIELD VARIATIONS IN COAXIAL STRUCTURES DUE TO PERTURBATION IN THE BOUNDARY CONDITIONS / [pt] APLICAÇÃO DE MÉTODO DE ELEMENTOS FINITOS NA ANÁLISE DE VARIAÇÕES DE CAMPO EM ESTRUTURAS COAXIAIS DEVIDO A PERTURBAÇÕES NAS CONDIÇÕES DE CONTORNOJANAINA CUNHA E SILVA ARTEAGA 20 May 2008 (has links)
[pt] Neste trabalho, o Método dos Elementos Finitos (MEF) é
aplicado na análise do campo eletromagnético em dutos
metálicos, utilizados no transporte de óleo e gás.
Instrumentos de inspeção podem utilizar esse método para
verificar a existência de anomalias do tipo amassamentos e
corrosões nos dutos, contribuindo para evitar vazamentos
que podem causar acidentes ecológicos. Atualmente, a
técnica de fuga de campo magnético (MFL) é a mais utilizada
para precisar a profundidade dos defeitos encontrados nos
dutos devido ao processo localizado de corrosão, onde a
espessura se reduz progressivamente até o inicio de um
eventual vazamento. Os picos dos sinais de fluxo magnético
defletido indicam a perda de massa metálica localizada. A
profundidade dos defeitos é correlacionada com a
amplitude dos picos. Neste trabalho, o algoritmo baseado em
Método de Elementos Finitos é utilizado para avaliar o
desempenho do MEF na localização da posição de defeitos na
parede do duto, através dos picos que aparecem nos
gráficos de campo magnético. Essa análise será feita
utilizando-se dois tipos de alimentação na entrada do duto:
(I) apenas propagação do modo fundamental TEM e (II)
diferença de potencial entre os cilindros externo e interno
do duto. Diversos tamanhos de deformação em posições
diferentes do duto serão analisados para se determinar para
quais tipos de problema o algoritmo é eficiente. / [en] In this work, the Finite Element Method (FEM) is applied in
the analysis of electromagnetic field generated by
instruments employed for inspection of pipelines for
natural gas-transmission. Inspection instruments can use
this method to evaluate anomalies as corrosion and kneading
in the superficies of the pipelines. This can avoid leaking
that can cause ecological accidents. The Magnetic Flux
Leakage (MFL) is the oldest and most commonly used in-line
inspection method for finding metal-loss regions in gas-
transmission pipelines. MFL inspections are typically used
to detect, locate, and characterize metal-loss and other
anomalies in natural gas-transmission pipelines. The
amplitude or magnitude of an MFL signal is strongly related
to defect depth. Alternatively, in this work, an instrument
that generates a TEM wave is explored for inspection of
pipelines. To evaluate the field distribution inside the
pipe, an algorithm based on Finite Element Method is used
to detect and locate an anomaly in the pipeline. It
is used two kinds of source in the duct´s port: (I) just
propagation in the FEM basic mode and (II) different
potential between the internal and external
cylindrical duct. Several sizes of anomalies will be
analyzed to identify for which kind of problems the
algorithm is useful.
|
84 |
[en] NUMERICAL AND EXPERIMENTAL INVESTIGATION OF THE MECHANICS OF ANEURYSMS IN ISOTROPIC RUBBER TUBES / [pt] INVESTIGAÇÕES NUMÉRICAS E EXPERIMENTAIS DA MECÂNICA DOS ANEURISMAS EM TUBOS ISOTRÓPICOS DE BORRACHALUCAS BOABAID IBRAHIM 26 January 2007 (has links)
[pt] Esta tese tem por objetivo investigar numérica e
experimentalmente a
mecânica da formação dos aneurismas na aorta. A parte
experimental foi realizada
no Laboratório de Membranas e Biomembranas utilizando-se
tubos cilíndricos de
látex sob pressão hidrostática e tubos de silicone com
geometria aproximada da
aorta. Foi investigada a pressão necessária à formação dos
aneurismas e o
comportamento do material ensaiado. A parte numérica foi
realizada por meio do
método dos elementos finitos através do programa ABAQUS.
Na análise
numérica foi validada a análise experimental e realizados
alguns estudos
paramétricos. / [en] The objective of this work was to investigate numerical
and experimentally
the mechanics of aortic aneurisms. The experimental part
was done in the
Laboratory of Membranes and Biomenbranes using latex
cylindrical tubes under
hydrostatic pressure and rubber tubes with the approximate
aortic geometry. The
required pressure to formation of aneurisms was
investigated as well the behavior
of the material during the experiments. The numerical part
was done with finite
element method with the ABAQUS program. In the numerical
analyses the
experimental analyses was validated and some parametric
studies was done.
|
85 |
Implementação do método das características na modelagem de problemas de convecção natural em cavidades cilíndricas /Jordam, Alice. January 2010 (has links)
Orientador: Vicente Luiz Scalon / Banca: Helio Aparecido Navarro / Banca: Sergio Rodrigues Fontes / Resumo: A fluidodinâmica computacinal (CFD) tem sido utilizada, estudadda e implementada ao longa das duas últimas décadas na solução dos mais diversos problemas de engenharia. O princípio básico desta ciência é a aplicação de métodos numéricos em problemas que envolvam mecânica dos fluidos. Nesse contexto, este trabalho utiliza essa técnica para analisar o comportamento de um fluido incompreensível, que se encontra numa cavidade cilíndrica fechada onde as faces inferior e superior são adiabáticas e as superfícies laterais se encontram em diferentes temperaturas. Os perfis de velocidade e temperatura resultantes - ocasionados pela convecção natural - serão avaliados em todo o domínio do problema. Existe uma série de técnicas para a solução de problemas envolvendo escoamentos, sendo as mais comuns as que se utilizam do "Esquema de Passo Fracionado" proposto por Chorin no final da década de 60. Dentre as diversas soluções que se utilizam desta técnica, este trabalho optou pelo uso do método das características e do algoritmo CBS de solução proposto por Zienkiewicz e Codina (1995). Para a implementação do algoritmo de solução do problema proposto foi realizada uma discretização geral através do método dos elementos finitos usando-se de uma malha formada por elementos bilineares. A solução foi obtida a partir de um ambiente matemático adequado, a GNU-Octave (2008). Os resultados foram analisados para diferentes razões de curvatura, números de Rayleigh e métodos de solução, sendo plotados para as suas diversas variáveis buscando descrever o comportamento do fenônemo / Abstract: The Computational Fluid Dynamics (CFD) has been used, studied and performed through the last two decades to solve the series of problems in Engineering. The most basic aim of this science is the appliance of numerical methods in cases that envolve fluid mechanics. In this context, this work uses this technic to analyze the behaviour of an incompressible fluid, which is found in a closed cylindrical cavity, a place where the inferior and superior surfaces are adiabatic and the lateral surfaces are shown in different temperatures. The resultant profiles of speed and temperature - induced by the free convection - are going to be appraised in all the dominion of the problem. There is a set of technics to solve the problems which involve the drainage, but the most usual are those which use the techic "Fractional Step Method" offered by Chorin in the final of 60s. Among the several solutions that are solved through this technic, this research used the characteristics method and of the CBS algorthm, offered by Zienkiewicz e Codina (1995). For the implementation of the algorithm, it was realized a general discretization through the finite elements method, making use of a loop formed by bilinear elements. The resolution was obtained from an adequated mathematical ambient, the GNU-Octave (2008). The results were analysed for different curvature ratios, Rayleigh numbers and methods of solution, being plotted for its different variables searching to describe the behavior of the phenomenon / Mestre
|
86 |
Estudo de problemas de escavação através da combinação elementos de contorno e elementos finitos / Study of excavation problems with a boundary element and finite element methods combinationKomatsu, José Sergio 15 September 1995 (has links)
Estuda-se uma combinação do método dos elementos finitos (MEF) com o método dos elementos de contorno (MEC) no acoplamento de uma estrutura reticulada em um domínio bidimensional. Para o caso em análise, os elementos uniaxiais são tratados através do MEF, enquanto que o MEC é utilizado na modelagem do meio contínuo que pode ser homogêneo ou não-homogêneo. Em problemas geomecânicos, é possível simular a sequência de escavação com as modificações estruturais necessárias. As equações do domínio bidimensional, com a influência das estruturas reticuladas, são agrupadas em um sistema que é resolvido pelo algoritmo proposto por Crotty. Utilizando-se o método dos elementos de contorno, a plasticidade do meio contínuo é analisada com um procedimento incremental e iterativo baseado no processo das tensões iniciais. / A combination of finite element method boundary element method (MEC) is studied (MEF) with the for the frame structure and two-dimensional domain. In the analyzed case, the uniaxial elements are treated by MEF, while the MEC be used to model the continuum media which can is homogeneous or not. In geomechanical problems, it is possible to simulate a sequence excavation with the necessary structural modifications. The equations of the two-dimensional domain with the influence of the frame structure, are assembled in the system which is solved by algorithm proposed by Crotty. Using the boundary element method, it is analyzed the plastic behavior of a continuum media by an incremental and iterative procedure, which is based in the initial stress process.
|
87 |
Avaliação das tensões geradas nos componentes protéticos de próteses obturadoras maxilares classe I, II e IV de Aramany por meio de análise de elementos finitos / Evaluation of stress generated in maxillar obturators prosthesis of class I, II and IV of Aramany in Finite Element AnalysisMiyashita, Elcio Ricardo 20 March 2013 (has links)
O tratamento cirúrgico do câncer bucal pode resultar em graves sequelas das estruturas anatômicas envolvidas na fisiologia da cavidade bucal. A extensão e a localização da ressecção cirúrgica, assim como a condição dental, são determinantes do planejamento da prótese obturadora. Este trabalho emprega a Análise de Elementos Finitos para o estudo do comportamento biomecânico de próteses obturadoras maxilares Classe I, II e IV de Aramany, em um modelo digital desenvolvido a partir de uma TC de um indivíduo adulto. O modelo tridimensional desenvolvido no programa Rhinoceros® versão 4.0 foi utilizado para a geração da malha de elementos finitos no programa Ansys Workbench 14.0®. Uma força de 120 N foi aplicada nas plataformas oclusal e incisal, correspondentes aos dentes posteriores e anteriores das próteses obturadoras, para observação das tensões geradas sobre as próteses obturadoras. Os resultados são apresentados sob a forma de análise qualitativa, indicada pelo estudo de Tensão Máxima Principal para as estruturas acrílicas da prótese e estudo de Von Mises para as estruturas metálicas da prótese, e análise quantitativa, representada por valores em Mega Pascal (MPa). Tanto nos carregamentos posteriores como nos carregamentos anteriores, foram observadas dissipações de cargas significativas localizadas na placa palatina próxima a ressecção e na junção do metal com o acrílico, sendo os maiores valores alcançados na placa de retenção metálica para retenção do acrílico localizada posteriormente na condição de carregamento posterior em ordem crescente na Classe II de Aramany com 149,17 MPa, Classe I de Aramany com 344,08 MPa e Classe IV de Aramany alcançando valor de 390,25 MPa. No acrílico da prótese obturadora, no caso de Classe I de Aramany em carregamento posterior, foi observada uma concentração de tensão principal na região correspondente ao limite lateral e anterior da junção da placa de retenção e do acrílico no valor de 19,796 MPa de tração. Na Classe IV de Aramany em carregamento posterior, foi observada uma concentração de tensão na região posterior da junção da placa de retenção e do acrílico no valor de 18,33 MPa. Na Classe II de Aramany, em carregamento posterior foi observada uma concentração de tensão na região anterior próxima ao limite da ressecção anterior no valor de 25,45 MPa de tração. Concluiu-se que os valores quantitativos alcançados neste estudo parecem não apresentar riscos de fratura sobre a armação metálica de liga de CrCo, porem podem apresentar riscos de falha na retenção metal/acrílico da prótese obturadora maxilar. A configuração de placa de recobrimento palatino proporciona uma melhor distribuição das tensões ao longo de sua superfície. / The surgical treatment of oral cancer may result in severe sequel of the anatomical structures involved in oral physiology. The extension and location of surgical resection, as well as the dental condition, determine the prosthetic planning of the obturator prostheses. This study employs the finite elements analysis to evaluate the biomechanical behavior of Aramany Class I,II andIV obturator prostheses, in a digital model developed from a computerized tomography of an adult individual. A tridimensional model was constructed applying the software Rhinoceros® version 4.0 and then it was used to develop the finite element mesh in the software Ansys Workbench 14.0 ®. A 120 N load was applied to the occlusal and incisal surfaces correspondent to the prosthetic teeth. The qualitative analysis was based on Principal Maximum stress and the quantitative analysis was expressed in Mega Pascal values. Results demonstrated that in the posterior loadings as in the anterior loadings were observed charge dissipation plate located in significant palate resection and near the junction of the metal with the acrylic and the highest values obtained in the metal retaining plate for retaining acrylic located further provided that subsequently loaded in ascending order in Class II Aramany with 149.17 MPa, Aramany with Class I and Class IV 344.08 MPa Aramany reaching value of 390.25 MPa. In acrylic prosthesis obturator in the case of Class I Aramany in subsequent loading, was observed a stress concentration in the region corresponding to the limit lateral and anterior junction of the retaining plate and the value of acrylic 19.796 MPa. Class IV for loading Aramany later, it was observed a stress concentration in the posterior region of the junction of the retaining plate of acrylic and the value of 18.336 MPa. In Class II Aramany in load was later observed a stress concentration in the anterior region near the boundary of the anterior resection worth 25.455 MPa. It was concluded that the quantitative values obtained in this study do not appear to present a risk of fracture of the metal frame CrCo alloy, but may present risks of failure in the metal retention / acrylic obturator prostesis . The configuration of palatal plate coating provides a better stress distribution along its surface.
|
88 |
Lógica difusa para controle não convencional de uma viga inteligente. / Fuzzy logic for non-conventional control of an intelligent beam.Nagamine, Renato Kazuki 02 August 2001 (has links)
Os avanços da indústria aeronáutica têm sido garantidos pelo emprego de tecnologias inovadoras. Para controle aeroelástico o conceito de estrutura inteligente tem ganho cada vez mais espaço. Uma estrutura inteligente é um sistema de controle estrutural onde elementos estruturais, sensores, atuadores, leis de controle, eletrônica associada e processamento de sinais estão altamente integradas garantindo aumento de desempenho. Os desenvolvimentos nesta área têm sido muito animadores e envolve uma série de disciplinas. Nesse contexto, o presente trabalho tem como meta estudar uma estrutura inteligente onde a lei de controle é representada através da lógica difusa. Esse método não convencional de controle tem proporcionado avanços no trato de sistemas complexos, não lineares e com parâmetros imprecisos e ambíguos. Um modelo em elementos finitos de uma viga inteligente com atuadores piezelétricos incorporados é desenvolvido. O modelo baseia-se nas hipóteses de viga Euler-Bernoulli e no princípio variacional eletromecânico. O modelo em elementos finitos é validado para garantir o uso no projeto do controlador não convencional. Estratégias de controle não convencional baseadas em dois tipos de metodologia difusa para controle, isto é, modelo de Mamdani e de Takagi-Sugeno-Kang, são estudas. O controlador difuso é aplicado para reduzir a resposta vibratória da viga inteligente quando submetida a distúrbios mecânicos externos. Um estudo comparativo das duas metodologias de controlador difuso é realizado e discutido. Os resultados satisfatórios alcançados mostram que o uso de controladores difusos para alívio de vibrações em vigas com atuadores piezelétricos é apropriado. Também se observou que o modelo Takagi-Sugeno-Kang é o que mais se ajustou as necessidades e requerimentos do problema. / The advances in aeronautical industry have been ensured by the application of novel technologies. For aeroelastic control the intelligent structure concept has become increasingly important. An intelligent structure is a structural control system in which structural elements, sensors, actuators, control laws, associated electronics, and signal processing are highly integrated providing better performance. The developments in this area have been encouraging and multi-disciplinar. In this context, this work aims studying control laws for intelligent structures represented through the fuzzy logic. This non-conventional method for control has provided advances in treating complex, non-linear, imprecise, and ambiguous systems.A finite element model of an intelligent beam with piezoelectric actuators is developed. The model is based in the assumptions of an Euler-Bernoulli beam and the electromechanic variational principle. The finite element model is validated to ensure its use in non-conventional control design. Non-conventional control strategies based on the two fuzzy control methodologies, that are, the Mamdani and Takagi-Sugeno-Kang models, are investigated. The fuzzy control is applied to reduce the vibratory response of the intelligent beam due to external mechanical disturbances. A comparison between the two fuzzy control methodologies is shown and discussed. The satisfactory results achieved by this work have shown that the application of fuzzy control for the alleviation of vibrations in beams with piezoelectric actuators is appropriate. It is also observed that the Takagi-Sugeno-Kang fuzzy model has presented a better performance when compared with the Mamdani one.
|
89 |
Eléments finis adaptatifs pour l'équation des ondes instationnaire / Adaptive finite elements for the time-dependent wave equationGorynina, Olga 22 February 2018 (has links)
La thèse porte sur l’analyse d’erreur a posteriori pour la résolution numérique de l’équation linéaire des ondes , discrétisée en temps par le schéma de Newmark et en espace par la méthode des éléments finis. Nous adoptons un choix particulier de paramètres pour le schéma de Newmark, notamment β = 1/4, γ = 1/2, qui assure que la méthode est conservative en énergie et d’ordre deux en temps. L’estimation d’erreur a posteriori, d’un ordre optimal en temps et en espace, est élaborée à partir de la discrétisation complète. L’erreur est mesurée dans une norme qui découle naturellement de la physique: H1 en espace et Linf en temps. Nous proposons d’abord un estimateur dit «à 3 points» qui fait intervenir la solution discrète en 3 points successifs du temps à chaque pas de temps. Cet estimateur fait appel à une approximation du Laplacien de la solution discrète qui doit être calculée à chaque pas de temps, en résolvant un problème auxiliaire d'éléments finis. Nous proposons ensuite un estimateur d’erreur alternatif qui permet d’éviter ces calculs supplémentaires: l’estimateur dit «à 5 points» puisqu’il met en jeu le schéma des différences finies d’ordre 4, qui fait intervenir la solution discrète en 5 points successifs du temps à chaque pas de temps. Nous démontrons que nos estimateurs en temps sont d’ordre optimal pour des solutions suffisamment lisses, sur des maillages quasi-uniformes en espace et uniformes en temps, en supposant que les conditions initiales soient discrétisées à l’aide de la projection elliptique. La trouvaille la plus intéressante de cette analyse est le rôle capitale de cette discrétisation : des discrétisations standards pour les conditions initiales, telles que l’interpolation nodale, peuvent être néfastes pour les estimateurs d’erreur en détruisant leur ordre de convergence, bien qu’elles fournissent des solutions numériques tout à fait acceptables. Des expériences numériques prouvent que nos estimateurs d’erreur sont d’ordre optimal en temps comme en espace, même dans les situations non couvertes par la théorie. En outre, notre analyse a posteriori s’étend au schéma de Newmark d’ordre deux plus général (γ = 1/2). Nous présentons des comparaisons numériques entre notre estimateur à 3 points généralisé et l’estimateur sur des grilles décalées, proposé par Georgoulis et al. Finalement, nous implémentons un algorithme adaptatif en temps et en espace basé sur notre estimateur d’erreur a posteriori à 3 points. Nous concluons par des expériences numériques qui montrent l’efficacité de l’algorithme adaptatif et révèlent l’importance de l’interpolation appropriée de la solution numérique d’un maillage à un autre, surtout vis à vis de l’optimalité de l’estimation d’erreur en temps. / This thesis focuses on the a posteriori error analysis for the linear second-order wave equation discretized by the second order Newmark scheme in time and the finite element method in space. We adopt the particular choice for the parameters in the Newmark scheme, namely β = 1/4, γ = 1/2, since it provides a conservative method with respect to the energy norm. We derive a posteriori error estimates of optimal order in time and space for the fully discrete wave equation. The error is measured in a physically natural norm: H1 in space, Linf in time. Numerical experiments demonstrate that our error estimators are of optimal order in space and time. The resulting estimator in time is referred to as the 3-point estimator since it contains the discrete solution at 3 points in time. The 3-point time error estimator contains the Laplacian of the discrete solution which should be computed via auxiliary finite element problems at each time step. We propose an alternative time error estimator that avoids these additional computations. The resulting estimator is referred to as the 5-point estimator since it contains the fourth order finite differences in time and thus involves the discrete solution at 5 points in time at each time step. We prove that our time estimators are of optimal order at least on sufficiently smooth solutions, quasi-uniform meshes in space and uniform meshes in time. The most interesting finding of this analysis is the crucial importance of the way in which the initial conditions are discretized: a straightforward discretization, such as the nodal interpolation, may ruin the error estimators while providing quite acceptable numerical solution. We also extend the a posteriori error analysis to the general second order Newmark scheme (γ = 1/2) and present numerical comparasion between the general 3-point time error estimator and the staggered grid error estimator proposed by Georgoulis et al. In addition, using obtained a posteriori error bounds, we implement an efficient adaptive algorithm in space and time. We conclude with numerical experiments that show that the manner of interpolation of the numerical solution from one mesh to another plays an important role for optimal behavior of the time error estimator and thus of the whole adaptive algorithm.
|
90 |
Avaliação das tensões transmitidas às estruturas de suporte, por meio da associação dos métodos fotoelástico e de elementos finitos, decorrentes das forças aplicadas sobre overdentures retidas em implantes com sistema de encaixe barr / Appraisal, by association of photoelastic and finite elements methods, of load transmission on support structures with overdentures retained in implants with bar-clip systemPigozzo, Mônica Nogueira 17 August 2010 (has links)
O objetivo desta pesquisa foi avaliar in vitro as tensões transmitidas às estruturas de suporte, por meio da associação dos métodos fotoelásticos e elementos finitos, decorrentes da aplicação de forças sobre overdentures retidas em implantes com sistema de encaixe barraclipe, particularmente quando duas posições dos implantes são simuladas. Assim, dois modelos fotoelásticos de mandíbulas foram confeccionados, com dois implantes cada um, posicionados na região interforaminal a 22 mm de distância: modelo 1, denominado AFIP, modelo fotoelástico com implantes paralelos e orientados verticalmente; e modelo 2, denominado AFII, modelo fotoelástico com implantes angulados 10º em relação à linha média da mandíbula. Sobre os implantes, em cada modelo, foi posicionada uma prótese tipo overdenture, com sistema de retenção barra-clipe. Para simulação da mucosa, foi interposta entre a base da prótese e o rebordo residual uma camada de 2 mm de silicone. Este mesmo modelo foi representado para uma análise tridimensional (3D) de elementos finitos. As estruturas foram consideradas homogêneas, isotrópicas e lineares, além de serem atribuídas as correspondentes propriedades mecânicas para cada uma delas. Foram simuladas quatro mandíbulas de elementos finitos, sendo modelo 1, denominado EFRIP, com implantes paralelos; modelo 2, denominado EFRII, com implantes angulados em 10º, ambos modelos (1 e 2) simulando o modelo mandibular de resina fotoelástica; modelo 3, denominado EFOIP, com implantes paralelos; e modelo 4, denominado EFOII, com implantes angulados em 10º, ambos os modelos (3 e 4), simulando os ossos cortical e trabecular. Foram aplicadas cargas de 0,5; 1,0; 1,5; 2,0 e 3 bar no método fotoelástico e 0,1 MPa no método de elementos finitos e em seguida as imagens foram fotografadas e analisadas. Os resultados mostraram que: não houve semelhança entre as áreas de tensão dos modelos de resina fotoelástica, quando avaliada a angulação dos implantes (modelo 1 AFIP - e modelo 2 AFII). Houve semelhança entre as áreas de tensão nos modelos de elementos finitos com mandíbula em resina (modelo 1 EFRIP - e modelo 2 EFRII), tendo ocorrido o mesmo entre os modelos de elementos finitos com simulação de osso (modelo 3 - EFOIP - e modelo 4 EFOII) quando avaliada a inclinação dos implantes. Quando comparados os diferentes modelos com a mesma angulação de implante, observou-se que houve semelhança na distribuição das áreas de tensão entre os modelos 1 (AFIP e EFRIP), com uma alta concentração de tensão nos ápices dos implantes. Porém, não houve relação entre os modelos de implantes inclinados, modelos 2 (AFII e EFRII). Da mesma forma, não houve relação entre os modelos de análise fotoelástica e elementos finitos com simulação de mandíbula em resina, independentemente da angulação dos implantes, comparado os modelos de elementos finitos com simulação de osso. Frente às limitações do presente trabalho, concluiu-se que a associação dos métodos de análise fotoelástica e de elementos finitos é de grande valia para a obtenção de informação em relação à biomecânica referente a esse tipo de prótese e sistema de retenção. Além disso, houve semelhança nas áreas de concentração de tensão, quando comparado os modelos de análise fotoelástica e de elementos finitos, quando as propriedades mecânicas inseridas no modelo de elementos finitos foram compatíveis com o modelo fotoelástico, pois, quando o modelo de elementos finitos simulou uma mandíbula com diferenciação entre osso cortical e medular, as áreas de concentração de tensão sofreram alterações. / The objective of this study was to appraise in vitro the load transmission in bar-clip retention system for overdenture when two implant positions are simulated, using for this purpose the association of photoelastic analysis and finite elements techniques. Therefore, two photoelastic mandible were manufactured with two implants each one, positioned in interforaminal region within 22 mm of distance: (i) model 1, named AFIP, photoelastic model with parallel and vertically oriented implants; and (ii) model 2, named AFII, photoelastic model with 10º angled implants in relation to mandible midline. A bar-clip retention system and an overdenture were fixed over both implants. To simulate oral mucosa were added a 2 mm silicon layer between overdenture prosthesis and photoelastic mandible. This same model was submitted to finite elements analysis in 3D, being attributed the correspondent mechanic properties to each structure analyzed and also being considered homogenous, isotropic and linear characteristics. Four finite elements mandible were simulated, (i) one with parallel implants and (ii) other with 10º angled implants, simulating photoelastic resin in mandible (models 1 and 2, named EFRIP and EFRII respectively); and a another couple of mandibles also with (iii) parallel and (iv) 10º angled implants, simulating cortical and cancellous bones in mandible (models 3 and 4, named EFOIP and EFOII respectively). 0.5; 1.0; 1.5; 2.0 and 3.0 bars loads were applied and correspondent images were analyzed. Results evidenced that there is no similarity in loaded areas between photoelastic resin models, when the implant angle is appraised in models 1 and 2 (AFIP and AFII). Notwithstanding, there are similarity between loaded areas in finite elements model with resin mandible, models 1 and 2 (EFRIP and EFRII), and the same have occurred with finite elements models with bone simulation, models 3 and 4 (EFOIP and EFOII), when appraised implant inclination. When different models are compared with the same implant angulations, was noted that there was similarity in the distribution of loaded areas between models 1 (AFIP and EFRIP), with high concentration of tension forces at the implant apex. Nevertheless, there was no relation between inclined models 2 (AFII and EFRII). At the same direction, there was no relation between photoelastic analysis and finite elements with resin model, with no influence of implant angulations, with finite element models with bone simulation. Considering the limitations of this study, it was concluded that association of photoelastic analysis and finite elements has a high potential value for obtaining information concerning the biomechanics of this kind of prosthesis and retention system; furthermore, there was similarity in areas where loads were applied, when compared photoelastic and finite elements models, provided that finite elements mechanic properties are compatible with photoelastic model, since concentration areas of tension forces suffer changes when finite elements model simulates a mandible with cortical and cancellous bones differentiation.
|
Page generated in 0.0298 seconds