A plot of the frequency-dependent wave propagation velocity for the case study problem with an overlocal plasticity model, with the elastic and local hardening wave speeds shown for reference (left). Stress histories using an overlocal plasticity model with a nonlocal length scale of 1m and a mesh resolution of 0.125m (right)
The following series of three articles (with common authors J. Burghardt and R. Brannon of the University of Utah) describes a state of insufficient experimental validation of conventional formulations of nonassociative plasticity (AKA nonassociated and non-normality). This work provides a confirmation that such models theoretically admit negative net work in closed strain cycles, but this simple prediction has never been validated or disproved in the laboratory!
- An early (mostly failed) attempt at experimental investigation of unvalidated plasticity assumptions (click to view),
- A simple case study confirming that nonassociativity can cause non-unique and unstable solutions to wave motion problems (click to view),
- An extensive study showing that features like rate dependence, hardening, etc. do not eliminate the instability and also showing that it is NOT related to conventional localization (click to view).
A. Sadeghirad, R. M. Brannon, and J. Burghardt
Three snapshots of the model with 248 particles in simulation of the radial expansion of a ring problem using: (a) CPDI method and (b) cpGIMP
A new algorithm is developed to improve the accuracy and efﬁciency of the material point method for problems involving extremely large tensile deformations and rotations. In the proposed procedure, particle domains are convected with the material motion more accurately than in the generalized interpolation material point method. This feature is crucial to eliminate instability in extension, which is a common shortcoming of most particle methods. Also, a novel alternative set of grid basis functions is proposed for efﬁciently calculating nodal force and consistent mass integrals on the grid. Speciﬁcally, by taking advantage of initially parallelogram-shaped particle domains, and treating the deformation gradient as constant over the particle domain, the convected particle domain is a reshaped parallelogram in the deformed conﬁguration. Accordingly, an alternative grid basis function over the particle domain is constructed by a standard 4-node ﬁnite element interpolation on the parallelogram. Effectiveness of the proposed modiﬁcations is demonstrated using several large deformation solid mechanics problems.
Below are shown comparisons of the exact and numerical solution for the vortex ring problem on a square domain.
MMS stands for “Method of Manufactured Solutions,” which is a rather sleazy sounding name for what is actually a respected and rigorous method of verifying that a finite element (or other) code is correctly solving the governing equations.
A simple introduction to MMS may be found on page 11 of The ASME guide for verification and validation in solid mechanics. The basic idea is to analytically determine forcing functions that would lead to a specific, presumably nontrivial, solution (of your choice) for the dependent variable of a differential equation. Then you would verify a numerical solver for that differential equation by running it using your analytically determined forcing function. The difference between the code’s prediction and your selected manufactured solution provides a quantitative measure of error.
Summary of the basic equations of mechanics, along with abbreviated derivations.
You may download the rest of the document here.