Scattering of flexural waves and dynamic stress concentrations in Mindlin's thick plates with a cutout

Using the complex variable method and conformal mapping, scattering of flexural waves and dynamic stress concentrations in Mindlin's thick plates with a cutout have been studied. The general solution of the stress problem of the thick plate satisfying the boundary conditions on the contour of cutouts is obtained. Applying the orthogonal function expansion technique, the dynamic stress problem can be reduced into the solution of a set of infinite algebraic equations. As examples, numerical results for the dynamic stress concentration factor in Mindlin's plates with a circular, elliptic cutout are graphically presented in sequence. The project supported by the National Natural Science Foundation of China     

含孔软铁磁材料Mindlin板中弹性波散射问题

基于考虑磁弹相互作用的Mindlin板弯曲波动方程，采用波函数展开法，分析研究 了含孔软铁磁材料Mindlin板中弹性波散射与动应力集中问题，给出了问题的分析 解和数值算例. 通过分析发现：磁感应强度对动弯矩集中系数和动剪力集中系数有 增加的作用，特别是在低频的情况下.     

弹性薄板弯曲问题的双互易杂交边界点法

基于一种板的修正变分泛函,将杂交边界点法与双互易法结合,用于薄板弯曲问题的分析。该方法将问题的解分为齐次方程的通解和非齐次的特解两部分,特解采用径向基函数插值得到,而通解则使用杂交边界点法求解。在杂交边界点法用于求解通解的列式过程中,边界变量采用移动最小二乘近似,域内变量则采用基本解插值。与有限元法相比,该方法仅需要边界上离散点的信息,无论插值还是积分都不需要网格,域内点仅用来插值非齐次项,因而仍是一种纯边界类型的无网格方法。数值算例表明,本文方法能以很少的计算自由度获得与其它方法同样的计算精度,且具有前后处理简单、收敛速度快等优点,适合于求解工程中各种薄板的弯曲问题。     

Dynamics stress concentrations in Ambartsumian's plate with a cutout

Based on the motion equations of flexural wave in Ambartsumian's plates including the effects of transverse shear deformations, by using perturbation method of small parameter, the scattering of flexural waves and dynamic stress concentrations in the plate with a cutout have been studied. The asymptotic solution of the dynamic stress problem is obtained. Numerical results for the dynamic stress concentration factor in Ambartsumian's plates with a circular cutour are graphically presented and discussed. The project supported by the National Natural Science Foundation of China.     

Multiple scattering of flexural waves in a semi-infinite thin plate with a cutout

The multiple scattering of flexural waves and dynamic stress concentration in a semi-infinite thin plate with a cutout are investigated, and the expressions of this problem are obtained. The analytical solutions of wave fields are expressed by employing the wave function expansion method and the expanded mode coefficients are solved by satisfying the boundary condition of the cutout. The image method is used to satisfy the traction free boundary condition of the plate. As an example, the numerical results of dynamic stress concentration factors are graphically presented and discussed. Numerical results show that the analytical results of the scattered waves and dynamic stress in semi-infinite plates are significantly different from those in infinite plates when the ratio of distance b/a is relatively little. In the region of low frequency and long wavelength, the maximum dynamic stress concentration factors occur on the illuminated side of the scattering body with θ = π, but not at the edge of the cutout with θ = π/2. As the incidence frequency increases (the wavelength becomes short), the dynamic stress on the illuminated side of the cutout decreases, however, the dynamic stress on the shadow side increases.     

Dual reciprocity hybrid radial boundary node method for Winkler and Pasternak foundation thin plate

An efficient dual reciprocity hybrid radial boundary node method is developed for the analysis of Winkler and Pasternak foundation thin plate, in which a hybrid displacement variational principle, radial point interpolation method (RPIM) and dual reciprocity method (DRM) are combined. Firstly, the hybrid displacement variational principle is developed, in which the domain variables are interpolated by two groups of symmetric fundamental solutions, while the boundary variables are interpolated by RPIM instead of the traditional moving least square, and the shape function obtained by RPIM satisfies the delta function property, so boundary conditions can be applied directly. Besides, DRM is exploited to evaluate the particular solutions of inhomogeneous terms, which can be used to transform the domain integrals arising from the inhomogeneous term into equivalent boundary integrals. Finally, some additional equations based on the DRM theory are proposed to overcome the problem that the boundary integral equations are not enough to solve all variables. This method has the advantages of both no element mesh of meshless method and dimensionality reduction of boundary element method. Numerical examples of Winkler and Pasternak foundation plates are given to illustrate that the present method is effective, accurate and it can be further expanded into practical engineering.     

