Solution of a multidimensional impact deformation problem for an elastic half-space with curved boundary on the basis of a modified ray method

A generalization of the method for constructing approximate solutions of boundary value problems of impact deformation dynamics in the form of ray expansions for two-dimensional plane deformation problems is presented. For each shock wave, the solution near its front is determined on the basis of ray coordinates consistent with this wave. The nonlinear divergence of curvilinear rays is taken into account. A mechanism of transformation from one ray coordinate system to another, which is crucially important in the ray method, is described. The developed technique is illustrated by solving the impact deformation problem for a half-space with boundary of nonzero curvature.     

Analytic Solution of an Inhomogeneous Kinetic Equation in the Problem of Second-Order Thermal Creep

An analytic solution of the problem of second-order thermal creep is obtained. A method for solving the half-space boundary value problem for an inhomogeneous linearized kinetic BGK equation forms the basis of the solution. The general solution of the input equation is constructed in the form of an expansion of the corresponding characteristic equation in terms of the eigenfunctions. Substitution of the solution in the boundary conditions leads to a Riemann boundary value problem. The unknown thermal creep velocity is found from the condition of solvability of the boundary value problem. The numerical analysis performed confirms the existence of negative thermophoresis (in the direction of the temperature gradient) for high-conductivity aerosol particles at low Knudsen numbers.     

Perturbation Method in the Problem of Compression-Shear Shock Load for a Nonlinear Elastic Half-Space

On the example of a one-dimensional nonstationary problem of oblique impact on the boundary of a nonlinear elastic isotropic half-space, the question of the manifestation of nonlinear deformation effects via basic evolution equations is studied. Much attention is given to the behavior of the solution behind the leading edge of a quasi-transverse shock wave. For particular cases of boundary conditions, it is shown that the onset region of the evolution equation of a quasi-transverse wave is preceded by a series of preliminary transitions to the intermediate internal problems of the small parameter method determined by the type of preliminary bulk deformation. This deformation consistently affects the distortion of the characteristic coordinates and the leading edge of the quasitransverse process. As a consequence, the transition to the evolution equation of quasi-transverse waves occurs with simultaneous change of all independent variables of the boundary value problem.     

NON-AXISYMMETRIC MIXED BOUNDARY VALUE PROBLEM FOR DYNAMIC INTERACTION BETWEEN SATURATED POROUS FOUNDATION AND RECTANGULAR PLATE$^{\bf 1)}$

在同一界面的不同区域具有多种边界条件, 称之为混合边界, 这是一个熟知的力学问题. 对这类问题进行精确分析时, 必须要进行混合边值问题的求解. 而对于一般的三维非轴对称情形, 混合边值问题的求解往往存在数学困难. 本文利用Hilbert定理和双重Fourier变换, 给出了一种求解三维非轴对称混合边值问题的解析方法, 利用该方法对具有混合透水边界的饱和多孔地基上矩形板的振动弯曲进行了解析研究(板与地基接触面为不透水边界, 其余为透水边界). 首先, 基于Kirchhoff理论和Biot多孔介质理论建立矩形板与饱和多孔地基的动力控制方程, 进行耦合求解. 针对板土接触面和非接触面的混合边值问题, 采用双重Fourier变换构造出两对二维对偶积分方程, 以接触应力和接触面孔隙压力为基本未知量, 用Jacobi正交多项式将未知量展开, 再利用Schmidt法对二维对偶积分方程完成求解, 最终推导出板土系统在动力作用下的位移和应力解析式. 通过将本文计算模型退化为单一弹性地基, 与已有研究结果进行对比, 验证了本文方法的正确性和有效性. 最后, 通过数值算例, 对饱和多孔地基上矩形板的动力响应及参数影响做出分析和讨论. 此外, 本文提出的解析法具有一般性, 可广泛应用于复杂接触问题和多场耦合问题的求解.     

Analysis of right closed thin-walled prismatic shells supported by stringers with geometric and physical nonlinearity taken into account

