Matching Items (10)
Filtering by
- All Subjects: Finite element method
- Creators: Neithalath, Narayanan
- Creators: Peralta, Pedro
Description
This study focuses on implementing probabilistic nature of material properties (Kevlar® 49) to the existing deterministic finite element analysis (FEA) of fabric based engine containment system through Monte Carlo simulations (MCS) and implementation of probabilistic analysis in engineering designs through Reliability Based Design Optimization (RBDO). First, the emphasis is on experimental data analysis focusing on probabilistic distribution models which characterize the randomness associated with the experimental data. The material properties of Kevlar® 49 are modeled using experimental data analysis and implemented along with an existing spiral modeling scheme (SMS) and user defined constitutive model (UMAT) for fabric based engine containment simulations in LS-DYNA. MCS of the model are performed to observe the failure pattern and exit velocities of the models. Then the solutions are compared with NASA experimental tests and deterministic results. MCS with probabilistic material data give a good prospective on results rather than a single deterministic simulation results. The next part of research is to implement the probabilistic material properties in engineering designs. The main aim of structural design is to obtain optimal solutions. In any case, in a deterministic optimization problem even though the structures are cost effective, it becomes highly unreliable if the uncertainty that may be associated with the system (material properties, loading etc.) is not represented or considered in the solution process. Reliable and optimal solution can be obtained by performing reliability optimization along with the deterministic optimization, which is RBDO. In RBDO problem formulation, in addition to structural performance constraints, reliability constraints are also considered. This part of research starts with introduction to reliability analysis such as first order reliability analysis, second order reliability analysis followed by simulation technique that are performed to obtain probability of failure and reliability of structures. Next, decoupled RBDO procedure is proposed with a new reliability analysis formulation with sensitivity analysis, which is performed to remove the highly reliable constraints in the RBDO, thereby reducing the computational time and function evaluations. Followed by implementation of the reliability analysis concepts and RBDO in finite element 2D truss problems and a planar beam problem are presented and discussed.
ContributorsDeivanayagam, Arumugam (Author) / Rajan, Subramaniam D. (Thesis advisor) / Mobasher, Barzin (Committee member) / Neithalath, Narayanan (Committee member) / Arizona State University (Publisher)
Created2012
Description
Shock loading is a complex phenomenon that can lead to failure mechanisms such as strain localization, void nucleation and growth, and eventually spall fracture. Studying incipient stages of spall damage is of paramount importance to accurately determine initiation sites in the material microstructure where damage will nucleate and grow and to formulate continuum models that account for the variability of the damage process due to microstructural heterogeneity. The length scale of damage with respect to that of the surrounding microstructure has proven to be a key aspect in determining sites of failure initiation. Correlations have been found between the damage sites and the surrounding microstructure to determine the preferred sites of spall damage, since it tends to localize at and around the regions of intrinsic defects such as grain boundaries and triple points. However, considerable amount of work still has to be done in this regard to determine the physics driving the damage at these intrinsic weak sites in the microstructure. The main focus of this research work is to understand the physical mechanisms behind the damage localization at these preferred sites. A crystal plasticity constitutive model is implemented with different damage criteria to study the effects of stress concentration and strain localization at the grain boundaries. A cohesive zone modeling technique is used to include the intrinsic strength of the grain boundaries in the simulations. The constitutive model is verified using single elements tests, calibrated using single crystal impact experiments and validated using bicrystal and multicrystal impact experiments. The results indicate that strain localization is the predominant driving force for damage initiation and evolution. The microstructural effects on theses damage sites are studied to attribute the extent of damage to microstructural features such as grain orientation, misorientation, Taylor factor and the grain boundary planes. The finite element simulations show good correlation with the experimental results and can be used as the preliminary step in developing accurate probabilistic models for damage nucleation.
ContributorsKrishnan, Kapil (Author) / Peralta, Pedro (Thesis advisor) / Mignolet, Marc (Committee member) / Sieradzki, Karl (Committee member) / Jiang, Hanqing (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2013
Description
Dwindling energy resources and associated environmental costs have resulted in a serious need to design and construct energy efficient buildings. One of the strategies to develop energy efficient structural materials is through the incorporation of phase change materials (PCM) in the host matrix. This research work presents details of a finite element-based framework that is used to study the thermal performance of structural precast concrete wall elements with and without a layer of phase change material. The simulation platform developed can be implemented for a wide variety of input parameters. In this study, two different locations in the continental United States, representing different ambient temperature conditions (corresponding to hot, cold and typical days of the year) are studied. Two different types of concrete - normal weight and lightweight, different PCM types, gypsum wallboard's with varying PCM percentages and different PCM layer thicknesses are also considered with an aim of understanding the energy flow across the wall member. Effect of changing PCM location and prolonged thermal loading are also studied. The temperature of the inside face of the wall and energy flow through the inside face of the wall, which determines the indoor HVAC energy consumption are used as the defining parameters. An ad-hoc optimization scheme is also implemented where the PCM thickness is fixed but its location and properties are varied. Numerical results show that energy savings are possible with small changes in baseline values, facilitating appropriate material design for desired characteristics.
ContributorsHembade, Lavannya Babanrao (Author) / Neithalath, Narayanan (Thesis advisor) / Rajan, Subramaniam D. (Thesis advisor) / Mobasher, Barzin (Committee member) / Arizona State University (Publisher)
Created2012
Description
Tall buildings are spreading across the globe at an ever-increasing rate (www.ctbuh.org). The global number of buildings 200m or more in height has risen from 286 to 602 in the last decade alone. The increasing complexity of building architecture poses unique challenges in the structural design of modern tall buildings. Hence, innovative structural systems need to be evaluated to create an economical design that satisfies multiple design criteria. Design using traditional trial-and-error approach can be extremely time-consuming and the resultant design uneconomical. Thus, there is a need for an efficient numerical optimization tool that can explore and generate several design alternatives in the preliminary design phase which can lead to a more desirable final design. In this study, we present the details of a tool that can be very useful in preliminary design optimization - finite element modeling, design optimization, translating design code requirements into components of the FE and design optimization models, and pre-and post-processing to verify the veracity of the model. Emphasis is placed on development and deployment of various FE models (static, modal and dynamic analyses; linear, beam and plate/shell finite elements), design optimization problem formulation (sizing, shape, topology and material selection optimization) and numerical optimization tools (gradient-based and evolutionary optimization methods) [Rajan, 2001]. The design optimization results of full scale three dimensional buildings subject to multiple design criteria including stress, serviceability and dynamic response are discussed.
ContributorsSirigiri, Mamatha (Author) / Rajan, Subramaniam D. (Thesis advisor) / Neithalath, Narayanan (Committee member) / Mobasher, Barzin (Committee member) / Arizona State University (Publisher)
Created2014
Description
Laminated composite materials are used in aerospace, civil and mechanical structural systems due to their superior material properties compared to the constituent materials as well as in comparison to traditional materials such as metals. Laminate structures are composed of multiple orthotropic material layers bonded together to form a single performing part. As such, the layup design of the material largely influences the structural performance. Optimization techniques such as the Genetic Algorithm (GA), Differential Evolution (DE), the Method of Feasible Directions (MFD), and others can be used to determine the optimal laminate composite material layup. In this thesis, sizing, shape and topology design optimization of laminated composites is carried out. Sizing optimization, such as the layer thickness, topology optimization, such as the layer orientation and material and the number of layers present, and shape optimization of the overall composite part contribute to the design optimization process of laminates. An optimization host program written in C++ has been developed to implement the optimization methodology of both population based and numerical gradient based methods. The performance of the composite structural system is evaluated through explicit finite element analysis of shell elements carried out using LS-DYNA. Results from numerical examples demonstrate that optimization design processes can significantly improve composite part performance through implementation of optimum material layup and part shape.
ContributorsMika, Krista (Author) / Rajan, Subramaniam D. (Thesis advisor) / Neithalath, Narayanan (Committee member) / Mobasher, Barzin (Committee member) / Arizona State University (Publisher)
Created2014
Description
Uranium Dioxide (UO2) is a significant nuclear fission fuel, which is widely used
in nuclear reactors. Understanding the influence of microstructure on thermo-mechanical behavior of UO2 is extremely important to predict its performance. In particular, evaluating mechanical properties, such as elasticity, plasticity and creep at sub-grain length scales is key to developing this understanding as well as building multi-scale models of fuel behavior with predicting capabilities. In this work, modeling techniques were developed to study effects of microstructure on Young’s modulus, which was selected as a key representative property that affects overall mechanical behavior, using experimental data obtained from micro-cantilever bending testing as benchmarks. Beam theory was firstly introduced to calculate Young's modulus of UO2 from the experimental data and then three-dimensional finite element models of the micro-cantilever beams were constructed to simulate bending tests in UO2 at room temperature. The influence of the pore distribution was studied to explain the discrepancy between predicted values and experimental results. Results indicate that results of tests are significantly affected by porosity given that both pore size and spacing in the samples are of the order of the micro-beam dimensions. Microstructure reconstruction was conducted with images collected from three-dimensional serial sectioning using focused ion beam (FIB) and electron backscattering diffraction (EBSD) and pore clusters were placed at different locations along the length of the beam. Results indicate that the presence of pore clusters close to the substrate, i.e., the clamp of the micro-cantilever beam, has the strongest effect on load-deflection behavior, leading to a reduction of stiffness that is the largest for any location of the pore cluster. Furthermore, it was also found from both numerical and i
analytical models that pore clusters located towards the middle of the span and close to the end of the beam only have a very small effect on the load-deflection behavior, and it is concluded that better estimates of Young's modulus can be obtained from micro- cantilever experiments by using microstructurally explicit models that account for porosity in about one half of the beam length close to the clamp. This, in turn, provides an avenue to simplify micro-scale experiments and their analysis.
in nuclear reactors. Understanding the influence of microstructure on thermo-mechanical behavior of UO2 is extremely important to predict its performance. In particular, evaluating mechanical properties, such as elasticity, plasticity and creep at sub-grain length scales is key to developing this understanding as well as building multi-scale models of fuel behavior with predicting capabilities. In this work, modeling techniques were developed to study effects of microstructure on Young’s modulus, which was selected as a key representative property that affects overall mechanical behavior, using experimental data obtained from micro-cantilever bending testing as benchmarks. Beam theory was firstly introduced to calculate Young's modulus of UO2 from the experimental data and then three-dimensional finite element models of the micro-cantilever beams were constructed to simulate bending tests in UO2 at room temperature. The influence of the pore distribution was studied to explain the discrepancy between predicted values and experimental results. Results indicate that results of tests are significantly affected by porosity given that both pore size and spacing in the samples are of the order of the micro-beam dimensions. Microstructure reconstruction was conducted with images collected from three-dimensional serial sectioning using focused ion beam (FIB) and electron backscattering diffraction (EBSD) and pore clusters were placed at different locations along the length of the beam. Results indicate that the presence of pore clusters close to the substrate, i.e., the clamp of the micro-cantilever beam, has the strongest effect on load-deflection behavior, leading to a reduction of stiffness that is the largest for any location of the pore cluster. Furthermore, it was also found from both numerical and i
analytical models that pore clusters located towards the middle of the span and close to the end of the beam only have a very small effect on the load-deflection behavior, and it is concluded that better estimates of Young's modulus can be obtained from micro- cantilever experiments by using microstructurally explicit models that account for porosity in about one half of the beam length close to the clamp. This, in turn, provides an avenue to simplify micro-scale experiments and their analysis.
ContributorsGong, Bowen (Author) / Peralta, Pedro (Thesis advisor) / Rajagopalan, Jagannathan (Committee member) / Solanki, Kiran (Committee member) / Arizona State University (Publisher)
Created2015
Description
This study employs a finite element method based modeling of cementitious composite microstructure to study the effect of presence of inclusions on the stress distribution and the constitutive response of the composite. A randomized periodic microstructure combined with periodic boundary conditions forms the base of the finite element models. Inclusion properties of quartz and light weight aggregates of size 600μm obtained from literature were made use of to study the effect of their material (including inclusion stiffness, stiffness of interfacial transition zone and matrix stiffening) and geometric properties (volume fraction of inclusion, particle size distribution of inclusion and thickness of the interfacial transition zone) on the composite. Traction-separation relationship was used to incorporate the effect of debonding at the interface of the matrix and the inclusion to study the effect on stress distribution in the microstructure. The stress distributions observed upon conducting a finite element analysis are caused due to the stiffness mismatch in both the quartz and the light weight aggregates as expected. The constitutive response of the composite microstructure is found to be in good conformance with semi-analytical models as well as experimental values. The effect of debonding throws up certain important observations on the stress distributions in the microstructure based on the stress concentrations and relaxations caused by the stiffness of the individual components of the microstructure. The study presented discusses the different micromechanical models employed, their applicability and suitability to correctly predict the composite constitutive response.
ContributorsMaroli, Amit (Author) / Neithalath, Narayanan (Thesis advisor) / Rajan, Subramanium (Committee member) / Mobasher, Barzin (Committee member) / Arizona State University (Publisher)
Created2016
Description
Composite materials are now beginning to provide uses hitherto reserved for metals in structural systems such as airframes and engine containment systems, wraps for repair and rehabilitation, and ballistic/blast mitigation systems. These structural systems are often subjected to impact loads and there is a pressing need for accurate prediction of deformation, damage and failure. There are numerous material models that have been developed to analyze the dynamic impact response of polymer matrix composites. However, there are key features that are missing in those models that prevent them from providing accurate predictive capabilities. In this dissertation, a general purpose orthotropic elasto-plastic computational constitutive material model has been developed to predict the response of composites subjected to high velocity impacts. The constitutive model is divided into three components – deformation model, damage model and failure model, with failure to be added at a later date. The deformation model generalizes the Tsai-Wu failure criteria and extends it using a strain-hardening-based orthotropic yield function with a non-associative flow rule. A strain equivalent formulation is utilized in the damage model that permits plastic and damage calculations to be uncoupled and capture the nonlinear unloading and local softening of the stress-strain response. A diagonal damage tensor is defined to account for the directionally dependent variation of damage. However, in composites it has been found that loading in one direction can lead to damage in multiple coordinate directions. To account for this phenomena, the terms in the damage matrix are semi-coupled such that the damage in a particular coordinate direction is a function of the stresses and plastic strains in all of the coordinate directions. The overall framework is driven by experimental tabulated temperature and rate-dependent stress-strain data as well as data that characterizes the damage matrix and failure. The developed theory has been implemented in a commercial explicit finite element analysis code, LS-DYNA®, as MAT213. Several verification and validation tests using a commonly available carbon-fiber composite, Toyobo’s T800/F3900, have been carried and the results show that the theory and implementation are efficient, robust and accurate.
ContributorsHoffarth, Canio (Author) / Rajan, Subramaniam D. (Thesis advisor) / Goldberg, Robert (Committee member) / Neithalath, Narayanan (Committee member) / Mobasher, Barzin (Committee member) / Liu, Yongming (Committee member) / Arizona State University (Publisher)
Created2016
Description
The focus of this investigation includes three aspects. First, the development of nonlinear reduced order modeling techniques for the prediction of the response of complex structures exhibiting "large" deformations, i.e. a geometrically nonlinear behavior, and modeled within a commercial finite element code. The present investigation builds on a general methodology, successfully validated in recent years on simpler panel structures, by developing a novel identification strategy of the reduced order model parameters, that enables the consideration of the large number of modes needed for complex structures, and by extending an automatic strategy for the selection of the basis functions used to represent accurately the displacement field. These novel developments are successfully validated on the nonlinear static and dynamic responses of a 9-bay panel structure modeled within Nastran. In addition, a multi-scale approach based on Component Mode Synthesis methods is explored. Second, an assessment of the predictive capabilities of nonlinear reduced order models for the prediction of the large displacement and stress fields of panels that have a geometric discontinuity; a flat panel with a notch was used for this assessment. It is demonstrated that the reduced order models of both virgin and notched panels provide a close match of the displacement field obtained from full finite element analyses of the notched panel for moderately large static and dynamic responses. In regards to stresses, it is found that the notched panel reduced order model leads to a close prediction of the stress distribution obtained on the notched panel as computed by the finite element model. Two enrichment techniques, based on superposition of the notch effects on the virgin panel stress field, are proposed to permit a close prediction of the stress distribution of the notched panel from the reduced order model of the virgin one. A very good prediction of the full finite element results is achieved with both enrichments for static and dynamic responses. Finally, computational challenges associated with the solution of the reduced order model equations are discussed. Two alternatives to reduce the computational time for the solution of these problems are explored.
ContributorsPerez, Ricardo Angel (Author) / Mignolet, Marc (Thesis advisor) / Oswald, Jay (Committee member) / Spottswood, Stephen (Committee member) / Peralta, Pedro (Committee member) / Jiang, Hanqing (Committee member) / Arizona State University (Publisher)
Created2012
Finite element analysis of silicon thin films on soft substrates as anodes for lithium ion batteries
Description
The wide-scale use of green technologies such as electric vehicles has been slowed due to insufficient means of storing enough portable energy. Therefore it is critical that efficient storage mediums be developed in order to transform abundant renewable energy into an on-demand source of power. Lithium (Li) ion batteries are seeing a stream of improvements as they are introduced into many consumer electronics, electric vehicles and aircraft, and medical devices. Li-ion batteries are well suited for portable applications because of their high energy-to-weight ratios, high energy densities, and reasonable life cycles. Current research into Li-ion batteries is focused on enhancing its energy density, and by changing the electrode materials, greater energy capacities can be realized. Silicon (Si) is a very attractive option because it has the highest known theoretical charge capacity. Current Si anodes, however, suffer from early capacity fading caused by pulverization from the stresses induced by large volumetric changes that occur during charging and discharging. An innovative system aimed at resolving this issue is being developed. This system incorporates a thin Si film bonded to an elastomeric substrate which is intended to provide the desired stress relief. Non-linear finite element simulations have shown that a significant amount of deformation can be accommodated until a critical threshold of Li concentration is reached; beyond which buckling is induced and a wavy structure appears. When compared to a similar system using rigid substrates where no buckling occurs, the stress is reduced by an order of magnitude, significantly prolonging the life of the Si anode. Thus the stress can be released at high Li-ion diffusion induced strains by buckling the Si thin film. Several aspects of this anode system have been analyzed including studying the effects of charge rate and thin film plasticity, and the results are compared with preliminary empirical measurements to show great promise. This study serves as the basis for a radical resolution to one of the few remaining barriers left in the development of high performing Si based electrodes for Li-ion batteries.
ContributorsShaffer, Joseph (Author) / Jiang, Hanqing (Thesis advisor) / Rajan, Subramaniam D. (Committee member) / Peralta, Pedro (Committee member) / Arizona State University (Publisher)
Created2011