非连续子母块体计算分析模块的开发.docxVIP

非连续子母块体计算分析模块的开发 1 理论和计算分析模块的验证 任何新的数值方法的准确性和有效性都必须通过具体的程序和软件进行验证。因此，项目是数值算法的一个重要环节。 在第Ⅰ部分提出的非连续子母块体理论模型的基础上 选取了3经典算例, 通过对算例的计算及分析, 对本文理论和计算分析模块的正确性及可行性进行了验证。其中, 斜面滑块算例验证了本文提出的非连续子母块体理论模型及计算分析模块在处理连续与非连续接触并存的问题时的正确性与计算精度;两块体相互作用算例验证了将修正的拉格朗日乘子法引入到非连续子母块体理论模型中处理接触问题的可行性与正确性;悬臂梁一端受集中力荷载的算例则验证了非连续子母块体理论在处理连续接触问题的可行性, 同时也表明本文的算法可能用于解决连续大变形问题。最后, 将本文的计算分析模块应用于某石拱桥受极限荷载作用破坏失稳的工程实例问题, 验证了非连续子母块体理论模型处理连续接触问题向非连续接触问题转化 (即岩体破坏问题) 的有效性, 并且将本文计算结果与其他算法的结果进行了对比讨论 2 初始接触模型 斜面滑块问题是非连续变形分析方法中的经典算例之一, 该算例由两个块体组成, 底部块体固定, 当滑动接触面的内摩擦角小于底部块体坡面角时, 上部块体会在重力作用下沿坡面加速下滑。该算例通过计算上部块体滑动的位移和速度来验证算法的精度。Yeung和MacLaughlin曾应用该算例测试其算法的精度, 其主要工作集中在通过改变几何参数和破坏模型来测试算法的静态及动态响应 对于一个受重力作用初始时静止于斜面上块体, 其位移随时间变化的解析方程为 式中:g为重力加速度;α为斜面的倾角;?为块体接触面的内摩擦角。 如图1的几何模型所示, 设底部块体为梯形块体, 底部边长为2 m高1 m斜面的倾角为30°, 块体的底部和右端均为固定约束;上部块体为矩形块体, 长为0.2 m, 宽0.1 m, 初始速度为0;上部块体与底部块体以非连续接触的形式联结在一起。将上部块体和底部块体均人为划分为若干个子块体, 子块体之间以连续接触的形式相互联结在一起。计算参数设置如下:块体密度为2 500 kg/m 从图2不难看出, 滑块随时间变化的位移计算结果与理论值吻合得很好, 误差均在0.01%以内。因此, 可以确定本文提出的非连续子母块体理论模型及其计算分析模块在计算多块体之间连续接触与非连续接触问题上是正确的, 并且具有很高的计算精度。 从图3中各个时步两块体的几何形态可以看出, 滑块在运动的过程中, 块体之间的接触关系是非常复杂的, 而且有连续和非连续接触两种情况, 对于这两种情况的处理是完全不同的。但从本文算例结果图中可以看出, 对于接触的动态寻找与释放是完全正确的。算例中子块体之间的位移协调也得到了很好的实现。 3 块体接触判断 两块体相互作用的算例是一个考虑两个块体在外力作用下相互作用的验证性实例 在两块体相互作用的算例中 (如图4 (a) 所示) , 底部块体的几何尺寸为1.0 m×1.0 m, 块体底部为固定约束;上部块体几何尺寸0.3 m×0.3 m, 左部边界中点受水平方向点荷载, 受力大小为10.0 kN, 顶部中点受垂直向下的点荷载, 受力大小为5.0 kN。将上部块体和底部块体均人为划分为若干个子块体, 块体内部的子块体之间以连续接触的形式相互联结在一起。计算参数:块体密度为2 500 kg/m 如图4可所示, 上部块体在垂直与水平两个方向的力作用下, 沿着底部块体的上表面滑动, 直至从底部块体的上表面滑出。可以看出, 虽然在计算中选取了较小的罚弹簧系数, 但在上部块体与底部块体发生相互作用的过程中, 两个块体没有明显地侵入现象发生。而应用罚弹簧法处理接触时, 如果罚弹簧系数取的过小, 计算中块体之间会出现很大的嵌入, 当嵌入量达到一定程度时, 同时还会造成接触判断的失效, 致使块体间的接触关系出现错误。这说明引入修正的拉格朗日乘子法对块体的接触进行处理是可行的, 同时还可以精确计算出块体之间的接触力。虽然引入修正的拉格朗日乘子法会增加对开-闭迭代的循环次数, 但计算结果证明, 通过少量增加求解的计算量来提高计算精度的做法是可取的。 4 静力问题解析解 对于悬臂梁一端受集中力荷载情况下的变形问题, Belytschko给出了静力问题理论解析解 如图5 (a) 所示, 悬臂梁长L=12 m, 高D=2.0 m, 取单位厚度1.0 m;材料弹性模量E=0.01 MPa, 泊松比μ=0.333 33, 密度ρ=0.68 kg m 5 荷载施加过程中裂缝的变化 图7为某石拱桥在中心偏左的位置受到垂直向下荷载作用下发生失稳的原型试验, Thavalingam应用DDA、离散元和非线性有限元3种方法对该石拱桥进行了分析, 不仅研究了