Calculating the forced response of cylinders and cylindrical shells using the wave and finite element method

The dynamic response of circular cylinders can be obtained analytically in very few (and simple) cases. For complicated (thick or anisotropic) circular cylinders, researchers often resort to the finite element (FE) method. This can lead to large models, especially at higher frequencies, which translates into high computational costs and memory requirements. In this paper, the response of axially homogenous circular cylinders (that can be arbitrarily complex through the thickness) is obtained using the wave and finite element (WFE) method. Here, the homogeneity of the cylinder around the circumference and along the axis are exploited to post-process the FE model of a small rectangular segment of the cylinder using periodic structure theory and obtain the wave characteristics of the cylinder. The full power of FE methods can be utilised to obtain the FE model of the small segment. Then, the forced response of the cylinder is posed as an inverse Fourier transform. However, since there are an integer number of wavelengths around the circumference of a closed circular cylinder, one of the integrals in the inverse Fourier transform becomes a simple summation, whereas the other can be resolved analytically using contour integration and the residue theorem. The result is a computationally efficient technique for obtaining the response to time harmonic, arbitrarily distributed loads of axially homogenous, circular cylinders with arbitrary complexity across the thickness.     

On the forced response of waveguides using the wave and finite element method

The forced response of waveguides subjected to time harmonic loading is treated. The approach starts with the wave and finite element (WFE) method where a segment of the waveguide is modeled using traditional finite element methods. The mass and stiffness matrices of the segment are used to formulate an eigenvalue problem whose solution yields the wave properties of the waveguide. The WFE formulation is used to obtain the response of the waveguide to a convected harmonic pressure (CHP). Since the Fourier transform of the response to a general excitation is a linear combination of the responses to CHPs, the response to a general excitation can be obtained via an inverse Fourier transform process. This is evaluated analytically using contour integration and the residue theorem. Hence, the approach presented herein enables the response of a waveguide to general loading to be found by: (a) modeling a segment of the waveguide using finite element methods and post-processing it to obtain the wave characteristics, (b) using Fourier transform and contour integration to obtain the wave amplitudes and (c) using the wave amplitudes to find the response at any point in the waveguide. Numerical examples are presented.     

Waves and energy in random elastic guided media through the stochastic wave finite element method

Energy propagation in random viscoelastic media is considered in this Letter. The forced response of uncertain waveguide subject to time harmonic loading is treated. This energy model is based on a spectral approach called the “Stochastic Wave Finite Element” (SWFE) method which is detailed in this Letter. Assuming that the random properties are spatially homogeneous in the media, the SWFE is a hybridization of the deterministic wave finite element and a parametric probabilistic approach. The proposed model is applicable in a wide frequency band with reduced time consumption. Numerical examples show the effectiveness of the proposed approach to predict the statistics of kinematic and quadratic variables of guided wave propagation. The results are compared to Monte Carlo simulations.     

Simulation of Lamb wave reflections at plate edges using the semi-analytical finite element method

In typical Lamb wave simulation practices, effects of plate edge reflections are often not considered in order to simplify the wave signal interpretations. Methods that are based on infinite plates such as the semi-analytical finite element method is effective in simulating Lamb waves as it excludes the effect of plate edges. However, the inclusion of plate edges in a finite plate could render this method inapplicable, especially for transient response simulations. Here, by applying the ratio of Lamb mode reflections at plate edges, and representing the reflection at plate edges using infinite plate solutions, the semi-analytical finite element method can be applied for transient response simulation, even when the plate is no longer infinite.     

PVDF水听器加速度响应特性的有限元分析

本文用有限元法分析了压电高聚物PVDF圆柱形水听器的材料参数与几何尺寸对其加速度响应的影响。针对理论分析,进行了水听器加速度响应实验,实验结果与理论分析基本相符。     

Viscoelastic analysis of multi-body systems using the finite element method

A study of the effect of viscoelastic material damping on the dynamic response of multibody systems, consisting of interconnected rigid, elastic and viscoelastic components, is presented. The motion of each elastic or viscoelastic body is identified by using three sets of modes: rigid body, reference and normal modes. Rigid body modes describe translation and large angular rotation of a body reference. Reference modes are the result of imposing the body-axis conditions. Normal modes define the deformation of the body relative to the body reference. Constraints between different components are formulated by using a set of non-linear algebraic equations that can be introduced to the dynamic formulation by using a Lagrange multiplier technique or can be utilized to eliminate dependent co-ordinates by partitioning the constraint Jacobian matrix. In developing the system equations of motion of the viscoelastic component, an assumption of a linear viscoelastic model is made. A Kelvin-Voigt model is employed, wherein the stress is assumed to be proportional to the strain and its time derivative. The formulation yields a constant damping matrix and the damping forces depend only on the local deformation; thus, no additional coupling between the reference and elastic co-ordinates appears in the formulation when considering the viscoelastic effects. It is demonstrated, by a numerical example, that the viscoelastic material damping can have a significant effect on the dynamic response of multibody systems.     

