A novel hybrid FS‐FEM/SEA for the analysis of vibro‐acoustic problems

Predicting the frequency response of a complex vibro‐acoustic system becomes extremely difficult in the mid‐frequency regime. In this work, a novel hybrid face‐based smoothed finite element method/statistical energy analysis (FS‐FEM/SEA) method is proposed, aiming to further improve the accuracy of ‘mid‐frequency’ predictions. According to this approach, the whole vibro‐acoustic system is divided into a combination of a plate subsystem with statistical behaviour and an acoustic cavity subsystem with deterministic behaviour. The plate subsystem is treated using the recently developed FS‐FEM, and the cavity subsystem is dealt with using the SEA. These two different types of subsystems can be coupled and interacted through the so‐called diffuse field reciprocity relation. The ensemble average response of the system is calculated, and the uncertainty is confined and treated in the SEA subsystems. The use of FS‐FEM ‘softens’ the well‐known ‘overly stiff’ behaviour in the standard FEM and reduces the inherent numerical dispersion error. The proposed FS‐FEM/SEA approach is verified and its features are examined by various numerical examples. Copyright © 2015 John Wiley & Sons, Ltd.     

Dispersion analysis of the gradient weighted finite element method for acoustic problems in one,two, and three dimensions

This paper reports a detailed analysis on the numerical dispersion error in solving one-, two-, and three-dimensional acoustic problems governed by the Helmholtz equation using the gradient weighted finite element method (GW-FEM) in comparison with the standard FEM and the modified methods presented in the literatures. The discretized system equations derived based on the gradient weighted operation corresponding to the considered method are first briefed. The discrete dispersion relationships relating the exact and numerical wave numbers defined in different dimensions are then formulated, which will be further used to investigate the dispersion effect mainly caused by the approximation of field variables. The influence of nondimensional wave number and wave propagation angle on the dispersion error is detailedly studied. Comparisons are made with the classical FEM and high-performance algorithms. Results of both theoretical and numerical experiments show that the present method can effectively reduce the pollution effect in computational acoustics owning to its crucial effectiveness in handing the dispersion error in the discrete numerical model.     

Selective smoothed finite element methods for extremely large deformation of anisotropic incompressible bio‐tissues

Two dynamic selective smoothed FEM (selective S‐FEM) are proposed for analysis of extremely large deformation of anisotropic incompressible bio‐tissues using the simplest four‐node tetrahedron elements. In the present two Selective S‐FEMs, the method that consists of the face‐based smoothed FEM (FS‐FEM) used for the deviatoric part of deformation and the node‐based smoothed FEM (NS‐FEM) used for the volumetric part is called FS/NS‐FEM; another method that replaces the deviatoric part of deformation in the first one by the edge‐based smoothed FEM (3D‐ES‐FEM) is call 3D‐ES/NS‐FEM. Both selective S‐FEMs can achieve outstanding accuracy, and stability of volumetric locking free. This is because the NS‐FEM offers an ‘overly‐soft’ feature (in contrast to the standard FEM ‘overly‐stiff’ model), which can be used to effectively mitigate the volumetric locking, and on the other hand, the 3D‐ES‐FEM and FS‐FEM produce close to exact stiffness for the discretized model leading to accurate solution. Numerical examples are presented to examine the performance of the selective S‐FEM methods, including soft bio‐tissues that may be isotropic, transversely isotropic, and anisotropic arterial layered materials. The present methods are found having good accuracy and performance. The examples also demonstrate that the proposed methods are very robust and possess remarkable capabilities of handling element distortion, which is very useful for simulating soft materials including bio‐tissues. Copyright © 2014 John Wiley & Sons, Ltd.     

Analysis of high‐order finite elements for convected wave propagation

In this paper, we examine the performance of high‐order finite element methods (FEM) for aeroacoustic propagation, based on the convected Helmholtz equation. A methodology is presented to measure the dispersion and amplitude errors of the p‐FEM, including non‐interpolating shape functions, such as ‘bubble’ shape functions. A series of simple test cases are also presented to support the results of the dispersion analysis. The main conclusion is that the properties of p‐FEM that make its strength for standard acoustics (e.g., exponential p‐convergence, low dispersion error) remain present for flow acoustics as well. However, the flow has a noticeable effect on the accuracy of the numerical solution, even when the change in wavelength due to the mean flow is accounted for, and an approximation of the dispersion error is proposed to describe the influence of the mean flow. Also discussed is the so‐called aliasing effect, which can reduce the accuracy of the solution in the case of downstream propagation. This can be avoided by an appropriate choice of mesh resolution. Copyright © 2013 John Wiley & Sons, Ltd.     

