For elastic-perfectly plastic solids (or structures) subjected to quasi-static cyclic loads, variational methods are presented for the direct eyâluation of the post-transient residual stresses, that is, the residual stresses in the structure at the end of the transient response phase, consequence of the plastic strains therein produced and crucial to predict the subsequent steady structural behaviour. The problem of the evaluation of the number of cycles spanned by the transient response is also...

We study the use of a GPU for the numerical approximation of the curvature dependent flows of graphs - the mean-curvature flow and the Willmore flow. Both problems are often applied in image processing where fast solvers are required. We approximate these problems using the complementary finite volume method combined with the method of lines. We obtain a system of ordinary differential equations which we solve by the Runge-Kutta-Merson solver. It is a robust solver with an automatic choice of the...

Let $K:=SO\left(2\right)A_1\cup SO\left(2\right)A_2\cdots SO\left(2\right)A_N$ where $A_1,A_2,\cdots ,A_N$ are matrices of non-zero determinant. We establish a sharp relation between the following two minimisation problems in two dimensions. Firstly the $N$-well problem with surface energy. Let $p\in \left[1,2\right]$, $\Omega \subset {ℝ}^2$ be a convex polytopal region. Define$$I_{ϵ}^p\left(u\right)={\int }_{\Omega }{d}^p\left(Du\left(z\right),K\right)+ϵ{\left|D^{2}u\left(z\right)\right|}^2\mathrm{d}{L}^{2}z$$and let $A_F$ denote the subspace of functions in $W^{2,2}\left(\Omega \right)$ that satisfy the affine boundary condition $Du=F$ on $\partial \Omega$ (in the sense of trace), where $F\notin K$. We consider the scaling (with respect to $ϵ$) of$$m_{ϵ}^p:=\underset{u\in A_F}{inf}{I}_{ϵ}^p\left(u\right).$$Secondly the finite element approximation to the $N$-well problem without surface...

Let $K:=SO\left(2\right)A_1\cup SO\left(2\right)A_2\cdots SO\left(2\right)A_N$ where $A_1,A_2,\cdots ,A_N$ are matrices of non-zero determinant. We establish a sharp relation between the following two minimisation problems in two dimensions. Firstly the N-well problem with surface energy. Let $p\in \left[1,2\right]$, $\Omega \subset {ℝ}^2$ be a convex polytopal region. Define $$I_{ϵ}^p\left(u\right)={\int }_{\Omega }{d}^p\left(Du\left(z\right),K\right)+ϵ{\left|D^{2}u\left(z\right)\right|}^2\mathrm{d}{L}^{2}z$$ and let AF denote the subspace of functions in $W^{2,2}\left(\Omega \right)$ that satisfy the affine boundary condition Du=F on $\partial \Omega$ (in the sense of trace), where $F\notin K$. We consider the scaling (with respect to ϵ) of $$m_{ϵ}^p:=\underset{u\in A_F}{inf}{I}_{ϵ}^p\left(u\right).$$ Secondly the finite element approximation to the N-well problem without...

For a general class of atomistic-to-continuum coupling methods, coupling multi-body interatomic potentials with a P1-finite element discretisation of Cauchy–Born nonlinear elasticity, this paper adresses the question whether patch test consistency (or, absence of ghost forces) implies a first-order error estimate. In two dimensions it is shown that this is indeed true under the following additional technical assumptions: (i) an energy consistency condition, (ii) locality of the interface correction,...

For a general class of atomistic-to-continuum coupling methods, coupling multi-body interatomic potentials with a P1-finite element discretisation of Cauchy–Born nonlinear elasticity, this paper adresses the question whether patch test consistency (or, absence of ghost forces) implies a first-order error estimate. In two dimensions it is shown that this is indeed true under the following additional technical assumptions: (i) an energy consistency condition, (ii) locality of the interface correction,...

In this paper we propose an original approach for the simulation of the time-dependent response of a floating elastic plate using the so-called Singularity Expansion Method. This method consists in computing an asymptotic behaviour for large time obtained by means of the Laplace transform by using the analytic continuation of the resolvent of the problem. This leads to represent the solution as the sum of a discrete superposition of exponentially damped oscillating motions associated to the poles...

We estimate the constant in the strengthened Cauchy-Bunyakowski-Schwarz inequality for hierarchical bilinear finite element spaces and elliptic partial differential equations with coefficients corresponding to anisotropy (orthotropy). It is shown that there is a nontrivial universal estimate, which does not depend on anisotropy. Moreover, this estimate is sharp and the same as for hierarchical linear finite element spaces.