Planar multiple-contact problems subject to unilateral and bilateral kinetic constraints with static Coulomb friction

This paper is concerned with contact problems. A planar multiple-contact problem subject to unilateral and bilateral kinetic constraints with static friction is studied using the complementarity method. First, this paper discusses the one-to-one correspondence of solutions of the contact problems of concern and of the corresponding complementarity models. An enhanced complementarity model is proposed by adding missed tangential acceleration constraints into previous complementarity models. Solutions of the proposed complementarity model and solutions of the contact problem are proven to exhibit one-to-one correspondence, which may not be guaranteed in the previous complementarity models. Then, this paper applies linear complementarity theory to investigate the properties of the solutions of the proposed complementarity model. For both unilaterally constrained contact problems and bilaterally constrained contact problems, the existence of solutions and boundedness of solutions are proven. Sufficient conditions for the uniqueness of solutions and finiteness of the number of solutions are also provided. Several numerical examples are given to show the non-uniqueness of solutions or the infiniteness of the number of solutions. Such phenomena demonstrate the non-smoothness of the contact problems discussed herein.     

The adhesive contact of viscoelastic spheres

We have formulated the restricted self-consistent model for the adhesive contact of linear viscoelastic spheres. This model is a generalization of both the Ting (J. Appl. Mech. 33 (1966) 845) approach to the viscoelastic contact of adhesionless spheres and the restricted self-consistent model for adhesive axisymmetric bodies. We also show how the model can be used in practice by giving a few examples of numerical solutions.     

弹塑性接触问题数值计算的二次规划算法

针对弹塑性接触问题所推得的数值求解式子,运用二次规划法具体设计了算法,该算法采用有限的基底交换运算就可得到收敛的数值解,具有较好的收敛性及较小的计算工作量.工程计算算例结果表明文章所提出的接触问题的求解方法是有效的,针对结构中接触问题所建立的数值计算模型能真实反映实际工作状况是可靠的.文章中还详细给出了接触单元相关矩阵和向量的具体推导形式.     

Complex variable formulation for non-slipping plane strain contact of two elastic solids in the presence of interface mismatch eigenstrain

In this paper, the problems of non-slipping contact, non-slipping adhesive contact, and non-slipping adhesive contact with a stretched substrate are sequentially studied under the plane strain theory. The main results are obtained as follows:(i) The explicit solutions for a kind of singular integrals frequently encountered in contact mechanics (and fracture mechanics) are derived, which enables a comprehensive analysis of non-slipping contacts. (ii) The non-slipping contact problems are formulated in terms of the Kolosov–Muskhelishvili complex potential formulae and their exact solutions are obtained in closed or explicit forms. The relative tangential displacement within a non-slipping contact is found in a compact form. (iii) The spatial derivative of this relative displacement will be referred to in this study as the interface mismatch eigenstrain. Taking into account the interface mismatch eigenstrain, a new non-slipping adhesive contact model is proposed and its solution is obtained. It is shown that the pull-off force and the half-width of the non-slipping adhesive contact are smaller than the corresponding solutions of the JKR model (Johnson et al., 1971). The maximum difference can reach 9% for pull-off force and 17% for pull-off width, respectively. In contrast, the new model may be more accurate in modeling the non-slipping adhesion. (iv) The non-slipping adhesions with a stretch strain (S-strain) imposed to one of contact counterparts are re-examined and the analytical solutions are obtained. The accurate analysis shows that under small values of the S-strain both the natural adhesive contact half-width and the pull-off force may be augmented, but for the larger S-strain values they are always reduced. It is also found that Dundurs’ parameter β may exert a considerable effect on the solution of the pull-off problem under the S-strain.These solutions may be used to study contacts at macro-, micro-, and nano-scales.     

考虑接触应力非线性分布的接触力元模式及其验证分析

在作者提出的非连续变形计算力学模型中,采用接触力元模型描述多体接触界面上的接触特性．由于这种模型中假定接触应力沿接触界面为线性分布,从而得到的接触界面应力分布往往出现跳跃等非光滑性特征,该文对此进行了改进,采用具有高阶光滑性的非线性函数建立了能够考虑界面上接触应力非线性分布的接触力元模式,以期合理地揭示多体系统中界面的接触特性．对某一典型算例进行了数值计算,通过与大型通用非线性有限元结构分析软件ABASQUS的计算结果对比,验证了所建议计算模型的合理性与有效性．两种方法计算得到的界面接触对上的接触力基本相同;而由于采用的应力分布模式的假定不同,接触应力有所差别,由于在该文计算模型中接触对上的接触应力是按照未知量直接求得的,因此按照所建议的非线性接触力元模式所得到的接触应力更为合理．     

