Search for: All records

Abstract This paper provides a comprehensive derivation and application of the nonlocal Nernst-Planck-Poisson (NNPP) system for accurate modeling of electrochemical corrosion with a focus on the biodegradation of magnesium-based implant materials under physiological conditions. The NNPP system extends and generalizes the peridynamic bi-material corrosion model by considering the transport of multiple ionic species due to electromigration. As in the peridynamic corrosion model, the NNPP system naturally accounts for moving boundaries due to the electrochemical dissolution of solid metallic materials in a liquid electrolyte as part of the dissolution process. In addition, we use the concept of a diffusive corrosion layer, which serves as an interface for constitutive corrosion modeling and provides an accurate representation of the kinetics with respect to the corrosion system under consideration. Through the NNPP model, we propose a corrosion modeling approach that incorporates diffusion, electromigration and reaction conditions in a single nonlocal framework. The validity of the NNPP-based corrosion model is illustrated by numerical simulations, including a one-dimensional example of pencil electrode corrosion and a three-dimensional simulation of a Mg-10Gd alloy bone implant screw decomposing in simulated body fluid. The numerical simulations correctly reproduce the corrosion patterns in agreement with macroscopic experimental corrosion data. Using numerical models of corrosion based on the NNPP system, a nonlocal approach to corrosion analysis is proposed, which reduces the gap between experimental observations and computational predictions, particularly in the development of biodegradable implant materials. 

Abstract We present an ordinary state-based peridynamic model in 2D and 3D consistent with rate-independent J2 plasticity with associated flow rule. The new contribution is the capability of the elastoplastic law to describe isotropic, kinematic and mixed hardening. The hardening formulations follow those available in the literature for classical elastoplasticity. The comparison between the results obtained with the peridynamic model and those obtained with a commercial FEM software shows that the two approaches are in good agreement. The extent of the plastic regions and von Mises stress computed with the new model for 2D and 3D examples match well those obtained with FEM-based solutions using ANSYS. 

Abstract We propose and prove several regularity criteria for the 2D and 3D Kuramoto–Sivashinsky equation, in both its scalar and vector forms. In particular, we examine integrability criteria for the regularity of solutions in terms of the scalar solution$$\phi $$ $\varphi$, the vector solution$$u\triangleq \nabla \phi $$ $u\triangleq \nabla \varphi$, as well as the divergence$$\text {div}(u)=\Delta \phi $$ $\text{div}\left(u\right)=\Delta \varphi$, and each component ofuand$$\nabla u$$ $\nabla u$. We also investigate these criteria computationally in the 2D case, and we include snapshots of solutions for several quantities of interest that arise in energy estimates. 

Abstract In this work, the essential work of fracture (EWF) method is introduced for a peridynamic (PD) material model to characterize fracture toughness of ductile materials. First, an analytical derivation for the path-independence of the PD J -integral is provided. Thereafter, the classical J -integral and PD J -integral are computed on a number of analytical crack problems, for subsequent investigation on how it performs under large scale yielding of thin sheets. To represent a highly nonlinear elastic behavior, a new adaptive bond stiffness calibration and a modified bond-damage model with gradual softening are proposed. The model is employed for two different materials: a lower-ductility bainitic-martensitic steel and a higher-ductility bainitic steel. Up to the start of the softening phase, the PD model recovers the experimentally obtained stress–strain response of both materials. Due to the high failure sensitivity on the presence of defects for the lower-ductility material, the PD model could not recover the experimentally obtained EWF values. For the higher-ductility bainitic material, the PD model was able to match very well the experimentally obtained EWF values. Moreover, the J -integral value obtained from the PD model, at the absolute maximum specimen load, matched the corresponding EWF value.