Some aspects of numerical simulation for Lamb wave propagation in composite laminates

Some aspects of numerical simulation of Lamb wave propagation in composite laminates using the finite element models with explicit dynamic analysis are addressed in this study. To correctly and efficiently describe the guided-wave excited/received by piezoelectric actuators/sensors, effective models of surface-bounded flat PZT disks based on effective force, moment and displacement are developed. Different finite element models for Lamb wave excitation, collection and propagation in isotropic plate and quasi-isotropic laminated composite are evaluated using continuum elements (3-D solid element) and structural elements (3-D shell element), to elaborate the validity and versatility of the proposed actuator/sensor models.     

Numerical analysis of the radiated fields of ultrasonic transducers

A finite element model is implemented for simulating the radiated fields of both planar and curved transducers in acoustic media. The approach is based on general finite element analysis developed for solving the governing equations of elastic wave propagation. The distributions of the wave fields are presented for both the nearfield and farfield regions of the transducer. Three excitation pulses with the same center frequency but different bandwidths are examined and the accuracy is indicated by a comparison of the simulation results for the axial and transverse fields with the analytical results for continuous wave excitation.     

成层介质中平面内自由波场的一维化时域算法

提出了一种弹性水平成层半空间中平面内波动斜入射时自由波场时域计算的一维化有限元方法。首先,基于弹性波在斜入射情形下的传播特点对计算区域进行自动虚拟网格划分。然后将集中质量有限元法和中心差分法相结合建立节点的二维运动方程组,并根据采用的离散化准则和显式有限元法的特点将其转化为一维方程组。求解此方程组,即得到自由场中竖向一列节点的运动。最后根据行波的传播规律得到全部自由波场。以P波为例给出了理论分析和数值算例。结果表明,该方法不仅简单实用,而且具有较高的精度和良好的稳定性。     

波在各向异性介质中传播规律的无网格法数值模拟

采用基于紧支距离基函数近似的配点型无网格方法对波在各向异性层状介质中的传播规律进行了数值模拟,得到了应力波的传播历程,并与冲击载荷作用下的有限元计算结果值进行了比较。该方法所得到的结果与有限元计算的结果吻合较好。说明该方法可以有效地模拟波在各向异性材料中的传播过程。     

Numerical modelling of elastic wave scattering in frequency domain by the partition of unity finite element method

In this paper, we investigate a numerical approach based on the partition of unity finite element method, for the time‐harmonic elastic wave equations. The aim of the proposed work is to accurately model two‐dimensional elastic wave problems with fewer elements, capable of containing many wavelengths per nodal spacing, and without refining the mesh at each frequency. The approximation of the displacement field is performed via the standard finite element shape functions, enriched by superimposing pressure and shear plane wave basis, which incorporate knowledge of the wave propagation. A variational framework able to handle mixed boundary conditions is described. Numerical examples dealing with the radiation and the scattering of elastic waves by a circular body are presented. The results show the performance of the proposed method in both accuracy and efficiency. Copyright © 2008 John Wiley & Sons, Ltd.     

Wave equation system for dilatation and torsion for piecewise continuous materials

Summary The wave equations governing the propagation and scatter of dilatation and torsion are developed from first principles for visco-elastic, piecewise continuous materials. The wave equations describe the conversion of dilatational to torsion waves in materials of spatially varying properties. The conditions for maximum excitation transfer at interfaces are identified. The physical insights provided by the mathematical formulation of the coupled dilatation and torsion wave equations are used to interpret results obtained from finite element numerical model calculations of vibrations in and around inclusions in composite materials excited by a sinusoidal pressure wave or pulse generated by noisy machinery.     

Stabilized nonlocal model for dispersive wave propagation in heterogeneous media

In this paper, a general-purpose computational model for dispersive wave propagation in heterogeneous media is developed. The model is based on the higher-order homogenization with multiple spatial and temporal scales and the C⁰-continuous mixed finite element approximation of the resulting nonlocal equations of motion. The proposed nonlocal Hamilton principle leads to the stable discrete system of equations independent of the mesh size, unit cell domain and the excitation frequency. The method has been validated for plane harmonic analysis and for transient wave motion insemi-infinite domain with various microstructures.This work was supported by the Sandia National Laboratories under Contract DE-AL04–94AL8500, the Office of Naval Research through grant number N00014–97–1-0687, and the Japan Society for the Promotion of Science under contract number Heisei11-nendo 06542.     

