[发明专利]一种基于能量耗散原理的挡土墙地震土压力时程计算方法

说明书

本发明提供了一种基于能量耗散原理的挡土墙地震土压力时程计算方法：建立计算模型；分别计算得到滑动土楔体的地震惯性力、重力、墙对滑动土楔体的支挡力；基于能量耗散原理分别计算得到滑动土楔体的外力所做的功率以及内能耗散；基于能量守恒原理建立平衡方程；计算得到地震持续作用下对应的挡土墙对滑动土楔体的支挡力的最大值时程，即为所求的地震土压力时程；绘制得到地震土压力时程曲线。本发明解决了传统地震土压力计算将滑裂面假定为斜直线，将地震作用简化为一个大小恒定、方向不变的惯性力，未考虑地震波在坡体中传播时对地震土压力影响，以及所得到的地震土压力大多为一个定值而失真，不能得出地震土压力如何随时间变化等一系列的问题。

技术领域

本发明属于地震土压力分析技术领域，尤其涉及一种基于能量耗散原理的挡土墙地震土压力时程计算方法。

背景技术

重力式挡土墙由于其形式简单、就地取材容易、施工方便等优点，已经成为目前最常用的一种挡土墙形式，其特点是以挡土墙自身重力来维持挡土墙在土压力作用下的稳定性，但目前关于该结构的抗震设计过于简单，对于设计的安全性和经济性缺乏切实有力的证据，缺乏关于重力式挡土墙在地震持续作用下土压力如何随时间变化的具体研究，因此需要对重力式挡土墙的地震土压力时程进行研究。

现有的关于地震条件下重力式挡土墙所受地震土压力的计算已有一些方法，例如：拟静力法、Newmark方法等。但已有的方法将地震荷载简化为一个大小恒定、方向不变的惯性力，未考虑地震波在坡体中传播时对地震土压力的影响，并且大多都假定滑裂面为直线以及所得到的地震土压力大多为一个定值而不能得出地震土压力随时间所发生的变化，这与实际工程不符。关于如何基于能量耗散原理通过功率的形式计算地震条件下曲线滑裂面重力式挡土墙地震土压力时程目前仍属空白。因此，亟待一种简便实用的基于能量耗散原理的重力式挡土墙地震土压力时程计算方法出现，从而更加精确可靠地进行重力式挡土墙的抗震设计，兼顾设计的安全性和经济性，对重力式挡土墙抗震减灾领域的发展具有重要意义。

发明内容

针对现有技术中的上述不足，本发明提供的一种基于能量耗散原理的挡土墙地震土压力时程计算方法解决了传统地震土压力计算将滑裂面假定为斜线，将地震作用简化为一个大小恒定、方向不变的惯性力，未考虑地震波在坡体中传播时对地震土压力的影响，以及所得到的地震土压力大多为一个定值而不能得出地震土压力如何随时间变化等一系列的问题。

为了达到以上目的，本发明采用的技术方案为：

本方案提供一种基于能量耗散原理的挡土墙地震土压力时程计算方法，包括如下步骤：

S1、根据挡土墙后土体曲线滑裂面建立计算模型；

S2、利用所述计算模型计算得到滑动土楔体的地震惯性力、重力以及挡土墙对滑动土楔体的支挡力P_a；

S3、根据所述地震惯性力、重力以及挡土墙对滑动土楔体的支挡力P_a，基于能量耗散原理计算得到滑动土楔体的外力所做功率，并基于能量耗散原理利用土体的黏聚力和内摩擦角计算得到滑动土楔体的内能耗散；

S4、根据所述滑动土楔体的外力所作功率和内能耗散建立平衡方程；

S5、根据所述平衡方程计算得到地震土压力时程；

S6、根据所述地震土压力时程绘制时程曲线，完成对重力式挡土墙地震土压力的计算。

进一步地，所述步骤S1中计算模型中摆线的切线斜率tanω的表达式如下：

其中，表示y对x的求导，x和y均表示曲线滑裂面的方程，θ表示生成摆线的圆的半径所经过的弧度，R表示土体作用在滑裂面上的反力。

再进一步地，所述步骤S2包括如下步骤：