CN115146565A - 一种基于内聚力算法的双曲拱坝坝基渗流评估方法及系统 - Google Patents
一种基于内聚力算法的双曲拱坝坝基渗流评估方法及系统 Download PDFInfo
- Publication number
- CN115146565A CN115146565A CN202211071029.XA CN202211071029A CN115146565A CN 115146565 A CN115146565 A CN 115146565A CN 202211071029 A CN202211071029 A CN 202211071029A CN 115146565 A CN115146565 A CN 115146565A
- Authority
- CN
- China
- Prior art keywords
- dam
- seepage
- arch dam
- hyperbolic
- hyperbolic arch
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种基于内聚力算法的双曲拱坝坝基渗流评估方法及系统,属于双曲拱坝坝基渗流监测技术领域。现有技术缺乏高效可靠的双曲拱坝坝基渗流计算方案。本发明的一种基于内聚力算法的双曲拱坝坝基渗流评估方法,构建基于内聚力算法的热渗流应力耦合有限元模型,用于对双曲拱坝的温度场、应力场和渗流场进行耦合计算,其采用内聚力算法模拟坝体和坝基交界面处的微裂缝和孔隙,能准确地掌握坝基渗流场的分布情况,实现双曲拱坝的结构安全状态的准确评估,无需进行坝体和坝基交界面的网格细化,可以有效降低网格数量,提高计算效率,并具有更好的计算稳定性,方案科学合理,为双曲拱坝坝基渗流监测和评估提供一个高效可靠的处理方法。
Description
技术领域
本发明涉及一种基于内聚力算法的双曲拱坝坝基渗流评估方法及系统,属于双曲拱坝坝基渗流监测技术领域。
背景技术
双曲拱坝是适合于修建在“U”形及“V”形河谷中的经济性与安全性都比较优越的坝型,同时又是挡水坝中结构最复杂的坝型,是高拱坝的发展趋势。在双曲拱坝的坝基岩体中,由于渗透荷载与岩体自重属同一数量级,可能会直接导致坝基失稳等危害。因此,对双曲拱坝的坝基渗流特性和渗流控制的研究就显得至关重要。
为了更加清楚地了解坝基各部位的扬压力和渗透力大小、岩体的抗渗透变形能力、以及岩体中渗流场和应力场的相互作用特性等,对渗流场进行计算就成了必不可少的工作。
同时坝基渗流场、压力梯度分布和渗流量等物理量对双曲拱坝的稳定性具有较大的影响,准确地掌握坝基渗流场的分布情况,有助于对双曲拱坝的结构安全状态进行评估。
国内外学者在大坝岩体渗流场和应力场计算方面进行了大量的实验和理论研究,总结了单一裂隙在应力作用下的理论模型,提出多种经验公式和间接公式,进而提出各种不同的渗流应力耦合模型并进行相应的数值理论分析。在数值理论分析方法中,最有效的一种就是有限元求解方法,如剩余流量法、初流量法、节点虚流量法、调整渗透系数法等。
但由于坝体和坝基交界面位置的地质构造复杂,造成计算方法中的高度非线性,以及渗流系数、岩体孔隙率等参数难以准确获取,对于坝基渗流场的仿真计算目前仍缺少有效的解决方案。
进一步,中国专利(公布号:CN 113447414 A)公开了一种土石坝绕坝及坝基渗流及渗流量实时监测装置及其使用方法,综合利用光纤监测技术以及电磁感应监测技术对反映土石坝绕坝及坝基渗流状态的多种参数进行监测,获得监测位置土石坝绕坝及坝基渗流状态的综合评价结果。另外,基于对土石坝坝内监测位置的长期实时监测数据,开展对土石坝绕坝及坝基渗流状态的综合分析,并根据参数监测结果对土石坝绕坝及坝基渗流状态进行预警。
上述方案结合光纤监测技术与电磁感应监测技术,实现对土石坝坝内监测点多种渗流参数的实时同步监测,但利用机械装置监测会存在较多的监测盲点,无法对坝体进行全面有效的监控;同时上述方案没有公开如何利用监测数据进行数据分析,得到坝基渗流值以及渗漏场的分布情况。
更进一步,上述方案以及现有技术,更多的是针对土石坝和重力坝的研究,对于双曲拱坝的坝基渗流场和应力场的耦合监测,以及在双曲拱坝运维期间的热渗流应力耦合的仿真计算仍缺乏高效可靠的处理方案,进而无法对双曲拱坝的结构安全状态进行准确评估,不利于推广使用。
发明内容
针对现有技术的缺陷,本发明的目的一在于提供一种通过构建基于内聚力算法的热渗流应力耦合有限元模型,并根据边界条件,得到双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况;进而根据双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况,对双曲拱坝的结构安全状态进行准确评估,方案科学合理,切实可行,适用于双曲拱坝坝基渗流监测和评估的基于内聚力算法的双曲拱坝坝基渗流评估方法。
本发明的目的二在于提供一种通过构建基于内聚力算法的热渗流应力耦合有限元模型,用于对双曲拱坝的温度场、应力场和渗流场进行耦合计算,其采用内聚力算法模拟坝体和坝基交界面处的微裂缝和孔隙,并通过带有孔隙水压力自由度的三维二次单元计算双曲拱坝坝基的渗流场,无需进行坝体和坝基交界面的网格细化,可以有效降低网格数量,提高计算效率,并具有更好的计算稳定性,方案科学合理,切实可行,能够对双曲拱坝坝基渗流进行有效监测和评估的基于内聚力算法的双曲拱坝坝基渗流评估方法。
本发明的目的三在于提供一种构建双曲拱坝设计参数存储模块、热渗流应力耦合有限元模块、双曲拱坝数据监测模块、边界条件处理模块、安全评估模块,对双曲拱坝的温度场、应力场和渗流场进行耦合计算,方案科学合理,切实可行,能够对双曲拱坝坝基渗流进行有效监测和评估的基于内聚力算法的双曲拱坝坝基渗流评估系统。
为实现上述目的之一,本发明的第一种技术方案为:
一种基于内聚力算法的双曲拱坝坝基渗流评估方法,
包括以下内容:
利用双曲拱坝监测数据生成热力学载荷数据和结构力学载荷数据;
获得边界条件;
将边界条件输入到预先构建的基于内聚力算法的热渗流应力耦合有限元模型中,得到双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况;
根据双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况,对双曲拱坝的结构安全状态进行评估。
本发明经过不断探索以及试验,通过构建基于内聚力算法的热渗流应力耦合有限元模型,并根据边界条件,得到双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况;进而根据双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况,能够对双曲拱坝的结构安全状态进行准确评估,方案科学合理,切实可行,利于推广使用,特别适用于双曲拱坝坝基的渗流监测和评估。
进一步,本发明使用内聚力算法进行双曲拱坝的渗流计算,内聚力算法可以在裂纹扩展过程中形成新的裂纹表面,在新裂纹面不连续的情况下,在数学模型上保持连续性条件,从而对裂纹的萌生、扩展建立统一的计算模型来描述;并通过构建热渗流应力耦合有限元模型可以进行双曲拱坝温度场、应力场和渗流场的耦合计算,能够同时考虑温度场、坝基应力场和渗流场之间的相互作用,进而可以完成双曲拱坝坝基的三维建模、边界条件设置、结果后处理和数据分析,从而能准确地掌握坝基渗流场的分布情况,实现双曲拱坝的结构安全状态的准确评估。
作为优选技术措施:
所述基于内聚力算法的热渗流应力耦合有限元模型,对双曲拱坝的温度场、应力场和渗流场进行耦合计算,其采用内聚力算法模拟坝体和坝基交界面处的微裂缝和孔隙,并通过带有孔隙水压力自由度的三维二次单元计算双曲拱坝坝基的渗流场。
所述三维二次单元包括四面体网格形状或/和六面体网格形状或/和三棱柱网格形状。
作为优选技术措施:
构建所述基于内聚力算法的热渗流应力耦合有限元模型,包括:
设置双曲拱坝拱圈的中心线参数方程,生成双曲拱坝中心线;
沿着双曲拱坝中心线,双曲拱坝法向厚度由一组随拱坝中心角变化的参数方程控制,所述参数方程为双曲拱坝中心线上任意点处的法向拱厚参数方程;
利用拱坝高程函数以及方程系数,计算中心线参数方程和法向拱厚参数方程中的未知参数,得到未知参数的参数值;
将参数值以及拱坝中心角带入中心线参数方程、拱厚参数方程,确定双曲拱坝坝体体型;
根据双曲拱坝坝体体型,通过参数曲线放样得到三维几何模型;
根据三维几何模型,通过坝体与坝基交界面处的几何拉伸自动生成坝基几何模型。
作为优选技术措施:
还包括构建三维几何模型,其包括:
根据双曲拱坝体型参数,采用具有六面体网格形状的三维二次单元,对双曲拱坝进行网格划分;
采用薄层内聚力单元生成坝体与坝基的交界面。
作为优选技术措施:
所述边界条件包括:
双曲拱坝温度场计算所需温度边界、上下游水位产生的静水压力和上游库底泥沙淤积压力;
双曲拱坝温度场计算所需温度边界包括坝体外部与空气接触部分温度边界,水位以下坝面温度边界,河谷、河岸在水位以下部分温度边界,河谷、河岸在水位以上部分温度边界,地基深部温度边界。
作为优选技术措施:
坝体外部与空气接触部分温度边界采用传热学第三类边界,在接触面产生对流换热;
水位以下坝面温度边界等效为第一类边界条件,相应温度为实测水温;
河谷、河岸在水位以下部分温度边界,按照第一类边界条件处理,取值为相应高程实测水温;
河谷、河岸在水位以上部分温度边界,将其转化为第一类边界条件,具体取值为每旬平均气温;
地基深部温度边界受地温影响,地基侧面按绝热边界处理;
上游库底泥沙淤积压力根据泥沙淤积压力标准值、淤积泥沙的浮容重、坝前泥沙淤积厚度和淤积泥沙内摩擦角进行计算。
作为优选技术措施:
基于内聚力算法的热渗流应力耦合有限元模型采用牛顿-拉夫森迭代方法进行求解,其包括温度场的计算、应力场的计算、渗流场的计算:
所述温度场采用热对流方程进行计算;
所述应力场采用内聚力单元内部的基于立方定律建立的平衡方程进行计算;
所述渗流场采用菲克扩散定律进行计算。
作为优选技术措施:
在双曲拱坝的结构安全状态的评估前,将双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况与已有的渗流监测数据,进行对比分析,以确定基于内聚力算法的热渗流应力耦合有限元模型仿真计算是否准确;
所述对比分析包括对不同测点的渗流量数据进行对比和残差分析。
为实现上述目的之一,本发明的第二种技术方案为:
一种基于内聚力算法的双曲拱坝坝基渗流评估方法,包括以下步骤:
步骤1:根据双曲拱坝体型参数,生成三维几何模型,构建基于内聚力算法的热渗流应力耦合有限元模型;
所述基于内聚力算法的热渗流应力耦合有限元模型,用于对双曲拱坝的温度场、应力场和渗流场进行耦合计算,其采用内聚力算法模拟坝体和坝基交界面处的微裂缝和孔隙,并通过带有孔隙水压力自由度的三维二次单元计算双曲拱坝坝基的渗流场;
步骤2:利用双曲拱坝监测数据生成热力学载荷数据和结构力学载荷数据;
根据热力学载荷数据和结构力学载荷数据构建步骤1中的基于内聚力算法的热渗流应力耦合有限元模型的边界条件;
步骤3:根据步骤2中的边界条件,求解基于内聚力算法的热渗流应力耦合有限元模型,得到双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况;
步骤4:利用步骤3中的双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况,实现双曲拱坝的结构安全状态的评估。
本发明经过不断探索以及试验,构建基于内聚力算法的热渗流应力耦合有限元模型,用于对双曲拱坝的温度场、应力场和渗流场进行耦合计算,其采用内聚力算法模拟坝体和坝基交界面处的微裂缝和孔隙,并通过带有孔隙水压力自由度的三维二次单元计算双曲拱坝坝基的渗流场,无需进行坝体和坝基交界面的网格细化,可以有效降低网格数量,提高计算效率,并具有更好的计算稳定性,方案科学合理,切实可行,为双曲拱坝坝基渗流监测和评估提供一个高效可靠的处理方法。
进一步,本发明使用内聚力算法进行双曲拱坝的渗流计算,内聚力算法可以在裂纹扩展过程中形成新的裂纹表面,在新裂纹面不连续的情况下,在数学模型上保持连续性条件,从而对裂纹的萌生、扩展建立统一的计算模型来描述;并结合基于内聚力算法的热渗流应力耦合有限元模型,进行双曲拱坝温度场、应力场和渗流场的耦合计算,能够同时考虑温度场、坝基应力场和渗流场之间的相互作用,并能完成三维建模、边界条件设置、结果后处理和数据分析的整个流程,从而能准确地掌握坝基渗流场的分布情况,实现双曲拱坝的结构安全状态的准确评估。
为实现上述目的之一,本发明的第三种技术方案为:
一种基于内聚力算法的双曲拱坝坝基渗流评估系统,
应用上述的一种基于内聚力算法的双曲拱坝坝基渗流评估方法;
其包括双曲拱坝设计参数存储模块、热渗流应力耦合有限元模块、双曲拱坝数据监测模块、边界条件处理模块、安全评估模块;
所述双曲拱坝设计参数存储模块,用于存储双曲拱坝体型参数;
所述热渗流应力耦合有限元模块,用于对双曲拱坝的温度场、应力场和渗流场进行耦合计算;
所述双曲拱坝数据监测模块,用于对双曲拱坝的状态进行监测;
所述边界条件处理模块,用于获得双曲拱坝的边界条件;
所述安全评估模块,用于对双曲拱坝的结构安全状态进行评估;
所述双曲拱坝设计参数存储模块将双曲拱坝体型参数输入到热渗流应力耦合有限元模块中;
所述双曲拱坝数据监测模块将监测数据传输到边界条件处理模块中,生成边界约束和载荷数据,并将边界约束和载荷数据输入到热渗流应力耦合有限元模块中;
所述热渗流应力耦合有限元模块根据双曲拱坝体型参数以及边界约束和载荷数据,对双曲拱坝的温度场、应力场和渗流场进行耦合计算,并将计算结果传输给安全评估模块;
所述安全评估模块根据热渗流应力耦合有限元模块的计算结果对双曲拱坝的结构安全状态进行评估。
本发明经过不断探索以及试验,构建双曲拱坝设计参数存储模块、热渗流应力耦合有限元模块、双曲拱坝数据监测模块、边界条件处理模块、安全评估模块,对双曲拱坝的温度场、应力场和渗流场进行耦合计算,方案科学合理,切实可行,为双曲拱坝坝基渗流监测和评估提供一个高效可靠的系统,进而能准确地掌握坝基渗流场的分布情况,实现双曲拱坝的结构安全状态的准确评估。
与现有技术相比,本发明具有以下有益效果:
本发明经过不断探索以及试验,通过构建基于内聚力算法的热渗流应力耦合有限元模型,并根据边界条件,得到双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况;进而根据双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况,能够对双曲拱坝的结构安全状态进行准确评估,方案科学合理,切实可行,利于推广使用,特别适用于双曲拱坝坝基的渗流监测和评估。
进一步,本发明构建基于内聚力算法的热渗流应力耦合有限元模型,用于对双曲拱坝的温度场、应力场和渗流场进行耦合计算,其采用内聚力算法模拟坝体和坝基交界面处的微裂缝和孔隙,并通过带有孔隙水压力自由度的三维二次单元计算双曲拱坝坝基的渗流场,无需进行坝体和坝基交界面的网格细化,可以有效降低网格数量,提高计算效率,并具有更好的计算稳定性,方案科学合理,切实可行,为双曲拱坝坝基渗流监测和评估提供一个高效可靠的处理方法。
再进一步,本发明使用内聚力算法进行双曲拱坝的渗流计算,内聚力算法可以在裂纹扩展过程中形成新的裂纹表面,在新裂纹面不连续的情况下,在数学模型上保持连续性条件,从而对裂纹的萌生、扩展建立统一的计算模型来描述;并结合基于内聚力算法的热渗流应力耦合有限元模型,进行双曲拱坝温度场、应力场和渗流场的耦合计算,能够同时考虑温度场、坝基应力场和渗流场之间的相互作用,并能完成三维建模、边界条件设置、结果后处理和数据分析的整个流程,从而能准确地掌握坝基渗流场的分布情况,实现双曲拱坝的结构安全状态的准确评估。
更进一步,本发明构建双曲拱坝设计参数存储模块、热渗流应力耦合有限元模块、双曲拱坝数据监测模块、边界条件处理模块、安全评估模块,对双曲拱坝的温度场、应力场和渗流场进行耦合计算,方案科学合理,切实可行,为双曲拱坝坝基渗流监测和评估提供一个高效可靠的系统,进而能准确地掌握坝基渗流场的分布情况,实现双曲拱坝的结构安全状态的准确评估。
附图说明
图1为本发明双曲拱坝坝基渗流评估方法的一种流程图;
图2为本发明生成的双曲拱坝中心线的一种示图;
图3为本发明生成的双曲拱坝上下游表面轮廓包络线的一种示图;
图4为本发明生成的双曲拱坝坝体的一种结构示图;
图5为本发明生成的坝基几何模型的一种结构示图;
图6为本发明双曲拱坝的一种结构示图;
图7为本发明内聚力单元网格的一种结构示图;
图8为本发明双曲拱坝温度场的一种示图;
图9为本发明双曲拱坝温度场的另一种示图;
图10为本发明渗流量计算结果的一种示图;
图11为本发明数据对比分析的一种示图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
相反,本发明涵盖任何由权利要求定义的在本发明的精髓和范围上做的替代、修改、等效方法以及方案。进一步,为了使公众对本发明有更好的了解,在下文对本发明的细节描述中,详尽描述了一些特定的细节部分。对本领域技术人员来说没有这些细节部分的描述也可以完全理解本发明。
本发明双曲拱坝坝基渗流评估方法的第一种具体实施例:
一种基于内聚力算法的双曲拱坝坝基渗流评估方法,包括以下内容:
利用双曲拱坝监测数据生成热力学载荷数据和结构力学载荷数据,获得边界条件;
将边界条件输入到预先构建的基于内聚力算法的热渗流应力耦合有限元模型中,得到双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况;
根据双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况,对双曲拱坝的结构安全状态进行评估。
本发明双曲拱坝坝基渗流评估方法的第二种具体实施例:
一种基于内聚力算法的双曲拱坝坝基渗流评估方法,包括以下步骤:
步骤1:根据双曲拱坝体型参数,生成三维几何模型,构建基于内聚力算法的热渗流应力耦合有限元模型;
所述基于内聚力算法的热渗流应力耦合有限元模型,用于对双曲拱坝的温度场、应力场和渗流场进行耦合计算,其采用内聚力算法模拟坝体和坝基交界面处的微裂缝和孔隙,并通过带有孔隙水压力自由度的三维二次单元计算双曲拱坝坝基的渗流场;
步骤2:利用双曲拱坝监测数据生成热力学载荷数据和结构力学载荷数据;
根据热力学载荷数据和结构力学载荷数据构建步骤1中的基于内聚力算法的热渗流应力耦合有限元模型的边界条件;
步骤3:根据步骤2中的边界条件,求解基于内聚力算法的热渗流应力耦合有限元模型,得到双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况;
步骤4:利用步骤3中的双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况,实现双曲拱坝的结构安全状态的评估。
本发明基于内聚力算法的热渗流应力耦合有限元模型的一种具体实施例::
所述步骤1中的基于内聚力算法的热渗流应力耦合有限元模型,其构建方法如下:
第一步,设置双曲拱坝拱圈的中心线参数方程,生成双曲拱坝中心线;
第二步,沿着第一步中的双曲拱坝中心线,双曲拱坝法向厚度由一组随拱坝中心角变化的参数方程控制,所述参数方程为双曲拱坝中心线上任意点 处的法向拱厚参数方程;
第三步,利用拱坝高程函数以及方程系数,计算第一步的中心线参数方程和第二步的法向拱厚参数方程中的未知参数,得到未知参数的参数值;
第四步,将第三步中的参数值以及拱坝中心角带入中心线参数方程、拱厚参数方程,确定双曲拱坝坝体体型;
第五步,根据第四步中双曲拱坝坝体体型,通过参数曲线放样得到三维几何模型;
第六步,根据第五步中三维几何模型,通过坝体与坝基交界面处的几何拉伸自动生成坝基几何模型。
本发明中心线参数方程的一种具体实施例::
所述中心线参数方程的计算公式如下:
X为拱圈中心线某点的横坐标,Y为拱圈中心线某点的纵坐标;
所述法向拱厚参数方程的计算公式如下:
所述拱坝高程函数的表达式如下:
其中,Z为拱坝高程;M为6个决定拱坝中心线方程的参数,分别为::拱冠梁上游
面的Y坐标,:拱冠梁处的拱圈厚度,:左拱端拱圈厚度,:右拱端拱圈厚度,:拱冠
梁处左岸侧的拱圈中心线曲率半径,:拱冠梁处右岸侧的拱圈中心线曲率半径;A/B/C/D
为方程系数。
本发明双曲拱坝的一种具体实施例::
所述双曲拱坝的坝基上下游宽度为1.5倍坝高,左右岸宽度为1.5倍坝高,坝基深度为2倍坝高;
同时对双曲拱坝的三维几何模型进行网格划分,整体采用20节点六面体网格形状的三维二次单元,坝体和坝基连接位置采用薄层内聚力单元生成。
本发明边界条件的一种具体实施例:
所述步骤2中,基于内聚力算法的热渗流应力耦合有限元模型的边界条件包括:
双曲拱坝温度场计算所需温度边界、上下游水位产生的静水压力和上游库底泥沙淤积压力;
双曲拱坝温度场计算所需温度边界包括坝体外部与空气接触部分温度边界,水位以下坝面温度边界,河谷、河岸在水位以下部分温度边界,河谷、河岸在水位以上部分温度边界,地基深部温度边界。
坝体外部与空气接触部分温度边界采用传热学第三类边界,在接触面产生对流换热按照下式计算:
水位以下坝面温度边界等效为第一类边界条件,相应温度为实测水温;
河谷、河岸在水位以下部分温度边界,按照第一类边界条件处理,取值为相应高程实测水温;
河谷、河岸在水位以上部分温度边界,将其转化为第一类边界条件,具体取值为每旬平均气温;
地基深部温度边界受地温影响,地基侧面按绝热边界处理,即:
上游库底泥沙淤积压力按照以下公式计算:
本发明模型求解的一种具体实施例:
所述步骤3中,基于内聚力算法的热渗流应力耦合有限元模型采用牛顿-拉夫森迭代方法进行求解,其包括温度场的计算、应力场的计算、渗流场的计算:
所述温度场采用热对流方程进行计算;
所述应力场采用内聚力单元内部的基于立方定律建立的平衡方程进行计算;
所述渗流场采用菲克扩散定律进行计算。
所述热对流方程的计算公式如下:
所述应力场的计算,包括以下步骤:
S1,通过节点位移插值得到裂隙宽度,其计算公式如下:
R为局部坐标系下的节点位移矩阵。
采用指数型内聚力张力位移关系,开裂过程中的断裂能控制方程为:
其中,为断裂过程中的总断裂能值,、分别为内聚力单元界面上的法向
与切向位移值,为纯法向开裂状态下界面完全开裂时的界面断裂能,、为法向与切
向界面开裂特征位移,即应力最大值点对应的位移值,q,r分别为:
内聚力单元界面上的各向应力由断裂能控制方程的偏导得到:
S2,根据S1中的裂隙宽度,并利用立方定律建立平衡方程,所述平衡方程的计算公式如下:
S3,采用牛顿-拉夫森迭代方法求解S2中的平衡方程,在求解过程中关闭矩阵奇异性检查;
所述菲克扩散定律的计算公式如下:
本发明判断模型是否准确的一种具体实施例:
所述步骤4中,在双曲拱坝的结构安全状态的评估前,将双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况与已有的渗流监测数据,进行对比分析,以确定基于内聚力算法的热渗流应力耦合有限元模型仿真计算是否准确;
所述对比分析包括对不同测点的渗流量数据进行对比和残差分析。
如图 1所示,本发明双曲拱坝坝基渗流评估方法的一种最佳具体实施例:
一种基于内聚力算法的双曲拱坝坝基渗流评估方法,包括如下步骤:
步骤1:根据双曲拱坝体型参数,生成三维几何模型,构建基于内聚力算法的热渗流应力耦合有限元模型。
所述基于内聚力算法的热渗流应力耦合有限元模型,其构建方法如下:
(1)根据以下双曲拱坝拱圈中心线参数方程,生成如图2所示的双曲拱坝中心线:
(2)沿着拱圈中心线,根据以下双曲拱坝法向厚度参数方程,生成如图3所示的双曲拱坝上下游表面轮廓的包络线:
(3)上述方程中的双曲拱坝体型参数是拱坝高程的函数,函数表达式如下:
其中,Z为拱坝高程,m;M为6个决定拱坝中心线方程的参数,分别为::拱冠梁上
游面的Y坐标(工程坐标系),:拱冠梁处的拱圈厚度,:左拱端拱圈厚度,:右拱端拱圈
厚度,:拱冠梁处左岸侧的拱圈中心线曲率半径,:拱冠梁处右岸侧的拱圈中心线曲率
半径,A/B/C/D为方程系数,具体取值如表1所示。
表1 双曲拱坝体型参数方程系数
(4)针对不同的拱坝中心线参数,具有一组不同的方程系数,通过方程系数和拱坝中心角的确定就可以确定双曲拱坝坝体体型,通过参数曲线放样得到三维几何模型。通过坝体与坝基交界面处的几何拉伸自动生成坝基几何模型,坝基上下游宽度为1.5倍坝高,左右岸宽度为1.5倍坝高,坝基深度为2倍坝高。生成的双曲拱坝坝体和坝基几何模型如图4、图5所示。
(5)对双曲拱坝的三维几何模型进行网格划分,整体采用20节点六面体网格形状的三维二次单元,坝体和坝基连接位置采用薄层内聚力单元。双曲拱坝和内聚力单元网格划分如图6、图7所示。
步骤2:读取双曲拱坝监测数据,设置所述的基于内聚力算法渗流应力耦合有限元模型的边界条件。所述双曲拱坝监测数据包括热力学载荷数据、结构力学载荷、渗流量。
所述边界条件包括:
双曲拱坝温度场计算所需温度边界(上游坝面不同高程温度边界、坝顶和下游坝面空气自然对流系数和环境温度边界)、上下游水位产生的静水压力和上游库底泥沙淤积压力。
所述坝顶和下游坝面空气自然对流系数和环境温度边界,其坝体外部与空气接触部分采用传热学第三类边界,在接触面产生对流换热按照下式计算:
表2 双曲拱坝坝体上游坝面不同高程温度边界
所述上游坝面不同高程温度边界(可参见表2),其水位以下坝面由于库水与混凝土表面放热系数很大,所以可以将其等效为第一类边界条件,相应温度为实测水温。河谷、河岸在水位以下部分,按照第一类边界条件处理,取值为相应高程实测水温。河谷、河岸在水位以上部分,将其转化为第一类边界条件,具体取值为每旬平均气温。地基深部温度主要受地温影响,地基侧面按绝热边界处理,即:
泥沙淤积压力按照以下公式计算:
步骤3:基于内聚力算法的热渗流应力耦合有限元模型的求解
(1)采用牛顿-拉夫森迭代方法,在求解过程中关闭矩阵奇异性检查;
(2)在基于内聚力算法的热渗流应力耦合有限单元内部,考虑立方定律建立平衡方程如下公式所示:
其他力学计算部分与传统内聚力单元一致。
(3)裂隙宽度由节点位移插值得到:
(4)在双曲拱坝坝基中计算渗流场,采用菲克扩散定律:
(5)温度场的计算采用热对流方程描述:
步骤4:根据步骤3的渗流仿真计算结果,进行渗流量、渗透压力等物理场分布显示,以及对不同测点的渗流量数据进行对比和残差分析。获得渗流量计算数据(图10),根据监测数据对渗流量进行数据对比分析,如图11所示。
本发明的内聚力算法是基于断裂力学裂纹萌生和扩展的一种技术。内聚力算法可以在裂纹扩展过程中形成新的裂纹表面,在新裂纹面不连续的情况下,在数学模型上保持连续性条件,从而对裂纹的萌生、扩展建立统一的计算模型来描述。使用内聚力算法进行双曲拱坝的渗流计算,无需进行坝体和坝基交界面的网格细化,可以有效降低网格数量,提高计算效率,并具有更好的计算稳定性。
本发明基于内聚力算法的双曲拱坝坝基渗流评估系统的一种具体实施例:
一种基于内聚力算法的双曲拱坝坝基渗流评估系统,应用上述的一种基于内聚力算法的双曲拱坝坝基渗流评估方法;其包括双曲拱坝设计参数存储模块、热渗流应力耦合有限元模块、双曲拱坝数据监测模块、边界条件处理模块、安全评估模块。
所述双曲拱坝设计参数存储模块,用于存储双曲拱坝体型参数;
所述热渗流应力耦合有限元模块,用于对双曲拱坝的温度场、应力场和渗流场进行耦合计算;
所述双曲拱坝数据监测模块,用于对双曲拱坝的状态进行监测;
所述边界条件处理模块,用于获得双曲拱坝的边界条件;
所述安全评估模块,用于对双曲拱坝的结构安全状态进行评估;
所述双曲拱坝设计参数存储模块将双曲拱坝体型参数输入到热渗流应力耦合有限元模块中;
所述双曲拱坝数据监测模块将监测数据传输到边界条件处理模块中,生成边界约束和载荷数据,并将边界约束和载荷数据输入到热渗流应力耦合有限元模块中;
所述热渗流应力耦合有限元模块根据双曲拱坝体型参数以及边界约束和载荷数据,对双曲拱坝的温度场、应力场和渗流场进行耦合计算,并将计算结果传输给安全评估模块;
所述安全评估模块根据热渗流应力耦合有限元模块的计算结果对双曲拱坝的结构安全状态进行评估。
应用本发明方法的一种设备实施例:
一种计算机设备,其包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现上述的一种基于内聚力算法的双曲拱坝坝基渗流评估方法。
应用本发明方法的一种计算机介质实施例:
一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述的一种基于内聚力算法的双曲拱坝坝基渗流评估方法。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。
Claims (10)
1.一种基于内聚力算法的双曲拱坝坝基渗流评估方法,其特征在于,
包括以下内容:
利用双曲拱坝监测数据生成热力学载荷数据和结构力学载荷数据,获得边界条件;
将边界条件输入到预先构建的基于内聚力算法的热渗流应力耦合有限元模型中,得到双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况;
根据双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况,对双曲拱坝的结构安全状态进行评估。
2.如权利要求1所述的一种基于内聚力算法的双曲拱坝坝基渗流评估方法,其特征在于,
所述基于内聚力算法的热渗流应力耦合有限元模型,对双曲拱坝的温度场、应力场和渗流场进行耦合计算,其采用内聚力算法模拟坝体和坝基交界面处的微裂缝和孔隙,并通过带有孔隙水压力自由度的三维二次单元计算双曲拱坝坝基的渗流场;
所述三维二次单元包括四面体网格形状或/和六面体网格形状或/和三棱柱网格形状。
3.如权利要求1所述的一种基于内聚力算法的双曲拱坝坝基渗流评估方法,其特征在于,
构建所述基于内聚力算法的热渗流应力耦合有限元模型,包括:
设置双曲拱坝拱圈的中心线参数方程,生成双曲拱坝中心线;
沿着双曲拱坝中心线,双曲拱坝法向厚度由一组随拱坝中心角变化的参数方程控制,所述参数方程为双曲拱坝中心线上任意点处的法向拱厚参数方程;
利用拱坝高程函数以及方程系数,计算中心线参数方程和法向拱厚参数方程中的未知参数,得到未知参数的参数值;
将参数值以及拱坝中心角带入中心线参数方程、拱厚参数方程,确定双曲拱坝坝体体型;
根据双曲拱坝坝体体型,通过参数曲线放样得到三维几何模型;
根据三维几何模型,通过坝体与坝基交界面处的几何拉伸自动生成坝基几何模型。
4.如权利要求1所述的一种基于内聚力算法的双曲拱坝坝基渗流评估方法,其特征在于,
还包括构建三维几何模型,其包括:
根据双曲拱坝体型参数,采用具有六面体网格形状的三维二次单元,对双曲拱坝进行网格划分;
采用薄层内聚力单元生成坝体与坝基的交界面。
5.如权利要求1所述的一种基于内聚力算法的双曲拱坝坝基渗流评估方法,其特征在于,
所述边界条件包括:
双曲拱坝温度场计算所需温度边界、上下游水位产生的静水压力和上游库底泥沙淤积压力;
双曲拱坝温度场计算所需温度边界包括坝体外部与空气接触部分温度边界,水位以下坝面温度边界,河谷、河岸在水位以下部分温度边界,河谷、河岸在水位以上部分温度边界,地基深部温度边界。
6.如权利要求5所述的一种基于内聚力算法的双曲拱坝坝基渗流评估方法,其特征在于,
坝体外部与空气接触部分温度边界采用传热学第三类边界,在接触面产生对流换热;
水位以下坝面温度边界等效为第一类边界条件,相应温度为实测水温;
河谷、河岸在水位以下部分温度边界,按照第一类边界条件处理,取值为相应高程实测水温;
河谷、河岸在水位以上部分温度边界,将其转化为第一类边界条件,具体取值为每旬平均气温;
地基深部温度边界受地温影响,地基侧面按绝热边界处理;
上游库底泥沙淤积压力根据泥沙淤积压力标准值、淤积泥沙的浮容重、坝前泥沙淤积厚度和淤积泥沙内摩擦角进行计算。
7.如权利要求1所述的一种基于内聚力算法的双曲拱坝坝基渗流评估方法,其特征在于,
基于内聚力算法的热渗流应力耦合有限元模型采用牛顿-拉夫森迭代方法进行求解,其包括温度场的计算、应力场的计算、渗流场的计算:
所述温度场采用热对流方程进行计算;
所述应力场采用内聚力单元内部的基于立方定律建立的平衡方程进行计算;
所述渗流场采用菲克扩散定律进行计算。
8.如权利要求1-7任一所述的一种基于内聚力算法的双曲拱坝坝基渗流评估方法,其特征在于,
在双曲拱坝的结构安全状态的评估前,将双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况与已有的渗流监测数据,进行对比分析,以确定基于内聚力算法的热渗流应力耦合有限元模型仿真计算是否准确;
所述对比分析包括对不同测点的渗流量数据进行对比和残差分析。
9.一种基于内聚力算法的双曲拱坝坝基渗流评估方法,其特征在于,
包括以下步骤:
步骤1:根据双曲拱坝体型参数,生成三维几何模型,构建基于内聚力算法的热渗流应力耦合有限元模型;
所述基于内聚力算法的热渗流应力耦合有限元模型,用于对双曲拱坝的温度场、应力场和渗流场进行耦合计算,其采用内聚力算法模拟坝体和坝基交界面处的微裂缝和孔隙,并通过带有孔隙水压力自由度的三维二次单元计算双曲拱坝坝基的渗流场;
步骤2:利用双曲拱坝监测数据生成热力学载荷数据和结构力学载荷数据;
根据热力学载荷数据和结构力学载荷数据构建步骤1中的基于内聚力算法的热渗流应力耦合有限元模型的边界条件;
步骤3:根据步骤2中的边界条件,求解基于内聚力算法的热渗流应力耦合有限元模型,得到双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况;
步骤4:利用步骤3中的双曲拱坝坝基渗流计算值以及坝基渗流场的分布情况,实现双曲拱坝的结构安全状态的评估。
10.一种基于内聚力算法的双曲拱坝坝基渗流评估系统,其特征在于,
应用如权利要求1-9任一所述的一种基于内聚力算法的双曲拱坝坝基渗流评估方法;
其包括双曲拱坝设计参数存储模块、热渗流应力耦合有限元模块、双曲拱坝数据监测模块、边界条件处理模块、安全评估模块;
所述双曲拱坝设计参数存储模块,用于存储双曲拱坝体型参数;
所述热渗流应力耦合有限元模块,用于对双曲拱坝的温度场、应力场和渗流场进行耦合计算;
所述双曲拱坝数据监测模块,用于对双曲拱坝的状态进行监测;
所述边界条件处理模块,用于获得双曲拱坝的边界条件;
所述安全评估模块,用于对双曲拱坝的结构安全状态进行评估;
所述双曲拱坝设计参数存储模块将双曲拱坝体型参数输入到热渗流应力耦合有限元模块中;
所述双曲拱坝数据监测模块将监测数据传输到边界条件处理模块中,生成边界约束和载荷数据,并将边界约束和载荷数据输入到热渗流应力耦合有限元模块中;
所述热渗流应力耦合有限元模块根据双曲拱坝体型参数以及边界约束和载荷数据,对双曲拱坝的温度场、应力场和渗流场进行耦合计算,并将计算结果传输给安全评估模块;
所述安全评估模块根据热渗流应力耦合有限元模块的计算结果对双曲拱坝的结构安全状态进行评估。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211071029.XA CN115146565B (zh) | 2022-09-02 | 2022-09-02 | 一种基于内聚力算法的双曲拱坝坝基渗流评估方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211071029.XA CN115146565B (zh) | 2022-09-02 | 2022-09-02 | 一种基于内聚力算法的双曲拱坝坝基渗流评估方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115146565A true CN115146565A (zh) | 2022-10-04 |
CN115146565B CN115146565B (zh) | 2022-11-25 |
Family
ID=83415942
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211071029.XA Active CN115146565B (zh) | 2022-09-02 | 2022-09-02 | 一种基于内聚力算法的双曲拱坝坝基渗流评估方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115146565B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104846779A (zh) * | 2015-05-22 | 2015-08-19 | 辽宁工程技术大学 | 一种堆筑不透水尾矿库坝体的方法 |
CN105354394A (zh) * | 2015-12-03 | 2016-02-24 | 中国电建集团西北勘测设计研究院有限公司 | 一种基于三维可视化的拱坝坝肩边坡稳定判断方法 |
CN106094011A (zh) * | 2016-06-30 | 2016-11-09 | 马克 | 双曲拱坝微震监测系统及方法 |
US20170039307A1 (en) * | 2015-08-07 | 2017-02-09 | Qrisq Analytics, LLC | Large scale analysis of catastrophic weather damage |
CN110673674A (zh) * | 2019-10-22 | 2020-01-10 | 浙江浙能华光潭水力发电有限公司 | 一种小型水电站全感知交互系统及方法 |
CN114330052A (zh) * | 2021-12-08 | 2022-04-12 | 中国水利水电科学研究院 | 一种谷幅收缩作用下拱坝真实变形特性分析方法 |
CN114547951A (zh) * | 2022-04-24 | 2022-05-27 | 浙江远算科技有限公司 | 一种基于数据同化的大坝状态预测方法及系统 |
-
2022
- 2022-09-02 CN CN202211071029.XA patent/CN115146565B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104846779A (zh) * | 2015-05-22 | 2015-08-19 | 辽宁工程技术大学 | 一种堆筑不透水尾矿库坝体的方法 |
US20170039307A1 (en) * | 2015-08-07 | 2017-02-09 | Qrisq Analytics, LLC | Large scale analysis of catastrophic weather damage |
CN105354394A (zh) * | 2015-12-03 | 2016-02-24 | 中国电建集团西北勘测设计研究院有限公司 | 一种基于三维可视化的拱坝坝肩边坡稳定判断方法 |
CN106094011A (zh) * | 2016-06-30 | 2016-11-09 | 马克 | 双曲拱坝微震监测系统及方法 |
CN110673674A (zh) * | 2019-10-22 | 2020-01-10 | 浙江浙能华光潭水力发电有限公司 | 一种小型水电站全感知交互系统及方法 |
CN114330052A (zh) * | 2021-12-08 | 2022-04-12 | 中国水利水电科学研究院 | 一种谷幅收缩作用下拱坝真实变形特性分析方法 |
CN114547951A (zh) * | 2022-04-24 | 2022-05-27 | 浙江远算科技有限公司 | 一种基于数据同化的大坝状态预测方法及系统 |
Non-Patent Citations (2)
Title |
---|
CHAOSHENG TANG等: "Influencing factors of geometrical structure of surface shrinkage cracks in clayey soils", 《ENGINEERING GEOLOGY 》 * |
赵俊宇: "库水位反复升降对均质土坝渗透稳定性影响研究", 《中国优秀硕士学位论文全文数据库(电子期刊)》 * |
Also Published As
Publication number | Publication date |
---|---|
CN115146565B (zh) | 2022-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Coupled thermo-hydro-mechanical modelling for geothermal doublet system with 3D fractal fracture | |
Wang et al. | A three-dimensional thermo-hydro-mechanical coupled model for enhanced geothermal systems (EGS) embedded with discrete fracture networks | |
Raymond et al. | Numerical analysis of thermal response tests with a groundwater flow and heat transfer model | |
Yao et al. | Effects of fracture density, roughness, and percolation of fracture network on heat-flow coupling in hot rock masses with embedded three-dimensional fracture network | |
Wu et al. | A coupled thermo-hydro-mechanical model for evaluating air leakage from an unlined compressed air energy storage cavern | |
Liu et al. | Numerical modeling of local scour and forces for submarine pipeline under surface waves | |
Zhou et al. | Flow–stress coupled permeability tensor for fractured rock masses | |
Zhou et al. | Numerically simulating the thermal behaviors in groundwater wells of groundwater heat pump | |
Kazemzadeh-Parsi et al. | Three dimensional smoothed fixed grid finite element method for the solution of unconfined seepage problems | |
Zhou et al. | An iterative method for evaluating air leakage from unlined compressed air energy storage (CAES) caverns | |
Li et al. | A numerical procedure for unsaturated seepage analysis in rock mass containing fracture networks and drainage holes | |
CN106202980A (zh) | 一种膨胀土在降雨入渗条件下增湿膨胀数值模拟方法 | |
Xue et al. | A fast numerical method and optimization of 3D discrete fracture network considering fracture aperture heterogeneity | |
Xiao et al. | Modelling of fractured horizontal wells with complex fracture network in natural gas hydrate reservoirs | |
CN112012712A (zh) | 一种嵌入式离散裂缝的注水生长缝数值模拟方法和装置 | |
Li et al. | Thermal deformation coordination analysis of CC-RCC combined dam structure during construction and operation periods | |
Chen et al. | Convective heat transfer of water flow in intersected rock fractures for enhanced geothermal extraction | |
Yan et al. | A three-dimensional thermal-hydro-mechanical coupling model for simulation of fracturing driven by multiphysics | |
Shi et al. | A modified embedded discrete-fracture model to study oil-water two-phase heat and mass transfer in the complex fracture network | |
Wang et al. | Influences of reservoir heterogeneity and anisotropy on CO 2 sequestration and heat extraction for CO 2-based enhanced geothermal system | |
Huo et al. | Porosity and permeability variations of a dam curtain during dissolution | |
Sha et al. | Modeling of Hydraulic Fracture of Concrete Gravity Dams by Stress‐Seepage‐Damage Coupling Model | |
CN115146565B (zh) | 一种基于内聚力算法的双曲拱坝坝基渗流评估方法及系统 | |
Nan et al. | Inversion of the Seepage Parameters of Earth/Rockfill Dams Considering the Coupling Effect of Seepage and Thermal Transfer | |
Sun et al. | Simulation of the nonplanar three-dimensional thermal cracking using the finite element-meshfree method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |