映射半解析边界元法解瞬态弹性动力问题的有关问题

提出了映射半解析边界元法解百轴对称瞬态弹性动力问题时存在解析方向不能正确反映波的传播现象的问题，产生该问题的原因和如何解决该问题的建议。     

基于TI介质模型的凸起地形对平面SH波的放大作用

在模拟地震波在介质中传播过程的研究中,正演模拟是一个重要的方面。该文采用间接边界元法（IBEM）对横观各向同性（TI）不同参数介质的凸起场地中入射SH波的散射问题进行了正演模拟,研究了TI凸起场地对SH波的放大效应。放大效应是通过地表位移与基岩露头位移幅值的比值来体现的。该间接边界元法结合了层状TI半空间精确动力刚度矩阵和均布斜线荷载动力格林函数,具有较高的精度。该文分别给出了TI介质凸起中SH波散射问题在频域和时域内的解答。频域内研究分析表明,土体TI介质参数的变化会对场地位移幅值的放大谱和空间分布产生显著的影响,即改变了凸起部分和土层的动力特性,使得两者动力相互作用发生改变。时域内研究表明,SH波在凸起周围的传播同时依赖于TI介质参数的具体取值和传播方向。时域位移幅值云图清晰地展示入射波、透射波、反射波和散射波的传播过程,凸起角点产生的散射波在TI介质中传播呈现"半椭圆"型,而在各向同性介质中呈现"半圆"型。     

饱和地基中单排桩远场被动隔振研究

针对均质饱和地基中单排桩远场被动隔振问题,基于饱和多孔介质的边界元法,建立了以薄层法(TLM)基本解作为动力Green函数的饱和土半解析边界元法,该方法可有效分析饱和半空间的动力问题,是解决土-结构动力相互作用问题的一种有效工具。在此基础上,根据排桩与周围土体的邻接条件,推导得到了单排桩对Rayleigh波散射的三维边界元方程;运用上述边界元方程,对单排桩的远场被动隔振问题进行了详细的参数分析。结果表明:单排桩能以一个整体对Rayleigh波进行隔振;增加单排桩的桩长、桩径和桩身剪切模量均能有效提高其隔振效果;减小桩间净距也能够有效提高屏障的隔振效果。     

边界元法在环境声学中的应用

边界元法是边界积分方程的数值解法 ,是随着计算机技术的发展而出现的。建立声学边界积分方程分两种方法 :直接法与间接法。本文介绍了边界元法在环境声学中的应用 ,如声屏障和不同情况下道路周围的声场分布、复杂气象条件对声传播的影响的问题等。由于边界元法是半解析半数值解法。在解边界积分方程时会遇到解的存在与唯一性问题。     

结构声辐射计算的边界无网格法

研究函数有限维逼近插值形函数的一般要求,介绍采用移动最小二乘构建无网格插值形函数的方法与步骤;通过配点法将Kirchhoff-Helmholtz边界积分方程离散为受边界条件约束的线性方程组;最后通过分块矩阵法求解约束方程组,得到离散后的声辐射传输模型数值表达式。在计算实例中,分别用边界无网格法和边界元法建立声辐射传输模型进行声场计算,计算声场值与解析值相对比的结果表明,由于边界无网格法插值形函数根据求解情况自行构建,因此更灵活,具有更高的插值和计算精度。     

离散高阶Higdon-like透射边界

基于比例边界有限元法(SBFEM)半离散思想和Higdon透射微分算子提出了一种用于模拟二维层状介质标量波传播的高效离散高阶Higdon-like透射边界。对无限介质边界进行迦辽金有限元离散后,描述标量波的偏微分方程转换为局部坐标系下半离散矩阵方程组;然后使用高阶Higdon透射算子和辅助变量,在时域内得到了一个阶数不超过2阶的离散高阶透射边界。透射边界是由一组常微分方程构成,可以采用通常的时步积分方法求解,它在截断边界上非局部,在时间域局部。算例表明:该文提出的透射边界的计算精度可以随着辅助变量的增加而提高,但计算量却呈线性化增加,因而计算效率较全局方法有了显著提高。另外,由于该文的边界条件是直接建立在离散节点上的,所以它很方便与近场有限单元法耦合。     

列车荷载下竖向非均匀地基波阻板主动隔振分析

对列车模型简化计算,基于薄层法研究波传播问题高效性与边界元法处理无限域问题精确性,建立薄层法基本解答作为格林函数边界元法。针对三维竖向非均匀地基模型,对WIB主动隔振体系进行参数研究。结果表明,在轨道交通荷载作用下,增加WIB剪切模量、厚度、宽度及选较小埋深均可增大WIB隔振体系隔振效果;非均匀系数对隔振效果略有影响。     

