Interacting indentors on a poroelastic half-space

This paper examines the interaction between two rigid circular indentors on a poroelastic half-space. The resulting mixed boundary value problem, when formulated in the Laplace transform domain, yields an infinite set of Fredholm integral equations. These integral equations are then solved for some special cases. Numerical results for the case of a single indentor show a good agreement with those obtained by using Heinrich and Desoyer's assumption. For the case in which the radius of one indentor reduces to zero (interaction between a rigid indentor and an externally placed load), the resulting equations are solved by a semi-inverse method to give analytical solutions for the resultant force and moment required to maintain the indentor with no normal displacement. When the indentor is subjected to an axial load but allowed to undergo an additional settlement and tilt, numerical results are presented to demonstrate the manner in which Poisson's ratio and the drainage boundary conditions influence the consolidation of the half-space. Numerical results are also given to illustrate the interaction between two identical indentors when ratio of the radius to the spatial distance between them is small.     

On magnetoelastic problems of a plane with an arbitrarily shaped hole under stress and displacement boundary conditions

An analytical method for the static plane problem of magnetoelasticityis developed for an infinite plane containing a hole of arbitraryshape under stress and displacement boundary conditions in aprimary uniform magnetic field. The magnetic field influencesthe elastic field by introducing a body force called the Lorentzponderomotive force in the equilibrium equations. The body forcecan be further described in a form relating with the electromagneticstress tensor. The complex variable method in conjunction withthe rational mapping function technique is used in the analysisfor both magnetic field and mechanical field. Governing equationsand boundary conditions are expressed in terms of complex functions.Complex magnetic potential and stress functions are obtainedusing Cauchy integrals for the paramagnetic and soft ferromagneticmaterials, respectively. The distributions of magnetic fieldand the stress components are shown for certain directions ofprimary magnetic fields in an infinite plane with a square hole,as an example. It is found that the stress distributions forthe two types of materials are identical despite the differenceof magnetic fields. The extreme cases of a free and a fixedhole reduced to a crack and a rigid fibre, respectively, arealso investigated. The stress intensity factors at the tipsof crack and rigid fibre are computed, and their variation forcertain directions of primary magnetic field is shown.     

On the linearized theory op flow around bodies. Method of force sources

The method of force sources is proposed for solving linear problems related to the interaction between rigid bodies, and fluids, or gases. Method is based on the introduction of perturbation force sources into equation of motion of fluid media. Boundary conditions at the rigid body surface make it possible to reduce the problem of hydrodynamic reactions to an integral equation defining the function of force sources. Method is illustrated by the solution of three simple problems in the field of acoustics, and of viscous, and compressible media flow around bodies.

In the linearized theory of flow around rigid bodies, as well as in acoustics, an important part of the sound wave generation analysis concerns the determination of hydrodynamic reactions of the medium on moving, pulsating, or oscillating bodies. Such reactions make themselves felt as constant, or variable mechanical forces, such as drag and lift, or in the case of sound wave emitters, as the wave resistance. Various methods had been proposed for the computation of such forces, as for example, in the monographs [1 to 6].

Here, a different approach to the problem of determination of surface forces exerted by liquids and gases on the rigid body is proposed. By resorting to the formalism of the generalized functions it is possible to introduce into the equations of motion of fluid media a perturbation source in the form volume density of forces exercised by the body on the gas. The distribution of surface tension entering into the expression of this force is selected in such a manner as to satisfy boundary conditions at the body surface. It becomes possible with the use of this device to reduce the problem of determination of forces acting on the body surface to the solution of certain Integral equations. The proposed method is in all respects completely analogous to the well-known method of sources and sinks [1 to 1]. Both methods reduce the problem of interaction between body and gas to the solution of Integral equations. The method of sources and sinks, however, leads to an integral equation which describes the distribution of fictitious sources and sinks in the volume of the body having the density of the medium, while the method of force sources yields an integral equation which directly defines the distribution of mechanical forces over the surface of the body (*).

