Shakedown in elastic contact problems with Coulomb friction

Elastic systems with frictional interfaces subjected to periodic loading are sometimes predicted to ‘shake down’ in the sense that frictional slip ceases after the first few loading cycles. The similarities in behaviour between such systems and monolithic bodies with elastic–plastic constitutive behaviour have prompted various authors to speculate that Melan’s theorem might apply to them – i.e., that the existence of a state of residual stress sufficient to prevent further slip is a sufficient condition for the system to shake down.In this paper, we prove this result for ‘complete’ contact problems in the discrete formulation (i) for systems with no coupling between relative tangential displacements at the interface and the corresponding normal contact tractions and (ii) for certain two-dimensional problems in which the friction coefficient at each node is less than a certain critical value. We also present counter-examples for all systems that do not fall into these categories, thus giving a definitive statement of the conditions under which Melan’s theorem can be used to predict whether such a system will shake down.     

Shakedown limits for an oscillating,elastic rolling contact with Coulomb friction

We examine experimentally and theoretically the effect of frictional shakedown of a three-dimensional elastic rolling contact. Small oscillations of the local normal forces lead to incremental sliding processes within the area of contact. Consequently, this causes a macroscopic slip motion of the two contacting bodies. If the oscillation amplitude is sufficiently small, the frictional slip ceases after the first few loading periods and a safe shakedown occurs. Otherwise the slip motion is continued and the contact fails.     

Shakedown analysis of shape memory alloy structures

Phase transformational shakedown of a structure refers to a status that plastic strains cease developing after a finite number of loading cycles, and subsequently the structure undergoes only elastic deformation and alternating phase transformations with limited magnitudes. Due to the intrinsic complexity in the constitutive relations of shape memory alloys (SMA), there is as yet a lack of effective methods for modeling the mechanical responses of SMA structures, especially when they develop both phase transformation and plastic deformation. This paper is devoted to present an algorithm for analyzing shakedown of SMA structures subjected to cyclic or varying loads within specified domains. Based on the phase transformation and plastic yield criteria of von Mises-type and their associated flow rules, a simplified three-dimensional phenomenological constitutive model is first formulated accounting for different regimes of elastic–plastic deformation and phase transformation. Different responses possible for SMA bodies exposed to varying loads are discussed. The classical Melan shakedown theorem is extended to determine a lower bound of loads for transformational shakedown of SMA bodies without necessity of a step-by-step analysis along the loading history. Finally, a simple example is given to illustrate the application of the present theory as well as some basic features of shakedown of SMA structures. It is interesting to find that phase transformation may either increase or decrease the load-bearing capacity of a structure, depending upon its constitutive relations, geometries and the loading mode.     

Shakedown reduced kinematic formulation,separated collapse modes,and numerical implementation

Our shakedown reduced kinematic formulation is developed to solve some typical plane stress problems, using finite element method. Whenever the comparisons are available, our results agree with the available ones in the literature. The advantage of our approach is its simplicity, computational effectiveness, and the separation of collapse modes for possible different treatments. Second-order cone programming developed for kinematic plastic limit analysis is effectively implemented to study the incremental plasticity collapse mode. The approach is ready to be used to solve general shakedown problems, including those for elastic–plastic kinematic hardening materials and under dynamic loading.     

A method for the evaluation of design limits for structural materials in a cyclic state of creep

The purpose of the present paper is to demonstrate how the minimum theorems proposed in an accompanying paper (Ponter and Boulbibane, 2002) can be utilised in the prediction of the deformation and life assessment of structures subjected to cyclic mechanical and thermal loadings. The developed method, which is based upon bounding theorems and an associate programming method, the Linear Matching method, takes into account the changes in residual stress field occurring within a cycle. Although the solution provided a bound on the inelastic work, it also appears that generally the displacements predicted by this solution are smaller than those that would be predicted by the rapid cycle solution. By way of illustration a simple non-linear viscous model is adopted and a number of solutions are presented involving a Bree plate problem subjected to cyclic histories of load and temperature. An elastic follow-up factor is identified as a key design parameter for high temperature dwell periods.     

A linear complementarity model for multibody systems with frictional unilateral and bilateral constraints

The Lagrange-I equations and measure differential equations for multibody systems with unilateral and bilateral constraints are constructed.For bilateral constraints,frictional forces and their impulses contain the products of the filled-in relay function induced by Coulomb friction and the absolute values of normal constraint reactions.With the time-stepping impulse-velocity scheme,the measure differential equations are discretized.The equations of horizontal linear complementarity problems(HLCPs),which are used to compute the impulses,are constructed by decomposing the absolute function and the filled-in relay function.These HLCP equations degenerate into equations of LCPs for frictional unilateral constraints,or HLCPs for frictional bilateral constraints.Finally,a numerical simulation for multibody systems with both unilateral and bilateral constraints is presented.     

