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.

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.

A hybrid method is presented for the analysis of layers, plates, and multi-layered systems consisting of isotropic and linear elastic materials. The problem is formulated for the general case of a multi-layered system using a total potential energy formulation and employing the layerwise laminate theory of Reddy. A one-dimensional finite element model is used for the analysis of the multi-layered system through its thickness, and integral Fourier transforms are used to obtain the exact solution for the in-plane problem. Explicit expressions are obtained for the fundamental solution of the typical infinite layer, which are applied in the two-dimensional boundary integral equation model to produce the integral representation of the solution. The boundary integral equation model is two-dimensional, displacement-based and assumes piecewise continuous distribution of the displacement components through the system's thickness. The developed model describes the three-dimensional displacement field, the stress field, the strains and the interlaminar stresses over the entire domain of the problem as continuous functions of the position. This detailed three-dimensional analysis is achieved by incorporating only contour integrals. The boundary integral equations are discretized using the boundary element method and a numerical model is developed for the single numerical layer (element). This model is extended to the case of a multilayered system by introducing appropriate continuity conditions at the interfaces between the layers (firmly bonded layers, or separation, slip and friction between the layers). Assembly of the element matrices yields the global system of equations, which can be solved via iterative techniques. In addition, numerical techniques are developed for the evaluation of the boundary and domain integrals involved in the construction of the element matrices. The singular boundary integrals are computed using a special coordinate transformation, along with a subdivision of the boundary element and a transformation of the Gauss points. The domain integrals (regular, singular or near-singular) are transformed to regular definite integrals along the boundary through a semi-analytical approach. The proposed method provides a simple, efficient, and versatile model for a three-dimensional analysis of thick plates or multilayered systems. It can also be used to study plates resting on elastic foundations or plates with internal supports. The proposed method can be applied in an obvious manner to anisotropic materials and vibration problems. / Ph. D.