竖向非均质饱和地基中埋置扭转荷载的动力响应

邹新军 贺琼 覃玉兰

摘   要：为探讨竖向非均质饱和地基中埋置简谐扭转荷载的动力响应，考虑地基为饱和半空间，并假定土体剪切模量随深度呈指数函数非线性分布，基于Biot固结理论与弹性动力学原理，建立饱和地基扭转振动的动力微分方程，引入边界条件并利用Hankel变换方法求解获得变换域内的剪应力、切向位移表达式，进而通过Hankel逆变换求得饱和地基内部的剪应力、切向位移解答，据此基于Mathematica编制出相应计算程序，并通过参数分析发现：地基土剪应力、切向位移沿径向呈现出明显的波动规律，曲线波动频率随荷载频率增大而增大;荷载作用面土体切向位移最大，剪应力发生突变，主要影响范围是荷载作用面上、下各两倍荷载作用半径（2a）区域;土体最大切向位移随荷载埋深h增大而减小，h = 2a时最大切向位移下降90%，h > 4a时，最大切向位移近似为零.

关键词：饱和地基;埋置扭转荷载;动力响应;剪切模量;Hankel变换

中图分类号：TU471.4                   文献标志码：A

Dynamic Response of Buried Torsional Load in Vertically

Non-Homogeneous Saturated Soil

ZOU Xinjun？覮，HE Qiong，QIN Yulan

（College of Civil Engineering，Hunan University，Changsha 410082，China）

Abstract：To discuss the dynamic response of buried harmonic torsional load in vertically non-homogeneous saturated soil，the soil shear modulus was assumed as a nonlinear distribution with the depth defined as an exponential function，and the dynamic differential equations of half space were established by using the Biot's consolidation theory and elastic-dynamic theory. The expressions of stress and tangential displacement in the Hankel transform domain were then acquired by solving the dynamic differential equations using the method of Hankel transform，and the true stress and tangential displacement can be obtained by Hankel inverse transformation as well. The corresponding calculation program by Mathematica was compiled based on the obtained solutions. A detailed parameter analysis completed by the program indicates that，the stress and tangential displacement of the soil show obvious fluctuations with the change of radius，and the frequency of fluctuant curves increases with the loading frequency. In addition，the maximum tangential displacement of soil occurs and there is a sharp change of stress on the loading surface. The influence range of the buried harmonic torsional load is about two times the action radius to the loading surface. Furthermore，the largest tangential displacement of the soil is negatively correlated with the depth of buried load，and it is reduced by 90% when the depth of buried load is 2 times the action radius. When the depth of buried load is greater than 4 times the action radius，the largest tangential displacement is approximately equal to zero.

Key words：saturated soil;buried torsional load;dynamic response;shear modulus;Hankel transform

彈性半空间的动力响应问题在岩土工程中具有举足轻重的地位，该问题在土-结构物的动力相互作用、地震工程、基础振动等领域均引起了国内外学者的关注.随着我国经济发展，近海工程建设如火如荼，如跨海大桥、海上风力发电机以及海上平台等，这些结构物受力较复杂，不仅要承受自然环境产生的风力、波浪力以及地震荷载，还要承受使用过程中的运转荷载，其不可避免地会受到扭转荷载的影响. 同时，实际工程中基础大多具有一定埋深，基础承受扭转动荷载时，周围土体亦会承受基础所传递的埋置扭转动荷载，因此研究弹性半空间中埋置扭转荷载的动力响应问题，有助于进一步掌握地基基础的承载变形特性.

