Wave propagation and dynamic analysis of smoothly graded heterogeneous continua using graded finite elements

The dynamic behavior of smoothly graded heterogeneous materials is investigated using the finite element method. The global variation of material properties (e.g., Young’s modulus, Poisson’s ratio and mass density) is treated at the element level using a generalized isoparametric formulation. Three classes of examples are presented to illustrate this approach and to investigate the influence of material inhomogeneity on the characteristics of wave propagation pattern and stress redistribution. First, a cantilever beam example is presented for verification purposes. Emphasis is placed on the comparison of numerical results with analytical ones, as well as modal analysis for beams with different material gradation profiles. Second, wave propagation patterns are explored for a fixed-free slender bar considering homogeneous, bi-material, tri-layered and smoothly graded materials (steel/alumina), which also provide further verification of the numerical procedures. Comparison of stress histories in these samples indicates that the smooth transition of material gradation considerably alleviates the stress discontinuity in the bi-material system (with sharp interface). Third, a three-point-bending epoxy/glass graded beam specimen is investigated for validation purposes. The beam is graded along the height direction. Stress evolution history at a location of interest is analyzed in detail, which not only reveals the dependence of stress evolution on material gradation direction, but also provides information predictive of potential material failure time for graded beams with different material gradation profiles. Jointly, these three classes of examples provide proper verification and validation for the present numerical techniques.     

成层半空间出平面自由波场的一维化时域算法

提出了一种计算出平面SH波斜入射时弹性水平成层半空间中自由波场时域计算的一维化有 限元方法. 在进行有限元网格划分时，竖向单元取满足有限元模拟精度的任意尺寸，水平向 网格尺寸由时间离散步长和水平视波速确定，并自动进行虚拟网格划分. 基底设置人工边界， 并将波动输入转化为等效荷载施加在边界节点上. 然后将集中质量有限元法和中心差分法相 结合建立节点运动方程，并将水平方向相邻节点的运动用该节点相邻时刻的运动表示，从而 将求解节点运动的二维方程组转化为一维方程组. 求解此方程组，即得到自由场中竖向一列 节点的运动. 最后根据行波传播的特点，可方便地确定全部自由波场. 理论分析和数值算例 表明，该方法具有较高的精度和良好的稳定性.     

Godunov‐type adaptive grid model of wave–current interaction at cuspate beaches

This paper presents a second‐order accurate Godunov‐type numerical scheme for depth‐ and period‐averaged wave–current interaction. A flux Jacobian is derived for the wave conservation equations and its eigensystem determined, enabling Roe's approximate Riemann solver to be used to evaluate convective fluxes. Dynamically adaptive quadtree grids are used to focus on local hydrodynamic features, where sharp gradients occur in the flow variables. Adaptation criteria based on depth‐averaged vorticity, wave‐height gradient, wave steepness and the magnitude of velocity gradients are found to produce accurate solutions for nearshore circulation at a half‐sinusoidal beach. However, the simultaneous combination of two or more separate criteria produces numerical instability and interference unless all criteria are satisfied for mesh depletion. Simulations of wave–current interaction at a multi‐cusped beach match laboratory data from the United Kingdom Coastal Research Facility (UKCRF). A parameter study demonstrates the sensitivity of nearshore flow patterns to changes in relative cusp height, angle of wave incidence, bed roughness, offshore wave height and assumed turbulent eddy viscosity. Only a small deviation from normal wave incidence is required to initiate a meandering longshore current. Nearshore circulation patterns are highly dependent on the offshore wave height. Reduction of the assumed eddy viscosity parameter causes the primary circulation cells for normally incident waves to increase in strength whilst producing rip‐like currents cutting diagonally across the surf zone. Copyright © 2004 John Wiley & Sons, Ltd.     

Experimental study and numerical simulation of mean longshore current for mild slope

