Contact Problems for Prestressed Elastic Bodies and Rigid and Elastic Punches

Contact problems for prestressed bodies and rigid and elastic punches are discussed. The influence of the prestresses on the contact characteristics is analyzed numerically     

Asymptotic Solutions of Boundary Value Problems for Third-Order Ordinary Differential Equations with Turning Points

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｓｉｎｇｕｌａｒｐｅｒｔｕｒｂａｔｉｏｎｓｏｆｂｏｕｎｄａｒｙｖａｌｕｅｐｒｏｂｌｅｍｓｆｏｒｏｒｄｉｎａｒｙｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓｗｉｔｈｎｏｔｕｒｎｉｎｇｐｏｉｎｔｈａｄｂｅｅｎｓｔｕｄｉｅｄｖｅｒｙｅａｒｌｙ ,ｂｕｔｔｈｅｓｔｕｄｙｏｆｂｏｕｎｄａｒｙｖａｌｕｅｐｒｏｂｌｅｍｓｗｉｔｈｔｕｒｎｉｎｇｐｏｉｎｔｐｒｏｃｅｅｄｅｄｒａｔｈｅｒｓｌｏｗｌｙｏｗｉｎｇｔｏｔｈｅｓｉｎｇｕｌａｒｐｏｉｎｔｏｆｔｈｅｄｅｇｅｎｅｒａｔｅｄｅｑｕａｔｉｏ…     

On the Asymptotic Solutions of Boundary Value Problems for a Class of Systems of Nonlinear Differential Equations (I)

A new method is applied to study the asymptotic behavior of solutions of boundary value problems for a class of systems of nonlinear differential equations . The asymptotic expansions of solutions are constructed, the remainders are estimated. The former works are improved and generalized.     

Asymptotic Research of Nonlinear Wave Processes in Saturated Porous Media

The problem of nonlinear wave dynamics of a fluid-saturated porous medium is investigated. The mathematical model proposed is based on the classical Frenkel--Biot--Nikolaevskiy theory concerning elastic wave propagation and includes mass, momentum, energy conservation laws, as well as rheological and thermodynamic relations. The model describes nonlinear, dispersive, and dissipative medium. To solve the system of differential equations, an asymptotic modified two-scales method is developed and a Cauchy problem for initial equations system is transformed to a Cauchy problem for nonlinear generalized Korteweg--de Vries--Burgers equation for modulated quick wave amplitudes and an inhomogeneous set of equations for slow background motion. Stationary solutions of the derived evolutionary equation that have been constructed numerically reflect different regimes of elastic wave attenuation: diffusive, oscillating, and soliton-like.     

Plane Elastic Boundary Value Problem Posed by Displacement and Force Orientations on a Closed Contour

A boundary value problem (BVP) of the plane elasticity posed in terms of the orientations of forces and displacements is considered. The main aim of the present paper is to investigate the solvability of BVPs of this kind. Firstly, analysis of two cases is performed: the case of a circle with special orientations of force and displacement vectors on the circumference and the case of an arbitrary contour with coaxial orientations of these vectors. The solutions obtained indicate that the problem can have a certain number of solutions or be unsolvable. Then the BVP is reduced to a boundary integral equation and its solvability is investigated for the general case of a smooth simple-connected closed contour. As a result, the number of linearly independent solutions is determined. This number only depends upon the angle between the force and displacement vectors. This revised version was published online in July 2006 with corrections to the Cover Date.     

Formulating Dynamic Multi-Rigid-Body Contact Problems with Friction as Solvable Linear Complementarity Problems

A linear complementarity formulation for dynamic multi-rigid-body contact problems with Coulomb friction is presented. The formulation, based on explicit Euler integration and polygonal approximation of the friction cone, is guaranteed to have a solution for any number of contacts and contact configuration. A model with the same property, based on the Poisson hypothesis, is formulated for impact problems with friction and nonzero restitution coefficients. An explicit Euler scheme based on these formulations is presented and is proved to have uniformly bounded velocities as the stepsize tends to zero for the Newton–Euler formulation in body co-ordinates.     

