CN112347687A - 一种自适应自由度电磁-温度多物理场耦合分析方法 - Google Patents
一种自适应自由度电磁-温度多物理场耦合分析方法 Download PDFInfo
- Publication number
- CN112347687A CN112347687A CN202011383160.0A CN202011383160A CN112347687A CN 112347687 A CN112347687 A CN 112347687A CN 202011383160 A CN202011383160 A CN 202011383160A CN 112347687 A CN112347687 A CN 112347687A
- Authority
- CN
- China
- Prior art keywords
- electromagnetic
- temperature
- freedom
- field
- physical
- 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.)
- Granted
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 64
- 230000008878 coupling Effects 0.000 title claims abstract description 55
- 238000010168 coupling process Methods 0.000 title claims abstract description 55
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 55
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 239000006185 dispersion Substances 0.000 claims abstract description 5
- 238000009795 derivation Methods 0.000 claims abstract description 4
- 230000005672 electromagnetic field Effects 0.000 claims description 29
- 230000003044 adaptive effect Effects 0.000 claims description 22
- 239000011159 matrix material Substances 0.000 claims description 20
- 238000000034 method Methods 0.000 claims description 15
- 230000004907 flux Effects 0.000 claims description 9
- 230000014509 gene expression Effects 0.000 claims description 8
- 239000000203 mixture Substances 0.000 claims description 6
- 239000002994 raw material Substances 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000012546 transfer Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 abstract description 10
- 238000010586 diagram Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 3
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000002940 Newton-Raphson method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Soft Magnetic Materials (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种自适应自由度电磁‑温度多物理场耦合分析方法,包括:建立电工装备几何模型及初始非结构网格离散;建立电磁‑温度多物理场弱耦合数值模型;基于非结构网格单元的电磁‑温度多物理场弱耦合数值模型分析有限元计算格式推导;计算求解电工装备内的电磁‑温度多物理场,并对数值解进行误差分析;调整各物理场的自由度,并根据调整结果重新求解直到数值解精度满足要求,完成分析。本发明在一套网格上实现对于两个物理场自由度的独立灵活调整,以较小的计算资源满足物理场对于离散的不同要求,避免进行实际操作的网格稀疏及后续的复杂操作,以及使用网格映射函数及其可能引入的误差,有效提高电磁‑温度耦合分析计算的计算效率。
Description
技术领域
本发明涉及电工装备及仿真计算的技术领域,尤其涉及一种自适应自由度电磁-温度多物理场耦合分析方法。
背景技术
随着电力、能源等产业的发展,高能量密度及高效的电工装备成为众多应用场合的追求目标,在相应设计中需综合考虑电磁场、温度场等多个物理场以实现对于材料的极限利用;作为优化设计中的重要一环,多物理场仿真分析方法是实现该目标的关键支撑。基于物理场独立数值求解构建的电磁-温度多物理场弱耦合分析方法是目前主要的研究和发展方向,电磁场和温度场分别采用有限差分、有限体积、有限元等方法进行求解,并通过网格映射程序实现物理场间耦合量的传递,然而传统弱耦合分析方法存在映射函数选择以及映射误差较大等问题。
发明内容
本部分的目的在于概述本发明的实施例的一些方面以及简要介绍一些较佳实施例。在本部分以及本申请的说明书摘要和发明名称中可能会做些简化或省略以避免使本部分、说明书摘要和发明名称的目的模糊,而这种简化或省略不能用于限制本发明的范围。
鉴于上述现有存在的问题,提出了本发明。
因此,本发明解决的技术问题是:传统弱耦合分析方法存在映射函数选择以及映射误差较大等问题。
为解决上述技术问题,本发明提供如下技术方案:一种自适应自由度电磁-温度多物理场耦合分析方法,其特性在于,包括:建立电工装备几何模型并进行离散,得到初始非结构网格;根据物理场控制方程及分析问题的边界条件,建立电磁-温度多物理场弱耦合数值模型;基于所述非结构网格单元的电磁-温度多物理场弱耦合数值模型分析有限元计算格式推导;计算求解所述电工装备内的电磁-温度多物理场,并对数值解进行误差分析;基于误差分析结果,调整各物理场的自由度,并根据调整结果重新求解直到数值解精度满足要求,完成分析。
作为本发明所述的自适应自由度电磁-温度多物理场耦合分析方法的一种优选方案,其中:完全离散后的电磁场控制方程包括,
作为本发明所述的自适应自由度电磁-温度多物理场耦合分析方法的一种优选方案,其中:完全离散后的温度场控制方程包括,
其中,T为温度,Sij为温度场刚度矩阵系数,Fi为温度场载荷,λ为热导率,N为权函数/插值函数,ρ为密度,c为比热容,Δt为时间步长,as为传热系数,Ω为定义域,QV为热源体密度,qg为热流,T0为环境温度,Γ2为第二类边界,Γ3为第三类边界。
作为本发明所述的自适应自由度电磁-温度多物理场耦合分析方法的一种优选方案,其中:待求解的电磁-温度多物理场方程的离散表达式包括,
作为本发明所述的自适应自由度电磁-温度多物理场耦合分析方法的一种优选方案,其中:待求解的电磁场和温度场的耦合关系的离散表达式包括,
作为本发明所述的自适应自由度电磁-温度多物理场耦合分析方法的一种优选方案,其中:所述误差分析包括基于磁场能量构建误差分析器对电磁场数值近似解进行评估,具体包括,
作为本发明所述的自适应自由度电磁-温度多物理场耦合分析方法的一种优选方案,其中:对温度场数值近似解的误差进行评估包括,
作为本发明所述的自适应自由度电磁-温度多物理场耦合分析方法的一种优选方案,其中:基于误差计算结果,对误差大区域进行加密,在两个物理场完成操作后进行网格同步;并将足够精确区域的自由度从有限元空间中移除,完成求解后恢复,其中待移除的自由度由周围主自由度通过插值函数表示,即从自由度的限制条件,具体包括,
其中,下标和上标分别表示自由度的索引和属性(m和s分别表示主、从),n表示主节点的数目,αi表示主节点i的系数。
作为本发明所述的自适应自由度电磁-温度多物理场耦合分析方法的一种优选方案,其中:所述主节点i的系数的正规化的权重函数包括,
其中,权重wi(x)包括,
系数γa,b包括,
其中,x表示从节点的坐标向量,k和j分别表示主节点附近主节点的数目和索引,V表示由从节点和其他三个下标标注节点组成四面体的体积,xa/xb表示节点a/b的坐标向量。
作为本发明所述的自适应自由度电磁-温度多物理场耦合分析方法的一种优选方案,其中:所述调整各物理场的自由度包括,
基于所构建的限制条件,首先在单元层面实现自由度的移除,所述单元特征方程包括,
Ku=f
其中,基于所述限制条件,未知数向量表示为:
所述单元特征方程进行如下变形后,即从所述单元层面实现有限元空间内自由度的移除:
本发明的有益效果:本发明在一套网格上实现对于两个物理场自由度的独立灵活调整,以较小的计算资源满足物理场对于离散的不同要求,避免进行实际的网格稀疏及后续的复杂操作以及使用网格映射函数及其可能引入的误差,有效提高电磁-温度耦合分析计算的计算效率。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。其中:
图1为本发明一个实施例提供的一种自适应自由度电磁-温度多物理场耦合分析方法的基本流程示意图;
图2为本发明一个实施例提供的一种自适应自由度电磁-温度多物理场耦合分析方法的四面体单元循环加密示意图;
图3为本发明一个实施例提供的一种自适应自由度电磁-温度多物理场耦合分析方法的TEAM Workshop问题7的几何结构;
图4为本发明一个实施例提供的一种自适应自由度电磁-温度多物理场耦合分析方法的线段A1-B1上的Z方向磁通密度分量图;
图5为本发明一个实施例提供的一种自适应自由度电磁-温度多物理场耦合分析方法的铝板在时刻25ms的温度分布图;
图6为本发明一个实施例提供的一种自适应自由度电磁-温度多物理场耦合分析方法的铝板在时刻25ms的涡流分布图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合说明书附图对本发明的具体实施方式做详细的说明,显然所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明的保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
其次,此处所称的“一个实施例”或“实施例”是指可包含于本发明至少一个实现方式中的特定特征、结构或特性。在本说明书中不同地方出现的“在一个实施例中”并非均指同一个实施例,也不是单独的或选择性的与其他实施例互相排斥的实施例。
本发明结合示意图进行详细描述,在详述本发明实施例时,为便于说明,表示器件结构的剖面图会不依一般比例作局部放大,而且所述示意图只是示例,其在此不应限制本发明保护的范围。此外,在实际制作中应包含长度、宽度及深度的三维空间尺寸。
同时在本发明的描述中,需要说明的是,术语中的“上、下、内和外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一、第二或第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
本发明中除非另有明确的规定和限定,术语“安装、相连、连接”应做广义理解,例如:可以是固定连接、可拆卸连接或一体式连接;同样可以是机械连接、电连接或直接连接,也可以通过中间媒介间接相连,也可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
实施例1
本发明在考虑电磁场和温度场离散要求差异的基础上,基于非结构网格,提出一种用于电磁-温度多物理场弱耦合分析的自适应自由度有限元方法,在一套网格上实现对于两个物理场自由度的独立灵活调整,以较小的计算资源满足物理场对于离散的不同要求,同时避免使用网格映射函数及其可能引入的误差;该发明中,自由度的增加通过网格加密实现,并在两个物理场网格独立加密后进行同步;自由度的移除仅通过主从技术在有限元空间实现,并在求解完成后恢复,进而避免进行实际的网格稀疏及后续的复杂操作。
参照图1~2,为本发明的一个实施例,提供了一种自适应自由度电磁-温度多物理场耦合分析方法,包括:
S1:建立电工装备几何模型并进行离散,得到初始非结构网格;
具体的,基于所分析的电工装备,建立几何模型并确定求解域边界;
利用开源网格剖分软件Gmsh对几何模型进行离散,生成初始网格。
S2:根据物理场控制方程及分析问题的边界条件,建立电磁-温度多物理场弱耦合数值模型;
具体的,根据电工装备的运行情况,确定各物理场的模型类型及耦合方式;
基于物理场控制方程,结合分析问题的边界条件,建立电磁-温度多物理场弱耦合数值模型,其中电磁场控制方程为:
温度场的控制方程为:
电磁场与温度场之间的双向耦合关系为:电磁场分析得出的损耗作为温度场分析中的热源;温度场分析得出的温度用于更新电磁场分析中的材料属性;可以通过以下方程表征:
Pcu=VJ2/σt
σt=σ0/(1.0+ασ×(T-T0))
μt=μ0(1.0+αμ×(T-T0))
其中,θ表示磁阻率,σ表示电导率,Js表示源电流密度,A表示矢量磁位,表示标量电位,T表示温度,ρ表示密度,c表示比热容,λ表示热导率,QV表示热源体密度,μ表示磁导率,α表示温度系数,T0表示参考温度,kh表示磁滞损耗系数,kc表示经典涡流系数,ke表示附加损耗系数,J表示电流密度,V表示单元体积,f表示频率,Bm表示磁通密度。
S3:基于非结构网格单元的电磁-温度多物理场弱耦合数值模型分析有限元计算格式推导;
具体的,首先采用伽辽金法对场控制方程进行处理:
定义矢量磁位、标量电位及温度的试探解为:
进一步的,在四面体单元特征矩阵建立过程中,单元基函数除用于表示场势函数的试探解外,还用作伽辽金法中的权函数,其可以写为:
基于由四面体单元确定的权函数及试探解,对控制方程进行离散,进而得出单元特征矩阵的表达式。
其中,完全离散后的电磁场控制方程包括,
完全离散后的温度场控制方程包括,
其中,T为温度,Sij为温度场刚度矩阵系数,Fi为温度场载荷,λ为热导率,N为权函数/插值函数,ρ为密度,c为比热容,Δt为时间步长,as为传热系数,Ω为定义域,QV为热源体密度,qg为热流,T0为环境温度,Γ2为第二类边界,Γ3为第三类边界。
更进一步的,基于单元局部索引与整体索引的对应关系,由单元特征矩阵组装为整体特征矩阵,推出待求解的离散方程表达式,其中电磁-温度多物理场方程的离散表达式为:
电磁场与温度场之间耦合关系的离散表达式为:
S4:计算求解电工装备内的电磁-温度多物理场,并对数值解进行误差分析;
利用牛顿-拉夫逊方法处理由电工材料引入的非线性,在把离散方程线性化后采用直接求解器进行求解;在求解完成后,采用误差分析器分别对两个物理场的近似解进行误差分析;由于多数位置的真实解无法得到,通过对近似解进行处理得到其估计值:
进而基于磁场能量构建误差分析器对电磁场近似解进行评估:
同理,对温度场近似解的误差进行评估:
S5:基于误差分析结果,调整各物理场的自由度,并根据调整结果重新求解直到数值解精度满足要求,完成分析;
基于误差计算结果,对误差大区域进行加密,在两个物理场完成操作后进行网格同步;并将足够精确区域的自由度从有限元空间中移除,完成求解后恢复,其中待移除的自由度由周围主自由度通过插值函数表示,即从自由度的限制条件,具体包括,
其中,下标和上标分别表示自由度的索引和属性(m和s分别表示主、从),n表示主节点的数目,αi表示主节点i的系数。
主节点i的系数的正规化的权重函数包括,
其中,权重wi(x)包括,
系数γa,b包括,
其中,x表示从节点的坐标向量,k和j分别表示主节点附近主节点的数目和索引,V表示由从节点和其他三个下标标注节点组成四面体的体积,xa/xb表示节点a/b的坐标向量。
在完成限制条件构建后,首先在单元层面实现自由度的移除,以温度场中一个节点2为从节点的网格为例,其单元特征方程为:
Ku=f
基于限制条件,未知数向量可以重写为:
原方程在进行如下变形后,即从单元层面实现有限元空间内自由度的移除,后续的整体特征矩阵组装与步骤S3一致:
重复步骤S4和步骤S5,直至数值解精度满足要求。
实施例2
参照图3~6,为本发明的另一个实施例,为对本方法中采用的技术效果加以验证说明,本实施例采用传统弱耦合分析方法与本发明方法进行对比测试,以科学论证的手段对比计算结果,以验证本方法所具有的真实效果。
该实施例为暂态电磁-温度耦合问题,自适应自由度有限元步进计算以0.5ms的时间步长计算到时间节点25ms,在自适应计算过程中生成了一个具有467202单元的网格,并且在每步计算中大约有17%的磁场自由度和29%的热场自由度被移出有限元空间,如图4~6所示,在时间节点25ms求得的磁通密度Z方向分量与测量结果吻合良好,且求得的温度分布与传统方案结果保持一致。对于实施例,所提方法在不影响计算精度的前提下,将计算量减小了18%左右,验证了本发明方法能够有效提高电磁-温度耦合分析的计算效率。
应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (10)
1.一种自适应自由度电磁-温度多物理场耦合分析方法,其特性在于,包括:
建立电工装备几何模型并进行离散,得到初始非结构网格;
根据物理场控制方程及分析问题的边界条件,建立电磁-温度多物理场弱耦合数值模型;
基于所述非结构网格单元的电磁-温度多物理场弱耦合数值模型分析有限元计算格式推导;
计算求解所述电工装备内的电磁-温度多物理场,并对数值解进行误差分析;
基于误差分析结果,调整各物理场的自由度,并根据调整结果重新求解直到数值解精度满足要求,完成分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011383160.0A CN112347687B (zh) | 2020-12-01 | 2020-12-01 | 一种自适应自由度电磁-温度多物理场耦合分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011383160.0A CN112347687B (zh) | 2020-12-01 | 2020-12-01 | 一种自适应自由度电磁-温度多物理场耦合分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112347687A true CN112347687A (zh) | 2021-02-09 |
CN112347687B CN112347687B (zh) | 2021-11-12 |
Family
ID=74427343
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011383160.0A Active CN112347687B (zh) | 2020-12-01 | 2020-12-01 | 一种自适应自由度电磁-温度多物理场耦合分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112347687B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113204887A (zh) * | 2021-05-17 | 2021-08-03 | 上海交通大学 | 一种电力装备的多物理场多参数反演方法 |
CN113591346A (zh) * | 2021-07-09 | 2021-11-02 | 上海大学 | 一种用于电磁场求解的自适应区域分解有限元方法 |
CN113722860A (zh) * | 2021-09-07 | 2021-11-30 | 上海交通大学 | 基于降阶模型的瞬态热力状态在线评估方法、设备及介质 |
CN113962086A (zh) * | 2021-10-22 | 2022-01-21 | 中国电子科技集团公司第二十六研究所 | 一种多物理场耦合的声表面波滤波器的计算方法 |
CN113962087A (zh) * | 2021-10-22 | 2022-01-21 | 中国电子科技集团公司第二十六研究所 | 一种耦合温度场的声表面波谐振器的计算方法 |
CN113962089A (zh) * | 2021-10-22 | 2022-01-21 | 中国电子科技集团公司第二十六研究所 | 一种基于精确理论解的声表面波滤波器的设计方法 |
CN114462264A (zh) * | 2021-12-14 | 2022-05-10 | 深圳供电局有限公司 | 一种基于多物理场耦合的电缆接头状态评估方法及装置 |
CN115933409A (zh) * | 2023-01-06 | 2023-04-07 | 青岛哈尔滨工程大学创新发展中心 | 一种基于多物理场弱耦合的自适应调整时间步长的方法 |
CN117236110A (zh) * | 2023-04-13 | 2023-12-15 | 安徽大学 | 基于电热多物理场随机fdtd的电磁辐射强度评估方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106295053A (zh) * | 2016-08-20 | 2017-01-04 | 三峡大学 | 一种基于自适应时间步长的瞬变电磁‑温度场耦合计算方法 |
KR20190103750A (ko) * | 2018-02-28 | 2019-09-05 | 조선대학교산학협력단 | 유한요소법을 이용한 원형단면 탄소섬유복합재(cfrp)의 정적압궤 특성 데이터 획득 방법 |
CN110427664A (zh) * | 2019-07-18 | 2019-11-08 | 大连理工大学 | 一种基于磁热耦合的永磁耦合器温度场分析方法 |
CN111881611A (zh) * | 2020-07-31 | 2020-11-03 | 珠海格力电器股份有限公司 | 电机物理场的仿真处理方法及装置 |
-
2020
- 2020-12-01 CN CN202011383160.0A patent/CN112347687B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106295053A (zh) * | 2016-08-20 | 2017-01-04 | 三峡大学 | 一种基于自适应时间步长的瞬变电磁‑温度场耦合计算方法 |
KR20190103750A (ko) * | 2018-02-28 | 2019-09-05 | 조선대학교산학협력단 | 유한요소법을 이용한 원형단면 탄소섬유복합재(cfrp)의 정적압궤 특성 데이터 획득 방법 |
CN110427664A (zh) * | 2019-07-18 | 2019-11-08 | 大连理工大学 | 一种基于磁热耦合的永磁耦合器温度场分析方法 |
CN111881611A (zh) * | 2020-07-31 | 2020-11-03 | 珠海格力电器股份有限公司 | 电机物理场的仿真处理方法及装置 |
Non-Patent Citations (2)
Title |
---|
YUNPENG ZHANG: "< An adaptive degrees-of-freedom finite element method for 3-D nonlinear magneto-thermal field analysis>", 《NUMERICAL HEAT TRANSFER,PART A》 * |
YUNPENG ZHANG: "<3-D Transient Magneto-Thermal Field Analysis Using Adaptive Degrees-of-Freedom Finite-Element Method>", 《IEEE TRANSACTIONS ON MAGNETICS》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113204887B (zh) * | 2021-05-17 | 2023-08-25 | 上海交通大学 | 一种电力装备的多物理场多参数反演方法 |
CN113204887A (zh) * | 2021-05-17 | 2021-08-03 | 上海交通大学 | 一种电力装备的多物理场多参数反演方法 |
CN113591346A (zh) * | 2021-07-09 | 2021-11-02 | 上海大学 | 一种用于电磁场求解的自适应区域分解有限元方法 |
CN113722860A (zh) * | 2021-09-07 | 2021-11-30 | 上海交通大学 | 基于降阶模型的瞬态热力状态在线评估方法、设备及介质 |
CN113722860B (zh) * | 2021-09-07 | 2024-02-23 | 上海交通大学 | 基于降阶模型的瞬态热力状态在线评估方法、设备及介质 |
CN113962086A (zh) * | 2021-10-22 | 2022-01-21 | 中国电子科技集团公司第二十六研究所 | 一种多物理场耦合的声表面波滤波器的计算方法 |
CN113962089A (zh) * | 2021-10-22 | 2022-01-21 | 中国电子科技集团公司第二十六研究所 | 一种基于精确理论解的声表面波滤波器的设计方法 |
CN113962087A (zh) * | 2021-10-22 | 2022-01-21 | 中国电子科技集团公司第二十六研究所 | 一种耦合温度场的声表面波谐振器的计算方法 |
CN114462264A (zh) * | 2021-12-14 | 2022-05-10 | 深圳供电局有限公司 | 一种基于多物理场耦合的电缆接头状态评估方法及装置 |
CN114462264B (zh) * | 2021-12-14 | 2024-10-08 | 深圳供电局有限公司 | 一种基于多物理场耦合的电缆接头状态评估方法及装置 |
CN115933409A (zh) * | 2023-01-06 | 2023-04-07 | 青岛哈尔滨工程大学创新发展中心 | 一种基于多物理场弱耦合的自适应调整时间步长的方法 |
CN115933409B (zh) * | 2023-01-06 | 2023-05-05 | 青岛哈尔滨工程大学创新发展中心 | 一种基于多物理场弱耦合的自适应调整时间步长的方法 |
CN117236110A (zh) * | 2023-04-13 | 2023-12-15 | 安徽大学 | 基于电热多物理场随机fdtd的电磁辐射强度评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112347687B (zh) | 2021-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112347687B (zh) | 一种自适应自由度电磁-温度多物理场耦合分析方法 | |
Chaniotis et al. | Model reduction in power systems using Krylov subspace methods | |
Dal et al. | Computational electro-chemo-mechanics of lithium-ion battery electrodes at finite strains | |
CN103413013B (zh) | 一种锂离子电池组系统内电池单体的原位热仿真分析方法 | |
Chai et al. | Fast ${\cal H} $-Matrix-Based Direct Integral Equation Solver With Reduced Computational Cost for Large-Scale Interconnect Extraction | |
Wang et al. | Spatial correlation-based incremental learning for spatiotemporal modeling of battery thermal process | |
Khalik et al. | Model simplifications and their impact on computational complexity for an electrochemistry-based battery modeling toolbox | |
Liu et al. | An efficient discontinuous Galerkin finite element method for highly accurate solution of Maxwell equations | |
CN109684740B (zh) | 一种基于混合网格及时间步长的电磁学多尺度计算方法 | |
CN114510892A (zh) | 一种模拟mosfet器件总剂量效应的高性能仿真方法 | |
CN113128159B (zh) | 混合流线迎风有限体积有限元方法、模型数值离散系统 | |
WO2024041642A1 (zh) | 时域热传导仿真方法及存储介质 | |
Dolean et al. | A domain decomposition method for solving the three-dimensional time-harmonic Maxwell equations discretized by discontinuous Galerkin methods | |
CN107908900A (zh) | 一种基于凸模型不确定性的连续体结构非概率可靠性拓扑优化方法 | |
CN113591346A (zh) | 一种用于电磁场求解的自适应区域分解有限元方法 | |
Gui et al. | Improved modal analyses using the novel quadrilateral overlapping elements | |
Lou et al. | A new explicit time-domain finite-element method based on element-level decomposition | |
Guo et al. | Nonlinear state-variable method for solving physics-based li-ion cell model with high-frequency inputs | |
CN112347638B (zh) | 一种基于对偶单元法的三维集成微系统电热耦合分析方法 | |
Abdi | Evolutionary topology optimization of continuum structures using X-FEM and isovalues of structural performance | |
CN105404751A (zh) | 基于热-力-电磁场网络统一的实现方法 | |
CN112836399A (zh) | 一种基于有限元算法的非线性接触热阻热分析求解方法 | |
JP4233513B2 (ja) | 解析装置、解析プログラム、および解析プログラムを記録したコンピュータ読取可能な記録媒体 | |
Xue et al. | Efficient transient thermal analysis based on spectral element time domain method with curvilinear hexahedrons | |
Liu et al. | Layered ${\cal H} $-Matrix Based Inverse and LU Algorithms for Fast Direct Finite-Element-Based Computation of Electromagnetic Problems |
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 |