In this paper, experimental and numerical studies were performed to investigate the characteristics of longshore current under two mild slopes, the results of which may complement the existing studies, which have mainly focused on steep slopes. The experimental results revealed that the average velocity distribution of the longshore current was significantly different under the two different mild slopes. Under the slope of 1:100, the distribution of the averaged longshore current velocity had a downward concave trend at the nearshore side, while under the slope of 1:40 the trend became convex. In addition, the analysis of the numerical results showed that the distribution of the averaged longshore current velocity was affected by the distribution of the wave height within the surf zone and the bottom friction equations, with the influence of the latter being more significant. For the slope of 1:100, the cross-shore variability of alongshore variability can be calculated using the flow-type bottom friction equation, while for the slope of 1:40, the wave-type bottom friction equation can be used instead. Finally, the secondary breaking wave heights for mild slopes and the distribution of wave set-up are also shown in the study.     

Application of Kalman Filter Finite Element Method and AIC

This paper presents a numerical simulation for application of the Kalman filter finite element method. The Kalman filter is employed frequently for the solution of time series analysis including observation and system noises. Applying the Kalman filter to the finite element method, the present method is capable of the estimation in time and space directions. In this method, the matrix generated by the finite element method is applied to the state transition matrix. Using the Kalman filter finite element method, the characteristics of both the Kalman filter and the finite element method can be strengthened. In this paper, the state transition matrix is based on the shallow water equations which are approximated by the finite element method. This method can estimate the tidal current not only in time but also in space directions.     

流固耦合地震波动问题的显式谱元模拟方法

流固耦合地震波动问题主要研究由流体和固体构成的复杂系统中地震波传播特性及其规律. 传统模拟方法中一般以声波方程、弹性波方程的数值解分别描述理想流体和弹性固体中的波动, 并实时地处理两种不同性质介质之间的相互耦合作用, 数值格式复杂且限制数值模拟精度与计算效率. 本文采用谱元法结合多次透射公式人工边界条件实现了一种流固耦合地震波动问题的高阶显式数值计算方法. 该方法利用了流固耦合问题统一计算框架,可将饱和多孔介质的Biot波动方程分别退化为理想流体的声波方程和弹性固体的弹性波方程. 通过P波垂直入射的水平成层理想流体-饱和多孔介质-弹性固体场地模型、P波斜入射的不规则层状界面以及任意形状界面的理想流体-饱和多孔介质-弹性固体场地模型等三个算例, 与传递函数法解析解以及集中质量有限元法计算结果进行对比分析, 证明了本文方法的正确性与有效性. 数值模拟结果表明, 本文方法相较传统有限元法可以少得多的节点数量获得更高的数值精度, 并且在较宽的频率范围内都能可靠地模拟出流固耦合系统的动力响应, 充分体现出本文方法兼顾高精度、计算效率和复杂场地建模灵活的特点.      

Two-dimensional wave run-up analysis by selective lumping finite element method

A new finite element technique for the analysis of wave run-up is presented in this paper. In this finite element approach, the movement of the shoreline is expressed by that of the nodal points at the wave front, and an auto mesh generation technique is effectively used. The present method is tested by the comparison with the experimental result of a channel with uniform slope, and two numerical examples are reported to show the efficiency of this method. As a final example, the tsunami run-up caused by the 1983 Nihonkai-Chubu earthquake is analysed and compared with actual records of the flooded area.     

Improved finite element forms for the shallow-water wave equations

This paper presents new finite element formulations of the shallow-water wave equations which use different basis functions for the velocity and height fields. These arrangements are analysed with the Fourier transform technique which was developed by Schoenstadt,¹ and they are also compared with other finite difference and finite element schemes. The new schemes are integrated in time for two initial states and compared with analytic solutions and numerical solutions from other schemes. The behaviour of the new forms is excellent and they are also convenient to apply in two dimensions with triangular elements.     

FINITE ELEMENT ANALYSIS FOR THE UNSTEADY NEARSHORE CIRCULATION DUE TO WAVE－CURRENT INTERACTION（II）──TWO－STEP EXPLICIT FINITE ELEMENT METHOD

ＦＩＮＩＴＥＥＬＥＭＥＮＴＡＮＡＬＹＳＩＳＦＯＲＴＨＥＵＮＳＴＥＡＤＹＮＥＡＲＳＨＯＲＥＣＩＲＣＵＬＡＴＩＯＮＤＵＥＴＯＷＡＶＥ－ＣＵＲＲＥＮＴＩＮＴＥＲＡＣＴＩＯＮ（ＩＩ）──ＴＷＯ－ＳＴＥＰＥＸＰＬＩＣＩＴＦＩＮＩＴＥＥＬＥＭＥＮＴＭＥＴＨＯＤＷ...     

