Matching Items (12)
Filtering by
- All Subjects: Finite element method
- Creators: Jiang, Hanqing
- Genre: Academic theses
Description
Woven fabric composite materials are widely used in the construction of aircraft engine fan containment systems, mostly due to their high strength to weight ratios and ease of implementation. The development of a predictive model for fan blade containment would provide great benefit to engine manufactures in shortened development cycle time, less risk in certification and fewer dollars lost to redesign/recertification cycles. A mechanistic user-defined material model subroutine has been developed at Arizona State University (ASU) that captures the behavioral response of these fabrics, namely Kevlar® 49, under ballistic loading. Previously developed finite element models used to validate the consistency of this material model neglected the effects of the physical constraints imposed on the test setup during ballistic testing performed at NASA Glenn Research Center (NASA GRC). Part of this research was to explore the effects of these boundary conditions on the results of the numerical simulations. These effects were found to be negligible in most instances. Other material models for woven fabrics are available in the LS-DYNA finite element code. One of these models, MAT234: MAT_VISCOELASTIC_LOOSE_FABRIC (Ivanov & Tabiei, 2004) was studied and implemented in the finite element simulations of ballistic testing associated with the FAA ASU research. The results from these models are compared to results obtained from the ASU UMAT as part of this research. The results indicate an underestimation in the energy absorption characteristics of the Kevlar 49 fabric containment systems. More investigation needs to be performed in the implementation of MAT234 for Kevlar 49 fabric. Static penetrator testing of Kevlar® 49 fabric was performed at ASU in conjunction with this research. These experiments are designed to mimic the type of loading experienced during fan blade out events. The resulting experimental strains were measured using a non-contact optical strain measurement system (ARAMIS).
ContributorsFein, Jonathan (Author) / Rajan, Subramaniam D. (Thesis advisor) / Mobasher, Barzin (Committee member) / Jiang, Hanqing (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
Concrete columns constitute the fundamental supports of buildings, bridges, and various other infrastructures, and their failure could lead to the collapse of the entire structure. As such, great effort goes into improving the fire resistance of such columns. In a time sensitive fire situation, a delay in the failure of critical load bearing structures can lead to an increase in time allowed for the evacuation of occupants, recovery of property, and access to the fire. Much work has been done in improving the structural performance of concrete including reducing column sizes and providing a safer structure. As a result, high-strength (HS) concrete has been developed to fulfill the needs of such improvements. HS concrete varies from normal-strength (NS) concrete in that it has a higher stiffness, lower permeability and larger durability. This, unfortunately, has resulted in poor performance under fire. The lower permeability allows for water vapor to build up causing HS concrete to suffer from explosive spalling under rapid heating. In addition, the coefficient of thermal expansion (CTE) of HS concrete is lower than that of NS concrete. In this study, the effects of introducing a region of crumb rubber concrete into a steel-reinforced concrete column were analyzed. The inclusion of crumb rubber concrete into a column will greatly increase the thermal resistivity of the overall column, leading to a reduction in core temperature as well as the rate at which the column is heated. Different cases were analyzed while varying the positioning of the crumb-rubber region to characterize the effect of position on the improvement of fire resistance. Computer simulated finite element analysis was used to calculate the temperature and strain distribution with time across the column's cross-sectional area with specific interest in the steel - concrete region. Of the several cases which were investigated, it was found that the improvement of time before failure ranged between 32 to 45 minutes.
ContributorsZiadeh, Bassam Mohammed (Author) / Phelan, Patrick (Thesis advisor) / Kaloush, Kamil (Thesis advisor) / Jiang, Hanqing (Committee member) / Arizona State University (Publisher)
Created2011
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
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
Origami and kirigami, the technique of generating three-dimensional (3D) structures from two-dimensional (2D) flat sheets, are now more and more involved in scientific and engineering fields. Therefore, the development of tools for their theoretical analysis becomes more and more important. Since much effort was paid on calculations based on pure mathematical consideration and only limited effort has been paid to include mechanical properties, the goal of my research is developing a method to analyze the mechanical behavior of origami and kirigami based structures. Mechanical characteristics, including nonlocal effect and fracture of the structures, as well as elasticity and plasticity of materials are studied. For calculation of relative simple structures and building of structures’ constitutive relations, analytical approaches were used. For more complex structures, finite element analysis (FEA), which is commonly applied as a numerical method for the analysis of solid structures, was utilized. The general study approach is not necessarily related to characteristic size of model. I believe the scale-independent method described here will pave a new way to understand the mechanical response of a variety of origami and kirigami based structures under given mechanical loading.
ContributorsLv, Cheng (Author) / Jiang, Hanqing (Thesis advisor) / Yu, Hongbin (Committee member) / Wang, Liping (Committee member) / Mignolet, Marc (Committee member) / Hildreth, Owen (Committee member) / Arizona State University (Publisher)
Created2016
Description
In this paper, at first, analytical formulation of J-integral for a non-local particle model (VCPM) using atomic scale finite element method is proposed for fracture analysis of 2D solids. A brief review of classical continuum-based J-integral and anon-local lattice particle method is given first. Following this, detailed derivation for the J-integral in discrete particle system is given using the energy equivalence and stress-tensor mapping between the continuum mechanics and lattice-particle system.With the help of atomistic finite element method, the J-integral is expressed as a summation of the corresponding terms in the particle system.
Secondly, a coupling algorithm between a non-local particle method (VCPM) and the classical finite element method (FEM) is discussed to gain the advantages of both methods for fracture analysis in large structures. In this algorithm, the discrete VCPM particle and the continuum FEM domains are solved within a unified theoretical framework. A transitional element technology is developed to smoothly link the 10-particles element with the traditional FEM elements to guaranty the continuity and consistency at the coupling interface. An explicit algorithm for static simulation is developed.
Finally, numerical examples are illustrated for the accuracy, convergence, and path-independence of the derived J-integral formulation. Discussions on the comparison with alternative estimation methods and potential application for fracture simulation are given. The accuracy and efficiency of the coupling algorithm are tested by several benchmark problems such as static crack simulation.
Secondly, a coupling algorithm between a non-local particle method (VCPM) and the classical finite element method (FEM) is discussed to gain the advantages of both methods for fracture analysis in large structures. In this algorithm, the discrete VCPM particle and the continuum FEM domains are solved within a unified theoretical framework. A transitional element technology is developed to smoothly link the 10-particles element with the traditional FEM elements to guaranty the continuity and consistency at the coupling interface. An explicit algorithm for static simulation is developed.
Finally, numerical examples are illustrated for the accuracy, convergence, and path-independence of the derived J-integral formulation. Discussions on the comparison with alternative estimation methods and potential application for fracture simulation are given. The accuracy and efficiency of the coupling algorithm are tested by several benchmark problems such as static crack simulation.
ContributorsZope, Jayesh (Author) / Liu, Yongming (Thesis advisor) / Oswald, Jay (Committee member) / Jiang, Hanqing (Committee member) / Arizona State University (Publisher)
Created2016
Description
Cohesive zone model is one of the most widely used model for fracture analysis, but still remains open ended field for research. The earlier works using the cohesive zone model and Extended finite element analysis (XFEM) have been briefly introduced followed by an elaborate elucidation of the same concepts. Cohesive zone model in conjugation with XFEM is used for analysis in static condition in order to check its applicability in failure analysis. A real time setup of pipeline failure due to impingement is analyzed along with a detailed parametric study to understand the influence of the prominent design variable. After verifying its good applicability, a creep model is built for analysis where the cohesive zone model with XFEM is used for a time dependent creep loading. The challenge in this simulation was to achieve coupled behavior of cracks initiation and propagation along with creep loading. By using Design of Experiment, the results from numerical simulation were used to build an equation for life prediction for creep loading condition. The work was further extended to account for fatigue damage accumulation for high cycle fatigue loading in cohesive elements. A model was conceived to account for damage due to fatigue loading along within cohesive zone model for cohesive elements in ABAQUS simulation software. The model was verified by comparing numerical modelling of Double cantilever beam under high cycle fatigue loading and experiment results from literature. The model was also applied to a major industrial problem of blistering in Cured-In-Plane liner pipelines and a demonstration of its failure is shown. In conclusion, various models built on cohesive zone to address static and time dependent loading with real time scenarios and future scope of work in this field is discussed.
ContributorsChandrasekhar, Vishal (Author) / Liu, Yongming (Thesis advisor) / Oswald, Jay (Committee member) / Jiang, Hanqing (Committee member) / Arizona State University (Publisher)
Created2016
Description
A method for modelling the interactions of dislocations with inclusions has been developed to analyse toughening mechanisms in alloys. This method is different from the superposition method in that infinite domain solutions and image stress fields are not superimposed. The method is based on the extended finite element method (XFEM) in which the dislocations are modelled according to the Volterra dislocation model. Interior discontinuities are introduced across dislocation glide planes using enrichment functions and the resulting boundary value problem is solved through the standard finite element variational approach. The level set method is used to describe the geometry of the dislocation glide planes without any explicit treatment of the interface geometry which provides a convenient and an appealing means for describing the dislocation. A method for estimating the Peach-Koehler force by the domain form of J-integral is considered. The convergence and accuracy of the method are studied for an edge dislocation interacting with a free surface where analytical solutions are available. The force converges to the exact solution at an optimal rate for linear finite elements. The applicability of the method to dislocation interactions with inclusions is illustrated with a system of Aluminium matrix containing Aluminium-copper precipitates. The effect of size, shape and orientation of the inclusions on an edge dislocation for a difference in stiffness and coefficient of thermal expansion of the inclusions and matrix is considered. The force on the dislocation due to a hard inclusion increased by 8% in approaching the sharp corners of a square inclusion than a circular inclusion of equal area. The dislocation experienced 24% more force in moving towards the edges of a square shaped inclusion than towards its centre. When the areas of the inclusions were halved, 30% less force was exerted on the dislocation. This method was used to analyse interfaces with mismatch strains. Introducing eigenstrains equal to 0.004 to the elastic mismatch increased the force by 15 times for a circular inclusion. The energy needed to move an edge dislocation through a domain filled with circular inclusions is 4% more than that needed for a domain with square shaped inclusions.
ContributorsVeeresh, Pawan (Author) / Oswald, Jay (Thesis advisor) / Jiang, Hanqing (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