基于快速精确数值重构技术的混凝土块冷却系统设计方法技术方案

本发明专利技术公开了一种基于快速精确数值重构技术的混凝土块冷却系统设计方法，将混凝土块中的水管冷却作用看作系统的负热源，考虑水泥水化热的影响，首先在所考察混凝土块的内部和边界配置若干测试点，获得测试点的温度值；建立热传导问题的等效控制方程，根据径向基函数配点法得到若干局部低阶矩阵；采用多尺度算法，计算径向基函数的形状参数；利用矩阵稀疏化算法，构造计算温度场所需的方程式；计算任意时刻任意内点的温度值；根据考察结构的内外最大温差，判断冷却系统方案的合理性。本发明专利技术简单高效，具有精度高，计算快，数学简单，程序简便的特点；切合混凝土块冷却水管设计需要数据重构技术快速、稳定、精确的特征要求。

Design method of concrete block cooling system based on fast and accurate numerical reconstruction technology

The invention discloses a design method of concrete block cooling system based on rapid and accurate numerical reconstruction. The cooling effect of water pipe in the concrete block is considered as the negative heat source of the system. Considering the influence of cement hydration heat, first of all the test points are arranged in the internal and boundary of the concrete block to obtain the temperature of the test point. The equivalent control equation of the heat conduction problem is established, and some local low order matrices are obtained according to the radial basis function matching point method. The shape parameters of the radial basis function are calculated by the multi-scale algorithm, and the equations for calculating the temperature places are constructed by using the matrix sparsity algorithm, and the temperature values of arbitrary points at any time are calculated. According to the maximum temperature difference between the inside and outside of the structure, the rationality of the cooling system scheme is judged. The invention has the characteristics of simple and efficient, high precision, fast calculation, simple mathematics and simple program. The design of concrete block cooling water pipe needs the characteristic requirements of rapid, stable and accurate data reconstruction.

