A 2.5-D dynamic model for a saturated porous medium. Part II: Boundary element method

The 3-D boundary integral equation is derived in terms of the reciprocal work theorem and used along with the 2.5-D Green’s function developed in Part I [Lu, J.F., Jeng, D.S., Williams, S., submitted for publication. A 2.5-D dynamic model for a saturated porous medium: Part I. Green’s function. Int. J. Solids Struct.] to develop the 2.5-D boundary integral equation for a saturated porous medium. The 2.5-D boundary integral equations for the wave scattering problem and the moving load problem are established. The Cauchy type singularity of the 2.5-D boundary integral equation is eliminated through introduction of an auxiliary problem and the treatment of the weakly singular kernel is also addressed. Discretisation of the 2.5-D boundary integral equation is achieved using boundary iso-parametric elements. The discrete wavenumber domain solution is obtained via the 2.5-D boundary element method, and the space domain solution is recovered using the inverse Fourier transform. To validate the new methodology, numerical results of this paper are compared with those obtained using an analytical approach; also, some numerical results and corresponding analysis are presented.     

三维动态裂纹问题的超奇异积分方程法

基于弹性材料的动态基本方程,结合广义Betti-Rayleigh互易等式与时域下的边界积分方程,推导得到时域下的超奇异积分方程组。引入Laplace域下的动态基本解,将经过主部分析的积分核函数分解为静态和动态部分,其中动态积分核不具有奇异性。在裂纹前沿附近单元,采用与理论分析一致的平方根位移模型。结合Lubich时间卷积实现拉氏变换,采用配置点法计算超奇异积分,获得问题的数值解。并针对椭圆裂纹算例编写Fortran程序,得到冲击荷载作用下张开型裂纹的动态应力强度因子变化规律,数值结果稳定且收敛速度快。     

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

结构粘弹性分析的线法

本文将以现代常微分方程求解器为支撑软件的线法推广应用于含时间变量的粘弹性力学问题。基基本思想是先用拦普拉斯积分变换把问题的控制偏微分方程组和定解条件进行交换,得到拉氏变换空间中的定常问题,然后用传统的线法借助于常微分方程求解器得到该变换空间中的数值解,再应用拉氏数值逆变换求得原时空域中的数值解。文中还给出方法的基本原理和实施过程。算例表明,采用拉普拉斯积分变换于线法分析结构粘弹性力学等含时间变量的     

Separation of principal stresses in dynamic photoelasticity

A new and effective method used to separate the transient principal stresses for dynamic photoelasticity is proposed. This is a hybrid method combining the optical method of dynamic caustics and the boundary element numerical method. Firstly, a modified Cranz-Schardin spark camera is used to record simultaneously the isochromatic fringe patterns of photoelasticity and the shadow spot patterns in the dynamic process. By means of the isochromatic fringe patterns, the difference between transient principal stresses in the whole domain and the principal stresses along the free boundary can be solved. In addition, the method of caustics is a very powerful technique for measuring the concentrative load. Then, the sum of the principal stresses is calculated by the boundary integral equation obtained from the Laplace integral transform of the wave equation. So, the transient principal stresses can be determined from the experimental and numerical results. As an example, the transient principal stresses in a polycarbonate disk under an impact load are resolved. Concurrently published in the Chinese Edition of Acta Mechanica Sinica, Vol. 26, No. 1, 1994     

Green quasifunction method for vibration of simply-supported thin polygonic plates on Pasternak foundation

A new numerical method—Green quasifunction is proposed.The idea of Green quasifunction method is clarified in detail by considering a vibration problem of simply-supported thin polygonic plates on Pasternak foundation.A Green quasifunction is established by using the fundamental solution and boundary equation of the problem. This function satisfies the homogeneous boundary condition of the problem.The mode shape differential equation of the vibration problem of simply-supported thin plates on Pasternak foundation is reduced to two simultaneous Fredholm integral equations of the second kind by Green formula.There are multiple choices for the normalized boundary equation.Based on a chosen normalized boundary equation,a new normalized boundary equation can be established such that the irregularity of the kernel of integral equations is overcome.Finally,natural frequency is obtained by the condition that there exists a nontrivial solution in the numerically discrete algebraic equations derived from the integral equations.Numerical results show high accuracy of the Green quasifunction method.     

求解饱和半空间上弹性圆板固结沉降的积分方程