We consider thin-walled right-angle closed prismatic shells with rigid contour of the transverse cross-section. Such shells underlie the schemes used in the analysis of various thin-walled spatial structures. The use of nonlinear physical and geometric relations in the computations permits numerically obtaining the strength margin of the corresponding structures. In the present paper, we propose methods for obtaining a boundary value problem and analyzing such shells with nonlinear factors taken into account; the problem is presented as a system of linear differential equations with variable coefficients. We show that, within the approach proposed, this boundary value problem has a fixed structure independent of the special form of nonlinearity. The entire variety of problems of static analysis of right-angle prismatic shells with nonlinear factors taken into account can be reduced to solving this boundary value problem. Methods for taking a specific nonlinearity into account are treated as various methods for obtaining expressions for the variable coefficients in the matrices of the boundary value problem. We present methods for solving this boundary value problem numerically; these methods are independent of the specific form of the nonlinearity.     

饱和多孔地基与矩形板动力相互作用的非轴对称混合边值问题

在同一界面的不同区域具有多种边界条件, 称之为混合边界, 这是一个熟知的力学问题. 对这类问题进行精确分析时, 必须要进行混合边值问题的求解. 而对于一般的三维非轴对称情形, 混合边值问题的求解往往存在数学困难. 本文利用Hilbert定理和双重Fourier变换, 给出了一种求解三维非轴对称混合边值问题的解析方法, 利用该方法对具有混合透水边界的饱和多孔地基上矩形板的振动弯曲进行了解析研究(板与地基接触面为不透水边界, 其余为透水边界). 首先, 基于Kirchhoff理论和Biot多孔介质理论建立矩形板与饱和多孔地基的动力控制方程, 进行耦合求解. 针对板土接触面和非接触面的混合边值问题, 采用双重Fourier变换构造出两对二维对偶积分方程, 以接触应力和接触面孔隙压力为基本未知量, 用Jacobi正交多项式将未知量展开, 再利用Schmidt法对二维对偶积分方程完成求解, 最终推导出板土系统在动力作用下的位移和应力解析式. 通过将本文计算模型退化为单一弹性地基, 与已有研究结果进行对比, 验证了本文方法的正确性和有效性. 最后, 通过数值算例, 对饱和多孔地基上矩形板的动力响应及参数影响做出分析和讨论. 此外, 本文提出的解析法具有一般性, 可广泛应用于复杂接触问题和多场耦合问题的求解.      

Nonlinear dynamical characteristics of piles under horizontal vibration

The pile-soil system is regarded as a visco-elastic half-space embedded pile. Based on the method of continuum mechanics, a nonlinear mathematical model of pile-soil interaction was established-a coupling nonlinear boundary value problem. Under the case of horizontal vibration, the nonlinearly dynamical characteristics of pile applying the axis force were studied in horizontal direction in frequency domain. The effects of parameters, especially the axis force on the stiffness were studied in detail. The numerical results suggest that it is possible that the pile applying an axis force will lose its stability. So,the effect of the axis force on the pile is considered.     

The general and basic boundary value problems of plane steady jet flows and the corresponding nonlinear equations

The general boundary value problem, including known plane steady jet flows of an ideal incompressible fluid, is formulated. The simplest problem retaining all the specific features of the general problem, known as the basic problem, is separated from the general problem. The solution of the basic problem is reduced to solving a non-linear integro-differential equation and also to solving nonlinear integral equations. Examples of flows whose determination is reduced to' solving the basic problem are cited.     

Seismic oscillations of an elastic rectangular structure protected by a viscoelastic stratum: in-plane problem

We study the problem of seismic protection of a rectangular elastic construction by mean of a viscoelastic stratum made of material of Kelvin–Voigt type. In particular we consider the in-plane problem. By using the standard boundary integral equations method (BIEs), the dynamic problem is reduced to a system of BIEs over the set of boundary lines. In order to reduce the dimension of the problem in its discrete form, for wave processes in the stratum and in the half-space foundation we use a special Green’s function satisfying the stress-free boundary conditions over the boundary of the half-space. A numerical algorithm is used to solve the problem after discretization and finally we discuss the physical meaning of the results and the efficiency of the seismic protection by the viscoelastic stratum.     

An analytical solution of an axisymmetric problem of deforming an isotropic half-space with an elastic fixed boundary under the action of a distributed load