We may note that the method of force sources had to a certain extent been already used in papers [6 and 7] for the determination of sound radiation by means of point-force sources.     

Numerical solution of an equilibrium problem for an elastic body with a thin delaminated rigid inclusion

Under consideration is a 2D-problem of elasticity theory for a body with a thin rigid inclusion. It is assumed that there is a delamination crack between the rigid inclusion and the elastic matrix. At the crack faces, the boundary conditions are set in the form of inequalities providing mutual nonpenetration of the crack faces. Some numerical method is proposed for solving the problem, based on domain decomposition and the Uzawa algorithm for solving variational inequalities.We give an example of numerical calculation by the finite element method.     

Solvability of a nonstationary problem of a rigid body motion in the flow of a viscous incompressible fluid in a pipe of arbitrary cross-section

The existence of a generalized weak solution is proved for the nonstationary problem of motion of a rigid body in the flow of a viscous incompressible fluid filling a cylindrical pipe of arbitrary cross-section. The fluid flow conforms to the Navier–Stokes equations and tends to the Poiseuille flow at infinity. The body moves in accordance with the laws of classical mechanics under the influence of the surrounding fluid and the gravity force directed along the cylinder. Collisions of the body with the boundary of the flow domain are not admitted and, by this reason, the problem is considered until the body approaches the boundary.     

The self-propulsion of a body with moving internal masses in a viscous fluid

An investigation of the characteristics of motion of a rigid body with variable internal mass distribution in a viscous fluid is carried out on the basis of a joint numerical solution of the Navier — Stokes equations and equations of motion for a rigid body. A nonstationary three-dimensional solution to the problem is found. The motion of a sphere and a drop-shaped body in a viscous fluid in a gravitational field, which is caused by the motion of internal material points, is explored. The possibility of self-propulsion of a body in an arbitrary given direction is shown.     

On the fall of a heavy rigid body in an ideal fluid

We consider a problem about the motion of a heavy rigid body in an unbounded volume of an ideal irrotational incompressible fluid. This problem generalizes a classical Kirchhoff problem describing the inertial motion of a rigid body in a fluid. We study different special statements of the problem: the plane motion and the motion of an axially symmetric body. In the general case of motion of a rigid body, we study the stability of partial solutions and point out limiting behaviors of the motion when the time increases infinitely. Using numerical computations on the plane of initial conditions, we construct domains corresponding to different types of the asymptotic behavior. We establish the fractal nature of the boundary separating these domains.     

Axisymmetrical vibrations of shells of revolution under abrupt loading

A frequency method is proposed for solving the problem of the vibrations of shells of revolution taking into account the energy dissipation under arbitrary force loading and on collision with a rigid obstacle. The Laplace transform is taken of the equation of the vibrations of a shell of revolution with non-zero initial conditions. For the inhomogeneous differential equation obtained, a variational method is used to solve the boundary-value problem, which consists of finding the Laplace-transformed boundary transverse and longitudinal forces and bending moments as functions of the boundary displacements. The equations of equilibrium of nodes, i.e. the corresponding equations of the finite-element method, are then compared, using results obtained earlier [1–4]. Amplitude-phase-frequency characteristics (APFCs) for the shell cross-sections selected are plotted. An inverse Laplace transformation is carried out using the clear relationship between the extreme points of the APFCs and the coefficients of the corresponding terms of the series in an expansion vibration modes [3]. In view of the fact that the proposed approach is approximate, numerical testing is used.     

The motion of a body with a cavity partly filled with a viscous liquid

