[发明专利]一种考虑气隙层影响的电磁吸附力计算方法有效
申请号: | 201710565851.4 | 申请日: | 2017-07-12 |
公开(公告)号: | CN107315890B | 公开(公告)日: | 2020-07-31 |
发明(设计)人: | 江鹏;张群 | 申请(专利权)人: | 英特工程仿真技术(大连)有限公司 |
主分类号: | G06F30/23 | 分类号: | G06F30/23 |
代理公司: | 大连东方专利代理有限责任公司 21212 | 代理人: | 王丹;李洪福 |
地址: | 116000 辽宁省大连市高*** | 国省代码: | 辽宁;21 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 考虑 气隙层 影响 电磁 吸附力 计算方法 | ||
1.一种考虑气隙层影响的电磁吸附力计算方法,其特征在于,包括如下步骤:
步骤1、建立邻近物体间含有气隙层的电磁场模型;
步骤2、建立所述气隙层两侧的邻近物体间所满足的物理关系,并将所建立的物理关系作为约束条件,引入到步骤1的电磁场模型所对应的控制方程中,得到含约束条件的控制方程对应的等效积分形式;
所述步骤2包括:
首先根据电磁场边界条件,建立区域Ωm、Ωs和Ωg间的电磁场量的关系式即气隙边界条件:
其中分别表示Γm面上沿Ωm外法向和切向的通量密度;分别表示Γs面上沿Ωs外法向和切向的通量密度;qg表示气隙层Ωg中沿Ωm外法向的通量密度;
其次采用电标量势或磁标量势来描述气隙边界条件(5a),则有:
在边界Γm处,表示为:
在边界Γs处,表示为:
在式(6)和(7)中,nm表示Γm面上的Ωm外法向,ns表示Γs面上的Ωs外法向,且它们的方向正好相反;αm和αs分别表示主域Ωm和从域Ωs的材料参数;和分别表示主域Ωm和从域Ωs的标量势;α0表示空气材料参数;d0表示气隙层的厚度;h0表示气隙层的等效电阻或磁阻;
最后分别在Ωm和Ωs中引入试探函数wm和ws,并应用Green公式,得到(1)、(6)、(7)的等效积分形式:
步骤3、采用有限元方法离散步骤2得到的所述等效积分形式,并得到含约束条件的控制方程所对应的有限元方程;
所述步骤3采用有限元方法进行离散,其包括:
步骤31、对区域Ωm和Ωs划分网格,并记Ni(i=1,2…n)为标准有限元空间的节点基函数,φ=[φ1,φ2…φn]T为标量势的节点自由度,并将求解域中标量势表示为节点基函数和节点自由度的乘积形式:
并将Ωm上的自由度拆分为两部分,即分别对应Ωm/Γm和Γm上的自由度分别为φim和将Ωs上的自由度拆分为两部分,即分别对应Ωs/Γs和Γs上的自由度分别φis和于是将φ表示为:
步骤32、将(10)和(11)式以及wm,s=Ni(i=1,2…n)代入到(8)和(9)中,得到含气隙边界条件的有限元方程即所述含约束条件的控制方程所对应的有限元方程:
将(12)式中统一表示成Kmm,则具有如下分量表达形式:
将(12)式中统一表示成Kss,则具有如下分量表达形式:
(13)和(14)式中,n2表示中的最大下标,
将(12)式等号右端中的载荷矢量fim、统一表示成fm,则具有如下分量表达形式:
将(12)式等号右端中的载荷矢量fis、统一表示成fs,则具有如下分量表达形式:
由于(12)式中是从(8)和(9)中Γm和Γs面上的边界积分得到,将其统一表示成NΓΓ,则具有如下分量表达形式:
(17)式中,n1是的最小下标;n3是的最大上标;
步骤4、求解步骤3所得到的有限元方程,根据所求得的有限元方程的数值解以及所建立的物理关系,得到与所述气隙层对应的电磁场分布;
步骤5、根据所述气隙层对应的电磁场分布,采用电磁力计算方法求得邻近物体间的电磁吸附力。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于英特工程仿真技术(大连)有限公司,未经英特工程仿真技术(大连)有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201710565851.4/1.html,转载请声明来源钻瓜专利网。