Absorbing boundary conditions for elastic waves in transversely isotropic media

In the numerical simulation of elastic wave propagation in the solid, it is essential to introduce absorbing boundary conditions to limit the large or unbounded domain of computation. In this paper, the absorbing boundaries for transversely isotropic media are composed of simple first-order partial differential operators, and each of the operators can perfectly absorb a plane wave outgoing at a certain angle. To test the absorbing ability, the reflection coefficient formulas for the quasi-P and quasi-S wave on the absorbing boundary are derived based on the potential functions theory of the elastic wave. Numerical examples show that the absorbing effect is good. The boundary conditions given here have a practical meaning.Supported by National Natural Science Foundation of China.     

Analytical thermo-elastodynamic solutions for a nonhomogeneous transversely isotropic hollow sphere

Summary  The spherically symmetric dynamic thermoelastic problem for a special nonhomogeneous transversely isotropic elastic hollow sphere is formulated by introduction of a dependent variable and separation of variables technique. The derived solution can be degenerated into that for a homogeneous transversely isotropic hollow sphere, a nonhomogeneous isotropic hollow sphere or a solid sphere. The present method, allow to avoid integral transforms, is suited for a hollow sphere of arbitrary thickness subjected to arbitrary spherical symmetric thermal and mechanical loads, and is convenient in dealing with different boundary conditions of dynamic thermoelasticity . The numerical calculation involved is easy to be performed and its results are also presented. Received 30 October 2001; accepted for publication 21 February 2002 The work was supported by the National Natural Science Foundation of China (No. 10172075 and No. 10002016)     

Free axisymmetric vibration of transversely isotropic laminated circular plates

Based on the three dimensional elastic theory, the state equation of the axisymmetric free vibration of transversely isotropic circular plates is established. Taking the advantage of finite Hankel transform, two exact solutions are derived for two boundary conditions, i. e. the rigid-slipping boundary and elastic simply supported boundary. Finally, numerical results are presented and compared with those of FEM. The project is supported by the National Natural Science Foundation of China and the Zhejiang Provincial Natural Science Foundation.     

DYNAMICAL ARTIFICIAL BOUNDARY FOR FLUID MEDIUM IN WAVE MOTION PROBLEMS

无限域流体介质的波动辐射效应是影响海域工程动力反应的重要因素,人工边界是实现此类开放系统近场波动问题数值分析的有效方法.基于位移格式的流体波动理论推导开放域流体介质的人工边界,分别给出一维、二维和三维空间中平面波、柱面波和球面波的流体介质动力人工边界条件,其中一维平面波动人工边界为经典的黏性边界,二维柱面波、三维球面波的人工边界处节点应力与节点速度和加速度成正比,可等效为由阻尼与质量系统构成的人工边界条件.讨论相应的数值模拟技术,给出流体介质动力人工边界在ANSYS软件平台的具体实现方法.近场流体介质动力反应问题的算例表明,所发展的流体动力人工边界对于轴对称波动与非轴对称波动在近场有限域截断处的透射吸收作用的模拟计算精度均较为良好,说明此流体介质人工边界具有较高的可靠性与实用性.所发展的流体介质动力人工边界可较为方便地与大型商用有限元软件结合,可为包括海域地形和海岛在内的海域工程的动力分析提供一定的方法借鉴.     

NUMERICAL SIMULATION OF ELASTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA

A new numerical simulation algorithm is presented for the elastic wave propagation in heterogeneous anisotropic media. We make discretization of the computational domain by using triangular and quadrangular grids. The scheme is based on integral equilibrium at each node to simulate the elastic wave propagation in heterogeneous anisotropic media. The scheme is very flexible in dealing with arbitrary surface topography, inner openings, liquid-solid boundaries and irregular interfaces. Moreover, the free-surface condition of complex geometrical boundaries can be satisfied naturally. This work is an extension of the grid method for the elastic wave propagation in heterogeneous isotropic media, and a quadrangular grid with low computational cost is also introduced. Project supported by the National Natural Science Foundation of China(19672016).     

层状横观各向同性饱和土的非轴对称动力响应

通过方位角的Fourier变换,将圆柱坐标系下横观各向同性饱和土的Biot非轴对称波动方 程转化为一组一阶常微分方程组. 然后基于径向Hankel变换,建立问题的状态方程;求解状态方程后,得到传递矩阵. 进而利用传递矩阵,结合饱和层状地基的边界条件、排水条件及层间接触和连续条件,求解 了任意震源力作用下层状横观各向同性饱和地基频域动力响应问题. 时域解可通过频率的Fourier积分得到.     

