Este trabalho apresenta a aplicação do Método de Otimização Topológica (MOT) considerando restrição de tensão mecânica em dois problemas de Engenharia: o projeto de estruturas mecânicas sujeitas a restrição de tensão e o projeto da distribuição de material em estruturas constituídas por Materiais com Gradação Funcional (MsGF). O MOT é um método numérico capaz de fornecer de forma automática o leiaute básico de uma estrutura mecânica para que esta atenda a um dado requisito de projeto, como o limite sobre a máxima tensão mecânica no componente. Os MsGF são materiais cujas propriedades variam gradualmente com a posição. Este gradiente de propriedades é obtido através da variação contínua da microestrutura formada por dois materiais diferentes. Neste trabalho o MOT foi implementado utilizando o modelo de material Solid Isotropic Microstructure with Penalization (SIMP) e o campo de densidades foi parametrizado utilizando a abordagem Aproximação Contínua da Distribuição de Material (ACDM). O modelo de material e utilizado em conjunto com um localizador de tensões, de modo a representar as tensões nas regiões com densidade intermediária. O projeto de estruturas tradicionais através do MOT possui dois problemas centrais aqui tratados: o fenômeno das topologias singulares, que consiste na incapacidade do algoritmo de otimização de retirar material de certas regiões da estrutura, onde a tensão mecânica supera o limite de tensão quando os valores da densidade tendem a zero, e o problema do grande número de restrições envolvidas, pois que a tensão mecânica é uma grandeza local e deve ser restrita em todos os pontos da estrutura. Para tratar o primeiro problema é utilizado o conceito de relaxação. Para o segundo são utilizadas duas abordagens: uma é a substituição das restrições locais por uma restrição global e a outra é a aplicação do Método do Lagrangeano Aumentado. Ambas foram implementadas e aplicadas para o projeto de estruturas planas e axissimétricas. No projeto da distribuição de material em estruturas constituídas por MsGF é utilizado um modelo de material baseado na interpolação dos limites de Hashin-Shtrikman. A partir deste modelo as tensões em cada fase são obtidas a partir das matrizes localizadoras de tensão. Para tratar o fenômeno das topologias singulares é proposto um índice estimativo de falha, baseado nas tensões de von Mises em cada fase da microestrutura, que evita tal problema. O grande número de restrições é tratado através da restrição global de tensão. Em ambos os problemas as formulações são apresentadas e sua eficiência é discutida através de exemplos numéricos. / This work presents the Topology Optimization Method (TOM) with stress constraint applied to two Engineering problems: the design of mechanical structures subjected to stress constraint and the design of material distribution in structures made of Functionally Graded Materials (FGMs). The TOM is a numerical method capable of synthesizing the basic layout of a mechanical structure accomplishing to a given design requirement, for example the maximum stress in the structure. The FGMs are materials with spatially varying properties, which are obtained through a continuum change of the microstructuremade of two different materials. In this work, the TOM was implemented with Solid Isotropic Microstructure with Penalization (SIMP) material model and the density field was parameterized with the Continuous Approximations of Material Distribution. To obtain the intermediate density stresses, the material model is applied together with a stress localization matrix. The design of mechanical structures through the TOM has two major problems: the singular topology phenomenon, which is characterized by the optimization algorithm impossibility of removing material from certain regions, where the stress overpasses the limiting stress when the density goes to zero, and the large number of constraints, once the stress is a local value that must be constrained everywhere in the structure. To deal with the first problem, it is applied the \"-realaxation concept, and for the second one two approaches are considered: one is to change the local stress constraint into a global stress constraint and the other is to apply the Augmented Lagrangian Method. Both approaches were implemented and applied to the design of plane and axisymmetric structures. In the design of material distribution in structures made of FGMs a material model based on Hashin-Shtrikman bounds is applied. From this model, stresses in each phase are obtained by the stress localization matrix. To deal with the singular topology phenomenon it is proposed a modified von Mises failure criteria index that avoids such problem. A global stress constraint is applied to deal with the large number of constraints. In both problems formulations are presented and their performance are discussed through numerical examples.

