摘要：为了探讨土一结构相互作用对煤矿采动损伤建筑地震动力破坏的影响，基于损伤力学和能量耗散理论，研究采动区土一基础一上部结构的协同作用，通过建立考虑土一结构相互作用的煤矿采动损伤建筑的动力学方程，重点分析了煤矿采动损害影响下的建筑物地震动力灾变演化过程。计算结果表明：土一结构的相互作用对煤矿采动建筑的地震动力响应影响较大，不考虑土一结构相互作用是偏于安全的；煤矿采动作用明显改变了建筑物的结构动力特性，建筑结构的薄弱层位置改变明显、塑性铰和层损伤分布规律发生改变，严重降低了建筑物的抗震性能，对于煤矿采动区建筑应当开展安全损伤评估工作，以保证矿区工程建设的安全性和可靠性。
关键词：土一结构相互作用；损伤力学；抗震性能；煤矿采动；地震动力；安全评估
中图分类号：T0352 文献标识码：A 文章编号：1000-0666（2015）02-0272-08
0 引言
我国拥有丰富的煤炭资源，煤炭资源位居世界首列。对地下的煤炭资源进行开采，采出空间周围的岩层失去支撑而逐渐向采空区内移动，导致上覆岩层原有的应力平衡状态被破坏，采动区的岩层发生应力重分布现象，采动区周围的岩体会发生移动变形、弯曲等破坏现象，导致地表发生沉陷破坏现象，进而对建筑物造成严重的损坏。我国80%以上的矿区位于地震多发带，采动区建筑物除了要承受开采沉陷的影响，还要承受地震灾害的影响（刘书贤等，2013）。因此，开展地震作用下的煤矿采动损伤建筑的抗震性能的研究工作非常有必要。
随着对煤矿采动建筑保护工作研究的不断深入，取得了丰硕的研究成果。周长海（2011）从不均匀沉降对建筑物的危害性冉发，以刚性地基假设条件，建立了层间高层钢筋混凝土框架结构模型，重点研究了沉降与地震耦合作用下结构构件的内力分布和变化规律；夏军武等（2007）建立了地基一独立基础一框架结构共同作用的理论计算模型，揭示了地表变形对建筑物附加变形和附加内力的关系；刘书贤等（2011）在分析地震作用下采动区建筑物抗震能力的基础上，提出了基于建筑隔震技术的的抗开采沉陷隔震保护装置。
以上研究对于地震作用下煤矿采动损伤建筑的研究人多是以刚性地基为假设的，忽略了土一结构相互作用（简称SSI）对建筑物抗震性能的影响。
考虑土结构相互作用后，地基的振动特性将对整个结构体系带来影响，同时土体的滤波效应和耗能改变了地震波的传播途径，对建筑物的动力效应影响较大。王海等（2012）通过研究发现当考虑土一结构的相互作用后，随着土体的剪切波速下降，结构的损伤指数会上升；曹青和张豪（2011）则认为土结构相互作用对风力发电机塔架结构的地震动力响应不容忽视；刘洁平（2009）提出了一种实用的高层建筑土一结构相互作用地震反应整体分析方法，研究了土一结构相互作用对位于厚和深厚地基土层上高层建筑地震反应的影响。
笔者针对煤矿采动和地震共同作用下建筑物损伤演化灾变研究的不足，从能量演化和损伤力学的角度分析煤矿采动损伤建筑的灾变演化过程，综合考虑土一结构相互作用对建筑物抗震性能的影响，深入探讨采动、地震对建筑物的不同致灾机理，重点研究了地震作用下考虑土一结构相互作用的煤矿采动损伤建筑的抗震性能。
1 土一结构相互作用影响下煤矿采动建筑的地震动力学方程
考虑到土体非线性对地震波传递的影响以及采动对土层间力学性能和相互作用的影响，土体本构关系采用Drucker-Prager理想弹塑性模型。地震作用下，建立考虑土一结构相互作用影响下建筑物的动力学计算模型如图1所示。
根据结构动力学理论（克拉夫，彭津，2006），可以建立考虑土结构相互作用影响下的煤矿采动建筑的地震动力学方程
2 煤矿采动建筑的地震动力损伤演化分析
煤矿采动引起的开采沉陷作用导致建筑物局部构件损伤，当地震发生时建筑物的损伤会加剧，由弹性阶段发展为弹塑性损伤阶段。本文采用牛荻涛和任立杰（1996）改进的钢筋混凝土结构的双参数地震损伤模型，来计算建筑结构的楼层损伤值：
2.3 建筑结构极限滞回耗能的计算
建筑结构的极限耗能是指结构在给定荷载幅值或变形幅值下循环至破坏的累积滞同耗能，具体计算如下（杨佑发，崔波，2004）：
3 地震作用下考虑土一结构相互作用的煤矿采动建筑物损伤演化致灾分析
3.1 有限元分析模型
某矿区地质条件良好，地基土体主要有砂土、砂岩、煤岩、泥岩，简化处理后各岩体和煤层力学性能参数如表1所示（刘刚，2011）。
该矿区（抗震设防烈度为7度）有现浇钢筋混凝土框架结构办公楼（图1），该楼共六层，底层层高为4.2 m，其余层高均为3.6 m，纵向4跨，跨度为4.5 m，横向2跨，跨度为6m。该建筑物总高度为4.2 m+5*3.6 m=22.2 m，总宽度为12m，总长度为18 m。柱、梁、楼板、基础均采用C30混凝土，弹性模量E=30 GPa，泊松比v=0.2，密度取2700 kg/m3。基础为筏板基础，尺寸为28 m*16 m*l m，为了更好的模拟土一结构相互作用效应对煤矿采动损伤建筑的影响，地基作用范围选120 m*60 m*11 m。
纵筋选用HRB400，根据《混凝土结构设计规范》（GB50010-2010）巾规定的框架结构最小配筋率和最大配筋率的要求，选取柱、梁板的配筋率如表2所示。
基于有限元分析软件ANSYS内嵌的APDL和FORTRAN语言进行二次开发，编制计算程序对建筑结构的梁、柱以及板巾钢筋进行处理，采用刚度EI等效方法调整弹性模量的整体式建模，整体式的单元刚度矩阵和材料本构矩阵为
经过上述处理后，不需考虑钢筋和混凝土分离的组合式模型，不仅缩短了计算时间，而且提高了计算精度，满足规范要求（何春林，邢静忠，2007）。 模拟地下煤层开采形成采空区的过程中，在土体和建筑物的自重作用下建筑物发生倾斜，其最人沉降值为0.038 m，平均沉降值2 mm/m。采动引起建筑物产生的附加变形和附加内力，局部结构构件产生初始损伤，如图3所示。
通过对比分析图4两种不同结构形式的自振频率可知：刚性地基下建筑结构的前3阶自振频率是考虑土一结构相互作用的1.605倍，由此可知地基土在一定程度上集中了建筑结构有效的地震动力响应，考虑土结构的相瓦作用后结构的自振频率大幅度的减小，即延长了结构的自振周期。由此可以判断：对地震作用下的煤矿采动损伤建筑的分析过程中需要考虑土结构的相互作用，而对地基土体采用刚性地基假设是不完善的。
根据《建筑抗震设计规范》（GB50011-2010）中的相关内容，选取的地震波要满足峰值、频谱特性、地震持时的要求，结合该办公楼所处工程地质条件及结构的自振周期，为了能够清晰的分析建筑物在采动与地震联合作用下的滞同耗能和损伤的变化过程，需要选取较长时间的强烈地震波。本文选用了20 s的El Centro地震波，满足包含最火幅值及持时5 T～10 T（T为建筑物的自振周期）的要求模拟烈度为7度时的罕遇地震作用，方向为南北向，调整后的最大加速度为220 cm・S-2，如图5所示。
3.2 数值模拟结果分析
地下煤炭开采的过程中，地表发生移动、曲率变形、不均匀沉降等现象，建筑物发生倾斜，导致局部构件刚度降低，发生损伤破坏（框架结构梁柱节点出现塑性铰）。地震发生时，煤矿采动损伤建筑在地震动力的作用下加剧了破坏程度（刘书贤等，2013，2014）。本文从损伤力学和能量耗散理论的角度重点考察了土一结构相互作用对煤矿采动损伤建筑的灾变影响。
分析图6可知：不同工况下的建筑物在地震动力作用下，煤矿采动作用对建筑物影响较大：地震发生2 s后，加速度迅速达到峰值，地震单独作用时为4.481m/S2，煤矿采动和地震联合作用时为6.266 m/s2，较地震单独作用时，结构顶层加速度被放大了1.397倍，这对建筑物的抗震是不利的。
分析图7a的楼层最人加速度可知：地震发生过程中底层和顶层的变化幅值最火，主要是地震发生初期短时间内输入建筑结构的能量急剧增大，部分结构构件的能量达到了损伤阀值极限，产生大量的损伤耗散能，致使梁端开裂。随着时间的增加，开裂区向梁柱节点处蔓延，形成塑性铰，降低了建筑物的底层刚度和强度，所以底层加速度出现增人趋势。说明地震作用下煤矿采动作用加剧了结构的动力反应，导致建筑物底部梁柱节点处形成塑性铰，构件强度减小和刚度退化，严重削弱了建筑物的抗震性能。
分析图7b、7c可知：地震发生后，随着层数的增加煤矿采动建筑物的最大楼层位移也增大，与未受采动影响相比，各层位移增加幅度明显，煤矿采动建筑的顶层最大位移较未采动增大了1.845倍。结合图7c可知：煤矿采动损害影响下建筑物的最人层问位移增大（其中底层和顶层的最人层间位移增幅最大），较末采动分别增大了2.23倍和2.67倍。煤矿采动作用对建筑物底层的抗震性能的影响最大，顶层次之，中间层最小，同时也明显改变了建筑物的薄弱层位置，有向下发展的趋势（末采动的建筑的最人层间位移发生在3～4层，煤矿采动后建筑的最大层间位移发生在2～3层），在煤矿采动建筑的损伤安全评估工作中必须予以重视。
由图7d可以看出：不考虑地震作用，煤矿采动引起的开采沉陷促使建筑物底部的剪力增大，增幅13.7 kN；建筑物上部的最大层间剪力变化不明显，此时建构处于弹性阶段。地震发生后建筑物的层间剪力迅速增大，受采动影响的建筑物底部最人层间剪力由143.12 kN增加到179.976 kN，增幅达到36.856 kN，是静力阶段增幅的2.69倍，接近或者超过构件的屈服强度，建筑物的底部损伤演化现象严重，且向薄弱层发展形成多处塑性铰，严重降低了建筑物的抗震性能。
分析图8不同工况下建筑物的楼层损伤值图可知：地震发生后，末受采动影响的建筑的楼层损伤集中在3～4层（这两层的层间位移最大），底部和顶部损伤相对较小，由此可以判断3～4层是建筑物的薄弱层，在结构的抗震设计分析时应予以重视；在煤矿采动损害影响下，建筑物的各楼层损伤值均出现明显增幅，且损伤最大值出现在底部，底层损伤急剧增大了1.643倍，说明建筑物底部已经严重破坏，有可能出现倒塌的危险。
对于煤矿采动区的拟建建筑物，采动和地震联合作用加大了结构的位移反应，在结构设计时应考虑这种不利影响可能产生的P-△效应，在地震中建筑物容易倾覆，因此需要提高建筑物的抗侧移能力，同时加强结构薄弱层的弹塑性变形验算。
4 研究结论
本文基于损伤力学和能量耗散理论，以有限元软件ANSYS为工具，建立煤矿采动与地震联合作用下的有限元模型，重点研究了土一结构相互作用对煤矿采动损伤建筑抗震性能的影响，通过分析层间加速度、层间位移、层间剪力、层间损伤得到如下研究成果：
（1）在煤矿采动建筑的地震动力灾变演化过程中，如果不考虑土一结构的相互作用，而采用刚性地基假设，使得计算结果过于保守、偏于安全。
（2）地震作用下煤矿采动作用加剧了结构的动力反应，导致建筑物底部梁柱节点处形成塑性铰，构件强度减小和刚度退化，严重削弱了建筑物的抗震性能。
（3）地震作用下煤矿采动损伤建筑出现多处塑性铰，且薄弱层位置向下发展，底层损伤相当严重，建筑物由局部构件损伤向结构整体破坏演化发展，增加了震害的影响，建议对煤矿采动建筑进行损伤安全评估时必须予以重视。
