二维层状地基波阻板隔振分析

为了研究层状地基中波阻板的隔振效果，基于薄层法在研究层状介质中波的传播问题的高效性、边界单元法处理无限域问题的精确性，结合两种方法的优点，本文采用以薄层法层状半空间基本解答作为格林函数的边界元法，分别对上软下硬和上硬下软两种层状半空间地基中波阻板的隔振效果进行分析。研究表明：增加波阻板的厚度、提高波阻板的弹性模量可以显著提高隔振效果；分层土体厚度和土性对于水平和竖向的位移振幅衰减系数有较大影响。     

A modified scaled boundary finite element method for three‐dimensional dynamic soil‐structure interaction in layered soil

This paper is devoted to the analysis of elastodynamic problems in 3D‐layered systems which are unbounded in the horizontal direction. For this purpose, a finite element model of the near field is coupled to a scaled boundary finite element model (SBFEM) of the far field. The SBFEM is originally based on describing the geometry of a half‐space or full‐space domain by scaling the geometry of the near field / far field interface using a radial coordinate. A modified form of the SBFEM for waves in a 2D layer is also available. None of these existing formulations can be used to describe a 3D‐layered medium. In this paper, a modified SBFEM for the analysis of 3D‐layered continua is derived. Based on the use of a scaling line instead of a scaling centre, a suitable scaled boundary transformation is proposed. The derivation of the corresponding scaled boundary finite element (SBFE) equations in displacement and stiffness is presented in detail. The latter is a nonlinear differential equation with respect to the radial coordinate, which has to be solved numerically for each excitation frequency considered in the analysis. Various numerical examples demonstrate the accuracy of the new method and its correct implementation. These include rigid circular and square foundations embedded in or resting on the surface of layered homogeneous or inhomogeneous 3D soil deposits over rigid bedrock. Hysteretic damping is assumed in some cases. The dynamic stiffness coefficients calculated using the proposed method are compared with analytical solutions or existing highly accurate numerical results. Copyright © 2011 John Wiley & Sons, Ltd.     

一种高阶精度人工边界条件:出平面外域波动问题

针对无限外域中的出平面波动问题,提出一种用于近场波动有限元分析的高阶精度人工边界条件。首先,采用变量分离法求解远场初边值问题,建立了时空全局的精确动力刚度人工边界条件;然后,发展了一种由有理函数近似和辅助变量实现构成的时间局部化方法,并将其应用于动力刚度人工边界条件,得到时间局部的高阶精度人工边界条件;最后,沿人工边界离散高阶精度人工边界条件,并将其与近场集中质量有限元方程耦合,形成对称的时间二阶常微分方程组,采用一种新的显式时间积分方法进行求解。数值算例表明:提出的高阶精度人工边界条件精确、高效、稳定并且容易在现有的有限元代码中实现。     

Explicit finite element artificial boundary scheme for transient scalar waves in two‐dimensional unbounded waveguide

To simulate the transient scalar wave propagation in a two‐dimensional unbounded waveguide, an explicit finite element artificial boundary scheme is proposed, which couples the standard dynamic finite element method for complex near field and a high‐order accurate artificial boundary condition (ABC) for simple far field. An exact dynamic‐stiffness ABC that is global in space and time is constructed. A temporal localization method is developed, which consists of the rational function approximation in the frequency domain and the auxiliary variable realization into time domain. This method is applied to the dynamic‐stiffness ABC to result in a high‐order accurate ABC that is local in time but global in space. By discretizing the high‐order accurate ABC along artificial boundary and coupling the result with the standard lumped‐mass finite element equation of near field, a coupled dynamic equation is obtained, which is a symmetric system of purely second‐order ordinary differential equations in time with the diagonal mass and non‐diagonal damping matrices. A new explicit time integration algorithm in structural dynamics is used to solve this equation. Numerical examples are given to demonstrate the effectiveness of the proposed scheme. Copyright © 2011 John Wiley & Sons, Ltd.     

模拟P波向无穷远域传播的一致边界

将人工边界设置在半无穷层单元和内部有限元区域的交界面上,建立了半无穷层单元的刚度矩阵后,得到了边界节点的动力平衡方程。任意给定激励圆频率,将边界节点系统的动力平衡方程转化为特征值方程。求解特征值方程得出边界节点系统的特征值和特征模态,利用模态叠加原理得到体现左半无穷层单元和右半无穷层单元对内部有限元区域作用的边界矩阵,这就是该文的一致边界。将其与内部有限元区域的刚度矩阵进行组装来模拟无穷远域介质对波的传播作用。最后用数值算例来说明一致边界的精确性和可行性。     