Modeling and analysis of rigid multibody systems with driving constraints and frictional translation joints

An approach is proposed for modeling and anal- yses of rigid multibody systems with frictional translation joints and driving constraints. The geometric constraints of translational joints with small clearance are treated as bilat- eral constraints by neglecting the impact between sliders and guides. Firstly, the normal forces acting on sliders, the driv- ing constraint forces （or moments） and the constraint forces of smooth revolute joints are all described by complementary conditions. The frictional contacts are characterized by a set- valued force law of Coulomb＇s dry friction. Combined with the theory of the horizontal linear complementarity problem （HLCP）, an event-driven scheme is used to detect the transi- tions of the contact situation between sliders and guides, and the stick-slip transitions of sliders, respectively. And then, all constraint forces in the system can be computed easily. Secondly, the dynamic equations of multibody systems are written at the acceleration-force level by the Lagrange multiplier technique, and the Baumgarte stabilization method is used to reduce the constraint drift. Finally, a numerical example is given to show some non-smooth dynamical behaviors of the studied system. The obtained results validate the feasibility of algorithm and the effect of constraint stabilization.     

Fixed Points of Multi-valued Maps and Static Coulomb Friction Problems

The existence of solutions to the Signorini problem with Coulomb friction is a long standing open question. We prove the existence of generalized solutions that satisfy the pointwise Coulomb friction conditions on the entire interface and the normal nonpenetration condition on the complement of a subset with arbitrarily small but possibly positive measure. Furthermore, the penetration itself can also be made arbitrarily small. Although “measure zero” instead of “arbitrarily small measure” would be needed to fully resolve the issue, these generalized solutions seem to be the closest answer available to date. Their existence is proved by a suitable application of Ky Fan's fixed point theorem for multi-valued maps. The same method can be used with a number of variants involving contact of two or more elastic bodies and possible debonding phenomena. This revised version was published online in August 2006 with corrections to the Cover Date.     

Shakedown analysis for a class of strengthening materials within the framework of gradient plasticity

The classical shakedown theory is extended to a class of perfectly plastic materials with strengthening effects (Hall–Petch effects). To this aim, a strain gradient plasticity model previously advanced by Polizzotto (2010) is used, whereby a featuring strengthening law provides the strengthening stress, i.e. the increase of the yield strength produced by plastic deformation, as a degree-zero homogeneous second-order differential form in the accumulated plastic strain with associated higher order boundary conditions. The extended static (Melan) and kinematic (Koiter) shakedown theorems are proved together with the related lower bound and upper bound theorems. The shakedown limit load problem is addressed and discussed in the present context, and its solution uniqueness shown out. A simple micro-scale structural system is considered as an illustrative example. The shakedown limit load is shown to increase with decreasing the structural size, which is a manifestation of the classical Hall–Petch effects in a context of cyclic loading.     

Shakedown theory for elastic plastic kinematic hardening bodies

Shakedown static and kinematic theorems for elastic–plastic (generally nonlinear) kinematic hardening solids are derived in classical (path-independence) spirit with new constructions. The generally plastic-deformation-history-dependent hardening curve is assumed to be limited by the initial yield stress and ultimate yield strength, and to obey a positive hysteresis postulate for closed plastic cycles, but else can be arbitrary and unspecified. The theorems reveal that the shakedown of structures is not affected by the particular form of the hardening curve, but just by the initial and ultimate yield stresses. While the ultimate yield strength is clearly defined macroscopically and attached to the incremental collapse mode with unbounded plastic deformations, the initial yield stress, which is responsible for the bounded cyclic plasticity collapse mode, should not be taken as the convenient one at a fixed amount of plastic deformation (0.2%), but is suggested to be taken as low as the fatigue limit to preserve the classical load-history-independence spirit of the shakedown theorems. Otherwise, for our pragmatic application purpose, it may be given empirical values between the low fatigue limit and high ultimate yield stresses according to particular loading processes considered, which may range anywhere between the high-cycle and low-cycle ones. The theorems appear as simple as those of Melan and Koiter for perfect plasticity but applied to the much larger class of more realistic kinematic hardening materials.     

薄板结构的安定性分析

本文关于薄板的分析采用样条函数插值及结构中的残余应力场用温度参数模拟，解决了薄板结构在交变载荷作用下的安定问题。由于温度参数的引用，使得研究的问题的自由度大幅度减少，文中还对薄板的Ｍｉｓｅｓ型屈服函数进行了线性化处理，它是安定定理的求解线性化的基础，算例表明本文所提出的方法及程序的可行和有效。     

