据国有重点煤矿的统计分析,全国“三下”(建筑物、水体、公路铁路)下压煤约137.9亿吨,其中建筑物下压煤就高达87.6亿吨(王金庄,郭增长,2002)。地下煤炭的开采在地下岩土层形成了一个无实物填充承载力的空间,破坏了岩土层的完整性,导致煤矿采动区上覆岩体发生变形、冒落、断裂、弯曲等一系列破坏,严重威胁了矿区的建筑安全。

目前国内专家学者已经对煤矿采动区建筑物保护问题开展了研究工作:于广明等(2004)利用有限元软件研究了地震作用下煤矿采动区建筑物的动力灾变问题; 滕永海和黄乐亭(1997)根据现有的实际测量资料库,指出了抗变形建筑物与其对应地表的下沉值在宏观上是基本一致的; 谭志祥和邓喀中(2006)为受采动影响的建筑物保护设计提供了科学依据,总结出了不同长高比条件下建筑物变形与地表变形之间的综合关系; 刘书贤等(2010,2013)提出了煤矿采动区建筑物抗震抗变形双重保护支座,对开采变形和地震耦合作用下的建筑物灾变机理进行了探讨。以上研究没有系统地从建筑物损伤的角度,深入探讨在开采变形和地震联合作用下建筑物的损伤演化灾变过程,缺乏对建筑物在开采变形和地震联合作用下损伤积累的研究(刘书贤,2014a,b,c)。

现有的研究对于煤矿开采变形与地震二者之间成灾机理的研究,更多的是应用“耦合”一词,由于煤矿采动区工程地质环境的复杂性,煤矿采动荷载与地震荷载对建筑物所产生的损伤并非简单的叠加,并且二者之间的成灾机理也不相同,所以“耦合”一词有失偏颇。笔者针对地下煤炭开采对建筑物抗震性能的扰动研究尚不完善的现状,基于损伤力学理论,建立煤矿采动损害影响下的建筑物动力损伤判据,通过建筑结构内力变化来探讨煤矿采动与地震联合作用下建筑物损伤演化灾变过程。

煤矿采空区的地表移动变形破坏了建筑物与地基间的初始平衡状态,在建筑物中产生了附加应力,导致了建筑物变形破坏。

定义张量(I-D)为由于损伤而造成的面元dA^*的缩减(李兆霞,2002),即

dA^*=(I-D)dA.(1)

式中:I为瞬时损伤; D为损伤张量。

煤矿采动对建筑物的损伤破坏表现在由于建筑自身重力的作用下的徐变损伤过程,而徐变过程中的总损伤包括瞬时损伤和徐变损伤,因此

D^·(σ,t)=D^·(σ,0)+D^·_cr(σ,t),(2)

式中:D^·(σ,t)为瞬时损伤; D^·_cr(σ,t)为徐变损伤; 即

D^·_cr(σ,t)=D_b(σ,t)+D_M(σ,t).(3)

一般的D^·_cr(σ,t)=F(D_cr,σ,t)。

当应力水平低于长期强度时,σ/f_c'<0.8徐变最终处于稳定阶段,对于给定的常应力,徐变速率和损伤值都趋于定值。这时,D_B=0,D_cr=D_B。

损伤变化率关于损伤是衰减的,并可认为是损伤和应力的显函数,因此D^·_cr=-γD_cr+f(σ).

边界条件为

D^·_cr(σ,0)=0,

D^·_cr(σ,t)〖JB<1*|〗_s=D_m(σ)-D(σ,0).(4)

D_m(σ)是在常应力下界面裂缝损伤D_m的最大值。

D^·_cr(σ,t)〖JB<1*|〗_s=[D_m(σ)-D(σ,0)](1-e^-rt).(5)

假设:混凝土是均质、弹性和各向同性的(李杰,杨卫忠,2009)并且徐变损伤是各向同性的,可选取标量D作为损伤的度量; 对于和时间相关的应力,有效应力相应地为

σ^～(t)=(σ(t))/(1-D(t)).(6)

将Lemaitre提出的应变等效假设推广到和时间相关的徐变问题,推导出弹性徐变耦合损伤的徐变损伤本构方程:

ε(t)=(σ^～(t))/(E(t))-∫_τ1^t()/(t)[1/(E(t))+c(t,τ)]σ^～(τ)dτ.(7)

由此得采动区建筑物的损伤本构关系为

ε(t)=(σ(t))/(E(t)[1-D(t)])-∫_τ1^t()/(t)[1/(E(t))+c(t,τ)]

(σ(t))/(1-D(t))dτ.(8)

式中:ε(t)为钢筋混凝土的应力; σ(t)为钢筋混凝土的应变; E(t)为弹性模量; c(t,τ)为龄期τ下的徐变度; [1-D(t)]=σ/(σ^～), σ^～为有效应力。

由公式(3)可知:煤矿采动损害影响下建筑物的次生损伤与其应力成反比(克拉夫,彭津,2006)。

耦合性即块间联系,是一种量度,泛指系统结构中各模块间相互联系紧密程度; 而联合则是一种新的数据类型,它是一种特殊形式的变量。由于煤矿采动区工程地质环境的复杂性,煤矿采动荷载与地震荷载对建筑物所产生的损伤并非简单的叠加,根据混凝土损伤本构方程,由于二者之间的成灾机理也不相同,所以“耦合”一词有失偏颇,“联合”一词更加准确。