本文采用解析方法分析了弹性圆板在饮和半空间上的固结沉降。考虑弹性圆板与饮和半空间的接触面上无摩擦力,且饱和半空间表面为全部透水的。运用Ｂｉｏｔ固结理论和积分方程技术,在Ｌａｐｌａｃｅ变换域上建立了弹性圆板固结沉降的对偶积分方程,并化此对偶积分方程为第二类Ｆｒｅｄｈｏｌｍ积分方程。通过对其核函数的有效数值发得到第二类Ｆｒｅｄｈｏｌｍ积分方程的解,再利用Ｌａｐａｃｅ反演技术获得弹性板在时间域中的固结沉     

Laplace transform-boundary element coupling method for viscous fluid-structure impact analysis

Based on linearized 2-D Navier-Stokes equation, a Laplace transform-boundary element coupling method for viscous fluid-structure impact analysis is proposed. Under assumption of incompressibility for the fluid, the corresponding equivalent boundary integral equation in terms of the potential function and stream function is first established by Lamb's transform in the Laplace transform domain. It enables us to solve impact water problems in frequency domain by the boundary element method, in which the effect of viscous flow on the dynamic response can be taken into account. Then a complete solution of the problem under consideration in time domain is obtained by means of Durbin's formulas for the numerical inversion of the Laplace transform. Finally, a practical example is given to confirm the validity of the present method. Project supported by the National Defence Foundation of Science & Technology of China (No. J14. 8. 1. JW0515).     

Free vibration analysis of a plate on foundation with completely free boundary by finite integral transform method

The theoretical solutions of eigenfrequencies and vibration modes of a rectangular thin plate on an elastic foundation with completely free boundary are derived by using a double finite cosine integral transform method. In the analysis procedure, the elastic foundation is regarded as a Winkler elastic foundation model. Because the basic dynamic elasticity equations of the thin plate on elastic foundation are only used, it is not needed to select the deformation function arbitrarily. Therefore, the solution developed in the present paper is more reasonable and more accurate. To prove the correctness of the solutions, numerical results obtained using the present solutions are compared with those in the literatures.     

Dynamic response of an elastic medium containing crack

A fundamental solution for an infinite elastic medium containing a penny-shaped crack subjected to dynamic torsional surface tractions is attempted. A double Laplace–Hankel integral transform with respect to time and space is applied both to motion equation and boundary conditions yielding dual integral equations. The solution of the derived dual integral equations is based on an analytic procedure using theorems of Bessel functions and ordinary differential equations. The dynamic displacements’ field is obtained by inversion of the corresponding Laplace–Hankel transformed variable. Results of a representative example for a crack subjected to pulse surface tractions are obtained and discussed.     

具有混合边界透水条件的多孔饱和半空间上刚性圆板的垂直振动

用积分变换和积分方程研究多孔饱和半空间上刚性圆板的垂直振动问题。首先应用逐次解耦方法求解多孔饱和固体的动力基本方程－Ｂｉｏｔ波动方程。然后考虑混合边界透水条件（半空间表面与圆板的接触面是不透水的,而其余表面是透水的）,建立子多孔饱和半空间上刚性圆板垂直振动的对偶积分方程,并化对偶积分方程为第二类Ｆｒｄｄｈｏｌｍ积分方程。     

Application of boundary integral equations for analyzing the dynamics of elastic,viscoelastic, and poroelastic bodies

Two approaches (classical and nonclassical) of the boundary integral equation method for solving three-dimensional dynamical boundary value problems of elasticity, viscoelasticity, and poroelasticity are considered. The boundary integral equation model is used for porous materials. The Kelvin–Voigt model and the weakly singular hereditary Abel kernel are used to describe the viscoelastic properties. Both approaches permit solving the dynamic problems exactly not only in the isotropic but also in the anisotropic case. The boundary integral equation solution scheme is constructed on the basis of the boundary element technique. The numerical results obtained by the classical and nonclassical approaches are compared.     

敷设多孔吸声材料声腔特征值分析的径向积分边界元法

由于Helmholtz方程的基本解是频率的函数,因此传统边界元法在处理声场特征值问题时具有天生的缺陷。本文采用Laplace方程基本解生成积分方程,通过径向积分法将在此过程中产生的域积分项转化为边界积分。此方法克服了传统边界元法系数矩阵对频率的依赖,同时克服了特解积分法对特解的依赖,并通过对表面声导纳的多项式逼近,将敷设多孔吸声材料声腔特征值问题转化为矩阵多项式,从而避免了复杂的非线性求解。通过数值算例验证了算法的有效性。     