Numerical simulation of ultrasonic wave propagation in anisotropic and attenuative solid materials

The axisymmetric elastodynamic finite element code developed is capable of predicting quantitatively accurate displacement fields for elastic wave propagation in isotropic and transversely isotropic materials. The numerical algorithm incorporates viscous damping by adding a time-dependent tensor to Hooke's law. Amplitude comparisons are made between the geometric attenuation in the far field and the corresponding finite element predictions to investigate the quality and validity of the code. Through-transmission experimental measurements made with a 1 MHz L-wave transducer attached to an aluminum sample support the code predictions. The algorithm successfully models geometric beam spreading dispersion and energy absorption due to viscous damping. This numerical model is a viable tool for the study of elastic wave propagation in nondestructive testing applications.     

Comparative study of different nonconserving time integrators for wave propagation in hyperelastic waveguides

This paper addresses the formulation and numerical efficiency of various numerical models of different nonconserving time integrators for studying wave propagation in nonlinear hyperelastic waveguides. The study includes different nonlinear finite element formulations based on standard Galerkin finite element model, time domain spectral finite element model, Taylor–Galerkin finite element model, generalized Galerkin finite element model and frequency domain spectral finite element model. A comparative study on the computational efficiency of these different models is made using a hyperelastic rod model, and the optimal computational scheme is identified. The identified scheme is then used to study the propagation of transverse and longitudinal waves in a Timoshenko beam with Murnaghan material nonlinearity.     

Finite element simulation and visualization of leaky Rayleigh wavesfor ultrasonic NDE

The generation and propagation properties of transient leaky Rayleigh waves are characterized by a two-dimensional finite element model. The displacement vector is used as the primary variable for the solid medium and a potential scalar, which is a replacement for the pressure, is taken as the fundamental variable for the fluid medium. The coupled system of finite element equations are solved in the time domain by direct integration through the central difference scheme. Three configurations are considered: the conversion of a Rayleigh surface wave into a leaky Rayleigh wave, a focused beam probing a fluid/solid interface at the Rayleigh angle, and the interaction of a defocused wave with the interface. The wave velocity in the fluid (water) is lower than the Rayleigh wave velocity in the solid (aluminum). The wave propagation profile in each case is predicted by the model. The finite element model proves to be an effective tool for surface acoustic device design and ultrasonic NDE     

弹性压应力波下直杆分又动力失稳特征值有限元分析

运用有限元特征值分析方法对弹性压应力波作用下直杆分叉动力失稳问题进行了研究。基于应力波理论和相邻平衡准则导出了直杆动力失稳时的有限元特征方程，把弹性直杆的动力失稳问题归结为特征值问题。通过引入直杆动力失稳时的波前约束条件实现了此类问题的有限元特征值解法。     

弹性压应力波下直杆分叉动力失稳特征值有限元分析

运用有限元特征值分析方法对弹性压应力波作用下直杆分叉动力失稳问题进行了研究。基于应力波理论和相邻平衡准则导出了直杆动力失稳时的有限元特征方程,把弹性直杆的动力失稳问题归结为特征值问题。通过引入直杆动力失稳时的波前约束条件实现了此类问题的有限元特征值解法。     

A continuum finite element for single-set jointed media

In response to the need for an advanced computational model for wave propagation in jointed-rock media a new finite element for jointed media with a single set of regularly spaced joints is developed. The element is a numerical implementation of the higher-order homogenization model recently proposed by Murakami and Hegemier. Due to the dispersive effects induced by regularly spaced joints, wave phenomena in jointed media are altered significantly. Therefore, in order to improve the interpretation of seismograms for accurate source identification, it is necessary to develop a higher-order continuum element. The accuracy and efficiency of the new element is investigated by applying it to wave-guide and wave-normal problems of a jointed half-space and by comparing the wave response with that of DYNA2D. The analyses by DYNA2D discretize explicitly the details of the joint microstructure, and are adopted as numerically exact measures for the assessment of the proposed finite element; good correlations were obtained. The validation study also confirmed the importance of wave dispersion for non-linear as well as linear joint responses. Finally, as a more practical application of the proposed element, the problem of a jointed full-space with a cylindrical cavity pressurized by step and pulse loadings was solved. Velocities at several observation points were compared with the numerically exact results of DYNA2D. Similar analyses carried out for elastic isotropic media predicted totally different velocity responses and confirmed the need for the proposed element.     

Finite element modeling of transient ultrasonic waves in linear viscoelastic media

