CN109408927B - 一种基于黑盒传输线模型的二维静磁场并行有限元加速方法 - Google Patents
一种基于黑盒传输线模型的二维静磁场并行有限元加速方法 Download PDFInfo
- Publication number
- CN109408927B CN109408927B CN201811192762.0A CN201811192762A CN109408927B CN 109408927 B CN109408927 B CN 109408927B CN 201811192762 A CN201811192762 A CN 201811192762A CN 109408927 B CN109408927 B CN 109408927B
- Authority
- CN
- China
- Prior art keywords
- solving
- iteration
- finite element
- nonlinear
- transmission line
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于黑盒传输线模型的二维非线性静磁场模型的有限元求解方法,所述方法包括如下步骤:一、确定待求解的变量以及求解域;二、建立rz直角坐标系;三、列出二维轴对称非线性静磁场中的控制方程和边界条件式并组成一微分方程组;四、对求解域进行分网;五、计算有限元的系数矩阵;六、建立等效的黑盒电路模型;七、在黑盒电路模型与线性电路之间***一段传输线线段;八、进行传输线法迭代;九、重复进行迭代,直到迭代结果收敛到固定误差,结束求解。本发明能够对所求问题并行地进行求解计算,并且能够用于复杂的有限元分网模型求解当中,从而解决牛顿迭代法求解有限元非线性问题时带来的求解时间长、效率低的问题。
Description
技术领域
本发明属于电磁场数值计算技术领域,涉及一种二维非线性静磁场模型的有限元求解方法,具体涉及一种将二维非线性静磁场的有限元问题转化为等价的黑盒电路问题,并借助传输线迭代法实现对二维非线性静磁场的有限元并行加速求解的方法。
背景技术
有限元法是工业领域中应用最广泛的数值计算方法,被诸多商用仿真软件采用。然而,随着求解模型的日益复杂化以及分网单元数目的不断增多,以传统的牛顿迭代法为核心的非线性有限元求解方法面临着求解耗时严重的问题,这直接关系到产品的仿真研发的速度和效率。
有限元问题的求解的核心在于求解线性方程组,而对于非线性问题来说,传统的牛顿迭代法每一步都要利用新的迭代结果重新生成有限元模型的全局矩阵,随着模型分网的不断增大,全局矩阵的维度不断变大,每一步矩阵的LU分解等消耗的时间会相应的增大,总体的求解时间可能随着分网的变密而成几何式增大。使用传输线迭代法能够避免重复的LU分解过程,求解时间。然而,传统的传输线迭代法只能用于一阶三角形单元的有限元模型当中,无法用于复杂的分网单元模型。因此,需要研究一种新的迭代方法,来提高传输线法的适用范围和应对复杂模型的能力,从而解决牛顿迭代法求解有限元非线性问题时带来的求解时间长、效率低的问题。
发明内容
本发明的目的是提供一种基于黑盒传输线模型的二维非线性静磁场模型的有限元求解方法,该方法能够对所求问题并行地进行求解计算,并且能够用于复杂的有限元分网模型求解当中,包括三角形单元和四边形单元等,从而解决牛顿迭代法求解有限元非线性问题时带来的求解时间长、效率低的问题。
本发明的目的是通过以下技术方案实现的:
一种基于黑盒传输线模型的二维非线性静磁场模型的有限元求解方法,包括如下步骤:
一、确定待求解的变量以及求解域,其中:所述待求解的变量为一二维轴对称非线性静磁场的磁势A,所述求解域为二维轴对称非线性静磁场所在的区域;
二、建立一个rz直角坐标系;
三、列出二维轴对称非线性静磁场中的控制方程和边界条件式并组成一微分方程组,其中:
所述控制方程为:
式中,J为电流密度变量,μ′=rμ,μ为三角单元的磁导率,A′=rA,A为磁势,r为横坐标,z为纵坐标;
所述边界条件式为:
Γ1:A′=0;
式中,m为求解边界上的法向量;
五、按照下式计算有限元的系数矩阵:
其中,每一项的计算表达式如下:
六、将步骤五当中的矩阵Ke看作电路的导纳矩阵,be看作电流源向量,建立等效的黑盒电路模型;
七、在黑盒电路模型与线性电路之间的连接导线任意位置***一段传输线线段;
八、进行传输线法迭代:
将电压信号入射到线性网络中,求解该电路,即方程组:(Ylinear+YTL)A′=b+2ViYTL,在整个迭代过程当中,矩阵(Ylinear+YTL)保持不变,在迭代的第1步执行一次矩阵的LU分解操作,后续步骤无需再次执行,求解完成后,计算反射电压Vr=A′-Vi;
式中,Ylinear为入射到线性网络内的线性电路的导纳矩阵;YTL为入射到线性网络内的传输线电路的导纳矩阵;A′为待求节点电压;Vr为反射到非线性电路的反射电压;Vi为入射到线性电路的入射电压;Ke为非线性单元的系数矩阵;为反射回非线性单元内的传输线的导纳矩阵;为反射回非线性单元内的反射电压;Je为雅克比矩阵;为第k次迭代时的雅克比矩阵;f为计算当中的中间变量;Ae′为单元内每个节点的磁势值。
九、重复进行步骤八的迭代,直到迭代结果收敛到固定误差,结束求解。
相比于现有技术,本发明具有如下优点:
1、采用有限元方法对静磁场进行离散建模,将有限元中的每个单元的单元系数矩阵等效为一个密封的黑盒电路模型,从而实现将有限元方程组等效为一个非线性的电路网络问题。
2、在非线性单元的黑盒模型与线性单元的网络之间***传输线线段,利用传输线的迭代规律,实现传输线迭代,从而实现静磁场问题的求解。
3、在传输线迭代过程中,每个非线性黑盒模型都是相互独立的,它的求解计算可以采用单独的计算资源(CPU或者GPU核心)进行计算,实现并行加速。
4、采用传输线法和黑盒模型能够进行有限元的迭代求解,能够处理多种复杂的有限元模型,例如三角形单元和四边形单元。
5、在迭代求解过程中,全局矩阵Y能够保持不变,在矩阵求解过程当中,采用LU分解法,只需要在计算的第一步进行LU分解,由于LU分解一般占用矩阵求解的95%左右的时间,使用这种方法将大幅减小每一步的迭代时间。
6、使用传输线法能够将非线性单元与线性求解区域进行隔离,实现并行加速效果。
附图说明
图1为黑盒电路模型;
图2为***传输线后的黑盒电路网络;
图3为电压入射到线性网络后的等效电路图;
图4为电压反射回非线性单元后的等效电路图;
图5为迭代求解流程图;
图6为某型号接触器的结构图;
图7为某型号接触器的剖面图;
图8为接触器简化模型;
图9为接触器计算模型;
图10为分网结果,(a)三角形单元,(b)四边形单元;
图11为采用传输线法求解得到不同状态下的磁场分布情况;
图12为单步计算时间加速效果。
具体实施方式
下面结合附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的范围,均应涵盖在本发明的保护范围中。
具体实施方式一:本实施方式提供了一种基于黑盒传输线模型的二维非线性静磁场模型的有限元求解方法,所述方法包括如下步骤:
一、确定待求解的变量以及求解域,待求解的变量为一二维轴对称非线性静磁场的磁势A,二维轴对称非线性静磁场由通电线圈中的电流产生,通电线圈周围的各元件均为铁磁材料,求解域为二维轴对称非线性静磁场所在的区域。
二、建立一个rz直角坐标系,r为横轴,z为纵轴。
三、列出二维轴对称非线性静磁场中的控制方程和边界条件式并组成一微分方程组,其控制方程为:其中,J为电流密度变量,μ′=rμ,μ为三角单元的磁导率,A′=rA,A为磁势,r为横坐标,z为纵坐标,边界条件为:Γ1:A′=0;m为求解边界上的法向量。
五、计算有限元的系数矩阵。根据伽辽金法,对控制方程进行加权积分,得到方程:
写成矩阵形式KeAe′=be,即为:
其中,每一项的计算表达式如下:
六、将步骤五当中的矩阵Ke看作电路的导纳矩阵,be看作电流源向量,建立等效的黑盒电路模型,如图1所示。
七、黑盒模型与线性电路之间是通过导线相连接的,在黑盒模型与线性电路之间的这些导线的任意位置***一段传输线线段,如图2所示。
八、进行传输线法迭代。
电压信号入射到线性网络中,如图3所示,求解该电路,即方程组:(Ylinear+YTL)A′=b+2ViYTL,在整个迭代过程当中,矩阵(Ylinear+YTL)保持不变,在迭代的第1步执行一次矩阵的LU分解操作,后续步骤无需再次执行,求解完成后,计算反射电压Vr=A′-Vi。
式中,Ylinear为图3当中线性电路的导纳矩阵;YTL为图3虚线框内传输线电路的导纳矩阵;A′为待求节点电压;Vr为反射到非线性电路的反射电压;Vi为入射到线性电路的入射电压;Ke为非线性单元的系数矩阵;为图4中传输线的导纳矩阵;为反射回单元内的反射电压;Je为雅克比矩阵;为第k次迭代时的雅克比矩阵。
不断地进行迭代,直到迭代结果收敛到固定误差,结束求解。求解流程图如图5。
具体实施方式二:本实施方式提供了一种基于黑盒传输线模型的二维非线性静磁场模型的有限元求解方法,所述方法具体实施步骤如下:
一、确定待求解的变量以及求解域,待求解的变量为一二维轴对称非线性静磁场的磁势A,二维轴对称非线性静磁场由通电线圈中的电流产生,通电线圈周围的各元件均为铁磁材料,求解域为二维轴对称非线性静磁场所在的区域。本实施方式针对某型号轴对称结构接触器进行计算,该模型的外观图如图6所示,剖面图如图7所示。对其进行简化,并建模,简化模型如图8所示,计算模型如图9所示。
二、建立一个rz直角坐标系,r为横轴,z为纵轴。
三、列出二维轴对称非线性静磁场中的控制方程和边界条件式并组成一微分方程组,其控制方程为:其中,J为电流密度变量,μ′=rμ,μ为三角单元的磁导率,A′=rA,A为磁势,r为横坐标,z为纵坐标,边界条件为:Γ1:A′=0;m为求解边界上的法向量。
四、采用分网程序对求解区域进行分网,将求解域离散为三角形单元或者四边形单元,分网结果如图10所示。在每一个有限元单元中,变量其中,n为有限元单元的节点数目,Nj为单元当中的形函数,为对应节点的变量值的大小。
五、计算有限元的系数矩阵。根据伽辽金法,对控制方程进行加权积分,得到方程:
写成矩阵形式KeAe′=be,即为:
其中,每一项的计算表达式如下:
六、将步骤五当中的矩阵Ke看作电路的导纳矩阵,be看作电流源向量,建立等效的黑盒电路模型,如图1所示。
七、在黑盒模型与线性电路之间***一段传输线线段,如图2所示。
八、进行传输线法迭代。
电压信号入射到线性网络中,如图3所示,求解该电路,即方程组:(Ylinear+YTL)A′=b+2ViYTL,然后,计算反射电压Vr=A′-Vi。
不断地进行迭代,直到迭代结果收敛到固定误差,结束求解。求解流程图如图5。
采用传输线法求解得到不同状态下的磁场分布情况,结果如图11所示。
表1不同的求解模型
case 1-5为三角形单元分网;case 6-10为四边形单元分网。
表2 N-R迭代法和本发明BB-TLM迭代法的误差对比
No.1-5为表1case 5模型当中任取五个点的误差对比;No.6-10为表1case 10模型任取五个点的误差对比。
对表1当中case10模型分别使用1、4、8、16、20个CPU核心进行加速,对比牛顿迭代法与传输线法的单步计算时间,由图12可见,传输线法的单步计算时间比牛顿迭代法快几倍。
对表1当中case10模型分别使用1、4、8、16、20个CPU核心进行加速,对比牛顿迭代法与传输线法的总计算时间,当CPU核心增大时,传输线法的加速效果比牛顿迭代法更加明显。
表3针对表1中的模型分别采用牛顿迭代法和传输线法求解(CPU核心为20)
由表3可知,当分网模型增大时,传输线法的计算时间比牛顿迭代法快。
Claims (1)
1.一种基于黑盒传输线模型的二维非线性静磁场模型的有限元求解方法,其特征在于所述方法包括如下步骤:
一、确定待求解的变量以及求解域,其中:所述待求解的变量为一二维轴对称非线性静磁场的磁势A,所述求解域为二维轴对称非线性静磁场所在的区域;
二、建立一个rz直角坐标系;
三、列出二维轴对称非线性静磁场中的控制方程和边界条件式并组成一微分方程组,其中:
所述控制方程为:
式中,J为电流密度变量,μ′=rμ,μ为三角单元的磁导率,A′=rA,A为磁势,r为横坐标,z为纵坐标;
所述边界条件式为:
Γ1:A′=0;
式中,m为求解边界上的法向量;
五、按照下式计算有限元的系数矩阵:
其中,每一项的计算表达式如下:
六、将步骤五当中的矩阵Ke看作电路的导纳矩阵,be看作电流源向量,建立等效的黑盒电路模型;
七、在黑盒电路模型与线性电路之间的连接导线任意位置***一段传输线线段;
八、进行传输线法迭代:
将电压信号入射到线性网络中,求解该电路,即方程组:(Ylinear+YTL)A′=b+2ViYTL,在整个迭代过程当中,矩阵(Ylinear+YTL)保持不变,在迭代的第1步执行一次矩阵的LU分解操作,后续步骤无需再次执行,求解完成后,计算反射电压Vr=A′-Vi;
将电压信号入射到非线性有限元单元中,求解该电路,即非线性方程组:采用牛顿迭代法进行求解,第k步的迭代公式如下:每个非线性单元内的方程组求解都单独放到独立的计算核心中求解,实现并行计算,求解完成后,再次计算入射到线性电路的入射电压
式中,Ylinear为入射到线性网络内的线性电路的导纳矩阵;YTL为入射到线性网络内的传输线电路的导纳矩阵;A′为待求节点电压;Vr为反射到非线性电路的反射电压;Vi为入射到线性电路的入射电压;Ke为非线性单元的系数矩阵;为反射回非线性单元内的传输线的导纳矩阵;为反射回非线性单元内的反射电压;为第k次迭代时的雅克比矩阵;
九、重复进行步骤八的迭代,直到迭代结果收敛到固定误差,结束求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811192762.0A CN109408927B (zh) | 2018-10-13 | 2018-10-13 | 一种基于黑盒传输线模型的二维静磁场并行有限元加速方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811192762.0A CN109408927B (zh) | 2018-10-13 | 2018-10-13 | 一种基于黑盒传输线模型的二维静磁场并行有限元加速方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109408927A CN109408927A (zh) | 2019-03-01 |
CN109408927B true CN109408927B (zh) | 2022-01-18 |
Family
ID=65467151
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811192762.0A Active CN109408927B (zh) | 2018-10-13 | 2018-10-13 | 一种基于黑盒传输线模型的二维静磁场并行有限元加速方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109408927B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112016232B (zh) * | 2020-08-31 | 2024-06-14 | 中国原子能科学研究院 | 一种撕裂有限元过程处理方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105631105A (zh) * | 2015-12-24 | 2016-06-01 | 北京航空航天大学 | 一种含外加场耦合的线缆等效电路建模方法 |
CN106649939A (zh) * | 2016-09-28 | 2017-05-10 | 哈尔滨工业大学 | 基于传输线迭代的2d轴对称非线性静磁场模型的求解方法 |
CN106991263A (zh) * | 2017-05-25 | 2017-07-28 | 清华大学 | 一种中低压配电变压器绕组雷电宽频模型的建模方法 |
CN107609274A (zh) * | 2017-09-14 | 2018-01-19 | 哈尔滨工业大学 | 基于传输线与级别调度法的二维静磁场并行有限元方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004109560A1 (ja) * | 2003-06-05 | 2004-12-16 | Fujitsu Limited | 回路解析と電磁波解析との連携による電磁界強度算出装置および算出方法 |
KR102066894B1 (ko) * | 2015-09-30 | 2020-01-16 | 오사카 유니버시티 | 전기회로에 있어서의 전위, 전류, 및 주변 전자계의 계산 방법, 및 프로그램 |
-
2018
- 2018-10-13 CN CN201811192762.0A patent/CN109408927B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105631105A (zh) * | 2015-12-24 | 2016-06-01 | 北京航空航天大学 | 一种含外加场耦合的线缆等效电路建模方法 |
CN106649939A (zh) * | 2016-09-28 | 2017-05-10 | 哈尔滨工业大学 | 基于传输线迭代的2d轴对称非线性静磁场模型的求解方法 |
CN106991263A (zh) * | 2017-05-25 | 2017-07-28 | 清华大学 | 一种中低压配电变压器绕组雷电宽频模型的建模方法 |
CN107609274A (zh) * | 2017-09-14 | 2018-01-19 | 哈尔滨工业大学 | 基于传输线与级别调度法的二维静磁场并行有限元方法 |
Non-Patent Citations (4)
Title |
---|
"Influence of electromagnetic pulse on electromagnetic relay by using equivalent circuit model";Li Bing等;《26th International Conference on Electrical Contacts (ICEC 2012)》;20120913;全文 * |
"传输线理论及电磁兼容计算的半解析方法研究";张旭锋;《中国博士学位论文全文数据库 信息科技辑》;20120715;全文 * |
"各向异性地层中可控源电磁法一维全参数反演及三维有限体积正演算法研究";周建美;《中国博士学位论文全文数据库 基础科学辑》;20140915;全文 * |
"基于巨磁阻抗效应的弱磁传感器技术研究";张雷;《中国博士学位论文全文数据库 信息科技辑》;20160115;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109408927A (zh) | 2019-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110110413B (zh) | 一种基于材料场缩减级数展开的结构拓扑优化方法 | |
JP7133894B2 (ja) | データに基づくインタラクティブ3dエクスペリエンス | |
RU2519331C2 (ru) | Способ вычисления физического значения, способ численного анализа, программа вычисления физического значения, программа численного анализа, устройство вычисления физического значения и устройство численного анализа | |
CN107526860B (zh) | 基于电场能建模技术的vlsi标准单元布局方法 | |
Hel-Or et al. | Relaxed parametric design with probabilistic constraints | |
CN110941021B (zh) | 基于网格点格架函数的重力异常及梯度异常的正演方法 | |
CN104361185A (zh) | 电缆虚拟设计用布线空间自动生成方法 | |
CN109408927B (zh) | 一种基于黑盒传输线模型的二维静磁场并行有限元加速方法 | |
WO2024007652A1 (zh) | 一种大型稀疏矩阵加速求解方法、***及存储介质 | |
CN108763777B (zh) | 基于泊松方程显式解的vlsi全局布局模型建立方法 | |
Brandi et al. | Improving stochastic dynamic programming on hydrothermal systems through an iterative process | |
CN111339688B (zh) | 基于大数据并行算法求解火箭仿真模型时域方程的方法 | |
CN116861822B (zh) | 基于笛卡尔网格的物面边界的处理方法及装置 | |
CN109960841B (zh) | 一种流体表面张力的仿真方法、终端设备及存储介质 | |
MOON et al. | Conservative treatment of boundary interfaces for overlaid grids andmulti-level grid adaptations | |
CN104216011B (zh) | 一种TTI介质稳定的qP波逆时偏移方法 | |
Zhou et al. | Reduced-order ℒ2–ℒ∞ filtering for singular systems: a linear matrix inequality approach | |
Burrows et al. | Learning overlap optimization for domain decomposition methods | |
Lin et al. | Turing instability and pattern formation of a fractional Hopfield reaction–diffusion neural network with transmission delay | |
Zhang et al. | Adaptive finite element analysis of elliptic problems based on bubble-type local mesh generation | |
Mavriplis | Viscous flow analysis using a parallel unstructured multigrid solver | |
CN113705969B (zh) | 一种圆柱坐标系下电流电荷分配方法及*** | |
US8676551B1 (en) | Multi-solver simulation of dynamic systems in a modeling environment | |
CN104809297A (zh) | 一种磁场-结构场耦合计算中的异型网格间的电磁力密度传递方法 | |
Sonar | Classical finite volume methods |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |