矩形移动荷载作用下路面-双层地基系统三维振动分析

考虑路面和地基之间的相互作用,建立了路面-双层地基的三维模型,将车辆荷载模拟成矩形移动荷载,利用Fourier变换方法对车辆荷载作用下路面-双层地基系统的三维振动问题进行了研究。假设地基为上部弹性土体和下部为饱和土体组成的双层混合结构,整个系统置于刚性基岩上。通过引入势函数,利用Lame分解理论和积分变换方法分别对弹性土层和饱和土层进行求解。在Fourier变换域内,联立路面和下卧双层地基系统的的运动方程,获得了车辆荷载作用下路面-双层地基系统三维振动的位移和孔压响应的积分形式解,并利用IFFT算法和自适应数值积分算法得到了数值计算结果。研究结果可为路基动力响应分析提供参考。     

公路交通荷载作用下分层地基的三维动响应分析

根据高速公路路基的实际情况,建立了三维分层地基模型,分别考虑了刚性路面和柔性路面的情况,利用传递矩阵法获得了置于刚性基岩上的层状地基和分层半无限地基的刚度矩阵,将车辆荷载简化为矩形均布移动荷载,利用Fourier变换方法得到了层状地基系统在交通荷载作用下的三维振动解,并用快速傅立叶逆变换方法得到了数值结果,与模型试验结果进行了比较,显示出了较好的一致性,同时分析了荷载大小、土层性状对响应结果的影响,得出了一些有意义的结论。研究结果可为路面施工设计及路基动力响应分析提供参考。     

双参数黏弹性地基上连续配筋混凝土路面振动参数分析

以连续配筋混凝土路面（CRCP）近年来应用比较广泛的路面结构形式为研究对象,采用考虑地基压缩系数和水平剪切系数的双参数地基模型,建立了考虑地基土体滞回阻尼的黏弹性地基上CRCP的振动微分方程,运用三角级数和Fourier变换得到了简谐、矩形均布荷载作用下路面竖向位移的解答,并利用Fourier逆变换得到了数值结果,较为全面的分析了荷载速度、频率、路面配筋率、板厚以及地基参数对板的动力响应的影响。研究结果可为路面动力响应分析及连续配筋混凝土路面的质量评价提供参考。     

基于传递-反射矩阵法的层状公路结构动力响应研究

利用传递-反射矩阵方法和叠加原理推导了双轮胎振动荷载作用下弹性层状公路结构的动力响应解,并利用Hankel数值逆变换技术进行了数值求解。弹性层状结构的传递-反射矩阵同时具有单个矩阵尺寸小、矩阵元素只含有负指数项2个优点,计算结果不受土层单元厚度与荷载频率大小的影响,可以大幅提高计算效率和精度。通过将数值计算与室内模型试验结果进行比较,验证了方法的正确性。在此基础上,分析了面层厚度及模量对层状公路结构动应力、动位移的影响程度和规律,得出了一些有意义的结论。研究结果可为公路路面、路基结构设计提供参考。     

软黏土在循环荷载作用下动力本构模型的研究

针对软黏土在循环荷载作用下的变形特性,即几乎不存在纯弹性变形阶段及应变趋于零时仍存在能量耗散,结合软土地铁车站振动台模型试验中对饱和软黏土进行的动三轴试验,利用边界面模型理论建立了软黏土的黏弹塑性动力本构模型。并进一步利用该模型对自由场振动台试验进行数值模拟计算,计算与试验的结果基本吻合。该本构模型可用于循环荷载作用下软黏土的动力响应分析。     

饱和多孔介质动力响应移动单元法分析

移动单元法在处理移动荷载下结构动力行为分析方面具有求解高效的优势,但目前针对饱和多孔介质动力响应的移动单元法的研究成果甚少。根据饱和多孔介质u-p格式动力控制方程,利用移动坐标系建立了饱和多孔介质瞬态及稳态动力控制方程的移动单元列式,通过编制相应的计算程序将计算结果与已有文献结果对比验证了算法的正确和有效性。基于移动单元法建立了移动荷载下饱和沥青路面-弹性基层系统计算模型,分析了移动荷载下该模型的瞬态动力响应规律,并与其稳态动力响应进行了对比分析,分析表明,其水动力特性较稳态响应呈现出明显的瞬态效应。基于稳态动力响应结果分析了荷载速度、排水边界、渗透系数对饱和沥青路面动力响应的影响规律,算例研究结果可以为分析水动力作用下沥青路面水稳定性功能损伤机制提供参考。     

