岩体粘弹塑性分析中的惯性力问题.pdfVIP

华中理I大学出版社 黪嚣_ 丛篓耋兰虹： !兰星盟!二!!!!二；?王!里 605 岩体粘弹塑性分析中的惯性力问题 王芝银牵云鹏 西安矿业荦院710054 摘要 本文探讨了在粘弹童性变形条件下，惯性力对总体基本方程的影响，建立了惯性力矩 阵、速度矩阵及位移矩阵的计算格式，给出了粘弹塑性变形条件下考虑和不考虑惯性力的比较． 关t词 粘弹童性．岩体．惯性力 1．引 言 流变特性是岩体的一种普遍存在的属性．在进行岩体工程地质力学行为的仿真模拟时，考 虑流变特性的常规作法有三种，即以现场观测或室内试验为基础的经验分析法，利用介质模型 建立微分本构模型的解析或数值分析法及以蠕变柔量或松弛模量为核函数的积分模式求解 法。无论是那种方法，无论是计算分析还是现场实测，通常所得到的位移分量随时间的变化规 律与应变随时间的变化规律有完全相似之处。如图1为某矿山巷道实铡位移与时间的变化曲 线LlJ，该曲线是巷道开挖后紧跟工作面作锚喷支护，然后及时设置测点进行观测得到的。图2 是由一般蠕变试验所获得的位移一时问曲线。 o f／d 图1某矿山巷道量测曲线 圈2试验位移时间曲线 由图1、图2可以看出，在变形开始阶段(oA段)，位移增加较快，位移速率迅速减小，则变 形体在此段具有负加速度；在稳定变形阶段(AB段)，位移速率趋于常数，加速度变为零值；图 1中当位移速率为零值时，工程岩体(围岩)是稳定的，流变变形属稳定变形，而对于图2所示 情况，经过位移速率的常数值之后，速率叉重新开始变化，这时加速度从零值向正值发展，且越 来越大，这种情况属不稳定变形。由以上分析可见，位移随时间变化的更一般的规律是负加速 度阶段、零加速度阶段和正加建度阶段。在零加速度阶段，当位移速率不为零值时，适时对工程 岩体进行加固，可以改变位移速率的变化趋势，强迫变形不出现正加速度(如图1)，使工程处 于稳定状态I而在零加速度阶段，位移速率不为零且没有适时采取加固措施，则必导致变形单 元体破坏，随着时间的持续，岩体从原小变形过渡到大变形，且破坏继续发展。这种破坏的范围 不断扩大和贯通，最终导致工程整体失稳。因此，在进行粘弹塑性分析时，位移加速度的存在是 ·368· 必然的。而在变形较大情况下，由于加速度所引起的惯性力对计算结果的影响是不可忽视的一 个重要方面。为此，本文对在粘弹塑性变形条件下，惯性力对总体基本方程的影响作了详细分 析探讨，推导了惯性力矩阵、速度矩阵和位移矩阵的计算格式，对粘弹塑性变形条件下考虑和 不考虑惯性力问题作了分析比较。 2．惯性力矩阵，速度矩阵 在粘弹塑性岩体的变形过程中，对于二维问题，单元e任惹一时刻t时的位移分量为 糍矧=[Ⅳ(剐)](㈨} (1) 式中，{u。(￡))为单元节点位移矩阵；[Nb，y)]为形函数矩阵。 设单元e内每单位面积质量为M，则每单位面积惯性力为 髓嬲卜M鬏搿；)=一M№神，竽 ㈣ 式中，正o，y，￡)，Lo，y，t)分别为单元c内任意点沿z方向和y方向的惯性力。 假设给定一时间步长为，时步开始时单元节点位移为{Ul(t)}，时步终了时，单元节点位 移为{U。O+血))，则由泰勒展开式，并忽略二阶以上微量有 (3) {以。+血))={阢(。}+出堕％≯+等掣 则驾掣一岳㈨(训一岳肌引 (4) 式中 (△仉(血)}为血时步内单元节点的位移增量，{y。(f)}为时步开始时单元节点的速 度，而在时步终了时单元节点的速度为 m(什训一限∽)+血掣