Application of the method of spatial characteristics to the solution of axially symmetric problems relating to the propagation of elastic waves

We present a difference scheme, based on the method of spatial characteristics, for solving axially symmetric dynamical problems of the theory of elasticity. Consideration is given to the possibility of solving a Cauchy problem and a problem for a solid or a hollow cylinder which takes boundary conditions into account. It is suggested that linear problems may be solved by this method. An example is given in which the parameters characterizing the stress-deformation state of a semiinfinite cylinder are calculated, the points of the end of the cylinder being given an initial axial velocity. The calculation of these parameters was carried out on the BÉSM-6 computer.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 101–109, July–August, 1971.     

Application of the angular superposition method to the contact problem on the compression of an elastic cylinder

The angular superposition method is used to construct an approximate solution of the contact problem on the compression of an elastic cylinder by two rigid plates. The solution thus obtained has a closed-form analytic expression and can be used in the entire domain of the cylinder cross-section. We analyze the absolute error, which takes the largest value near the points of contact between the plates and the cylinder, where the boundary conditions are discontinuous. According to the von Mises criterion, when moving into the depth of the cylinder from the contact site along the symmetry axis, the second invariant J ₂ of the stress deviator tensor first decreases and then, after attaining a minimum, increases and attains the largest value at a small depth, which agrees with Johnson’s photoelastic experiments and Dinnik’s computations. We present the graphs of the displacement and normal stress distributions over the contact site, the dependence of the compressing force on the displacements of rigid plates, and the dependence of the invariant J ₂ on the coordinate along the symmetry axis. If 640 computation points are chosen on the cylinder boundary and the Hertz law for the normal pressure on the contact site is used, then the error in the approximate solution near the endpoint of the contact site is approximately 55%, and if the proposed two-parameter normal law is used, then the error is of the order of 4%. On the free lateral surface of the cylinder boundary, we find the critical pointM*, which separates the cylinder contraction and extension parts.The contact problems are the most difficult problems, and their solution is complicated by the discontinuous boundary conditions [1–5]. In [6], the contact problem is solved by the Fourier method, which can be used only for bodies of classical shapes. In such cases, the problem can be reduced to solving coupled integral equations [7]. The interaction between the bandage and a cylindrical body is considered in [2, 6, 7]. In [8], the possibility of using the finite element method is investigated in the case of contact problems for a differential wheel with roughness of the contacting surfaces taken into account. In [9, 10], the method of homogeneous solutions is used to consider contact problems for a finite-dimensional elastic cylinder loaded on its end surfaces. Note that only error estimates are given in the literature cited above; the absolute error over the entire domain of the elastic body is not studied, although this is one of the important characteristics of the obtained approximate solution. A sufficiently complete survey of the literature in the field of contact interactions of elastic bodies is given in [3–5].In what follows, we propose to solve contact problems by the angular superposition method [11]. This method can be used for bodies of nonclassical shapes, which can be multiply connected, and the friction on the contact site can be taken into account. In the present paper, as a first example of applied character, we show how this method can be used in the simplest case. The multiple connectedness and the curvilinearity of the shape of the body, as well as taking into account the friction on the boundary, do not create new essential difficulties in this method.     

A hybrid immersed‐boundary and multi‐block lattice Boltzmann method for simulating fluid and moving‐boundaries interactions

A numerical method is developed for modelling the interactions between incompressible viscous fluid and moving boundaries. The principle of this method is introducing the immersed‐boundary concept in the framework of the lattice Boltzmann method, and improving the accuracy and efficiency of the simulation by refining the mesh near moving boundaries. Besides elastic boundary with a constitutive law, the method can also efficiently simulate solid moving‐boundary interacting with fluid by employing the direct forcing technique. The method is validated by the simulations of flow past a circular cylinder, two cylinders moving with respect to each other and flow around a hovering wing. The versatility of the method is demonstrated by the numerical studies including elastic filament flapping in the wake of a cylinder and fish‐like bodies swimming in quiescent fluid. Copyright © 2006 John Wiley & Sons, Ltd.     

An explicit formulation for the evolution of nonlinear surface waves interacting with a submerged body

An explicit formulation to study nonlinear waves interacting with a submerged body in an ideal fluid of infinite depth is presented. The formulation allows one to decompose the nonlinear wave–body interaction problem into body and free‐surface problems. After the decomposition, the body problem satisfies a modified body boundary condition in an unbounded fluid domain, while the free‐surface problem satisfies modified nonlinear free‐surface boundary conditions. It is then shown that the nonlinear free‐surface problem can be further reduced to a closed system of two nonlinear evolution equations expanded in infinite series for the free‐surface elevation and the velocity potential at the free surface. For numerical experiments, the body problem is solved using a distribution of singularities along the body surface and the system of evolution equations, truncated at third order in wave steepness, is then solved using a pseudo‐spectral method based on the fast Fourier transform. A circular cylinder translating steadily near the free surface is considered and it is found that our numerical solutions show excellent agreement with the fully nonlinear solution using a boundary integral method. We further validate our solutions for a submerged circular cylinder oscillating vertically or fixed under incoming nonlinear waves with other analytical and numerical results. Copyright © 2007 John Wiley & Sons, Ltd.     