爆破中隔振孔对Rayleigh波的隔振效果分析

以理论分析为基础,研究了隔振孔对Rayleigh波的主动隔振计算流程,在弹性体半解析边界元法的基础上,推导出弹性半空间地层介质中隔振孔减振的边界元方程,并用MATLAB编制相关计算程序,计算并分析了弹性半空间介质中单排隔振孔的尺寸等参数对减振效果的影响规律并用实际案例进行验证。结果表明:随着隔振孔半径的增大,有效减振区域面积逐渐增大,减振效果逐渐上升,并且隔振孔中心线位置处的减振效果最佳。当采用单排隔振孔进行隔振时,隔振孔的半径对整体的减振效果有重要影响,为求获得良好的减振效果以及扩大减振区域面积,应尽可能使用孔径较大的隔振孔,隔振孔孔距应控制在一定范围内,并尽可能降低。     

爆破中隔振孔对Rayleigh波的隔振效果分析

以理论分析为基础,研究了隔振孔对Rayleigh波的主动隔振计算流程,在弹性体半解析边界元法的基础上,推导出弹性半空间地层介质中隔振孔减振的边界元方程,并用MATLAB编制相关计算程序,计算并分析了弹性半空间介质中单排隔振孔的尺寸等参数对减振效果的影响规律并用实际案例进行验证。结果表明:随着隔振孔半径的增大,有效减振区域面积逐渐增大,减振效果逐渐上升,并且隔振孔中心线位置处的减振效果最佳。当采用单排隔振孔进行隔振时,隔振孔的半径对整体的减振效果有重要影响,为求获得良好的减振效果以及扩大减振区域面积,应尽可能使用孔径较大的隔振孔,隔振孔孔距应控制在一定范围内,并尽可能降低。     

各向异性弹性层中SH波传播的边界元方法

本文利用坐标变换下的映象法求出了底部为刚性边界的各向异性弹性层中SH波传播问题的格林函数,并由此建立了求解分层介质问题的边界元方法。文中给出了在稳态线荷载作用下带有圆孔的弹性层中反应的计算实例。     

无限域波动问题的有限元模型

本文基于半解析波动方程,推导了无限域粘弹性人工边界单元。在有限元波动模型中,统一了节点计算稳定性问题;提出了合理的人工边界反射系数公式,为有限元模型提供了理论基础。SH波数值分析表明本文有限元模型具有多向波动透射能力。     

A frequency-domain approach for modelling transient elastodynamics using scaled boundary finite element method

This study develops a frequency-domain method for modelling general transient linear-elastic dynamic problems using the semi-analytical scaled boundary finite element method (SBFEM). This approach first uses the newly-developed analytical Frobenius solution to the governing equilibrium equation system in the frequency domain to calculate complex frequency-response functions (CFRFs). This is followed by a fast Fourier transform (FFT) of the transient load and a subsequent inverse FFT of the CFRFs to obtain time histories of structural responses. A set of wave propagation and structural dynamics problems, subjected to various load forms such as Heaviside step load, triangular blast load and ramped wind load, are modelled using the new approach. Due to the semi-analytical nature of the SBFEM, each problem is successfully modelled using a very small number of degrees of freedom. The numerical results agree very well with the analytical solutions and the results from detailed finite element analyses.     

钢轨垂直振动模态导波频散曲线、波结构及检测应用*

频散曲线和波结构是钢轨导波检测时设计检测方法和分析回波波形的基础。首先介绍了半解析有限元计算任意截面弹性波导中导波传播的频散曲线和波结构的基本方法。然后应用该方法求解了U60型钢轨中垂直振动模态的导波频散曲线和波结构,并讨论了频散曲线和波结构在钢轨轨头缺陷导波检测中的应用。最后采用模态力锤实验及时频联合分析的方法验证了频散曲线数值计算结果。钢轨轨头横向裂纹垂直振动模态导波检测的有限元数值模拟和实验表明该模式导波是检测钢轨轨头缺陷的有效检测模式波型。     

A finite element Laplace transform solution technique for the wave equation

A technique is described for the solution of the wave equation with time dependent boundary conditions. The finite element solution accompanied by the numerical Laplace inversion process seems to be an efficient procedure to treat such problems. The programming involved is straightforward in the sense that numerical Laplace inversion routines can be directly used as a time integration procedure after obtaining standard finite element differential equation solutions in the transformed domain. Some results are presented for one- and two- dimensional applications, such as wave propagation in longitudinal bars and wave propagation in harbours.     

多边形比例边界有限元非线性高效分析方法