弹性半空间二维出平面自由波场的一维化时域算法

提出了一种计算出平面SH波斜入射时弹性半空间自由波场时域计算的一维化有限元方法。首先利用Snell定律确定平面波沿水平方向的传播规律,在用有限元法对弹性半空间进行离散化时,竖向单元尺寸根据波动有限元模拟精度要求确定,而水平向有限元网格尺寸根据水平向波的传播规律和采用的离散时间步长确定,使得有限元离散模型中任意节点的运动可以用水平向相邻节点的运动表示,从而将二维有限元节点运动方程组化为一维的形式。求解此一维方程组,可得到弹性半空间中一列节点的运动,再根据行波的传播规律,可确定全空间自由波场。理论分析和数值算例表明,该方法具有较高的精度和良好的稳定性。     

Experimental investigation of slosh dynamics of liquid-filled containers

This paper is concerned with the experimental studies on the sloshing response of liquid-filled containers. A three-dimensional finite element analysis is carried out for the numerical simulation of this problem. The effects of sloshing are computed in the time domain using Newmark's time integration scheme. A simple experimental setup is designed and fabricated in-house to conduct experiments for measuring some of the basic parameters of sloshing. A sensor device is especially developed to record the free-surface wave heights. Each wave height sensor is a capacitance probe that detects the change in level of liquid (water) precisely with no time lag. The sensors are used in conjunction with a signal-processing unit in which the capacitance values are transduced to a voltage signal between 0 V and 10 V. These wave height sensors simultaneously record the slosh wave height near the periphery of the container wall from 16 predetermined locations to give the free-surface profiles of liquid at desired time steps. The experimental results are compared with those obtained from the present theoretical analysis, and good agreements are observed.     

Numerical modeling of PZT-induced Lamb wave-based crack detection in plate-like structures

This paper presents the numerical modeling and simulations of PZT-induced Lamb wave propagation in plate-like structures by using the spectral finite element method. A novel spectral plate finite element, which can efficiently model the three-dimensional (3D) behavior of Lamb waves, is proposed. In the formulation, linear displacement distributions in the thickness direction are assumed for both the PZT layer and the base plate. A way to avoid the thickness locking is proposed and used in the formulations. Two examples, one for the validation of the proposed two-dimensional (2D) spectral finite element and the other for the demonstration of crack detection in plates, are presented and discussed. The contact between the two faces of crack is considered. Numerical results show that (1) only the anti-symmetric mode is prone to thickness locking thus remedy should be made only on this part, (2) the proposed 2D spectral finite element can adequately model the Lamb wave propagation in plate-like structures and the complex scattering for the crack, and (3) crack location can be well determined by a PZT-induced Lamb wave-based diagnosis algorithm.     

Longshore boundary conditions for numerical wave models

The problem of truncating nearshore finite element wave models is addressed. Incorrect treatment of the artificial boundaries of the model will cause spurious wave reflections. Three methods for dealing with these boundaries: application of constraints, use of the Smith condition and longshore dampers, are proposed. Numerical results show the dampers to be the best method.     

Generalized fourier analyses of the advection–diffusion equation—Part II: two‐dimensional domains

Part I of this work presents a detailed multi‐methods comparison of the spatial errors associated with the one‐dimensional finite difference, finite element and finite volume semi‐discretizations of the scalar advection–diffusion equation. In Part II we extend the analysis to two‐dimensional domains and also consider the effects of wave propagation direction and grid aspect ratio on the phase speed, and the discrete and artificial diffusivities. The observed dependence of dispersive and diffusive behaviour on propagation direction makes comparison of methods more difficult relative to the one‐dimensional results. For this reason, integrated (over propagation direction and wave number) error and anisotropy metrics are introduced to facilitate comparison among the various methods. With respect to these metrics, the consistent mass Galerkin and consistent mass control‐volume finite element methods, and their streamline upwind derivatives, exhibit comparable accuracy, and generally out‐perform their lumped mass counterparts and finite‐difference based schemes. While this work can only be considered a first step in a comprehensive multi‐methods analysis and comparison, it serves to identify some of the relative strengths and weaknesses of multiple numerical methods in a common mathematical framework. Published in 2004 by John Wiley & Sons, Ltd.     

