Free vibration analysis of an elliptical plate with eccentric elliptical cut-outs

The elaborated collocation multipole method is employed to obtain a semi-analytical solution, involving proper products of angular and radial Mathieu functions, for the free flexural vibrations of a fully clamped thin elastic plate of elliptical planform containing multiple elliptical cutouts of arbitrary size, location, and orientation. The problem boundary conditions are satisfied by uniformly collocating points on the boundaries, and exactly calculating the normal derivative of plate displacement at the collocation points through use of appropriate directional derivative in each coordinate system. The multipole expansion is truncated to yield a coupled algebraic linear system of equations that is then solved for the nontrivial eigensolutions. Extensive numerical simulations present the first three calculated natural frequencies and the associated deformed mode shapes of an elliptical plate with elliptical/circular cutouts, for a wide range of plate/cutout aspect ratios, and cutout location/orientation parameters. The accuracy of solutions is checked through appropriate convergence studies, and the validity of results is established with the aid of a commercial finite element package as well as by comparison with the data available in the existing literature.     

A wideband fast multipole boundary element method for half-space/plane-symmetric acoustic wave problems

This paper presents a novel wideband fast multipole boundary element approach to 3D half-space/planesymmetric acoustic wave problems.The half-space fundamental solution is employed in the boundary integral equations so that the tree structure required in the fast multipole algorithm is constructed for the boundary elements in the real domain only.Moreover,a set of symmetric relations between the multipole expansion coefficients of the real and image domains are derived,and the half-space fundamental solution is modified for the purpose of applying such relations to avoid calculating,translating and saving the multipole/local expansion coefficients of the image domain.The wideband adaptive multilevel fast multipole algorithm associated with the iterative solver GMRES is employed so that the present method is accurate and efficient for both lowand high-frequency acoustic wave problems.As for exterior acoustic problems,the Burton-Miller method is adopted to tackle the fictitious eigenfrequency problem involved in the conventional boundary integral equation method.Details on the implementation of the present method are described,and numerical examples are given to demonstrate its accuracy and efficiency.     

A semi-analytical approach for the non-linear large deflection analysis of laminated rectangular plates under general out-of-plane loading

A semi-analytical approach for the geometrically non-linear analysis of rectangular laminated plates with general inplane and out-of-plane boundary conditions under a general distribution of out-of-plane loads is developed. The analysis is based on the elastic thin plate theory with geometrically non-linear von Kármán strains. The solution of the non-linear partial differential equations is reduced to an iterative sequential solution of non-linear ordinary differential equations using the multi-term extended Kantorovich method. The efficiency, accuracy, and convergence of the proposed method are examined through a comparison with other semi-analytical methods and with finite element analyses. The capabilities of the approach and its applicability to the non-linear large deflection analysis of plate structures are demonstrated through various numerical examples. Emphasis is placed on combinations of lamination, boundary, and loading conditions that cannot be analyzed using alternative semi-analytical methods.     

Analytical investigation for true and spurious eigensolutions of multiply-connected membranes containing elliptical boundaries using the dual BIEM

It is well known that the boundary element method may induce spurious eigenvalues while solving eigenvalue problems. The finding that spurious eigenvalues depend on the geometry of inner boundary and the approach utilized has been revealed analytically and numerically in the literature. However, all the related efforts were focused on eigenproblems involving circular boundaries. On the other hand, the extension to elliptical boundaries seems not straightforward and lacks of attention. Accordingly, this paper performs an analytical investigation of spurious eigenvalues for a confocal elliptical membrane using boundary integral equation methods (BIEM) in conjunction with separable kernels and eigenfunction expansion. To analytically study this eigenproblem, the elliptic coordinates and Mathieu functions are adopted. The fundamental solution is expanded into the separable kernel by using the elliptic coordinates and the boundary densities are expanded by using the eigenfunction expansion. The Jacobian terms may exist in the separable kernel, boundary density and boundary contour integration and they can cancel each other out. Therefore, the orthogonal relations are reserved in the boundary contour integration. In this way, a similar finding about the mechanism of spurious eigenvalues is found and agrees with those corresponding to the annular case. To verify this finding, the boundary element method and the commercial finite-element code ABAQUS are also utilized to provide eigensolutions, respectively, for comparisons. Good agreement is observed from comparisons. Based on the adaptive observer system, the present approach can deal with eigenproblems containing circular and elliptical boundaries at the same time in a semi-analytical manner. By using the BIEM, it is found that spurious eigenvalues are the zeros of the modified Mathieu functions which depend on the inner elliptical boundary and the integral formulation. Finally, several methods including the CHIEF method, the SVD updating technique and the Burton & Miller method are employed to filter out the spurious eigenvalues, respectively. In addition, the efficiency of the CHIEF method is better than those of the SVD updating technique and the Burton & Miller approach, since not only hypersingularity is avoided but also computation effort is saved.     