Many papers are concerned with the dynamics of a rigid body with a cavity filled with liquid (see the bibliography in [1]). The present paper deals with the motion of a rigid body having a cavity partly filled with a viscous incompressible liquid, and having a free surface. The shape of the cavity is arbitrary. The problem is considered in a linear formulation. The oscillations of the body with respect to its center of inertia and the motion of the liquid in the cavity are assumed small. The viscosity of the liquid is considered low. The solution of the problem of the oscillations of a body with a cavity partly filled with an ideal liquid is used as an initial approximation [1 to 6]. The viscosity is taken into consideration by the boundary layer method used before in similar problems [1 and 7 to 10). General equations are derived for the dynamics of a body filled with a liquid, for an arbitrary form of cavity. The coefficients of those integro-differential equations depend only on the solution of the problem of the oscillations of a body with a cavity of the given form filled with an ideal liquid. Since the corresponding problem has been solved for cavities of many forms [1 to 6, 11 and 12] in the case of an ideal liquid, the determination of the characteristic coefficients is reduced to the evaluation of quadratures. Several particular cases of motion are considered.     

带弹性附件充液矩形贮箱俯仰运动动态响应

首先建立了俯仰运动矩形贮箱刚-液-弹耦合系统在外力矩作用下的耦合动力学模型,给出满足边界条件的速度势函数和液面波高的级数表达式,采用伽辽金法离散,将动力学模型转化为常微分方程组,得到刚-液-弹耦合系统的固有频率,给出简单的近似表达式,分析了转动中心距静液面不同位置时刚-液-弹耦合系统各阶固有频率的变化规律,系统转动中心距静液面较近时,耦合后液体反对称模态和刚体的固有频率对比耦合前减小,较远时则增大,最后进行数值验证,比较分析了液体和弹性体对刚体姿态的影响．     

Modelling two-dimensional flow past arbitrary cylindrical bodies using boundary element formulations

This paper presents an efficient method of solving Queen's linearized equations for steady plane flow of an incompressible, viscous Newtonian fluid past a cylindrical body of arbitrary cross-section. The numerical solution technique is the well known direct boundary element method. Use of a fundamental solution of Oseen's equations, the ‘Oseenlet’, allows the problem to be reduced to boundary integrals and numerical solution then only requires boundary discretization. The formulation and solution method are validated by computing the net forces acting on a single circular cylinder, two equal but separated circular cylinders and a single elliptic cylinder, and comparing these with other published results. A boundary element representation of the full Navier-Stokes equations is also used to evaluate the drag acting on a single circular cylinder by matching with the numerical Oseen solution in the far field.     

Potential flow past a porous body with a core of different permeability

It is well known that a uniform flow past a non-permeable rigid body does not exert a total force upon the surface of the body, however this is not the case when the body is permeable. Power et. al. (1984, 1986) first solved the problem of uniform potential flow past a two-dimensional permeable circular cylinder, with constant permeability, and found that the exterior flow exerts a drag force upon the surface of the cylinder independent of its size and secondly the problem when the uniform potential flow past a porous sphere, with constant permeability, in this case the exterior flow exerts a drag force on the sphere which is linearly dependent on the radius of the sphere. Here we will present the solution of two problems, a uniform potential flow past a porous circular cylinder and past a porous sphere, for each case the porous body is composed of two materials with different permeabilities. In both cases the total force exerted by the exterior flow upon the body is dependent on the thickness of the porous materials, and in the limit when the two permeabilities are equal, the previous results, circular cylinder and sphere, with constant permeability, are recovered. Atlhough, the mathematics involved in the solution of the present problem is simple, due to the nice boundary geometry of the bodies, the final expression for the total force found in each case is quite interesting on the way it depends on the permeability relation, in particular, in the limiting cases of a porous body with solid or hollow core.     

Approximate solutions and error estimates for a Stokes boundary value problem

The aim of this paper is the numerical treatment of a boundary value problem for the system of Stokes’ equations. For this we extend the method of approximate approximations to boundary value problems. This method was introduced by Maz’ya (DFG-Kolloquium des DFG-Forschungsschwerpunktes Randelementmethoden, 1991) and has been used until now for the approximation of smooth functions defined on the whole space and for the approximation of volume potentials. In the present paper we develop an approximation procedure for the solution of the interior Dirichlet problem for the system of Stokes’ equations in two dimensions. The procedure is based on potential theoretical considerations in connection with a boundary integral equations method and consists of three approximation steps as follows. In a first step the unknown source density in the potential representation of the solution is replaced by approximate approximations. In a second step the decay behavior of the generating functions is used to gain a suitable approximation for the potential kernel, and in a third step Nyström’s method leads to a linear algebraic system for the approximate source density. For every step a convergence analysis is established and corresponding error estimates are given.     

