Matching Items (12)
Filtering by
- All Subjects: Mechanics
- Creators: Chawla, Nikhilesh
- Creators: Rajan, Subramaniam D.
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
As one of the most promising materials for high capacity electrode in next generation of lithium ion batteries, silicon has attracted a great deal of attention in recent years. Advanced characterization techniques and atomic simulations helped to depict that the lithiation/delithiation of silicon electrode involves processes including large volume change (anisotropic for the initial lithiation of crystal silicon), plastic flow or softening of material dependent on composition, electrochemically driven phase transformation between solid states, anisotropic or isotropic migration of atomic sharp interface, and mass diffusion of lithium atoms. Motivated by the promising prospect of the application and underlying interesting physics, mechanics coupled with multi-physics of silicon electrodes in lithium ion batteries is studied in this dissertation. For silicon electrodes with large size, diffusion controlled kinetics is assumed, and the coupled large deformation and mass transportation is studied. For crystal silicon with small size, interface controlled kinetics is assumed, and anisotropic interface reaction is studied, with a geometry design principle proposed. As a preliminary experimental validation, enhanced lithiation and fracture behavior of silicon pillars via atomic layer coatings and geometry design is studied, with results supporting the geometry design principle we proposed based on our simulations. Through the work documented here, a consistent description and understanding of the behavior of silicon electrode is given at continuum level and some insights for the future development of the silicon electrode are provided.
ContributorsAn, Yonghao (Author) / Jiang, Hanqing (Thesis advisor) / Chawla, Nikhilesh (Committee member) / Phelan, Patrick (Committee member) / Wang, Yinming (Committee member) / Yu, Hongyu (Committee member) / Arizona State University (Publisher)
Created2014
Description
Pb-free solder joints are commonly used as interconnects in semiconductor packaging. One of the major defects affecting the mechanical performance of solder joints are reflow pores that form during processing. These pores exhibit significant variability in size and distribution, and understanding the effects of pore geometry on failure is an important reliability concern. In this thesis, the pore microstructures of solder joint samples and the localized plastic deformation around individual pores was characterized in 3D using lab scale X-ray Microtomography. To observe the deformation of a solder joint in 3D, a solder joint was imaged with Microtomography after reflow and then deformed in shear in several loading steps with additional tomography data taken between each. The 3D tomography datasets were then segmented using the 3D Livewire technique into regions corresponding to solder and pores, and used to generate 3D models of the joint at each strain value using Mimics software. The extent of deformation of individual pores in the joint as a function of strain was quantified using sphericity measurements, and correlated with the observed cracking in the joint. In addition, the error inherent in the data acquisition and 3D modeling process was also quantified. The progression of damage observed with X-ray Microtomography was then used to validate the deformation and failure predicted by a Finite Element (FE) simulation. The FE model was based on the as-reflowed tomography data, and incorporated a ductile damage failure model to simulate fracture. Using the measured sphericity change and cracking information obtained from the tomography data, the FE model is shown to correctly capture the broad plastic deformation and strain localization seen in the actual joint, as well as the crack propagation. Lastly, Digital Image Correlation was investigated as a method of obtaining improved local strain measurements in 3D. This technique measures the displacement of the inherent microstructural features of the joint, and can give localized strain measurements that can be directly comparable to that predicted by modeling. The technique is demonstrated in 2D on Pb-Sn solder, and example 3D data is presented for future analysis.
ContributorsPadilla, Erick (Author) / Chawla, Nikhilesh (Thesis advisor) / Alford, Terry (Committee member) / Krause, Stephen (Committee member) / Arizona State University (Publisher)
Created2012
Description
The objective of this research is to develop robust, accurate, and adaptive algorithms in the framework of the extended finite element method (XFEM) for fracture analysis of highly heterogeneous materials with complex internal geometries. A key contribution of this work is the creation of novel methods designed to automate the incorporation of high-resolution data, e.g. from X-ray tomography, that can be used to better interpret the enormous volume of data generated in modern in-situ experimental testing. Thus new algorithms were developed for automating analysis of complex microstructures characterized by segmented tomographic images.
A centrality-based geometry segmentation algorithm was developed to accurately identify discrete inclusions and particles in composite materials where limitations in imaging resolution leads to spurious connections between particles in close contact.To allow for this algorithm to successfully segment geometry independently of particle size and shape, a relative centrality metric was defined to allow for a threshold centrality criterion for removal of voxels that spuriously connect distinct geometries.
To automate incorporation of microstructural information from high-resolution images, two methods were developed that initialize signed distance fields on adaptively-refined finite element meshes. The first method utilizes a level set evolution equation that is directly solved on the finite element mesh through Galerkins method. The evolution equation is formulated to produce a signed distance field that matches geometry defined by a set of voxels segmented from tomographic images. The method achieves optimal convergence for the order of elements used. In a second approach, the fast marching method is employed to initialize a distance field on a uniform grid which is then projected by least squares onto a finite element mesh. This latter approach is shown to be superior in speed and accuracy.
Lastly, extended finite element method simulations are performed for the analysis of particle fracture in metal matrix composites with realistic particle geometries initialized from X-ray tomographic data. In the simulations, particles fracture probabilistically through a Weibull strength distribution. The model is verified through comparisons with the experimentally-measured stress-strain response of the material as well as analysis of the fracture. Further, simulations are then performed to analyze the effect of mesh sensitivity, the effect of fracture of particles on their neighbors, and the role of a particles shape on its fracture probability.
A centrality-based geometry segmentation algorithm was developed to accurately identify discrete inclusions and particles in composite materials where limitations in imaging resolution leads to spurious connections between particles in close contact.To allow for this algorithm to successfully segment geometry independently of particle size and shape, a relative centrality metric was defined to allow for a threshold centrality criterion for removal of voxels that spuriously connect distinct geometries.
To automate incorporation of microstructural information from high-resolution images, two methods were developed that initialize signed distance fields on adaptively-refined finite element meshes. The first method utilizes a level set evolution equation that is directly solved on the finite element mesh through Galerkins method. The evolution equation is formulated to produce a signed distance field that matches geometry defined by a set of voxels segmented from tomographic images. The method achieves optimal convergence for the order of elements used. In a second approach, the fast marching method is employed to initialize a distance field on a uniform grid which is then projected by least squares onto a finite element mesh. This latter approach is shown to be superior in speed and accuracy.
Lastly, extended finite element method simulations are performed for the analysis of particle fracture in metal matrix composites with realistic particle geometries initialized from X-ray tomographic data. In the simulations, particles fracture probabilistically through a Weibull strength distribution. The model is verified through comparisons with the experimentally-measured stress-strain response of the material as well as analysis of the fracture. Further, simulations are then performed to analyze the effect of mesh sensitivity, the effect of fracture of particles on their neighbors, and the role of a particles shape on its fracture probability.
ContributorsYuan, Rui (Author) / Oswald, Jay (Thesis advisor) / Chawla, Nikhilesh (Committee member) / Liu, Yongming (Committee member) / Solanki, Kiran (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2015
Description
Hydrogen embrittlement (HE) is a phenomenon that affects both the physical and chemical properties of several intrinsically ductile metals. Consequently, understanding the mechanisms behind HE has been of particular interest in both experimental and modeling research. Discrepancies between experimental observations and modeling results have led to various proposals for HE mechanisms. Therefore, to gain insights into HE mechanisms in iron, this dissertation aims to investigate several key issues involving HE such as: a) the incipient crack tip events; b) the cohesive strength of grain boundaries (GBs); c) the dislocation-GB interactions and d) the dislocation mobility.
The crack tip, which presents a preferential trap site for hydrogen segregation, was examined using atomistic methods and the continuum based Rice-Thompson criterion as sufficient concentration of hydrogen can alter the crack tip deformation mechanism. Results suggest that there is a plausible co-existence of the adsorption induced dislocation emission and hydrogen enhanced decohesion mechanisms. In the case of GB-hydrogen interaction, we observed that the segregation of hydrogen along the interface leads to a reduction in cohesive strength resulting in intergranular failure. A methodology was further developed to quantify the role of the GB structure on this behavior.
GBs play a fundamental role in determining the strengthening mechanisms acting as an impediment to the dislocation motion; however, the presence of an unsurmountable barrier for a dislocation can generate slip localization that could further lead to intergranular crack initiation. It was found that the presence of hydrogen increases the strain energy stored within the GB which could lead to a transition in failure mode. Finally, in the case of body centered cubic metals, understanding the complex screw dislocation motion is critical to the development of an accurate continuum description of the plastic behavior. Further, the presence of hydrogen has been shown to drastically alter the plastic deformation, but the precise role of hydrogen is still unclear. Thus, the role of hydrogen on the dislocation mobility was examined using density functional theory and atomistic simulations. Overall, this dissertation provides a novel atomic-scale understanding of the HE mechanism and development of multiscale tools for future endeavors.
The crack tip, which presents a preferential trap site for hydrogen segregation, was examined using atomistic methods and the continuum based Rice-Thompson criterion as sufficient concentration of hydrogen can alter the crack tip deformation mechanism. Results suggest that there is a plausible co-existence of the adsorption induced dislocation emission and hydrogen enhanced decohesion mechanisms. In the case of GB-hydrogen interaction, we observed that the segregation of hydrogen along the interface leads to a reduction in cohesive strength resulting in intergranular failure. A methodology was further developed to quantify the role of the GB structure on this behavior.
GBs play a fundamental role in determining the strengthening mechanisms acting as an impediment to the dislocation motion; however, the presence of an unsurmountable barrier for a dislocation can generate slip localization that could further lead to intergranular crack initiation. It was found that the presence of hydrogen increases the strain energy stored within the GB which could lead to a transition in failure mode. Finally, in the case of body centered cubic metals, understanding the complex screw dislocation motion is critical to the development of an accurate continuum description of the plastic behavior. Further, the presence of hydrogen has been shown to drastically alter the plastic deformation, but the precise role of hydrogen is still unclear. Thus, the role of hydrogen on the dislocation mobility was examined using density functional theory and atomistic simulations. Overall, this dissertation provides a novel atomic-scale understanding of the HE mechanism and development of multiscale tools for future endeavors.
ContributorsAdlakha, Ilaksh (Author) / Solanki, Kiran (Thesis advisor) / Mignolet, Marc (Committee member) / Chawla, Nikhilesh (Committee member) / Jiang, Hanqing (Committee member) / Liu, Yongming (Committee member) / Arizona State University (Publisher)
Created2015
Corrosion and corrosion-fatigue behavior of 7075 aluminum alloys studied by in situ X-ray tomography
Description
7XXX Aluminum alloys have high strength to weight ratio and low cost. They are used in many critical structural applications including automotive and aerospace components. These applications frequently subject the alloys to static and cyclic loading in service. Additionally, the alloys are often subjected to aggressive corrosive environments such as saltwater spray. These chemical and mechanical exposures have been known to cause premature failure in critical applications. Hence, the microstructural behavior of the alloys under combined chemical attack and mechanical loading must be characterized further. Most studies to date have analyzed the microstructure of the 7XXX alloys using two dimensional (2D) techniques. While 2D studies yield valuable insights about the properties of the alloys, they do not provide sufficiently accurate results because the microstructure is three dimensional and hence its response to external stimuli is also three dimensional (3D). Relevant features of the alloys include the grains, subgrains, intermetallic inclusion particles, and intermetallic precipitate particles. The effects of microstructural features on corrosion pitting and corrosion fatigue of aluminum alloys has primarily been studied using 2D techniques such as scanning electron microscopy (SEM) surface analysis along with post-mortem SEM fracture surface analysis to estimate the corrosion pit size and fatigue crack initiation site. These studies often limited the corrosion-fatigue testing to samples in air or specialized solutions, because samples tested in NaCl solution typically have fracture surfaces covered in corrosion product. Recent technological advancements allow observation of the microstructure, corrosion and crack behavior of aluminum alloys in solution in three dimensions over time (4D). In situ synchrotron X-Ray microtomography was used to analyze the corrosion and cracking behavior of the alloy in four dimensions to elucidate crack initiation at corrosion pits for samples of multiple aging conditions and impurity concentrations. Additionally, chemical reactions between the 3.5 wt% NaCl solution and the crack surfaces were quantified by observing the evolution of hydrogen bubbles from the crack. The effects of the impurity particles and age-hardening particles on the corrosion and fatigue properties were examined in 4D.
ContributorsStannard, Tyler (Author) / Chawla, Nikhilesh (Thesis advisor) / Solanki, Kiran N (Committee member) / Goswami, Ramasis (Committee member) / Liu, Yongming (Committee member) / Arizona State University (Publisher)
Created2017
Description
Phase change materials (PCMs) are combined sensible-and-latent thermal energy storage materials that can be used to store and dissipate energy in the form of heat. PCMs incorporated into wall-element systems have been well-studied with respect to energy efficiency of building envelopes. New applications of PCMs in infrastructural concrete, e.g., for mitigating early-age cracking and freeze-and-thaw induced damage, have also been proposed. Hence, the focus of this dissertation is to develop a detailed understanding of the physic-chemical and thermo-mechanical characteristics of cementitious systems and novel coating systems for wall-elements containing PCM. The initial phase of this work assesses the influence of interface properties and inter-inclusion interactions between microencapsulated PCM, macroencapsulated PCM, and the cementitious matrix. The fact that these inclusions within the composites are by themselves heterogeneous, and contain multiple components necessitate careful application of models to predict the thermal properties. The next phase observes the influence of PCM inclusions on the fracture and fatigue behavior of PCM-cementitious composites. The compliant nature of the inclusion creates less variability in the fatigue life for these composites subjected to cyclic loading. The incorporation of small amounts of PCM is found to slightly improve the fracture properties compared to PCM free cementitious composites. Inelastic deformations at the crack-tip in the direction of crack opening are influenced by the microscale PCM inclusions. After initial laboratory characterization of the microstructure and evaluation of the thermo-mechanical performance of these systems, field scale applicability and performance were evaluated. Wireless temperature and strain sensors for smart monitoring were embedded within a conventional portland cement concrete pavement (PCCP) and a thermal control smart concrete pavement (TCSCP) containing PCM. The TCSCP exhibited enhanced thermal performance over multiple heating and cooling cycles. PCCP showed significant shrinkage behavior as a result of compressive strains in the reinforcement that were twice that of the TCSCP. For building applications, novel PCM-composites coatings were developed to improve and extend the thermal efficiency. These coatings demonstrated a delay in temperature by up to four hours and were found to be more cost-effective than traditional building insulating materials.
The results of this work prove the feasibility of PCMs as a temperature-regulating technology. Not only do PCMs reduce and control the temperature within cementitious systems without affecting the rate of early property development but they can also be used as an auto-adaptive technology capable of improving the thermal performance of building envelopes.
The results of this work prove the feasibility of PCMs as a temperature-regulating technology. Not only do PCMs reduce and control the temperature within cementitious systems without affecting the rate of early property development but they can also be used as an auto-adaptive technology capable of improving the thermal performance of building envelopes.
ContributorsAguayo, Matthew Joseph (Author) / Neithalath, Narayanan (Thesis advisor) / Rajan, Subramaniam D. (Committee member) / Mobasher, Barzin (Committee member) / Underwood, Benjamin (Committee member) / Liu, Yongming (Committee member) / Arizona State University (Publisher)
Created2018
Description
In order to verify the dispersive nature of transverse displacement in a beam, a deep understanding of the governing partial differential equation is developed. Using the finite element method and Newmark’s method, along with Fourier transforms and other methods, the aim is to obtain consistent results across each numerical technique. An analytical solution is also analyzed for the Euler-Bernoulli beam in order to gain confidence in the numerical techniques when used for more advance beam theories that do not have a known analytical solution. Three different beam theories are analyzed in this report: The Euler-Bernoulli beam theory, Rayleigh beam theory and Timoshenko beam theory. A comparison of the results show the difference between each theory and the advantages of using a more advanced beam theory for higher frequency vibrations.
ContributorsTschetter, Ryan William (Author) / Hjelmstad, Keith D. (Thesis advisor) / Rajan, Subramaniam D. (Committee member) / Mobasher, Barzin (Committee member) / Arizona State University (Publisher)
Created2016
Description
This study explores the possibility of two matrices containing metallic particulates to act as smart materials by sensing of strain due to the presence of the conducting particles in the matrix. The first matrix is a regular Portland cement-based one while the second is a novel iron-based, carbonated binder developed at ASU. Four different iron replacement percentages by volume (10%, 20%, 30% and 40%) in a Portland cement matrix were selected, whereas the best performing iron carbonate matrix developed was used. Electrical impedance spectroscopy was used to obtain the characteristic Nyquist plot before and after application of flexural load. Electrical circuit models were used to extract the changes in electrical properties under application of load. Strain sensing behavior was evaluated with respect to application of different stress levels and varying replacement levels of the inclusion. A similar approach was used to study the strain sensing capabilities of novel iron carbonate binder. It was observed that the strain sensing efficiency increased with increasing iron percentage and the resistivity increased with increase in load (or applied stress) for both the matrices. It is also found that the iron carbonate binder is more efficient in strain sensing as it had a higher gage factor when compared to the OPC matrix containing metallic inclusions.
Analytical equations (Maxwell) were used to extract frequency dependent electrical conductivity and permittivity of the cement paste (or the host matrix), interface, inclusion (iron) and voids to develop a generic electro-mechanical coupling model to for the strain sensing behavior. COMSOL Multiphysics 5.2a was used as finite element analysis software to develop the model. A MATLAB formulation was used to generate the microstructure with different volume fractions of inclusions. Material properties were assigned (the frequency dependent electrical parameters) and the coupled structural and electrical physics interface in COMSOL was used to model the strain sensing response. The experimental change in resistance matched well with the simulated values, indicating the applicability of the model to predict the strain sensing response of particulate composite systems.
Analytical equations (Maxwell) were used to extract frequency dependent electrical conductivity and permittivity of the cement paste (or the host matrix), interface, inclusion (iron) and voids to develop a generic electro-mechanical coupling model to for the strain sensing behavior. COMSOL Multiphysics 5.2a was used as finite element analysis software to develop the model. A MATLAB formulation was used to generate the microstructure with different volume fractions of inclusions. Material properties were assigned (the frequency dependent electrical parameters) and the coupled structural and electrical physics interface in COMSOL was used to model the strain sensing response. The experimental change in resistance matched well with the simulated values, indicating the applicability of the model to predict the strain sensing response of particulate composite systems.
ContributorsChowdhury, Swaptik (Author) / Neithalath, Narayanan (Thesis advisor) / Rajan, Subramaniam D. (Committee member) / Hoover, Christian G (Committee member) / Arizona State University (Publisher)
Created2017
Description
Ordered buckling of stiff films on elastomeric substrates has many applications in the field of stretchable electronics. Mechanics plays a very important role in such systems. A full three dimensional finite element analysis studying the pattern of wrinkles formed on a stiff film bonded to a compliant substrate under the action of a compressive force has been widely studied. For thin films, this wrinkling pattern is usually sinusoidal, and for wide films the pattern depends on loading conditions. The present study establishes a relationship between the effect of the load applied at an angle to the stiff film. A systematic experimental and analytical study of these systems has been presented in the present study. The study is performed for two different loading conditions, one with the compressive force applied parallel to the film and the other with an angle included between the application of the force and the alignment of the stiff film. A geometric model closely resembling the experimental specimen studied is created and a three dimensional finite element analysis is carried out using ABAQUS (Version 6.7). The objective of the finite element simulations is to validate the results of the experimental study to be corresponding to the minimum total energy of the system. It also helps to establish a relation between the parameters of the buckling profile and the parameters (elastic and dimensional parameters) of the system. Two methods of non-linear analysis namely, the Newton-Raphson method and Arc-Length method are used. It is found that the Arc-Length method is the most cost effective in terms of total simulation time for large models (higher number of elements).The convergence of the results is affected by a variety of factors like the dimensional parameters of the substrate, mesh density of the model, length of the substrate and the film, the angle included. For narrow silicon films the buckling profile is observed to be sinusoidal and perpendicular to the direction of the silicon film. As the angle increases in wider stiff films the buckling profile is seen to transit from being perpendicular to the direction of the film to being perpendicular to the direction of the application of the pre-stress. This study improves and expands the application of the stiff film buckling to an angled loading condition.
ContributorsKondagari, Swathi Sri (Author) / Jiang, Hanqing (Thesis advisor) / Yu, Hongyu (Committee member) / Rajan, Subramaniam D. (Committee member) / Arizona State University (Publisher)
Created2010