An analytical solution to the axisymmetric problem on the action of a distributed load on an isotropic half-space when the load is given by a function dependent on the radial coordinate is obtained. The surface of the half-space is elastically fixed outside the circular domain of load application, the shear stresses are absent along the entire boundary, and the stresses vanish at infinity. At the boundary and inside the elastic half-space, the solutions are represented by the formulas for the stress tensor components and for the displacement vector components.     

Mechanics of indentation for piezoelectric thin films on elastic substrate

Frictionless normal indentation problem of rigid flat-ended cylindrical, conical and spherical indenters on piezoelectric film, which is either in frictionless contact with or perfectly bonded to an elastic half-space (substrate), is investigated. Both conducting and insulating indenters are considered. With Hankel transform, the general solutions of the homogeneous governing equations for the piezoelectric layer and the elastic half-space are presented. Using the boundary conditions for a vertical point force or a point electric charge, and the boundary conditions on the film/substrate interface, the Green’s functions can be obtained by solving sets of simultaneous linear algebraic equations. The solution of the indentation problem is obtained by integrating these Green’s functions over the contact area with unknown surface tractions or electric charge distribution, which will be determined from the boundary conditions on the contact surface between the indenter and the film. The solution is expressed in terms of dual integral equations that are converted to a Fredholm integral equation of the second kind and solved numerically. Numerical examples are also presented. The comparison between two film/substrate bonding conditions is made. It shows that the indentation rigidity of the film/substrate system is lower when the film is in frictionless contact with the substrate. The effects of the Young’s modulus and Poisson’s ratio of the elastic substrate, indenter electrical condition and indenter prescribed electric potential on the indentation responses are presented.     

半空间弹性波双参数强非线性逆散射

半空间弹性波强散射源逆问题是一个实用而又难度很大的强非线性问题,而经典的Born迭代方法只适用于弱非线性问题．基于DBI方法（Distorted-Born Iterative方法）,提出了一种求解弹性波强非线性逆散射问题的迭代方法．在数值模拟运算时利用矩量法进行离散处理,并采用正则化原则避免求解病态矩阵方程。计算结果表明该方法具有较快的收敛速度及较高的精度。     

Third-order nonlinear singularly perturbed boundary value problem

ＩｎｔｒｏｄｕｃｔｉｏｎＳｉｎｃｅ 1 958,Ｌｅｖｓｉｏｎａｎｄｓｅｖｅｒａｌｏｔｈｅｒａｕｔｈｏｒｓ[1]- [4 ]ｈａｖｅｍａｄｅａｓｕｂｓｔａｎｔｉａｌｓｔｕｄｙｆｏｒｅｘｉｓｔｅｎｃｅａｎｄｅｓｔｉｍａｔｅｏｆｓｏｌｕｔｉｏｎｓｆｏｒｔｈｅｔｈｉｒｄ＿ｏｒｄｅｒｓｉｎｇｕｌａｒｌｙｐｅｒｔｕｒｂｅｄｂｏｕｎｄａｒｙｖａｌｕｅｐｒｏｂｌｅｍｗｉｔｈｔｗｏｐｏｉｎｔｂｏｕｎｄａｒｙｃｏｎｄｉｔｉｏｎｓｘ( 0 ) =Ａ ,ｘ′( 0 ) =Ｂ ,ｘ′( 1 ) =Ｃ ,ｏｒｓｌｉｇｈｔｌｙｅｘｔｅｎｓｉｖｅｌｙｌｉｎｅａｒｂｏ…     

非线性大变形问题边界元分析的自校正方法

针对用增量法求解非线性方程解的漂移问题，在非线性问题边界元法计算中建立了自我校正方法，对在拖带坐标上建立的增量形式的基本方程，引入Ｌａｎｇｒａｎｇｅ校正因子，以全量形式的基本方程作为其辅助方程，在此基础上导出含校正项的边界积分方程，边界元自我校正方法的建立有效地保证了在非线性问题的计算中最终收敛在其解附近，提高了计算精度和运算效率。     

Solution of the optimal turn problem for a spherically symmetric rigid body with arbitrary boundary conditions in the class of generalized conical motions

