Crack nucleation in a peridynamic solid

A condition for the emergence of a discontinuity in an elastic peridynamic body is proposed, resulting in a material stability condition for crack nucleation. The condition is derived by determining whether a small discontinuity in displacement, superposed on a possibly large deformation, grows over time. Stability is shown to be determined by the sign of the eigenvalues of a tensor field that depends only on the linearized material properties. This condition for nucleation of a discontinuity in displacement can be interpreted in terms of the dynamic stability of plane waves with very short wavelength. A numerical example illustrates that cracks in a peridynamic body form spontaneously as the body is loaded

Positive approximations of the inverse of fractional powers of SPD M-matrices

This study is motivated by the recent development in the fractional calculus and its applications. During last few years, several different techniques are proposed to localize the nonlocal fractional diffusion operator. They are based on transformation of the original problem to a local elliptic or pseudoparabolic problem, or to an integral representation of the solution, thus increasing the dimension of the computational domain. More recently, an alternative approach aimed at reducing the computational complexity was developed. The linear algebraic system $\cal A^\alpha \bf u=\bf f$, $0< \alpha <1$ is considered, where $\cal A$ is a properly normalized (scalded) symmetric and positive definite matrix obtained from finite element or finite difference approximation of second order elliptic problems in $\Omega\subset\mathbb{R}^d$, $d=1,2,3$. The method is based on best uniform rational approximations (BURA) of the function $t^{\beta-\alpha}$ for $0 < t \le 1$ and natural $\beta$. The maximum principles are among the major qualitative properties of linear elliptic operators/PDEs. In many studies and applications, it is important that such properties are preserved by the selected numerical solution method. In this paper we present and analyze the properties of positive approximations of $\cal A^{-\alpha}$ obtained by the BURA technique. Sufficient conditions for positiveness are proven, complemented by sharp error estimates. The theoretical results are supported by representative numerical tests

Towards an Efficient Finite Element Method for the Integral Fractional Laplacian on Polygonal Domains

We explore the connection between fractional order partial differential equations in two or more spatial dimensions with boundary integral operators to develop techniques that enable one to efficiently tackle the integral fractional Laplacian. In particular, we develop techniques for the treatment of the dense stiffness matrix including the computation of the entries, the efficient assembly and storage of a sparse approximation and the efficient solution of the resulting equations. The main idea consists of generalising proven techniques for the treatment of boundary integral equations to general fractional orders. Importantly, the approximation does not make any strong assumptions on the shape of the underlying domain and does not rely on any special structure of the matrix that could be exploited by fast transforms. We demonstrate the flexibility and performance of this approach in a couple of two-dimensional numerical examples

A peridynamic based machine learning model for one-dimensional and two-dimensional structures

With the rapid growth of available data and computing resources, using data-driven models is a potential approach in many scientific disciplines and engineering. However, for complex physical phenomena that have limited data, the data-driven models are lacking robustness and fail to provide good predictions. Theory-guided data science is the recent technology that can take advantage of both physics-driven and data-driven models. This study presents a novel peridynamics based machine learning model for one and two-dimensional structures. The linear relationships between the displacement of a material point and displacements of its family members and applied forces are obtained for the machine learning model by using linear regression. The numerical procedure for coupling the peridynamic model and the machine learning model is also provided. The numerical procedure for coupling the peridynamic model and the machine learning model is also provided. The accuracy of the coupled model is verified by considering various examples of a one-dimensional bar and two-dimensional plate. To further demonstrate the capabilities of the coupled model, damage prediction for a plate with a pre-existing crack, a two-dimensional representation of a three-point bending test, and a plate subjected to dynamic load are simulated

Derivation of dual horizon state-based peridynamics formulation based on euler-lagrange equation

The numerical solution of peridynamics equations is usually done by using uniform spatial discretisation. Although implementation of uniform discretisation is straightforward, it can increase computational time significantly for certain problems. Instead, non-uniform discretisation can be utilised and different discretisation sizes can be used at different parts of the solution domain. Moreover, the peridynamic length scale parameter, horizon, can also vary throughout the solution domain. Such a scenario requires extra attention since conservation laws must be satisfied. To deal with these issues, dual-horizon peridynamics was introduced so that both non-uniform discretisation and variable horizon sizes can be utilised. In this study, dual-horizon peridynamics formulation is derived by using Euler–Lagrange equation for state-based peridynamics. Moreover, application of boundary conditions and determination of surface correction factors are also explained. Finally, the current formulation is verified by considering two benchmark problems including plate under tension and vibration of a plate

Positive correlation between Merkel cell polyomavirus viral load and capsid-specific antibody titer

Merkel cell polyomavirus (MCPyV or MCV) is the first polyomavirus to be clearly implicated as a causal agent underlying a human cancer, Merkel cell carcinoma (MCC). Infection with MCPyV is common in the general population, and a majority of adults shed MCPyV from the surface of their skin. In this study, we quantitated MCPyV DNA in skin swab specimens from healthy volunteers sampled at different anatomical sites over time periods ranging from 3 months to 4 years. The volunteers were also tested using a serological assay that detects antibodies specific for the MCPyV virion. There was a positive correlation between MCPyV virion-specific antibody titers and viral load at all anatomical sites tested (dorsal portion of the hands, forehead, and buttocks) (Spearman’s r 0.644, P < 0.0001). The study results are consistent with previous findings suggesting that the skin is primary site of chronic MCPyV infection in healthy adults and suggest that the magnitude of an individual’s seroresponsiveness against the MCPyV virion generally reflects the overall MCPyV DNA load across wide areas of the skin. In light of previous reports indicating a correlation between MCC and strong MCPyV-specific seroresponsiveness, this model suggests that poorly controlled chronic MCPyV infection might be a risk factor in the development of MCC