一种电力电子集成模块冷板液流通道的拓扑优化设计方法技术

一种电力电子集成模块冷板液流通道的拓扑优化设计方法，首先要建立有限元模型并设定生长初始点及完成自适应生长的参数初始化，然后拟定冷板液流通道的自适应生长策略，给出确定冷板液流通道生长方向及分叉的方法，最后进行冷板液流通道的逐级、迭代生长；本发明专利技术可直接获得优化的冷板液流通道布局，同时摆脱了基结构节点的束缚，所以也更接近于最优结果，既不必苛求设计者具有大量的设计经验，也不再需要重复模拟、改进、再模拟的工作，明显提高了工作效率与设计质量。

Topology optimization design method for cold plate liquid flow channel of power electronic integration module

An integrated power electronic module cold plate liquid flow channel topology optimization method, we must first establish the finite element model and set the initial growth point and complete the initialization parameters of adaptive growth, and then develop adaptive cold plate liquid flow channel growth strategy, determination method of cold plate liquid flow channel growth and bifurcation direction are given, finally step by step, iterative growth of cold plate flow channels; the invention can directly obtain the cold plate liquid flow channel layout optimization, and gets rid of the base structure of the node, so it is more close to the optimal result, not only the critical design has extensive design experience, but also no longer need to repeat the simulation, and then improved the simulation work, significantly improve the working efficiency and quality of design.

