Dynamic vibrations of a damageable viscoelastic beam in contact with two stops

A model for the material damage, due to dynamic vibrations of a Kelvin‐Voigt viscoelastic beam whose tip is constrained to move between two stops, is presented and numerically analyzed. The contact of the free tip with the stops is described by the normal compliance condition. The evolution of damage of the beam's material, which measures the reduction of its load carrying capacity, is modeled with a parabolic inclusion. The existence of the unique local solution is stated. A numerical algorithm is presented, in which spatially it is approximated by finite elements, and the time derivatives are discretized with the Euler scheme. Error estimates are derived for sufficiently regular solutions, and four numerical simulations are shown. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Thermomechanical behaviour of a damageable beam in contact with two stops

A model for the thermomechanical behaviour of a beam which allows for the general evolution of material damage is presented and investigated. One end of the beam is fixed while the other is constrained to move between two stops. The contact of the free tip with the stops is modelled by the normal compliance condition. The thermal interaction between the stops and the free tip is described by a heat exchange condition where the heat transfer coefficient is a general function of the gaps between the tip and the stops. The effects on the mechanical properties of the material due to crack expansion are described by a damage field, which measures the decrease in the load-bearing capacity of the material. The damage evolves as a constrained diffusion process in which the microcracks that develop may grow or disappear. The mathematical model consists of a coupled system of energy--elasticity equations together with a nonlinear parabolic inclusion for the damage field. The existence of a local solution is established using truncation, penalization, and a priori estimates.     

A frictional contact problem with normal damped response for elastic-visco-plastic materials

We study an evolution problem which describes the dynamic contact of an elastic-visco-plastic body with a foundation. We model the contact with normal damped response and a local friction law. A damage of the material caused by elastic deformation is taken into account, its evolution is described by an inclusion of parabolic type. We derive variational formulation for the model which is in the form of a system involving the displacement field, the stress field and the damage field. We prove the existence and uniqueness result of the weak solution. The proof is based on arguments of evolution equations with monotone operators, a classical existence and uniqueness result on parabolic inequalities and fixed point.     

Existence and Regularity for Dynamic Viscoelastic Adhesive Contact with Damage

A model for the dynamic process of frictionless adhesive contact between a viscoelastic body and a reactive foundation, which takes into account the damage of the material resulting from tension or compression, is presented. Contact is described by the normal compliance condition. Material damage is modelled by the damage field, which measures the pointwise fractional decrease in the load-carrying capacity of the material, and its evolution is described by a differential inclusion. The model allows for different damage rates caused by tension or compression. The adhesion is modelled by the bonding field, which measures the fraction of active bonds on the contact surface. The existence of the unique weak solution is established using the theory of set-valued pseudomonotone operators introduced by Kuttler and Shillor (1999). Additional regularity of the solution is obtained when the problem data is more regular and satisfies appropriate compatibility conditions.     

Variational and numerical analysis of a quasistatic viscoelastic problem with normal compliance,friction and damage

We consider a model for quasistatic frictional contact between a viscoelastic body and a foundation. The material constitutive relation is assumed to be nonlinear. The mechanical damage of the material, caused by excessive stress or strain, is described by the damage function, the evolution of which is determined by a parabolic inclusion. The contact is modeled with the normal compliance condition and the associated version of Coulomb's law of dry friction. We derive a variational formulation for the problem and prove the existence of its unique weak solution. We then study a fully discrete scheme for the numerical solutions of the problem and obtain error estimates on the approximate solutions.     

Numerical analysis of a frictional contact problem for viscoelastic materials with long-term memory

We consider a mathematical model which describes the bilateral contact between a deformable body and an obstacle. The process is quasistatic, the material is assumed to be viscoelastic with long memory and the friction is modeled with Tresca’s law. The problem has a unique weak solution. Here we study spatially semi-discrete and fully discrete schemes using finite differences and finite elements. We show the convergence of the schemes under the basic solution regularity and we derive order error estimates. Finally, we present an algorithm for the numerical realization and simulations for a two-dimensional test problem.     

One-dimensional dynamic thermoviscoelastic contact with damage

The problem of thermoviscoelastic dynamic contact between a rod and a rigid obstacle, when the material damage is taken into account, is modeled and analyzed. The contact is modeled by the normal compliance condition and the stress-strain constitutive equation is of Kelvin-Voigt type. The damage, which describes the reduction of the load carrying capacity of the rod, evolves because of the opening of microcracks as a result of tension or compression. When the damage reaches a critical value at a point on the rod the material cannot carry any load and the system breaks down. Mathematically, this is expressed by the quenching of the solution. The existence of a local weak solution is established using penalization and a priori estimates.     