Contact problems for a circular plate with sliding fixation at the end face

Two mixed elasticity problems of punch indentation into a circular plate placed without clearance in a rigid cylindrical holder with smooth walls are considered. In the first problem, the plate lies without friction on a rigid base, and in the second problem, the plate is rigidly fixed to the base. The problems are solved by a method that was developed for bodies of finite dimensions and is based on the properties of closed systems of orthogonal functions. Each of the problems is reduced to two integral equations, namely, a Volterra integral equation of the first kind for the contact pressure function and a Fredholm integral equation of the first kind for the derivatives of the displacement of the plate upper surface outside the punch. The displacement function is sought as the sum of a trigonometric series and a power function with a root singularity. After truncation, the obtained illposed system of linear algebraic equation has a stable solution. A method for solving Volterra integral equations is given. The contact pressure distribution function and the dimensionless indentation force are determined. Examples of calculation of the plate interaction with the plane punch are given. Contact problems were earlier studied for a rectangle and a circular plate with a stress-free end both without taking account of their fixation [1, 2] and with regard for their fixation [3, 4]. The solution method described here was used to study the interaction of elastic hollow cylinder of finite length with a rigid bandage and a rigid insert [5, 6]. Other papers dealing with contact problems for bodies of finite dimensions, in particular, for a circular plate, should also be mentioned. In these papers, the problems under study were solved by the method of homogeneous solutions [7, 8] and by the method of coupled series-equations [9].     

Design of Plane Thermoelastic Composite Constructions with Uniformly Stressed Reinforcement

The problem of uniformly stressed reinforcement of plane composite constructions under thermoforce loading is formulated. An asymptotic analysis of the corresponding boundary–value problem is performed. Based on this analysis, it is shown that the problem may have two solutions due to the significant nonlinearity of static boundary conditions. An iterative method for solving the problem is proposed. Particular analytical and numerical solutions are analyzed, and the level of influence of the thermal action on uniformly stressed reinforced constructions is studied.     

Existence of Weak Solutions for the Unsteady Interaction of a Viscous Fluid with an Elastic Plate

The purpose of this work is to study the existence of solutions for an unsteady fluid-structure interaction problem. We consider a three-dimensional viscous incompressible fluid governed by the Navier–Stokes equations, interacting with a flexible elastic plate located on one part of the fluid boundary. The fluid domain evolves according to the structure’s displacement, itself resulting from the fluid force. We prove the existence of at least one weak solution as long as the structure does not touch the fixed part of the fluid boundary. The same result holds also for a two-dimensional fluid interacting with a one-dimensional membrane.     

Three-dimensional thin viscous shock layer in the absence of planes of symmetry in the flow

The basic laws of viscous homogeneous gas flow at high supersonic speeds past smooth blunt bodies with a permeable surface are investigated within the framework of the thin viscous shock layer model. An efficient numerical method of solving these equations, which makes it possible to consider cases of flow past bodies at angles of attack and slip, when there are no planes of symmetry in the flow, is proposed. Some results of calculating the flow past a triaxial ellipsoid with an axial ratio of 103n73 at angles of attack =0–45° and slip angles =0–45° over a broad interval of Reynolds numbers are presented as an example. The effect of the principal determining parameters of the problem on the flow structure in the shock layer and the surface friction and heat transfer coefficients is analyzed. An expression for calculating the heat fluxes to the impermeable surface of smooth blunt bodies in a supersonic homogeneous viscous gas flow over a broad interval of Reynolds numbers is proposed on the basis of the solutions obtained and the results of other authors.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 150–158, March–April, 1989.     

Three-dimensional laminar boundary layer on a permeable surface in the neighborhood of a symmetry plane

Analytical and numerical methods are used to investigate a three-dimensional laminar boundary layer near symmetry planes of blunt bodies in supersonic gas flows. In the first approximation of an integral method of successive approximation an analytic solution to the problem is obtained that is valid for an impermeable surface, for small values of the blowing parameter, and arbitrary values of the suction parameter. An asymptotic solution is obtained for large values of the blowing or suction parameters in the case when the velocity vector of the blown gas makes an acute angle with the velocity vector of the external flow on the surface of the body. Some results are given of the numerical solution of the problem for bodies of different shapes and a wide range of angles of attack and blowing and suction parameters. The analytic and numerical solutions are compared and the region of applicability of the analytic expressions is estimated. On the basis of the solutions obtained in the present work and that of other authors, a formula is proposed for calculating the heat fluxes to a perfectly catalytic surface at a symmetry plane of blunt bodies in a supersonic flow of dissociated and ionized air at different angles of attack. Flow near symmetry planes on an impermeable surface or for weak blowing was considered earlier in the framework of the theory of a laminar boundary layer in [1–4]. An asymptotic solution to the equations of a three-dimensional boundary layer in the case of strong normal blowing or suction is given in [5, 6].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 37–48, September–October, 1980.     

