基于热-力-电磁场网络统一的实现方法
【技术领域】
[0001] 本发明涉及热力电磁分析领域，具体地，涉及基于热-力-电磁场网络统一的实现 方法。
【背景技术】
[0002] 目前市面上求解电磁感应电流、电压的电磁软件，如ANS0FT;计算流体及热变形 的软件，如MP;它们都是基于不同的网络划分处理模型，前者采用大多基于六面体网格处 理模型，后者采用大多基于三角形网格处理模型。因此，两者无法做到网格以及文件格式统 一，这也就意味着无法实现热力场分布位置与电磁场分布位置关系的一一对应，右侧带来 模型仿真分析结果不够准确、分析判断过程时间长等问题。
[0003]由于热力分析软件所解决和分析的问题是基于由三角形网格构成的四面体模型， 采用该网格可较为准确的描绘模型的结构尺寸及外观，且较为准确的实现热辐射计算，但 在电磁场辐射下其模型存在发散不收敛的情况，导致计算迭代次数很长，需要高性能计算 机处理，势必带来硬件成本增加。而电磁分析软件通常采用的是基于六面体网格处理模型， 其电磁场辐射收敛性较高且计算量小、效率高，但它在热力分析模型中不能描述非规则几 何边缘效应所引起的弹射，进而带来热辐射计算不准确的问题。
[0004] 因此，如何对一个分析模型进行合理的有限元网格划分，即根据不同的模型能自 动找到一种计算量小、热辐射与电磁场计算都比较准确的网格划分方法。最终形成热力场 与电磁场统一的网格模型以及如何实现热力场与电磁场分析软件之间的网格文件格式转 换是本次发明所要解决的核心问题。
【发明内容】
[0005] 针对现有技术中的缺陷，本发明的目的是提供一种基于热-力-电磁场网络统一 的实现方法。
[0006] 根据本发明提供的基于热-力-电磁场网络统一的实现方法，包括如下步骤：
[0007] 步骤1 :利用热力分析软件对模型进行有限元网格划分，从而得到模型的热力耦 合网格；
[0008] 步骤2 :将模型的热力耦合网格文件格式转换成电磁软件能够识别的文件格式；
[0009] 步骤3:将转换后的热力耦合网格文件输入到电磁软件中，分析得到模型的热力 场分布数据，从而得到基于温度场分布及变形的电磁场数据分布图。
[0010] 优选地，所述步骤1包括：
[0011] 步骤1. 1 :开启热力分析软件，设置模块的存储路径，启动网格模块；
[0012] 步骤1. 2 :选择几何模型，设置该模型的网格划分条件和边界条件；
[0013] 步骤1. 3:对所述模型结构进行分析，对模型内部结构采用六面体网格划分，对模 型的边缘结构采用棱柱网格划分；
[0014] 步骤1. 4:通过热力分析软件中的程序针对模型进行热力耦合计算，得到模型的 热力耦合网格。
[0015] 优选地，所述步骤L3包括：
[0016] 步骤1. 3. 1 :针对模型的结构，将模型内部划分成用笛卡尔坐标值所表示的子区 域，并从Fluidyn求解模型中生成模型内部的初始六面体网格；
[0017] 步骤1. 3. 2 :利用累计弦长参数化法对初始六面体网格的每条边进行参数化；计 算公式如下：
[0019] 式中：Γ?η?ι表示最小边长长度，ΛXm表示X方向上增量，Λym表示y方向上增量， Λzm表示z方向上增量；
[0020] 步骤1. 3. 3 :利用拉格朗日插值求解模型的边界函数；计算公式如下：
[0022] 式中：f(r)表示模型的边界函数，r表示网格边长，r。表示边网格边长初始设定 值，A表示按增量递增第i次数值，rn表示按增量递增第η次数值；
[0023] 步骤1. 3. 4 :根据得到的模型的边界函数，运用孔斯线性混和插值公式求得所述 模型内部子区域、几何边界及内部边界的笛卡尔坐标值采用无量纲大于1的rst表示三维 空间方向；计算公式如下：
[0024] x(r,s,t) = (1-s) (1-t)fi(r) + (1-s)tf2 (r)+stf3 (r)+s(1-t)f4 (r)
[0025] + (1-r) (1-t)f5 (s) + (l~r)tf6 (s) +rtf7 (s)+r(1-t)f8 (s) + (1-r) (1-s)f9 (t)
[0026] + (1-r)sf10 (t) +rsfn (t)+r(1-s)f12 (r)+c(r,s,t)
[0027] 式中：x(r,s,t)表示曲面方程，r表示空间r方向，s表示s方向，t表示t方向， f\(r)表示积累参数弦长边界1的值，f2(r)积累参数弦长边界2的值，以此类推，c(r，s，t) 表示表示其他曲面的曲线方程；
[0028] 其中，
[0029] c(r,s,t) = -2[f(l~r) (l~s) (l-t)x(0,0,0) + (1-r) (l-s)tx(O, 0, 1) + (1-r) s(l~t)x(0, 1, 0)
[0030] +(l_r)stx(0, 1,l)+r(l-s)x(l，0，0) + (r(l-s)tx(l，0,l)+rs(l_t) x(1, 1,0)+rstx(1, 1, 1)]
[0031] 式中：c(r,s,t)表示；同理求得x(r,s,t)、y(r,s,t)、z(r,s,t)的值；
[0032] 步骤 1. 3. 5 :利用步骤 1. 3. 4 同理求得x(r,s,t)、y(r,s,t)、z(r,s,t)的值，将坐 标值所在的节点相连，得到模型新划分的有限元网格及网格坐标。
[0033] 优选地，所述步骤2包括：
[0034] 步骤2. 1 :以模型新划分的有限元网格为基础，确定求解域的物理性质和几何区 域；
[0035] 步骤2. 2 :实现求解域的离散化，并确定状态变量边界条件；
[0036] 步骤2. 3 :在新划分的网格节点上增加力学自由度数，即通过单元刚度矩阵数学 计算方法求得各节点的自由度数；
[0037] 步骤2. 4 :通过热力网格特征计算，并采用参数映射方法进行网格对应，得到电磁 分析软件能够调用的数据格式。
[0038] 优选地，所述步骤2. 3中的单元刚度矩阵的计算公式如下：
[0039] [K] =JΩ(e) [B]T [C] [B]det[J]Dω;
[0040] 式中：[Κ]表示刚度；[Β]表示应变位移矩阵，所述矩阵由单元形状函数偏导数构 成；[C]表示应力应变矩阵，所述由材料特性参数构成；[J]表示Jacobi变换矩阵，即将各 单元的自然坐标变换为相同的母单元局部坐标；Ω(ε)表示单元域，[Β]τ表示矩阵的转置； det[J]表示扭转力，D表示弹性力，ω表示矫正参赛；
[0041] 应用Gauss数值求积分法，得到如下计算公式：
[0043] 式中：1?]表示刚度矩阵，心表示刚度，w⑴表示矫正参赛，det[J]⑴表示变换 矩阵，[1?表示几何矩阵转置，[(：](1)表示应变矩阵，^^(^表示几何矩阵；其中i，j= 1，2,…，m,m是单元的节点数，1 = 1，2,…，n，为η个数值积分点的权系数；
[0044] 则对一个线性结构，其单元刚度矩阵就具有如下所示的分块形式如下：
[0045]
[0046] 式中：[K(e)]表示刚度矩阵值，[KJ表示mn段刚度矩阵；即通过计算矩阵的下三 角部分，即式中的i多j部分，求解得到电磁分析软件能够调用的数据格式。
[0047] 与现有技术相比，本发明具有如下的有益效果：
[0048] 1、本发明提供的方法是基于温度场分布及变形的电磁场数据分布图，采用此电磁 场仿真手段比传统独立热-力-磁场分析方法能够更加精确、高效。
[0049] 2、本发明提供的方法能够对电磁场、电路和系统全集成的设计环境做到精确的细 节考虑，确保了系统性能，减低设计风险。
[0050] 3、本发明提供的方法中通过对被分析结构单元采取内部结构与边缘结构不同的 网格划分方法，即内部结构采用六面体网格划分，结构单元边缘采用棱柱网格划分，能够有 效提高计算精度、减少计算规模。
【附图说明】
[0051] 通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、 目的和优点将会变得更明显：
[0052] 图1为本发明提供的基于热-力-电磁场网络统一的实现方法的流程示意图；
[0053] 图2为本发明中应用的笛卡尔坐标系示意图。
【具体实施方式】
[0054] 下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术 人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术 人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明 的保护范围。
[0055] 本发明基于机箱热效应计算平台工作环境下实现的一种统一的热力耦合网格 及电磁场网格的方法，并最终得出基于温度场分布及变形的电磁场数据分布图，帮助在 多物理场条件下相关产品设计、制造中的问题预测以及完成准确的仿真分析，真正实现 热-电-力-磁多物理场数值计算模拟分析，更快更准确的指导产品结构、电磁场、电路设 计和系统全集成化的设计环境能够在设计时精确考虑细节的电磁场效应，同时对已有产品 的电磁兼容性、环境适应性等问题能够精确定位。
[0056] 根据本发明提供的基于热-力-电磁场网络统一的实现方法，包括如下步骤：
[0057] 步骤1 :利用热力分析软件对模型进行有限元网格划分，从而得到模型的热力耦 合网格；
[0058] 步骤2 :将模型的热力耦合网格文件格式转换成电磁软件能够识别的文件格式；
[0059] 步骤3:将转换后的热力耦合网格文件输入到电磁软件中，分析得到模型的热力 场分布数据，从而得到基于温度场分布及变形的电磁场数据分布图。
[0060] 优选地，所述步骤1包括：
[0061] 步骤1. 1:开启热力分析软件，设置模块的存储路径，启动网格模块；
[0062] 步骤1. 2 :选择几何模型，设置该模型的网格划分条件和边界条件；
[0063] 步骤1. 3 :对所述模型结构进行分析，对模型内部结构采用六面体网格划分，对模 型的边缘结构采用棱柱网格划分；
[0064] 步骤1. 4 :通过热力分析软件中的程序针对模型进行热力耦合计算，得到模型的 热力耦合网格。
[0065] 所述步骤1. 3包括：
[0066] 步骤1. 3. 1 :针对模型的结构，将模型内部划分成用笛卡尔坐标值所表示的子区 域，并从Fluidyn求解模型中生成模型内部的初始六面体网格；
[0067] 步骤1. 3. 2 :利用累计弦长参数化法对初始六面体网格的每条边进行参数化；计 算公式如下：
[0069] 式中：Γ?η?ι表示最小边长长度，ΛXm表示X方向上增量，Λym表示y方向上增量， Λzm表示z方向上增量；
[0070] 步骤1. 3. 3 :利用拉格朗日插值求解模型的边界函数；计算公式如下：
[0072] 式中：f(r)表示模型的边界函数，r表示网格边长，r。表示边网格边长初始设定 值，A表示按增量递增第i次数值，rn表示按增量递增第η次数值；
[0073] 步骤1. 3. 4 :根据得到的模型的边界函数，运用孔斯线性混和插值公式求得所述 模型内部子区域、几何边界及内部边界的笛卡尔坐标值采用无量纲大于1的rst表示三维 空间方向；计算公式如下：
[0074] x(r,s,t) = (1-s) (1-t)fi(r) + (1-s)tf2 (r)+stf3 (r)+s(1-t)f4 (r