考虑粗糙结合面弹塑性接触的黏滑摩擦建模

提出一种同时考虑粗糙面上微凸体弹性变形和塑性接触的切向黏滑摩擦建模方法。采用Hertz弹性理论和Mindlin解描述弹性接触微凸体的切向载荷和相对变形的关系;采用AF(Abbott-Firstone)塑性理论和Fujimoto模型描述塑性接触微凸体切向载荷和相对变形的关系。再利用GW(Greenwood-Williamson)模型统计分析方法建立粗糙表面切向载荷和相对变形之间的关系。将模型与仅考虑微凸体弹性接触情况的模型进行对比,并研究了不同塑性指数对切向载荷和相对变形关系的影响。结果表明：与完全弹性接触模型相比,本文模型引入了塑性接触理论,能够更好地描述粗糙表面切向载荷和相对变形关系,并且考虑不同接触条件下弹性变形微凸体和塑性变形微凸体对切向接触载荷的贡献,在微滑移阶段,主要由弹性接触变形影响,而在进入宏观滑移阶段之后,切向行为主要由塑性变形影响。界面切向载荷由黏着和滑移接触作用共同决定,随着切向变形的增加,滑移接触力逐渐增加,而黏着接触力先增加后减少,反映了界面由微滑移逐渐向宏滑移演化的过程。随着塑性指数的增加,粗糙面上发生塑性接触的微凸体数目逐渐增加,切向黏滑行为主要受到塑性接触特征的控制。     

薄壁细胞受微吸管与探针共同作用接触模型及数值模拟

建立了以典型的薄壁球型植物细胞为原型的细胞、微吸管及探针接触模型.模型的细胞壁采用封闭球形薄膜,其本构关系为体积不可压超弹性,膜球内充满有压流体以模拟细胞质.应用轴对称几何非线性方法得出了基本微分方程组,并应用龙格-库塔法进行了求解;同时,应用流固耦合有限元法进行了数值模拟以资比较.两种方法得出了较为一致的变形和应力分...     

Symmetry Analysis and Numerical Modelling of Invasion by Malignant Tumour Tissue

We develop a model for the early stages of malignant tumour invasion dueto random motility, cellular proliferation, proteolysis and haptotaxis.At early times in the absence of tumour cell diffusion, a compressedtumour layer is evident. Transient protease production-decay dynamicsand diffusion, must be present in order for the tumour concentrationpeak to be smoothed to realistic levels.We demonstrate that invasion profiles asymptotically evolve totravelling wave solutions and that kink-like profiles, previouslythought to be due to contact inhibition and haptotaxis, can equally beexplained by cellular diffusion with a decreasing nonlinear diffusivity.As well as generalising the model and examining its robustness, afull Lie Symmetry classification is carried out.     

两种线接触非稳态部分弹流模型数值解的比较

本文利用一种快速求解弹流问题完全数值方法，求解了包含挤压项的非稳态部分弹流的Ｐａｔｉｒ和Ｃｈｅｎｇ的平均流动模型以及吴承伟的郑林庆的接触因子模型，研究了微峰接触动态变化特点以及这两种部分膜型的数值解异同，为更加直观、简单的部分弹流问题的接触因子模型的工程应用提供了理论依据。     

Spherical indentation of freestanding circular thin films in the membrane regime

We present theoretical and experimental results to describe the mechanics of indentation of a clamped circular membrane with a frictionless spherical indenter. Analytical expressions and numerical simulations are presented for the relationships between contact radius, finite indentation strains (and stresses), pre-stretch, loads and deflection. These closed-form solutions are contrasted with point-load models that neglect the contact size (i.e. classical Schwerin-type solutions), and lead to important differences in the indentation strain and load-deflection response. The accuracy of these closed form expressions is illustrated by comparisons with detailed numerical results and experiments on thin elastomer films. We show that the closed-form solutions can be used to extract mechanical properties from indentation testing of freestanding films, with important implications for developing new tests on nanoscale films and/or compliant materials such as polymers and biological substances.     