Lamb[1]最先对弹性半空间的动力学问题进行研究，其分别考虑了半空间表面作用和内部埋置的法向线荷载与点荷载等四种典型荷载. Chao[2]研究了弹性半空间表面受水平、竖向荷载时的动力响应问题.王贻荪[3]研究了竖向集中简谐荷载作用于弹性半空间表面的问题，采用拉普拉斯变换方法获得了特定条件下半空间表面位移的精确解. Reissner和Sagoci[4]首次研究了弹性半空间表面基础受扭转荷载作用下的响应问题. Rahman[5]在Reissner的理论基础上进一步对弹性半空间内部埋置刚性板的扭转振动问题进行了研究. 自Biot[6] 提出饱和多孔介质波的传播理论，并推导出土体的多维固结方程以来，不少学者开始在此理论基础上研究饱和半空间的动力学问题. Philippacopoulos[7]考虑多孔弹性半空间中埋置点源荷载，利用傅立叶变换与汉克尔变换方法求解地基土的位移. 陈胜立等[8]对饱和地基中埋置简谐竖向集中荷载下的动力响应问题进行研究，分析了饱和地基土表面竖向位移的变化规律. Chen等[9]研究饱和地基中埋置任意荷载的三维动力响应问题，求解了弹性半空间中的应力、位移分量. 张智卿等[10]研究了均质非饱和滞回阻尼土层中弹性支撑桩的扭转振动. 陈刚等[11]对半空间饱和土体中埋置扭转简谐荷载的稳态响应问题进行研究，分析了土体参数对响应的影响规律.

以上研究均視地基为各向同性弹性半空间，但实际上地基土在沉积过程中会受到各个方向有效应力的影响，具有一定的各向异性，主要有径向非均质性和竖向非均质性两种情况. Naggar[12]研究了径向非均匀地基中土体的竖向、扭转动力响应. Zhang等[13]考虑地基土的径向非均匀性（将桩周土沿径向分为n圈层，每个圈层为剪切模量不同的均质各向同性弹性体），研究弹性半空间中单桩的扭转振动响应问题. Awojobi等[14]考虑弹性半空间的竖向非均质性，假定土体剪切模量随深度线性变化，对半空间的平面应变与轴对称问题进行了研究. 邹新军等[15]基于桩侧土双折线模型，探讨了地基土分层时单桩的受扭性状. 王国才等[16]对层状地基中单桩的扭转振动问题进行研究，采用积分变换及Muki 虚拟桩的方法求解并分析了桩土参数对动力响应的影响. Wroth[17]等认为土壤沉积时的地质环境和荷载历史随着时间而发生变化，故假定土体剪切模量随深度呈非线性分布更符合实际. Rajapakse[18]假定土体剪切模量随深度分别呈幂函数与指数函数非线性分布，采用经典变分定理的方法来分析荷载传递问题. 邹新军等[19]考虑地基土的剪切模量随深度非线性分布，对单桩静力受扭问题进行了弹塑性分析，但相应的动力扭转分析仍属鲜见.

天然地基土通常在竖向呈现非均匀性，其剪切模量、渗透系数及孔隙率等均随深度发生变化，为简化问题，本文主要考虑剪切模量随深度的非线性变化，即假定地基土的剪切模量随深度呈指数函数非线性分布，基于Biot固结理论与弹性动力学原理，建立饱和地基土中作用简谐扭转荷载下的动力微分方程，引入边界条件并利用Hankel变换与逆变换求解获得饱和地基土中的剪应力与切向位移解答，据此基于Mathematica编制出相应的计算程序，进一步进行参数分析并讨论土体参数对地基土中埋置简谐扭转荷载动力响应的影响规律.

1   基本扭转动力微分方程的建立

弹性半空间中埋置简谐扭转荷载问题可在柱坐标系下建立简化分析计算模型，如图1所示，其中：弹性半空间地基内部z = h处作用埋置简谐扭矩T0eiωt （ω为频率，i = （-1）1/2），荷载作用范围是半径为a的圆形区域.

许多学者在研究扭转问题时，为进一步简化问题均做出一定假定，如张智卿[10]、Naggar[12]、Rajapakse[18]等，为建立问题的基本方程，本文参照已有文献特做出如下假定：

3） 地基土的扭转振动为小变形，只考虑土骨架切向位移uθ（r，z，t）eiωt和孔隙水相对于土骨架的切向位移wθ（r，z，t）eiωt，为书写方便下文省略简谐因子eiωt.

4） 假定简谐扭转荷载作用范围内，扭转剪应力F（r）与半径成正比例关系，则荷载作用面处扭转剪应力分布函数为：

2   基本扭转动力微分方程的求解

为求解剪切模量呈指数函数非线性分布的地基土扭转动力响应，将剪切模量分布式（1）代入饱和地基土的基本扭转动力微分方程式（9）中并化简：

3   边界条件与待定系数求解