Exact solutions for free vibration of transversely isotropic piezoelectric circular plates

By employing the general solution for coupled three-dimensional equations of a transversely isotropic piezoelectric body, this paper investigates the free vibration of a circular plate made of piezoelectric material. Three-dimensional exact solutions are then obtained under two specified boundary conditions, which can be used for both axisymmetric and non-axisymmetric cases. Numerical results are finally presented. The project supported by the National Natural Science Foundation of China (No. 19872060)     

Non-reflecting boundary conditions applicable to general purpose CFD simulators

In simulations of propagating blast waves the effects of artificial reflections at open boundaries can seriously degrade the accuracy of the computations. In this paper, a boundary condition based on a local approximation by a plane traveling wave is presented. The method yields small artificial reflections at open boundaries. The derivation and the theory behind these so-called plane-wave boundary conditions are presented. The method is conceptually simple and is easy to implement in two and three dimensions. These non-reflecting boundary conditions are employed in the three-dimensional computational fluid dynamics (CFD) solver FLACS, capable of simulating gas explosions and blast-wave propagation in complex geometries. Several examples involving propagating waves in one and two dimensions, shock tube and an example of a simulation of a propagating blast wave generated by an explosion in a compressor module are shown. The numerical simulations show that artificial reflections due to the boundary conditions employed are negligible. © 1998 John Wiley & Sons, Ltd.     

Representations for transversely hemitropic and transversely isotropic stress-strain relations

The sets of polynomial stress-strain relations for elastic points which are transversely hemitropic and transversely isotropic are presented as projections of free algebraic modules having 20 and 10 generators, respectively. Complete sets of relations for the projections are presented which allow the sets of interest to be identified as free submodules having 12 and 6 generators, respectively. The results are established using the Cartan decomposition of the representation of the adjoint action of the two-dimensional rotation and orthogonal groups on the space of three-by-three symmetric matrices. The results are compared to known representations for nonlinear transversely isotropic stress-strain relations and for linear, transversely hemitropic and transversely isotropic ones.Work supported in part by National Science Foundation Grant INT-9106519.     

The 3-D stress intensity factor for a half plane crack in a transversely isotropic solid due to the motion of loads on the crack faces

Three-dimensional analysis of a half plane crack in a transversely isotropic solid is performed. The crack is subjected to a pair of normal point loads moving in a direction perpendicular to the crack edge on its faces. Transform methods are used to reduce the boundary value problem to a single integral equation that can be solved by the Wiener-Hopf technique. The Cagniard-de Hoop method is employed to invert the transforms. An exact expression is derived for the mode I stress intensity factor as a function of time and position along the crack edge. Some features of the solution are discussed through numerical results. The project supported by the Guangdong Provincial Natural Science Foundation and the Science Foundation of Shantou University     

A penny-shaped crack in a transversely isotropic piezoelectric solid: Modes II and III problems

An exact analysis of the modes II and III problems of a pennyshaped crack in a transversely isotropic piezoelectric medium is performed in this paper. The potential theory method is employed based on the general solution of three-dimensional piezoelasticity and the four harmonics involved are represented by one complex potential. Previous results in potential theory are then utilized to obtain the exact solution that is expressed in terms of elementary functions. Comparison is made between the current results with those published and good agreement is obtained. The project supported by the National Natural Science Foundation of China (No. 19872060)     

横观各向同性饱和层状地基的三维稳态动力响应

首先引入状态向量,将直角坐标系下横观各向同性饱和土的Biot波动方程转化为一组状态方程,然后基于双重Fourier变换,求解了状态方程,得到传递矩阵.进而利用传递矩阵,并结合饱和地基的边界条件、排水条件及层间接触和连续条件,求解了横观各向同性饱和层状地基的稳态动力响应问题.     

POINT FORCE SOLUTION FOR A TRANSVERSELY ISOTROPIC ELASTIC LAYER