The objective of this study is to determine stress intensity factors (SIFs) for a crack in a radially graded FGM layer on a substrate. Functionally graded coating with an edge crack perpendicular to the interface and a homogeneous substrate are bonded together. In order to make the problem analytically tractable, geometry is modeled as an FGM strip attached to a homogeneous layer. Introducing the elastic foundation underneath the homogeneous layer, an FGM coating on a thin walled cylinder can be modeled. At first, governing equations are obtained from stress displacement and equilibrium equations. Then using an assumed form of solution in terms of Fourier Transforms for displacements and applying the boundary conditions, a singular integral equation is obtained for the mode-I problem. Solving this singular integral equation numerically, stress intensity factors are obtained as functions of crack length, strip thicknesses and inhomogeneity parameter.

Este trabalho apresenta a aplicação do Método de Otimização Topológica (MOT) considerando restrição de tensão mecânica em dois problemas de Engenharia: o projeto de estruturas mecânicas sujeitas a restrição de tensão e o projeto da distribuição de material em estruturas constituídas por Materiais com Gradação Funcional (MsGF). O MOT é um método numérico capaz de fornecer de forma automática o leiaute básico de uma estrutura mecânica para que esta atenda a um dado requisito de projeto, como o limite sobre a máxima tensão mecânica no componente. Os MsGF são materiais cujas propriedades variam gradualmente com a posição. Este gradiente de propriedades é obtido através da variação contínua da microestrutura formada por dois materiais diferentes. Neste trabalho o MOT foi implementado utilizando o modelo de material Solid Isotropic Microstructure with Penalization (SIMP) e o campo de densidades foi parametrizado utilizando a abordagem Aproximação Contínua da Distribuição de Material (ACDM). O modelo de material e utilizado em conjunto com um localizador de tensões, de modo a representar as tensões nas regiões com densidade intermediária. O projeto de estruturas tradicionais através do MOT possui dois problemas centrais aqui tratados: o fenômeno das topologias singulares, que consiste na incapacidade do algoritmo de otimização de retirar material de certas regiões da estrutura, onde a tensão mecânica supera o limite de tensão quando os valores da densidade tendem a zero, e o problema do grande número de restrições envolvidas, pois que a tensão mecânica é uma grandeza local e deve ser restrita em todos os pontos da estrutura. Para tratar o primeiro problema é utilizado o conceito de relaxação. Para o segundo são utilizadas duas abordagens: uma é a substituição das restrições locais por uma restrição global e a outra é a aplicação do Método do Lagrangeano Aumentado. Ambas foram implementadas e aplicadas para o projeto de estruturas planas e axissimétricas. No projeto da distribuição de material em estruturas constituídas por MsGF é utilizado um modelo de material baseado na interpolação dos limites de Hashin-Shtrikman. A partir deste modelo as tensões em cada fase são obtidas a partir das matrizes localizadoras de tensão. Para tratar o fenômeno das topologias singulares é proposto um índice estimativo de falha, baseado nas tensões de von Mises em cada fase da microestrutura, que evita tal problema. O grande número de restrições é tratado através da restrição global de tensão. Em ambos os problemas as formulações são apresentadas e sua eficiência é discutida através de exemplos numéricos. / This work presents the Topology Optimization Method (TOM) with stress constraint applied to two Engineering problems: the design of mechanical structures subjected to stress constraint and the design of material distribution in structures made of Functionally Graded Materials (FGMs). The TOM is a numerical method capable of synthesizing the basic layout of a mechanical structure accomplishing to a given design requirement, for example the maximum stress in the structure. The FGMs are materials with spatially varying properties, which are obtained through a continuum change of the microstructuremade of two different materials. In this work, the TOM was implemented with Solid Isotropic Microstructure with Penalization (SIMP) material model and the density field was parameterized with the Continuous Approximations of Material Distribution. To obtain the intermediate density stresses, the material model is applied together with a stress localization matrix. The design of mechanical structures through the TOM has two major problems: the singular topology phenomenon, which is characterized by the optimization algorithm impossibility of removing material from certain regions, where the stress overpasses the limiting stress when the density goes to zero, and the large number of constraints, once the stress is a local value that must be constrained everywhere in the structure. To deal with the first problem, it is applied the \"-realaxation concept, and for the second one two approaches are considered: one is to change the local stress constraint into a global stress constraint and the other is to apply the Augmented Lagrangian Method. Both approaches were implemented and applied to the design of plane and axisymmetric structures. In the design of material distribution in structures made of FGMs a material model based on Hashin-Shtrikman bounds is applied. From this model, stresses in each phase are obtained by the stress localization matrix. To deal with the singular topology phenomenon it is proposed a modified von Mises failure criteria index that avoids such problem. A global stress constraint is applied to deal with the large number of constraints. In both problems formulations are presented and their performance are discussed through numerical examples.