为了求解所得解答式（16）与式（17）中的待定系数A、B，首先需要确定相应的边界条件. 由图1所示弹性半空间中埋置简谐扭转荷载的简化计算分析模型可知，荷载作用面将半空间分为区域（1）与区域（2）两部分，根据基本扭转动力微分方程的解式（16）以及土体的剪应力分布式（17）可得区域（1）、（2）的切向位移、剪应力表达式分别为：

将待定系数表达式（24）代入到饱和地基土的切向位移式（18）与剪应力式（19）中可获得Hankel变化域内的剪应力、切向位移解答，并将该解答进行Hanke逆变换即可得到饱和地基土的应力位移.

根据上述获得的地基土应力位移解答，采用Mathematica编制出Hankel逆变换的数值积分程序，用以计算埋置简谐扭转荷载作用下竖向非均质饱和地基中的最终剪应力与切向位移.

4   结果验证与参数分析

4.1   结果验证

为验证本文推导结果的正确性，令α→0，则土体剪切模量Gs不随深度z变化，此时问题退化为均质饱和半空间中埋置扭转荷载作用的动力响应问题，与陈刚等[11]的理论解进行对比（如图2），计算参数为：地表处剪切模量μ = 9.4 × 106 Pa，孔隙率n = 0.4，土颗粒密度ρs = 2 650 kg/m3，水的密度ρw = 1 000 kg/m3，土体动力渗透系数kd = 10-7 m/s，扭转荷载频率f = 1、埋深 h = 8，计算深度 z = 8.8.

从图2可看出，本文退化结果与已有理论解吻合较好，从而验证了本文解的正确性.

4.2   参数分析

为了分析影响饱和地基扭转动力响应的主要参数及其规律，下面分别探讨简谐扭转荷载的频率与埋深、土体剪切模量非均匀系数及地基土的动力渗透系数等参数对地基扭转动力响应的影响.

以饱和黏性土地基为例[11]，其物理力学参数为：地表处剪切模量μ = 9.4 × 106 Pa，孔隙率n = 0.4，土颗粒密度ρs = 2 650 kg/m3，水的密度ρw = 1 000 kg/m3. 算例中计算得到的应力、位移皆取幅值并除以k即

4.2.1   简谐扭转荷载无量纲频率 f的影响

为探讨简谐扭转荷载无量纲频率f的影响，分别取f为1.00、1.25、1.50、1.75、2.00、2.25、2.50及2.75，荷载作用半径a = 1 m，埋深 h = 1，土体剪切模量非均匀系数α = 1，土体动力渗透系数kd = 10-7 m/s，根据编制的计算程序分别计算深度z1 = 0.5与 z2 = 1.5处土体的剪应力以及切向位移，获得如图3、图4所示的土体切向位移、剪应力随无量纲半径的变化曲线.

由图3、图4可看出：饱和地基土的剪应力与切向位移沿径向呈现出明显的波动变化规律，随着简谐扭转荷载无量纲频率f的增大，曲线波动频率逐渐增大.地基土的切向位移幅值随荷载频率的增大先增后减，即存在某一临界荷载频率，使得切向位移幅值最大，与共振理论相一致.总体趋势上，地基土的应力幅值随荷载频率的增大而增大.在同一深度处，当r > a时，地基土的应力、位移幅值随着到荷载中心距离的增大而减小，与地基中应力波传播时能量逐渐降低的规律一致.

4.2.2   简谐扭转荷载埋深h的影响

为探讨简谐扭转荷载埋深的影响，分别取 h为0.5、1.0、2.0及4.0，荷载频率 f = 1，其余参数同上，根据所编制的计算程序分别计算r > a处土体不同深度处的剪应力和切向位移，获得如图5、图6所示的地基土切向位移、剪应力随深度变化曲线.

由图5、图6可看出：荷载作用面以上，土的剪应力、切向位移随荷载埋深h的增大而增大，荷载作用面以下，土体的剪应力与切向位移随荷载埋深h的增大而减小，荷载作用面土体的剪应力发生突变，切向位移达到峰值.当z - h≥2a时，土体的剪应力、切向位移近似为零，说明埋置荷载主要影响其作用面上、下各两倍荷载作用半径（2a）区域.

根据所编制的程序分别计算荷载频率 f为0.1、0.5、1.0及1.5时，土体不同深度处的剪应力与切向位移，获得如图7所示的土体最大切向位移与荷载埋深h之间的关系曲线.