Linear viscoelasticity offers a minimal framework within which to construct a causal model for wave propagation in absorptive media. Viscoelastic media are often described as media with `fading memory,' that is, the present state of stress is dependent on the present strain and the complete time history of strain convolved with appropriate time-dependent shear and bulk stress relaxation moduli. An axisymmetric, displacement-based finite element method for modeling pulsed ultrasonic waves in linear, homogeneous, and isotropic (LHI) viscoelastic media is developed that does not require storage of the complete time history of displacement at every node. This is accomplished by modeling stress relaxation moduli as discrete or continuous spectra of decaying exponentials and relaxation times. Details of the construction and computation of the time-dependent stiffness matrix are presented. As an application of the finite element method, a finite number of exponentials (amplitudes and relaxation times) are employed to represent a typical model for a continuous relaxation spectrum. It is demonstrated that a small number of discrete exponentials are required to model ultrasonic wave propagation of a typical band-limited pulse in a model material accurately. Previous work has shown this model to be consistent with other analytic models for wave propagation in viscoelastic media     

Finite element analysis of elastic–plastic materials displaying mixed hardening

The problem of non-radial loading (in the stress space) of structures of elastic–plastic metals is discussed. A concept of combined isotropic and kinematic hardening is proposed. For a material displaying such mixed hardening, constitutive equations are derived and incorporated in a finite element program for the analysis of elastic–plastic 2D-structures. When determining the numerical solution of the non-linear equations, the program enables the investigator to select an optimal combination of step by step and iteration procedures. The computed results are compared with experimental results and with results from calculations based on the overlay model, Reference 25.     

声表面波换能器激励的有限元仿真

采用有限元法分析了声表面波换能器电极上的激励问题。从声场波动方程、麦克斯韦方程以及压电本构方程出发,利用哈密顿原理,推导了在压电介质中声表面波有限元方程,然后采用Newmark法对有限元方程进行时域变换。分析了换能器电极上的静态电荷分布和动态电荷分布。对压电介质中声表面波振动振幅进行计算并分析了质点振动振幅随深度的变化情况。     

基于谱元法的裂Lamb纹梁波传播特性研究

首次提出用一种无质量弹簧元来模拟含横向裂纹梁的轴弯耦合效应,并结合谱元法分析含裂纹梁内Lamb波的传播特性.由卡氏定理和断裂力学方法推导弹簧元的刚度,以此构建裂纹处的平衡条件和位移协调条件,建立损伤谱元模型.通过和传统的有限元模型进行比较,表明在显著提高计算效率的同时,所提出模型在分析结构固有特性和Lamb波传播特性上都具有较高的精度.在所提出模型的基础上又推导出基于谱元法的能量计算公式,通过裂纹处的能量守恒再次验证损伤模型的正确性和有效性,同时研究结果表明裂纹处转化生成的Lamb 波各模态能量随裂纹深度的变化具有单调性,该结论可以为结构健康监测中定量识别裂纹提供实用依据.     

Three-dimensional finite element analysis of finite deformation micromorphic linear isotropic elasticity

A finite deformation micromorphic materially linear isotropic elastic model is formulated and implemented for three dimensional finite element analysis. The model is based on the kinematics, balance equations and thermodynamic equations proposed by Eringen and Suhubi (1964). The constitutive equations are calculated in the reference configuration, and the resulting stresses are mapped to the current configuration. The balance of linear momentum and the balance of first moment of momentum are linearized to construct the consistent tangent for three dimensional finite element implementation for solution by the Newton–Raphson method. Three dimensional numerical examples are analyzed to demonstrate preliminarily the implementation.     

基于谱元法的复合材料裂纹梁Lamb波传播特性研究

基于谱元法提出了一种弹簧元来模拟复合材料梁由于横向裂纹导致的轴弯耦合效应,分析复合材料裂纹梁中Lamb波的传播特性。由断裂力学的相关知识求得弹簧元的刚度,建立复合材料裂纹梁的损伤谱元模型。通过模拟复合材料裂纹梁内Lamb波传播,并和传统的有限元结果进行比较,验证了所提出模型的可行性和有效性。推导了频域内Lamb波各模态的能量计算公式,裂纹处的能量守恒证明了所提出模型的正确性,同时计算表明复合材料梁中裂纹处反射与透射的Lamb波各模态能量随着裂纹深度的变化规律具有单调性,结论可以为定量识别复合材料梁裂纹提供实用依据。