A boundary element method for calculating elastic waves scattered by axisymmetric inclusion

A Boundary Element Method (BEM) is described to compute the scattering of elastic waves by an axisymmetric inclusion in an infinite elastic medium. The boundary loads applied to the inclusion is expanded in terms of Fourier series in an infinite space. The boundary integral equation is solved in the general direction of the axisymmetric inclusion by BEM. The problem of the 3-D scattering of elastic waves is reduced to a 1-Done. According to the geometric features of the axisymmetric in clusion the ring shell elements are adopted in this method. A comparison is made with other BEM methods. The numerical results show this method can reduce the amount of calculation and enhance the speed of convergence. Supported by Foundation of Ph. D Program of State Education Commission of China     

薄板小波有限元理论及其应用

利用样条小波尺度函数构造了常用的三角形和矩形薄板单元的位移函数,得到了利用小波函数表示的形函数。采用合理的局部坐标,对单元进行压缩,使单元在局部坐标区间上有其值,成功地推导出了分域的三角形和矩形薄板小波有限元列式。在此基础上,提出了弹性地基薄板的小波有限元求解方法。通过两个算例对薄板的挠度和弯矩进行了计算,数值结果表明,求解结果具有收敛快、精度高的特点。     

Zienkiewicz薄板单元的数值流形法实现

Zienkiewicz薄板单元是最早出现的板单元,因节点参数简单,所以在70年代之前曾起过重要作用;但同时它又是非协调的,仅对由三组间距相等的平行直线所形成的网格才能满足分片检验,这限制了这一单元的应用.因此,先后有不少学者对其进行过改进,使其能对任意形状的网格收敛.考虑到数值流形元法(NMM)总是用最佳质量的数学网格来进行逼近,本文尝试用NMM来解决Zienkiewicz单元的网格依赖问题.基于板的最小势能原理,提出了适用于NMM的混合变分提法,然后利用有解析解的椭圆板问题对方法进行了验证.本文所提出的方法,还可被用于其他类型的非协调单元的改造.     

各向异性平板开孔动应力集中问题的研究

采用各向异性平板弯曲波动理论及摄动方法，对正交各向异性平板开孔弯曲波的散射及动应力集中问题进行了分析研究，得到了此种平板稳态弯曲波动问题的渐近形式的分析解。同时采用保角映射技术，为求解正交各向异性平板开孔弹性波的散射及动应力集中问题提供了一种统一规范的方法。     

Research on the companion solution for a thin plate in the meshless local boundary integral equation method

The meshless local boundary integral equation method is a currently developed numerical method, which combines the advantageous features of Galerkin finite element method ( GFEM ), boundary element method (BEM) and element free Galerkin method (EFGM), and is a truly meshless method possessing wide prospects in engineeringapplications. The companion solution and all the other formulas required in the meshless local boundary integral equation for a thin plate were presented, in order to make this method apply to solve the thin plate problem.     

Multiple scattering of flexural waves by random configurations of inclusions in thin plates

Flexural waves are scattered by inclusions in a thin plate. For a single inclusion of arbitrary shape, reciprocity relations are obtained connecting coefficients in circular multipole expansions. Then, a formula for the effective wavenumber in a random arrangement of identical circular inclusions is derived, using the Lax quasi-crystalline approximation.     

A boundary element method for calculating the SH waves scattering from arbitrarily shaped holes in an anisotropic medium

The scattering problem of elastic wave by arbitrarily shaped cavities in an infinite anisotropic medium is investigated by the boundary integral equation (BIE) method. The formulations of BIE are derived with the help of generalized Green's formula. The discretization of BIE is based upon constant elements. After confirmation of the accuracy of the present method, some numerical examples are given for various cavities in a full space, in which an isotropic body with a circular cylinder hole is used for comparison and good agreement is observed. It has been proved that the method developed in this paper is effective.     

螺栓松动边界下纤维增强复合薄板的振动测试方法

