Matching Items (15)
Filtering by
- All Subjects: Mechanics
- Creators: Jiang, Hanqing
Description
Pb-free solders are used as interconnects in various levels of micro-electronic packaging. Reliability of these interconnects is very critical for the performance of the package. One of the main factors affecting the reliability of solder joints is the presence of porosity which is introduced during processing of the joints. In this thesis, the effect of such porosity on the deformation behavior and eventual failure of the joints is studied using Finite Element (FE) modeling technique. A 3D model obtained by reconstruction of x-ray tomographic image data is used as input for FE analysis to simulate shear deformation and eventual failure of the joint using ductile damage model. The modeling was done in ABAQUS (v 6.10). The FE model predictions are validated with experimental results by comparing the deformation of the pores and the crack path as predicted by the model with the experimentally observed deformation and failure pattern. To understand the influence of size, shape, and distribution of pores on the mechanical behavior of the joint four different solder joints with varying degrees of porosity are modeled using the validated FE model. The validation technique mentioned above enables comparison of the simulated and actual deformation only. A more robust way of validating the FE model would be to compare the strain distribution in the joint as predicted by the model and as observed experimentally. In this study, to enable visualization of the experimental strain for the 3D microstructure obtained from tomography, a three dimensional digital image correlation (3D DIC) code has been implemented in MATLAB (MathWorks Inc). This developed 3D DIC code can be used as another tool to verify the numerical model predictions. The capability of the developed code in measuring local displacement and strain is demonstrated by considering a test case.
ContributorsJakkali, Vaidehi (Author) / Chawla, Nikhilesh K (Thesis advisor) / Jiang, Hanqing (Committee member) / Solanki, Kiran (Committee member) / Arizona State University (Publisher)
Created2011
Description
Advanced composites are being widely used in aerospace applications due to their high stiffness, strength and energy absorption capabilities. However, the assurance of structural reliability is a critical issue because a damage event will compromise the integrity of composite structures and lead to ultimate failure. In this dissertation a novel homogenization based multiscale modeling framework using semi-analytical micromechanics is presented to simulate the response of textile composites. The novelty of this approach lies in the three scale homogenization/localization framework bridging between the constituent (micro), the fiber tow scale (meso), weave scale (macro), and the global response. The multiscale framework, named Multiscale Generalized Method of Cells (MSGMC), continuously bridges between the micro to the global scale as opposed to approaches that are top-down and bottom-up. This framework is fully generalized and capable of modeling several different weave and braids without reformulation. Particular emphasis in this dissertation is placed on modeling the nonlinearity and failure of both polymer matrix and ceramic matrix composites.
ContributorsLiu, Guang (Author) / Chattopadhyay, Aditi (Thesis advisor) / Mignolet, Marc (Committee member) / Jiang, Hanqing (Committee member) / Li, Jian (Committee member) / Rajadas, John (Committee member) / Arizona State University (Publisher)
Created2011
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
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
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
Description
The study of deflagration to detonation transition (DDT) in explosives is of prime importance with regards to insensitive munitions (IM). Critical damage owing to thermal or shock stimuli could translate to significant loss of life and material. The present study models detonation and deflagration of a commonly used granular explosive: cyclotetramethylene-tetranitramine, HMX. A robust literature review is followed by computational modeling of gas gun and DDT tube test data using the Sandia National Lab three-dimensional multi-material Eulerian hydrocode CTH. This dissertation proposes new computational practices and models that aid in predicting shock stimulus IM response. CTH was first used to model experimental data sets of DDT tubes from both Naval Surface Weapons Center and Los Alamos National Laboratory which were initiated by pyrogenic material and a piston, respectively. Analytical verification was performed, where possible, for detonation via empirical based equations at the Chapman Jouguet state with errors below 2.1%, and deflagration via pressure dependent burn rate equations. CTH simulations include inert, history variable reactive burn and Arrhenius models. The results are in excellent agreement with published HMX detonation velocities. Novel additions include accurate simulation of the pyrogenic material BKNO3 and the inclusion of porosity in energetic materials. The treatment of compaction is especially important in modeling precursory hotspots, caused by hydrodynamic collapse of void regions or grain interactions, prior to DDT of granular explosives. The CTH compaction model of HMX was verified within 11% error via a five pronged validation approach using gas gun data and employed use of a newly generated set of P-α parameters for granular HMX in a Mie-Gruneisen Equation of State. Next, the additions of compaction were extended to a volumetric surface burning model of HMX and compare well to a set of empirical burn rates. Lastly, the compendium of detonation and deflagration models was applied to the aforementioned DDT tubes and demonstrate working functionalities of all models, albeit at the expense of significant computational resources. A robust hydrocode methodology is proposed to make use of the deflagration, compaction and detonation models as a means to predict IM response to shock stimulus of granular explosive materials.
ContributorsMahon, Kelly Susan (Author) / Lee, Taewoo (Thesis advisor) / Herrmann, Marcus (Committee member) / Chen, Kangping (Committee member) / Jiao, Yang (Committee member) / Huang, Huei-Ping (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
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
Polymer matrix composites (PMCs) are attractive structural materials due to their high stiffness to low weight ratio. However, unidirectional PMCs have low shear strength and failure can occur along kink bands that develop on compression due to plastic microbuckling that carry strains large enough to induce nonlinear matrix deformation. Reviewing the literature, a large fraction of the existing work is for uniaxial compression, and the effects of stress gradients, such as those present during bending, have not been as well explored, and these effects are bound to make difference in terms of kink band nucleation and growth. Furthermore, reports on experimental measurements of strain fields leading to and developing inside these bands in the presence of stress gradients are also scarce and need to be addressed to gain a full understanding of their behavior when UDCs are used under bending and other spatially complex stress states.
In a light to bridge the aforementioned gaps, the primary focus of this work is to understand mechanisms for kink band evolution under an influence of stress-gradients induced during bending. Digital image correlation (DIC) is used to measure strains inside and around the kink bands during 3-point bending of samples with 0°/90° stacking made of Ultra-High Molecular Weight Polyethylene Fibers. Measurements indicate bands nucleate at the compression side and propagate into the sample carrying a mixture of large shear and normal strains (~33%), while also decreasing its bending stiffness. Failure was produced by a combination of plastic microbuckling and axial splitting. The microstructure of the kink bands was studied and used in a microstructurally explicit finite element model (FEM) to analyze stresses and strains at ply level in the samples during kink band evolution, using cohesive zone elements to represent the interfaces between plies. Cohesive element properties were deduced by a combination of delamination, fracture and three-point bending tests used to calibrate the FEMs. Modeling results show that the band morphology is sensitive to the shear and opening properties of the interfaces between the plies.
In a light to bridge the aforementioned gaps, the primary focus of this work is to understand mechanisms for kink band evolution under an influence of stress-gradients induced during bending. Digital image correlation (DIC) is used to measure strains inside and around the kink bands during 3-point bending of samples with 0°/90° stacking made of Ultra-High Molecular Weight Polyethylene Fibers. Measurements indicate bands nucleate at the compression side and propagate into the sample carrying a mixture of large shear and normal strains (~33%), while also decreasing its bending stiffness. Failure was produced by a combination of plastic microbuckling and axial splitting. The microstructure of the kink bands was studied and used in a microstructurally explicit finite element model (FEM) to analyze stresses and strains at ply level in the samples during kink band evolution, using cohesive zone elements to represent the interfaces between plies. Cohesive element properties were deduced by a combination of delamination, fracture and three-point bending tests used to calibrate the FEMs. Modeling results show that the band morphology is sensitive to the shear and opening properties of the interfaces between the plies.
ContributorsPatel, Jay K (Author) / Peralta, Pedro D (Thesis advisor) / Oswald, Jay (Committee member) / Jiang, Hanqing (Committee member) / Solanki, Kiran (Committee member) / Ayyar, Adarsh (Committee member) / Arizona State University (Publisher)
Created2016
Description
A previously developed small time scale fatigue crack growth model is improved, modified and extended with an emphasis on creating the simplest models that maintain the desired level of accuracy for a variety of materials. The model provides a means of estimating load sequence effects by continuously updating the crack opening stress every cycle, in a simplified manner. One of the significant phenomena of the crack opening stress under negative stress ratio is the residual tensile stress induced by the applied compressive stress. A modified coefficient is introduced to determine the extent to which residual stress impact the crack closure and is observed to vary for different materials. Several other literature models for crack closure under constant loading are also reviewed and compared with the proposed model. The modified model is then shown to predict several sets of published test results under constant loading for a variety of materials.
The crack opening stress is formalized as a function of the plastic zone sizes at the crack tip and the current crack length, which provided a means of approximation, accounting for both acceleration and retardation effects in a simplified manner. A sensitivity parameter is introduced to modify the enlarged plastic zone due to overload, to better fit the delay cycles with the test data and is observed to vary for different materials. Furthermore, the interaction effect induced by the combination of overload and underload sequence is modeled by depleting the compressive plastic zone due to an overload with the tensile plastic zone due to an underload. A qualitative analysis showed the simulation capacity of the small time scale model under different load types. A good agreement between prediction and test data for several irregular load types proved the applicability of the small time scale model under variable amplitude loading.
The crack opening stress is formalized as a function of the plastic zone sizes at the crack tip and the current crack length, which provided a means of approximation, accounting for both acceleration and retardation effects in a simplified manner. A sensitivity parameter is introduced to modify the enlarged plastic zone due to overload, to better fit the delay cycles with the test data and is observed to vary for different materials. Furthermore, the interaction effect induced by the combination of overload and underload sequence is modeled by depleting the compressive plastic zone due to an overload with the tensile plastic zone due to an underload. A qualitative analysis showed the simulation capacity of the small time scale model under different load types. A good agreement between prediction and test data for several irregular load types proved the applicability of the small time scale model under variable amplitude loading.
ContributorsVenkatesan, Karthik Rajan (Author) / Liu, Yongming (Thesis advisor) / Oswald, Jay (Committee member) / Jiang, Hanqing (Committee member) / Arizona State University (Publisher)
Created2016