ＰＯＩＮＴＦＯＲＣＥＳＯＬＵＴＩＯＮＦＯＲＡＴＲＡＮＳＶＥＲＳＥＬＹＩＳＯＴＲＯＰＩＣＥＬＡＳＴＩＣＬＡＹＥＲＰＯＩＮＴＦＯＲＣＥＳＯＬＵＴＩＯＮＦＯＲＡＴＲＡＮＳＶＥＲＳＥＬＹＩＳＯＴＲＯＰＩＣＥＬＡＳＴＩＣＬＡＹＥＲ￥ＤｉｎｇＨａｏｊｉａｎｇ（丁...     

横观各向同性介质中弹性波的吸收边界条件

在数值求解固体中的弹性波动问题时，常需引入吸收边界条件来限制大范围或无边界的求解区域，使数值计算得以顺利进行。本文通过合成简单的一阶偏微分算子，给出了横观各向同性介质中弹性波的吸收边界条件，其中每个单一的算子均可完全吸收沿某一角度出射的平面波。文中还基于弹性波的势函数理论，导出了准Ｐ波和准Ｓ波在吸收边界处的反射系数公式，用以检验其吸收能力。本文所给出的吸收条件，形式简单，且算例表明吸收效果良好，因     

Circular crack in a transversely isotropic piezoelectric space under point forces and point charges

In this paper, two kinds of circular crack including external circular crack and penny-shaped crack in a transversely isotropic piezoelectric space are considered. Firstly, we obtain the solution to the problem of an external circular crack in a transversely isotropic piezoelectric space subjected to antisymmetric normal point forces and point charges. Based on this, the solution of one-sided loading of an external circular crack is constructed. Secondly, the real shape of an external circular crack and the opening displacement of a penny-shaped crack under an arbitrary point force and point charge are further obtained. At last, the results are presented in a graphical form. The project supported by the National Natural Science Foundation of China (19872060 and 69982009) and the Postdoctoral Foundation of China     

The numerical solutions of Green's functions for transversely isotropic elastic strata

ＩｎｔｒｏｄｕｃｔｉｏｎＷｈｅｎａｎａｌｙｚｉｎｇｔｈｅｄｙｎａｍｉｃｉｎｔｅｒａｃｔｉｏｎｏｆａｓｏｉｌ＿ｓｔｒｕｃｔｕｒｅｓｙｓｔｅｍ,ｔｈｅｋｅｙｓｔｅｐｉｓｔｏｃａｌｃｕｌａｔｅｔｈｅｄｙｎａｍｉｃｓｔｉｆｆｎｅｓｓｍａｔｒｉｘｏｆｔｈｅｓｏ...     

各向同性和横观各向同性材料的统一点力解

本文根据横观各向同性弹性力学的通解获得了无限体的点力解，由它可以直接退化到各向同性情形的Ｋｅｌｖｉｎ解，利用这个点力解编制的边界元法程序，适用于横观各向同性材料也适用于各向同性材料，因此是真正的统一点力解。还用边界元法计算了两个数值例题。     

3-D DYNAMIC RESPONSE OF TRANSVERSELY ISOTROPIC SATURATED SOILS

Introduction Thedynamicanalysisofsaturatedporoelasticmediaisofbroadmeaninginmanyaspectssuchasseismology,earthquakeengineering,soilmechanicsandgeophysics,etc.Thetransverselyisotropicsaturatedsoilsissolid fluidcoupledtwophasemediainwhichthe soilskeletondisp…     

横观各向同性圆柱体端部问题的三维弹性理论解

从位移的通解出发，用分离变量法得到横观各向同性圆柱体的位移和应力的特征函数展开式，并把位移势函数的解用付里叶积分的形式表示。利用留数运算，该积分解可以转换成类似于特征函数的展开式。通过混合端部边界问题，得到与特征函数解成双正交关系的另一组函数。利用这种双正交关系，可以处理不同的端部边界问题。     

Thermoelastic interaction with energy dissipation in a transversely isotropic thin circular disc

The distribution of stresses due to step input of temperature on the boundaries of a homogeneous transversely isotropic circular disc is investigated by applying Laplace transform technique in the context of generalized theories of thermo-elasticity. The inverse of the transformed solution is carried out by applying a method of Bellman et al. The stresses are computed numerically and presented graphically in a number of figures. A comparison of the results for different theories (CTE, CCTE, TRDTE(GL), TEWED(GN)) and the effect of anisotropy on the stresses are also presented. When the material is isotropic and outer radius of the disc tends to infinity, the corresponding result agrees with that of existing literature.