Slip-shakedown analysis of a system of circular beams in frictional contact

In automotive components, the cumulative microslip phenomenon is often observed for engine assemblies. This phenomenon results in an accumulation of the relative slips in a preferred tangential direction on the contact interface of two solids under cyclic loadings. A significant relative displacement may occur and leads to the assembly failure. In particular, a global rotation of the bearing shell may result from this mechanism of cumulated slips in conrod big end systems. To discuss this rotation problem, a model of two circular beams in frictional contact and submitted to a periodical rotating load is considered here. The aim is to give some simplified estimates of the critical rotation load based on a slip-shakedown analysis. The discussion holds for Tresca friction and can be extended to Coulomb friction under the assumption of small coupling. The static and kinematic slip-shakedown approaches are discussed. The obtained analytical results are shown to be in agreement with the finite element computations.     

组合旋转壳体的安定上限分析

对旋转组合壳体的各个子区的位移场进行双三次Ｂ样条插值，同时对子区交界处的广义力向量以单三次Ｂ样条展开，推导出以广义向量表示的组合壳体的平衡方程，然后从Ｂｅｔｉ互换定理和虚功原理出发推导出离散的壳体结构中可能的塑性应变率循环，从而给出了组合壳体安定上限分析的有限元列式，把组合壳体的安定上限分析归结为求解一个数学规划问题，文末给出了算例     

薄板结构的安定优化设计

本文通过对薄板的Ｂ样条函数插值及结构中残余应力场的温度参数模拟，研究了薄板安定分析的样条有限元列式，在分析结构安定分析的数学规划列式的基础上，解决了薄板结构在给定可变载荷作用下的安定最小厚度设计问题。由于温度参数的引用，使得研究的问题的自由度大幅度减少。算例表明本文所提出的方法及所编程序的可行和有效。     

The shakedown load boundary of an elastic-perfectly plastic structure

In the hypothesis of small displacements and combined time-variable/steady loads, the geometrical-mechanical properties of the shakedown load boundary are investigated. It is shown that, in the load space, the shakedown load boundary plays the role of yield surface, and that a certain plastic strain accumulation vector—characterizing some impending inadaptation collapse mechanism—obeys the normality rule, whereas a specific form of the maximum plastic work theorem constitutes an effective tool for the evaluation of the shakedown limit load corresponding to a specified inadaptation collapse mode. The equations governing the state of the structure at the shakedown limit are provided and the related collapse mechanism is shown to specify the shape of the steady-state response of the structure to a periodic load enveloping the load domain with an intensity slightly above the shakedown limit.

---

Sommario Nell'ipotesi di piccoli spostamenti e di plasticità perfetta, si studiano le proprietà geometrico-meccaniche della frontiera del dominio dei carichi di adattamento. Si mostra come, nello spazio dei carichi, la suddetta frontiera giuoca il ruolo di superficie di plasticizzazione ed un particolare vettore di deformazione plastica accumulata—che caratterizza il meccanismo di non adattamento—segue la legge della normalità, mentre una forma specifica del teorema del massimo lavoro plastico costituisce un valido strumento per la valutazione del carico limite corrispondente ad un dato meccanismo di non adattamento. Si forniscono altresì le equazioni che governano lo stato della struttura al limite di adattamento, mostrando la capacità del relativo meccanismo di collasso a rappresentare la risposta a lungo termine della struttura a dei carichi periodici inviluppanti il dominio dei carichi con una intensità un poco al di sopra del limite di adattamento.

     

Dynamics of spatial structure-varying rigid multibody systems

Summary Couplings in machines and mechanisms exhibiting backlash and friction phenomena can be modeled as multibody systems with unilateral constraints and Coulomb friction. The structure of the differential-algebraic equations describing the system depends on the state of the constraints. The contact forces occurring at active constraints are taken into account in the equations of motion as Lagrange multipliers. Additionally, the kinematic conditions of all active constraints are formulated on the acceleration level. Contact and friction laws are sufficient conditions for state transitions of active constraints, and are represented by nonsmooth characteristics. Several formulations, like the linear complementarity problem, and two different nonlinear systems of equations are presented together with their solution method. The theory is applied to a mechanical system containing three-dimensional and coupled unilateral constraints with friction. Received 14 May 1998; accepted for publication 5 January 1999     

On oscillations of a frictional pendulum

Conditions are established under which a standard limit cycle occurs in the system under consideration, or the trajectory closes under the influence of a stagnation domain. It is pointed out that when the solution falls into the stagnation domain it makes no sense to use the asymptotic method because of a large error __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 2, pp. 104–112, February 2006.