On the Equivalence between the Additive Hypo-Elasto-Plasticity and Multiplicative Hyper-Elasto-Plasticity Models and Adaptive Propagation of Discontinuities

On the Equivalence between the Additive Hypo-Elasto-Plasticity and Multiplicative HyperElasto-Plasticity Models and Adaptive Propagation of Discontinuities Yang Jiao Ductile and brittle failure of solids are closely related to their plastic and fracture behavior, respectively. The two most common energy dissipation mechanisms in solids possess distinct kinematic characteristics, i.e. large strain and discontinuous displacement, both of which pose challenges to reliable, efficient numerical simulation of material failure in engineering structures. This dissertation addresses the reliability and efficiency issues associated with the kinematic characteristics of plasticity and fracture. At first, studies are conducted to understand the relation between two well recognized large strain plasticity models that enjoy widespread popularity in numerical simulation of plastic behavior of solids. These two models, termed the additive hypo-elasto-plasticity and multiplicative hyper-elasto-plasticity models, respectively, are regarded as two distinct strategies for extending the classical infinitesimal deformation plasticity theory into the large strain regime. One of the most recent variants of the additive models, which features the logarithmic stress rate, is shown to give rise to nonphysical energy dissipation during elastic unloading. A simple modification to the logarithmic stress rate is accordingly made to resolve such a physical inconsistency. This results in the additive hypo-elasto-plasticity models based on the kinetic logarithmic stress rate in which energy dissipation-free elastic response is produced whenever plastic flow is absent. It is then proved that for isotropic materials the multiplicative hyper-elasto-plasticity models coincide with the additive ones if a newly discovered objective stress rate is adopted. Such an objective stress rate, termed the modified kinetic logarithmic rate, reduces to the kinetic logarithmic rate in the absence of strain-induced anisotropy which is characterized as kinematic hardening in the present dissertation. In the second part of the dissertation, the computational complexity of finite element analysis of the onset and propagation of interface cracks in layered materials is addressed. The study is conducted in the context of laminated composites in which interface fracture (delamination) is a dominant failure mode. In order to eliminate the complexities of remeshing for constant initiation and propagation of delamination, two hierarchical approaches, the extended finite element method (XFEM) and the s-version of the finite element method (s-method) are studied in terms of their effectiveness in representing displacement discontinuity across delaminated interfaces. With one single layer of 20-node serendipity solid elements resolving delamination-free response of the layered materials, it is proved that the delamination representations based on the s-method and the XFEM result in the same discretization space as the conventional non-hierarchical ply-by-ply approach which employs one layer of solid elements for each ply as well as double nodes on delaminated interfaces. Delamination indicators based on the s-method representation of delamination are then proposed to detect the onset and propagation of delamination. An adaptive methodology is accordingly developed in which the s-method displacement field enrichment for delamination is adaptively added to interface areas with high likelihood of delamination. Numerical examples show that the computational cost of the adaptive s-method is significantly lower than that incurred by the conventional ply-by-ply approach despite the fact that the two approaches produce practically identical results.