A unified Krylov-Bogoliubov method for solving second-order non-linear systems

A unified theory is presented for obtaining the transient response of second-order non-linear systems by the Krylov-Bogoliubov method. The method is a generalization of Bogoliubov's asymptotic method and covers all three cases when the roots of the corresponding linear equation are real, complex conjugate, or pure imaginary. It is shown that by suitable substitution for the roots in the general result, that the solution corresponding to each of the three cases can be obtained. The solution for the equation governing the motion of a simple pendulum with and without damping derived from the general solution reduces to that obtained by Popov's [4] method.     

Interface crack extension at any constant speed in orthotropic or transversely isotropic bimaterials––I. General exact solutions

A semi-infinite crack along the interface of two dissimilar half-spaces extends under in-plane loading. Each half-space belongs to a class of orthotropic or transversely isotropic elastic materials, the crack can extend at any constant speed, and all six possible relations between the four body wave speeds are considered. A steady dynamic situation is treated, and exact full displacement fields derived. A key step is a factorization that produces, despite anisotropy, simple solution forms and compact crack speed-dependent functions that exhibit the Rayleigh and Stoneley speeds as roots. These roots are calculated for various representative bimaterials.Closed-form crack opening displacement gradient and interface stress fields are also derived from a general set of coupled singular integral equations. The equation eigenvalues can, depending on crack speed, be complex/imaginary conjugates, purely real, or zero. This suggests possibilities observed in other studies: oscillations and square-root singular behavior at the crack edge, non-singular behavior, singular behavior not of square-root order, and the radiation of displacement gradient discontinuities at crack speeds beyond the purely sub-sonic range.These possibilities are explored further in terms of two important special cases in Part II of this study [Int. J. Solids Struct., 39, 1183–1198].     

Impact of a plane die on an elastic orthotropic half-plane

To study the process of impact of a rigid body on the surface of an elastic body made of a composite material, we consider a nonstationary dynamic contact problem about the impact of a plane rigid die on an elastic orthotropic half-plane. The problem is reduced to solving an integral equation of the first kind for the Laplace transform of the contact stresses under the die base. An approximate solution of the integral equation is constructed with the use of a special approximation to the symbol of the kernel of the integral equation in the complex plane. The inverse Laplace transform of the solution results in determining the scalar contact stress field on the die base, the force exerted by the die on the elastic medium, and the vertical displacement field of the free surface of the orthotropic medium out side the die. The solutions thus obtained permit studying specific features of the process of die penetration into an orthotropic medium and the strain properties of the medium.     

加权残数配点法解正交各向异性板的积分方程

本文推导了一般各向异性板弯曲的积分方程,运用加权残数配点法求解了正交各向异性板弯曲的积发方程,本文将部分配点取在边界上,另一部分配点取在域外,只用关于找度的基本积分方程,而不用关于转角的补充积分方程,简化了方程求解和计算程序,由于正交各向异性板没有争析形式的、实用的基本解,本文提出了两种新的近似基本解;加权双三角级数;广义各向同性板解析形式的基本解和加权双三角级数的叠加,算例表明,本文提出的解法和近似基本解适用于各类边界条件的正交各向异性板,具有简单、可靠、精度高等优点。     

Dynamic behavior of rectangular crack in three-dimensional orthotropic elastic medium by means of non-local theory

The dynamic behavior of a rectangular crack in a three-dimensional(3D)orthotropic elastic medium is investigated under a harmonic stress wave based on the non-local theory. The two-dimensional(2D) Fourier transform is applied, and the mixedboundary value problems are converted into three pairs of dual integral equations with the unknown variables being the displacement jumps across the crack surfaces. The effects of the geometric shape of the rectangular crack, the circular frequency of the incident waves,and the lattice parameter of the orthotropic elastic medium on the dynamic stress field near the crack edges are analyzed. The present solution exhibits no stress singularity at the rectangular crack edges, and the dynamic stress field near the rectangular crack edges is finite.     

Weight functions for interface cracks in dissimilar anisotropic materials

Bueckner‘s work conjugate integral customarily adopted for linear elastic materials is established for an interface crack in dissimilar anisotropic materials. The difficulties in separating Stroh‘s six complex arguments involved in the integral for the dissimilar materials are overcome and then the explicit function representations of the integral are given and studied in detail. It is found that the pseudo-orthogonal properties of the eigenfunction expansion form (EEF) for a crack presented previously in isotropic elastic cases, in isotopic bimaterial cases, and in orthotropic cases are also valid in the present dissimilar arbitrary anisotropic cases. The relation between Bueckner‘s work conjugate integral and the J-integral in these cases is obtained by introducing a complementary stressdisplacement state. Finally, some useful path-independent integrals and weight functions are proposed for calculating the crack tip parameters such as the stress intensity factors.     

