饱和多孔介质动力及断裂分析的多尺度有限元和近场动力学方法

发布时间：2020-05-19 06:32

【摘要】：多孔介质材料作为一种由固体颗粒骨架、大量微小孔隙以及孔隙中的液相和气相组成的组合体广泛存在于自然环境和工程应用中。尤其是土体和岩石,作为一种普遍和价廉的材料而被大量应用到房屋、路隧、桥梁和水坝等土木结构中,其中主要涉及到的力学问题有固结、动力响应和断裂破坏等。多孔介质是一种由多相物质共存的组合体,其涉及到的问题又属于多物理场耦合问题。如对饱和多孔介质的力学分析,不仅需要考虑固体颗粒骨架的变形对孔隙流体流动的作用,而且需要考虑孔隙流体的流动对固体颗粒骨架变形产生的影响。同时,由于多孔介质材料具有强烈的非均质特性,采用有限元方法往往需要精细的有限元网格才能准确反映材料的微观特性和满足计算精度要求。这使得数值分析计算量很大,从而导致计算效率较低甚至计算的失败。另外,对于多孔介质复杂的断裂问题,由于需要在考虑间断不连续的同时还要考虑多物理场之间的相互耦合作用,这也使得采用有限元方法的计算效率较低甚至很难实现。因此,结合基于连续介质局部理论的多尺度有限元方法的高效性和基于非局部理论的近场动力学方法处理间断不连续问题的优势,发展高效的新型的数值算法对于理解和研究土体和岩石等多孔介质在力学方面的特性和机制以及在工程应用中的设计、稳定性和安全性的检测和评估等都具有重要的意义。本文基于多尺度有限元和近场动力学方法的基础理论,对单相固体断裂问题、饱和多孔介质液固耦合固结和动力响应及水力压裂等问题的多尺度有限元方法和近场动力学方法展开研究,主要研究包括以下五个部分:第一、针对均质和非均质饱和多孔介质液固耦合动力问题,提出了扩展多尺度有限元方法。不同于静力问题的扩展多尺度有限元方法,对于动力问题,先对每个单胞在时间和空间上进行有限元离散,然后利用单胞的等效静力平衡方程来进行升尺度和降尺度计算。其中,采用单胞的等效刚度矩阵和给定的线性边界条件来求解等效数值基函数,因此该数值基函数不仅能反映单胞的静力和动力惯性特征,而且还能反映固体骨架和孔隙流体之间的相互耦合作用。其次,给出了求解等效解数值基函数的一般化公式。提出了改进的降尺度计算公式和采用了多节点单胞技术来提高计算精度。最后,通过均质和非均质饱和多孔介质动力问题的算例验证了算法的正确性和高效性。第二、针对饱和多孔介质液固耦合固结和动力问题,提出了近场动力学方法。首先,根据饱和多孔介质有效应力原理和近场动力学基本理论,建立饱和多孔介质液固耦合固结和动力问题的近场动力学模型,并给出相应的控制方程。然后,根据一阶泰勒展开技术对控制方程进行线性化处理,得到近场动力学方法隐式的求解列式。最后,通过典型的饱和多孔介质固结和动力问题的数值算例验证了算法的正确性和有效性。第三、针对单相固体裂纹动态扩展问题,提出了耦合有限元和近场动力学方法。首先,基于键型近场动力学模型,利用一阶泰勒展开技术、Newmark积分方法和Newton-Raphson方法推导得到动力问题的近场动力学方法等效的隐式增量方程。其次,将整个求解区域分解成有限元列式求解子域和近场动力学列式求解子域,并根据采用基于网格节点和质点直接进行耦合的策略得到系统耦合的等效的隐式增量方程。其中,前者主要为包含边界条件或无裂纹扩展的区域,而后者主要为包含存在裂纹和裂纹扩展的区域。然后,针对此非线性问题,给出位移/载荷控制的增量迭代求解算法和相应的计算流程。最后,通过数值算例验证了耦合方法的正确性和有效性。第四、针对单相固体裂纹准静态扩展问题,提出了耦合扩展多尺度有限元和近场动力学方法。首先,基于常规态型近场动力学模型,利用泰勒展开的线性化处理方法推导得到相应的隐式增量方程。其次,将整个求解域划分为扩展多尺度有限元列式求解子域和近场动力学列式求解子域两个部分。在耦合区域中,采用扩展多尺度有限元单胞的数值基函数来建立扩展多尺度有限元单胞的节点和近场动力学的质点之间的相互约束关系。同时,采用拉格朗日乘子法将此约束条件引入到耦合的应变能密度函数中,并利用变分原理得到耦合的控制方程。然后,在近场动力学裂纹扩展分析中,采用双线性软化的材料模型来描述材料从损伤到破坏的过程,并同时给出相应的位移/载荷控制的增量迭代求解算法。最后,通过数值算例验证了耦合方法的正确性和高效性。最后,针对饱和多孔介质水力压裂问题,发展了耦合扩展多尺度有限元和近场动力学方法。首先,提出了饱和多孔介质水力压裂问题的近场动力学模型,即固体骨架运动方程、孔隙流体渗流方程和裂隙内流体运动方程。采用基于泰勒展开技术的线性化方法、Newmark方法和Newton-Raphson方法,推导得到相应控制方程的增量形式的等效静力平衡方程。其次,基于饱和多孔介质动力问题的连续介质方程,推导得到了饱和多孔介质动力问题的扩展多尺度有限元方法的增量形式的等效静力平衡方程。然后,在耦合策略中,采用拉格朗日乘子法将耦合区域中单胞的节点和质点之间的相互约束关系引入到耦合的应变能密度函数中,再利用变分原理推导得到系统耦合的等效静力平衡方程。同时,针对该非线性问题,给出位移/载荷控制的增量迭代求解算法。最后,采用水力压裂问题的数值算例验证了提出的耦合方法的正确性和有效性。
【图文】：

