SUMMARYNodal integration can be applied to the Galerkin weak form to yield a particle-type method where stress and material history are located exclusively at the nodes and can be employed when using meshless or finite element shape functions. This particle feature of nodal integration is desirable for large deformation settings because it avoids the remapping or advection of the state variables required in other methods. To a lesser degree, nodal integration can be desirable because it relies on fewer stress point evaluations than most other methods. In this work, aspects regarding stability, consistency, efficiency and explicit time integration are explored within the context of nodal integration. Both small and large deformation numerical examples are provided.

SUMMARYThis work introduces the weighted radial basis collocation method for boundary value problems. We first show that the employment of least-squares functional with quadrature rules constitutes an approximation of the direct collocation method. Standard radial basis collocation method, however, yields a larger solution error near boundaries. The residuals in the least-squares functional associated with domain and boundary can be better balanced if the boundary collocation equations are properly weighted. The error analysis shows unbalanced errors between domain, Neumann boundary, and Dirichlet boundary least-squares terms. A weighted least-squares functional and the corresponding weighted radial basis collocation method are then proposed for correction of unbalanced errors. It is shown that the proposed method with properly selected weights significantly enhances the numerical solution accuracy and convergence rates.

No abstract

Existing and emerging methods in computational mechanics are rarely validated against problems with an unknown outcome. For this reason, Sandia National Laboratories, in partnership with US National Science Foundation and Naval Surface Warfare Center Carderock Division, launched a computational challenge in mid-summer, 2012. Researchers and engineers were invited to predict crack initiation and propagation in a simple but novel geometry fabricated from a common off-the-shelf commercial engineering alloy. The goal of this international Sandia Fracture Challenge was to benchmark the capabilities for the prediction of deformation and damage evolution associated with ductile tearing in structural metals, including physics models, computational methods, and numerical implementations currently available in the computational fracture community. Thirteen teams participated, reporting blind predictions for the outcome of the Challenge. The simulations and experiments were performed independently and kept confidential. The methElectronic supplementary material The online version of this article (doi:10.1007/s10704-013-9904-6) contains supplementary material, which is available to authorized users.Sandia National Laboratories, Albuquerque, NM, USA e-mail: blboyce@sandia.gov ods for fracture prediction taken by the thirteen teams ranged from very simple engineering calculations to complicated multiscale simulations. The wide variation in modeling results showed a striking lack of consistency across research groups in addressing problems of ductile fracture. While some methods were more successful than others, it is clear that the problem of ductile fracture prediction continues to be challenging. Specific areas of deficiency have been identified through this effort. Also, the effort has underscored the need for additional blind prediction-based assessments.