裂隙岩体不稳定温度场的复合单元算法研究

基于复合单元法，结合三维热传导-对流方程和“充填模型”，提出了裂隙岩体不稳定温度场的复合单元模型。该模型前处理简便快捷，计算网格生成时无需考虑裂隙的存在，网格剖分不受限制，随后利用复合单元前处理程序，依据裂隙的位置和方位将其自动离散在单元内。对常规热传导-对流方程进行自伴随性调整，应用变分原理，推导出裂隙岩体不稳定温度场的复合单元算法，该算法可分别计算出岩块子单元和裂隙的温度值，且可真实反映裂隙中水流与相邻岩块间的热能量交换规律。将复合单元数值模型计算的不稳定温度场结果与相应的实测数据进行对比分析可知，数值计算结果与实测数据基本一致，验证了裂隙岩体不稳定温度场复合单元算法的可靠性与有效性。算例分析表明，裂隙中水流与相邻岩块间有明显的热传导和热对流作用。