由图7可看出：土的最大切向位移随荷载埋深h的增大而减小，荷载埋深小于两倍荷载作用半径（2a）时，曲线急剧下降，h = a时最大切向位移下降75%，h = 2a时最大切向位移下降90%，h > 4a时，饱和地基土的最大切向位移近似为零.

4.2.3   剪切模量非均匀系数α的影响

为探讨地基土剪切模量分布形式的影响，分别取剪切模量指数函数式（1）中的非均匀系数α为0.1、0.5、1.0及1.5，荷载频率 f = 1，埋深 h = 1，其余参数同上，计算获得如图8、图9所示深度 z1 = 0.5处地基土的剪应力、切向位移随半径的分布曲线.

由图8、图9可看出：随地基土剪切模量非均匀系数α的增大，土中剪应力、切向位移逐渐减小，系数α每增加0.5个单位，剪应力峰值下降约23%，切向位移峰值下降约40%，这是因为系数α越小，地基土越接近均质状态，其变形越大，说明相同荷载条件下，剪切模量随深度呈指数函数分布的非均质地基比均质地基变形小.

4.2.4   地基土动力渗透系数kd的影响

为探讨饱和地基土动力渗透系数的影响，分别取kd为10-2 m/s、10-5 m/s及10-7 m/s，α = 1，荷载频率 f = 1，埋深 h = 1，其余参数同上，计算 z1 = 0.5处土体的切向位移，获得如图10所示地基土切向位移与无量纲半径的变化曲线.

.

5   结   论

假定土体剪切模量随深度呈指数函数非线性分布，基于Biot固结理论与动力学原理，建立竖向非均质饱和地基内部作用简谐扭转荷载时的动力微分方程，通过Hankel变换与逆变换求得土体的剪应力与切向位移解答，据此基于Mathematica编制出相应计算程序，通过参数分析获得土体参数对地基扭转动力响应的影响规律，主要结论如下：

1）埋置简谐扭转荷载作用下，土体的剪应力、切向位移均沿径向呈现出明显波动规律，波动频率随荷载频率的增大而增大，切向位移幅值随荷载频率的增大先增后减，即存在某一临界荷载频率，使得切向位移幅值最大.

2）对比分析表明，随着土体剪切模量指数分布函数中非均勻系数α增大，土体的剪应力、切向位移逐渐减小，系数α每增大0.5个单位，剪应力峰值下降约23%，切向位移峰值下降约40%，说明相同荷载条件下，剪切模量随深度呈指数函数分布的非均质地基比均质地基变形小.

3）埋置简谐扭转荷载作用下，土体的切向位移自地面开始沿深度表现为先增加后减小的变化趋势，荷载作用面处，土的切向位移最大，剪应力发生突变，荷载主要影响范围为其作用面上、下各两倍荷载作用半径（2a），且随荷载埋深h的增大，土体最大切向位移逐渐减小，h = a时最大切向位移值下降75%，h = 2a时最大切向位移值下降约90%，h > 4a时最大切向位移值近似为零.

4）地基土动力渗透系数kd的变化对土体应力以及切向位移几乎没有影响.

诚然，上述解答暂仅针对竖向非均质（土体剪切模量随深度呈指数函数非线性分布）饱和地基中埋置简谐扭转荷载的动力响应问题，对于地基为非饱和、双向非均质或分层等更为复杂的情况，有待进一步的深入研究.

参考文献

[1]    LAMB H. On the propagation of tremors over the surface of an elastic solid[J]. Philosophical Transactions of the Royal Society of London，1904，203：1—42.

[2]    CHAO C C. Dynamical response of an elastic half-space to tangential surface loadings[J]. Journal of Applied Mechanics，1960，27（3）：559—567.

[3]    王贻荪.半无限体表面在竖向集中谐和力作用下表面竖向位移的精确解[J]. 力学学报，1980，16（4）：63—68.

WANG Y S. Exact solution of surface vertical displacement of half-space under vertical concentrated harmonic force[J]. Chinese Journal of Theoretical and Applied Mechanics，1980，16（4）：63—68. （In Chinese）

[4]    REISSNER E，SAGOCI H F. Forced torsional oscillations of elastic half-space. I[J]. Journal of Applied Physics，1944，15（9）：655—662.

