Matching Items (31)
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
An integrated experimental and numerical investigation for laser-generated optoacoustic wave propagation in structural materials is performed. First, a multi-physics simulation model is proposed to simulate the pulsed laser as a point heat source which hits the surface of an aluminum sheet. The pulsed laser source can generate a localized heating on the surface of the plate and induce an in-plane stress wave. ANSYS – a finite element analysis software – is used to build the 3D model and a coupled thermal-mechanical simulation is performed in which the heat flux is determined by an empirical laser-heat conversion relationship. The displacement and stress field-histories are obtained to get the time of arrival and wave propagation speed of the stress wave. The effect of an added point mass is investigated in detail to observe the local material perturbation and remote wave signals. Following this, the experimental investigation of optoacoustic wave is also performed. A new experimental setup and control is developed and assembled in-house. Various laser firing parameters are investigated experimentally and the optimal combination is used for the experimental testing. Matrix design for different testing conditions is also proposed to include the effect of wave path, sampling procedure, and local point mass on the optoacoustic wave propagation. The developed numerical simulation results are validated with experimental observations. It is shown that the proposed experimental setup can offer a potential fast scanning method for damage detection (local property change) for plate-like structural component.
ContributorsLiu, Chen (Author) / Liu, Yongming (Thesis advisor) / Wang, Liping (Committee member) / Jiao, Yang (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
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
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
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
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