A spectrally formulated finite element for wave propagation analysis in functionally graded beams

In this paper, spectral finite element method is employed to analyse the wave propagation behavior in a functionally graded (FG) beam subjected to high frequency impulse loading, which can be either thermal or mechanical. A new spectrally formulated element that has three degrees of freedom per node (based upon the first order shear deformation theory) is developed, which has an exact dynamic stiffness matrix, obtained by exactly solving the homogeneous part of the governing equations in the frequency domain. The element takes into account the variation of thermal and mechanical properties along its depth, which can be modeled either by explicit distribution law like the power law and the exponential law or by rule of mixture as used in composite. Ability of the element in capturing the essential wave propagation behavior other than predicting the propagating shear mode (which appears only at high frequency and is present only in higher order beam theories), is demonstrated. Propagation of stress wave and smoothing of depthwise stress distribution with time is presented. Dependence of cut-off frequency and maximum stress gradient on material properties and FG material (FGM) content is studied. The results are compared with the 2D plane stress FE and 1D Beam FE formulation. The versatility of the method is further demonstrated through the response of FG beam due to short duration highly transient temperature loading.     

Boundary-type finite element method for wave propagation analysis

The boundary-type finite element method has been investigated and applied to the Helmholz and mild-slope equations. Four types of interpolation function are examined based on trigonometric function series. Three-node triangular, four-node quadrilateral, six-node triangular and eight-node quadrilateral elements are tested; these are all non-conforming elements. Three types of numerical example show that the three-node triangular and four-node quadrilateral elements are useful for practical analysis.     

A wavelet deconvolution method for impact force identification

The inverse problem of solving for impact force history using experimentally measured structural responses tends to be ill conditioned. A computationally efficient deconvolution method with similarities to Fourier analysis and wavelet analysis is introduced. Force reconstructions obtained using measured acceleration responses from beam and plate models are used to verify the method.     

裂纹扩展过程中线性内聚力模型计算的半解析有限元法

提出了求解基于线性内聚力模型的平面裂纹扩展问题的半解析有限元法,利用弹性平面扇形域哈密顿体系的方程,通过分离变量法及共轭辛本征函数向量展开法,推导了一个环形和一个圆形奇异超级解析单元列式,组装这两个超级单元能准确地描述裂纹表面作用有双线性内聚力的平面裂纹尖端场。将该解析元与有限元相结合,构成半解析的有限元法,可求解任意几何形状和载荷的基于线性内聚力模型的平面裂纹扩展问题。典型算例的计算结果表明本文方法简单有效,具有令人满意的精度。     

PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

IntroductionIt is known that standard finite element procedure is unable to simulate the wavepropagation with high oscillations or gradients in space in the media with reasonableefficiency and accuracy due to the nature of polynomial interpolation approxi…     

Accurate 3-D finite element simulation of elastic wave propagation with the combination of explicit and implicit time-integration methods

One of the big issues in finite element solutions of wave propagation problems is the presence of spurious high-frequency oscillations that may lead to divergent results at mesh refinement. The paper deals with the extension of the new two-stage time-integration technique developed in our previous papers to the solution of wave propagation problems with explicit time-integration methods.The explicit central difference method is used for accurate time-integration of the semi-discrete system of elastodynamics at the stage of basic computations and allows spurious high-frequency oscillations. To filter these oscillations, pre- or/and post-processing (the filtering stage) is applied using a few time increments of the implicit time-continuous Galerkin method with large numerical dissipation.A special calibration procedure is used for the selection of the minimum necessary amount of numerical dissipation (in terms of a time increment) at the filtering stage. In contrast to existing approaches that use a time-integration method with the same dissipation (or artificial viscosity) for all time increments, the new technique yields accurate and non-oscillatory results for wave propagation problems without interaction between user and computer code. The solutions of 3-D wave propagation and impact problems show the effectiveness of the new approach.     

A finite wavelet domain method for wave propagation analysis in thick laminated composite and sandwich plates

