O propósito principal do presente trabalho é estabelecer diretrizes para o uso de elementos finitos sólidos na modelagem de estruturas delgadas submetidas à flexão. Por meio de um conjunto de problemas estudou-se o travamento de elementos finitos da elasticidade bi e tridimensionais, quando usados para modelar estruturas delgadas submetidas à flexão. Foram usados elementos planos de 4, 8 e 9 nós e elementos tridimensionais de 4, 10, 11, 8, 20 e 27 nós, considerando elementos não distorcidos e distorcidos. A análise dos resultados nesses modelos permitiram entender melhor o comportamento e o uso de elementos sólidos em estruturas delgadas. / The main objective of the present work is to establish guidelines for the use of solid finite elements in the modeling of thin structures submitted to bending. By means of a set of problems the locking of elements three-dimensional and two-dimensional elasticity was studied, when used to model problems thin structures subjected to bending. We use 2-D displacements-based elements of 4, 8 and 9 nodes and 3-D displacement based of 4, 10, 11, 8, 20 and 27 nodes, considering undistorted and distorted elements. The analysis of the results of these models lead to the understanding of the behavior and the use of solid elements in thin structures.

O propósito principal do presente trabalho é estabelecer diretrizes para o uso de elementos finitos sólidos na modelagem de estruturas delgadas submetidas à flexão. Por meio de um conjunto de problemas estudou-se o travamento de elementos finitos da elasticidade bi e tridimensionais, quando usados para modelar estruturas delgadas submetidas à flexão. Foram usados elementos planos de 4, 8 e 9 nós e elementos tridimensionais de 4, 10, 11, 8, 20 e 27 nós, considerando elementos não distorcidos e distorcidos. A análise dos resultados nesses modelos permitiram entender melhor o comportamento e o uso de elementos sólidos em estruturas delgadas. / The main objective of the present work is to establish guidelines for the use of solid finite elements in the modeling of thin structures submitted to bending. By means of a set of problems the locking of elements three-dimensional and two-dimensional elasticity was studied, when used to model problems thin structures subjected to bending. We use 2-D displacements-based elements of 4, 8 and 9 nodes and 3-D displacement based of 4, 10, 11, 8, 20 and 27 nodes, considering undistorted and distorted elements. The analysis of the results of these models lead to the understanding of the behavior and the use of solid elements in thin structures.

Ce travail consiste à développer un logiciel métier basé sur la modélisation thermomécanique linéaire de l'échangeur tout en intégrant les techniques d'homogénéisation. Une méthodologie de modélisation de l'échangeur constitué par un empilement de différentes ondes et tôles brasées, est donc adoptée. Cette méthodologie suppose que le comportement global de chaque couche d'ondes et tôles, est encadré par deux comportements limites déterminés par des approches mécanique périodique (HMP) et cinématique périodique (HCP). Ces deux approches sont mises en œuvre pour l’application numérique tout en étudiant le chargement interne global dû à la température et à la pression. Un outil d'homogénéisation (HomPass) est ensuite développé afin de déterminer automatiquement les comportements équivalents à chaque onde et tôle brasées. Cela contribue au développement de l'outil métier final (SiTEME) dédié à l'étude thermomécanique globale de l'échangeur. / This work consists of developing a software tool based on the linear thermomechanical modeling of the heat exchanger using homogenization techniques. A methodology for modeling the heat exchanger constituted by stacking of different brazed fins and sheets is adopted. This methodology assumes that the global behavior of every layer of fins and sheets is found between two behavior limits determined by periodic mechanical approach (HMP) and periodic kinematical approach (HCP). These techniques are implemented for numerical application while studying global loading due to the temperature and pressure internal loads. A homogenization tool (HomPass) is then developed in order to calculate automatically equivalent behaviors to each layer of brazed fins and sheets. That contributes to the development of the final software tool (SiTEME) dedicated to the global thermomechanical study of the heat exchanger.