【技术实现步骤摘要】
基于快速精确数值重构技术的混凝土块冷却系统设计方法
本专利技术涉及含冷却水管的混凝土结构降温问题，具体涉及一种基于快速精确数值重构技术的混凝土块冷却系统设计方法。
技术介绍
热裂纹是大型混凝土结构中经常发生的一种现象。由于混凝土在凝固过程中产生大量的水化热，加之其自身导热性能差，容易产生热聚集，使得结构内部温度迅速上升。由于混凝土结构的弹性模量不大，此时徐变较大，温升引起压应力。于此同时，如果结构外表面散热不加以控制，那么外表面的温度将快速下降，当温度梯度变得很陡时，温度应力随之增大，最终导致严重的温度裂纹。因此需要设计合理的冷却系统，如外表面保温控制、水管布置、水管口径、管内通水量、及通水温度等，使得结构内外温度满足工程要求。于是，通过数值手段来实现快速预测混凝土结构的温度场，具有很高的理论和工程价值。传统上处理温度场动态数据重构，一般有有限元方法、边界元方法、有限差分法等。但作为传统区域型网格方法，有限元在处理如复杂的混凝土结构、含复杂水管冷却系统问题的数据重构时，往往会出现网格划分困难、计算速度慢、计算荷载大等难以克服的问题。而边界元方法作为边界型网格方法，虽然在一定程度上克服了有限元方法需要划分区域网格的缺点，提高了计算效率，但由于在计算过程中需要处理大量奇异与近奇异积分，严重影响了数据的重构速度(见文献1.ChengAH-D,ChengDT.Heritageandearlyhistoryoftheboundaryelementmethod.EngineeringAnalysiswithBoundaryElements2005；29(3):268-302.)。
技术实现思路
本专利技术的目的在于提供一种简单、高效、无网格的基于径向基函数配点法的温度场动态数据重构方法，节约重构时间，提高冷却系统的方案设计效率，解决现有技术的混凝土块冷却系统温度场计算过程复杂、计算荷载大、计算时间长的问题。本专利技术采用如下技术方案，基于快速精确数值重构技术的混凝土块冷却系统设计方法，包括以下步骤：(1)在所考察混凝土块的内部和边界均匀配置若干测试点，获得测试点的温度值；(2)建立热传导问题的等效控制方程，根据径向基函数配点法得到若干局部低阶矩阵；(3)采用多尺度算法，计算径向基函数的形状参数；(4)利用矩阵稀疏化算法，构造计算温度场所需的方程式；(5)计算任意时刻任意内部的测试点的温度值；(6)根据考察混凝土结构的内部测试点与边界测试点的最大温差，判断水管冷却系统设计方案是否合理。优选地，在所述步骤(1)中，初始时刻在所考察混凝土块的内部均匀配置Nip个测试点，获得初始时刻测试点的温度值Ti,i＝1,...,Nip；在所考察混凝土块的边界配置Nbp个测试点，随着热传导的进行，获得不同时刻边界所有测试点的温度值Ti,i＝Nip+1,...,Ncp，Ncp＝Nip+Nbp，Ncp表示所有测试点的个数。优选地，在所述步骤(2)中，采用的热传导问题的等效控制方程为：其中，Ω表示所考察物质的区域，x为空间坐标，τ为坐标x所对应的时刻，φ代表混凝土块和水管间的温度传递关系，θ表示绝热温升，a表示热传导系数，T表示温度场函数，表示边界条件，T0表示初值条件，Tw表示水管通水温度，Δ表示Laplace算子，Γ表示区域边界；热传导问题所用的径向基函数为：其中，Φ(r)表示多重二次曲面MQ径向基函数，c表示形状参数，r表示两两测试点间的距离；以每个测试点为中心，取Nsp个配置点构造局部区域，再利用径向基函数配点法在每个局部构造如下低阶矩阵：其中表示局部区域s上所有测试点的温度近似值，Nsp为局部区域上的配置点数，Φs是局部径向基函数矩阵αs是待定系数为径向基函数表示权函数，表示临近测试点的配置点；温度梯度由下式求得：其中Δ为拉普拉斯算子，表示空间方向上的二阶导数。式(1)、式(3)和式(4)共同构成局部径向基函数配点法的差值矩阵。4、根据权利要求1所述的基于快速精确数值重构技术的混凝土块冷却系统设计方法，其特征在于，在所述步骤(3)中，计算径向基函数的形状参数具体为：由多尺度算法推导出如下关系式：由上式求出每个局部区域上的形状参数c，用于步骤(2)构造径向基函数差值矩阵。优选地，在所述步骤(4)中采用矩阵稀疏化算法，通过定义脉冲函数，将每个局部区域构造的低阶矩阵拓展成如下全域矩阵其中xk表示配置点，k是配置点的编号，k＝1,2,…,Ncp，Ncp表示配置点总数，Ωm表示所考察的第m个局部区域，为测试点与其临近的配置点所构造的行向量，Φm是第m个局部区域的径向基函数矩阵优选地，在所述步骤(5)中，计算在时刻τ任意点x的温度场T具体为：将公式(3)、(4)、(6)代入步骤(1)中的控制方程中，并结合隐式差分格式得出：其中，Tn+1代表当前时刻的温度场，Tn代表延迟时刻的温度场，θn+1表示当前时刻的绝热温升，φn+1表示当前时刻混凝土与水管间的导热关系式，τ表示时刻，dτ表示时间步长，Nip表示内部测试点数量，Nbp表示边界测试点数量，Ncp＝Nip+Nbp，测试点数量等于配置点数量，Mm表示重构的全域矩阵的第m行向量，Em表示单位向量，表示当前时刻测试点的温度场近似值，dτ为时间步长，公式(8)和(9)为任意时刻任意测试点上温度值的计算公式。优选地，在所述步骤(6)中，根据混凝土结构的任一内部测试点与任一外部测试点的最大温差不大于设定值，即任意时刻的温度场需满足最大温差小于设定值，判断混凝土块冷却系统设计方案为合理，否则为不合理。优选地，设定值为25℃。优选地，用于调整混凝土块冷却系统方案的参数包括：外表面保温控制、水管布置、水管口径、管内通水量和通水温度。局部径向基函数配点法，将考察区域分成若干局部区域，这样可以更好的描述某些特定区域温度梯度比较大的情况；在每个局部区域构造低阶矩阵，使得求解规模大大降低，并且通过矩阵重构，最终生成含有大量零元素的矩阵便于矩阵稀疏化；通过矩阵稀疏化技术，大幅度提高计算效率，实现快速精确运算。专利技术所达到的有益效果：本专利技术是基于快速精确数值重构技术的混凝土块冷却系统的设计方法，具有精度高，计算快，数学简单，程序简便的特点，切合混凝土块冷却水管设计需要数据重构技术快速、稳定、精确的特征要求。采用局部径向基函数配点法，将考察区域分成若干局部区域，更好地描述混凝土块冷却水管附近温度骤变的现象，实现精确求解；引入多尺度技术能够快速确定径向基函数的形状参数；采用矩阵稀疏化算法，实现混凝土块的温度场分布的快速精确求解。附图说明图1为本专利技术的基于快速精确数值重构技术的混凝土块冷却系统设计方法流程图；图2为含五根水管的混凝土块截面热传导问题边界和内部数据采集点布置示意图；图3为实施例1的四个测试点上的温度值随时间的变化曲线；图4为实施例1所提方法与有限元方法在四个测试点上的数值结果偏差曲线；图5为实施例2中τ＝1d时三个测试线上的温度值；图6为实施例2所提方法与有限元方法的数值绝对误差比较。具体实施方式下面根据附图并结合实施例对本专利技术的技术方案作进一步阐述。图1为基于快速精确数值重构技术的混凝土块冷却系统设计方法的流程图，具体步骤如下：(1)初始时刻在所考察混凝土结构的内部配置Nip个测试点，获得初始时刻这些测试点的温度本文档来自技高网...