The Wavelet Finite Element Method (WFEM) involves combining the versatile wavelet analysis with the classical Finite Element Method (FEM) by utilizing the wavelet scaling functions as interpolating functions; providing an alternative to the conventional polynomial interpolation functions used in classical FEM. Wavelet analysis as a tool applied in WFEM has grown in popularity over the past decade and a half and the WFEM has demonstrated potential prowess to overcome some difficulties and limitations of FEM. This is particular for problems with regions of the solution domain where the gradient of the field variables are expected to vary fast or suddenly, leading to higher computational costs and/or inaccurate results. The properties of some of the various wavelet families such as compact support, multiresolution analysis (MRA), vanishing moments and the “two-scale” relations, make the use of wavelets in WFEM advantageous, particularly in the analysis of problems with strong nonlinearities, singularities and material property variations present. The wavelet based finite elements (WFEs) presented in this study, conceptually based on previous works, are constructed using the Daubechies and B-spline wavelet on the interval (BSWI) wavelet families. These two wavelet families possess the desired properties of multiresolution, compact support, the “two scale” relations and vanishing moments. The rod, beam and planar bar WFEs are used to study structural static and dynamic problems (moving load) via numerical examples. The dynamic analysis of functionally graded materials (FGMs) is further carried out through a new modified wavelet based finite element formulation using the Daubechies and BSWI wavelets, tailored for such classes of composite materials that have their properties varying spatially. Consequently, a modified algorithm of the multiscale Daubechies connection coefficients used in the formulation of the FGM elemental matrices and load vectors in wavelet space is presented and implemented in the formulation of the WFEs. The approach allows for the computation of the integral of the products of the Daubechies functions, and/or their derivatives, for different Daubechies function orders. The effects of varying the material distribution of a functionally graded (FG) beam on the natural frequency and dynamic response when subjected to a moving load for different velocity profiles are analysed. The dynamic responses of a FG beam resting on a viscoelastic foundation are also analysed for different material distributions, velocity and viscous damping profiles. The approximate solutions of the WFEM converge to the exact solution when the order and/or multiresolution scale of the WFE are increased. The results demonstrate that the Daubechies and B-spline based WFE solutions are highly accurate and require less number of elements than FEM due to the multiresolution property of WFEM. Furthermore, the applied moving load velocities and viscous damping influence the effects of varying the material distribution of FG beams on the dynamic response. Additional aspects of WFEM such as, the effect of altering the layout of the WFE and selection of the order of wavelet families to analyse static problems, are also presented in this study.

No description available.

