一种有限元求解轧制过程温度场的集中热容矩阵方法，通过将有限元线性方程组中热容矩阵的同行或同列元素相加代替对角线元素进行计算，在不影响计算效率的情况下有效地克服了传统有限元方法计算瞬态温度场时产生的振荡现象，保证了计算的稳定性；该方法计算的中厚板轧制过程中的钢板表面温度值与实测温度的结果对比表明，该方法保证了计算瞬态温度场时具有较高的计算精度。
【技术实现步骤摘要】
本专利技术属于轧制
，特别涉及一种有限元求解轧制过程温度场的集中热容 矩阵方法。
技术介绍
在现代工程
中，有限元法作为一种有效的数值分析方法已广泛应用于结构、力 和热等的分析过程中；近年来，许多研究人员采用有限元法分析热轧和淬火等过程中瞬态温 度场。在用有限元法分析瞬态温度场时，通常采用的方法是在空间上进行结构离散、在时间 上采用有限差分格式进行求解。该方法从初始时刻的温度场开始，每隔一个时间步长进行一 次迭代计算，进而求出各个时刻的温度场；其优点是节省计算机内存，但是往往产生时间和空间上的振荡现象，影响了计算的稳定性和精度。对于上述振荡现象，研究人员从差分格式、网格划分方式、加权变系数、迭代变步长、在实数子空间上求解和用Norsette法代替差分格式计算等多方面进行了深入的研究，得出了 一些降低或消除振荡的方法，但效果都不太理想。
技术实现思路
本专利技术的目的就是通过采用集中热容矩阵的方法克服上述采用有限元法求解瞬态温度场 时常常产生的时间和空间上的振荡现象，保证计算的稳定性，提高的计算精度。实现本专利技术巨的的技术解决方案是① 采集轧制过程数据，包括轧制参数，材料热物性参数，单元划分信息。 轧制参数轧制时间，轧件宽度，轧件厚度，初始时间，初始温度，轧件周围介质温度，时间步长。 材料热物性参数热传导系数，黑度，比热，密度。 单元划分信息宽度单元数和厚度单元数② 根据单元划分数据、轧件宽度和厚度尺寸建立有限元分析模型(见图1)，然后进行单 元节点编号、确定换热边界和计算节点坐标。单元和节点编号沿厚度方向和宽度方向逐渐增加，图1中，i为单元编号，j为节点编号，H为厚度，W为宽度。确定边界条件AB和AD边界绝热，BC和CD边界换热。以A点坐标为零，计算各节点坐标，在宽度方向和厚度方向上单元均匀划分。③根据不同轧制过程，确定边界换热系数力。在板带轧制过程中，不同轧制阶段，具有不同换热系数计算模型。热轧板带 在空冷过程中，其表面换热方式主要为辐射和自然对流，辐射系数表述为<formula>formula see original document page 5</formula> (i)其中/B 为辐射系数，cr为Stefan-Boltzman常数，cr = 5,67xl(T8r/(m2 .〖4); s为黑度系数，￡与温度的关系式为e = 0.125(771000)2 - 0.38(771000) + 1.1 。热轧板带在高压水除磷过程中，主要换热方式为强迫对流和侧面辐射，辐射计算方法同上，对流系数表达式为<formula>formula see original document page 5</formula> (2)其中w (i:/min.附为水流密度；r(《)板带表面温度。在轧制过程中，板带与轧辊之间接触换热是主要热损失方式。接触换热系数与氧化铁皮 厚度、氧化铁皮导热率和轧制压力有关。接触换热系数表达式为<formula>formula see original document page 5</formula> (3)式中凡(M尸)-轧制压力 利用有限元基本原理，计算四边形等参单元的形函数iV、 B矩阵和雅克 比矩阵J。 以二维热传导基本方程为基础，利用欧拉方程建立等效泛函，确定温度场求解的系统 方程。以热力学第一定律为依据建立无内热源强度的二维热传导的微分方程为<formula>formula see original document page 5</formula>式中[&]-温度刚度矩阵，fc]=i(k(e)]+[At)D; [A]-变温矩阵，[&]=i;k3w];常数项列式，(p} = it(e)}; {T^-温度列式；五-单元总数；上标e表示每个单元。<formula>formula see original document page 5</formula>对每个单元来说，刚度矩阵、变温矩阵和常数项可以通过下式求解:'<formula>formula see original document page 6</formula>r-瞬时温度(a:) p-材料密度(&/w3)c-材料比热 f-时间(s)h热传导系数利用欧拉方程将二维热传导问题方程(4)变为等效泛涵:<formula>formula see original document page 6</formula>(5)根据热传导问题的变分原理，对泛函式(5)求一阶偏导数并置零，得到温度求解的系统方程:(6)式中其中:<formula>formula see original document page 6</formula>yo-材料密度(紐/w3);c-材料比热(々(％.尺))；h-换热系数； 7V-形函数；,'，y'-节点编号；◎利用二点向后差分格式，将系统方程转化为瞬态温度场求解的线性方程组。将系统方 程(6)中的温度对时间偏导数表示为二点向后差分格式<formula>formula see original document page 7</formula> (7) 将时间向后差分格式(4)带入系统方程式(3)得到温度场求解的线性方程组<formula>formula see original document page 7</formula>(8)⑦集中热容矩阵式(8)中，[K」是由各单元质量阵和比热的乘积组装成的矩阵，称为协调质量热容矩阵(简 称热容矩阵)，该矩阵为n (n为节点总数)阶对称带状矩阵。由于并非只有对角元素为零， 由(1)式中各方程可以看出流入某节点的热量不仅与该节点的温度变化有关还与其周围节点 的温度变化有关，这是不合理的。对于该问题一个合理的解决办法是将单元的质量集中到各节点上去。 一种可行的做法是 将热容矩阵的同行或同列元素相加代替对角线元素，新的热容矩阵只有对角线元素有值，其 余元素均为零，如(3)式所示。<formula>formula see original document page 7</formula>集中热容矩阵后(1)式中各方程的物理意义是十分明确且合理的，即流入某节点的热量 与该节点温度的变化相平衡。⑧ 求解线性方程组，获得温度场⑨ 根据迭代次数判断程序计算是否结束，如果迭代次数大于设定迭代次数则退出程序， 否则继续迭代计算。计算流程如图2所示。本专利技术在传统算法的基础上通过采用集中热容矩阵的方法，在不影响计算效率的情况下 有效地克服了振荡现象，能够保证了计算的稳定性，提高计算精度。附图说明图1本专利技术方法有限元分析模型图，图2本专利技术方法单元划分图，图3本专利技术方法计算流程图，图4实施本专利技术方法之前温度场的空间分布图，图5实施本专利技术方法之前温度随时间变化图，图6实施本专利技术方法之后温度场的空间分布图，图7实施本专利技术方法之后温度随时间变化图。图中i为单元编号，j为节点编号，H为厚度，本文档来自技高网...
【技术保护点】
一种有限元求解轧制过程温度场的集中热容矩阵方法，其特征在于该方法包括以下步骤：（１）采集轧制过程数据，包括：轧制参数，材料热物性参数，单元划分信息轧制参数：初始时间，轧制时间，轧件宽度，轧件厚度，初始温度，轧件周围介质温度， 时间步长材料热物性参数：热传导系数，黑度，比热，密度单元划分信息：宽度单元数和厚度单元数（２）建立有限元分析模型，进行单元节点编号、确定换热边界和计算节点坐标（３）根据不同轧制过程，确定边界换热系数ｈ 热轧板带，在空冷过程中，辐射系数表述为：ＨＲ＝σ.ε.（Ｔ＋Ｔ↓［ａｉｒ］）（Ｔ↑［２］＋Ｔ↓［ａｉｒ］↑［２］）式中：ＨＲ为辐射系数，σ＝５．６７×１０↑［－８］Ｗ／（ｍ↑［２］.Ｋ↑［４］）ε＝０．１２５（Ｔ／１ ０００）↑［２］－０．３８（Ｔ／１０００）＋１．１热轧板带在高压水除磷过程中，辐射系数ＨＲ，对流系数表达式为：ＨＣ↓［Ｗ］＝１２４．７×ｗ↑［０．６６３］×１０↑［－０．００１４７（Ｔ－２７３．１６）］式中：ｗ为水流 密度，Ｔ板带表面温度在轧制过程中，接触换热系数表达式为：ＩＨＴＣ＝６９５ｐ↓［ｍ］－３４４００（Ｗ／ｍ↑［２］Ｋ）式中：ｐ↓［ｍ］－轧制压力（４）利用有限元基本原理，计算四边形等参单元的形函数Ｎ、Ｂ矩阵和雅克 比矩阵Ｊ（５）以二维热传导基本方程为基础，利用欧拉方程建立等效泛函，确定温度场求解的系统方程以热力学第一定律为依据建立无内热源强度的二维热传导的微分方程为：ｋ（＊↑［２］Ｔ／＊ｘ↑［２］＋＊↑［２］Ｔ／＊ｙ↑［２］） －ρｃ＊Ｔ／＊ｔ＝０其中：Ｔ－瞬时温度，ρ－材料密度，ｃ－材料比热，ｔ－时间ｋ－热传导系数利用欧拉方程将二维热传导问题方程（４）变为等效泛涵：Ｉ＝１／２∫∫↓［Ｓ］［ｋ［（＊Ｔ／＊ｘ）↑［２］＋（＊Ｔ ／＊ｙ）↑［２］］＋２ρｃ＊Ｔ／＊ｔＴ］ｄＳ＋１／２∫↓［ｌ］ｈ（Ｔ－Ｔ↓［∞］）ｄｌ根据热传导问题的变分原理，对泛函式求一阶偏导数并置零，得到温度求解的系统方程［Ｋ↓［Ｔ］］｛Ｔ｝＋［Ｋ↓［３］］｛＊Ｔ／＊ｔ｝＝｛ｐ｝ 式中：［Ｋ↓［Ｔ］］－温度刚度矩阵［Ｋ↓［３］］－变温矩阵｛ｐ｝－常数项列式ｋ－热传导系数ρ－材料密度ｃ－材料比热ｈ－换热系数Ｎ－...
【技术特征摘要】
1、一种有限元求解轧制过程温度场的集中热容矩阵方法，其特征在于该方法包括以下步骤(1)采集轧制过程数据，包括轧制参数，材料热物性参数，单元划分信息轧制参数初始时间，轧制时间，轧件宽度，轧件厚度，初始温度，轧件周围介质温度，时间步长材料热物性参数热传导系数，黑度，比热，密度单元划分信息宽度单元数和厚度单元数(2)建立有限元分析模型，进行单元节点编号、确定换热边界和计算节点坐标(3)根据不同轧制过程，确定边界换热系数h热轧板带，在空冷过程中，辐射系数表述为HR＝σ·ε·(T+Tair)(T2+Tair2)式中HR为辐射系数，σ＝5.67×10-8W/(m2...
【专利技术属性】
技术研发人员：刘刚，李长生，刘相华，
申请(专利权)人：东北大学，
类型：发明
国别省市：89[中国|沈阳]
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。