Boundary Integral Equations in Dynamic Problems for Elastic Plates

Initial-boundary value problems with Dirichlet and Neumann conditions arising in the theory of bending of plates with transverse shear deformation are reduced to time-dependent boundary integral equations by means of layer potentials. The solvability of these equations is then investigated in Sobolev-type spaces. This revised version was published online in June 2006 with corrections to the Cover Date.     

On the Form of the Contact Reaction in a Solid Circular Plate Simply Supported Along Two Antipodal Edge Arcs and Deflected by a Central Transverse Concentrated Force

A static, purely flexural mechanical analysis is presented for a Kirchhoff solid circular plate, deflected by a transverse central force, and bilaterally supported along two antipodal periphery arcs, the remaining part of the boundary being free. Two kinds of contact reactions are considered, namely the case of distributed reaction force alone, and the situation in which the distributed force is added to a distributed couple of properly selected profile. For both cases this plate problem is formulated in terms of an integral equation of the Prandtl type, coupled with two constraint conditions. The existence of solutions in an appropriate scaled weighted Sobolev space is discussed, and the behaviour of the solution at the endpoints of the support is exhibited. This revised version was published online in June 2006 with corrections to the Cover Date.     

Adjoint Equation-Based Methods for Control Problems in Incompressible, Viscous Flows

A review of adjoint equation-based methodologies for viscous,incompressible flow control and optimization problems is given and illustrated by a drag minimization example. A number of approaches to ameliorating the high storage and CPU costs associated with straightforward implementations of adjoint equation based methodologies are discussed. Other issues, including the relative merits of the differentiate-then-discretize and discretize-then-differentiate approaches to deriving discrete adjoint equations, the incorporation of side constraints into adjoint equation-based methodologies, and inaccuracies that occur due to differentiations at the boundary, are also discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

一种求解平面接触应力问题的光弹性—数值组合方法

本文提出一种求解平面接触应力最大值的光弹性—数值计算组合方法,仅从光弹性等差线上采集实验数据,然后用牛顿—拉斐逊技术与最小二乘法相结合计算出结果。文中推导了控制方程,阐述了方法,研究了实验数据取值区域、取值个数以及接触长度测试误差等因素的影响。用两类平面接触问题作为应用实例,结果符合工程要求的精度。     

On the Solution of Mixed Problems in Linear Anti-plane Piezoelectricity

In the present paper we consider interior and exterior mixed boundary value problems of anti-plane shear in the static theory of linear piezoelectricity. Using the boundary integral equation method we reduce the problems to systems of singular integral equations with discontinuous coefficients to which the classical Nöether’s theorems on existence of the solution can be applied. This allows us to establish well-posedness results and to obtain integral solutions of the corresponding mixed boundary value problems for a rather general class of piezoelectric materials. Mathematics Subject Classifications (2000) 45E05, 45F15, 74F15.     

Hertz接触问题的解的唯一性条件

讨论了接触面为圆面的Ｈｅｒｔｚ接触问题。若压力分布是轴对称的，则该接触问题的解必是唯一的。且在上述条件下，该接触问题的积分方程可化为两个推广的Ａｂｅｌ积分方程组，此方程组的解便给出此接触问题的解。     

正交异性板条的渐近分析

本文阐明由三维薄弹性体的渐近分析导出各级精度板壳理论的基本方法。将多尺度分析用于板的内部区域和边界层区域导出应力,应变和位移等物理量的不同的无量纲小厚度参数ε的渐近展开式。与工程的方法不同,推导仅基于ε→0的渐近分析,对板的变形不做任何假定。给出正交异性板的平面应变渐近展开的具体列式,仅有一些常数待定。结果表明,内域解渐近级数的首项正是熟知的Kirchhoff板理论解。本文为内域解和边界层解的渐近匹配,从而正确表述圣维南原理并用以建立各级渐近解的边界条件进而求解的研究做好了准备。     

Analytical Solution for Sliding Rounded-Edge Contact