[5]    RAHMAN M. The Reissner-Sagoci problem for a half-space under buried torsional forces[J]. International Journal of Solids & Structures，2000，37（8）：1119—1132.

[6]    BIOT M A. Theory of propagation of elastic waves in a fluid-saturated porous solid. I. Low-frequency range[J]. Journal of the Acoustical Society of America，1956，28（2）：179—191.

[7]   PHILIPPACOPOULOS A J. Buried point source in a poroelastic half-space[J]. Journal of Engineering Mechanics，1997，123（8）：860—869.

[8]    陈胜立，张建民，陈龙珠. 饱和地基中埋置点源荷载的动力Green函数[J]. 岩土工程学报，2001，23（4）：423—426.

CHEN S L，ZHANG J M，CHEN L Z.  Dynamic Green′s functions of saturated soils subjected to the internal excitation[J]. Chinese Journal of Geotechnical Engineering，2001，23（4）：423—426.（In Chinese）

[9]   CHEN S L，CEN L Z，PAN E. Three-dimensional time-harmonic green′s functions of saturated soil under buried loading[J]. Soil Dynamics & Earthquake Engineering，2007，27（5）：448—462.

[10]  张智卿，王奎华，谢康和，等. 非饱和土中弹性支承桩的扭转振动响应分析[J]. 计算力学学报，2008，25（4）：552—556.

ZHANG Z Q，WANG K H，XIE K H，et al. Torsional vibration of a pile with elastic bottom boundaries embedded in unsaturated soil[J]. Chinese Journal of Computational Mechanics，2008，25（4）：552—556. （In Chinese）

[11]  陳刚，徐长节，吴大志，等. 半空间饱和土在埋置扭转振动荷载作用下的响应[J]. 浙江大学学报（工学版），2005，39（10）：1566—1570.

CHEN G，XU C J，WU D Z，et al. Dynamic response of buried torsional oscillations in saturated half-space[J]. Journal of Zhejiang University （Engineering Science），2005，39（10）：1566—1570. （In Chinese）

[12]  NAGGAR M H E. Vertical and torsional soil reactions for radially inhomogeneous soil layer[J]. Structural Engineering & Mechanics，2000，10（4）：299—312.

[13]  ZHANG Z，PAN E. Dynamic torsional response of an elastic pile in a radially inhomogeneous soil[J]. Soil Dynamics & Earthquake Engineering，2017，99：35—43.

[14]  AWOJOBI A O，GIBSON R E. Plane strain and axially symmetric problems of a linearly nonhomogeneous elastic half-space[J]. Quarterly Journal of Mechanics & Applied Mathematics，1973，26（3）：285—302.

[15]  鄒新军，徐洞斌，王亚雄. 成层地基中单桩受扭弹塑性分析[J]. 湖南大学学报（自然科学版），2014，41（9）：72—78.

ZOU X J，XU D B，WANG Y X. Elastic-plastic analysis of torsional single pile in layered subsoil[J]. Journal of Hunan University （Natural Sciences），2014，41（9）：72—78. （In Chinese）

[16]  王国才，丁翠红，王哲，等. 层状地基中单桩扭转振动特性分析[J]. 计算力学学报，2009，26（5）：715—721.

WANG G C，DING C H，WANG Z，et al. Study on the torsional oscillations of single pile embedded in layered foundation[J]. Chinese Journal of Computational Mechanics，2009，26（5）：715—721. （In Chinese）

[17]  WROTH C P，RANDOLPH M F，HOULSBY G T，et al. A review of the engineering properties of soils with particular reference to shear modulus[R]. London：CUED/D-Soils TR，Cambridge University，1979.

[18]  RAJAPAKSE R K N D. A torsion load transfer problem for a class of non-homogeneous elastic solids[J]. International Journal of Solids and Structures，1988，24（2）：139—151.

[19]  邹新军，徐洞斌，王亚雄，等. 非均质地基中单桩受扭弹塑性分析[J]. 土木工程学报，2015，48（11）：103—110.

ZOU X J，XU D B，WANG Y X. Torsional elasto-plastic analysis of single piles in heterogeneous ground[J]. China Civil Engineering Journal，2015，48（11）：103—110. （In Chinese）