A quasi‐static crack growth simulation based on the singular ES‐FEM

In the edge‐based smoothed finite element method (ES‐FEM), one needs only the assumed displacement values (not the derivatives) on the boundary of the edge‐based smoothing domains to compute the stiffness matrix of the system. Adopting this important feature, a five‐node crack‐tip element is employed in this paper to produce a proper stress singularity near the crack tip based on a basic mesh of linear triangular elements that can be generated automatically for problems with complicated geometries. The singular ES‐FEM is then formulated and used to simulate the crack propagation in various settings, using a largely coarse mesh with a few layers of fine mesh near the crack tip. The results demonstrate that the singular ES‐FEM is much more accurate than X‐FEM and the existing FEM. Moreover, the excellent agreement between numerical results and the reference observations shows that the singular ES‐FEM offers an efficient and high‐quality solution for crack propagation problems. Copyright © 2011 John Wiley & Sons, Ltd.     

Efficient implementation of high‐order finite elements for Helmholtz problems

Computational modeling remains key to the acoustic design of various applications, but it is constrained by the cost of solving large Helmholtz problems at high frequencies. This paper presents an efficient implementation of the high‐order finite element method (FEM) for tackling large‐scale engineering problems arising in acoustics. A key feature of the proposed method is the ability to select automatically the order of interpolation in each element so as to obtain a target accuracy while minimizing the cost. This is achieved using a simple local a priori error indicator. For simulations involving several frequencies, the use of hierarchical shape functions leads to an efficient strategy to accelerate the assembly of the finite element model. The intrinsic performance of the high‐order FEM for 3D Helmholtz problem is assessed, and an error indicator is devised to select the polynomial order in each element. A realistic 3D application is presented in detail to demonstrate the reduction in computational costs and the robustness of the a priori error indicator. For this test case, the proposed method accelerates the simulation by an order of magnitude and requires less than a quarter of the memory needed by the standard FEM. © 2016 The Authors. International Journal for Numerical Methods in Engineering published by John Wiley & Sons Ltd.     

Dispersion free analysis of acoustic problems using the alpha finite element method

The classical finite element method (FEM) fails to provide accurate results to the Helmholtz equation with large wave numbers due to the well-known “pollution error” caused by the numerical dispersion, i.e. the numerical wave number is always smaller than the exact one. This dispersion error is essentially rooted at the “overly-stiff” feature of the FEM model. In this paper, an alpha finite element method (α-FEM) is then formulated for the acoustic problems by combining the “smaller wave number” model of FEM and the “larger wave number” model of NS-FEM through a scaling factor ${a\in [0,1]}The classical finite element method (FEM) fails to provide accurate results to the Helmholtz equation with large wave numbers due to the well-known “pollution error” caused by the numerical dispersion, i.e. the numerical wave number is always smaller than the exact one. This dispersion error is essentially rooted at the “overly-stiff” feature of the FEM model. In this paper, an alpha finite element method (α-FEM) is then formulated for the acoustic problems by combining the “smaller wave number” model of FEM and the “larger wave number” model of NS-FEM through a scaling factor a ? [0,1]{a\in [0,1]} . The motivation for this combined approach is essentially from the features of “overly-stiff” FEM model and “overly-soft” NS-FEM model, and accurate solutions can be obtained by tuning the α-FEM model. A technique is proposed to determine a particular alpha with which the α-FEM model can possess a very “close-to-exact” stiffness, which can effectively reduce the dispersion error leading to dispersion free solutions for acoustic problems. Theoretical and numerical studies shall demonstrate the excellent properties of the present α-FEM.     

One‐dimensional dispersion analysis for the element‐free Galerkin method for the Helmholtz equation