Transient solutions of Markov processes and generalized continued fractions

In 1974 J. A. Murphy and M. R. O'Donohoe numerically approximatedthe minimal solution of the Kolmogorov forward equation forthe generalized birth and death process by use of continuedfractions. This paper generalizes this approach by suggestingan algorithm for q-matrices of lower band structure (n, 1).This is achieved by analogy with generalized continued fractions.Applications involving q-matrices of this type include, forexample, many types of queueing systems with batch processingor birth–death–catastrophe population processesin biology.     

Damageable contact between an elastic body and a rigid foundation

In this work, the contact problem between an elastic body and a rigid obstacle is studied, including the development of material damage which results from internal compression or tension. The variational problem is formulated as a first-kind variational inequality for the displacements coupled with a parabolic partial differential equation for the damage field. The existence of a unique local weak solution is stated. Then, a fully discrete scheme is introduced using the finite element method to approximate the spatial variable and an Euler scheme to discretize the time derivatives. Error estimates are derived on the approximate solutions, from which the linear convergence of the algorithm is deduced under suitable regularity conditions. Finally, three two-dimensional numerical simulations are performed to demonstrate the accuracy and the behaviour of the scheme.     

A dynamic viscoelastic contact problem with normal compliance,finite penetration and nonmonotone slip rate dependent friction

We consider a mathematical model which describes the dynamic evolution of a viscoelastic body in frictional contact with an obstacle. The contact is modelled with normal compliance and unilateral constraint, associated to a rate slip-dependent version of Coulomb’s law of dry friction. In order to approximate the contact conditions, we consider a regularized problem wherein the contact is modelled by a standard normal compliance condition without finite penetrations. For each problem, we derive a variational formulation and an existence result of the weak solution of the regularized problem is obtained. Next, we prove the convergence of the weak solution of the regularized problem to the weak solution of the initial nonregularized problem. Then, we introduce a fully discrete approximation of the variational problem based on a finite element method and on a second order time integration scheme. The solution of the resulting nonsmooth and nonconvex frictional contact problems is presented, based on approximation by a sequence of nonsmooth convex programming problems. Finally, some numerical simulations are provided in order to illustrate both the behaviour of the solution related to the frictional contact conditions and the convergence result.     

Shape and topology optimization of contact problems by level set methods

This paper deals with the numerical solution of a topology and shape optimization problems of an elastic body in unilateral contact with a rigid foundation. The contact problem with the prescribed friction is considered. The structural optimization problem consists in finding such shape of the boundary of the domain occupied by the body that the normal contact stress along the contact boundary of the body is minimized. In the paper shape as well as topological derivatives formulae of the cost functional are provided using material derivative and asymptotic expansion methods, respectively. These derivatives are employed to formulate necessary optimality condition for simultaneous shape and topology optimization. Level set based numerical algorithm for the solution of the shape optimization problem is proposed. Level set method is used to describe the position of the boundary of the body and its evolution on a fixed mesh. This evolution is governed by Hamilton – Jacobi equation. The speed vector field driving the propagation of the boundary of the body is given by the shape derivative of a cost functional with respect to the free boundary. Numerical examples are provided. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Multi-symplectic Runge-Kutta-type methods for Hamiltonian wave equations

^** Email: kzhang{at}math.cuhk.edu.hk The non-linear wave equation is taken as a model problem forthe investigation. Different multi-symplectic reformulationsof the equation are discussed. Multi-symplectic Runge–Kuttamethods and multi-symplectic partitioned Runge–Kutta methodsare explored based on these different reformulations. Some popularand efficient multi-symplectic schemes are collected and constructed.Stability analyses are performed for these schemes.     

Singularities of Hyperbolic Gauss Maps

In this paper we adopt the hyperboloid in Minkowski space asthe model of hyperbolic space. We define the hyperbolic Gaussmap and the hyperbolic Gauss indicatrix of a hypersurface inhyperbolic space. The hyperbolic Gauss map has been introducedby Ch. Epstein [J. Reine Angew. Math. 372 (1986) 96–135]in the Poincaré ball model, which is very useful forthe study of constant mean curvature surfaces. However, it isvery hard to perform the calculation because it has an intrinsicform. Here, we give an extrinsic definition and we study thesingularities. In the study of the singularities of the hyperbolicGauss map (indicatrix), we find that the hyperbolic Gauss indicatrixis much easier to calculate. We introduce the notion of hyperbolicGauss–Kronecker curvature whose zero sets correspond tothe singular set of the hyperbolic Gauss map (indicatrix). Wealso develop a local differential geometry of hypersurfacesconcerning their contact with hyperhorospheres. 2000 MathematicalSubject Classification: 53A25, 53A05, 58C27.     

