近年来,在“一带一路”战略以及十四五布局规划引领下,大体量复杂岩体工程项目逐渐增多,随之而来的技术难度增大。针对大尺度范围岩体工程的特点以及岩体研究现状,基于等效连续介质理论,从本构模型及其有限元程序开发两方面对含低序次裂隙系统岩体中最为典型的饱和多组节理岩体进行了研究。不仅完善了岩体本构理论体系,还为该类岩体工程数值分析及设计提供了一种新的计算模型,利于分析优化比对。首先,明确基于等效连续介质理论研究的前提是多组节理岩体可等效为多孔连续介质。通过对岩体单元的渗透性参数和力学性参数进行类常量和类张量判定,检验岩体REV的存在性,以此确认岩体是否可以等效为多孔连续介质,并以REV尺寸对应的岩体属性参数作为等效连续性参数。通过算例,展示了多组节理岩体等效连续性判定步骤,证明节理间距较密集均匀时,此类包含低序次裂隙系统的岩体可视作多孔连续介质进行研究。其次,克服以往饱和多孔介质本构理论中不能合理考虑固相基质变形的缺点。基于具有普适性的机械功守恒原理和工程混合物理论,通过推导和分析具有普适性质的饱和多孔介质体积变形功守恒方程,提出小应变条件下考虑基质压缩性的饱和多孔介质双应力原理,即Terzaghi有效应力唯一决定固相体积分数应变;固相基质压力唯一决定固相基质体应变;Terzaghi有效应力与固相基质压力共同决定固相体应变。双应力原理为可等效为多孔连续介质的饱和岩体本构模型的建立奠定了理论基础。再次,建立了泊松比随主应变变化的完整岩块损伤本构模型以及刚度随节理变形和节理应力变化的动态节理本构模型,讨论了节理变形对完整岩块的耦合损伤作用,通过整合获得多组节理岩体的非线性本构模型。利用混合物理论和双应力原理,探究了饱和多孔介质应力-渗流耦合作用,并综合多组节理岩体非线性本构模型特点,建立了饱和多组节理岩体本构模型。最后,鉴于大部分可由用户自定义材料的有限元软件(如ABAQUS的UMAT子程序)在计算岩土力学问题时,采用的本构模型都是基于有效应力原理而不考虑固流两相基质压缩变形。故为了运用双应力原理而开发新型有限元程序,由饱和多组节理岩体本构模型推导出岩体应力-渗流耦合基本方程,并利用Galerkin加权余量法和二元函数分部积分法进行处理,得到有限元方程和单元刚度矩阵,进行有限元程序开发。程序考虑了完整岩块和节理水压缩变形、岩体非线性以及应力-渗流耦合等因素,并被用以模拟深海岩体隧道开挖。数值分析结果表明,由提出的饱和岩体本构模型建立的有限元程序可以较真实地模拟非线性饱和岩体的工程响应,为设计施工提供借鉴。