The standard finite element method (FEM) is unreliable to compute approximate solutions of the Helmholtz equation for high wave numbers due to the dispersion, unless highly refined meshes are used, leading to unacceptable resolution times. The paper presents an application of the element‐free Galerkin method (EFG) and focuses on the dispersion analysis in one dimension. It shows that, if the basis contains the solution of the homogenized Helmholtz equation, it is possible to eliminate the dispersion in a very natural way while it is not the case for the finite element methods. For the general case, it also shows that it is possible to choose the parameters of the method in order to minimize the dispersion. Finally, theoretical developments are validated by numerical experiments showing that, for the same distribution of nodes, the element‐free Galerkin method solution is much more accurate than the finite element one. Copyright © 2000 John Wiley & Sons, Ltd.     

Stability and accuracy of finite element methods for flow acoustics. I: general theory and application to one‐dimensional propagation

The dispersion properties of finite element models for aeroacoustic propagation based on the convected scalar Helmholtz equation and on the Galbrun equation are examined. The current study focusses on the effect of the mean flow on the dispersion and amplitude errors present in the discrete numerical solutions. A general two‐dimensional dispersion analysis is presented for the discrete problem on a regular unbounded mesh, and results are presented for the particular case of one‐dimensional acoustic propagation in which the wave direction is aligned with the mean flow. The magnitude and sign of the mean flow is shown to have a significant effect on the accuracy of the numerical schemes. Quadratic Helmholtz elements in particular are shown to be much less effective for downstream—as opposed to upstream—propagation, even when the effect of wave shortening or elongation due to the mean flow is taken into account. These trends are also observed in solutions obtained for simple test problems on finite meshes. A similar analysis of two‐dimensional propagation is presented in an accompanying article. Copyright © 2005 John Wiley & Sons, Ltd.     

F‐bar aided edge‐based smoothed finite element method using tetrahedral elements for finite deformation analysisof nearly incompressible solids

A new smoothed finite element method (S‐FEM) with tetrahedral elements for finite strain analysis of nearly incompressible solids is proposed. The proposed method is basically a combination of the F‐bar method and edge‐based S‐FEM with tetrahedral elements (ES‐FEM‐T4) and is named ‘F‐barES‐FEM‐T4’. F‐barES‐FEM‐T4 inherits the accuracy and shear locking‐free property of ES‐FEM‐T4. At the same time, it also inherits the volumetric locking‐free property of the F‐bar method. The isovolumetric part of the deformation gradient ( F ^iso) is derived from the F of ES‐FEM‐T4, whereas the volumetric part ( F ^vol) is derived from the cyclic smoothing of J(=det( F )) between elements and nodes. Some demonstration analyses confirm that F‐barES‐FEM‐T4 with a sufficient number of cyclic smoothings suppresses the pressure oscillation in nearly incompressible materials successfully with no increase in DOF. Moreover, they reveal that our method is capable of relaxing the corner locking issue arising at the corner in the cylinder barreling analysis. Copyright © 2016 John Wiley & Sons, Ltd.     

光滑有限元的声学研究：时域和频域分析

摘　要：在使用有限元进行声场的数值模拟中,存在着两个主要误差,一个是数值方法中常规的插值误差,另外一个是计算声学中所特有的耗散误差(dispersion error),后者则是影响声学模拟仿真置信度的最重要因素。产生耗散误差的本质原因是由于有限元的数值模型刚度“偏硬”造成的。为了控制耗散误差,最重要的是使数值模型更好的反映真实模型。本文采用了一种基于边光滑的有限元方法(ES-FEM)来对声场的时域和频域进行数值模拟研究。该方法只采用对复杂问题域适应性很强的三角形网格,通过引进基于边的广义梯度光滑技术,能够使得有限元系统得到适当的“软化”。关于时域和频域的算例表明了在使用同样网格的情况下,本方法在声学模拟中的精度都要比有限元模型的高。     

Upper and lower bounds for natural frequencies: A property of the smoothed finite element methods