【技术保护点】
1.一种基于快速精确数值重构技术的混凝土块冷却系统设计方法，其特征在于，包括以下步骤：(1)在所考察混凝土块的内部和边界均匀配置若干测试点，获得测试点的温度值；(2)建立热传导问题的等效控制方程，根据径向基函数配点法得到若干局部低阶矩阵；(3)采用多尺度算法，计算径向基函数的形状参数；(4)利用矩阵稀疏化算法，构造计算温度场所需的方程式；(5)计算任意时刻任意内部的测试点的温度值；(6)根据考察混凝土结构的内部测试点与边界测试点的最大温差，判断水管冷却系统设计方案是否合理。

【技术特征摘要】
1.一种基于快速精确数值重构技术的混凝土块冷却系统设计方法，其特征在于，包括以下步骤：(1)在所考察混凝土块的内部和边界均匀配置若干测试点，获得测试点的温度值；(2)建立热传导问题的等效控制方程，根据径向基函数配点法得到若干局部低阶矩阵；(3)采用多尺度算法，计算径向基函数的形状参数；(4)利用矩阵稀疏化算法，构造计算温度场所需的方程式；(5)计算任意时刻任意内部的测试点的温度值；(6)根据考察混凝土结构的内部测试点与边界测试点的最大温差，判断水管冷却系统设计方案是否合理。2.根据权利要求1所述的基于快速精确数值重构技术的混凝土块冷却系统设计方法，其特征在于，在所述步骤(1)中，初始时刻在所考察混凝土块的内部均匀配置Nip个测试点，获得初始时刻测试点的温度值Ti,i＝1,...,Nip；在所考察混凝土块的边界配置Nbp个测试点，随着热传导的进行，获得不同时刻边界所有测试点的温度值Ti,i＝Nip+1,...,Ncp，Ncp＝Nip+Nbp，Ncp表示所有测试点的个数。3.根据权利要求1所述的基于快速精确数值重构技术的混凝土块冷却系统设计方法，其特征在于，在所述步骤(2)中，采用的热传导问题的等效控制方程为：其中，Ω表示所考察物质的区域，x为空间坐标，τ为坐标x所对应的时刻，φ代表混凝土块和水管间的温度传递关系，θ表示绝热温升，a表示热传导系数，T表示温度场函数，表示边界条件，T0表示初值条件，Tw表示水管通水温度，Δ表示Laplace算子，Γ表示区域边界；热传导问题所用的径向基函数为：其中，Φ(r)表示多重二次曲面MQ径向基函数，c表示形状参数，r表示两两测试点间的距离；以每个测试点为中心，取Nsp个配置点构造局部区域，再利用径向基函数配点法在每个局部构造如下低阶矩阵：其中表示局部区域s上所有测试点的温度近似值，Nsp为局部区域上的配置点数，Φs是局部径向基函数矩阵αs是待定系数为径向基函数，表示权函数，表示临近测试点的配置点；温度梯度由下式求得：其中Δ为拉普拉斯算子，表示空间方向上的二阶导数。式(1)、式(3)和式(4)共同构成局部径向基函数配点法的差值矩阵。4.根据权利要求1所述的基于快速精确数值重构技术的混凝土块冷却系统设计方法，其特征...

【专利技术属性】
技术研发人员：林继，洪永兴，陈文，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏,32