岩体损伤增量数值流形方法.pdfVIP

_l乜 下 空 间 Vol19 第19卷第5期 No：! ：罂竺!!里 g!墼!竺!!竺!!!竺 ：：————竺些竺 岩体损伤增量数值流形方法+ 王芝银， 剩云鹏， 王思敬2 杨志法2 (1．西安矿业学簏西安7160542．中罔科学院地质研览所北京 100029) ¨膏要】 本文基于盘值浣形方法的基蕾拉话，建立了考虑岩俸扭曲分析的数值流彤方法增量理论噩其平衡方程， 并时边墟存在圈孤滑动面时的变彤、滑葫过程进行了模拟舟析口 I关量词】 岩体攫伤敷值漉彤方法_II量理论 1 引言 在岩体工程中，随着静态或准静忽的过程的不断发展，岩体系统自发形成时空有序的状 态，或达到新的平衡，或发生非平衡非线性现象，这种非平衡导致系统向动态发展。此外。岩体 变形与破坏规律受岩体结构的控制…。不同结构岩体，具有不同的力学行为。在同一系统中，岩 体的结构类型可能有多种，是一种连续与非连续介质的混合体，而且由于工程地质力学的长期 作用，各种不连续面相互贯通，形成可以导致工程岩体整体失稳的潜在界面。对于这种情况，既 有向动态发展的趋势，又是连续介质与非连续介质的耦合体。90年代由石根华发展的数值流 形方法应用有限覆盖技术，吸收有限元与DDA各自的优点，通过在分析域内各物理覆盖上建 立一般覆盖函数和加权求和形成总体位移函数，从而把连续和非连续变形力学问题统一到这 种方法之中，使得其具有解决动静交叉和睦续与非连续介质耦合问题的能力。纵观目前国内外 已发表的有关文献”。’，在数值流形方法中对岩体损伤分析还未涉及，同时，基本理论和方程 均基于全量格式，与实际岩体变形过程和分阶段施工增量加卸载的情况不相一致。本文视岩体 结构体的各种缺陷(如微裂隙、微裂纹和孔隙等)为初始损伤，利用数值流形方法的基本理论． 建立了考虑岩体材料损伤的流形元增量理论，推导出了相应的计算格式，并针对岩土工程实际 对数值流形方法中的分布力固定点处理停出了台理的改进，最后给出了应用算例。 2 位移、应变及应力 在数值流形方法中，总体位移函数可逋过物理覆盖函数来形成，其一般形式为 ，a(x，y)、 l I=【T【。}(x，ly)]{D{e)} (1) I v(x，y)J 式中 【T{。)(x，Y)】=[T(叭】)(x，Y)…R(，)(牝，Y)] (2) 幸中科院工程地质力学开敢研究实验室基金责助项日和煤炭基金(96采10408)资助项目内容 999年第5期 王芝银等：岩体损伤增量数值流形方法+ ’岛0…￡，0， ’Lm(x，，)=ffu。：。：·0 ￡，…0 ￡。1 矗=W。“l(x，Y)(S《(x，Y)i=l，…，q；j=1，…，171 (3) D小)}’ (4) {D【“1=(D。…Dew…+ {Dtml={U．1v，l…… v．。v．。17…… (5) 其中q为单元e上的物理覆盖数；W。¨)为覆盖i上位移函数对单元e的贡献权函数；in／1是物理 11‘o‘I_l-··卜 覆盖i上所定义的位移函数的项数；“Ⅱ，U。为覆盖i上第J个沿x向和Y向的位移变量；Stj(x， Y)为覆盖i上的基本函数。 由式(2)可给出从t时刻到t+At时刻(或某一增量荷载)的位移增量为 ’‰Au(x’(x,vyj：}=n“xty川△D【e)l (6)