提出了基于预埋压力传感器的量化测试方法,研究了螺栓松动边界对纤维增强复合薄板振动特性的影响。首先,自主设计并开发了带有预埋压力传感器的螺栓松动边界下复合薄板的振动测试系统,并详细介绍了系统各个部件的组成和功能;然后,归纳出一套合理、规范的松动边界下复合薄板的振动测试流程,并对HF10碳纤维/树脂复合薄板进行了实际测试。结果表明:随着螺栓松动程度的不断增加,复合薄板的固有频率逐渐降低,模态振型的节线位置也发生了不同程度的变化,但其阻尼结果呈现先增大后减小的趋势;而共振和非共振响应呈现先减小后增大的趋势。     

弯曲波对含多圆孔薄板的散射与动应力集中

采用复变函数法和多极坐标方法,研究了弯曲波对含有多圆孔薄板的散射问题。通过板的弯曲波动方程和内力方程的推导,求出在入射弯曲波条件下该问题的一般解的函数逼近序列和边界条件的表达式。用展开正交函数的方法将待解的问题归结为对一组无穷代数方程组的求解。最后,给出了含3圆孔薄板的孔边动应力集中系数的结果,并分析了孔间距和波数对动应力分布的影响。     

一种适用于简支多边形薄板问题的有限单元法

简支多边形薄板双调和定解问题可分解为二个互不耦合的Ｐｏｉｓｓｏｎ方程的定解问题，通过比拟，可用二维平面应力（应变）单元进行有限元分析。数值算例表明，本文采用的方法在较粗的网格下即可得到较高的精度，而且具有速度快、占用计算机内存少等优点。     

Local Petrov-Galerkin method for a thin plate

The meshless local Petrov-Galerkin (MLPG) method for solving the bendingproblem of the thin plate were presented and discussed. The method used the moving least-squares approximation to interpolate the solution variables, and employed a local symmetricweak form. The present method was a truly meshless one as it did not need a finite elementor boundary element mesh, either for purpose of interpolation of the solution, or for theintegration of the energy. All integrals could be easily evaluated over regularly shapeddomains ( in general, spheres in three-dimensional problems ) and their boundaries. Theessential boundary conditions were enforced by the penalty method. Several numericalexamples were presented to illustrate the implementation and performance of the presentmethod. The numerical examples presented show that high accuracy can be achieved forarbitrary grid geometries for clamped and simply-supported edge conditions. No postprocessing procedure is required to computer the strain and stress, since the originalsolution from the present method, using the moving least squares approximation, is already smooth enough.     

Natural convection in porous media—dual reciprocity boundary element method solution of the Darcy model

This paper describes the solution of a steady state natural convection problem in porous media by the dual reciprocity boundary element method (DRBEM). The boundary element method (BEM) for the coupled set of mass, momentum, and energy equations in two dimensions is structured by the fundamental solution of the Laplace equation. The dual reciprocity method is based on augmented scaled thin plate splines. Numerical examples include convergence studies with different mesh size, uniform and non‐uniform mesh arrangement, and constant and linear boundary field discretizations for differentially heated rectangular cavity problems at filtration with Rayleigh numbers of Ra*=25, 50, and 100 and aspect ratios of A=1/2, 1, and 2. The solution is assessed by comparison with reference results of the fine mesh finite volume method (FVM). Copyright © 2000 John Wiley & Sons, Ltd.     

Scattering of water waves by inclined thin plate submerged in finite-depth water

Summary  The problem of water wave scattering by an inclined thin plate submerged in water of uniform finite depth is investigated here under the assumption of irrotational motion and linear theory. A hypersingular integral equation formulation of the problem is obtained by an appropriate use of Green's integral theorem followed by utilization of the boundary condition on the plate. This hypersingular integral equation involves the discontinuity in the potential function across the plate, which is approximated by a finite series involving Chebyshev polynomials. The coefficients of this finite series are obtained numerically by collocation method. The quantities of physical interest, namely the reflection and transmission coefficients, force and moment acting on the plate per unit width, are then obtained numerically for different values of various parameters, and are depicted graphically against the wavenumber. Effects of finite-depth water, angle of inclination of the plate with the vertical over the deep water and vertical plate results for these quantities are shown. It is observed that the deep-water results effectively hold good if the depth of the mid-point of the submerged plate below the free surface is of the order of one-tenth of the depth of the bottom. Received 30 November 2000; accepted for publication 26 June 2001