Initial value problem for the dynamic system of anisotropic elasticity

The main object of this paper is the Cauchy problem for the dynamic system of anisotropic elasticity. Existence and uniqueness theorems of weak and smooth solutions of this problem are established by the reduction of the original elasticity system into a symmetric hyperbolic system of the first order. The numerical method of the Cauchy problem solving for anisotropic elastic system with polynomial data is obtained and its correctness is established. The simulations of the numerical solutions are presented.     

基于混合变量条形传递函数方法的圆柱壳屈曲分析

提出了一种分析交各向异性圆柱壳和阶梯圆柱壳稳定性问题的混合变量条形传递函数方法。首先基于Fluegge薄壳理论，通过定义广义位移变量和对应的广义力变量，建立了圆柱壳混合变量能量泛函；然后通过引入条形单元，定义混合状态变量和采用传递函数方法对超级壳单元求解，得到具有多种边界条件圆柱壳屈曲问题的半解析解；最后通过位移连续和力平衡条件，可以得到阶梯圆柱壳屈曲问题的解。理论解推导过程表明此方法在引入边界条件和进行阶梯圆柱壳求解时非常方便。算例分析的结果验证了本方法的正确性。     

Aspects of self-similar diffraction gas flows and their numerical simulation

A large number of papers has been devoted to the investigation of the interaction of a plane shock wave with bodies of various geometric shapes, and they have been generalized and classified for a stationary body in [1, 2]. Separate results of experimental and theoretical investigations of the interaction of a shock wave with a wedge, cone, sphere, and cylinder moving with supersonic velocities are contained in [3–9]. Analysis of the available results shows that the features of the unsteady gas flows formed in this case largely depend on the nature of the boundary-value problem that arises for the system of differential gas dynamic equations. The question of the wave structure of the unsteady gas flow and the accuracy of the obtained solution is central to the numerical investigation of the present class of problems. The most characteristic types of unsteady self-similar gas flows that arise on the interaction of a plane shock wave with bodies of a wedge or convex corner type are calculated on the basis of an explicit numerical continuous calculation method of the second order of accuracy. The accuracy of the numerical solutions is discussed on the basis of a comparison with the experimental data. The case of the interaction of a shock wave with the rarefaction wave that arises in a supersonic flow past a convex corner is considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 146–152, July–August, 1986.     

Poromechanical response of a finite elastic cylinder under axisymmetric loading

Drained or undrained cylindrical specimens under axisymmetric loading are commonly used in laboratory testing of soils and rocks. Poroelastic cylindrical elements are also encountered in applications related to bioengineering and advanced materials. This paper presents an analytical solution for an axisymmetrically-loaded solid poroelastic cylinder of finite length with permeable (drained) or impermeable (undrained) hydraulic boundary conditions. The general solutions are derived by first applying Laplace transforms with respect to the time and then solving the resulting governing equations in terms of Fourier–Bessel series, which involve trigonometric and hyperbolic functions with respect to the z-coordinate and Bessel functions with respect to the r-coordinate. Several time-dependent boundary-value problems are solved to demonstrate the application of the general solution to practical situations. Accuracy of the numerical solution is confirmed by comparing with the existing solutions for the limiting cases of a finite elastic cylinder and a poroelastic cylinder under plane strain conditions. Selected numerical results are presented for different cylinder aspect ratios, loading and hydraulic boundary conditions to demonstrate the key features of the coupled poroelastic response.     

Interaction between an acoustic wave and a disc

A numerical solution is obtained to the problem concerning a pressure measurement at the boundary between an ideal compressible fluid and a solid wall. It is assumedthat the fluid occupies a semiinfinite cylinder with a rigid bottom into which an elastic disc is inserted and heldfirmly around its edges. Motion is produced by a pressure wave originating at infinity. A finite-difference grid for this application is described and the results of actual calculations are shown.Deceased.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 84–91, January–February, 1972.The authors thank L. M. Flitman for reviewing the work.     

One-dimensional and two-dimensional analogs for three-dimensional viscous flows in the neighborhood of the plane of symmetry of a blunt body