An efficient numerical method is developed for the simulation of three dimensional transient dynamic response in thick laminated composite and sandwich plate structures involving very high frequencies and wave numbers. The proposed method incorporates Daubechies wavelet scaling functions for the interpolation of the in-plane displacements with a Galerkin formulation. It further explores the orthonormality and compact support of wavelet scaling functions to produce near diagonal consistent mass matrices and banded stiffness matrices. Hence, an uncoupled equivalent discrete spatial dynamic system is formulated, synthesized and rapidly solved in the wavelet domain using an explicit time integration scheme. The in-plane wavelet interpolation is further combined with an efficient high order layerwise laminate plate theory, that implements Hermite cubic splines for the through-the-thickness approximation of displacement fields. Numerical results are presented on the prediction of guided waves in laminated and thick sandwich composite plates and compared with respective solutions obtained by analytical, semi-analytical and time domain spectral element models. The method yielded higher convergence rates and substantial reductions in computational effort compared to respective time domain spectral finite elements.     

Free vibration and wave propagation analysis of uniform and tapered rotating beams using spectrally formulated finite elements

This paper presents a formulation of an approximate spectral element for uniform and tapered rotating Euler–Bernoulli beams. The formulation takes into account the varying centrifugal force, mass and bending stiffness. The dynamic stiffness matrix is constructed using the weak form of the governing differential equation in the frequency domain, where two different interpolating functions for the transverse displacement are used for the element formulation. Both free vibration and wave propagation analysis is performed using the formulated elements. The studies show that the formulated element predicts results, that compare well with the solution available in the literature, at a fraction of the computational effort. In addition, for wave propagation analysis, the element shows superior convergence.     

固体中短波传播单位分解有限元法的解析积分

针对固体中短波传播数值模拟的单位分解有限元法中单元矩阵积分的被积函数的强烈振荡特性,应用直角坐标系下标准有限元形函数和单元内的波动方向知识提出了一种单元矩阵的解析积分方案。它对于平面三,六,四,八和九节点的直边单位分解有限单元是完全解析的,对于与这些单元相应的曲边单元则是半解析的。数值结果显示所提出的积分方案在计算效率上比高斯-勒让德积分有大幅度提高。     

A depth‐integrated non‐hydrostatic finite element model for wave propagation

This paper presents a two‐dimensional finite element model for simulating dynamic propagation of weakly dispersive waves. Shallow water equations including extra non‐hydrostatic pressure terms and a depth‐integrated vertical momentum equation are solved with linear distributions assumed in the vertical direction for the non‐hydrostatic pressure and the vertical velocity. The model is developed based on the platform of a finite element model, CCHE2D. A physically bounded upwind scheme for the advection term discretization is developed, and the quasi second‐order differential operators of this scheme result in no oscillation and little numerical diffusion. The depth‐integrated non‐hydrostatic wave model is solved semi‐implicitly: the provisional flow velocity is first implicitly solved using the shallow water equations; the non‐hydrostatic pressure, which is implicitly obtained by ensuring a divergence‐free velocity field, is used to correct the provisional velocity, and finally the depth‐integrated continuity equation is explicitly solved to satisfy global mass conservation. The developed wave model is verified by an analytical solution and validated by laboratory experiments, and the computed results show that the wave model can properly handle linear and nonlinear dispersive waves, wave shoaling, diffraction, refraction and focusing. Copyright © 2013 John Wiley & Sons, Ltd.     

水下声波传播过程的时域间断Galerkin有限元方法

本文构建了声压波动方程的改进时域间断Galerkin有限元方法.传统时域连续有限元方法在计算高梯度、强间断特征水中声波传播问题时往往会出现虚假数值振荡现象,这些数值振荡会影响正常波动的计算精度.为了解决这一问题,本文通过引入人工阻尼的方式构建了改进的时域间断Galerkin有限元方法,并针对具有高梯度、强间断特征的多障...     

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

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

Shell element based model for wave propagation analysis in multi-wall carbon nanotubes

A spectrally formulated finite element is developed to study very high frequency elastic waves in carbon nanotubes (CNTs). A multi-walled nanotube (MWNT) is modelled as an assemblage of shell elements connected throughout their length by distributed springs, whose stiffness is governed by the van der Waals force acting between the nanotubes. The spectral element is formulated using the recently developed strategy based on singular value decomposition (SVD) and polynomial eigenvalue problem (PEP). The element can model a MWNT with any number of walls. Studies are carried out to investigate the effect of the number of walls on the spectrum and dispersion relation. The importance of shell element based model over the beam model is established. The zone of validity of the previously developed beam model is also investigated. It is shown that the shell model is required to capture the symmetric Lamb wave modes. It is also shown through numerical examples that the developed element efficiently captures the response of MWNT for Tera-hertz level frequency loading.     