Stationary vortex sheets in a stirring flow

Stationary vortex sheets in a two-dimensional stirring flow may be approximated by arrays of stationary point vortices arranged along the support of the sheets. These vortices lie at the roots of a polynomial that satisfies a generalized Lamé differential equation; the polynomial itself (not the roots) determines the complex potential and stream function. In this paper, sufficient conditions for the stirring flow are found so that the differential equation has two independent polynomial solutions with simple closed-form expressions, analogous to hypergeometric polynomials. The corresponding point vortex array then depends on a complex parameter that controls the location of the sheet, so that it may pass through any selected point. Stationary sheets in a periodic flow are approximated by the same method.     

Method for determining power-type complex singularities in solutions of singular integral equations with generalized kernels and complex conjugate unknowns

We develop a method for determining power-type complex singularities of solutions for a class of one-dimensional singular integral equations with generalized kernels and complex conjugate unknown functions. By analyzing the characteristic part of a singular integral equation, we reduce the problem of determining the solution singularity exponents at the ends of the integration interval to two independent transcendental equations for these exponents. We show that the distribution of admissible singularity exponents is of continuous character. We present numerical results for a two-dimensional elasticity problem whose mathematical statement leads to a singular integral equation of the class under study. We also reveal the drawbacks of one classical approach to the determination of stress field singularities.     

Application of Betti's reciprocal work theorem to the construction of the hypersingular integral equation of a plane crack in three-dimensional elasticity

The two-dimensional hypersingular integral equation of a plane crack under an arbitrary normal pressure distribution inside an infinite, three-dimensional, isotropic elastic medium is rederived here by application of the classical Betti's reciprocal work theorem. This approach is a very simple one and, moreover, it illustrates the usefulness of the aforementioned theorem for the derivation of hypersingular integral equations of elasticity problems.     

Extraordinary acoustic transmission,symmetry, blaschke products and resonators

The phenomenon of extraordinary acoustic transmission ( eat) in a resonator, which has recently been investigated experimentally, is studied theoretically. It is shown that the combination of a single propagating mode and a symmetry orthogonal to the direction of propagation for a resonator leads to eat. This is accomplished by decomposing the problem using symmetry, the Blaschke product and the properties of functions of a single complex variable which have modulus one on the real axis. The conditions of a resonator requires that the solution has singularities in the analytic extension to complex frequencies (resonances) and it is precisely near these resonances that we observe eat. The condition of a Blaschke product requires that there is a zero at the complex conjugate of the singularity and eat occurs when the solution on the real axis passes between these complex conjugate pairs of poles and zeros. A detailed numerical study of the problem is conducted and we show that once the single mode of propagation or the symmetry is broken then eat (at least perfect transmission) no longer holds generally.     

The stability of the interface between two bodies compressed along interface cracks. 3. Exact solutions for the combined case of equal and unequal roots

The problem of the stability of the interface between two bodies is considered for the case where several plane cracks are located in the interface, and the bodies are compressed along them (along the interface of two different materials). The study is carried out for a plane problem by using the three-dimensional linearized theory of stability of deformable bodies. Complex variables and potentials of the above-mentioned linearized theory are used. This problem is reduced to the problem of linear conjugation of two analytical functions of complex variable. The exact solution of the above-mentioned problem is derived for the case where the basic equation has unequal roots for the first material and equal roots for the second material. In earlier authors' publications, the exact solutions were obtained for the cases where both materials have either equal or unequal roots. Some mechanical effects are analyzed for the general formulation of the problem (elastic, elastoplastic compressible and incompressible isotropic and orthotropic bodies). It is pointed out that, in accordance with the exact solutions, the main result and conclusions have a general form for the above-mentioned cases of roots. S.P. Timoshenko Institute of Mechanics, National Academy of Sciences of Ukraine, Kiev. Translated from Prikladnaya Mekhanika, Vol. 36, No. 6, pp. 67–77, June, 2000     

Three-Dimensional Compositional Elastic Field in Cubic Anisotropic Solid Films on Compliant Substrates

