考虑Gauss移动热源瞬态效应的传热结构拓扑优化.pdf

应用数学和力学编委会，：考虑 移动热源瞬态效应的传热结构拓扑优化周崇伟，赵清海，陈建良，时高松（青岛大学 机电工程学院，山东 青岛；青岛大学 电动汽车智能化动力集成技术国家地方联合工程研究中心，山东 青岛）摘要：针对热源位置随时间发生变化的结构热传导问题，考虑 移动热源进行瞬态热传导拓扑优化设计 分别以整个时间历程内传热结构散热弱度最小化与区域温度最大值最小化为设计目标，体积分数为约束条件，采用伴随变量法推导目标函数与约束条件的敏度信息，借助移动渐进线法更新设计变量，研究了不同 热源移动路径与移动速度对拓扑优化结果的影响 结果表明，瞬态拓扑结构相较于稳态结果具有明显时变性，同时最佳传热构型受到热源加热时间和移动速度及路径的多重影响 关 键 词：瞬态效应；移动热源；拓扑优化；散热弱度；传热结构中图分类号：；文献标志码：，（，；，）：，应用数学和力学 卷 期 年 月 ，收稿日期：；修订日期：基金项目：国家自然科学基金（）作者简介：周崇伟（），男，硕士生（：）；赵清海（），男，副教授，博士（通讯作者：）；陈建良（），男，硕士（：）；时高松（），男，硕士（：）引用格式：周崇伟，赵清海，陈建良，时高松 考虑 移动热源瞬态效应的传热结构拓扑优化 应用数学和力学，（）：；引 言移动热源散热问题在工程应用中普遍存在，如汽车制动盘摩擦生热和金属材料的热处理等 考虑移动热源条件下传热结构的瞬态热效应，设计高效散热拓扑结构，将移动热源的热量迅速地传递到环境中，从而有效提高结构的散热性能，具有重要研究意义 分析了移动热源在半无限大物体中引起的热传导问题，无限时移动热源周围的温度场表现为稳态 等考虑热传导系数因素，在较宽的热源速度范围内研究了运动热源所产生的温度场，给出了不同坐标系中，不同速度情况下，温度场随时间变化的情形，并给出了不同热源速度下的温度场分布 对于稳态热传导问题，左孔天等通过对热传导结构进行数学建模，设计出结构优化的最佳散热路径 对拓扑优化方法系统进行了比较，并对拓扑优化方法在热传导问题中的实际工程应用现状进行了综述 陈拥平等以变密度法为基础，建立了以结构散热弱度最小化为目标的稳态热传导拓扑优化模型 赵清海等提出了一种周期性约束下的稳态热传导拓扑优化设计方法 张晖等针对热源随结构变化的稳态热传导结构优化问题，考虑热载荷的拓扑相关性进行了优化设计 龙凯等基于独立联系映射法，解决了多工况下的稳态热传导拓扑优化问题 闫浩等提出了一种密度指数函数插值的有序多材料性能近似模型，实现了三类热边界条件下复杂热环境的多材料传热结构拓扑优化 等基于无网格的广义有限差分法和具有惩罚插值模型的固体各向同性微观结构，提出了稳态散热结构的拓扑优化方法 上述文献都是基于稳态热传导开展研究的，然而移动热源温度场本质上为瞬态，其温度分布大都随时间和空间变化 因此，将 方法扩展应用到瞬态拓扑优化，考虑了瞬态效应对拓扑优化结果的影响 吴书豪等基于区域温度控制函数构建瞬态热传导拓扑优化数学模型，开展了瞬态热传导结构拓扑优化方法研究 李信卿等搭建了基于有序有理近似材料属性法的多材料插值模型，提出了一种周期性多材料瞬态热传导拓扑优化设计方法 等对瞬态热传导设计的多重温度准则和三相拓扑优化进行了探究，考虑了瞬态效应的热传导拓扑优化过程中热柔顺度最大值最小化的问题 上述考虑瞬态效应的拓扑优化设计，极大地推动了考虑热源时变下传热结构设计优化理论的发展 综上可知，散热结构拓扑优化设计缺乏考虑移动热源瞬态效应对拓扑优化结果的影响 考虑移动热源的瞬态效应时，热源随时空变化增加了传热结构优化的复杂性 因此，本文基于瞬态热传导优化设计问题，建立了以变密度法为基础的拓扑优化模型 针对 热源模型进行仿真模拟，分别考虑了以传热结构区域温度最大值最小化与最小散热弱度为设计目标，通过伴随变量法对目标函数与约束条件进行了灵敏度推导，并基于移动渐近线法对设计变量进行了更新迭代 最后通过数值算例验证了所提方法的必要性与有效性 瞬态热传导效应分析热源的瞬态效应涉及到结构热传导随时间变化，需要建立复杂的瞬态热传导模型 这要求对瞬态热传导过程进行准确建模，包括考虑瞬态热传导方程、边界条件和初始条件等 尤其在移动热源进行传热结构优化时考虑热源瞬态效应，其随时间空间均发生变化 考虑瞬态响应的传热结构考虑瞬态效应的传热结构如图 所示，由热源（）加热的设计域应用了两种边界条件，一个是恒温边界，另一个是绝热边界，设计域内部采用高热传导属性的材料填充，设计域需要离散成有限数量的单元，非设计域为低导热系数基板结构 采用有限元法，温度场的瞬态热传导数学模型可描述为 ，（）式中 和 分别为热容矩阵和热传导矩阵；为节点温度向量；为移动热源的热载荷向量 表示节点温度向量 对时间 的导数，其可表示为第 期 周崇伟，等：考虑 移动热源瞬态效应的传热结构拓扑优化 （）（）利用 时间差分法求解结构温度和载荷向量：（），（）（）（）联立式（）（）并代入式（）中，得到温度场迭代公式：（）（）（），（）式中 为算法参数，时，式中向前差分，时，式中向后差分；为第 时刻的热载荷向量 图 瞬态效应传热材料结构 热源模型分析工程中常用的热源模型主要有 热源、的解析热源模型、椭球形热源及多热源组合的复合热源模型 热源具有中心温度高而周围温度递减的特点，例如制动器摩擦生热及金属热处理等，温度分布符合正态分布，模型如图 所示 图 热源模型 图 中，为热源中心，为热源有效加热半径 热源加热过程中任意一点的表面热流密度可确定为 （），（）式中 为外部热源有效加热范围内半径为 处的表面热流大小；为加热中心处的最大热流，在 热源有效范围内 ；为热源的集中系数，取值为 基于有限元思想，将设计域划分为网格状，将 热源施加到单元节点，计算热源作用区域中所有单元节点的热流密度，具体步骤如下：确定热源中心节点的位置 计算单元各节点到热源中心的距离 将 与热源有效加热半径 比较，若 不大于（如图 中 四个节点），则按照 模型计算该表面的热流密度；若 大于（如图 中节点），则热流密度赋为 应 用 数 学 和 力 学 年 第 卷 判断热源在移动过程中满足规律的节点并赋予符合 热源规律的热流密度 如图 所示，热源在移动过程中，热源作用在基于有限元网格划分后的节点上，相比正态分布的 热源，点热源作用在单个节点上 将移动热源的路径在时间和空间上进行提前规划，使其在单位时间上从指定节点间依次移动 使热源从设计域网格顶点 到达顶点，分别设计对角线移动路径以及区域 形移动路径 图 热源热流密度及模拟路径 经典 解析热源模型，其主要被简化为集中热源模型，按照加热材料及尺寸形状等方式将其简化为点状、线状、面状三种热源模式 其热源温度高度集中，如激光加热、焊接等，热源施加在结构表面时，可以将热源看作一个点热源 点热源形成的温度场其解析式为 （），（）式中 为热源所提供的热流量；为材料热扩散率；表示到热源点的距离；为导热介质的比热容；为热源供热时间 热源仿真通过对 热源模型移动过程进行仿真，规划移动路径，观察仿真过程中热源温度分布规律 以顶点为散热边界条件，四周为绝热材料，热源运动模型为 ，（）式中 为移动路径总节点；为热源移动速度，表示单位时间经过节点的快慢；为运动时间步长，它由总时间 和时间划分段数 决定 令 热源沿对角线以每秒经过一个节点的速度移动，取 ，取移动过程中 个位置点观察热源规律，分别取热源在开始、中间位置以及最终位置 如图 所示，热源匀速缓慢移动过程中符合 分布规律，在初始位置可以明显观察到中间高、四周低的温度分布规律，移动过程中留下明显的“彗尾”状温度，到达终点后，移动路径温度向散热边界方向降低 图 热源移动过程仿真 注 为了解释图中的颜色，读者可以参考本文的电子网页版本，后同 第 期 周崇伟，等：考虑 移动热源瞬态效应的传热结构拓扑优化 构建瞬态热传导拓扑优化模型基于变密度理论，考虑移动热源瞬态热传导问题建立优化模型：，：（，），）（）（）（），（），（）式中 为设计域单元相对密度，即设计变量；为设计域内单元格总数目；为各单元体积；为设计域的体积比；为设计域的原始体积；和分别为单元相对密度的两个极限值，为防止计算过程中热传导矩阵奇异，设置 ；（，），）表示目标函数 针对目标函数进行分析，考虑如下两种优化目标：目标 散热弱度（，），）（）（）；（）目标 区域温度（，），）（），（）式中 为列向量，在目标节点处为，其余为 灵敏度推导基于伴随变量法求解散热弱度最小化的目标函数 引入 因子，构建 函数：（）（）（）（）（）（）求解新目标函数的灵敏度，对任意设计变量 求导：（）（）（）（）（）（）（）（）（）（）将关于温度的偏导数进行分部积分计算，可得（）（）（）考虑 ，进行回代，所以上式经过化简为（）（）（）（）（）（）（）忽略载荷设计相关性影响，消除温度灵敏度，则目标函数灵敏度方程为（），（）（）（）（）（），（）其中，伴随向量 计算公式为 ，（）在整个时间历程上以区域内的最高温度最小化为目标，最高温度在瞬态下非连续，导致优化问题的求解变得更加困难，因此用聚合函数 重新定义最高温度目标函数为应 用 数 学 和 力 学 年 第 卷（）（）（）；（）定义空间函数为（），（）式中 为聚合参数，当 时，表示温控区域总节点数量 结合目标函数与约束条件构建优化准则公式，定义 函数，为 因子：（）（）（）（）目标函数 对设计变量 进行灵敏度推导：（）（）（）（）（），（）其中 （）（）（）（），（）（）（）对温度偏导数进行分部积分，考虑|同时不考虑载荷设计相关性，消除温度灵敏度项 忽略热载荷与设计变量关系，目标函数为 （）（）（）（）（）（）计算 因子：（）（）（），（）其中 ，（）（）（）（）（）（）算法求解过程考虑瞬态热传导响应，基于 材料插值模型，分析移动热源结构拓扑优化设计研究，具体流程如下：定义散热结构设计区域、材料属性、边界条件、热载荷工作时间 初始化材料设计变量值、温度边界和热负荷移动路径，将设计域进行离散网格划分，获得单元热容矩阵和热传导矩阵 计算整体热容矩阵和热传导矩阵，利用向前差分法得到温度场和载荷向量 计算两个目标函数对于设计变量的导数，得到灵敏度并进行灵敏度过滤 利用移动渐进线法对设计变量进行更新 判断优化设计结果是否收敛，若收敛，则输出拓扑构型；否则返回步骤，重复步骤，直至满足收敛 第 期 周崇伟，等：考虑 移动热源瞬态效应的传热结构拓扑优化 数 值 算 例考虑 移动热源在结构中瞬态热传导的问题，针对不同 热源移动路径，给出数值算例 定义传热结构背景区域的热传导系数 与热容系数 分别为 （），（）；在结构区域中施加热源并布置高热传导材料，导热系数为 （），热容系数为 （）对角线路径算例如图 所示的结构优化设计域，设计域选取为 的正方形结构，并将其离散成 个平面四边形单元 热源中心点强度为，在固定热源设计域中，热源位于设计域的中心位置，四顶点为 点，四边是绝热边界，在不同加热时间为 ，和 下，研究热传导瞬态效应对结构散热拓扑优化的影响规律（）对角线路径设计域（）中心固定热源设计域（）（）图 结构优化设计域 移动热源沿对角线从左上顶点移动至右下顶点，设计域的余下顶点温度恒为 ，其他边界为绝热材料，过滤半径选取 图 中，基于稳态传热的拓扑优化模型只考虑热量传递，忽略了设计域温度不断被上升的情况 随着加热时间不断增加，结构整体温度场并非呈现稳态不变，其温度在逐渐升高，同时拓扑构型受热源瞬态效应的影响变化明显 当加热时间足够大时，考虑瞬态热传导的传热结构拓扑构型趋向于稳态热传导拓扑结果 图 稳态热传导与瞬态热传导的拓扑优化结果 针对不同的优化目标（目标为散热弱度，目标为区域温度）对传热结构进行优化 作为对比，将集中热源模型与 热源使其沿对角线路径移动，传热结构优化结果列于表 考虑移动集中热源和 移动热源所优化的拓扑结果沿其对角线路径对称分布，且不同优化目标所优化的散热结构材料分布有明显不同 集中热源其拓扑构型受其移动路径的影响，分布更加集中，而 热源热传导路径效果更加明显 应 用 数 学 和 力 学 年 第 卷控制 移动热源加热时间，以热源移动速度 为变量，观察传热结构拓扑构型，可以清楚得到其运动过程中受到瞬态热传导影响的变化规律，如图 所示 表 不同热源拓扑优化结果 图 稳态热传导与移动热源不同速度的优化结果 图 中考虑了移动热源瞬态效应对传热结构的优化结果的影响，随加热时间的不同，其瞬态效应明显，可与图 中固定热源的传热结构瞬态热传导进行对比 移动热源速度较为缓慢时，其传热材料的拓扑构型瞬态效应明显，沿移动路径向恒温点传热，材料聚集于热源移动路径周围 一定时间内，当热源在路径上处于高速移动时，可以将其视为固定在路径上的稳态热源 热源沿路径移动，热源速度的变化导致整体结构的温度并非稳态温度场，因此，非稳态的温度变化必须要考虑移动热源在瞬态效应下对传热结构的热传导拓扑优化设计 图 为 热源移动速度下，传热结构的目标函数随迭代步数的变化 在迭代过程中，目标函数能快速收敛到最优解，不同速度的移动 热源对于传热结构优化迭代过程会有一定影响，产生轻微波动，但整体趋势平缓，证明该方法具有较好的稳定性 考虑设计域体积对于传热结构瞬态热传导效果的影响，确定热源移动速度及时间，分别以散热弱度最小化和温度最大值最小化为目标进行讨论 将传热结构材料体积分数约束 分别设置为，研究移动热源随材料体积变化对拓扑构型的影响规律 由图 可得出其优化后的材料最高温度变化趋势 两个优化目标随体积分数变化的温度变化规律，其温度变化趋势呈现一致性，即随着传热材料的增加，移动热源下的整体结构最高温度下降明显 目标在对角第 期 周崇伟，等：考虑 移动热源瞬态效应的传热结构拓扑优化线路径随材料体积变化的最高温度小于目标，而优化目标的温度变化幅度明显比目标更迅速（）目标（）目标（）（）图 不同热源速度下的目标函数随迭代步数变化曲线 图 最高温度变化 表 材料体积占比不同下拓扑优化结果 如表 所示，随着设计域体积分数的递增，两个目标优化后的拓扑构型材料的分布具有一定的规律 以散热弱度最小化为优化目标，材料占比较小时，其传热结构不够稳定，受移动热源瞬态效应影响较大，随体积应 用 数 学 和 力 学 年 第 卷分数的增加，优化后的拓扑结构趋于稳定，材料向热源移动路径集中，同时向 点延伸 区域温度最大值最小化的优化目标，其优化后的材料分布逐渐向对角线移动路径周围集中靠拢，即从加热路径向零温度点延伸的变化，材料占比足够大时，优化后的构型边界清晰、结构合理 方形区域路径针对区域路径移动热源的结构散热优化，设计域几何尺寸同样取 ，将其离散为 的四节点方形单元 热源移动区域路径如图 所示，令其设计域四顶点为温度恒定为 的恒温材料，其余材料为隔温绝热材料 图 区域 形路径结构设计域 如图 所示，给定移动热源速度 ，在离散单元中取中心 的移动范围，以热源加热时间为变量，考虑 移动热源对传热结构的拓扑结果具有明显的瞬时性，其材料分布随着加热时间的变化边界更加清晰，结构更加合理 图 形路径移动热源的拓扑构型 移动热源在离散的结构单元内部区域进行移动，规划热源移动路径范围，使其分别在单元中心区域以，的方形单元范围内，移动速度为 经过 单元，利用奇偶性改变运动方向进行 形路径移动，控制迭代次数直至达到优化结构的传热路径最佳 优化拓扑构型列于表 根据优化结果显示，在针对目标进行优化时，同一速度 移动热源的在不同移动范围内的运动对传热结构优化结果具有一定影响，运动范围相对整体结构较小时，其热源可近似为固定点热源，其拓扑构型会受热源瞬态效应的影响，构型分枝较多，其温度分布受到移动热源运动终点的影响呈现非对称性 然而，拓扑结果呈现一定规律性，在热源移动区域较小时，进行传热优化的结构其主散热结构仍然是向角点延伸，随着热源移动区域的增大及加热时间的增长，角点延伸区域材料分布越来越多，最终呈现 形结构 随着运动范围占据整体结构的比例上升，散热材料集中分布在运动路径并向恒温点发散，其拓扑构型呈现出 形，中间温度较高区域也随之变大，最高温度点位于热源运动终点处并偏向与位置最近的恒温点散热 由区域温度最大值最小化为目标的拓扑结果可知，当移动热源小范围运动时，其拓扑构型呈现 形，温第 期 周崇伟，等：考虑 移动热源瞬态效应的传热结构拓扑优化度分布大体与拓扑构型重合，中心移动范围扩大，其温度分布不再沿其拓扑构型，受运动终点的位置影响较大 表 移动热源不同移动范围的拓扑构型 对比不同目标函数下的拓扑构型及温度云图，在相同热源移动区域下，目标温度幅值明显高于目标 移动区域较小（）时，热源作用时间短，因此热源瞬态效应明显 目标优化结构传热路径较为复杂，路径尚未连接恒温点，整体设计域温度偏高，热源最高温度高达 ；目标优化后的构型形成了更为高效的传热路径，设计域温度最大温度仅为 移动区域的增加使优化后地路径有效地连通了恒温顶点，传热材料分布更加集中，温度得到有效降低 相比目标优化后的温度分布，随着移动路径的变化，其目标整体结构最高温度更低，以区域温度最大值最小化为目标的拓扑构型相比以最小散热弱度为目标的优化，在考虑移动热源的瞬态情况下具有更佳的散热效果 结 论）本文以对角线路径为例，对 热源模型移动过程中温度分布进行模拟，验证了将 移动热源作为热载荷应用于结构优化的有效性和可行性）基于变密度法建立了考虑移动 热源瞬态效应的传热结构拓扑优化模型，分析了不同优化目标下的拓扑构型 算例表明，材料比例小于 时，不同目标的温度赋值变化幅度较大，超过后最高温度趋向稳定降低 调整散热材料体积占比，可以得到传热更佳的拓扑构型，验证了该方法的有效性）优化结果表明，基于控制变量的思想，不同 热源移动速度、移动路径范围以及材料体积占比优应 用 数 学 和 力 学 年 第 卷化的拓扑构型均会受到热源瞬态效应的影响 温度云图显示目标优化后的结构整体温度明显低于目标的优化结构，最高温度仅为目标的 左右 不同优化目标的传热结构优化结果存在一定的差异性和规律性，可为将来工程应用中考虑移动热源的传热结构拓扑优化设计研究提供理论基础 参考文献（）：，匡震邦 运动热源产生的温度场的有限元分析 应用数学和力学，（）：（，（）：（）左孔天，陈立平，张云清，等 用拓扑优化方法进行热传导散热体的结构优化设计 机械工程学报，（）：（，（）：（），：陈拥平，高亮，肖蜜 基于变密度法的散热结构拓扑优化设计 计算机集成制造系统，（）：（，（）：（）赵清海，张洪信，蒋荣超，等 考虑拓扑相关热载荷的散热结构多相材料拓扑优化设计 中国机械工程，（）：（，（）：（）张晖，刘书田，张雄 拓扑相关热载荷作用下稳态热传导结构拓扑优化 中国机械工程，（）：（，（）：（）龙凯，左正兴 稳态热传导下的连续体结构拓扑优化 中国机械工程，（）：（，（）：（）闫浩，吴晓明 基于 模型的多材料传热结构拓扑优化 航空动力学报，（）：（，（）：（），（）：，（）：吴书豪，张永存，刘书田 考虑瞬态效应的承载隔热多功能结构拓扑优化 计算力学学报，（）：（，（）：（），：李信卿，赵清海，龙凯，等 考虑瞬态效应的周期性多材料传热结构拓扑优化 航空学报，（）：（，（）：（）第 期 周崇伟，等：考虑 移动热源瞬态效应的传热结构拓扑优化 ，（）：，（）：，（）：黄飞，马永斌 移动热源作用下基于分数阶应变的三维弹性体热机响应 应用数学和力学，（）：（，（）：（）周焕林，严俊，余波，等 识别含热源瞬态热传导问题的热扩散系数 应用数学和力学，（）：（，（）：（）昌俊康，段宝岩 连续体结构拓扑优化的一种改进变密度法及其应用 计算力学学报，（）：（，（）：（）乔赫廷，张永存，刘书田 散热结构拓扑优化目标函数的讨论 中国机械工程，（）：（，（）：（），（）：黄河，高佳徐，任智彬，等 内三角管式快速蓄放热单元的肋片拓扑优化 应用数学和力学，（）：（，（）：（）应 用 数 学 和 力 学 年 第 卷