286无单元伽辽金方法在触地结构传热中及应用.docVIP

无单元伽辽金方法在触地结构传热中的应用 湖南大学土木工程学院 李向前 陈友明☆ 摘要 本文介绍了无单元伽辽金方法的基本原理及在二维稳态触地围护结构传热中的应用。结果表明：无单元伽辽金方法建模只需要节点数据，不需要单元及单元的相关性数据，大大减少了建模初期的准备工作，计算精度高，很适合处理复杂的多维传热问题。 关键词 触地围护结构传热；无单元伽辽金方法 Abstract A meshless element-free Galerkin method (EFGM) is introduced and applied to two-dimensional steady-state ground-coupled heat transfer problems. The construction of the EFGM approximation requires only nodal data—no element structure is needed. The numerical integration of the Galerkin weak form needs only simple integration cells—no element connectivity is needed. The EFGM reduced considerably the preparation of the model. The EFGM results agree very well with the finite element method results. keywords: Ground-coupled heat transfer; Element free Galerkin method 引言 节能减排是我国的能源政策。目前建筑地上部分的围护结构保温（寒冷地区）或隔热（炎热地区）材料的性能越来越好，地上部分围护结构的节能效率越来越高，以至于建筑地下部分围护结构的节能效率直接关系到整个建筑的节能效率，因此建筑地下结构的节能问题不能再忽视了。地下结构节能效率的关键在于充分地分析地下结构的传热问题。而处理这类具有二维和三维特性并且边界条件复杂的传热问题，数值方法是最佳的选择。Davies [1] 利用二维的有限差分方法把地下传热问题整合进了建筑能耗仿真软件；Deru [2] 发展了一个二维的热湿传递有限元模型，这个模型详细地处理了土壤表面的热湿平衡和土壤结冰对触地围护结构热传递的影响；Mihalakakou 等[3] 利用瞬态的有限容积方法预测了建筑地下部分不同深度土壤的温度。以上所提到的有限差分法、有限元法、有限容积法有一个通用的特性——都是以单元为基础的数值方法。然而把复杂问题离散为一个个小单元紧密相关的线性问题是一个非常耗时而且很容易出错的过程。为了避免这个繁琐的需要考虑单元之间关联性的离散过程，本文引进了无单元伽辽金方法[4]。无单元伽辽金方法形函数的构建只需要节点数据，不需要单元数据，也不需要单元之间的相关性数据。实例表明无单元伽辽金方法计算精度高，稳定性好，建模简单，非常适合处理复杂的二维和三维传热问题。 无单元伽辽金方法 无单元伽辽金方法最主要的特点是形函数的构建是通过移动最小二乘方法来完成的。移动最小二乘方法由三个部分组成：权函数、基函数和系数。根据移动最小二乘原理，全域内温度场的近似表达式[5]为： 式中，是向量坐标 是一组基函数，m 是基函数的项数。二维空间的线性基 ；是一组随变化的系数 全域内任一点的局部近似式为 取二次函数 式中，n 是邻域内使权函数大于0的节点数, 是节点I 的参数。本文采用的是四次样条权函数 式中，r 是正则化半径。 对求极值可得 式中 把式代入式，可得移动最小二乘近似式 变分原理和方程离散 各向同性均匀介质的稳态热传导方程： 本质边界条件： 自然边界条件： 式中，k是导热系数，h是对流换热系数，是本质边界上的已知温度，是环境空气温度，n 是自然边界上的单位外法向向量。 采用拉格朗日乘子方法[5]处理本质边界条件，式-的变分为 令，得 结果 图1方形地下室的几何形状和边界条件 本文算例的几何形状如图1所示，地下室地板的半宽长度a为5m，深度b为4m，遥远边界距离F为10m，深层土壤边界深度D为12m，室内对流换热系数hi为5 W/(m2 K)，室外对流换热系数he为7 W/(m2 K)，室内空气温度Ti为22 oC，室外地表面温度Te 为26 oC，深层土壤边界温度Tdg为18 oC。地下室周围的土壤简单假定是一个层状的土壤，每一层的土壤是各向同性均匀的，一层土壤的导热系数k1为2 W/(m