Using experimental data obtained from standard fracture test configurations, theoretical and numerical tools are developed to mathematically describe non-self-similar progression of cracks without specifying an initial crack. A cohesive-decohesive zone model, similar to the cohesive zone model known in the fracture mechanics literature as the Dugdale-Barenblatt model, is adopted to represent the degradation of the material ahead of the crack tip. This model unifies strength-based crack initiation and fracture-mechanics-based crack progression. The cohesive-decohesive zone model is implemented with an interfacial surface material that consists of an upper and a lower surface that are connected by a continuous distribution of normal and tangential nonlinear elastic springs that act to resist either Mode I opening, Mode II sliding, Mode III sliding, or a mixed mode. The initiation of fracture is determined by the interfacial strength and the progression of the crack is determined by the critical energy release rate. The adhesive is idealized with an interfacial surface material to predict interfacial fracture. The interfacial surface material is positioned within the bulk material to predict discrete cohesive cracks. The interfacial surface material is implemented through an interface element, which is incorporated in ABAQUS using the user defined element (UEL) option. A procedure is established to formulate a rate dependent model based on experiments carried out on compact tension test specimens. The rate dependent model is incorporated into the interface element approach to capture the unstable crack growth observed in experiments under quasi-static loading conditions. The compact tension test gives the variation of the fracture toughness with the rate of loading, this information is processed and a relationship between the fracture toughness and the rate of the opening displacement is established. The cohesive-decohesive zone model is implemented through a material model to be used in an explicit code (LS-DYNA). Dynamic simulations of the standard test configurations for Mode I (Double Cantilever Beam) and Mode II (End Load Split) are carried out using the explicit code. Verification of these coupon tests leads to the crash analysis of realistic structures like the square composite tube. Analyses of bonded and unbonded square tubes are presented. These tubes shows a very uncharacteristic failure mode: the composite material disintegrates on impact, and this has been captured in the analysis. Disadvantages of the interface element approach are well documented in the literature. An alternative method, known as the Extended Finite Element Method (XFEM), is implemented here through an eight-noded quadrilateral plane strain element. The method, based on the partition-of-unity, is used to study simple test configuration like the three-point bend problem and a double cantilever beam. Functionally graded materials are also simulated and the results are compared to the experimental results available in the literature. / Ph. D.

The goal of this research was to create a microstructure graded glass-ceramic and investigate the resulting properties as a function of crystallization processing. The desired glass-ceramic was a lithium disilicate material that has a crystallization gradient across the sample, leading to functionally graded properties as a result of the microstructure gradient. Samples were prepared by melting and pouring glass at 1400°C, annealing at 400°C for 48 hours, and nucleating at 480°C for 2 hours. To ensure that crystallization would not occur homogeneously throughout the sample, a temperature gradient was imposed during crystallization. Samples were crystallized on a self-constructed resistance wire furnace that was open to air. Several crystallization processing parameters were tested, including high temperature for a short time and low temperature for a long time. Samples were ground and polished to 0.25 microns before characterization methods were performed. Scanning electron microscopy (SEM) showed the microstructure transition across the sample cross section, with crystals present on the crystalline side and only nuclei present on the glassy side. Raman spectroscopy showed a transformation of the characteristic spectra across the sample cross section, with defined, high-intensity peaks on the crystalline side and broad, low-intensity peaks on the glassy side. Microhardness showed a slight transition in hardness values across the sample cross section, however the variability was too great to draw any conclusions. The characterization methods showed that the desired material was created and the resulting properties were a function of the crystallization processing parameters. / Master of Science

In this study the plane linear elastic problem of a functionally graded layer which contains periodic cracks is considered. The main objective of this study is to determine the thermal stress intensity factors for edge cracks. In order to find an analytic solution, Young&rsquo / s modulus and thermal conductivity are assumed to be varying exponentially across the thickness, whereas Poisson ratio and thermal diffusivity are taken as constant. First, one dimensional transient and steady state conduction problems are solved (heat flux being across the thickness) to determine the temperature distribution and the thermal stresses in a crack free layer. Then, the thermal stress distributions at the locations of the cracks are applied as crack surface tractions in the elasticity problem to find the stress intensity factors. By defining an appropriate auxiliary variable, elasticity problem is reduced to a singular integral equation, which is solved numerically. The influence of such parameters as the grading, crack length and crack period on the stress intensity factors is investigated.

No description available.

The main objective of this study is to evaluate mixed mode stress intensity factors for inclined embedded cracks in functionally graded materials. Fracture analysis of inclined cracks requires the calculation of both Mode I and Mode II stress intensity factors ( I K , II K ). In this study, k J -integral is used to calculate I K and II K . Equivalent domain integral approach is utilized to evaluate the k J - integral around the crack tip. The present study aims at developing a finite element model to study inclined crack problems in graded media under thermomechanical loading. A two dimensional finite element model is developed for inclined cracks located in a functionally graded medium. Structural and thermal problems are solved using two dimensional finite elements namely 8- noded triangles. Material properties are sampled directly at the integration points of the elements, as required by the numerical integral evaluation. The main results of the study are the stress intensity factors at the crack tip for functionally graded materials subjected to thermomechanical loading.