We consider the problem of time- and energy consumption-optimal turn of a rigid body with spherical mass distribution under arbitrary boundary conditions on the angular position and angular velocity of the rigid body. The optimal turn problem is modified in the class of generalized conical motions, which allows one to obtain closed-form solutions for equations of motion with arbitrary constants. Thus, solving the optimal control boundary value problem is reduced to solving a system of nonlinear algebraic equations for the constants. Numerical examples are considered to illustrate the proximity between the solutions of the traditional and modified problems of optimal turn of a rigid body.     

分层介质中Love波的扩展W-W算法

应用精细积分法(PIM)和扩展Wittrick-Williams(W-W)算法求解横观各向同性分层半空间中的Love波问题.Love波对应于波数-频率域线性常微分方程的本征值问题.精细积分法是求解线性常微分方程两端边值问题和初值问题的高精度算法.利用本征值计数技术,扩展W-W算法可以不遗漏地找到所有本征值.因此,文中使用的方法可以得到计算机精度意义下的精确解.     

Action of an elliptic punch on a transversally isotropic half-space

The 3D contact problem on the action of a punch elliptic in horizontal projection on a transversally isotropic elastic half-space is considered for the case in which the isotropy planes are perpendicular to the boundary of the half-space. The elliptic contact region is assumed to be given (the punch has sharp edges). The integral equation of the contact problem is obtained. The elastic rigidity of the half-space boundary characterized by the normal displacement under the action of a given lumped force significantly depends on the chosen direction on this boundary. In this connection, the following two cases of location of the ellipse of contact are considered: it can be elongated along the first or the second axis of Cartesian coordinate system on the body boundary. Exact solutions are obtained for a punch with base shaped as an elliptic paraboloid, and these solutions are used to carry out the computations for various versions of the five elastic constants. The structure of the exact solution is found for a punch with polynomial base, and a method for determining the solution is proposed.     

Contact problem for an orthotropic half-space

Numerical and analytical solutions of the 3D contact problem of elasticity on the penetration of a rigid punch into an orthotropic half-space are obtained disregarding the friction forces.A numericalmethod ofHammerstein-type nonlinear boundary integral equations was used in the case of unknown contact region, which permits determining the contact region and the pressure in this region. The exact solution of the contact problem for a punch shaped as an elliptic paraboloid was used to debug the program of the numerical method. The structure of the exact solution of the problem of indentation of an elliptic punch with polynomial base was determined. The computations were performed for various materials in the case of the penetration of an elliptic or conical punch.     

Stress disturbances caused by the inhomogeneity in an elastic half-space subjected to contact loading

Inhomogeneities can increase localized stress and cause microstructural alterations to initiate fatigue failures in rolling elements under cyclic contact loading. To study the stress disturbances created by the inhomogeneity, a two-dimensional contact stress analysis is presented for a cylindrical indenter sliding on an elastic half-space containing an inhomogeneity of arbitrary shape. The boundary element method is used to analyze the contact problem, where actual contact boundary, contact pressure as well as tractions and displacements at inhomogeneity–substrate interface are determined by solving a set of integral equations numerically. Numerical results are presented to investigate effects and the stress disturbances caused by the inhomogeneity with various locations, sizes and material properties of inhomogeneity. The results also show that hard inclusions are more detrimental than soft deformable particles in rolling contact elements.     

A boundary-field equation method for a nonlinear exterior elasticity problem in the plane

This paper presents a new method for solving a nonlinear exterior boundary value problem arising in two-dimensional elasto-plasticity. The procedure is based on the introduction of a sufficiently large circle that divides the exterior domain into a bounded region and an unbounded one. This allows us to consider the Dirichlet-Neumann mapping on the circle, which provides an explicit formula for the stress in terms of the displacement by using an appropriate infinite Fourier series. In this way we can reduce the original problem to an equivalent nonlinear boundary value problem on the bounded domain with a natural boundary condition on the circle. Hence, the resulting weak formulation includes boundary and field terms, which yields the so called boundary-field equation method. Next, we employ the finite Fourier series to obtain a sequence of approximating nonlinear problems from which the actual Galerkin schemes are derived. Finally, we apply some tools from monotone operators to prove existence, uniqueness and approximation results, including Cea type error estimates for the corresponding discrete solutions.