We propose an unconditionally stable semi-implicit time discretization of the phase field crystal evolution. It is based on splitting the underlying energy into convex and concave parts and then performing H-1 gradient descent steps implicitly for the former and explicitly for the latter. The splitting is effected in such a way that the resulting equations are linear in each time step and allow an extremely simple implementation and efficient solution. We provide the associated stability and error...

In order to study the effects of earthquakes on tombstones and monumental columns in recent years the dynamical analysis of rigid blocks subjected to ground accelerations has interested many researchers. Mainly, the rocking motion has been investigated and many numerical difficulties have been pointed out in such analysis [1-2-3-4]. Some computational advantages can be obtained by modelling the bonding between two blocks or between block and foundation by means of an elastic layer of Winkler's springs...

A new very high-order finite volume method to solve problems with harmonic and biharmonic operators for onedimensional geometries is proposed. The main ingredient is polynomial reconstruction based on local interpolations of mean values providing accurate approximations of the solution up to the sixth-order accuracy. First developed with the harmonic operator, an extension for the biharmonic operator is obtained, which allows designing a very high-order finite volume scheme where the solution is...

We present here a theoretical study of eigenmodes in axisymmetric elastic layers. The mathematical modelling allows us to bring this problem to a spectral study of a sequence of unbounded self-adjoint operators $A_n$, $n\in \mathbb{N}$, in a suitable Hilbert space. We show that the essential spectrum of $A_n$ is an interval of type $[\gamma ,+\infty [$ and that, under certain conditions on the coefficients of the medium, the discrete spectrum is non empty.

We present here a theoretical study of eigenmodes in axisymmetric elastic layers. The mathematical modelling allows us to bring this problem to a spectral study of a sequence of unbounded self-adjoint operators An, $n\in \mathbb{N}$, in a suitable Hilbert space. We show that the essential spectrum of An is an interval of type $[\gamma ,+\infty [$ and that, under certain conditions on the coefficients of the medium, the discrete spectrum is non empty.

The purpose of this paper is to provide a priori error estimates on the approximation of contact conditions in the framework of the eXtended Finite-Element Method (XFEM) for two dimensional elastic bodies. This method allows to perform finite-element computations on cracked domains by using meshes of the non-cracked domain. We consider a stabilized Lagrange multiplier method whose particularity is that no discrete inf-sup condition is needed in the convergence analysis. The contact condition is...

The purpose of this paper is to provide a priori error estimates on the approximation of contact conditions in the framework of the eXtended Finite-Element Method (XFEM) for two dimensional elastic bodies. This method allows to perform finite-element computations on cracked domains by using meshes of the non-cracked domain. We consider a stabilized Lagrange multiplier method whose particularity is that no discrete inf-sup condition is needed in the convergence analysis. The contact condition is...

The purpose of this paper is to provide a priori error estimates on the approximation of contact conditions in the framework of the eXtended Finite-Element Method (XFEM) for two dimensional elastic bodies. This method allows to perform finite-element computations on cracked domains by using meshes of the non-cracked domain. We consider a stabilized Lagrange multiplier method whose particularity is that no discrete inf-sup condition is needed in the convergence analysis. The contact condition is...

In this work we consider a stabilized Lagrange multiplier method in order to approximate the Coulomb frictional contact model in linear elastostatics. The particularity of the method is that no discrete inf-sup condition is needed. We study the existence and the uniqueness of solution of the discrete problem.

The aim of this paper is to study a quasistatic unilateral contact problem between an elastic body and a foundation. The constitutive law is nonlinear and the contact is modelled with a normal compliance condition associated to a unilateral constraint and Coulomb's friction law. The adhesion between contact surfaces is taken into account and is modelled with a surface variable, the bonding field, whose evolution is described by a first-order differential equation. We establish a variational formulation...

The operator $L_0:D_{L_0}\subset H\to H$, $L_{0}u=\frac{1}{r}\frac{d}{dr}\left\{r\frac{d}{dr}\left[\frac{1}{r}\frac{d}{dr}\left(r\frac{du}{dr}\right)\right]\right\}$, $D_{L_0}=\{u\in C^4\left([0,R]\right),u^{\text{'}}\left(0\right)=u^{\text{'}\text{'}\text{'}\text{'}}\left(0\right)=0,u\left(R\right)=u^{\text{'}}\left(R\right)=0\}$, $H=L_{2,r}(0,R)$ is shown to be essentially self-adjoint, positive definite with a compact resolvent. The conditions on $L_0$ (in fact, on a general symmetric operator) are given so as to justify the application of the Fourier method for solving the problems of the types $L_{0}u=g$ and $u_{tt}+L_{0}u=g$, respectively.

In this paper we present a unified approach to obtain integral representation formulas for describing the propagation of bending waves in infinite plates. The general anisotropic case is included and both new and well-known formulas are obtained in special cases (e.g. the classical Boussinesq formula). The formulas we have derived have been compared with experimental data and the coincidence is very good in all cases.

In this report, a control method for the stabilization of periodic orbits for a class of one- and two-dimensional discrete-time systems that are topologically conjugate to symbolic dynamical systems is proposed and applied to a population model in an ecosystem and the Smale horseshoe map. A periodic orbit is assigned as a target by giving a sequence in which symbols have periodicity. As a consequence, it is shown that any periodic orbits can be globally stabilized by using arbitrarily small control...

We study preconditioning techniques for discontinuous Galerkin discretizations of isotropic linear elasticity problems in primal (displacement) formulation. We propose subspace correction methods based on a splitting of the vector valued piecewise linear discontinuous finite element space, that are optimal with respect to the mesh size and the Lamé parameters. The pure displacement, the mixed and the traction free problems are discussed in detail. We present a convergence analysis of the proposed...