Formulação de elemento finito posicional para modelagem numérica de pórticos planos constituídos por compósitos laminados: uma abordagem não linear geométrica baseada na teoria Layerwise
Data
2017-11-15
Autores
Título da Revista
ISSN da Revista
Título de Volume
Editor
Biblioteca Digital de Teses e Dissertações da USP
Universidade de São Paulo
Escola de Engenharia de São Carlos
Universidade de São Paulo
Escola de Engenharia de São Carlos
Resumo
Descrição
A análise de compósitos laminados apresenta grandes desafios, pois, diferentemente dos materiais isotrópicos homogêneos, os compósitos laminados são constituídos de materiais heterogêneos e anisotrópicos. Além disso, as distribuições de tensões interlaminares obtidas com as formulações convencionais são descontínuas e imprecisas. Sua melhoria, portanto, é imprescindível para buscar e modelar critérios de falha relacionados às estruturas formadas por compósitos laminados. Diante disso, este trabalho se concentrou no desenvolvimento e implementação computacional de um elemento finito posicional de pórtico plano laminado cuja cinemática é descrita ao longo da espessura do laminado de acordo com a teoria Layerwise. A formulação do elemento considera a não linearidade geométrica, originada pela ocorrência de grandes deslocamentos e rotações, e admite deformações moderadas, em função da lei constitutiva de Saint-Venant-Kirchhoff. O desenvolvimento deste trabalho se iniciou com uma preparação teórica sobre mecânica dos sólidos deformáveis e métodos numéricos para que fossem adquiridos os subsídios teóricos necessários ao desenvolvimento de códigos computacionais, à interpretação dos resultados e à tomada de decisões quando das análises numéricas. A formulação desenvolvida é Lagrangiana total com emprego do método dos elementos finitos baseado em posições. Inicialmente o elemento finito posicional de pórtico plano homogêneo é proposto, uma vez que sua cinemática possibilita uma expansão natural para o caso laminado. Os graus de liberdade são compostos por posições nodais e por vetores generalizados que representam o giro e a variação na altura da seção transversal. A eficiência do elemento é constatada através de análises realizadas em problemas de pórtico sujeitos a grandes deslocamentos e rotações. Os resultados obtidos apresentaram excelente concordância com soluções numéricas e analíticas disponíveis na literatura. Uma expansão natural da cinemática é empregada na formulação do elemento laminado. Os graus de liberdade do elemento são as posições nodais e as componentes de vetores generalizados associados às seções transversais de cada lâmina. Dessa forma, as lâminas têm liberdade para variação de espessura e giro independente das demais, mas com as posições compatibilizadas nas interfaces. Os resultados de análises numéricas realizadas em vários exemplos demonstram a eficiência da formulação proposta, pois as distribuições de deslocamentos e tensões ao longo da espessura do laminado apresentaram excelente concordância com as obtidas a partir de análises numéricas utilizando um elemento finito bidimensional em uma discretização bastante refinada. Os exemplos analisados contemplam problemas com seção laminada fina ou espessa.
The analysis of laminated composites presents challenges because, unlike homogeneous isotropic materials, the laminated composites are made up of heterogeneous and anisotropic materials. Moreover, the distribution of interlaminar stresses obtained with conventional formulations are discontinuous and inaccurate. His improvement is therefore essential to check and modeling failure criteria related to structures formed by laminates. Thus, this work focused on developing and computational implementation of a positional finite element of laminated plane frame whose kinematics is described throughout the thickness of the laminate according to Layerwise theory. The formulation element considers the geometric nonlinearity, caused by the occurrence of large displacements and rotations, and admits moderate deformation, in the constitutive law function of Saint-Venant-Kirchhoff. The development of this work began with a theoretical preparation on mechanics of deformable solids and numerical methods for the acquired of the theoretical support needed for the development of computational codes, interpretation of results and decision-making when of the numerical analyzes. The developed formulation is total Lagrangian with use of the finite element method based on positions. Initially the positional finite element of homogeneous plane frame is proposed, since their kinematic enables a natural expansion for the laminate case. The degrees of freedom are composed of nodal positions and generalized vectors representing the spin and the variation in the height of the cross section. The efficiency of the element is verified through analyzes performed in frame problems subject to large displacements and rotations. The results showed excellent agreement with numerical and analytical solutions available in the literature. A natural expansion of the kinematics is used in the formulation of the laminate element. The degrees of freedom of the element are the nodal positions and components of the generalized vectors associated to cross-sections of each lamina. Thus, the laminas are free for the thickness variation and for independent spin, but with the positions matched in the interfaces. The results of numerical analysis performed in various examples show the effectiveness of the proposed formulation, since the distributions of displacements and stresses through the thickness of the laminate agreed well with those obtained from numerical analysis using a discretization with two-dimensional finite elements in a very refined. The examples discussed include problems with thin or thick laminated section.
The analysis of laminated composites presents challenges because, unlike homogeneous isotropic materials, the laminated composites are made up of heterogeneous and anisotropic materials. Moreover, the distribution of interlaminar stresses obtained with conventional formulations are discontinuous and inaccurate. His improvement is therefore essential to check and modeling failure criteria related to structures formed by laminates. Thus, this work focused on developing and computational implementation of a positional finite element of laminated plane frame whose kinematics is described throughout the thickness of the laminate according to Layerwise theory. The formulation element considers the geometric nonlinearity, caused by the occurrence of large displacements and rotations, and admits moderate deformation, in the constitutive law function of Saint-Venant-Kirchhoff. The development of this work began with a theoretical preparation on mechanics of deformable solids and numerical methods for the acquired of the theoretical support needed for the development of computational codes, interpretation of results and decision-making when of the numerical analyzes. The developed formulation is total Lagrangian with use of the finite element method based on positions. Initially the positional finite element of homogeneous plane frame is proposed, since their kinematic enables a natural expansion for the laminate case. The degrees of freedom are composed of nodal positions and generalized vectors representing the spin and the variation in the height of the cross section. The efficiency of the element is verified through analyzes performed in frame problems subject to large displacements and rotations. The results showed excellent agreement with numerical and analytical solutions available in the literature. A natural expansion of the kinematics is used in the formulation of the laminate element. The degrees of freedom of the element are the nodal positions and components of the generalized vectors associated to cross-sections of each lamina. Thus, the laminas are free for the thickness variation and for independent spin, but with the positions matched in the interfaces. The results of numerical analysis performed in various examples show the effectiveness of the proposed formulation, since the distributions of displacements and stresses through the thickness of the laminate agreed well with those obtained from numerical analysis using a discretization with two-dimensional finite elements in a very refined. The examples discussed include problems with thin or thick laminated section.
Palavras-chave
Teoria Layerwise, Pórtico plano, Não linearidade geométrica, Elemento finito posicional, Compósito laminado, Geometric nonlinearity, Layerwise theory, Plane frame, Positional finite element, Composite laminate