Characterization of silicon nanophotonic devices using the finite element method

The Poynting vector and the full-vectorial H and E-field profiles are considered for use in nanoscale silicon waveguides in this article. This paper reveals that the mode profile of a circular silicon nanowire is not circular and also has a strong axial field component. From the analysis, the characteristics of single mode operation and the vector field profiles of both circular and planar silicon nanowires are presented. The modal birefringence of rectangular silicon nanowires and power density in low-index region of a slot-type waveguide and designs of a compact polarization rotator are also presented in this work.     

Hybrid finite volume/ finite element method for radiative heat transfer in graded index media

The rays propagate along curved path determined by the Fermat principle in the graded index medium. The radiative transfer equation in graded index medium (GRTE) contains two specific redistribution terms (with partial derivatives to the angular coordinates) accounting for the effect of the curved ray path. In this paper, the hybrid finite volume with finite element method (hybrid FVM/FEM) (P.J. Coelho, J. Quant. Spectrosc. Radiat. Transf., vol. 93, pp. 89–101, 2005) is extended to solve the radiative heat transfer in two-dimensional absorbing-emitting-scattering graded index media, in which the spatial discretization is carried out using a FVM, while the angular discretization is by a FEM. The FEM angular discretization is demonstrated to be preferable in dealing with the redistribution terms in the GRTE. Two stiff matrix assembly schemes of the angular FEM discretization, namely, the traditional assembly approach and a new spherical assembly approach (assembly on the unit sphere of the solid angular space), are discussed. The spherical assembly scheme is demonstrated to give better results than the traditional assembly approach. The predicted heat flux distributions and temperature distributions in radiative equilibrium are determined by the proposed method and compared with the results available in other references. The proposed hybrid FVM/FEM method can predict the radiative heat transfer in absorbing-emitting-scattering graded index medium with good accuracy.     

应用声有限元素法模拟坐席吸声低谷效应

本文给出了应用声有限元素法模拟坐席吸声低谷效应的计算模型。该模型首先在时间域计算房间脉冲响应,再通过傅立叶变换得到声压频谱特性。将有限元分析结果与已有的模型实验及大厅实测数据相比较,证明该方法是研究坐席吸声低谷效应的有效手段。     

反射声波成像测井的有限元模拟

应用有限元法对新兴的反射声波成像测井进行了数值模拟研究，计算了井旁附近存在不同倾角声阻抗不连续界面对声波的反射，使用的井内激励源为幅度加权相控线阵声波辐射器。从计算出的声全波波形上可以清楚地看到沿井壁传播的折射波和来自井旁界面的反射纵波。应用偏移叠加等数据处理技术对数值模拟波形进行了处理，从偏移剖面上反演出的反射界面的尺寸、位置和倾角等几何特性与数值模拟输入参数基本一致。     

Boundary element method for calculation of elastic wave transmission in two-dimensional phononic crystals

A boundary element method (BEM) is presented to compute the transmission spectra of two-dimensional (2-D) phononic crystals of a square lattice which are finite along the x-direction and infinite along the y-direction. The cross sections of the scatterers may be circular or square. For a periodic cell, the boundary integral equations of the matrix and the scatterers are formulated. Substituting the periodic boundary conditions and the interface continuity conditions, a linear equation set is formed, from which the elastic wave transmission can be obtained. From the transmission spectra, the band gaps can be identified, which are compared with the band structures of the corresponding infinite systems. It is shown that generally the transmission spectra completely correspond to the band structures. In addition, the accuracy and the efficiency of the boundary element method are analyzed and discussed.     

Simulating finger phenomena in porous media with a moving finite element method

