浅埋盾构隧道地基弹簧刚度的求解方法

董正方　王君杰　赵东晓　王文彪

摘要：把浅埋盾构隧道简化为半无限平面的挖孔问题，采用平面弹性理论的复变函数方法，推导了反应位移法中浅埋盾构隧道地基弹簧刚度的解析公式。设定不同的埋深和土体泊松比，计算了隧道周边的压缩和剪切地基弹簧刚度，讨论了隧道埋深和土体泊松比对刚度的影响，以及压缩和剪切地基弹簧刚度的关系。结果表明：浅埋隧道与深埋隧道的地基弹簧刚度存在差异；隧道埋深和土体泊松比对地基弹簧刚度的大小和分布规律影响很大；浅埋隧道的剪切与压缩地基弹簧刚度之比沿隧道周边变化。

关键词：盾构隧道；复变函数；地基弹簧刚度；反应位移法

中图分类号：TU311.4文献标志码：A文章编号：16744764（2013）06002805

20世纪70年代，日本学者提出了盾构隧道抗震设计的反应位移法[1]。该法需要计算隧道周围的地基弹簧刚度，并且地基弹簧刚度值的变化对隧道结构内力有很大影响[23]，然而地基弹簧刚度的求解至今未有满意的成果，因此研究地基弹簧刚度的求解显得十分重要。

地基弹簧刚度的求解有3种方法：理论推导、试验测得、数值计算[4]。其中后2种方法费时费力，而理论推导一般会得到解析或半解析公式，使用简单、快捷，因此大量学者在理论推导方面开展了大量研究[5]。理论推导大都简化为弹性理论的无限平面或半无限平面中挖孔的边界求解问题，但以前的研究要么不是针对反应位移法、要么就是不考虑隧道的埋深[69]。Verruijt[5]提出了利用复变函数理论求解半无限平面的圆孔问题，可以给出特殊边界条件的应力场和位移场，但没有给出地基弹簧刚度的解析公式。本文在Verruijt研究成果的基础上，利用复变函数理论，给出求解浅埋盾构隧道的径向压缩和剪切地基弹簧刚度的解析公式，并讨论二者的关系。1问题的描述

盾构隧道的反应位移法的地基弹簧分为压缩弹簧和剪切弹簧，其刚度的求解一般根据局部变形理论，利用一点上的力除以该力引起的变形得到。在弹性理论中，就是求解隧道周边应力边界条件下的应力场和位移场，然后得到弹簧刚度。深埋隧道的问题可简化为无限平面中的挖孔问题，浅埋隧道则可简化为半无限平面中的挖孔问题。其中，压缩弹簧刚度可通过隧道周边径向应力q边界条件得到，剪切弹簧刚度可通过隧道周边切向应力τ边界条件得到，如图1。

可以看出：深埋时，隧道地基弹簧沿周边均匀分布；当rh值较小时（例如≤0.01），基本上跟深埋情况类似。随着rh的增大，即系数α逐渐增大，隧道埋深逐渐减小，隧道的地基弹簧沿周边分布发生变化；隧道顶部弹簧刚度比深埋时变小，底部弹簧刚度比深埋时变大，这是由于埋深的减小造成隧道孔周变形发生变化造成的。从隧道顶部到底部压缩弹簧刚度先变小后变大，剪切弹簧刚度逐渐变大（除了埋深很浅的情况）。弹簧刚度最大值一般出现在隧道底部，压缩弹簧刚度的最大值可达深埋时的1.5倍，剪切弹簧刚度的最大值可达深埋时的1.2倍；压缩弹簧最小值一般出现在正负60°处，最小可达深埋时的0.1倍，剪切弹簧最小值一般出现在顶部，最小可达深埋时的0.4倍。泊松比越大，地基弹簧刚度最大值越小，说明土质条件越差，地基弹簧刚度越小。

为了讨论某点处剪切弹簧刚度和压缩弹簧刚度的关系，可以把剪切弹簧刚度和压缩弹簧刚度相除，得到的关系如图5。

从图5可看出：深埋时，地基弹簧的剪切刚度和压缩刚度相等。浅埋时，地基弹簧的剪切刚度和压缩刚度不总是相等；靠近隧道顶部和底部的区域，剪切刚度小于压缩刚度，其余区域情况相反。除了埋深较小的情况（例如rh≥0.6），剪切刚度与压缩刚度相差都在50%以内。土体泊松比对剪切与压缩地基弹簧刚度比的影响不大。

5结论

通过保角映射把无限半平面的挖孔问题转变为ζ平面的环域问题，利用复变函数的理论，推导出浅埋隧道的地基弹簧刚度的解析公式，并得出以下结论：

1）浅埋隧道的地基弹簧刚度相当于在深埋隧道地基弹簧刚度上乘以一个系数，沿隧道周边的分布不再均匀。

2）随着隧道埋深逐渐减小，隧道顶部弹簧刚度比深埋时变小，底部弹簧刚度比深埋时变大；泊松比越大，地基弹簧刚度最大值越小。

3）浅埋隧道剪切与压缩地基弹簧刚度不总是相等。

参考文献：

[1]Tateishi A.Underground structures using static finite element method [J].Journal of Structural Engineering/ Earthquake Engineering， JSCE，2005，22（1）：4153.

[2]刘如山，胡少卿，石宏彬.地下结构抗震计算中拟静力法的地震荷载施加方法研究[J].岩土工程学报，2007，29（2）：237242.

[3]刘晶波，王文晖，赵冬冬.地下结构横截面抗震设计分析方法综述[J].施工技术，2010，39（6）：9197.

[4]Sadrekarimi J， Akbarzad M.Comparative study of methods of determination of coefficient of subgrade reaction [J]. Electronic Journal of Geotechnical Engineering，2009，14（E）：114.

[5]Verruijt A.Deformations of an elastic half plane with a circular cavity [J]. Journal of Solids Structures， 1998，35（21）：27952804.

[6]Bobet A.Analytical solution for shallow tunnels in saturated ground [J].Journal of Engineering Mechanics，2001，127（12）：12581266.

[7]Park K H.Elastic solution for tunnelinginduced ground movement in clays [J].International Journal of Geomechanics，2004，4（4）：310318.

[8]刘齐建.软土地铁建筑结构抗震设计计算理论的研究[D].上海：同济大学，2005.

[9]涂忠仁，孙钧，蔡晓鸿.海底隧道围岩抗力系数计算方法研究[J].岩土工程学报，2006，28（8）：10021007.

[10]谷音，刘晶波，杜义欣.三维一致粘弹性人工边界及等效粘弹性边界单元[J].工程力学，2007，24（12）：3137.

[11]小泉·淳.盾构隧道的抗震研究及算例[M].张稳军，袁大军.译.北京：中国建筑工业出版社，2009.

[12]黄茂松，曹杰.隧道地震响应简化分析与动力离心试验验证[J].岩石力学与工程学报，2010，29（2）：271280.

[13]England A H.Complex variable methods in elasticity[M].New York： Dover Publications，2003：2425.

[14]Wang L Z，Li L L，Lyu X J.Complex variable solutions for tunnelinginduced ground movement [J].International Journal of Geomechanics， ASCE，2009，9（2）：6372.

[15]吕爱钟，张路青.地下隧洞力学分析的复变函数方法[M].北京：科学出版社，2007.

（编辑王秀玲）