水工隧洞岩体脆塑性损伤数值分析方法研究.pdfVIP

现代隧道技术 2006年增刊(140—143) 水工隧洞岩体脆塑性损伤数值分析方法研究 肖 明 田华付颖 闫军南 (水资源与水电工程国家重点实验室，武汉大学水电学院，武汉430072) 摘要对于一些呈脆性特征破坏的水工隧洞，由于岩体本身存在的微裂隙，微空洞等缺陷．以及应变软化的 性质，采用经典的弹塑性理论去求解，已经偏离了工程实际。文章将脆塑性损伤力学应用于地下隧洞开挖分析，并 考虑岩体由于应变软化导致的应力跌落现象，采用该理论较好地模拟了岩体开裂破坏特性，在模拟脆性岩体水工 隧洞开挖时，与工程实际情况更接近，为呈脆性特征破坏的岩体洞室开挖提供了一种有效的数值分析方法。 ‘ 关键词臆塑性 损伤破坏 数值分析 水工鹱洞 生和发展的过程，实际上也是围岩的塑性应变的累 1前言 积过程和能量累积过程。因此当围岩进入损伤后， 对于岩性较硬的岩体。在其应力达到峰值强度 围岩的物理参数和塑性屈服面都随损伤参数的增加 时，往往会有一个应力突发性的跌落过程，这种带有 而减小。为此，我们可将损伤变量D看作是一种塑 脆性破坏特征的现象一旦发生，将会显著影响隧洞 性内变量，进入损伤状态的微分应力增量，则应按下 的围岩稳定，甚至导致结构失效，给地下巷道，地下 式确定： 洞室结构的稳定带来很大的危害”J。怎样用合理 的力学模型来模拟其破坏的过程，是国内外学者共 d一；=筹+曲，箸dD (3) 同关心的问题。本文应用脆塑性损伤力学的方法来 根据损伤的应力分式(1)．对其求偏微分得： 研究岩体缺陷给洞室开挖稳定带来的影响，并给出 了应力跌落的求法，得出了较为符合工程实际的结 果。该方法较好地模拟了应力的跌落过程和损伤破 将式(2)和式(4)代人式(3)，由此得出损伤状 坏的发展过程，为硬岩的损伤破坏迭代提供了一种 态下微分应力增量表达式： 有效的分析方法。 2岩体损伤的基本方程 (5) 对于卸载工况，结构材料符合广义弹性虎克定 Frantziskt和Desai在描述岩石和混凝土结构材 料出现裂隙时‘”，曾提出由于结构受力后产生微裂 律，材料的损伤变量D．是一个常量，则dD=0，因此 损伤状态下的微分应力为： 隙，在微裂隙区应力将被释放和减小，从而形成应力 损伤区，应力损伤可用一个材料损伤系数D来表 (6) 示： 对于加载工况，围岩的损伤破坏是一个不可逆 的破坏过程。随着损伤破坏程度的增加，损伤系数 h岛 (1) qD=(1一D)卿+尹D D也在不断增大。根据围岩的破坏特性，通常采用 式中：口i是符合弹塑性本构关系的应力值，用 围岩的应变率来描述围岩损伤的破坏程度。由于结 微分方程来表示则为： 构受静水压力时不产生破坏，结构损伤破坏主要是 由于塑性应变的偏张量，则三维损伤破坏的演化方 d唧=[％】d8=([D。]-[D，])dei(2) [D。]为弹塑性应力矩阵。 程可以用下式表示： 围岩的损伤破坏