交通荷载作用下公路路基工作区深度研究

根据路面不平整实测资料,获得了汽车动荷载的计算模型和参数,并利用此荷载模式来分析公路结构的动力响应。考虑到公路结构的三维分层性状,利用Fourier变换技术和求解层状结构的精确刚度矩阵法,研究了公路结构在移动汽车荷载作用下的动力响应问题。为了方便工程应用,利用Odemark厚度和模量的当量转换公式,将公路结构简化为由路面、路基及地基组成的3层体系,并以此为基础分析了路基动应力的衰减规律,同时重点研究了汽车轴载大小、当量路面厚度及路基模量对路基工作区深度的影响,最后建议了路基工作区深度的定量表达式,研究结果可为公路路基设计提供参考。     

移动荷载下高铁路基段振动加速度频谱衰减特性

采用精细化和多尺度建模技术,建立了设计速度为350 km/h的双线高速铁路轨道-路基-地基非线性耦合系统的真三维动力分析模型,将地基土和路基填筑材料的非线性纳入分析模型中,采用三维黏弹性静-动力统一人工边界技术对地基无限域进行模拟,采用动接触算法模拟底座板底面和基床表层表面之间的相互作用,考虑移动荷载作用前路基-地基中客观的静应力状态对后续动力计算的影响,依托大规模并行计算技术和单元生死技术,模拟了地基的初始应力场生成、路基结构和轨道系统的施工过程,在此基础上模拟了与8辆编组的某型动车组轮对空间位置相对应的、作用于钢轨顶部的压力荷载的移动过程。基于分析结果,归纳了移动荷载作用下高速铁路轨道-路基-地基系统中振动加速度频谱的衰减特性。     

横观各向同性层状地基动应力响应的求解与分析

动力响应问题的求解对于地基在外荷载作用下引起的弹性波动问题研究有重要的意义。本文提出了一种求解横观各向同性层状地基在施加时间简谐荷载作用下任意点的应力响应的算法。此算法利用傅里叶变换将广义平面应变问题频率-空间域的动力方程转化到频率-波数域内,结合对偶变量的引入,利用高精度的精细积分算法对状态方程进行求解,在得到频率-波数域内的位移响应的基础上,利用傅里叶逆变换得到任意点的动应力响应。简谐荷载不仅可以施加在地基表面,而且可以施加在地基内部。对比算例验证了本文算法的准确性,同时对地基各向异性特性、激励频率和阻尼比对动应力响应的影响进行了参数分析,为工程实际提供可靠的数值依据。     

移动荷载下黏弹性半空间体上双层板的动力响应

采用黏弹性半空间体上无限大双层板模型模拟路面结构,通过三重Fourier积分变换得到单位脉冲荷载作用下路面动力响应的Green函数。根据线性系统的叠加原理,利用广义Duhamel积分推导出移动荷载作用下路面动力响应的解析解。采用自适应Simpson法编制了计算奇异、振荡函数的广义积分计算程序,完成了路面动力响应从波数-频率域到时间-空间域的转换。结合算例,对移动荷载作用下路面的振动规律进行研究,讨论上、下层板厚度和板材料的弹性模量对路面动力响应的影响,明确了路面结构的振动特性,研究结果可为路面结构的设计、施工提供理论指导。     

Uncoupled analysis of stabilizing piles in weathered slopes

This paper describes a simplified numerical approach for analyzing the slope/pile system subjected to lateral soil movements. The lateral one-row pile response above and below the critical surface is computed by using load transfer approach. The response of groups was analyzed by developing interaction factors obtained from a three-dimensional nonlinear finite element study. An uncoupled analysis was performed for stabilizing piles in slope in which the pile response and slope stability are considered separately. The non-linear characteristics of the soil–pile interaction in the stabilizing piles are modeled by hyperbolic load transfer curves. The Bishop's simplified method of slope stability analysis is extended to incorporate the soil-pile interaction and evaluate the safety factor of the reinforced slope. Numerical study is performed to illustrate the major influencing parameters on the pile-slope stability problem. Through comparative studies, it has been found that the factor of safety in slope is much more conservative for an uncoupled analysis than for a coupled analysis based on three-dimensional finite element analysis.     