O objetivo precípuo deste estudo é a implementação computacional de modelos constitutivos elásticos e elastoplásticos para materiais com gradação funcional em regime de grandes deslocamentos e elevadas deformações. Para simular numericamente um problema estrutural, são empregados aqui elementos finitos sólidos (tetraédrico e hexaédrico) com ordem de aproximação polinomial qualquer. Grandezas da Mecânica Não Linear do Contínuo, como deformação e tensão, são utilizadas na formulação deste estudo. Para reproduzir os grandes deslocamentos, é empregada a análise não linear geométrica. A descrição adotada aqui é a Lagrangiana total, e o equilíbrio da estrutura é expresso pelo Princípio da Mínima Energia Potencial Total. Com relação à resposta elástica do material, são usadas leis constitutivas hiperelásticas, nas quais a relação tensão-deformação é obtida a partir de um potencial escalar. O comportamento elastoplástico do material é definido pela decomposição da deformação nas parcelas elástica e plástica, pelo critério de plastificação de von-Mises, pela lei de fluxo associativa, pelas condições de consistência e de complementaridade, pelo parâmetro de encruamento isotrópico e pelo tensor das tensões inversas, relacionado ao encruamento cinemático. Duas formulações elastoplásticas são utilizadas aqui: a de Green-Naghdi, na qual a deformação é decomposta de forma aditiva; e a hiperelastoplástica, em que o gradiente é decomposto de forma multiplicativa. É empregado também o conceito de material com gradação funcional (GF), a qual é definida como a variação gradual (contínua e suave) das propriedades constitutivas do material. A solução numérica do equilíbrio de forças é feita via método iterativo de Newton-Raphson. Para satisfazer o critério de plastificação, são utilizadas as estratégias de previsão elástica, e de correção plástica via algoritmos de retorno. Basicamente foram desenvolvidos cinco programas computacionais: o gerador automático das funções de forma; o gerador de malhas de elementos finitos sólidos; o código para análise de materiais em regime elástico; o código para análise de materiais em regime elastoplástico; e o programa de pós-processamento. Além desses, o aluno teve contato com os programas EPIM3D e DD3IMP ao longo de seu estágio de doutorado na Universidade de Coimbra (Portugal). Os programas EPIM3D e DD3IMP são empregados para analisar, respectivamente, materiais em regime elastoplástico, e processos de conformação de metais. Para o problema da barra sob tração uniaxial uniforme, são descritas equações e soluções analíticas para materiais homogêneos e com GF em regime elastoplástico. Para reduzir o tempo de simulação, foi empregada a programação em paralelo. De acordo com os resultados das simulações numéricas, as principais conclusões são: o refinamento da malha de elementos finitos melhora a precisão dos resultados para materiais em regimes elástico e elastoplástico; as formulações elastoplásticas de Green-Naghdi e hiperelastoplástica parecem ser equivalentes para pequenas deformações; a formulação hiperelastoplástica é equivalente ao modelo mecânico dos programas EPIM3D e DD3IMP para materiais em regime de pequenas deformações elásticas; foram constatados ganhos significativos, em termos de tempo de simulação, com a paralelização dos códigos computacionais de análise estrutural; e os programas desenvolvidos são capazes de simular - com precisão - problemas complexos, como a membrana de Cook e o cilindro fino transversalmente tracionado. / The main objective of this study is the computational implementation of elastic and elastoplastic constitutive models for functionally graded materials in large deformation regime. In order to numerically simulate a structural problem, the finite elements used are solids (tetrahedric and hexahedric) of any order of approximation. Entities from Nonlinear Continnum Mechanics, as strain and stress, are used in the present formulation. To reproduce the finite displacements, the geometrically nonlinear analysis is employed. The description adopted here is the total Lagrangian, and the structural equilibrium is expressed by means of the Principal of Minimum Total Potential Energy. Regarding the elastic material response, hyperelastic constitutive laws are used, in which the stress-strain relation is obtained from a scalar potential. The elastoplastic material behavior is defined by the strain decomposition in the elastic and plastic parts, by the von-Mises yield criterion, by the associative flow law, by the consistency and complementarity conditions, by the isotropic hardening parameter, and by the backstress tensor, related to the kinematic hardening. Two elastoplastic formulations are used here: the Green-Naghdi one, in which the strain is additively decomposed; and the hyperelastoplasticiy, in which the gradient is multiplicatively decomposed. The concept of functionally graded (FG) material, in which the constitutive properties vary gradually (continuous and smoothly), is also used. The numerical solution of the forces equilibrium is obtained via Newton-Raphson iterative procedure. In order to satisfy the yield criterion, the strategies of elastic prediction and plastic correction (via return algorithms) are used. Basically, five computer codes have been developed: the automatic shape functions generator; the solid mesh generator; the code for analysis of materials in the elastic regime; the code for analysis of materials in the elastoplastic regime; and the post-processor. Besides these, the student had contact with the programs EPIM3D and DD3IMP during his doctoral stage in the University of Coimbra (Portugal). The programs EPIM3D and DD3IMP are employed to analyze, respectively, materials in the elastoplastic regime, and sheet-metal forming processes. For the problem of the bar under uniform uniaxial tension, equations and analytical solutions are described for homogeneous and FG materials. To reduce the simulation time, the parallel programming has been employed. According to the numerical simulation results, the main conclusions are: the results accuracy is improved with mesh refinement for materials in the elastic and elastoplastic regimes; the Green-Naghdi elastoplastic formulation and the hyperelastoplasticity appear to be equivalent for small strains; the hyperelastoplastic formulation is equivalent to the mechanical model of the programs EPIM3D and DD3IMP for materials the small elastic strains regime; simulation time reduction has been obtained with the parallelization of the computer codes for structural analysis; the developed programs are capable of simulating, precisely, complex problems, such as the Cook\'s membrane and the pulled thin cylinder.