比例边界有限元法作为一种高精度的半解析数值求解方法,特别适合于求解无限域与应力奇异性等问题,多边形比例边界单元在模拟裂纹扩展过程、处理局部网格重剖分等方面相较于有限单元法具有明显优势。目前,比例边界有限元法更多关注的是线弹性问题的求解,而非线性比例边界单元的研究则处于起步阶段。该文将高效的隔离非线性有限元法用于比例边界单元的非线性分析,提出了一种高效的隔离非线性比例边界有限元法。该方法认为每个边界线单元覆盖的区域为相互独立的扇形子单元,其形函数以及应变-位移矩阵可通过半解析的弹性解获得;每个扇形区的非线性应变场通过设置非线性应变插值点来表达,引入非线性本构关系即可实现多边形比例边界单元高效非线性分析。多边形比例边界单元的刚度通过集成每个扇形子单元的刚度获取,扇形子单元的刚度可采用高斯积分方案进行求解,其精度保持不变。由于引入了较多的非线性应变插值点,舒尔补矩阵维数较大,该文采用Woodbury近似法对隔离非线性比例边界单元的控制方程进行求解。该方法对大规模非线性问题的计算具有较高的计算效率,数值算例验证了算法的正确性以及高效性,将该方法进行推广,对实际工程分析具有重要意义。     

A dynamic infinite element for three-dimensional infinite-domain wave problems

A three-dimensional dynamic infinite element which satisfies the following requirements: (1) displacement compatibility on the interface between finite and infinite elements; (2) definition of the wave propagation and amplitude attenuation behaviours in the infinite element using wave propagation functions; (3) convergence of the generalized integrals related to mass and stiffness matrices of the infinite element: and (4) displacement continuity along the common boundary of neighbouring infinite elements in the case of simulating multiple material layers or multiple wave numbers within the foundation, is presented in this paper. Since P-waves, S-waves and R-waves in the foundation can be simulated Simultaneously in the present infinite element, the seismic response of an arch-dam-foundation system, especially a thin double-curvature arch-dam-foundation system where the boundary element loses its competitive capacity with the finite element, can be economically calculated by coupling this infinite element with conventional finite elements. The good accuracy obtained using the present infinite element and finite element coupling model to simulate foundation wave problems has been proven by comparing the current numerical results with previous analytical results.     

“Causal” shape functions in the time domain boundary element method

Since failing to respect the causality condition has been identified as one of the main sources of inaccuracies in the time domain boundary element method for elastodynamics and scalar wave propagation problems, in this contribution new shape functions are investigated, which permit a more accurate simulation of the continuous propagation of wave fronts. The performance of these shape functions in 2D scalar wave propagation problems is tested both for the potential (displacement) and for the time gradient (velocity) equations. Analytical time integrations are developed and numerical results are presented.     

基于水平集算法的扩展比例边界有限元法研究

扩展比例边界有限元法在裂纹贯穿单元采用Heaviside阶跃函数描述裂纹面两侧的不连续位移,在裂尖则采用半解析的比例边界有限元描述奇异应力场。该方法具有无需预先知道裂尖渐进场的形式,无需采用特殊的数值积分技术直接生成裂尖刚度阵,对多种应力奇异类型可根据定义直接求解广义应力强度因子的特点。该文将扩展比例边界有限元法与水平集方法相结合,进一步发展了扩展比例边界有限元法,并将其应用于解决裂纹扩展的问题。在数值算例中,通过编写完整的MATLAB分析计算程序,求解了单边缺口的三点弯曲梁和四点剪切梁的裂纹扩展问题,计算结果显示扩展比例边界有限元法能有效地预测裂纹轨迹和荷载-位移曲线。通过参数敏感性分析,还可得出该方法具有较低的网格依赖性,且对裂纹扩展步长不敏感。     

Modelling crack propagation in reinforced concrete using a hybrid finite element–scaled boundary finite element method

A previously developed hybrid finite element–scaled boundary finite element method (FEM–SBFEM) is extended to model multiple cohesive crack propagation in reinforced concrete. This hybrid method can efficiently extract accurate stress intensity factors from the semi-analytical solutions of SBFEM and is also flexible in remeshing multiple cracks. Crack propagation in the concrete bulk is modelled by automatically inserted cohesive interface elements with nonlinear softening laws. The concrete–reinforcement interaction is also modelled by cohesive interface elements. The bond shear stress–slip relation of CEB-FIP Model Code 90 and an empirical confining stress–crack opening relation are used to characterise slip and split failure at the concrete–reinforcement interface, respectively. Three RC beams were simulated. The numerical results agreed well with both experimental and numerical results available in the literature. Parametric studies demonstrated the importance of modelling both slip and split failure mechanisms at the concrete–reinforcement interface.