[发明专利]基于改进三重互易边界元法的瞬态温度计算方法

摘要

本发明公开了一种基于改进三重互易边界元法的瞬态温度计算方法，包括以下步骤：采用倍数时刻方案确定计算时刻序列；基于Laplace变换得到瞬态热传导问题的频域方程，并基于分析时刻t和Laplace数值逆变换方法确定频域方程的频域参数；基于修正亥姆霍兹基本解得到频域方程的边界‑域积分方程；基于改进三重互易公式将边界‑域积分方程包含的域积分转化为等效边界积分，得到纯边界积分方程；基于边界元法理论表面网格与域内插值点，求解频域边界积分方程；对频域解进行Laplace逆变换，求解时域温度分布。本发明提出了低自由度表达的改进三重互易法计算公式，大大地减少了三重互易法的计算时间，同时节约了存储空间。

主权项

1.一种基于改进三重互易边界元法的瞬态温度计算方法，包括以下步骤：S1：基于计算时长TT和最小时间步长Δt，采用倍数时刻方案确定计算时刻序列；S2：基于Laplace变换得到瞬态热传导问题的频域方程，并基于分析时刻t和Laplace数值逆变换方法确定频域方程的频域参数；S3：基于修正亥姆霍兹基本解得到频域方程的边界‑域积分方程；S4：基于改进三重互易公式将边界‑域积分方程包含的域积分转化为等效边界积分，得到纯边界积分方程；S5：基于边界元法理论表面网格与域内插值点，求解频域边界积分方程；S6：对频域边界积分方程的解进行Laplace逆变换，求解时域温度分布，得到计算时刻序列中所有时刻的温度值；S7：判断所得所有温度值中相邻时刻的最大温度变化值是否小于Δumax，Δumax为相邻时刻间允许最大温度变化值，若否，则返回步骤S1，基于倍数时刻方案进行新的计算时刻插入，重复S2‑S6步骤计算新插入时刻的温度值，直至满足所有相邻时刻的最大温度变化值小于Δumax。