Finite Difference Method for the Acoustic Radiation of an Elastic Plate Excited by a Turbulent Boundary Layer: A Spectral Domain Solution

A finite difference method is developed to study, on a two-dimensional model, the acoustic pressure radiated when a thin elastic plate, clamped at its boundaries, is excited by a turbulent boundary layer. Consider a homogeneous thin elastic plate clamped at its boundaries and extended to infinity by a plane, perfectly rigid, baffle. This plate closes a rectangular cavity. Both the cavity and the outside domain contain a perfect fluid. The fluid in the cavity is at rest. The fluid in the outside domain moves in the direction parallel to the system plate/baffle with a constant speed. A turbulent boundary layer develops at the interface baffle/plate. The wall pressure fluctuations in this boundary layer generates a vibration of the plate and an acoustic radiation in the two fluid domains. Modeling the wall pressure fluctuations spectrum in a turbulent boundary layer developed over a vibrating surface is a very complex and unresolved task. Ducan and Sirkis [1] proposed a model for the two-way interactions between a membrane and a turbulent flow of fluid. The excitation of the membrane is modeled by a potential flow randomly perturbed. This potential flow is modified by the displacement of the membrane. Howe [2] proposed a model for the turbulent wall pressure fluctuations power spectrum over an elastomeric material. The model presented in this article is based on a hypothesis of one-way interaction between the flow and the structure: the flow generates wall pressure fluctuations which are at the origin of the vibration of the plate, but the vibration of the plate does not modify the characteristics of the flow. A finite difference scheme that incorporates the vibration of the plate and the acoustic pressure inside the fluid cavity has been developed and coupled with a boundary element method that ensures the outside domain coupling. In this paper, we focus on the resolution of the coupled vibration/interior acoustic problem. We compare the results obtained with three numerical methods: (a) a finite difference representation for both the plate displacement and the acoustic pressure inside the cavity; (b) a coupled method involving a finite difference representation for the displacement of the plate and a boundary element method for the interior acoustic pressure; (c) a boundary element method for both the vibration of the plate and the interior acoustic pressure. A comparison of the numerical results obtained with two models of turbulent wall pressure fluctuations spectrums - the Corcos model [3] and the Chase model [4] - is proposed. A difference of 20 dB is found in the vibro-acoustic response of the structure. In [3], this difference is explained by calculating a wavenumber transfer function of the plate. In [6], coupled beam-cavity modes for similar geometry are calculated by the finite difference method. This revised version was published online in July 2006 with corrections to the Cover Date.     

Structural–acoustic sensitivity analysis of radiated sound power using a finite element/ discontinuous fast multipole boundary element scheme

The complete interaction between the structural domain and the acoustic domain needs to be considered in many engineering problems, especially for the acoustic analysis concerning thin structures immersed in water. This study employs the finite element method to model the structural parts and the fast multipole boundary element method to model the exterior acoustic domain. Discontinuous higher‐order boundary elements are developed for the acoustic domain to achieve higher accuracy in the coupling analysis. Structural–acoustic design sensitivity analysis can provide insights into the effects of design variables on radiated acoustic performance and thus is important to the structural–acoustic design and optimization processes. This study is the first to formulate equations for sound power sensitivity on structural surfaces based on an adjoint operator approach and equations for sound power sensitivity on arbitrary closed surfaces around the radiator based on the direct differentiation approach. The design variables include fluid density, structural density, Poisson's ratio, Young's modulus, and structural shape/size. A numerical example is presented to demonstrate the accuracy and validity of the proposed algorithm. Different types of coupled continuous and discontinuous boundary elements with finite elements are used for the numerical solution, and the performances of the different types of finite element/continuous and discontinuous boundary element coupling are presented and compared in detail. Copyright © 2016 John Wiley & Sons, Ltd.     

