岩体损伤流变模型及有限元分析.pdfVIP

·828· (工程力学)增刊 1999年 岩体损伤流变模型及其有限元分析 朱珍德 胡定 (四川大学术利工程系．成都610065) 提要 本文运用流变力学与损伤理论．提出了一个损伤演化方程来描述岩体力学性质随时间不断 弱化的现象．从而给出相应的有限元解。 关麓词 流变，扭伤，本掏模型 1前言 考虑岩体流变特性的工程问题粘弹塑性有限元分析已有许多学者做过不少工作，取得 了许多有用的成果【I剖，但是，仍有许多同题值得进一步探讨。譬如，地下工程中巷道开挖 后由于围岩应力重力布，往往存在一个松动圈，该区域微裂纹张开、扩展，岩体力学性质也随 之改变。因此，采用粘弹塑性有限元分析巷道稳定性时，如何考察岩体力学性质随时间不断 恶化这一不利因素，将具有重要的工程意义。 近年来发展起来的损伤理论，为解决这个问题提供了新途径。KachanovL¨M建立丁 第一个蠕变损伤模型¨o；文献”J从微裂纹模型出发得到粘弹性材料的损伤本构方程； Kawamoto¨1提出的二阶损伤张量引用于岩体损伤问题，提出了节理裂隙岩体的损伤张量。 本人基于前人的工作．提出了一个损伤演化方程来描述岩体流变过程中裂隙随时间不 断扩展而导致岩体力学性质不断劣化的现象。最后将损伤变量引入枯弹性本构方程．并对 某地下巷道进行流变一损伤有限元分析，计算结果是可靠的。 2岩体损伤本构方程 岩体裂隙的演化，必将反映到其宏观物性参数和本构关系的变化。由于初始地应力场 与开挖干扰力的作用，岩体裂隙扩展、贯通使其剐度降低粘滞性态改变等劣化现象。因而可 通过考虑损伤对物性参数的影响。建立流变损伤本构方程。 设无损时的本构方程为： E=F(口，E，f，T，…) (1) 桎青年科学基金资助项旨 朱珍静．男，1964．6出生．博士舌．副救授 (工程力学)增刊 J999年 式中：d，e一分别为应力和应变列阵；￡一时间； E=(EJ，E2，Ej。…E。)7一物性参数列阵；T一温度。 损伤产生时，将对物性参数产生影响。物性参数E。(扣J，2。一”)的变化可通过下式 加以考虑： 雷，=似(口，E，￡。T，…) (i=J，2，…H) (2) 称重。为E。的有效形式，对E而言有： E：(E (3J J，E2，…E。) 式(2)中口是包含损伤变量在内的内变量列阵；红则体现了损伤对不同的材料、不同物理参 数的影响，也可由E到重的变化，考察损伤变量的变化。 如将式(3)代入式(1)，可得到一个考虑损伤的本构方程： e=F(口，E，f，T，…) (41 以三体的Kelvin—Voignt模型(见图1)为例， 其无损本构方程为： m !+鲁：=?巨+墨学!】(5) 考虑损伤流变本构关系时，可将本构参数的有效形 二￡3 式代入式(5)得： 1+鲁：=孕巨+墨铲纠㈣ 刚刚咖～浪型 式中：D是仅与泊松比P有关的系数矩阵；西是D的有效形式，是泊松比p有效形式；的矩 7 阵。 3损伤演化方程 文献[7)由热动力学理论出发，对粘弹性本构方程作如下推导： 巷+哦磐=Q。 (7) o