Wear resistance of bodies protected by a thin composite coating

The mechanical contact interaction of bodies with a thin composite coating is investigated with account of wear. The thermal effects are not considered. The coating is modeled by a thin plate. Between the body and the coating is an interlayer, which is modeled by a Winkler body with one modulus of subgrade reaction. Under the action of a rigid stamp on the coating, the process of abrasive wear proceeds. The contact interaction of the coating with the base is described by using the model of an intermediate layer. To determine the stress-strain state of the coating, equations of the generalized theory of plates including the shear strains and the compression of normal are utilized. For the contact wear problem formulated, the basic integral equation with a Fredholm-type kernel is derived, and its solution algorithm is proposed. Numerical results are presented. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 42, No. 3, pp. 319–330, May–June, 2006.     

A regularized damage-plasticity model: micromorphic formulation and numerical aspects

A gradient-extended damage-plasticity model is discussed which is based on a micromorphic approach according to Forest [1]. Damage and plasticity are treated as independent but strongly coupled dissipative phenomena by considering separate yield and damage loading functions to describe the onset of plastic flow and / or damage evolution. A numerical benchmark test conducted in the study reveals that the model is able to essentially cure the well-known mesh-dependence issue which is known from finite element simulations involving conventional (i. e. ‘local’) damage material models. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A combined BIE-FE method for the Stokes equations

A numerical algorithm for the biharmonic equation in domainswith piecewise smooth boundaries is presented. It is intendedfor problems describing the Stokes flow in the situations whereone has corners or cusps formed by parts of the domain boundaryand, due to the nature of the boundary conditions on these partsof the boundary, these regions have a global effect on the shapeof the whole domain and hence have to be resolved with sufficientaccuracy. The algorithm combines the boundary integral equationmethod for the main part of the flow domain and the finite-elementmethod which is used to resolve the corner/cusp regions. Twoparts of the solution are matched along a numerical ‘internalinterface’ or, as a variant, two interfaces, and theyare determined simultaneously by inverting a combined matrixin the course of iterations. The algorithm is illustrated byconsidering the flow configuration of ‘curtain coating’,a flow where a sheet of liquid impinges onto a moving solidsubstrate, which is particularly sensitive to what happens inthe corner region formed, physically, by the free surface andthe solid boundary. The ‘moving contact line problem’is addressed in the framework of an earlier developed interfaceformation model which treats the dynamic contact angle as partof the solution, as opposed to it being a prescribed functionof the contact line speed, as in the so-called ‘slip models’.     

Quasistatic evolution of damage in an elastic body with random inputs

Abstract

A model for the quasistatic evolution of the motion of an elastic body which is subject to material damage is presented and analyzed. The model takes the form of an elliptic system for the displacements coupled with a parabolic inclusion for the damage field. In both the system and the inclusion the coefficient and input functions that are present are assumed to be stochastic processes dependent on a random variable. The existence of a weak solution to the model is established using a sequence of approximate problems and passage to a limit. Moreover, the weak solution is shown to be product measurable.     

Global existence for a nonlocal model for adhesive contact

Abstract

In this paper, we address the analytical investigation into a model for adhesive contact introduced in a paper by Freddi and Fremond, which includes nonlocal sources of damage on the contact surface, such as the elongation. The resulting PDE system features various nonlinearities rendering the unilateral contact conditions, the physical constraints on the internal variables, as well as the contributions related to the nonlocal forces. For the associated initial-boundary value problem, we obtain a global-in-time existence result by proving the existence of a local solution via a suitable approximation procedure and then by extending the local solution to a global one by a nonstandard prolongation argument.     

Quadratic Pareto optimal control of parabolic equation with state-control constraints and an infinite number of variables

A distributed Pareto optimal control problem for the parabolicoperator with an infinite number of variables is considered.The performance index has an integral form. Constraints on controlsand on states are imposed. To obtain optimality conditions forthe Neumann problem, the generalization of the Dubovitskii–MilyutinTheorem given by WALCZAK, S. (1984a) Folia Mathematica, 1, 187–196and (1984b) J. Optimiz. Theory Appl., 42, 561–582, wasapplied.     

On the formation of acceleration and reaction-diffusion wavefronts in autocatalytic-type reaction-diffusion systems

^** Email: d.j.needham{at}reading.ac.uk We consider generalization of the theory for the evolution ofreaction–diffusion and accelerating wavefronts in KPP-typesystems as developed in Needham (2004, Proc. R. Soc. Lond. A,460, 1921–1934) (DN). These generalizations allow forthe removal of a number of technical restrictions imposed inthe paper of DN.