Diagonal form fast multipole boundary element method for 2D acoustic problems based on Burton-Miller boundary integral equation formulation and its applications

This paper describes formulation and implementation of the fast multipole boundary element method (FMBEM) for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the FMBEM are described in details. They are moment calculation, moment to moment translation, moment to local translation, and local to local translation. A data structure for the quad-tree construction is proposed which can facilitate implementation. An analytical moment expression is derived, which is more accurate, stable, and efficient than direct numerical computation. Numerical examples are presented to demonstrate the accuracy and efficiency of the FMBEM, and radiation of a 2D vibration rail mode is simulated using the FMBEM.     

THE NO-SLIP BOUNDARY CONDITION IN FINITE DIFFERENCE APPROXIMATIONS

A finite difference method for the Navier–Stokes equations in vorticity –streamfunction formulation is proposed to resolve the difficulty of the lack of a vorticity boundary condition at a no-slip boundary. It is particularly suitable for flows in regions with complicated geometries. Convergence with second-order accuracy in vorticity and velocity is established. In numerical experiments the convergence rates agree with theoretical predictions. Test results for the two-dimensional driven cavity problem and for the flow in expansion and contraction channels are given.     

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

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

Moving least squares reconstruction for sharp interface immersed boundary methods

We propose a new approach for reconstructing velocity boundary conditions in sharp-inerface immersed boundary (IB) methods based on the moving least squares (MLS) interpolation method. The MLS is employed to not only reconstruct velocity boundary conditions but also to calculate the pressure and velocity gradients in the vicinity of the immersed body, which are required in fluid structure interaction problems to obtain the force exerted by the fluid on the structure. To extend the method to arbitrarily complex geometries with nonconvex shaped boundaries, the visibility method is combined with the MLS method. The performance of the proposed curvilinear IB MLS (CURVIB-MLS) is demonstrated by systematic grid-refinement studies for two- and three-dimensional tests and compared with the standard CURVIB method employing standard wall-normal interpolation for reconstructing boundary conditions. The test problems are flow in a lid-driven cavity with a sphere, uniform flow over a sphere, flow on a NACA0018 airfoil at incidence, and vortex-induced vibration of an elastically-mounted cylinder. We show that the CURVIB-MLS formulation yields a method that is easier to implement in complex geometries and exhibits higher accuracy and rate of convergence relative to the standard CURVIB method. The MLS approach is also shown to dramatically improve the accuracy of calculating the pressure and viscous forces imparted by the flow on the body and improve the overall accuracy of FSI simulations. Finally, the CURVIB-MLS approach is able to qualitatively capture on relatively coarse grids important features of complex separated flows that the standard CURVIB method is able to capture only on finer grids.     

Variety of acoustic streaming in 2D resonant channels

Acoustic streaming in 2D resonant channels with uniform or non-uniform cross-sections is studied within this work. An inertial force as well as a vibrating boundary are assumed for driving the acoustic field. The method of successive approximations is employed to derive linear equations for calculation of primary acoustic and time-averaged secondary fields including the radiation pressure and the mass transport velocity. The model equations have a standard form which allows their numerical integration using a universal solver; in this case, COMSOL Multiphysics was employed. As this software is based on the finite element method, it is simple and straightforward to perform the calculations with moderate computational costs even for complex geometries, which makes the proposed approach an operative tool for study of acoustic streaming. The numerical results are validated for the case of a rectangular channel by comparison with previously published analytical results; an excellent agreement is found. The numerical results show that the acoustic streaming can be quite complex even in rectangular channels and its structure depends on the manner of driving. Examples of acoustic streaming in wedged and elliptical channels are given to demonstrate a strong dependence of the acoustic streaming structure on the resonator shape.     