Transient dynamic response of multilayered saturated media subjected to impulsive loadings

This paper presents a stable and efficient method for calculating the transient solution of layered saturated media subjected to impulsive loadings by means of the analytical layer element method. Starting with the field equations based on Biot's linear theory for porous, fluid‐saturated media, and the seepage continuity equation, an analytical layer element for a single layer is established by applying Laplace‐Hankel integral transform. The global stiffness matrix in the transform domain for a layered saturated half‐space subjected to a transient circular patch loading is obtained by assembling the layer elements of each layer. The displacements in the time domain are derived by Laplace‐Hankel inverse transform of the global stiffness matrix. Numerical examples are conducted to verify the accuracy of the method and to demonstrate the influences of type of transient loading, buried depth of loading, permeability, and stratification of materials on the transient response of the multilayered saturated poroelastic media.     

Stochastic finite element methods for the seismic response of soils

Some of the available stochastic finite element methods are adapted and evaluated for the analyses of response of soils with uncertain properties subjected to earthquake induced random ground motion. In this study, the dynamic response of a soil mass, with finite element discretization, is formulated in the frequency domain. The spectral density function of the response variables are obtained from which the evaluation of the root-mean-squared and the most probable extreme values of the response are made. The material non-linearities are incorporated by using strain compatible moduli and damping of soils using an equivalent linear model for stress–strain behaviour of soils and an iterative solution of the response. The spatial variability of the shear modulus is described through a random field model and the earthquake included motion is treated as a stochastic process. The available formulations of direct Monte-Carlo simulation, first-order perturbation method, a spectral decomposition method with Neumann expansion and a spectral decomposition method with Polynomial Chaos are used to develop stochastic finite element analyses of the seismic response of soils. The numerical results from these approaches are compared with respect to their accuracy and computational efficiency. © 1998 John Wiley & Sons Ltd.     

任意荷载下双面半透水边界分数阶导数黏弹性饱和土层一维固结

基于Terzaghi一维固结理论,分析了考虑半透水边界条件的分数阶导数黏弹性饱和土层在随时间变化的任意荷载作用下一维固结问题。首先,应用Laplace变换联立求解饱和土层一维固结微分方程和分数阶Kelvin-Voigt黏弹性本构方程,推导出有效应力和沉降在Laplace变换域内的解析解,采用Crump方法进行Laplace逆变换,得到了时间域内的半解析解。然后将本文得到的半解析解分别退化为半透水边界条件下基于黏弹性假设的一维固结半解析解和双面透水边界条件下基于分数阶黏弹性假设的一维固结半解析解,结果与已有文献的半解析解相同,验证了本研究所提出解的可靠性。最后通过算例分别考察了半透水边界参数、分数阶黏弹性模型参数和荷载参数对饱和土层固结沉降的影响。研究表明,半透水边界条件参数、分数阶次与黏滞系数主要影响饱和土层固结的发展快慢,而饱和土层的最终沉降量主要受到土层压缩模量的影响;另外,饱和土层的固结规律与外荷载变化规律一致。     

Viscoelastic analysis of piles and pile groups

The axial and lateral load–displacement response of single piles and pile groups during the secondary consolidation (creep) phase of deformation is obtained by idealizing the soil as a linear viscoelastic material. The boundary element method is used to model the soil halfspace and the piles, which are assumed to behave elastically, are idealized as beam-columns. By using the correspondence principle, the problem is solved in the Laplace transform domain and the final solution is obtained by numerical inversion of the transformed data. The salient features of the viscoelastic response are identified in several non-dimensional plots which encompass the spectrum of geometrical and material parameters commonly encountered in practice. The results of the analysis compare favourably with the limited data available in the literature.     

Hybrid FE‐DQ for dynamic analysis of multilayered half‐space subjected to concentrated point impulsive loading