Node‐based smoothed finite element method (NS‐FEM) using triangular type of elements has been found capable to produce upper bound solutions (to the exact solutions) for force driving static solid mechanics problems due to its monotonic ‘soft’ behavior. This paper aims to formulate an NS‐FEM for lower bounds of the natural frequencies for free vibration problems. To make the NS‐FEM temporally stable, an α‐FEM is devised by combining the compatible and smoothed strain fields in a partition of unity fashion controlled by α∈[0, 1], so that both the properties of stiff FEM and the monotonically soft NS‐FEM models can be properly combined for a desired purpose. For temporally stabilizing NS‐FEM, α is chosen small so that it acts like a ‘regularization parameter’ making the NS‐FEM stable, but still with sufficient softness ensuring lower bounds for natural frequency solution. Our numerical studies demonstrate that (1) using a proper α, the spurious non‐zero energy modes can be removed and the NS‐FEM becomes temporally stable; (2) the stabilized NS‐FEM becomes a general approach for solids to obtain lower bounds to the exact natural frequencies over the whole spectrum; (3) α‐FEM can even be tuned for obtaining nearly exact natural frequencies. Copyright © 2010 John Wiley & Sons, Ltd.     

Dispersion analysis and error estimation of Galerkin finite element methods for the Helmholtz equation

When applying numerical methods for the computation of stationary waves from the Helmholtz equation, one obtains ‘numerical waves’ that are dispersive also in non-dispersive media. The numerical wave displays a phase velocity that depends on the parameter k of the Helmholtz equation. In dispersion analysis, the phase difference between the exact and the numerical solutions is investigated. In this paper, the authors' recent result on the phase difference for one-dimensional problems is numerically evaluated and discussed in the context of other work directed to this topic. It is then shown that previous error estimates in H¹-norm are of nondispersive character but hold for medium or high wavenumber on extremely refined mesh only. On the other hand, recently proven error estimates for constant resolution contain a pollution term. With certain assumptions on the exact solution, this term is of the order of the phase difference. Thus a link is established between the results of dispersion analysis and the results of numerical analysis. Throughout the paper, the presentation and discussion of theoretical results is accompanied by numerical evaluation of several model problems. Special attention is given to the performance of the Galerkin method with a higher order of polynomial approximation p(h-p-version).     

Element‐wise a posteriori estimates based on hierarchical bases for non‐linear parabolic problems

We show that the issue of a posteriori estimate the errors in the numerical simulation of non‐linear parabolic equations can be reduced to a posteriori estimate the errors in the approximation of an elliptic problem with the right‐hand side depending on known data of the problem and the computed numerical solution. A procedure to obtain local error estimates for the p version of the finite element method by solving small discrete elliptic problems with right‐hand side the residual of the p‐FEM solution is introduced. The boundary conditions are inherited by those of the space of hierarchical bases to which the error estimator belongs. We prove that the error in the numerical solution can be reduced by adding the estimators that behave as a locally defined correction to the computed approximation. When the error being estimated is that of a elliptic problem constant free local lower bounds are obtained. The local error estimation procedure is applied to non‐linear parabolic differential equations in several space dimensions. Some numerical experiments for both the elliptic and the non‐linear parabolic cases are provided. Copyright © 2005 John Wiley & Sons, Ltd.     

声学数值计算的分区光滑径向点插值无网格法

针对标准的有限元法分析声学问题时由于数值色散导致高波数计算结果不可靠问题,将分区光滑径向点插值法(cell-based smoothed radial point interpolation method, CS-RPIM)应用到二维声学分析中,推导了分区光滑径向点插值法分析二维声学问题的原理公式。该方法将问题域划分为三角形背景单元,每个单元进一步分成若干个光滑域,对每个光滑域进行声压梯度光滑处理,运用光滑Galerkin弱形式构造系统方程,并按有限元中方法施加必要的边界条件。CS-RPIM提供了合适的模型硬度,能有效降低色散效应,提高计算精度。对管道和二维轿车声学问题的数值分析结果表明,与标准有限元法相比,CS-RPIM具有更高的精度和准确度,在高波数计算时这种优势特别明显。     

Smoothed finite element method for topology optimization involving incompressible materials

It is well known that the finite element method (FEM) suffers severely from the volumetric locking problem for incompressible materials in topology optimization owing to its numerical ‘overly stiff’ property. In this article, two typical smoothed FEMs with a certain softened effect, namely the node-based smoothed finite element method (NS-FEM) and the cell-based smoothed finite element method, are formulated to model the compressible and incompressible materials for topology optimization. Numerical examples have demonstrated that the NS-FEM with an ‘overly soft’ property is fairly effective in tackling the volumetric locking problem in topology optimization when both compressible and incompressible materials are involved.     