双周期加筋微穿孔板的振动响应及声透射的研究

本文主要研究了水下无穷大双周期加筋微穿孔薄板,在平面声波斜入射下的振动响应和声透射,并提出了一种半解析半数值的计算方法。利用微穿孔板的声阻抗以及薄板表面的振速边界条件,建立了加筋穿孔薄板的振动方程,并根据傅立叶变换及空间波数法将振动位移表达为波数分量的迭加形式。采用数值计算的方法对波数分量进行求解并通过傅里叶逆变换,最终得到了双周期加筋穿孔薄板的振动响应及透射系数。通过与Takahashi穿孔板声压结果的对比,证明了本方法的正确性。在算例中,分析了加强筋及穿孔率对薄板结构的振动和声透射的影响。     

快速多极虚边界元法对含圆孔薄板有效弹性模量的模拟分析

针对虚边界元法,引入快速多极展开和广义极小残值法(GMRES)的思想,以形成快速多极虚边界元法的求解思想,并将此方法用于含圆孔薄板有效弹性模量的模拟分析.由于本文方法采用了"源点"多极展开和"场点"局部展开的组合处理方案,从而使得原问题方程组求解的计算耗时量和储存量降至与所求问题的计算自由度数成线性比例.本文工作的研究目的在于:提高虚边界元法在普通台式机上的运算能力和拓宽虚边界元法对大规模复杂问题的求解(或数值模拟).文中给出了均布圆孔的正方形薄板和之字形分布圆孔薄板二个算例,以验证该方法的可行性,计算精度和计算效率.     

A reined global-local approach for evaluation of singular stress ield based on scaled boundary inite element method

A reined global-local approach based on the scaled boundary inite element method(SBFEM) is proposed to improve the accuracy of predicted singular stress ield. The proposed approach is carried out in conjunction with two steps. First, the entire structure is analyzed by employing an arbitrary numerical method. Then, the interested region, which contains stress singularity, is re-solved using the SBFEM by placing the scaling center right at the singular stress point with the boundary conditions evaluated from the irst step imposed along the whole boundary including the side-faces. Beneiting from the semi-analytical nature of the SBFEM, the singular stress ield can be predicted accurately without highly reined meshes. It provides the FEM or other numerical methods with a rather simple and convenient way to improve the accuracy of stress analysis. Numerical examples validate the effectiveness of the proposed approach in dealing with various kinds of problems.     

广义有限差分法在含阻抗边界空腔声学分析中的应用

声学分析在噪声控制、室内隔音等工程计算中有着重要的作用. 由于现实生活中的声学模型往往伴随着吸声材料, 因此分析含阻抗边界条件的声学问题显得十分必要. 广义有限差分法是一种新型区域型无网格数值离散方法, 该方法基于多元函数泰勒级数展开式和加权最小二乘拟合, 将控制方程中未知参量的各阶偏导数表示为相邻节点函数值的线性组合. 本文首次将广义有限差分法应用于含阻抗边界条件空腔声学问题的分析中, 建立了空腔声场问题的广义有限差分法数值离散格式. 与传统算法相比, 所建立的数值模型具有无需网格剖分和数值积分、计算精度高、适用于大规模声学分析等优点. 通过具有解析解的经典算例, 研究了总节点数目和局部支撑点数目对数值结果的影响, 得到了最大计算频率与节点间距之间关系的经验公式. 此外, 将广义有限差分法应用于无解析解的二维和三维复杂声学模型, 并与COMSOL Multiphysics软件所得的有限元结果进行了比较分析. 数值实验表明, 该算法是一种高效、精确、稳定、收敛的数值模拟方法, 在含阻抗边界空腔声学分析中具有广阔的应用前景.      

边界节点法的计算稳定性研究