煤矿采动区的开采变形是随着开采过程逐渐发展的,是长期缓慢发展的过程,不仅发生在煤炭开采过程中,而且在开采完成后仍然不断发展,对煤矿采动区的土地及建筑物损害严重; 煤矿采动区的地震则是一个瞬间剧烈破坏的过程,在短时间内能量急剧释放,地表的移动变形强烈,当地震结束时地表移动变形几乎完全消失,对建筑物的损伤破坏也仅存在于地震发生的过程中。

通过有限元数值计算,所得到的不同工况下建筑物的地震动力响应如图3所示。工况一为不考虑煤矿采动作用仅考虑地震荷载,工况二为考虑煤矿采动与地震联合作用。

分析图3a可知:工况1建筑物层间位移变化曲线比较平滑,顶层最大位移为0.045 m; 工况2建筑物层间位移变化曲线比较层间位移变化急剧,顶层最大值达到0.075 m,是工况1的1.67倍。由此可以判断:煤矿采动作用对建筑物所产生的次生损伤对其地震动力响应影响较大,地震发生时,建筑物的次生损伤进一步发展演化。地震作用下煤矿采动作用加剧了结构的动力反应,严重降低了建筑物的抗震安全性能。

分析图3b可知:工况1建筑物底层的最大加速度为1.225 m·s^-2,顶层的最大加速度为7.163 m·s^-2; 工况2建筑物的底层最大加速度为4.788 m·s^-2,顶层的顶层为17.163 m·s^-2; 分别是工况1底层的加速度3.91倍,顶层的加速度2.5倍。对于煤矿采动与地震联合作用下建筑的最大加速度曲线在五层六层发生突变现象(加速度变化较大),说明在煤矿采动损害影响下建筑物产生损伤,其强度和刚度都有一定程度的退化,地震发生时建筑物鞭梢效应明显,导致其动力响应较大。

分析图3c可知:不考虑煤矿采动对建筑物所产生的次生损伤,地震作用下建筑物的最大层间剪力曲线变化较为光滑均匀; 而煤矿采动与地震联合作用下建筑物的最大层间剪力曲线变化趋势非常陡峭,而且其转折点非常明显(转折点主要集中在楼层的第二层和第五层),说明这两层为建筑物的结构薄弱处,煤矿采动作用对建筑物的抗震能力影响较大,需要根据煤矿采动对建筑物所产生的次生损伤予以控制,开展建筑物的煤矿采动损伤评估以保证其安全性。

分析图3d可知:在两种不同荷载工况影响下的建筑物最大层间轴力变化曲线均较为光滑,最大轴力在底层达到最大(与最大剪力相似),说明地震作用下建筑物的底层位置(与建筑物是否发生次生损伤无关)需要予以重视,同时由于地震作用力是按照建筑物的刚度进行分配的,建筑物底层一旦发生损伤破坏,建筑物的质量中心和刚度中心不一致会导致建筑物产生扭转效应,加剧

建筑物的地震动力响应,加剧建筑的倒塌破坏。

分析图3e可知:随着建筑物高度的增加,建筑物的弯矩峰值整体上呈现出减小的趋势,但是在建筑物的第二层却达到最大。结合建筑的层间剪力峰值、层间轴力峰值、位移峰值和加速度峰值可知:在建筑物的第二层处为该结构的特殊位置(薄弱层位置)。在煤矿采动损害影响下,薄弱

图3 不同工况下建筑物的地震动力响应(a)、最大加速度(b)、最大层间剪力(c)、最大层间轴力(d)、最大层间弯矩(e)
Fig.3 The seismic dynamic response(a),the biggest acceleration(b),the biggest interlaminar stress(c),the biggest intertaminar axial force(d),the biggest interlaminar bending

层位置极易出现大量的塑性铰,尤其是在建筑物底部梁柱节点处形成塑性铰,构件强度减小、刚度退化,严重削弱了建筑物的抗震性能,在煤矿采动建筑的损伤安全评估工作中必须对结构的薄弱层位置予以重视,并采取合理的抗震加固措施。

本文从损伤角度出发,基于损伤结构理论,利用有限元计算软件ANSYS建立了煤炭开采区建筑物在地震作用和采动地震联合作用下的计算模型,得到各工况下的最大加速度,层间剪力轴力等结果,分析得到如下研究成果:

(1)由于煤矿采动区工程地质环境的复杂性,煤矿采动荷载与地震荷载对建筑物所产生的损伤并非简单的叠加,并且二者之间的成灾机理也不相同,所以“耦合”一词有失偏颇,“联合作用”能够更好的描述煤矿采动与地震对建筑的破坏。

(2)基于损伤因子所建立的煤矿采动区建筑物的损伤本构关系为

ε(t)=(σ(t))/(E(t)[1-D(t)])-∫^t∫_τ1/(t)[1/(E(t))+c(t,τ)](σ(t))/(1-D(t))dτ.

该本构关系较好的反应了煤矿采动损害影响下建筑物的次生损伤与其应力成反比,揭示了煤矿采动作用降低建筑物抗震性能的本质原因。

(3)地震作用下煤矿采动对建筑物所产生的次生损伤不断发展演化,导致建筑物的薄弱层位置形成塑性铰,严重削弱了建筑物的抗震性能,同时煤矿采动也加剧了建筑的地震动力响应,对与煤矿采动建筑必须对煤矿开采变形与地震荷载的成灾机理,以及二者的联合作用成灾机制予以重视,开展煤矿采动建筑物抗开采变形隔震保护体系相关的研究工作。