Three-dimensional fundamental solutions for transient heat transfer by conduction in an unbounded medium,half-space,slab and layered media

This paper presents analytical Green's functions for the transient heat transfer phenomena by conduction, for an unbounded medium, half-space, slab and layered formation when subjected to a point heat source. The transient heat responses generated by a spherical heat source are computed as Bessel integrals, following the transformations proposed by Sommerfeld [Sommerfeld A. Mechanics of deformable bodies. New York: Academic Press; 1950; Ewing WM, Jardetzky WS, Press F. Elastic waves in layered media. New York: McGraw-Hill; 1957]. The integrals can be modelled as discrete summations, assuming a set of sources equally spaced along the vertical direction. The expressions presented here allow the heat field inside a layered formation to be computed without fully discretizing the interior domain or boundary interfaces.The final Green's functions describe the conduction phenomenon throughout the domain, for a half-space and a slab. They can be expressed as the sum of the heat source and the surface terms. The surface terms need to satisfy the boundary conditions at the surfaces, which can be of two types: null normal fluxes or null temperatures. The Green's functions for a layered formation are obtained by adding the heat source terms and a set of surface terms, generated within each solid layer and at each interface. These surface terms are defined so as to guarantee the required boundary conditions, which are: continuity of temperatures and normal heat fluxes between layers.This formulation is verified by comparing the frequency responses obtained from the proposed approach with those where a double-space Fourier transformation along the horizontal directions [Tadeu A, António J, Simões N. 2.5D Green's functions in the frequency domain for heat conduction problems in unbounded, half-space, slab and layered media. CMES: Computer Model Eng Sci 2004;6(1):43–58] is used. In addition, time domain solutions were compared with the analytical solutions that are known for the case of an unbounded medium, a half-space and a slab.     

竖向成层介质中标量波传播问题的高精度人工边界条件

为了实现含竖向成层介质以及表面不规则地形场地中标量波传播问题的高效且高精度求解,该文基于连分式展开和扩展的一致边界,建立了一种频域下折线形高精度人工边界条件。通过在每个竖向地层内引入独立的斜角坐标变换,新的人工边界条件可以用于多起伏地表地形条件。新的折线形人工边界在频域下推导,仅含有连分式阶数一个待定实参数,用于调整计算精度,该参数不随外行波的频率和传播角度改变。人工边界条件可以与内域有限元方程无缝耦合,应用简单方便。由于新边界条件的高精度,内域尺寸可以取较小甚至可以直接将人工边界加在结构周围或者地表,从而极大提高计算效率。通过典型数值算例,将人工边界计算模型与有限元大模型的解进行了对比分析,验证了该文提出的折线形人工边界条件的有效性和高精度。     

地震作用下水-轴对称柱体相互作用的子结构分析方法

针对水-轴对称柱体动力相互作用问题,提出了一种地震作用下水-结构相互作用的时域子结构分析方法.基于三维不可压缩水体的波动方程和边界条件,利用分离变量法将其转换为环向解析、竖向和径向数值的二维模型;基于比例边界有限元推导了截断边界处无限域水体的动力刚度方程,并将水体内域有限元方程和人工边界处的动水压力进行耦合,从而得到结...     

一种高精度圆柱形人工边界条件:水-柱体相互作用问题

针对水-柱体动力相互作用问题,提出一种用于模拟无限域水体的圆柱形高精度时域人工边界条件。首先,基于三维可压缩水体的波动方程和边界条件,采用分离变量法建立了时空全局的精确人工边界条件;然后,将其动力刚度表示为外域模型和波导模型人工边界条件动力刚度的嵌套形式;之后,应用时间局部化方法得到时间局部的高精度人工边界条件;最后,离散高精度人工边界条件,并将其与近场有限元方程耦合,形成一种能够采用显式时间积分方法求解的时间二阶常微分方程组。数值算例表明:提出的三维圆柱形高精度人工边界条件精确、高效、稳定。     

Green's functions for transversely isotropic piezoelectric multilayered half-spaces

