1、煤矿采动与地震联合作用下建筑物损伤演化致灾分析摘要：针对煤矿采动区建筑物的地震损伤演化过程尚不明晰的问题，考虑到煤矿采动与地震联合成灾机制的复杂性，基于损伤力学理论，采用地震工程学和开采学理论，建立煤矿采动损害影响下的建筑物动力损伤判据，探讨煤矿采动与地震联合作用下建筑物损伤演化灾变过程，通过有限数值计算分析煤矿采动建筑物地震动力破坏的损伤成灾机制。结果表明：“联合作用”能够更好的描述煤矿采动与地震对建筑的破坏作用，煤矿采动损害影响下建筑物的次生损伤与其应力成反比，煤矿采动加剧了建筑的地震动力响应，所产生的次生损伤不断发展演化，导致建筑物的薄弱层位置形成塑性铰，严重削弱了建筑物的抗震性能。 关2、键词：成灾机制；开采；地震作用；损伤演化；采动损伤中图分类号：TV441 文献标识码：A 文章编号：1000-0666（2015）01-0155-060 引言据国有重点煤矿的统计分析，全国“三下”（建筑物、水体、公路铁路）下压煤约137.9亿吨，其中建筑物下压煤就高达87.6亿吨（工金庄，郭增长，2002）。地下煤炭的开采在地下岩土层形成了一个无实物填充承载力的空间，破坏了岩土层的完整性，导致煤矿采动区上覆岩体发生变形、冒落、断裂、弯曲等一系列破坏，严重威胁了矿区的建筑安全。目前国内专家学者已经对煤矿采动区建筑物保护问题开展了研究工作：于广明等（ 2004）利用有限元软件研究了地震作用下煤矿采3、动区建筑物的动力灾变问题；滕永海和黄乐亭（ 1997）根据现有的实际测量资料库，指出了抗变形建筑物与其对应地表的下沉值在宏观上是基本一致的；谭志祥和邓喀中（2006）为受采动影响的建筑物保护设计提供了科学依据，总结出了不同长高比条件下建筑物变形与地表变形之间的综合关系；刘书贤等（2010，2013）提出了煤矿采动区建筑物抗震抗变形双重保护支座，对开采变形和地震耦合作用下的建筑物灾变机理进行了探讨。以上研究没有系统地从建筑物损伤的角度，深入探讨在开采变形和地震联合作用下建筑物的损伤演化灾变过程，缺乏对建筑物在开采变形和地震联合作用下损伤积累的研究（刘书贤，2014a，b，c）。现有的研究对于煤矿4、开采变形与地震二者之间成灾机理的研究，更多的是应用“耦合”一词，由于煤矿采动区工程地质环境的复杂性，煤矿采动荷载与地震荷载对建筑物所产生的损伤并非简单的叠加，并且二者之间的成灾机理也不相同，所以“耦合”一词有失偏颇。笔者针对地下煤炭开采对建筑物抗震性能的扰动研究尚不完善的现状，基于损伤力学理论，建立煤矿采动损害影响下的建筑物动力损伤判据，通过建筑结构内力变化来探讨煤矿采动与地震联合作用下建筑物损伤演化灾变过程。1 基于损伤因子的煤矿采动区建筑物的损伤本构关系煤矿采空区的地表移动变形破坏了建筑物与地基间的初始平衡状态，在建筑物中产生了附加应力，导致了建筑物变形破坏。式中：I为瞬时损伤；D为损伤张5、量。煤矿采动对建筑物的损伤破坏表现在由于建筑自身重力的作用下的徐变损伤过程，而徐变过程中的总损伤包括瞬时损伤和徐变损伤，因此由公式（3）可知：煤矿采动损害影响下建筑物的次生损伤与其应力成反比（克拉夫，彭津，2006）。耦合性即块间联系，是一种量度，泛指系统结构中各模块间相互联系紧密程度；而联合则是一种新的数据类型，它是一种特殊形式的变量。由于煤矿采动区工程地质环境的复杂性，煤矿采动荷载与地震荷载对建筑物所产生的损伤并非简单的叠加，根据混凝土损伤本构方程，由于二者之间的成灾机理也不相同，所以“耦合”一词有失偏颇，“联合”一词更加准确。2 地震作用下煤矿采动区建筑物损伤演化的有限元分析模型2.1 6、有限元模型的选取为了探讨地震发生时煤矿采动建筑的次生损伤累积演化过程（杨伟，欧进萍，2009；张春礼，2011），重点分析煤矿采动与地震联合作用下建筑物的损伤演化成灾机制，尤其是建筑物从弹性阶段逐步发展到塑性铰出现的塑性损伤破坏过程，故选取唐山某矿区一框架结构的办公楼（余红发等，2006）。唐山某矿区地质条件良好，抗震设防烈度为度，该现浇钢筋混凝土框架结构办公楼，楼层层高统一为3.3 m，纵向4跨，横向2跨，跨度均为6m，该建筑物总高19.8 m，长24 m，宽12 m。柱子截面尺寸均为600 mm600 mm，梁截面均为300 mm650 mm，楼板厚度为120 mm，梁柱均采用C45混凝土7、，弹性模量E=32.5 GPa，泊松比v=0.2，密度为2450 kg/m?（图1）。该矿区由于地下煤炭开采致使地表发生移动和变形（张永波，2006），根据现场监测站的实际监测数据，由于采煤导致办公楼产生了2 mm/m的倾斜（颜荣贵，1995）。2.2 地震波的选取根据建筑抗震设计规范（GB50011-2010）中的相关内容，要求选取地震波要满足峰值、频谱特性、地震持时的要求，结合该办公楼所处区域条件及结构的自振周期，为了能够分析建筑物在采动与地震联合作用下的耗散能量的变化过程，需要选取较长时间的强烈地震波。本文选用了15 s的唐山地震波，其地震等级为7级，方向为南北向，最大加速度为0. 248、 m.s-2（图2）。3 有限元数值计算结果分析煤矿采动区的开采变形是随着开采过程逐渐发展的，是长期缓慢发展的过程，不仅发生在煤炭开采过程中，而且在开采完成后仍然不断发展，对煤矿采动区的土地及建筑物损害严重；煤矿采动区的地震则是一个瞬间剧烈破坏的过程，在短时间内能量急剧释放，地表的移动变形强烈，当地震结束时地表移动变形几乎完全消失，对建筑物的损伤破坏也仅存在于地震发生的过程中。通过有限元数值计算，所得到的不同工况下建筑物的地震动力响应如图3所示。工况一为不考虑煤矿采动作用仅考虑地震荷载，工况二为考虑煤矿采动与地震联合作用。 分析图3a可知：工况1建筑物层间位移变化曲线比较平滑，顶层最大位移为09、.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；分别是工况l底层的加速度3.91倍，顶层的加速度2.5倍。对于煤矿采动与地震联合作用10、下建筑的最大加速度曲线在五层六层发生突变现象（加速度变化较大），说明在煤矿采动损害影响下建筑物产生损伤，其强度和刚度都有一定程度的退化，地震发生时建筑物鞭梢效应明显，导致其动力响应较大。分析图3c可知：不考虑煤矿采动对建筑物所产生的次生损伤，地震作用下建筑物的最大层间剪力曲线变化较为光滑均匀；而煤矿采动与地震联合作用下建筑物的最大层间剪力曲线变化趋势非常陡峭，而且其转折点非常明显（转折点主要集中在楼层的第二层和第五层），说明这两层为建筑物的结构薄弱处，煤矿采动作用对建筑物的抗震能力影响较大，需要根据煤矿采动对建筑物所产生的次生损伤予以控制，开展建筑物的煤矿采动损伤评估以保证其安全性。分析图3d11、可知：在两种不同荷载工况影响下的建筑物最大层间轴力变化曲线均较为光滑，最大轴力在底层达到最大（与最大剪力相似），说明地震作用下建筑物的底层位置（与建筑物是否发生次生损伤无关）需要予以重视，同时由于地震作用力是按照建筑物的刚度进行分配的，建筑物底层一旦发生损伤破坏，建筑物的质量中心和刚度中心不一致会导致建筑物产生扭转效应，加剧建筑物的地震动力响应，加剧建筑的倒塌破坏。分析图3e可知：随着建筑物高度的增加，建筑物的弯矩峰值整体上呈现出减小的趋势，但是在建筑物的第二层却达到最大。结合建筑的层间剪力峰值、层间轴力峰值、位移峰值和加速度峰值可知：在建筑物的第二层处为该结构的特殊位置（薄弱层位置）。在煤矿12、采动损害影响下，薄弱层位置极易出现大量的塑性铰，尤其是在建筑物底部梁柱节点处形成塑性铰，构件强度减小、刚度退化，严重削弱了建筑物的抗震性能，在煤矿采动建筑的损伤安全评估工作中必须对结构的薄弱层位置予以重视，并采取合理的抗震加固措施。4 结论本文从损伤角度出发，基于损伤结构理论，利用有限元计算软件ANSYS建立了煤炭开采区建筑物在地震作用和采动地震联合作用下的计算模型，得到各工况下的最大加速度，层间剪力轴力等结果，分析得到如下研究成果：（1）由于煤矿采动区工程地质环境的复杂性，煤矿采动荷载与地震荷载对建筑物所产生的损伤并非简单的叠加，并且二者之间的成灾机理也不相同，所以“耦合”一词有失偏颇，“联合作用”能够更好的描述煤矿采动与地震对建筑的破坏。（2）基于损伤因子所建立的煤矿采动区建筑物的损伤本构关系为 该本构关系较好的反应了煤矿采动损害影响下建筑物的次生损伤与其应力成反比，揭示了煤矿采动作用降低建筑物抗震性能的本质原因。（3）地震作用下煤矿采动对建筑物所产生的次生损伤不断发展演化，导致建筑物的薄弱层位置形成塑性铰，严重削弱了建筑物的抗震性能，同时煤矿采动也加剧了建筑的地震动力响应，对与煤矿采动建筑必须对煤矿开采变形与地震荷载的成灾机理，以及二者的联合作用成灾机制予以重视，开展煤矿采动建筑物抗开采变形隔震保护体系相关的研究工作。