An approximate method of determining the heat transfer and friction stress in three-dimensional flow problems using the two-dimensional and one-dimensional solutions is proposed. This method is applicable over a wide range of Reynolds numbers — from low to high. On the basis of a theoretical analysis of the approximate analytic solution of the equations of a three-dimensional viscous shock layer it is shown that the problem of determining the heat flux in the neighborhood of the plane of symmetry of bodies inclined to the flow at an angle of attack can be reduced, firstly, to the problem of determining that quantity for an axisymmetric body and, secondly, to the problem of determining the heat transfer to an axisymmetric stagnation point. On the basis of an analysis of the results of a numerical solution of the problem it is shown that corresponding analogs can also be used for the friction stress. The accuracy of the similarity relations established is estimated by solving the problem by a finite-difference method. A similarity relation of the same kind was previously obtained in [1] for a double-curvature stagnation point.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 117–122, January–February, 1990.     

Elastoplastic state of flexible cylindrical shells with two circular holes

The elastoplastic state of thin cylindrical shells with two equal circular holes is analyzed with allowance made for finite deflections. The shells are made of an isotropic homogeneous material. The load is internal pressure of given intensity. The distribution of stresses along the hole boundary and in the stress concentration zone (when holes are closely spaced) is analyzed by solving doubly nonlinear boundary-value problems. The results obtained are compared with the solutions that allow either for physical nonlinearity (plastic strains) or geometrical nonlinearity (finite deflections) and with the numerical solution of the linearly elastic problem. The stresses near the holes are analyzed for different distances between the holes and nonlinear factors.Translated from Prikladnaya Mekhanika, Vol. 40, No. 10, pp. 107–112, October 2004.     

A Pressure Field in an Acoustic Medium Containing a Cylindrical Solid and a Sphere Oscillating in a Predetermined Manner

A problem on the interaction of a spherical body oscillating in a predetermined fashion and a rigid cylinder is formulated. The bodies do not intersect, are immersed into an ideal compressible liquid, and their centers are in one plane. The solution is based on the possibility of representing the partial solution of the Helmholtz equation, written in cylindrical coordinates, in terms of partial solutions in spherical coordinates, and vice versa. An infinite system of linear algebraic equations is obtained by satisfying the boundary conditions on the sphere and cylinder surfaces. The system is intended for determining the coefficients of the expansion of the velocity potential into a series in terms of spherical and trigonometric functions. The system obtained is solved by the reduction method. The appropriateness of this method is substantiated. The hydrodynamic characteristics of the liquid surrounding the spherical and cylindrical bodies are determined. A comparison is made with the problem on a sphere oscillating in an infinite incompressible liquid that contains also a cylinder and in a compressible liquid that contains nothing more. Two types of motion of the sphere — pulsation and oscillation — are considered     

Influence of Anisotropy on the Response Characteristics of Finite Cylinders Under Free Vibrations

A numerical–analytical approach is proposed to solve a problem on the free vibrations of cylindrical bodies. The approach is based on three-dimensional elastic theory and the semianalytic finite-element method. The free vibrations of isotropic and anisotropic solid cylinders of finite length are examined. It is studied how boundary conditions and mechanical and geometrical parameters affect the distribution of dynamic properties. The efficiency of the approach proposed is tested by comparing results produced by different approaches     

Boundary element calculation of shallow water waves

A boundary element method is proposed for studying periodic shallow water problems. The numerical model is based on the shallow water equation. The key feature of this method is that the boundary integral equations are derived using the weighted residual method and the fundamental solutions for shallow water wave problems are obtained by solving the simultaneous singular equations. The accuracy of this method is studied for the wave reflection problem in a rectangular tank. As a result of this test, it has been shown that the number of element divisions and the distribution of nodes are significant to the accuracy. For numerical examples of external problems, the wave diffraction problems due to single cylindrical, double cylindrical and plate obstructions are analysed and compared with the exact and other numerical solutions. Relatively accurate solutions are obtained.     

Thermoelastic Contact of a Shroud Ring and a Cylinder under Unsteady Friction-Induced Heat Release

Mathematical formulation is performed and a solution is found for a quasi-static thermoelastic problem of contact interaction of an elastic shroud ring and a hollow circular cylinder inserted into this ring, which are compressed by a load varied along the axis of the system, under the condition of an unloaded contact over the ring surface or over the circumference contour. The radial displacements of the contact surface of the shroud ring are approximated by displacements of the surface of a long circular hollow cylinder. Unsteady friction-induced heat release caused by the action of friction forces owing to shroud ring rotation over the cylinder with a time-dependent low angular velocity is taken into account. The problem is reduced to a system of integral equations whose structure is determined by the form of thermophysical contact conditions. A numerical algorithm of the solution is proposed, and the influence of the problem parameters on the contact pressure and temperature distributions is considered. Based on an analysis of results, a conclusion is made that the character of axial variation of the compressing load has a significant effect on the distribution of contact pressure in describing the kinematic condition of interaction of bodies in accordance with Hertz’s theory.__________Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 4, pp. 161–178, July– August, 2005.