A hybrid finite element method and differential quadrature method (DQM) is developed to estimate the dynamic response of two‐dimensional multilayered half‐spaces subjected to impulsive point loading. Nonreflecting absorbing boundary conditions consist of appropriate springs, and dampers are considered. The capabilities of the finite element method for solving boundary value problems with general domain, loading and systematic boundary treatment are combined with accurate and stable time marching capabilities of the DQM to develop an accurate and efficient numerical technique. The capability, efficiency, robustness and convergence of the DQM for solving the dynamic problem are demonstrated through numerical simulations of various half‐spaces with different time increments and layer arrangement. Also, comparison study when using Newmark's time integration scheme for the same problem is done. It can be concluded that the DQM as an unconditionally stable method is suitable for solving such a problem. Also, parametric study is performed to show the effect of the absorbing boundary conditions on the dynamic response. Copyright © 2012 John Wiley & Sons, Ltd.     

爆炸荷载作用下深埋圆形隧洞的饱和黏弹性土-衬砌系统动力响应

假定土骨架服从标准线性固体黏弹性本构关系,研究了深埋圆形隧洞的饱和黏弹性土-弹性衬砌耦合系统在轴对称爆炸作用下的瞬态动力响应。首先,基于饱和土的Biot模型和衬砌的弹性理论,通过引入势函数和Laplace变换,利用弹性衬砌和饱和黏弹性土界面处的连续性条件以及边界条件,得到饱和黏弹性土体和弹性衬砌位移、应力和孔隙水压力等在Laplace变换域中的解析解。其次,利用Laplace数值Crump逆变换得到耦合系统在时间域的动力响应,数值分析了不同土体模型下土体-衬砌耦合系统的径向位移和环向应力以及土体孔隙水压力等。结果表明：对不同土体模型的土体-衬砌耦合系统,其在爆炸载荷作用下的动力响应性态基本一致,但动力响应的振动周期和幅值等具有明显的差异。同时,对于饱和黏弹性土-弹性衬砌系统,土体黏性参数对土体径向位移和孔隙水压力有明显的影响,但对土体环向应力影响较小。     

结构-地基动力相互作用的时域模型

结构-地基动力相互作用是结构地震响应分析及安全评估的一个非常重要课题。基于比例边界有限元法,提出了一种新颖的结构-地基动力相互作用的时域模型,即采用比例边界有限元子域模拟近场有限域部分,采用高阶透射边界模拟远场无限域部分。通过采用连分式展开和引入辅助变量,有限域的动力方程采用高阶的静力刚度矩阵和质量矩阵表示。高阶透射边界精确满足无限远处的辐射边界条件,具有全局精确、时间局部和收敛速度快等优点。它是基于改进的连分式法求解无限域动力刚度矩阵而建立的,在时域里表示为一阶常微分方程组。通过联立有限域和无限域的运动方程,建立了结构-地基相互作用的标准动力学方程,采用Newmark法可直接求解。3个算例结果表明,该算法在时域里比黏弹性边界更精确、有效。     

An operator‐split ALE model for large deformation analysis of geomaterials

Analysis of large deformation of geomaterials subjected to time‐varying load poses a very difficult problem for the geotechnical profession. Conventional finite element schemes using the updated Lagrangian formulation may suffer from serious numerical difficulties when the deformation of geomaterials is significantly large such that the discretized elements are severely distorted. In this paper, an operator‐split arbitrary Lagrangian–Eulerian (ALE) finite element model is proposed for large deformation analysis of a soil mass subjected to either static or dynamic loading, where the soil is modelled as a saturated porous material with solid–fluid coupling and strong material non‐linearity. Each time step of the operator‐split ALE algorithm consists of a Lagrangian step and an Eulerian step. In the Lagrangian step, the equilibrium equation and continuity equation of the saturated soil are solved by the updated Lagrangian method. In the Eulerian step, mesh smoothing is performed for the deformed body and the state variables obtained in the updated Lagrangian step are then transferred to the new mesh system. The accuracy and efficiency of the proposed ALE method are verified by comparison of its results with the results produced by an analytical solution for one‐dimensional finite elastic consolidation of a soil column and with the results from the small strain finite element analysis and the updated Lagrangian analysis. Its performance is further illustrated by simulation of a complex problem involving the transient response of an embankment subjected to earthquake loading. Copyright © 2007 John Wiley & Sons, Ltd.