This paper presents Green's functions for transversely isotropic piezoelectric and layered half-spaces. The surface of the half-space can be under general boundary conditions and a point source (point-force/point-charge) can be applied to the layered structure at any location. The Green's functions are obtained in terms of two systems of vector functions, combined with the propagator-matrix method. The most noticeable feature is that the homogeneous solution and propagator matrix are independent of the choice of the system of vector functions, and can therefore be treated in a unified manner. Since the physical-domain Green's functions involve improper integrals of Bessel functions, an adaptive Gauss-quadrature approach is applied to accelerate the convergence of the numerical integral. Typical numerical examples are presented for four different half-space models, and for both the spring-like and general traction-free boundary conditions. While the four half-space models are used to illustrate the effect of material stacking sequence and anisotropy, the spring-like boundary condition is chosen to show the effect of the spring constant on the Green's function solutions. In particular, it is observed that, when the spring constant is relatively large, the response curve can be completely different to that when it is small or when it is equal to zero, with the latter corresponding to the traction-free boundary condition.     

A spectral finite element for axial-flexural-shear coupled wave propagation analysis in lengthwise graded beam

A new spectral element (SE) is formulated to analyse wave propagation in anisotropic inhomogeneous beam. The inhomogeneity is considered in the longitudinal direction. Due to this particular pattern of inhomogeneity, the governing partial differential equations (PDEs) have variable coefficients and an exact solution for arbitrary variation of material properties, even in frequency domain, is not possible to obtain. However, it is shown in this work that for exponential variation of material properties, the equations can be solved exactly in frequency domain, when the same parameter governs the variation of elastic moduli and density. The SE is formed using this exact solution as interpolating polynomial. As a result a single element can replace hundreds of finite elements (FEs), which are essential for all wave propagation analysis and also for accurate representation of the inhomogeneity. The developed element is used for eliciting several advantages of the gradation, including mode selection, mode blockage and smoothening of stress waves.     

An optimal high‐order non‐reflecting finite element scheme for wave scattering problems

A new finite element scheme is proposed for the numerical solution of time‐harmonic wave scattering problems in unbounded domains. The infinite domain in truncated via an artificial boundary ?? which encloses a finite computational domain Ω. On ?? a local high‐order non‐reflecting boundary condition (NRBC) is applied which is constructed to be optimal in a certain sense. This NRBC is implemented in a special way, by using auxiliary variables along the boundary ??, so that it involves no high‐order derivatives regardless of its order. The order of the scheme is simply an input parameter, and it may be arbitrarily high. This leads to a symmetric finite element formulation where standard C⁰ finite elements are used in Ω. The performance of the method is demonstrated via numerical examples, and it is compared to other NRBC‐based schemes. The method is shown to be highly accurate and stable, and to lead to a well‐conditioned matrix problem. Copyright © 2002 John Wiley & Sons, Ltd.     

SV波超临界角斜入射时层状地基地震动输入在ABAQUS中的实现

有限元数值模拟中,超临界角斜入射SV波作用下层状地基的地震动输入是一个亟待解决的问题。该文采用层状地基频域内精确的动力刚度矩阵(即频域刚度矩阵法)推导了SV波任意角度斜入射下的地震动输入等效节点力计算公式,通过ABAQUS有限元软件模拟SV波斜入射下均匀半空间、层状地基的地震波场,探讨了基于频域刚度矩阵法的层状地基任意角度斜入射地震动输入应用于ABAQUS有限元软件的有效性和准确性。结果表明,采用频域刚度矩阵法可以在ABAQUS中实现层状地基任意角度斜入射SV波地震动输入,且方法具有很高的计算精度,尤其对于SV波超临界角入射情况,有限元数值模拟很好地再现了均匀半空间地表质点的椭圆型运动轨迹和层状地基的行波特点。在此基础上进一步将频域刚度矩阵法与等效线性化方法相结合,解决了考虑土体非线性的层状地基在任意角度平面波入射条件下的地震动输入问题。     

层状饱和半空间中沉积谷地对斜入射平面P1波的三维散射

在作者给出层状饱和场地三维精确动力刚度矩阵和层状饱和半空间中移动荷载动力格林函数基础上,采用间接边界元方法在频域内求解了层状流体饱和场地中沉积谷地对斜入射平面P1波的三维散射问题。该方法的特点在于虚拟移动均布荷载和斜线孔隙水压可以直接施加在沉积与层状饱和半空间交界面而不存在奇异性。该文通过与已有结果的比较验证了方法的正确性,并以均匀饱和半空间和弹性基岩上单一饱和土层中沉积谷地为例进行了数值计算分析。研究表明,沉积谷地对平面P1波的三维散射与二维散射之间存在本质差别,入射角度、孔隙率、饱和土层刚度和饱和土层厚度等参数对沉积谷地附近动力响应有着显著影响。     

一种新型的时域动力无穷单元