杆杆型冲击拉伸试验装置的二维轴对称弹塑性有限元分析

对带有哑铃状圆柱形试件的杆杆型冲击拉伸试验装置,建立了含有多个物理和几何间断面的二维轴对称弹塑性有限元模型,采用动态增量非线性有限元程序ＡＤＩ－ＮＡ,进行了数值模拟模拟分析,揭示了应力波在试验系统中的传播规律,在弹性框架内初步论证了杆杆型冲击拉伸试验装置测试原理的有效性。     

Finite element analysis for the unsteady nearshore circulation due to wave-current interaction (II) — Two-step explicit finite element method

In this paper, the two-step explicit finite element analysis for the numerical model of the unsteady nearshore circulation proposed in·Ref. [1] and its realization of Fortran program are presented. A circulation has been, clearly shown in the calculated wave-current velocity field, and it is in good agreement with observations. First Received March 27, 1987     

成层饱和介质平面波斜入射问题的一维化时域方法

地震波斜入射下自由场的输入是大型结构抗震分析中亟待解决的问题之一,尤其是成层饱和多孔介质自由场问题,由于问题的复杂性,目前研究甚少. 本文基于Biot提出的饱和多孔介质动力方程,建立了一种新的求解平面波斜入射下基岩上覆饱和多孔介质成层场地自由场分析的一维化时域计算方法. 该方法首先根据Snell定律将饱和多孔介质二维空间问题转化为一维时域问题,通过对深度方向的有限元离散,得到饱和多孔介质波动问题的一维化有限元方程,然后采用单相弹性介质精确人工边界条件模拟基岩半空间的波动辐射和输入特征,通过考虑基岩与饱和多孔介质间透水或不透水边界条件以及不同饱和多孔介质交界面边界条件,形成基岩上覆成层饱和介质系统的整体有限元方程,最后采用中心差分法与Newmark平均加速度近似格式相结合的方法对时间进行离散,得到节点的动力时程的显式表达. 典型场地的地震反应分析表明,本文方法的计算结果与传递矩阵法结合傅里叶变换的计算结果完全吻合,证明了其有效性.      

Finite element method for shallow water equation including open boundary condition

A method to deal with an open boundary condition in the analysis of water surface waves, the tide, etc. by means of the finite element method is proposed in this paper. The present method has two important features relating to the treatment of the open boundary condition. The first feature is to consider the non-reflective virtual boundary condition which has been developed in the numerical wave analysis method. The incident wave conditions without spurious reflected waves can be imposed at the open boundary. The second feature is to identify the amplitude of the components of incident waves in terms of observed water elevations in the field of standing waves. This can be done as a parameter identification based on an optimization problem by applying the conjugate gradient method. The applicability of this method to wave propagation problems is verified by several numerical computations.     

基于不同时间步长时域非结构有限体积法模拟声-弹性耦合问题

介绍一种改进的时域非结构有限体积法(FVM),并将其应用于声-弹性耦合问题。在流体与固体介质中分别求解声波动方程与弹性波方程,根据交界面上的力平衡与质点振速连续条件考虑二者的相互作用。同时考虑双线性四边形单元的线性变化项及常数项,并结合常应变三角形单元处理混合网格问题。分别对三角形单元和四边形单元进行色散分析,给出声波动方程的稳定性条件。在不同介质中采用不同时间步长,提高计算效率。求解弹性波问题、声-弹性耦合问题,结果表明,改进后的方法求解声-弹性耦合问题是有效和准确的,并且具有良好的数值稳定性。     

The analysis of unsteady incompressible flows by a three-step finite element method

This paper describes a three-step finite element method and its applications to unsteady incompressible fluid flows. Stability analysis of the one-dimensional pure convection equation shows that this method has third-order accuracy and an extended numerical stability domain in comparison with the Lax--Wendroff finite element method. The method is cost-effective for incompressible flows because it permits less frequent updates of the pressure field with good accuracy. In contrast with the Taylor-Galerkin method, the present method does not contain any new higher-order derivatives, which makes it suitable for solving non-linear multidimensional problems and flows with complicated boundary conditions. The three-step finite element method has been used to simulate unsteady incompressible flows. The numerical results obtained are in good agreement with those in the literature.