《电磁场流场温度场耦合计算中的异型网格映射方法.pdf》由会员分享，可在线阅读，更多相关《电磁场流场温度场耦合计算中的异型网格映射方法.pdf（9页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 103617367 A (43)申请公布日 2014.03.05 CN 103617367 A (21)申请号 201310652717.X (22)申请日 2013.12.06 G06F 19/00(2011.01) (71)申请人 三峡大学 地址 443002 湖北省宜昌市大学路 8 号 (72)发明人 张宇娇 李杰 黄雄峰 (74)专利代理机构 宜昌市三峡专利事务所 42103 代理人 彭娅 (54) 发明名称 电磁场 - 流场 - 温度场耦合计算中的异型网 格映射方法 (57) 摘要 一种电磁场 - 流场 - 温度场耦合计算中的异 型网格映射方法， 该方法包括。
2、以下步骤 ： 1） 建立 流场 - 温度场计算模型 ； 2） 按读取各模型离散单 元的节点坐标值 ； 3） 再次建立研究对象在电磁场 计算中的几何模型， 将模型离散为三角形有限元 网格， 进行电磁场数值计算 ； 4） 将功率损耗映射 到温度场中的各单元 ； 5） 得到映射修正系数 ； 6） 将修正后的各单元功率损耗密度作为载荷加载到 温度场中进行求解计算， 得到温度结果 ； 从而得 到温度场中各单元功率损耗。本发明提供的电磁 场 - 流场 - 温度场耦合计算中的异型网格映射方 法， 可实现不同物理场采用不同类型网格之间的 数据传递， 大大减少了网格数量， 提高了计算时间 及精度。 (51)In。
3、t.Cl. 权利要求书 1 页 说明书 6 页 附图 1 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书1页 说明书6页 附图1页 (10)申请公布号 CN 103617367 A CN 103617367 A 1/1 页 2 1. 一种电磁场 - 流场 - 温度场耦合计算中的异型网格映射方法， 其特征在于该方法包 括以下步骤 ： 1） 建立流场 - 温度场计算模型， 将模型离散为四边形有限元网格 ； 2） 按逆时针顺序读取各模型离散单元的节点坐标值 ； 3） 再次建立研究对象在电磁场计算中的几何模型， 根据电磁场理论中磁场分布特点， 将模型离散为三角形有限元网格，。
4、 进行电磁场数值计算 ； 4） 将功率损耗映射到温度场中的各单元 ； 5） 将映射后的总功率与电磁场计算结果中导体总功率进行对比， 得到映射修正系数 ； 6） 将修正后的各单元功率损耗密度作为载荷加载到温度场中进行求解计算， 得到温度 结果 ； 从而得到温度场中各单元功率损耗。 2. 根据权利要求 1 所述的电磁场 - 流场 - 温度场耦合计算中的异型网格映射方法， 其 特征在于 ： 步骤 1） 中， 建立研究对象在流场 - 温度场计算中的几何模型， 并根据流体力学理 论中流体运动特点， 将模型离散为四边形有限元网格。 3. 根据权利要求 1 所述的电磁场 - 流场 - 温度场耦合计算中的异型。
5、网格映射方法， 其 特征在于 ： 步骤 4） 中， 采用高斯积分和等参单元矩阵变换方法， 将功率损耗密度映射到流 场 - 温度场模型的各单元上。 4. 根据权利要求 1 所述的电磁场 - 流场 - 温度场耦合计算中的异型网格映射方法， 其 特征在于 ： 步骤 6） 中， 通过修正系数， 对流场 - 温度场模型的各网格映射后的功率损耗进行 修正， 并作为载荷代入能量方程， 施加边界条件， 对连续性方程、 Navier-Stokes 方程及能量 守恒方程进行数值求解后， 得到温度结果。 权 利 要 求 书 CN 103617367 A 2 1/6 页 3 电磁场 - 流场 - 温度场耦合计算中的异。
6、型网格映射方法 技术领域 0001 本发明涉及一种网格映射方法， 尤其是一种电磁场 - 流场 - 温度场耦合计算中的 异型网格映射方法。 背景技术 0002 多物理场耦合系统中最突出的矛盾是各物理场对网格离散的不同特征要求和网 格间的数据映射。 传统的处理方法分为两种， 在分别进行各物理场求解时， 若根据各自特征 方式离散， 建立不同的有限元模型， 耦合时则无法直接一一对应地读取单元或节点计算结 果， 通常采用平均载荷加载方式， 然而尤其是对于时域瞬态多物理场耦合问题的实时求解 （特别是含有媒质运动、 流体流动等） ， 这种加载方式显然无法满足工程设计需要 ； 而若建立 统一的有限元模型， 则。
7、必须同时兼顾各物理场对网格形状和精度的要求， 导致网格数量极 大增加， 对计算机硬件的要求更为苛刻。 0003 对于多物理场分析的弱耦合方法， 在分别进行各物理场求解时， 若根据各自特征 方式离散， 建立不同的有限元模型， 耦合时则无法直接一一对应地读取单元或节点计算结 果 ( 如图 1 所示 )， 通常采用平均载荷加载方式， 然而尤其是对于时域瞬态多物理场耦合问 题的实时求解 （特别是含有媒质运动、 流体流动等） ， 这种加载方式显然无法满足工程设计 需要 ； 而若建立统一的有限元模型， 则必须同时兼顾各物理场对网格形状和精度的要求， 导 致网格数量极大增加， 对计算机硬件的要求更为苛刻。要。
8、更好地实现多物理场耦合分析的 弱耦合求解， 首先必须解决耦合计算中的网格数据映射问题， 鉴于目前传统的处理方法难 以满足计算精度及计算效率的要求， 多物理场耦合计算网格映射急需突破传统的方法， 寻 求新的理论及方法的指导。 0004 不匹配网格的数据传递 (Data Transfer between Non-matching Meshes) 是电 磁场数值计算中的重要问题， 它允许将求解域分解为多个子域， 在各区域以最适合子域特 征的方式离散， 而在各区域交界面上， 边界节点不要求逐点匹配 （如图 2 所示） ， 网格数据在 交界面上实现映射。近年来， 国内外研究者对不匹配网格的数据传递的研究。
9、取得了一些成 果 ： 0005 (1)1996年由德国Fraunhofer研究院SCAI中心提出， 主要应用于航空航天领域飞 机设计中流固耦合计算的不同区域不匹配网格间的数据传递 ； 0006 (2)2004 年美国伊利诺斯大学的 Xiangmin Jiao 和 Michael T.Heath 对流固耦合 问题中的不匹配网格映射从数值方法上提出了新的思路， 其中包括逐点内外插值法、 面积 加权平均法， Mortar 元法等 ； 0007 (3)2006年德国弗里德里希 亚历山大埃尔兰根纽伦堡大学的Faisal Ahmed在硕 士论文中提出了解决流固耦合问题中的不匹配网格数据映射的数值方法 ； 。
10、0008 (4)2008 年日本物理化学研究所的 Yoshifumi Okamoto 将不匹配网格的映射应用 到电机的电磁 - 热耦合分析中 ； 0009 (5)2011 年武汉大学刘守豹博士在处理运动导体涡流场数值计算时采用的 说 明 书 CN 103617367 A 3 2/6 页 4 Mortar 元法也可以解决边界节点不匹配问题 ； 0010 (6)2011 年全球最大的有限元分析软件公司之一的美国 ANSYS 在最新的版本中推 出了适用于电磁场、 流场、 温度场耦合计算中各场之间异型网格映射的功能。 0011 不匹配网格的数据传递方法和原理为多物理场耦合数值计算中异型网格映射提 供了。
11、一种新的思路， 可以推广到全模型的网格映射。因此， 需要面向电磁装置的温升分析， 探索出适用于多物理场耦合计算的异型网格映射新理论与方法。 发明内容 0012 本发明所要解决的技术问题是提供一种电磁场-流场-温度场耦合计算中的异型 网格映射方法， 提出异型网格映射方法来处理耦合场分析中整体模型的不一致网格间数据 传递问题， 解决不同物理场对网格剖分精度以及形状不同需求的问题， 可实现不同物理场 采用不同类型网格之间的数据传递， 大大减少了网格数量， 提高了计算时间及精度。 0013 为解决上述技术问题， 本发明所采用的技术方案是 ： 一种电磁场-流场-温度场耦 合计算中的异型网格映射方法， 该。
12、方法包括以下步骤 ： 0014 1） 建立流场 - 温度场计算模型， 将模型离散为四边形有限元网格 ； 0015 2） 按逆时针顺序读取各模型离散单元的节点坐标值 ； 0016 3） 再次建立研究对象在电磁场计算中的几何模型， 根据电磁场理论中磁场分布特 点， 将模型离散为三角形有限元网格， 进行电磁场数值计算 ； 0017 4） 将功率损耗映射到温度场中的各单元 ； 0018 5） 将映射后的总功率与电磁场计算结果中导体总功率进行对比， 得到映射修正系 数 ； 0019 6） 将修正后的各单元功率损耗密度作为载荷加载到温度场中进行求解计算， 得到 温度结果 ； 0020 从而得到温度场中各单。
13、元功率损耗。 0021 步骤 1） 中， 建立研究对象在流场 - 温度场计算中的几何模型， 并根据流体力学理 论中流体运动特点， 将模型离散为四边形有限元网格。 0022 步骤 4） 中， 采用高斯积分和等参单元矩阵变换方法， 将功率损耗密度映射到流 场 - 温度场模型的各单元上。 0023 步骤 6） 中， 通过修正系数， 对流场 - 温度场模型的各网格映射后的功率损耗进行 修正， 并作为载荷代入能量方程， 施加边界条件， 对连续性方程、 Navier-Stokes 方程及能量 守恒方程进行数值求解后， 得到温度结果。 0024 本发明提供的电磁场 - 流场 - 温度场耦合计算中的异型网格映。
14、射方法， 允许在不 同物理场求解时， 模型以最适合控制方程求解特征的方式离散， 所有单元和节点不要求逐 一匹配， 而是通过高斯 - 勒让德积分、 等参单元矩阵变换及数据插值等方式来实现不同网 格间的数据传递 ； 在针对电磁装置的电磁场、 流场、 温度场耦合计算中， 提出异型网格映射 方法， 可解决不同物理场对网格剖分精度以及形状不同需求的问题， 可实现不同物理场采 用不同类型网格之间的数据传递， 大大减少了网格数量， 提高了计算时间及精度。 0025 除用于解决电磁 - 流体 - 温度的耦合场计算， 还可以推广至含有矢量力的数据传 递的电磁 - 结构的耦合场数值计算。 说 明 书 CN 103。
15、617367 A 4 3/6 页 5 附图说明 0026 下面结合附图和实施例对本发明作进一步说明 ： 0027 图 1 为本发明方法的流程图 ； 0028 图 2 是步骤 2） 采用的四边形单元的示意图 ； 0029 图 3 是步骤 3） 得到的三角形单元的示意图。 具体实施方式 0030 如图 1 所示， 本发明方法由以下步骤组成 ： 0031 1） 建立研究对象在流场 - 温度场计算中的几何模型 （即场域模型） ， 并根据流体力 学理论中流体运动特点， 将模型离散为四边形有限元网格 ； 0032 2）按逆时针顺序读取模型中各离散网格的各节点坐标值 (x,y)， 保存于文件 A （A.tx。
16、t） 中 ； 0033 3） 再次建立研究对象在电磁场计算中的几何模型 （即场域模型） ， 根据电磁场理论 中磁场分布特点， 将模型离散为三角形有限元网格 ； 0034 4） 对以下电磁场方程 （1） ，（2） 进行数值计算 ： 0035 0036 0037 V1是产生涡流的区域， V2是施加电流的区域。 0038 5） 读取文件 A 中各节点坐标值， 采用高斯积分和等参单元矩阵变换方法， 将功率 损耗密度映射到流场 - 温度场模型的各单元上 ； 0039 四边形单元如图 2 所示， 、 可以看作是四边形单元的自然坐标， 自然坐标系的 原点一般取在单元重心上。根据形函数的特征， 可以写出四边形。
17、单元的形函数为 ： 0040 0041 根据等参元将直角坐标转换为自然坐标， 结合高斯 - 勒让德积分公式， 单元功率 损耗可以表示为 ： 说 明 书 CN 103617367 A 5 4/6 页 6 0042 0043 高斯积分点的位置及其权系数如表 1 所示 ： 0044 表 1 高斯积分点与权重 0045 0046 0047 公式 (4) 和 (5) 中， i 和 j 为高斯积分点， wi和 wj为各高斯积分点的权系数， 为电 导率， J(xi,yj) 为各高斯点在电磁场计算模型中对应的电流密度， Ja为雅克比矩阵， 即 0048 0049 其中 s 为节点编号， 按逆时针顺序读取。 说。
18、 明 书 CN 103617367 A 6 5/6 页 7 0050 0051 将不同的高斯积分点数值和每个单元四个节点坐标代入式 (6) 和 (7)， 即可以求 出雅克比矩阵。由于 J=E， 而场向量如果使用一次单元进行数值计算， 则 无法通过插值得到各高斯积分点处不同的电流密度值， 因此在涡流场计算中， 使用二次单 元六节点三角形单元进行离散， 三角形单元如图 3 所示 ； 0052 用面积坐标 （L1， L2， L3） 表示此二次单元的插值函数。其中三个角点的插值函数可 表示为 ： 0053 Ni (2Li-1)Li (i 1， 2， 3) (8) 0054 三个边中点表示为 ： 005。
19、5 0056 求解涡流场方程之后， 将流场计算有限元模型中四边形单元所有高斯点的坐标转 换到电磁场计算有限元模型中的坐标位置， 通过插值可得到各点的电流密度及功率损耗， 代入式 (4) 可得到温度场中各单元功率损耗。 0057 6） 将映射后的总功率与电磁场计算结果中导体总功率进行对比， 得到映射修正系 数 ； 0058 0059 7） 通过修正系数， 对流场 - 温度场模型的各网格映射后的功率损耗进行修正， 并 作为载荷代入能量方程， 施加边界条件， 对连续性方程 （11） 、 Navier-Stokes 方程 （12） 及能 量守恒方程 （13） 进行数值求解。 0060 0061 说 明 书 CN 103617367 A 7 6/6 页 8 0062 0063 得到温度结果 ； 0064 从而得到温度场中各单元功率损耗。 说 明 书 CN 103617367 A 8 1/1 页 9 图 1 图 2 图 3 说 明 书 附 图 CN 103617367 A 9 。