Rolling contact of a rigid sphere/sliding of a spherical indenter upon a viscoelastic half-space containing an ellipsoidal inhomogeneity

In this paper, the frictionless rolling contact problem between a rigid sphere and a viscoelastic half-space containing one elastic inhomogeneity is solved. The problem is equivalent to the frictionless sliding of a spherical tip over a viscoelastic body. The inhomogeneity may be of spherical or ellipsoidal shape, the later being of any orientation relatively to the contact surface. The model presented here is three dimensional and based on semi-analytical methods. In order to take into account the viscoelastic aspect of the problem, contact equations are discretized in the spatial and temporal dimensions. The frictionless rolling of the sphere, assumed rigid here for the sake of simplicity, is taken into account by translating the subsurface viscoelastic fields related to the contact problem. Eshelby's formalism is applied at each step of the temporal discretization to account for the effect of the inhomogeneity on the contact pressure distribution, subsurface stresses, rolling friction and the resulting torque. A Conjugate Gradient Method and the Fast Fourier Transforms are used to reduce the computation cost. The model is validated by a finite element model of a rigid sphere rolling upon a homogeneous vciscoelastic half-space, as well as through comparison with reference solutions from the literature. A parametric analysis of the effect of elastic properties and geometrical features of the inhomogeneity is performed. Transient and steady-state solutions are obtained. Numerical results about the contact pressure distribution, the deformed surface geometry, the apparent friction coefficient as well as subsurface stresses are presented, with or without heterogeneous inclusion.     

自适应无网格热弹塑性接触模型研究

提出一种自适应无网格热弹塑性接触求解模型,求解接触问题的线性规划-增量初应力法与基于应变能梯度的自适应无网格法相结合,给出了模型计算理论和算法实现.通过圆柱体与弹塑性平面热弹塑性接触算例对模型进行验证.对是否考虑材料应变硬化,是否考虑摩擦力和热输入,是否考虑材料屈服强度温度相关等情况的两种算例进行了讨论.结果表明,该模型能有效地求解考虑不同情况下的热弹塑性接触问题,在较真实地模拟接触状况的同时,具有较高的计算精度和计算效率.     

Normal Compliance Contact Models with Finite Interpenetration

We study contact problems with contact models of normal compliance type, where the compliance function tends to infinity for a given finite interpenetration. Such models are physically more realistic than standard normal compliance models, where the interpenetration is not restricted, because the interpenetration is typically justified by deformations of microscopic asperities of the surface; therefore it should not exceed a certain value that corresponds to a complete flattening of the asperities. The model can be interpreted as intermediate between the usual normal compliance models and the unilateral contact condition of Signorini type. For the static problem without friction, we prove the existence and uniqueness of solutions and establish the equivalence to an optimization problem. For the static problem with Coulomb friction, we show the existence of a solution. The analysis is based on an approximation of the problems by standard normal compliance models, the derivation of regularity results for these auxiliary problems in Sobolev spaces of fractional order by a special translation technique, and suitable limit procedures.     

A unified treatment of axisymmetric adhesive contact problems using the harmonic potential function method

A unified treatment of axisymmetric adhesive contact problems is provided using the harmonic potential function method for axisymmetric elasticity problems advanced by Green, Keer, Barber and others. The harmonic function adopted in the current analysis is the one that was introduced by Jin et al. (2008) to solve an external crack problem. It is demonstrated that the harmonic potential function method offers a simpler and more consistent way to treat non-adhesive and adhesive contact problems. By using this method and the principle of superposition, a general solution is derived for the adhesive contact problem involving an axisymmetric rigid punch of arbitrary shape and an adhesive interaction force distribution of any profile. This solution provides analytical expressions for all non-zero displacement and stress components on the contact surface, unlike existing ones. In addition, the newly derived solution is able to link existing solutions/models for axisymmetric non-adhesive and adhesive contact problems and to reveal the connections and differences among these solutions/models individually obtained using different methods at various times. Specifically, it is shown that Sneddon’s solution for the axisymmetric punch problem, Boussinesq’s solution for the flat-ended cylindrical punch problem, the Hertz solution for the spherical punch problem, the JKR model, the DMT model, the M-D model, and the M-D-n model can all be explicitly recovered by the current general solution.     

