Explore the vast range of titles available.

A stabilized enhanced strain finite element procedure for poromechanics is fully integrated with an elasto-plastic cap model to simulate the hydro-mechanical interactions of fluid-infiltrating porous rocks with associative and non-associative plastic flow. We present a quantitative analysis on how macroscopic plastic volumetric response caused by pore collapse and grain rearrangement affects the seepage of pore fluid, and vice versa. Results of finite element simulations imply that the dissipation of excess pore pressure may significantly affect the stress path and thus alter the volumetric plastic responses.

The Sandia Fracture Challenges provide a forum for the mechanics community to assess its ability to predict ductile fracture through a blind, round-robin format where mechanicians are challenged to predict the deformation and failure of an arbitrary geometry given experimental calibration data. The Third Challenge (SFC3) required participants to predict fracture in an additively manufactured (AM) 316L stainless steel bar containing through holes and internal cavities that could not have been conventionally machined. The volunteer participants were provided extensive data including tension and notched tensions tests of 316L specimens built on the same build-plate as the Challenge geometry, micro-CT scans of the Challenge specimens and geometric measurements of the feature based on the scans, electron backscatter diffraction (EBSD) information on grain texture, and post-test fractography of the calibration specimens. Surprisingly, the global behavior of the SFC3 geometry specimens had modest variability despite being made of AM metal, with all of the SFC3 geometry specimens failing under the same failure mode. This is attributed to the large stress concentrations from the holes overwhelming the stochastic local influence of the AM voids and surface roughness. The teams were asked to predict a number of quantities of interest in the response based on global and local measures that were compared to experimental data, based partly on Digital Image Correlation (DIC) measurements of surface displacements and strains, including predictions of variability in the resulting fracture response, as the basis for assessment of the predictive capabilities of the modeling and simulation strategies. Twenty-one teams submitted predictions obtained from a variety of methods: the finite element method (FEM) or the mesh-free, peridynamic method; solvers with explicit time integration, implicit time integration, or quasi-statics; fracture methods including element deletion, peridynamics with bond damage, XFEM, damage (stiffness degradation), and adaptive remeshing. These predictions utilized many different material models: plasticity models including J2 plasticity or Hill yield with isotropic hardening, mixed Swift-Voce hardening, kinematic hardening, or custom hardening curves; fracture criteria including GTN model, Hosford-Coulomb, triaxiality-dependent strain, critical fracture energy, damage-based model, critical void volume fraction, and Johnson-Cook model; and damage evolution models including damage accumulation and evolution, crack band model, fracture energy, displacement value threshold, incremental stress triaxiality, Cocks-Ashby void growth, and void nucleation, growth, and coalescence. Teams used various combinations of calibration data from tensile specimens, the notched tensile specimens, and literature data. A detailed comparison of results based of these different methods is presented in this paper to suggest a set of best practices for modeling ductile fracture in situations like the SFC3 AM-material problem. All blind predictions identified the nominal crack path and initiation location correctly. The SFC3 participants generally fared better in their global predictions of deformation and failure than the participants in the previous Challenges, suggesting the relative maturity of the models used and adoption of best practices from previous Challenges. This paper provides detailed analyses of the results, including discussion of the utility of the provided data, challenges of the experimental-numerical comparison, defects in the AM material, and human factors.

We propose a variational procedure for the recovery of internal variables, in effect extending them from integration points to the entire domain. The objective is to perform the recovery with minimum error and at the same time guarantee that the internal variables remain in their admissible spaces. The minimization of the error is achieved by a three-field finite element formulation. The fields in the formulation are the deformation mapping, the target or mapped internal variables and a Lagrange multiplier that enforces the equality between the source and target internal variables. This formulation leads to an L 2 projection that minimizes the distance between the source and target internal variables as measured in the L 2 norm of the internal variable space. To ensure that the target internal variables remain in their original space, their interpolation is performed by recourse to Lie groups, which allows for direct polynomial interpolation of the corresponding Lie algebras by means of the logarithmic map. Once the Lie algebras are interpolated, the mapped variables are recovered by the exponential map, thus guaranteeing that they remain in the appropriate space.

The third Sandia Fracture Challenge highlighted the geometric and material uncertainties introduced by modern additive manufacturing techniques. Tasked with the challenge of predicting failure of a complex additively-manufactured geometry made of 316L stainless steel, we combined a rigorous material calibration scheme with a number of statistical assessments of problem uncertainties. Specifically, we used optimization techniques to calibrate a rate-dependent and anisotropic Hill plasticity model to represent material deformation coupled with a damage model driven by void growth and nucleation. Through targeted simulation studies we assessed the influence of internal voids and surface flaws on the specimens of interest in the challenge which guided our material modeling choices. Employing the Kolmogorov–Smirnov test statistic, we developed a representative suite of simulations to account for the geometric variability of test specimens and the variability introduced by material parameter uncertainty. This approach allowed the team to successfully predict the failure mode of the experimental test population as well as the global response with a high degree of accuracy.

This dissertation presents a formulation of incompatibility based strain gradient plasticity utilizing discontinuous Galerkin finite element methods. Foundations of the classical theory of plasticity are laid out including a discussion of computational implementation and a series of numerical examples. A gradient plasticity constitutive theory is developed based on micromechanical arguments emanating from dislocation theory generalized into a continuum, tensorial treatment. The variational statement of the gradient plasticity equations utilizes concepts from discontinuous Galerkin methods to account for the continuity requirements dictated by the theory. Implementation of the method and solution procedures are discussed and numerical examples are presented showing the well established size effect for materials at small scales and mesh independence for a boundary value problem exhibiting localization due to softening.

The advent of fabrication techniques such as additive manufacturing has focused attention on the considerable variability of material response due to defects and other microstructural aspects. This variability motivates the development of an enhanced design methodology that incorporates inherent material variability to provide robust predictions of performance. In this work, we develop plasticity models capable of representing the distribution of mechanical responses observed in experiments using traditional plasticity models of the mean response and recently developed uncertainty quantification (UQ) techniques. We demonstrate that the new method provides predictive realizations that are superior to more traditional ones, and how these UQ techniques can be used in model selection and assessing the quality of calibrated physical parameters.

We use machine learning (ML) to infer stress and plastic flow rules using data from repre- sentative polycrystalline simulations. In particular, we use so-called deep (multilayer) neural networks (NN) to represent the two response functions. The ML process does not choose ap- propriate inputs or outputs, rather it is trained on selected inputs and output. Likewise, its discrimination of features is crucially connected to the chosen input-output map. Hence, we draw upon classical constitutive modeling to select inputs and enforce well-accepted symmetries and other properties. With these developments, we enable rapid model building in real-time with experiments, and guide data collection and feature discovery.