Reliable analysis of the local stress fields in the vicinity of sliding frictional contacts of engineering components is a pre-requisite for the reliable assessment of structural integrity. In this paper we present the use of Muskhelishvili potentials to derive an analytical solution for a semi-infinite punch with a rounded edge pressed against a half-pane substrate, and a general numerical solution for a kind of Cauchy integral involved in the analytical contact solution. Using these solutions, the effect of the friction coefficient on the normal traction distribution is investigated. Numerical results show that the peak normal traction value is altered in comparison with the frictionless case solution, but this variation is mild (less than 5%), provided the friction coefficient does not exceed about 0.6. Finally, the adaptation of the analytical result to the solution of practical contact problems is addressed.     

NOVEL REGULARIZED BOUNDARY INTEGRAL EQUATIONS FOR POTENTIAL PLANE PROBLEMS

The universal practices have been centralizing on the research of regulariza-tion to the direct boundary integal equations (DBIEs). The character is elimination of singularities by using the simple solutions. However, up to now the research of regular ization to the first kind integral equations for plane potential problems has never been found in previous literatures. The presentation is mainly devoted to the research on the regularization of the singular boundary integral equations with indirect unknowns. A novel view and idea is presented herein, in which the regularized boundary integral equations with indirect unknowns without including the Cauchy principal value (CPV) and Hadamard-finite-part (HFP) integrals are established for the plane potential problems. With some numerical results, it is shown that the better accuracy and higher efficiency, especially on the boundary, can be achieved by the present system.     

Contact problems involving beams

Elastic contact problems involving Euler–Bernoulli beams or Kirchhoff plates generally involve concentrated contact forces. Linear elasticity (e.g. finite element) solutions of the same problems show that finite contact regions are actually developed, but these regions have dimensions that are typically of the order of the beam thickness. Thus if beam theory is appropriate for a given structural problem, the local elasticity fields can be explored by asymptotic methods and will have fairly general (problem independent) characteristics. Here we show that the extent of the contact region is a fixed ratio of the beam thickness which is independent of the concentrated load predicted by the beam theory, and that the distribution of contact pressure in this region has a universal form, which is well approximated by a simple algebraic expression.     

A Unified Approach to Dynamic Contact Problems in Viscoelasticity

In this paper we consider mathematical models describing dynamic viscoelastic contact problems with the Kelvin–Voigt constitutive law and subdifferential boundary conditions. We treat evolution hemivariational inequalities which are weak formulations of contact problems. We establish the existence of solutions to hemivariational inequalities with different types of non-monotone multivalued boundary relations. These results are consequences of an existence theorem for second order evolution inclusions. In a particular case we deliver sufficient conditions under which the solution to a hemivariational inequality is unique. Finally, applications to several unilateral and bilateral problems in contact mechanics are provided.*Research supported in part by the State Committee for Scientific Research of the Republic of Poland (KBN) under Grants no. 2 P03A 003 25 and 4 T07A 027 26.     

THE EXACT INTEGRAL EQUATION OF HERTZ’S CONTACT PROBLEM

This paper presents the exact integral equation of Hertz‘s contact problem,which isobtained by taking into account the horizontal displacement of points in the contactedsurfaces due to pressure.     

A Comparative Study of Least-squares,SUPG and Galerkin Methods for Convection Problems

Least-squares, SUPG and Galerkin methods are compared for a model convection problem in the context of smooth and discontinuous solutions. Our experiments focus on the impact of boundary condition implementation and grid orientation upon these methods, as well as their relative performance for both types of solutions.     

The Effect of Freezing and Discretization to the Asymptotic Stability of Relative Equilibria

In this paper we prove nonlinear stability results for the numerical approximation of relative equilibria of equivariant parabolic partial differential equations in one space dimension. Relative equilibria are solutions which are equilibria in an appropriately comoving frame and occur frequently in systems with underlying symmetry. By transforming the PDE into a corresponding PDAE via a freezing ansatz [2] the relative equilibrium can be analyzed as a stationary solution of the PDAE. The main result is the fact that nonlinear stability properties are inherited by the numerical approximation with finite differences on a finite equidistant grid with appropriate boundary conditions. This is a generalization of the results in [14] and is illustrated by numerical computations for the quintic complex Ginzburg Landau equation.