O objetivo precípuo deste estudo é a implementação computacional de modelos constitutivos elásticos e elastoplásticos para materiais com gradação funcional em regime de grandes deslocamentos e elevadas deformações. Para simular numericamente um problema estrutural, são empregados aqui elementos finitos sólidos (tetraédrico e hexaédrico) com ordem de aproximação polinomial qualquer. Grandezas da Mecânica Não Linear do Contínuo, como deformação e tensão, são utilizadas na formulação deste estudo. Para reproduzir os grandes deslocamentos, é empregada a análise não linear geométrica. A descrição adotada aqui é a Lagrangiana total, e o equilíbrio da estrutura é expresso pelo Princípio da Mínima Energia Potencial Total. Com relação à resposta elástica do material, são usadas leis constitutivas hiperelásticas, nas quais a relação tensão-deformação é obtida a partir de um potencial escalar. O comportamento elastoplástico do material é definido pela decomposição da deformação nas parcelas elástica e plástica, pelo critério de plastificação de von-Mises, pela lei de fluxo associativa, pelas condições de consistência e de complementaridade, pelo parâmetro de encruamento isotrópico e pelo tensor das tensões inversas, relacionado ao encruamento cinemático. Duas formulações elastoplásticas são utilizadas aqui: a de Green-Naghdi, na qual a deformação é decomposta de forma aditiva; e a hiperelastoplástica, em que o gradiente é decomposto de forma multiplicativa. É empregado também o conceito de material com gradação funcional (GF), a qual é definida como a variação gradual (contínua e suave) das propriedades constitutivas do material. A solução numérica do equilíbrio de forças é feita via método iterativo de Newton-Raphson. Para satisfazer o critério de plastificação, são utilizadas as estratégias de previsão elástica, e de correção plástica via algoritmos de retorno. Basicamente foram desenvolvidos cinco programas computacionais: o gerador automático das funções de forma; o gerador de malhas de elementos finitos sólidos; o código para análise de materiais em regime elástico; o código para análise de materiais em regime elastoplástico; e o programa de pós-processamento. Além desses, o aluno teve contato com os programas EPIM3D e DD3IMP ao longo de seu estágio de doutorado na Universidade de Coimbra (Portugal). Os programas EPIM3D e DD3IMP são empregados para analisar, respectivamente, materiais em regime elastoplástico, e processos de conformação de metais. Para o problema da barra sob tração uniaxial uniforme, são descritas equações e soluções analíticas para materiais homogêneos e com GF em regime elastoplástico. Para reduzir o tempo de simulação, foi empregada a programação em paralelo. De acordo com os resultados das simulações numéricas, as principais conclusões são: o refinamento da malha de elementos finitos melhora a precisão dos resultados para materiais em regimes elástico e elastoplástico; as formulações elastoplásticas de Green-Naghdi e hiperelastoplástica parecem ser equivalentes para pequenas deformações; a formulação hiperelastoplástica é equivalente ao modelo mecânico dos programas EPIM3D e DD3IMP para materiais em regime de pequenas deformações elásticas; foram constatados ganhos significativos, em termos de tempo de simulação, com a paralelização dos códigos computacionais de análise estrutural; e os programas desenvolvidos são capazes de simular - com precisão - problemas complexos, como a membrana de Cook e o cilindro fino transversalmente tracionado. / The main objective of this study is the computational implementation of elastic and elastoplastic constitutive models for functionally graded materials in large deformation regime. In order to numerically simulate a structural problem, the finite elements used are solids (tetrahedric and hexahedric) of any order of approximation. Entities from Nonlinear Continnum Mechanics, as strain and stress, are used in the present formulation. To reproduce the finite displacements, the geometrically nonlinear analysis is employed. The description adopted here is the total Lagrangian, and the structural equilibrium is expressed by means of the Principal of Minimum Total Potential Energy. Regarding the elastic material response, hyperelastic constitutive laws are used, in which the stress-strain relation is obtained from a scalar potential. The elastoplastic material behavior is defined by the strain decomposition in the elastic and plastic parts, by the von-Mises yield criterion, by the associative flow law, by the consistency and complementarity conditions, by the isotropic hardening parameter, and by the backstress tensor, related to the kinematic hardening. Two elastoplastic formulations are used here: the Green-Naghdi one, in which the strain is additively decomposed; and the hyperelastoplasticiy, in which the gradient is multiplicatively decomposed. The concept of functionally graded (FG) material, in which the constitutive properties vary gradually (continuous and smoothly), is also used. The numerical solution of the forces equilibrium is obtained via Newton-Raphson iterative procedure. In order to satisfy the yield criterion, the strategies of elastic prediction and plastic correction (via return algorithms) are used. Basically, five computer codes have been developed: the automatic shape functions generator; the solid mesh generator; the code for analysis of materials in the elastic regime; the code for analysis of materials in the elastoplastic regime; and the post-processor. Besides these, the student had contact with the programs EPIM3D and DD3IMP during his doctoral stage in the University of Coimbra (Portugal). The programs EPIM3D and DD3IMP are employed to analyze, respectively, materials in the elastoplastic regime, and sheet-metal forming processes. For the problem of the bar under uniform uniaxial tension, equations and analytical solutions are described for homogeneous and FG materials. To reduce the simulation time, the parallel programming has been employed. According to the numerical simulation results, the main conclusions are: the results accuracy is improved with mesh refinement for materials in the elastic and elastoplastic regimes; the Green-Naghdi elastoplastic formulation and the hyperelastoplasticity appear to be equivalent for small strains; the hyperelastoplastic formulation is equivalent to the mechanical model of the programs EPIM3D and DD3IMP for materials the small elastic strains regime; simulation time reduction has been obtained with the parallelization of the computer codes for structural analysis; the developed programs are capable of simulating, precisely, complex problems, such as the Cook\'s membrane and the pulled thin cylinder.