Adhesion between elastic cylinders based on the double-Hertz model

A cohesive zone model for two-dimensional adhesive contact between elastic cylinders is developed by extending the double-Hertz model of Greenwood and Johnson (1998). In this model, the adhesive force within the cohesive zone is described by the difference between two Hertzian pressure distributions of different contact widths. Closed-form analytical solutions are obtained for the interfacial traction, deformation field and the equilibrium relation among applied load, contact half-width and the size of cohesive zone. Based on these results, a complete transition between the JKR and the Hertz type contact models is captured by defining a dimensionless transition parameter μ, which governs the range of applicability of different models. The proposed model and the corresponding analytical results can serve as an alternative cohesive zone solution to the two-dimensional adhesive cylindrical contact.     

A discrete contact model for complex arbitrary-shaped convex geometries

The shape of particles has a significant influence on the behavior of suspensions, as the particle-fluid, particle-particle, and particle-wall interactions depend on it. However, the simultaneous consideration of complex particle shapes and four-way coupling remains a major challenge. This is mainly due to a lack of suitable contact models. Contact models for complex shapes have been proposed in literature, and most limit the accuracy of the particle-fluid interaction. For this reason, this paper presents a novel contact model for complex convex particle shapes for use with partially saturated methods, in which we propose to obtain necessary contact properties, such as the indentation depth, by a discretization of the contact area. The goal of the proposed model is to enable comprehensive and accurate studies of particulate flows, especially with high volume fractions, that lead to new insights and contribute to the improvement of existing industrial processes. To ensure correctness and sustainability, we validate the model extensively by studying cases with and without fluid. In the latter case, we use the homogenized lattice Boltzmann method. The provided investigations show a great agreement of the proposed discrete contact model with analytical solutions and the literature.     

Analytical Solutions for Multicomponent,Two-Phase Flow in Porous Media with Double Contact Discontinuities

This article presents the first instance of a double contact discontinuity in analytical solutions for multicomponent, two-phase flow in porous media. We use a three-component system with constant equilibrium ratios and fixed injection and initial conditions, to demonstrate this structure. This wave structure occurs for two-phase injection compositions. Such conditions were not considered previously in the development of analytical solutions for compositional flows. We demonstrate the stability of the double contact discontinuity in terms of the Liu entropy condition and also show that the resulting solution is continuously dependent on initial data. Extensions to four-component and systems with adsorption are presented, demonstrating the more widespread occurrence of this wave structure in multicomponent, two-phase flow systems. The developments in this article provide the building blocks for the development of a complete Riemann solver for general initial and injection conditions.     

Navier–Stokes solutions for steady parallel-sided pendent rivulets

We investigate exact solutions of the Navier–Stokes equations for steady rectilinear pendent rivulets running under inclined surfaces. First we show how to find exact solutions for sessile or hanging rivulets for any profile of the substrate (transversally to the direction of flow) and with no restrictions on the contact angles. The free surface is a cylindrical meniscus whose shape is determined by the static equilibrium between gravity and surface tension, by the shape of the solid surface, and by the contact angles on both contact lines. Given this, the velocity field can be obtained by integrating numerically a Poisson equation. We then perform a systematic study of rivulets hanging below an inclined plane, computing some of their global properties, and discussing their stability.     

Matched asymptotic expansions for twisted elastic knots: A self-contact problem with non-trivial contact topology

We derive solutions of the Kirchhoff equations for a knot tied on an infinitely long elastic rod subjected to combined tension and twist, and held at both endpoints at infinity. We consider the case of simple (trefoil) and double (cinquefoil) knots; other knot topologies can be investigated similarly. The rod model is based on Hookean elasticity but is geometrically nonlinear. The problem is formulated as a nonlinear self-contact problem with unknown contact regions. It is solved by means of matched asymptotic expansions in the limit of a loose knot. We obtain a family of equilibrium solutions depending on a single loading parameter (proportional to applied twisting moment divided by square root of pulling force), which are asymptotically valid in the limit of a loose knot, ε→0. Without any a priori assumption, we derive the topology of the contact set, which consists of an interval of contact flanked by two isolated points of contacts. We study the influence of the applied twist on the equilibrium.