针对频域动力无穷元不便与时域有限元结合的缺陷,采用时域的波传函数,推导出一种新型的时域动力无穷单元,并采用UPFs工具对ANSYS进行二次开发,将其嵌入到ANSYS当中。然后采用经典算例对所提出的新型动力无穷单元进行验证,并将其结果与远置边界、粘性边界、粘弹性边界的结果进行比较。验证结果表明,所提出的新型时域动力无穷单元是成功的,并可方便地实现与时域的有限单元结合。     

Fundamental solutions for transient heat transfer by conduction and convection in an unbounded, half-space, slab and layered media in the frequency domain

Analytical Green's functions in the frequency domain are presented for the three-dimensional diffusion equation in an unbounded, half-space, slab and layered media. These proposed expressions take into account the conduction and convection phenomena, assuming that the system is subjected to spatially sinusoidal harmonic heat line sources and do not require any type of discretization of the space domain. The application of time and spatial Fourier transforms along the two horizontal directions allows the solution of the three-dimensional time convection-diffusion equation for a heat point source to be obtained as a summation of one-dimensional responses. The problem is recast in the time domain by means of inverse Fourier transforms using complex frequencies in order to avoid aliasing phenomenon. Further, no restriction is placed on the source time dependence, since the static response is obtained by limiting the frequency to zero and the high frequency contribution to the response is small.

The proposed functions have been verified against analytical time domain solutions, known for the case of an unbounded medium, and the Boundary Element Method solutions for the case of the half-space, slab and layered media.     

Transient analysis of three-dimensional dynamic interaction between multilayered soil and rigid foundation

The study of dynamic soil-structure interaction is significant to civil engineering applications, such as machine foundation vibration, traffic-induced vibration, and seismic dynamic response. The scaled boundary finite element method (SBFEM) is a semi-analytical algorithm, which is used to solve the dynamic response of a three-dimensional infinite soil. It can automatically satisfy the radiation boundary condition at infinity. Based on the dynamic stiffness matrix equation obtained by the modified SBFEM, a continued fraction algorithm is proposed to solve the dynamic stiffness matrix of layered soil in the frequency-domain. Then, the SBFEM was coupled with the finite element method (FEM) at the interface to solve the dynamic stiffness matrices of the rigid surface/buried foundation. Finally, the mixed-variable algorithm was used to solve the three-dimensional transient dynamic response of the foundation in the time domain. Numerical examples were performed to verify the accuracy of the proposed algorithm in solving the dynamic stiffness matrix of the infinite domain in the frequency domain and the dynamic transient displacement response of the foundation in the time domain. Compared with the previous numerical integration technique, the dynamic stiffness matrix in the frequency domain calculated by using the proposed algorithm has higher accuracy and higher efficiency.     

基于有限元与边界元耦合法的波浪力分析

利用有限元与边界元耦合法对三维无界区域中直立圆柱所受的波浪力进行进行计算,把整个求解区域分成内域或外域两部分,在内域采用有限元法,对外域采用边界元法,数值计算的结果与理论解吻合良好,表明该方法有效。     

Eddy-current interaction of a long coil with a slot in a conductive plate

We describe a truncated-domain method for calculating eddy currents in a plate with a long flaw. The plate is modeled as a conductive half-space and the flaw is a long slot with a rectangular cross section. A long two-dimensional (2-D) coil carrying an alternating current is aligned parallel to the slot. The coil impedance variation with frequency is determined for an arbitrary coil location. The electromagnetic field due to a long coil above a conductive half-space can be expressed as integrals of trigonometric functions. For a half-space with a long slot, however, additional boundary conditions must be satisfied at the slot walls. The truncated-domain method makes this possible by recasting the problem in a finite domain; as a result, the Fourier integral is replaced by a series. The domain can be made arbitrarily large, thereby yielding results that are numerically as close to the infinite domain solution as desired. We have used the truncated domain approach to study both eddy-current flaw interactions and edge effects in the limiting case of a very wide and deep slot. We confirmed the theoretical predictions by comparing them with results of a 2-D finite element calculation and of experiments.     

3-D boundary element–finite element method for the dynamic analysis of piled buildings

A boundary element–finite element model is presented for the three-dimensional dynamic analysis of piled buildings in the frequency domain. Piles are modelled as compressible Euler–Bernoulli beams founded on a linear, isotropic, viscoelastic, zoned-homogeneous, unbounded layered soil, while multi-storey buildings are assumed to be comprised of vertical compressible piers and rigid slabs. Soil–foundation–structure interaction is rigorously taken into account with an affordable number of degrees of freedom. The code allows the direct analysis of multiple piled buildings, so that the influence of other constructions can be taken into account in the analysis of a certain element. The formulation is outlined before presenting validation results and an application example.