焦油砂,多孔介质,毛细血管网,页岩

１．１研究背景和意义逡逑多孔介质材料，如页岩⑴、毛细血管网络［２］、焦油砂［３］、活性炭［４ｉ、大豆蛋白多孔水逡逑凝胶Ｗ等（如图１．１所示），作为一种由固体颗粒骨架、大量微小孔隙以及孔隙中的液逡逑相和气相组成的组合体广泛存在于自然环境和工程应用中。尤其是土壤和岩石，，作为一逡逑种常见和价廉的材料与人们的生产生活息息相关，比如自然和环境工程中涉及到的房屋逡逑和桥梁的地基、水利水坝、深埋管道的掩体和渣土堆填体等。多孔介质是一种由多相物逡逑质共存的组合体，主要包括固相、液相和气相。如果多孔介质的孔隙中为干的，则称其逡逑为干的多孔介质？，如果多孔介质的孔隙中为饱和的流体相，则称其为饱和多孔介质？，如逡逑果多孔介质的孔隙中为流体相和气相共存，则称其为非饱和多孔介质。针对于饱和多孔逡逑介质，不仅需要考虑固体颗粒骨架和孔隙流体之间的相互耦合作用关系，即固体颗粒骨逡逑架的变形将直接影响孔隙流体的流动状态，而且孔隙流体的流动对固体颗粒骨架的变形逡逑也会产生影响。因此

核燃料,核电站,日本,事件

失时便发生液化现象。另外，在大力倡导保护环境的时代背景下，如何有效防止被工业逡逑废水废物、生活垃圾、农业化肥和农药等污染的地表水渗入到地下也是现在必须面对和逡逑急需解决的问题。如图１．２所示，在２０１１年日本福岛核电站在地震和海啸的作用下发生逡逑核泄露事故，其中被核燃料污染的水渗入地下对当地环境造成极大伤害，使得灾区２０公逡逑里范围内的百姓不得不离^u家园。从本质上来讲，以上提及的地基沉降、土壤液化和土逡逑壤渗流等问题都属于多孔介质的固结和动力等力学问题。因此，多孔介质固结和动力力逡逑学行为特征和机制的研宄，对于合理、有效和科学地解决上述实际工程问题具有重要指逡逑导意义。逡逑图１．２日本福岛核电站核燃料泄露事件ｗ逡逑Ｆｉｇ．邋１．２邋Ｎｕｃｌｅａｒ邋ｆｕｅｌ邋ｌｅａｋ邋ａｔ邋Ｆｕｋｕｓｈｉｍａ邋Ｎｕｃｌｅａｒ邋Ｐｏｗｅｒ邋Ｐｌａｎｔ邋ｉｎ邋Ｊａｐａｎ１６】逡逑－２－逡逑
【学位授予单位】：大连理工大学
【学位级别】：博士
【学位授予年份】：2019
【分类号】：TU43

【参考文献】