边界节点法利用满足控制方程的非奇异通解作为基函数,半解析边界数值离散偏微分方程,具有精度高、收敛快、易编程等优点,是一种纯无网格配点方法.但是在求解具体问题时,随着节点数的增加,边界节点法经常得到严重病态的插值矩阵.本文利用有效条件数评价边界节点法求解Helmholtz问题线性方程组的计算稳定性;然后利用三种正则化方法处理其病态的线性方程组,并与高斯消元法比较计算精度和收敛性.通过数值实验,本文研究了有效条件数、误差和正则化方法之间的关系.     

Nonreflecting boundary conditions based on nonlinear multidimensional characteristics

Nonlinear characteristic boundary conditions based on nonlinear multidimensional characteristics are proposed for 2‐ and 3‐D compressible Navier–Stokes equations with/without scalar transport equations. This approach is consistent with the flow physics and transport properties. Based on the theory of characteristics, which is a rigorous mathematical technique, multidimensional flows can be decomposed into acoustic, entropy, and vorticity waves. Nonreflecting boundary conditions are derived by setting corresponding characteristic variables of incoming waves to zero and by partially damping the source terms of the incoming acoustic waves. In order to obtain the resulting optimal damping coefficient, analysis is performed for problems of pure acoustic plane wave propagation and arbitrary flows. The proposed boundary conditions are tested on two benchmark problems: cylindrical acoustic wave propagation and the wake flow behind a cylinder with strong periodic vortex convected out of the computational domain. This new approach substantially minimizes the spurious wave reflections of pressure, density, temperature, and velocity as well as vorticity from the artificial boundaries, where strong multidimensional flow effects exist. The numerical simulations yield accurate results, confirm the optimal damping coefficient obtained from analysis, and verify that the method substantially improves the 1‐D characteristics‐based nonreflecting boundary conditions for complex multidimensional flows. Copyright © 2009 John Wiley & Sons, Ltd.     

内部压力作用下矩形板中源于椭圆孔的分支裂纹应力强度因子的一种数值分析

应用一种边界元方法来研究内部压力作用下矩形板中源于椭圆孔的分支裂纹。该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者最近提出的裂尖位移不连续单元构成。在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界。本数值结果进一步证实这种数值方法对计算有限大板中复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示裂纹体几何对应力强度因子的影响。     

Semi-analytical solutions to one-dimensional consolidation for unsaturated soils with semi-permeable drainage boundary

The semi-analytical solutions to Fredlund and Hasan's one-dimensional(1 D)consolidation for unsaturated soils with a semi-permeable drainage boundary are presented. Two variables are introduced to transform the two coupled governing equations of pore-water and pore-air pressures into an equivalent set of partial differential equations(PDFs), which are easily solved by the Laplace transform method. Then, the pore-water pressure, pore-air pressure, and soil settlement are obtained in the Laplace domain. The Crump method is adopted to perform the inverse Laplace transform in order to obtain the semi-analytical solutions in the time domain. It is shown that the proposed solutions are more applicable to various types of boundary conditions and agree well with the existing solutions from the literature. Several numerical examples are provided to investigate the consolidation behavior of an unsaturated single-layer soil with single, double, mixed, and semi-permeable drainage boundaries. The changes in the pore-air and pore-water pressures and the soil settlement with the time factor at different values of the semi-permeable drainage boundary parameters are illustrated. In addition, parametric studies are conducted on the pore-air and pore-water pressures at different ratios(the air permeability coefficient to the water permeability coefficient) and depths.     

Contact Transform for the Biharmonic Equation Applicable to Plane Stress Elastoplastic Elliptical Hole Problems

An analytical technique is developed that reduces the unknown elastic-plastic boundary of a linear elastic-perfectly plastic material containing an elliptical hole under tensile plane stress loading conditions into an equivalent mathematical problem with known boundaries. This mathematical transformation may facilitate this problem’s solution by either analytical or numerical means. Yield is assumed to occur in this analysis under the Tresca yield criterion. An example elastic-plastic problem illustrating this method is drawn from existing literature in the form of a perturbation solution for an elliptical hole derived by a series expansion about a circular boundary.