The uniform heating and bending of a two-layer plate

The problem of the uniform heating of a two-layer plate is solved. The transversely isotropic layer considered (a soft plate) is in ideal contact with a rigid isotropic thin elastically deformed layer. The ends of the plate are load-free. A boundary layer of the soft plate (a thin contact layer) is introduced, which enables the boundary conditions on the ends of the plate to be formulated in such a way that the problem has a bounded smooth solution [1]. The two-layer plate, generally speaking, is bounded along the axis perpendicular to the axes directed along the length and thickness of the plate. The resultant force and the resultant moment, applied to the end transverse sections, are equal to zero. The exact solution of the temperature problem is sought using the equations of the theory of elasticity. The plane problem of the bending of a two-layer plate acted upon by a uniformly distributed pressure applied to the side surface of an anisotropic layer is solved by a similar method. The ends of the rigid isotropic layer are clamped.     

非对称载荷作用的外部圆形裂纹问题

使用边界积分方程方法,研究了三维无限弹性体中受非对称载荷作用的外部圆形裂纹问题。通过使用Fourier级数和超几何函数,将问题的二维边界奇异积分方程简化为Abel型方程,获得了一般非对称载荷作用的外部圆形裂纹问题的应力强度因子精确解,比用Hankel变换法得到的结果更为一般。结果表明:边界积分方程法在解析分析方面还有很大的潜力。     

Method of integral equations in the scalar problem of diffraction on a system consisting of a “soft” and a “hard” screen and an inhomogeneous body

We consider the scalar problem on the diffraction of a plane wave on a system of two screens with boundary conditions of the first and the second kind and a solid inhomogeneous body in the semiclassical setting. The original boundary value problem for the Helmholtz equation is reduced to a system of singular integral equations over the body domain and the screen surfaces. We prove the equivalence of the integral and differential statements of the problem, the solvability of the system of integral equations in Sobolev spaces, and the smoothness of its solutions. To solve the integral equations approximately, we use the Bubnov-Galerkin method; we introduce basis functions on the body and the screens and prove the consistency and convergence of the numerical method.     

Existence of strong solutions for the problem of a rigid-fluid system

This Note is devoted to the study of a fluid–rigid body interaction problem. The motion of the fluid is modelled by the Navier–Stokes equations, written in an unknown bounded domain depending on the displacement of the rigid body. Our main result yields the existence and uniqueness of strong solutions, which are global provided that the rigid body does not touch the boundary. To cite this article: T. Takahashi, C. R. Acad. Sci. Paris, Ser. I 336 (2003).     

Existence of solutions for the equations

We introduce a concept of weak solution for a boundary value problem modelling the interactive motion of a coupled system consisting in a rigid body immersed in a viscous fluid. The fluid, and the solid are contained in a fixed open bounded set of R3. The motion of the fluid is governed by the incompresible Navier-Stokes equations and the standard conservation's laws of linear, and angular momentum rules the dynamics of the rigid body. The time variation of the fluid's domain (due to the motion of the rigid body) is not known apriori, so we deal with a free boundary value problem. Our main theorem asserts the existence of at least one weak solution for this problem. The result is global in time provided that the rigid body does not touch the boundary     

Numerical methods for some inverse problems of heart electrophysiology

We consider numerical methods for solving inverse problems that arise in heart electrophysiology. The first inverse problem is the Cauchy problem for the Laplace equation. Its solution algorithm is based on the Tikhonov regularization method and the method of boundary integral equations. The second inverse problem is the problem of finding the discontinuity surface of the coefficient of conductivity of a medium on the basis of the potential and its normal derivative given on the exterior surface. For its numerical solution, we suggest a method based on the method of boundary integral equations and the assumption on a special representation of the unknown surface.     

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)