A two-dimensional Laplace transform technique is used to find an analytical solution for the three-dimensional elastic field in solid films attached to a compliant substrate. The stress field arises from the misfit between the film and substrate and also from the non-uniform composition field in the film itself. Both the film and substrate are assumed to have cubic symmetry. The solution method can also be applied to materials with general orthotropic symmetry. Mathematics Subject Classifications (2000) 74F05, 74G10, 74K35, 74E10, 44A30     

The problem of the simple smooth crack in an infinite anisotropic elastic medium

The problem of the simple smooth curvilinear crack in an infinite anisotropic elastic medium under conditions of generalized plane stress or plane strain and under the supposition that the plane of the problem is a plane of elastic symmetry of the anisotropic medium is reduced to a complex Cauchy-type singular integral equation along the crack together with a condition of single-valuedness of displacements around the crack by using the complex potentials technique. Application to the case of a straight crack is also given.     

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法对二维对偶积分方程完成求解, 最终推导出板土系统在动力作用下的位移和应力解析式. 通过将本文计算模型退化为单一弹性地基, 与已有研究结果进行对比, 验证了本文方法的正确性和有效性. 最后, 通过数值算例, 对饱和多孔地基上矩形板的动力响应及参数影响做出分析和讨论. 此外, 本文提出的解析法具有一般性, 可广泛应用于复杂接触问题和多场耦合问题的求解.     

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

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

Research on an orthogonal relationship for orthotropic elasticity

IntroductionAsymplecticsystematicmethodology[1- 3]forelasticitywasestablishedbyZhongWan_xie .Hepresentedcreativelythedualvectorsandthesymplecticorthogonalrelationshipandopenedaworkplatformparalleledtothetraditionalelasticity[4 - 9].AnewdualvectorandanewdualdifferentialmatrixLwerepresentedforasymplecticsystematicmethodologyfortwo_dimensionalelasticityandaneworthogonalrelationshipwasdiscoveredforisotropicplaneproblems[4 ]byLuoJian_hui.Theneworthogonalrelationshipisgeneralizedfororthotropicelas…     

横观各向同性弹性半空间地基上正交异性矩形中厚板弯曲解析解

对横观各向同性体通解进行双重傅里叶变换,获得了直角坐标系下横观各向同性弹性半空间地基受任意竖向荷载作用下的位移积分变换解;在此基础上建立了板与地基的变形协调方程,并与三个广义位移变量描述的弹性地基上四边自由正交各向异性矩形中厚板的弯曲控制方程相结合,用三角级数法,得出横观各向同性弹性半空间地基上四边自由正交异性矩形中厚板受任意竖向荷载作用的弯曲解析解。相关算例分析表明,本文方法是有效的。     

复杂载荷三维裂纹分析双重边界元法

提出可用于高温、高转速状态下的热动力机械三维含裂构件热弹性分析方法——双重边界元法.首先建立了考虑温度及离心载荷的双重边界积分方程组,并对边界积分方程组的选取及适用范围进行了讨论。然后提出角非快调元模型离散技术。接着提出超奇异积分方程分析去除奇异性方法及数值积分技术.数值算例表明计算结果与有关权函数解十分吻合,说明了用双重边界元法计算复杂载荷条件下三维应力强度因子的有效性.还讨论了有关热应力强度因子权函数解的适用范围.     

Stress singularities for three-dimensional corners using the boundary integral equation method

The boundary integral equation method is developed to study three-dimensional asymptotic singular stress fields at vertices of a pyramidal notch or inclusion in an isotropic elastic space. Two-dimensional boundary integral equations are used for the infinite body with pyramidal notches and inclusions when either stresses or displacements are specified on its surface. Applying the Mellin integral transformation reduces the problem to one-dimensional singular integral equations over a closed, piece-wise smooth line. Using quadrature formulas for regular and singular integrals with Hilbert and logarithmic kernels, these integral equations are reduced to a homogeneous system of linear algebraic equations. Setting its determinant to zero provides a characteristic equation for the determination of the stress singularity power. Numerical results are obtained and compared against known eigenvalues from the literature for an infinite region with a conical notch or inclusion, for a Fichera vertex, and for a half-space with a wedge-shaped notch or inclusion.     

充液层合压电圆柱壳的自由振动

推导得到轴向极化的圆柱型正交各向异性压电弹性力学的三维状态方程，采用细分近似方法，得到了状态方程的解，并建立了圆柱壳内外表面边界量之间的传递关系，分析了内充可压缩流体的层合压电圆柱壳的自由振动问题，给出了频率方程的精确形式，并作了具体计算。