There are several engineering applications in which the assumptions of homogenization and scale separation may be violated, in particular, for metallic structures constructed through additive manufacturing. Instead of resorting to direct numerical simulation of the macroscale system with an embedded fine scale, an alternative approach is to use an approximate macroscale constitutive model, but then estimate the model-form error using a posteriori error estimation techniques and subsequently adapt the macroscale model to reduce the error for a given boundary value problem and quantity of interest. Here, we investigate this approach to multiscale analysis in solids with unseparated scales using the example of an additively manufactured metallic structure consisting of a polycrystalline microstructure that is neither periodic nor statistically homogeneous. As a first step to the general nonlinear case, we focus here on linear elasticity in which each grain within the polycrystal is linear elastic but anisotropic.
An improved electrical contact resistance (ECR) model for elastic rough electrode contact is proposed, incorporating the effects of asperity interactions and temperature rise by frictional and joule heating. The analytical simulation results show that the ECR decreases steeply at the beginning of the contact between Al and Cu. However, it becomes stabilized after reaching a specific contact force. It is also found that the longer elapsed sliding contact time, the higher ECR due to the increase in electrical resistivity of electrode materials by the frictional temperature rise at the interface. The effects of surface roughness parameters on ECR are studied through the 32 full-factorial design-of-experiment analysis. Based on the two representative roughness parameters, i.e., root-mean-square (rms) roughness and asperity radius, their individual and coupled effects on the saturated ECR are examined. The saturated ECR increases with the rms roughness for a rough machined surface condition, but it is hardly affected by the asperity radius. On the other hand, the saturated ECR increases with both the rms roughness and the asperity radius under a smooth thin film surface condition.
The 2020 Nonlinear Mechanics and Dynamics (NOMAD) Research Institute was successfully held from June 15 to July 30, 2020. NOMAD brings together participants with diverse technical backgrounds to work in small teams to cultivate new ideas and approaches in engineering mechanics and dynamics research. NOMAD provides an opportunity for researchers – especially early career researchers - to develop lasting collaborations that go beyond what can be established from the limited interactions at their institutions or at annual conferences. A total of 11 students participated in the seven-week long program held virtually due to the COVID-19 health pandemic. The students collaborated on one of four research projects that were developed by various mentors from Sandia National Laboratories, the University of New Mexico, and other academic and research institutions. In addition to the research activities, the students attended weekly technical seminars, various virtual tours, and socialized at virtual gatherings. At the end of the summer, the students gave a final technical presentation on their research findings. Many of the research discoveries made at NOMAD 2020 are published as proceedings at technical conferences and have direct alignment with the critical mission work performed at Sandia.
This study employs nonlinear ultrasonic techniques to track microstructural changes in additively manufactured metals. The second harmonic generation technique based on the transmission of Rayleigh surface waves is used to measure the acoustic nonlinearity parameter, β. Stainless steel specimens are made through three procedures: traditional wrought manufacturing, laser-powder bed fusion, and laser engineered net shaping. The β parameter is measured through successive steps of an annealing heat treatment intended to decrease dislocation density. Dislocation density is known to be sensitive to manufacturing variables. In agreement with fundamental material models for the dislocation-acoustic nonlinearity relationship in the second harmonic generation, β drops in each specimen throughout the heat treatment before recrystallization. Geometrically necessary dislocations (GNDs) are measured from electron back-scatter diffraction as a quantitative indicator of dislocations; average GND density and β are found to have a statistical correlation coefficient of 0.852 showing the sensitivity of β to dislocations in additively manufactured metals. Moreover, β shows an excellent correlation with hardness, which is a measure of the macroscopic effect of dislocations.
We develop a generalized stress inversion technique (or the generalized inversion method) capable of recovering stresses in linear elastic bodies subjected to arbitrary cuts. Specifically, given a set of displacement measurements found experimentally from digital image correlation (DIC), we formulate a stress estimation inverse problem as a partial differential equation-constrained optimization problem. We use gradient-based optimization methods, and we accordingly derive the necessary gradient and Hessian information in a matrix-free form to allow for parallel, large-scale operations. By using a combination of finite elements, DIC, and a matrix-free optimization framework, the generalized inversion method can be used on any arbitrary geometry, provided that the DIC camera can view a sufficient part of the surface. We present numerical simulations and experiments, and we demonstrate that the generalized inversion method can be applied to estimate residual stress.
A stable and nodally integrated meshfree formulation for modeling shock waves in fluids is developed. The reproducing kernel approximation is employed to discretize the conservation equations for compressible flow, and a flux vector splitting approach is applied to allow proper numerical treatments for the advection and pressure parts, respectively, based on the characteristics of each flux term. To capture the essential shock physics in fluids, including the Rankine–Hugoniot jump conditions and the entropy condition, local Riemann enrichment is introduced under the stabilized conforming nodal integration (SCNI) framework. Meanwhile, numerical instabilities associated with the advection flux are eliminated by adopting a modified upwind scheme. To further enhance accuracy, a MUSCL-type method is introduced in conjunction with an oscillation limiter to avoid Gibbs phenomenon and ensure monotonic piecewise linear reconstruction in the smooth region. The present meshfree formulation is free from tunable artificial parameters and is capable of capturing shock and rarefaction waves without over/undershoots. Several numerical examples are analyzed to demonstrate the effectiveness of the proposed MUSCL-SCNI approach in meshfree modeling of complex shock phenomena, including shock diffraction, shock–vortex interaction, and high energy explosion processes.
Meshfree methods for solid mechanics have been in development since the early 1990’s. Initial motivations included alleviation of the burden of mesh creation and the desire to overcome the limitations of traditional mesh-based discretizations for extreme deformation applications. Here, the accuracy and robustness of both meshfree and mesh-based Lagrangian discretizations are compared using manufactured extreme deformation fields. For the meshfree discretizations, both moving least squares and maximum entropy are considered. Quantitative error and convergence results are presented for the best approximation in the H1 norm.
General polyhedral discretizations offer several advantages over classical approaches consisting of standard tetrahedra and hexahedra. These include increased flexibility and robustness in the meshing of geometrically complex domains and higher-quality solutions for both finite element and finite volume schemes. Currently, the use of general polyhedra is hampered by the lack of general-purpose polyhedral meshing algorithms and software. One approach for generating polyhedral meshes is the use of tetrahedral subdivisions and dual-cell aggregation. In this approach, each tetrahedron of an existing tetrahedral mesh is subdivided using one of several subdivision schemes. Polyhedral-dual cells may then be formed and formulated as finite elements with shape functions obtained through the use of generalized barycentric coordinates. We explore the use of dual-cell discretizations for applications in nonlinear solid mechanics using a displacement-based finite element formulation. Verification examples are presented that yield optimal rates of convergence. Accuracy of the methodology is demonstrated via several nonlinear examples that include large deformation and plasticity.
Additive Manufacturing (AM) offers the opportunity to transform design, manufacturing, and qualification with its unique capabilities. AM is a disruptive technology, allowing the capability to simultaneously create part and material while tightly controlling and monitoring the manufacturing process at the voxel level, with the inherent flexibility and agility in printing layer-by-layer. AM enables the possibility of measuring critical material and part parameters during manufacturing, thus changing the way we collect data, assess performance, and accept or qualify parts. It provides an opportunity to shift from the current iterative design-build-test qualification paradigm using traditional manufacturing processes to design-by-predictivity where requirements are addressed concurrently and rapidly. The new qualification paradigm driven by AM provides the opportunity to predict performance probabilistically, to optimally control the manufacturing process, and to implement accelerated cycles of learning. Exploiting these capabilities to realize a new uncertainty quantification-driven qualification that is rapid, flexible, and practical is the focus of this paper.
Residual stress is a common result of manufacturing processes, but it is one that is often overlooked in design and qualification activities. There are many reasons for this oversight, such as lack of observable indicators and difficulty in measurement. Traditional relaxation-based measurement methods use some type of material removal to cause surface displacements, which can then be used to solve for the residual stresses relieved by the removal. While widely used, these methods may offer only individual stress components or may be limited by part or cut geometry requirements. Diffraction-based methods, such as X-ray or neutron, offer non-destructive results but require access to a radiation source. With the goal of producing a more flexible solution, this LDRD developed a generalized residual stress inversion technique that can recover residual stresses released by all traction components on a cut surface, with much greater freedom in part geometry and cut location. The developed method has been successfully demonstrated on both synthetic and experimental data. The project also investigated dislocation density quantification using nonlinear ultrasound, residual stress measurement using Electronic Speckle Pattern Interferometry Hole Drilling, and validation of residual stress predictions in Additive Manufacturing process models.
Crystal plasticity-finite element method (CP-FEM) is now widely used to understand the mechanical response of polycrystalline materials. However, quantitative mesh convergence tests and verification of the necessary size of polycrystalline representative volume elements (RVE) are often overlooked in CP-FEM simulations. Mesh convergence studies in CP-FEM models are more challenging compared to conventional finite element analysis (FEA) as they are not only computationally expensive but also require explicit discretization of individual grains using many finite elements. Resolving each grains within a polycrystalline domain complicates mesh convergence study since mesh convergence is strongly affected by the initial crystal orientations of grains and local loading conditions. In this work, large-scale CP-FEM simulations of single crystals and polycrystals are conducted to study mesh sensitivity in CP-FEM models. Various factors that may affect the mesh convergence in CP-FEM simulations, such as initial textures, hardening models and boundary conditions are investigated. In addition, the total number of grains required to obtain adequate RVE is investigated. This work provides a list of guidelines for mesh convergence and RVE generation in CP-FEM modeling.
We describe an approach to predict failure in a complex, additively-manufactured stainless steel part as defined by the third Sandia Fracture Challenge. A viscoplastic internal state variable constitutive model was calibrated to fit experimental tension curves in order to capture plasticity, necking, and damage evolution leading to failure. Defects such as gas porosity and lack of fusion voids were represented by overlaying a synthetic porosity distribution onto the finite element mesh and computing the elementwise ratio between pore volume and element volume to initialize the damage internal state variables. These void volume fraction values were then used in a damage formulation accounting for growth of these existing voids, while new voids were allowed to nucleate based on a nucleation rule. Blind predictions of failure are compared to experimental results. The comparisons indicate that crack initiation and propagation were correctly predicted, and that an initial porosity field superimposed as higher initial damage may provide a path forward for capturing material strength uncertainty. The latter conclusion was supported by predicted crack face tortuosity beyond the usual mesh sensitivity and variability in predicted strain to failure; however, it bears further inquiry and a more conclusive result is pending compressive testing of challenge-built coupons to de-convolute materials behavior from the geometric influence of significant porosity.