【技术实现步骤摘要】
一种电力电子集成模块冷板液流通道的拓扑优化设计方法
本专利技术属于液流通道布局设计领域，具体涉及一种电力电子集成模块冷板液流通道的拓扑优化设计方法。技术背景冷板液流通道的布局设计对大功率集成化电力电子设备提高性能有着十分重要的作用；在传统的工业生产中，设计者一般会根据自身经验并结合车间中成熟的生产工艺，设计出冷板液流通道的布局，然后利用分析软件仿真或制作样品测试，找出设计中欠缺的地方，再来更正冷板液流通道的布局。这种过于依赖设计人员经验与直觉的冷板液流通道的布局设计方法非常缺少理论依据，需要经过大量实验及修正，不仅浪费材料、达不到目标，而且设计效率低，增加了企业的成本；并且从采用的冷板液流通道的布局构形看，大多数采用的是平行或蛇形布局，虽然形式简单，易于加工制造，但是与实际的生热情况不一定匹配，难以保证布局设计的合理性。
技术实现思路
为了克服上述现有设计方式的缺点，本专利技术提供一种电力电子集成模块冷板液流通道的拓扑优化设计方法，提高了设计的效率与质量。为了达到上述目标，本专利技术采取的技术方案为：一种电力电子集成模块冷板液流通道的拓扑优化设计方法，包括以下步骤：1)冷板液流通道自适应生长的初始化：1.1)生长区域的构建：根据实际待优化的散热问题划定设计域，建立设计域的有限元模型，将该有限元模型称为基结构；根据设计域内热量的产生和传导情况，对基结构施加热载荷边界条件；1.2)生长参数的初始化：对冷板液流通道自适应生长的相关参数进行初始化，参数包括流体圆管的内外半径r1与r2，管内流体的导热系数λl，管壁自身的导热系数λs，流体散热圆管的用量体积上限Vmax，冷板液流通道的生长步长L；除以上参数外，还需对冷板液流通道的生长点进行初始化；根据给定的热载荷边界条件，设定一个或多个初始生长点；冷板液流通道由初始生长点开始生长，并且随着生长，生长点的位置不断更新；2)冷板液流通道的自适应生长策略：2.1)冷板液流通道热刚度矩阵Ke的构造：冷板液流通道中的流体与通道自身在布局设计时作为整体进行考虑，采用等效法，冷板液流通道的散热效果用固体高导热材料来代替，以此来实现冷板液流通道热刚度矩阵的构造；冷板液流通道的等效导热系数λe：以高导热材料的导热系数来描述冷板液流通道的导热能力，即高导热材料的导热系数即为冷板液流通道的等效导热系数，表达式为：其中：r1是流体圆管内壁半径；r2是流体圆管外壁半径；λs是固体管壁的导热系数；h是对流换热系数，为表征冷板液流通道流固界面处热量交换强度的参数，h的取值与管内流体的导热系数、流速及流固界面的粗糙度有关，即：h＝g(λl，vmax，R)(2)其中：λl是通道内流体的导热系数；vmax是通道横截面内流体流速的最大值；R是表征冷板液流通道内流固界面处的粗糙度的参数；假设在圆管内流动的流体遵循泊肃叶定律，则简化了通道内流体流动的控制方程；其中，Q为流量体积流率，ΔP为圆管两端的压降，η为流体的粘滞系数，L为通道的轴向长度，r为通道的半径；流场速度在轴线上取得最大值，假设液流通道轴线方向的压力梯度恒定，即为常数，则：进而冷板液流通道的热刚度矩阵为：其中，Le是流体液流通道的轴向长度；2.2)生长方向的确定：利用有限单元法计算设计域在热载荷边界条件下的温度场，稳态热传导的有限元方程表达为：KTb＝F(6)冷板液流通道与基结构作为整体进行温度场计算，因此K为包含了基结构及液流通道的整体热刚度矩阵；Tb和F分别是基结构的节点温度向量和节点载荷向量；应用连续温度场插值的数值处理方法，将液流流体通道对基结构的影响通过基结构的节点温度插值表达，突破了基结构对冷板液流通道布局的束缚,具体描述如下：设计域内任意一点的温度由基结构的节点温度插值得到，如下式所示：其中，s是设计域内任意一点的位置坐标，wi是第i个插值函数，αi是相应的权重系数，N是基结构的节点总数；利用式(7)插值得到的基结构节点的温度值应等于其真实值，即：其中，是基结构上第k个节点的位置坐标，wi是该节点的第i个插值函数，αi是相应的权重系数，N是基结构的节点总数；利用式(7)计算某一冷板液流通道单元的节点温度向量为：Tc＝[αiwi(sj)]2X1(i＝1,2…N；j＝1,2)(9)其中，sj是液流通道单元上第j个节点的位置坐标，wi是该节点的第i个插值函数，αi是相应的权重系数，N是基结构的节点总数；由式(8)、(9)得：Tc＝CTb(10)其中，矩阵C为：考虑仅有冷板液流通道的热传导控制方程的有限元形式：KeTc＝Fc(12)其中，Ke由步骤2.1)中计算得到；将式(10)代入式(12)，并在等式两边左乘矩阵CT得：CTKeCTb＝CTFc(13)由式(13)得冷板液流通道在全局坐标系下的等效热传导矩阵：Kceq＝CTKeC(14)那么，包含冷板液流通道和基结构的整体热刚度矩阵K为：K＝CTKeC+Kb(15)其中，Kb是基结构的刚度矩阵；通过连续温度场插值的数值处理方法，将冷板液流通道的节点温度表示了为基结构节点温度的加权；每个基结构节点的温度所占的权重与基结构节点和冷板液流通道节点的相对位置及所选的插值函数有关；将冷板液流通道的等效热刚度矩阵转换为全局坐标下的等效热刚度矩阵，从而实现了结构整体热刚度矩阵的组装；以结构整体的散热弱度作为确定冷板液流通道生长方向的目标函数，结构整体最小的散热弱度表明结构具有最大的散热强度与可散热性，结构整体的散热弱度表达式为：D＝TbTKTb(16)冷板液流通道的生长方向的确定方法为：以生长点为起始点分别向周围各个方向生长出固定长度L的冷板液流通道，然后计算冷板液流通道在不同方向时结构整体的散热弱度，并将散热弱度取得最小值时的生长方向确定为冷板液流通道的最终生长方向，并将确定的生长方向以变量θ标识；同时更新基结构的热刚度矩阵Kb；新生长出的冷板液流通道的均按照步骤2.3)判断其是否具有分叉能力；如果该新生长的通道具有分叉能力，则将按照步骤2.3)更新下一步生长的相关参数；否则将直接以该新生长的通道的终止点作为下一步生长的起始点；2.3)分叉参数的确定：当步生长迭代中计算每段冷板液流通道生长前后结构整体的散热弱度的减少量ΔD；如果某一生长点生长出冷板液流通道前后的散热弱度减小量大于某一阀值Dv，则这一生长点被认为具有分叉能力；否则，这一生长点被认为不具有分叉能力；具有分叉能力的生长点按照植物叶脉分主、次脉的形式分级生长，主脉在分叉点分为主脉与次脉；2.3.1)分叉中主次脉半径的确定方式：为了保证通道内的流体在分叉处的流阻最小，分叉前母支与分叉后两个子支的内径满足Murray’s法则，即r30＝r31+r32，其中r0表示分叉前母支内半径，r1表示分叉后主脉内半径；r2表示分叉后次脉内半径；设定r1与r2的比值为定值，不小于1；2.3.2)分叉中主次脉生长角度的确定：当步生长迭代中生长出的冷板液流通道生长如果其生长点被认为具有分叉能力，称为分叉点；主脉在分叉点分为主脉与次脉，已经长出的液流通道作为主脉，其通道布置方向不作调整，仍记为θ，只对其半径进行调整；以分叉点为起始点预置次脉，预置的次脉轴向长记为L’；为使得分叉处流体液流通道的流动阻力最小化，预置的次脉和主脉之间的夹角需要满足一定条件，且该夹本文档来自技高网...

