基于ANSYSLS―DYNA矿区地下巷道三维动力响应分析.docVIP

基于ANSYSLS―DYNA矿区地下巷道三维动力响应分析 　　摘要：为研究煤矿地下巷道在地震作用下的动力响应，考虑巷道结构的损伤效应，基于显式算法的土结构体系动力学方程，建立ANSYS/LS-DYNA三维数值模型，并在模型指定位置设置观测点。通过观测点的动力响应特征研究巷道的动态响应。结果表明：巷道在地震作用下，半拱形巷道应力主要分布在拱帮和墙角立板、底板处，且拱帮和墙角处塑性体积应变明显，巷道墙角处的损伤明显大于其他各部，是巷道容易发生破坏的部位。 　　关键词：地下巷道；动力响应；有限元分析；地震作用；ANSYS/LS-DYNA 　　中图分类号：TD32 文献标识码：A 文章编号：1000-0666（2016）01-0022-06 　　0 引言 　　近年来，随着我国煤炭资源开采的增多，我国煤矿区形成了大量形式各异、大小功能不同、纵横交错的地下巷道。我国80%的矿区处于地震多发区，由于缺少完善的地下结构抗震设计规范，巷道在地震作用下的稳定性值得我们关注，因此研究巷道的地震动力响应特征及其变化规律十分必要。孙超等（2009）已经对地下结构的动力稳定性问题做了大量研究：刘晶波等（2005）对地铁盾构隧道的地震响应进行了数值分析；刘书贤等（2013）研究了煤矿采动与地震耦合作用下建筑物灾变机理；章泳健等（2010）基于LS-DY-NA研究了隧道的初应力计算方法；姜耀东等（2005）分析了放炮震动引发煤矿巷道动力失稳的机理；高富强等（2009）分析了动力扰动下巷道围岩力学响应。 　　对于地下结构的研究比较多的是对地铁等浅埋结构的研究，但对于地震作用下巷道响应的研究较少，且考虑因素单一，对于巷道结构在地震中的损伤（赵宝友，2009）、岩层分布及其与巷道的接触和波动效应未给予充分的考虑。本文针对煤矿开采地下巷道在初始应力场和地震共同作用下巷道的地震响应研究的不足，基于ANSYS/LS-DYNA的显式动力有限元模型，考虑巷道在地震过程中的损伤及巷道与围岩的接触效应，在模型边界设置无反射边界，对煤矿采煤地下巷道进行了三维数值地震动响应分析研究。 　　1 巷道围岩-结构体系动力有限元分析 　　由于巷道结构与围岩之间存在相互作用，在分析巷道围岩体系结构的同时还要对结构围岩进行分析。围岩作为一种半无限介质，在数值计算中一方面主要是通过引入人工边界（如粘性边界、粘弹性边界等）模拟连续介质的辐射阻尼，保证散射波能量从有限区域穿过人工边界（何伟，2011）。另一方面是模拟围岩的非线性，本文考虑围岩材料非线性和接触非线性，围岩采用D-P理想弹塑性模型，就可以很好地描述围岩的非线性行为。 　　1.1 巷道围岩-结构体系动力方程 　　根据达朗贝尔原理可以建立该体系地震动作用下任意时刻有限元平衡方程（克拉夫，彭津，2006）： 　　Mü+Cu+Ku=P. （1）式中，M为质量矩阵，C为阻尼矩阵，K为刚度矩阵，P是节点荷载向量；u、u、ü分别是体系的位移、速度和加速度。该运动方程一般将整个地震过程分成若干微小时间段h，假设体系在微小时间段内体系是线性变化的，然后采用直接积分法解得，LS-DYNA采用显式中心差分法，即： 　　阻尼采用Rayleigh阻尼，它是将结构整体阻尼矩阵C表示成整体质量矩阵M和整体刚度矩阵K的线性组合：C=αM+βK，式中比例常数α和β阻尼系数可以由振型阻尼比计算得到，在LS-DYNA中通过关键字：*DAMPING设置（ISCT，2007）。 　　2 地震作用下巷道-围岩动力失稳演化分析 　　2.1 有限元分析模型 　　建立三维有限元模型，模型尺寸为y向610m，X向660m，Z向300m；岩层和巷道物理力学性能如表1，岩层分布如表2。巷道为半拱形巷道，拱直径4m，巷道的纵向长度为200m，衬砌厚度为400mm。围岩和巷道衬砌使用solid164体单元，人工边界采用无反射边界。岩层材料采用理想弹塑性DP材料，巷道结构采用C50混凝土，LS-DYNA中通过*MAT_072113材料模型实现巷道损伤本构。 　　2.2 地震波的选取与调整 　　国内外研究表明，选取地震波需满足强度、频谱特性、波的持时的要求，不同地震波参数得出的地震响应可能相差明显，因而选择合理的地震波参数进行动力分析是必要的。本文选用了15s的El Centro地震波，最大加速度调整为150cm/s2，如图2所示，地震波在基岩底部沿X水平输入。 　　2.3 数值模拟结果分析 　　在巷道100m处截取截面，并在截面不同部位设置观测点，如图3所示。 　　2.3.1 巷道应力分析 　　图4为巷道应力的典型分布图，从图中可知，巷道的等效应力和S1主要分布在墙角和拱帮偏下，S3主要分布两侧拱腰处；最大剪应力在墙角处最大，拱帮偏下处次之；这说明巷道的