The non-equilibrium Richards equation is solved using a moving finite element method in this paper. The governing equation is discretized spatially with a standard finite element method, and temporally with second-order Runge–Kutta schemes. A strategy of the mesh movement is based on the work by Li et al. [R.Li, T.Tang, P.W. Zhang, A moving mesh finite element algorithm for singular problems in two and three space dimensions, Journal of Computational Physics, 177 (2002) 365–393]. A Beckett and Mackenzie type monitor function is adopted. To obtain high quality meshes around the wetting front, a smoothing method which is based on the diffusive mechanism is used. With the moving mesh technique, high mesh quality and high numerical accuracy are obtained successfully. The numerical convergence and the advantage of the algorithm are demonstrated by a series of numerical experiments.     

二维多流体域耦合声场的光滑有限元解法

针对声学有限元分析中四节点等参单元计算精度低,对网格质量敏感的问题,将光滑有限元法引入到多流体域耦合声场的数值分析中,提出了二维多流体域耦合声场的光滑有限元解法。该方法在Helmholtz控制方程与多流体域耦合界面的声压/质点法向速度连续条件的基础上,得到二维多流体耦合声场的离散控制方程,并采用光滑有限元的分区光滑技术将声学梯度矩阵形函数导数的域内积分转换形函数的域边界积分,避免了雅克比矩阵的计算。以管道二维多流体域耦合内声场为数值分析算例,研究结果表明,与标准有限元相比,对单元尺寸较大或扭曲严重的四边形网格模型,光滑有限元的计算精度更高。因此光滑有限元能很好地应用于大尺寸单元或扭曲严重的网格模型下二维多流体域耦合声场的预测,具有良好的工程应用前景。     

A nodal discontinuous Galerkin finite element method for nonlinear elastic wave propagation

A nodal discontinuous Galerkin finite element method (DG-FEM) to solve the linear and nonlinear elastic wave equation in heterogeneous media with arbitrary high order accuracy in space on unstructured triangular or quadrilateral meshes is presented. This DG-FEM method combines the geometrical flexibility of the finite element method, and the high parallelization potentiality and strongly nonlinear wave phenomena simulation capability of the finite volume method, required for nonlinear elastodynamics simulations. In order to facilitate the implementation based on a numerical scheme developed for electromagnetic applications, the equations of nonlinear elastodynamics have been written in a conservative form. The adopted formalism allows the introduction of different kinds of elastic nonlinearities, such as the classical quadratic and cubic nonlinearities, or the quadratic hysteretic nonlinearities. Absorbing layers perfectly matched to the calculation domain of the nearly perfectly matched layers type have been introduced to simulate, when needed, semi-infinite or infinite media. The developed DG-FEM scheme has been verified by means of a comparison with analytical solutions and numerical results already published in the literature for simple geometrical configurations: Lamb's problem and plane wave nonlinear propagation.     

Calculation of waveguides by the finite element method using the Bunch-Kaufman procedure

An algorithm for the numerical solution of the problem for eigenvalues in a waveguide in the complete vector statement using the finite-element method and the Bunch-Kaufman procedure for the factorization of the matrix of the derived set of linear algebraic equations was built and implemented.     

Three-dimensional radiative transfer modeling using the control volume finite element method

In the present study, a three-dimensional algorithm for the treatment of radiative heat transfer in emitting, absorbing and scattering media is developed. The approach is based on the utilization of control volume finite element method (CVFEM) which, to the knowledge of the authors, is applied at the first time to 3D radiative heat transfer in participating media. The accuracy of the present algorithm is tested by comparing its predictions to other published works. Comparisons show that CVFEM produces good results. Moreover, this approach permits compatibility with other numerical methods used for computational fluids mechanics problems.     

Solitary wave solutions of selective nonlinear diffusion-reaction equations using homogeneous balance method

An auto-Bäcklund transformation derived in the homogeneous balance method is employed to obtain several new exact solutions of certain kinds of nonlinear diffusion-reaction (D-R) equations. These equations arise in a variety of problems in physical, chemical, biological, social and ecological sciences.     

激光激发黏弹表面波有限元数值模拟

研究黏弹性材料中激光激发的Rayleigh波的传播特征. 考虑到黏弹性材料的黏性特征,在频域内建立黏弹性材料中激光激发Rayleigh波的有限元数值模型. 在验证有限元频域数值模型正确性的基础上,模拟脉冲激光作用在黏弹性材料上激发出Rayleigh波,进而讨论激光激发的黏弹Rayleigh波的传播特征,并比较黏弹性材料与弹性材料中激光激发的Rayleigh波差异,同时分析了材料的黏性劲度参量变化对Rayleigh波特征的影响. 关键词： 表面波 激光超声 有限元方法 黏弹性