【技术保护点】
一种电力电子集成模块冷板液流通道的拓扑优化设计方法，其特征在于，包括以下步骤：1)冷板液流通道自适应生长的初始化：1.1)生长区域的构建：根据实际待优化的散热问题划定设计域，建立设计域的有限元模型，将该有限元模型称为基结构；根据设计域内热量的产生和传导情况，对基结构施加热载荷边界条件；1.2)生长参数的初始化：对冷板液流通道自适应生长的相关参数进行初始化，参数包括流体圆管的内外半径r

【技术特征摘要】
1.一种电力电子集成模块冷板液流通道的拓扑优化设计方法，其特征在于，包括以下步骤：1)冷板液流通道自适应生长的初始化：1.1)生长区域的构建：根据实际待优化的散热问题划定设计域，建立设计域的有限元模型，将该有限元模型称为基结构；根据设计域内热量的产生和传导情况，对基结构施加热载荷边界条件；1.2)生长参数的初始化：对冷板液流通道自适应生长的相关参数进行初始化，参数包括流体圆管的内外半径r1与r2，管内流体的导热系数λl，管壁自身的导热系数λs，流体散热圆管的用量体积上限Vmax，冷板液流通道的生长步长L；除以上参数外，还需对冷板液流通道的生长点进行初始化；根据给定的热载荷边界条件，设定一个或多个初始生长点；冷板液流通道由初始生长点开始生长，并且随着生长，生长点的位置不断更新；2)冷板液流通道的自适应生长策略：2.1)冷板液流通道热刚度矩阵Ke的构造：冷板液流通道中的流体与通道自身在布局设计时作为整体进行考虑，采用等效法，冷板液流通道的散热效果用固体高导热材料来代替，以此来实现冷板液流通道热刚度矩阵的构造；冷板液流通道的等效导热系数λe：以高导热材料的导热系数来描述冷板液流通道的导热能力，即高导热材料的导热系数即为冷板液流通道的等效导热系数，表达式为：其中：r1是流体圆管内壁半径；r2是流体圆管外壁半径；λs是固体管壁的导热系数；h是对流换热系数，为表征冷板液流通道流固界面处热量交换强度的参数，h的取值与管内流体的导热系数、流速及流固界面的粗糙度有关，即：h＝g(λl，vmax，R)(2)其中：λl是通道内流体的导热系数；vmax是通道横截面内流体流速的最大值；R是表征冷板液流通道内流固界面处的粗糙度的参数；假设在圆管内流动的流体遵循泊肃叶定律，则简化了通道内流体流动的控制方程；其中，Q为流量体积流率，ΔP为圆管两端的压降，η为流体的粘滞系数，L为通道的轴向长度，r为通道的半径；流场速度在轴线上取得最大值，假设液流通道轴线方向的压力梯度恒定，即为常数，则：进而冷板液流通道的热刚度矩阵为：其中，Le是流体液流通道的轴向长度；2.2)生长方向的确定：利用有限单元法计算设计域在热载荷边界条件下的温度场，稳态热传导的有限元方程表达为：KTb＝F(6)冷板液流通道与基结构作为整体进行温度场计算，因此K为包含了基结构及液流通道的整体热刚度矩阵；Tb和F分别是基结构的节点温度向量和节点载荷向量；应用连续温度场插值的数值处理方法，将液流流体通道对基结构的影响通过基结构的节点温度插值表达，突破了基结构对冷板液流通道布局的束缚,具体描述如下：设计域内任意一点的温度由基结构的节点温度插值得到，如下式所示：其中，s是设计域内任意一点的位置坐标，wi是第i个插值函数，αi是相应的权重系数，N是基结构的节点总数；利用式(7)插值得到的基结构节点的温度值应等于其真实值，即：其中，是基结构上第k个节点的位置坐标，wi是该节点的第i个插值函数，αi是相应的权重系数，N是基结构的节点总数；利用式(7)计算某一冷板液流通道单元的节点温度向量为：Tc＝[αiwi(sj)]2X1(i＝1,2…N；j＝1,2)(9)其中，sj是液流通道单元上第j个节点的位置坐标，wi是该节点的第i个插值函数，αi是相应的权重系数，N是基结构的节点总数；由式(8)、(9)得：Tc＝CTb(10)其中，矩阵C为：考虑仅有冷板液流通道的热传导控制方程的有限元形式：KeTc＝Fc(12)其...

【专利技术属性】
技术研发人员：李宝童，洪军，刘國光，葛柳华，闫素娜，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：陕西,61