黏弹性体界面裂纹的冲击响应

研究两半无限大黏弹性体界面Griffith裂纹在反平面剪切突出载荷下,裂纹尖端动应力强度因子的时间响应,首先,运用积分变换方法将黏弹性混合黑社会问题化成变换域上的对偶积分方程,通过引入裂纹位错密度函数进一步化成Cauchy型奇异积分方程,运用分片连续函数法数值求解奇异积分方程,得到变换域内的动应力强度因子,再用Laplace积分变换数值反演方法,将变换域的解反演到时间域内,最终求得动应力强度因子的时间响应,并对黏弹性参数的影响进行分析。     

Transient dynamic responses of a cracked solid subjected to in-plane loadings

Transient dynamic responses of an elastic cracked solid subjected to in-plane surface loadings are investigated in this study. Two vertical cracks, a surface-breaking crack and a sub-surface crack, are considered. The frequency responses of the plane strain problem are calculated by the computational mechanics combining the finite element method with the boundary integral equation. The finite element method is used for the near-field enclosing the crack, while the boundary integral equation is applied for the far-field to satisfy the Sommerfeld radiation condition. The transient responses are then obtained using fast Fourier transform. Surface displacements, crack opening displacements, and dynamic stress intensity factors are presented to show the significant effects of the cracks. The interaction between the elastic waves and the cracks as well as the mode conversion phenomena can be observed and understood through the numerical simulations.     

变厚度板弯曲的边界元分析法

由于变厚度板弯曲问题的控制分方程复杂，直接求解其基本解推导边界积分方程建立边界元分析法较为困难，本文通过引入等效荷载，等效刚度，将此问题的控制微分方程化成与普通薄板弯曲问题基本方程相同的形式，利用求解通板弯曲问题的边界元迭代求解，建立了分析变厚度板弯曲问题的蛤法，算例表明本方法理正确，精度良好。     

Waves generated by an inclined-plate wave generator

This paper describes the characteristics of small-amplitude waves generated by a sinusoidally oscillating, inclined paddle-type wavemaker operating in a constant-depth channel. Two-dimensional, linearized potential flow is assumed. A semi-analytical method, the boundary collocation method, is used to establish the relationship between wave amplitude and paddle stroke. The numerical results are compared with the numerical results of the boundary integral equation method. It is found that the boundary collocation method is simpler and more flexible to implement and faster to compute. In addition, the numerical results are in reasonably good agreement with the laboratory experimental data. For the vertical wavemaker, which is a special case of the inclined wavemaker, an analytical series solution can be found. By using the boundary collocation method and the boundary integral equation method to solve the vertical wavemaker problem and comparing the results with the analytical series solution, it is found that the boundary collocation method yields a solution which is much more accurate than that from the boundary integral equation method. Finally, the relationships between wave amplitude and paddle stroke are established for different inclinations of the paddle-type wavemaker, based on the boundary collocation method.     

Quasi-Green’s function method for free vibration of clamped thin plates on Winkler foundation

The quasi-Green’s function method is used to solve the free vibration problem of clamped thin plates on the Winkler foundation. Quasi-Green’s function is established by the fundamental solution and the boundary equation of the problem. The function satisfies the homogeneous boundary condition of the problem. The mode-shape differential equation of the free vibration problem of clamped thin plates on the Winkler foundation is reduced to the Fredholm integral equation of the second kind by Green’s formula. The irregularity of the kernel of the integral equation is overcome by choosing a suitable form of the normalized boundary equation. The numerical results show the high accuracy of the proposed method.     

The Flow Analysis of Viscoelastic Fluid with Fractional Order Derivative in Horizontal Well

The fractional calculus approach is introduced into the seepage mechanics. A three-dimensional relaxation model of viscoelastic fluid is built. The models based on four boundary conditions of exact solution in Laplace space for some unsteady flows in an infinite reservoir is obtained by using the Laplace transform and Fourier sine and cosine integral transform. The pressure transient behavior of non-Newtonian viscoelastic fluid is studied by using Stehfest method of the numerical Laplace transform inversion and Gauss–Laguerre numerical integral formulae. The viscoelastic fluid is very sensitive to the order of the fractional derivative. The change rules of pressure are discussed when the parameters of the models change. The plots of type pressure curves are given, and the results can be provided to theoretical basis and well-test method for oil field.