Spectral element based model for wave propagation analysis in multi-wall carbon nanotubes

A spectrally formulated finite element is developed to study elastic waves in carbon nanotubes (CNT), where the frequency content of the exciting signal is at terahertz level. A multi-walled nanotube (MWNT) is modelled as an assemblage of Euler–Bernoulli beams connected throughout their length by distributed springs, whose stiffness is governed by the van der Waals force acting between the nanotubes. The spectral element is developed using the recently developed formulation strategy based on the solution of polynomial eigenvalue problem (PEP). A single element can model a MWNT with any number of walls. Studies are carried out to investigate the effect of the number of walls on the spectrum and dispersion relation. Effect of the number of walls on the frequency response function is investigated. Response of MWNT for terahertz level loading is analyzed for broad-band shear pulse.     

Acoustic wave propagation in inhomogeneous,layered waveguides based on modal expansions and hp-FEM

A new model is presented for harmonic wave propagation and scattering problems in non-uniform, stratified waveguides, governed by the Helmholtz equation. The method is based on a modal expansion, obtained by utilizing cross-section basis defined through the solution of vertical eigenvalue problems along the waveguide. The latter local basis is enhanced by including additional modes accounting for the effects of inhomogeneous boundaries and/or interfaces. The additional modes provide implicit summation of the slowly convergent part of the local-mode series, rendering the remaining part to be fast convergent, increasing the efficiency of the method, especially in long-range propagation applications. Using the enhanced representation, in conjunction with an energy-type variational principle, a coupled-mode system of equations is derived for the determination of the unknown modal-amplitude functions. In the case of multilayered environments, h$h$- and p$p$-FEM have been applied for the solution of both the local vertical eigenvalue problems and the resulting coupled mode system, exhibiting robustness and good rates of convergence. Numerical examples are presented in simple acoustic propagation problems, illustrating the role and significance of the additional mode(s) and the efficiency of the present model, that can be naturally extended to treat propagation and scattering problems in more complex 3D waveguides.     

A wavelet finite-difference method for numerical simulation of wave propagation in fluid-saturated porous media

In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.     

Closed-form solution for shock wave propagation in density-graded cellular material under impact

Density-graded cellular materials have tremendous potential in structural applications where impact resistance is required. Cellular materials subjected to high impact loading result in a compaction type deformation, usually modeled using continuum-based shock theory. The resulting governing differential equation of the shock model is nonlinear, and the density gradient further complicates the problem. Earlier studies have employed numerical methods to obtain the solution. In this study, an analytical closed-form solution is proposed to predict the response of density-graded cellular materials subjected to a rigid body impact. Solutions for the velocity of the impinging rigid body mass, energy absorption capacity of the cellular material, and the incident stress are obtained for a single shock propagation. The results obtained are in excellent agreement with the existing numerical solutions found in the literature. The proposed analytical solution can be potentially used for parametric studies and for effectively designing graded structures to mitigate impact.     

High-amplitude elastic solitary wave propagation in 1-D granular chains with preconditioned beads: Experiments and theoretical analysis

Elastic solitary waves resulting from Hertzian contact in one-dimensional (1-D) granular chains have demonstrated promising properties for wave tailoring such as amplitude-dependent wave speed and acoustic band gap zones. However, as load increases, plasticity or other material nonlinearities significantly affect the contact behavior between particles and hence alter the elastic solitary wave formation. This restricts the possible exploitation of solitary wave properties to relatively low load levels (up to a few hundred Newtons). In this work, a method, which we term preconditioning, based on contact pre-yielding is implemented to increase the contact force elastic limit of metallic beads in contact and consequently enhance the ability of 1-D granular chains to sustain high-amplitude elastic solitary waves. Theoretical analyses of single particle deformation and of wave propagation in a 1-D chain under different preconditioning levels are presented, while a complementary experimental setup was developed to demonstrate such behavior in practice. The experimental results show that 1-D granular chains with preconditioned beads can sustain high amplitude (up to several kN peak force) solitary waves. The solitary wave speed is affected by both the wave amplitude and the preconditioning level, while the wave spatial wavelength is still close to 5 times the preconditioned bead size. Comparison between the theoretical and experimental results shows that the current theory can capture the effect of preconditioning level on the solitary wave speed.