基于热-力-电磁场网络统一的实现方法技术

本发明专利技术提供了一种基于热-力-电磁场网络统一的实现方法，包括步骤1：利用热力分析软件对模型进行有限元网格划分，从而得到模型的热力耦合网格；步骤2：将模型的热力耦合网格文件格式转换成电磁软件能够识别的文件格式；步骤3：将转换后的热力耦合网格文件输入到电磁软件中，分析得到模型的热力场分布数据，从而得到基于温度场分布及变形的电磁场数据分布图。本发明专利技术使得热-力-电磁场网格统一，能够更加精确、高效地仿真和验证设计方案，确保了系统性能，减低设计风险。

【技术实现步骤摘要】

本专利技术涉及热力电磁分析领域，具体地，涉及基于热-力-电磁场网络统一的实现 方法。
技术介绍
目前市面上求解电磁感应电流、电压的电磁软件，如ANS0FT;计算流体及热变形 的软件，如MP;它们都是基于不同的网络划分处理模型，前者采用大多基于六面体网格处 理模型，后者采用大多基于三角形网格处理模型。因此，两者无法做到网格以及文件格式统 一，这也就意味着无法实现热力场分布位置与电磁场分布位置关系的一一对应，右侧带来 模型仿真分析结果不够准确、分析判断过程时间长等问题。由于热力分析软件所解决和分析的问题是基于由三角形网格构成的四面体模型， 采用该网格可较为准确的描绘模型的结构尺寸及外观，且较为准确的实现热辐射计算，但 在电磁场辐射下其模型存在发散不收敛的情况，导致计算迭代次数很长，需要高性能计算 机处理，势必带来硬件成本增加。而电磁分析软件通常采用的是基于六面体网格处理模型， 其电磁场辐射收敛性较高且计算量小、效率高，但它在热力分析模型中不能描述非规则几 何边缘效应所引起的弹射，进而带来热辐射计算不准确的问题。 因此，如何对一个分析模型进行合理的有限元网格划分，即根据不同的模型能自 动找到一种计算量小、热辐射与电磁场计算都比较准确的网格划分方法。最终形成热力场 与电磁场统一的网格模型以及如何实现热力场与电磁场分析软件之间的网格文件格式转 换是本次专利技术所要解决的核心问题。
技术实现思路
针对现有技术中的缺陷，本专利技术的目的是提供一种基于热-力-电磁场网络统一 的实现方法。 根据本专利技术提供的，包括如下步骤： 步骤1 :利用热力分析软件对模型进行有限元网格划分，从而得到模型的热力耦 合网格； 步骤2 :将模型的热力耦合网格文件格式转换成电磁软件能够识别的文件格式； 步骤3:将转换后的热力耦合网格文件输入到电磁软件中，分析得到模型的热力 场分布数据，从而得到基于温度场分布及变形的电磁场数据分布图。 优选地，所述步骤1包括： 步骤1. 1 :开启热力分析软件，设置模块的存储路径，启动网格模块； 步骤1. 2 :选择几何模型，设置该模型的网格划分条件和边界条件； 步骤1. 3:对所述模型结构进行分析，对模型内部结构采用六面体网格划分，对模 型的边缘结构采用棱柱网格划分； 步骤1. 4:通过热力分析软件中的程序针对模型进行热力耦合计算，得到模型的 热力耦合网格。 优选地，所述步骤L3包括： 步骤1. 3. 1 :针对模型的结构，将模型内部划分成用笛卡尔坐标值所表示的子区 域，并从Fluidyn求解模型中生成模型内部的初始六面体网格； 步骤1. 3. 2 :利用累计弦长参数化法对初始六面体网格的每条边进行参数化；计 算公式如下： 式中：Γ?η?ι表示最小边长长度，ΛXm表示X方向上增量，Λym表示y方向上增量， Λzm表示z方向上增量； 步骤1. 3. 3 :利用拉格朗日插值求解模型的边界函数；计算公式如下： 式中：f(r)表示模型的边界函数，r表示网格边长，r。表示边网格边长初始设定 值，A表示按增量递增第i次数值，rn表示按增量递增第η次数值； 步骤1. 3. 4 :根据得到的模型的边界函数，运用孔斯线性混和插值公式求得所述 模型内部子区域、几何边界及内部边界的笛卡尔坐标值采用无量纲大于1的rst表示三维 空间方向；计算公式如下： x(r,s,t) = (1-s) (1-t)fi(r) + (1-s)tf2 (r)+stf3 (r)+s(1-t)f4 (r) + (1-r) (1-t)f5 (s) + (l~r)tf6 (s) +rtf7 (s)+r(1-t)f8 (s) + (1-r) (1-s)f9 (t) + (1-r)sf10 (t) +rsfn (t)+r(1-s)f12 (r)+c(r,s,t) 式中：x(r,s,t)表示曲面方程，r表示空间r方向，s表示s方向，t表示t方向， f\(r)表示积累参数弦长边界1的值，f2(r)积累参数弦长边界2的值，以此类推，c(r，s，t) 表示表示其他曲面的曲线方程； 其中， 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) +(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)] 式中：c(r,s,t)表示；同理求得x(r,s,t)、y(r,s,t)、z(r,s,t)的值； 步骤 1. 3. 5 :利用步骤 1. 3. 4 同理求得x(r,s,t)、y(r,s,t)、z(r,s,t)的值，将坐 标值所在的节点相连，得到模型新划分的有限元网格及网格坐标。 优选地，所述步骤2包括： 步骤2. 1 :以模型新划分的有限元网格为基础，确定求解域的物理性质和几何区 域； 步骤2. 2 :实现求解域的离散化，并确定状态变量边界条件； 步骤2. 3 :在新划分的网格节点上增加力学自由度数，即通过单元刚度矩阵数学 计算方法求得各节点的自由度数； 步骤2. 4 :通过热力网格特征计算，并采用参数映射方法进行网格对应，得到电磁 分析软件能够调用的数据格式。 优选地，所述步骤2. 3中的单元刚度矩阵的计算公式如下： =JΩ(e) T detDω; 式中：表示刚度；表示应变位移矩阵，所述矩阵由单元形状函数偏导数构 成；表示应力应变矩阵，所述由材料特性参数构成；表示Jacobi变换矩阵，即将各 单元的自然坐标变换为相同的母单元局部坐标；Ω(ε)表示单元域，τ表示矩阵的转置； det表示扭转力，D表示弹性力，ω表示矫正参赛； 应用Gauss数值求积分法，得到如下计算公式： 式中：1?]表示刚度矩阵，心表示刚度，w⑴表示矫正参赛，det⑴表示变换 矩阵，(1)表示应变矩阵，^^(^表示几何矩阵；其中i，j= 1，2,…，m,m是单元的节点数，1 = 1，2,…，n，为η个数值积分点的权系数； 则对一个线性结构，其单元刚度矩阵就具有如下所示的分块形式如下： 式中：表示刚度矩阵值，[KJ表示mn段刚度矩阵；即通过计算矩阵的下三 角部分，即式中的i多j部分，求解得到电磁分析软件能够调用的数据格式。 与现有技术相比，本专利技术具有如下的有益效果： 1、本专利技术提供的方法是基于温度场分布及变形的电磁场数据分布图，采用此电磁 场仿真手段比传统独立热-力-磁场分析方法能够更加精确、高效。 2、本专利技术提供的方法能够对电磁场、电路和系统全集成的设计环境做到精确的细 节考虑，确保了系统性能，减低设计风险。 3、本专利技术提供的方法中通过对被分析结构单元采取内部结构与边缘结构不同的 网格划分方法，即内部结构采用六面体网格划分，结构单元边缘采用棱柱网格划分，能够有 效提高计算精度、减少计算规模。【附图说明】 通过阅读参照以下附图对非限制性实施例所作的详细描述，本专利技术的其它特征、 目的和优点将会变得更明显： 图1为本专利技术提供的的流程示意图； 图2为本专利技术中应用的笛卡尔坐标系示意图。【具体实施方式】 下面结合具体实施例对本专利技术进行详细说明。以下实施例将有助于本领域的本文档来自技高网...

【技术保护点】
一种基于热‑力‑电磁场网络统一的实现方法，其特征在于，包括如下步骤：步骤1：利用热力分析软件对模型进行有限元网格划分，从而得到模型的热力耦合网格；步骤2：将模型的热力耦合网格文件格式转换成电磁软件能够识别的文件格式；步骤3：将转换后的热力耦合网格文件输入到电磁软件中，分析得到模型的热力场分布数据，从而得到基于温度场分布及变形的电磁场数据分布图。

【技术特征摘要】

【专利技术属性】
技术研发人员：方明，龚东磊，陈洁，段晓峰，
申请(专利权)人：中国电子科技集团公司第三十二研究所，
类型：发明
国别省市：上海;31