Stability and accuracy of finite element methods for flow acoustics. II: Two‐dimensional effects

This is the second of two articles that focus on the dispersion properties of finite element models for acoustic propagation on mean flows. We consider finite element methods based on linear potential theory in which the acoustic disturbance is modelled by the convected Helmholtz equation, and also those based on a mixed Galbrun formulation in which acoustic pressure and Lagrangian displacement are used as discrete variables. The current paper focuses on the effects of numerical anisotropy which are associated with the orientation of the propagating wave to the mean flow and to the grid axes. Conditions which produce aliasing error in the Helmholtz formulation are of particular interest. The 9‐noded Lagrangian element is shown to be superior to the more commonly used 8‐noded serendipity element. In the case of the Galbrun elements, the current analysis indicates that isotropic meshes generally reduce numerical error of triangular elements and that higher order mixed quadrilaterals are generally less effective than an equivalent mesh of lower order triangles. Copyright © 2005 John Wiley & Sons, Ltd.     

An edge‐based smoothed finite element method for 3D analysis of solid mechanics problems

The edge‐based smoothed finite element method (ES‐FEM) was proposed recently in Liu, Nguyen‐Thoi, and Lam to improve the accuracy of the FEM for 2D problems. This method belongs to the wider family of the smoothed FEM for which smoothing cells are defined to perform the numerical integration over the domain. Later, the face‐based smoothed FEM (FS‐FEM) was proposed to generalize the ES‐FEM to 3D problems. According to this method, the smoothing cells are centered along the faces of the tetrahedrons of the mesh. In the present paper, an alternative method for the extension of the ES‐FEM to 3D is investigated. This method is based on an underlying mesh composed of tetrahedrons, and the approximation of the field variables is associated with the tetrahedral elements; however, in contrast to the FS‐FEM, the smoothing cells of the proposed ES‐FEM are centered along the edges of the tetrahedrons of the mesh. From selected numerical benchmark problems, it is observed that the ES‐FEM is characterized by a higher accuracy and improved computational efficiency as compared with linear tetrahedral elements and to the FS‐FEM for a given number of degrees of freedom. Copyright © 2013 John Wiley & Sons, Ltd.     

Overview of the discontinuous enrichment method,the ultra‐weak variational formulation,and the partition of unity method for acoustic scattering in the medium frequency regime and performance comparisons

The efficient finite element discretization of the Helmholtz equation becomes challenging in the medium frequency regime because of numerical dispersion, or what is often referred to in the literature as the pollution effect. A number of FEMs with plane wave basis functions have been proposed to alleviate this effect, and improve on the unsatisfactory preasymptotic convergence of the polynomial FEM. These include the partition of unity method, the ultra‐weak variational formulation, and the discontinuous enrichment method. A previous comparative study of the performance of such methods focused on the first two aforementioned methods only. By contrast, this paper provides an overview of all three methods and compares several aspects of their performance for an acoustic scattering benchmark problem in the medium frequency regime. It is found that the discontinuous enrichment method outperforms both the partition of unity method and the ultra‐weak variational formulation by a significant margin. Copyright © 2011 John Wiley & Sons, Ltd.     

Dispersion and pollution of the FEM solution for the Helmholtz equation in one,two and three dimensions

For high wave numbers, the Helmholtz equation suffers the so‐called ‘pollution effect’. This effect is directly related to the dispersion. A method to measure the dispersion on any numerical method related to the classical Galerkin FEM is presented. This method does not require to compute the numerical solution of the problem and is extremely fast. Numerical results on the classical Galerkin FEM (p‐method) is compared to modified methods presented in the literature. A study of the influence of the topology triangles is also carried out. The efficiency of the different methods is compared. The numerical results in two of the mesh and for square elements show that the high order elements control the dispersion well. The most effective modified method is the QSFEM [1,2] but it is also very complicated in the general setting. The residual‐free bubble [3,4] is effective in one dimension but not in higher dimensions. The least‐square method [1,5] approach lowers the dispersion but relatively little. The results for triangular meshes show that the best topology is the ‘criss‐cross’ pattern. Copyright © 1999 John Wiley & Sons, Ltd.