第五章 3D静态磁场分析(标量法)
ANSYS电磁场分析报告指南设计
ANSYS电磁场分析指南（共17章）ANSYS电磁场分析指南第一章磁场分析概述：ANSYS电磁场分析指南第二章 2-D静态磁场分析：ANSYS电磁场分析指南第三章２-Ｄ谐波（ＡＣ）磁场分析：ANSYS电磁场分析指南第四章２-Ｄ瞬态磁场分析：ANSYS电磁场分析指南第五章３-Ｄ静态磁场分析（标量法）：ANSYS电磁场分析指南第六章３-Ｄ静态磁场分析（棱边元方法）：ANSYS电磁场分析指南第七章３-Ｄ谐波磁场分析（棱边单元法）：ANSYS电磁场分析指南第八章３-D瞬态磁场分析（棱边单元法）：ANSYS电磁场分析指南第九章 3-D静态、谐波和瞬态分析（节点法）：ANSYS电磁场分析指南第十章高频电磁场分析：ANSYS电磁场分析指南第十一章磁宏：ANSYS电磁场分析指南第十二章远场单元：ANSYS电磁场分析指南第十三章电场分析：ANSYS电磁场分析指南第十四章静电场分析（h方法）：ANSYS电磁场分析指南第十五章静电场分析（P方法）：ANSYS电磁场分析指南第十六章电路分析：ANSYS电磁场分析指南第十七章其它分析选项和求解方法：第一章磁场分析概述1.1磁场分析对象利用ANSYS/Emag或ANSYS/Multiphysics模块中的电磁场分析功能，ANSYS可分析计算下列的设备中的电磁场，如：·电力发电机·磁带及磁盘驱动器·变压器·波导·螺线管传动器·谐振腔·电动机·连接器·磁成像系统·天线辐射·图像显示设备传感器·滤波器·回旋加速器在一般电磁场分析中关心的典型的物理量为：·磁通密度·能量损耗·磁场强度·磁漏·磁力及磁矩· S-参数·阻抗·品质因子Q·电感·回波损耗·涡流·本征频率存在电流、永磁体和外加场都会激励起需要分析的磁场。
ANSYS热耦合-学习笔记
＝＝＝＝＝【热力耦合分析单元简介】＝＝＝＝＝＝SOLID5－三维耦合场实体具有三维磁场、温度场、电场、压电场和结构场之间有限耦合的功能。
本单元由8个节点定义，每个节点有6个自由度。
在静态磁场分析中，可以使用标量势公式（对于简化的RSP，微分的DSP，通用的GSP）。
在结构和压电分析中，具有大变形的应力钢化功能。
与其相似的耦合场单元有PLANE13、SOLID62和SOLID98。
INFIN9－二维无限边界用于模拟一个二维无界问题的开放边界。
具有两个节点，每个节点上带有磁向量势或温度自由度。
所依附的单元类型可以为PLANE13和PLANE53磁单元，或PLANE55和PLANE77和PLANE35热单元。
使用磁自由度（ＡＺ）时，分析可以是线性的也可以是非线性的，静态的或动态的。
使用热自由度时，只能进行线性稳态分析。
PLANE13－二维耦合场实体具有二维磁场、温度场、电场和结构场之间有限耦合的功能。
由4个节点定义，每个节点可达到4个自由度。
具有非线性磁场功能，可用于模拟B－H曲线和永久磁铁去磁曲线。
具有大变形和应力钢化功能。
当用于纯结构分析时，具有大变形功能，相似的耦合场单元有SOLID5、SOLID98和SOLID62。
LINK31－辐射线单元用于模拟空间两点间辐射热流率的单轴单元。
每个节点有一个自由度。
可用于二维（平面或轴对称）或三维的、稳态的或瞬态的热分析问题。
允许形状因子和面积分别乘以温度的经验公式是有效的。
发射率可与温度相关。
如果包含热辐射单元的模型还需要进行结构分析，辐射单元应当被一个等效的或（空）结构单元所代替。
LINK32－二维传导杆用于两节点间热传导的单轴单元。
该单元每个节点只有一个温度自由度。
可用于二维（平面或轴对称）稳态或瞬态的热分析问题。
如果包含热传导杆单元的模型还需进行结构分析，该单元可被一个等效的结构单元所代替。
LINK33－三维传导杆用于节点间热传导的单轴单元。
该单元每个节点只有一个温度自由度。
静态电磁场
2. 电位的表达式 对于连续的体分布电荷，由
R r r
E(r )
1
4π
V
(r) R3
RdV
1
4π
V
(r)( 1 )dV
R
[ 1
4π
V
(r)( 1 )dV ]
R
故得
(r )
1
4π
V
(r)dV C
R
1
R
( R) R3
同理得，面电荷的电位： (r ) 1
4π
S
S (r R3
) dS
C
线电荷的电位：
应用高斯定理可得到内外导体间任一点的电场强度为
E()
e
l 2π
内外导体间的电位差
U
b a
E( ) ed
l 2π
b 1d a
ab
同轴线
l ln(b / a)
2π
故得同轴线单位长度的电容为
C1
l
U
2π
ln(b / a)
(F/m)
＊ 2. 部份电容
在多导体系统中，任何两个导体间的电压都要受到其余导体 上的电荷的影响。因此，研究多导体系统时，必须把电容的 概念加以推广，引入部分电容的概念。
由于所以介质片所受到的力有将其拉进电容器的趋势此题也可用式来计算q不变设极板上保持总电荷q不变则由此可得由于同样得到32本节内容321恒定电场的基本方程和边界条件322恒定电场与静电场的比拟323漏电导可知导体中若存在恒定电流则必有维持该电流的电场虽然导体中产生电场的电荷作定向运动但导体中的电荷分布是一种不随时间变化的恒定分布这种恒定分布电荷产生的电场称为恒定电场
Δl
2
2
n
静态电磁场及其边值问题的解chap3
ϕ ( P) = ∫
∞
P
v v E ⋅d l
（以无限远处为零电位做参考，任意点P的电位表示） 以无限远处为零电位做参考，任意点P的电位表示）
2、静电位的微分方程
v E = −∇ϕ ⇒ D = ε E = −ε∇ϕ
∇ ⋅ D = ρ ⇒ ∇⋅ ( −ε∇ϕ ) = ρ ⇒ ∇⋅ ( ∇ϕ ) = − ρ ⇒ ∇2ϕ = − ρ
ρS = 0
∂ϕ1 ∂ϕ2 ε1 =ε2 ∂n ∂n
导体
∂ϕ ε =−ρS ∂n
【例3.1.1】 求电偶极子 p = qdl 的电位 ϕ ( r ) 3.1.1】
当
z
+q d
r+ r− = r
P ( r,θ ,ϕ )
r >> d 1 1 1 ≈ + 2 d cosθ r+ r r
因此
ϕ=
θ
−q
解：取如图所示坐标系，场点 P ( r,θ ,ϕ ) 取如图所示坐标系， 的电位等于两个点电荷电位的叠加
a a
Cl =
ρl
U
=
D−a ln a
πε 0
≈
πε 0
ln( D / a)
ρl 1 ρl D − a 1 ( + )dx = ln 2πε 0 x D − x πε 0 a
F /m
【例3.1.5】同轴线的内导体的半径为a，外导 3.1.5】同轴线的内导体的半径为a 体的半径为b 体的半径为b，内外导体间填充介电常数为 ε 的均匀电介质，试求同轴线单位长度的电容。 的均匀电介质，试求同轴线单位长度的电容。
电 位 的 泊 松 方 程
ε
ε
若空间电荷分布为零， 若空间电荷分布为零，则有 ∇2ϕ = 0
ANSYS电磁场分析指南
ANSYS电磁场分析指南（共17章）ANSYS电磁场分析指南第一章磁场分析概述：ANSYS电磁场分析指南第二章2-D静态磁场分析：ANSYS电磁场分析指南第三章２-Ｄ谐波（ＡＣ）磁场分析：ANSYS电磁场分析指南第四章２-Ｄ瞬态磁场分析：ANSYS电磁场分析指南第五章３-Ｄ静态磁场分析（标量法）：ANSYS电磁场分析指南第六章３-Ｄ静态磁场分析（棱边元方法）：ANSYS电磁场分析指南第七章３-Ｄ谐波磁场分析（棱边单元法）：ANSYS电磁场分析指南第八章３-D瞬态磁场分析（棱边单元法）：ANSYS电磁场分析指南第九章3-D静态、谐波和瞬态分析（节点法）：ANSYS电磁场分析指南第十章高频电磁场分析：ANSYS电磁场分析指南第十一章磁宏：ANSYS电磁场分析指南第十二章远场单元：ANSYS电磁场分析指南第十三章电场分析：ANSYS电磁场分析指南第十四章静电场分析（h方法）：ANSYS电磁场分析指南第十五章静电场分析（P方法）：ANSYS电磁场分析指南第十六章电路分析：ANSYS电磁场分析指南第十七章其它分析选项和求解方法：第一章磁场分析概述1.1磁场分析对象利用ANSYS/Emag或ANSYS/Multiphysics模块中的电磁场分析功能，ANSYS可分析计算下列的设备中的电磁场，如：·电力发电机·磁带及磁盘驱动器·变压器·波导·螺线管传动器·谐振腔·电动机·连接器·磁成像系统·天线辐射·图像显示设备传感器·滤波器·回旋加速器在一般电磁场分析中关心的典型的物理量为：·磁通密度·能量损耗·磁场强度·磁漏·磁力及磁矩·S-参数·阻抗·品质因子Q·电感·回波损耗·涡流·本征频率存在电流、永磁体和外加场都会激励起需要分析的磁场。
电磁场与电磁波课件第5章 静态场的边值问题
1 2 ，
然后进行 证明.同样可得出结论,其解唯一.
设φ1φ2是同一有源区域的边值问题
2 的解。 | f1 ( S )
即在区域V内，φ1和φ2满泊松方程，即
1 2 2
2
在V的边界S上，φ1和φ2满足同样的边界条件， 即
5.3.1 导体平面镜像
设在无限大导体平面(z=0)附近有一点电荷与平面距离为z=h 。 若导体平面接地，则导体平面电位为零，如图所示。求上半 空间中的电场。 分析：上半空间任一点 P处的电位，应等于点 电荷q和无限大导体平 板上感应的负电荷产生 的的电位总和。因此， 上半空间的电位问题可 表示为 ：
2
C (常数)
0
1 2
C 0
5.3 镜像法
实质:是以一个或几个等效电荷代替边界的影响，将原来具有边
界的非均匀空间变成无限大的均匀自由空间，从而使计算过程 大为简化。
依据：惟一性定理。等效电荷的引入必须维持原来的边界 条件不变。这些等效电荷通常处于镜像位置，因此称为镜 像电荷，而这种方法称为镜像法。
2 A ( A) A J
人为规定
A 0
这个规定被称为库仑规范
于是有
2 A J
此式即为矢量磁位的泊松方程。
在没有电流的区域有J 0
2 A0
此式即为矢量磁位的拉普拉斯方程。 (2) 磁场的标量位函数 在没有电流分布的区域内，恒定磁场的基本方程变为 H 0 B 0 这样，在无源区域内，磁场也成了无旋场，具有位场的性 质，因此，象静电场一样，我们可以引入一个标量函数， 即标量磁位函数
磁场有限元分析
电机磁场分析主要关注电机的磁通密度、磁通路径、磁阻、涡流和磁力线分布等参数。通过有限元方 法，可以模拟电机的磁场分布和变化，从而优化电机设计，提高电机的功率密度、效率和使用寿命。
磁悬浮系统分析
总结词
磁悬浮系统分析是磁场有限元分析的重要应用之一，通过分析磁悬浮系统的磁场分布和作用力，可以优化磁悬浮 系统的控制和稳定性。
磁场有限元分析
contents
目录
• 引言 • 磁场有限元分析的基本原理 • 磁场有限元分析的实现过程 • 磁场有限元分析的应用案例 • 磁场有限元分析的挑战与展望 • 参考文献
01 引言
背景介绍
磁场有限元分析是计算电磁场问题的一种数值方法，通过将 连续的磁场分布离散化为有限个小的单元，利用数学模型和 物理定律建立每个单元的方程，然后通过求解这些方程得到 磁场的近似解。
在磁场有限元分析中，复杂的边 界条件（如开域、闭域、周期性 边界等）需要特殊处理，以确保 求解的准确性和可靠性。
材料属性
不同材料的磁导率、磁化强度等 属性可能存在较大的差异，需要 在模型中准确描述，以便更准确 地模拟磁场分布和磁力作用。
多物理场耦合的磁场有限元分析
耦合方式
多物理场耦合的磁场有限元分析需要 考虑磁场与其他物理场（如电场、流 体场等）之间的相互作用和耦合效应， 需要采用适当的耦合方式进行建模和 分析。
结果后处理
结果可视化
将计算结果以图形或图像的形式呈现出来，便 于观察和分析。
结果评估
对计算结果进行评估，判断其准确性和可靠性。
结果优化
根据需要对计算结果进行优化处理，如滤波、平滑等。
04 磁场有限元分析的应用案 例
电机磁场分析
总结词
第五章 有限元法-1-泛函与变分
设待求变分问题（5-4）的解答（极值函数）为 y=y(x) （5-7）
因y是x的函数，但讨论的是y的变化
设想函数y从极值解（5-7）稍稍变动到y+dy，并把变分dy改记为：eh(x)，
e是一个任意给定的微量实参数（实变量）；
h(x)是定义于区间[x1，x2]，且满足齐次边界条件的任意选定的可微函数，即有： h(x1)=h(x2)=0。
15
与多元函数的极值问题相对应，在几何、力学上的求解泛 函极值的问题。 最速降线问题。
研究当质点从定点A自由下滑到定点B时，为使滑行时间最短，试 求质点应沿着怎样形状的光滑轨道y=y(x)下滑。 取A点为坐标原点，y轴竖直向下（图5-1）。
则沿曲线y=y(x)滑行线段ds所需的时间为
16
18
在最速下降问题，在端点x1和x2给定的无数个函数之中， y ( x) 仅有一个函数 能使式（ 5-2a）中的定积分达到极小 y ( x) 值函数，这一函数 被称为极值函数。 所谓变分问题就在于寻求使泛函达到极值的该极值函数， 即分析研究泛函的极值问题。 物理学各分支都存在有相应的变分问题（变分原理），例 如
因此
式中
26
故可得
简写为
将上式与式（5-6）相比较，只相差一个数值因子e。
27
故（5-8）等价于变分方程
也即
（线性主部）
利用分部积分，根据变分与微分顺序可以互换的原理，即 dy’=(dy)’，得
28
在变分问题中，变分dy在端点保持为零
于是，必要条件（5-12）成为
电流域标量处理法在横向磁场电机三维磁场分析中的应用
Cartesian coordinates 设 A 磁位是连续的，沿线圈右侧面磁位差为
∫ ∆ψ = l Jsdl = H1 • y
（21）
2.2 圆柱坐标系下电流区域处理方法 2.2.1 圆柱坐标系下任意截面形状的电流区域处
理方法
图 4 所示，Ω2 为圆柱坐标系下任意形状的电流 区域在 r−θ 平面内的截面。各量参考方向如图所示。
（22）
eθ —— θ 向单位矢量 设 C 为电流区域边界上任意一点，其半径为 r，
其法向与 A 点半径方向夹角为α，在该点设置面磁
荷密度
ρs=µH1·sinα 同时，该点设置面电流密度
（23）
Js= H1·cosα
（24）
在电流区域边界同样存在磁位差。设 A 点磁位
差∆ψ =0，沿边界逆时针方向不同位置的磁位差为
• ab +
J
•a•(y
− b)
y ∈[b,b + c] （16）
BC 面
∆ψ P
=
J
•a•b+
1 2
J
•a•c−
1 2
J
• tanα
•
x2
x ∈[a,0]
三维磁场有限元分析技术
第39卷第2期原子能科学技术Vol.39,No.2 2005年3月Atomic Energy Science and TechnologyMar.2005100Me V 紧凑型回旋加速器主磁铁的三维磁场有限元分析技术钟俊晴,张天爵,杨建俊,储诚节(中国原子能科学研究院核技术应用研究所,北京 102413)摘要:100MeV 紧凑型回旋加速器主磁铁的几何结构十分复杂,但为了形成加速器束流动力学所要求的磁场分布,本文对初步设计的磁铁进行必要的简化。
综合采用各种适当的三维有限元网格剖分技术,对该磁铁的磁场进行数值分析,计算精度满足加速器物理设计的要求。
关键词:紧凑型回旋加速器;有限元法;剖分;双标量位法中图分类号:TL542 文献标识码:A 文章编号:100026931(2005)022*******3D Finite Element Analysis Method for Main Magnet Design of 100Me V Compact CyclotronZHON G J un 2qing ,ZHAN G Tian 2jue ,YAN G Jian 2jun ,C HU Cheng 2jie(China I nstitute of A tomic Energ y ,P.O.B ox 27523,B ei j ing 102413,China )Abstract : To get t he magnetic field distribution for beam dynamics calculation in t he 100MeV co mpact cyclot ron feat ured wit h a complex geomet ric magnet ,we simplified t he magnet struct ure and used various 3D finite element meshing met hods in t he initial design stage.After t he result s comparison by different meshing met hods ,we analysed numerically t he t hree dimentional magnetic field and found t hat t his 3D finite element calculation p rocess could meet t he requirement of t he p hysics design for 100MeV com 2pact cyclot ron.K ey w ords :compact cyclot ron ;finite element met hod ;meshing ;double scalar potential收稿日期:2004211225基金项目:国家自然科学基金资助项目(10125218)作者简介:钟俊晴(1978—),男,江西瑞金人,实习研究员,加速器专业 串列加速器升级工程主要由以下几部分组成:100MeV 回旋加速器、在线同位素分离器、超导直线增能器以及现有HI 213串列加速器注入器的改造。
3静态场分析法
∫ J ⋅ dS = 0
S
或
∇⋅ J = 0
电流恒定时，电荷分布不随时间变化， ◇ 电流恒定时，电荷分布不随时间变化，恒定电场同静电场 具有相同的性质， 具有相同的性质，因此
∫ E ⋅ dl = 0
l
或
∇× E = 0
P.9
一、恒定电场的电位
与静电场的讨论类似， 与静电场的讨论类似，由 ∇× E = 0 可引入恒定电场的电位函数
E = −∇ϕ
由
J = σE J =γ E
∇⋅J = 0
2 ∇ϕ =0
二、恒定电场的边界条件
将恒定电场的基本方程
∫ J ⋅ dS = 0、∫ E ⋅ dl = 0 分别用于二种不同导电媒
S l
质的分界面上，与推导静电场边界条件方法类似，可导出恒定电场的边界条件。 质的分界面上，与推导静电场边界条件方法类似，可导出恒定电场的边界条件。
1 ρl dl ,ϕ = ∫ R +C 4πε0 l
P.5
电位的边界条件： 电位的边界条件： 对于电位 ϕ 由
E1t = E2t
∂ϕ1 ∂ϕ 2 ⇒− =− ∂τ ∂τ
⇒ ϕ1 = ϕ2
由 D1n − D2 n = ρ s
⇒ ε1 E1n − ε 2 E2 n = ρ s
∂ϕ1 ∂ϕ2 − ε2 = −ρs ⇒ ε1 ∂n ∂n
ϕ1 = Aln r + B
P.4
∫
El ⋅ dl
对于点电荷的电场， ◇ 对于点电荷的电场， 处的电位为零， 若取 R P → ∞ 处的电位为零，则
q ϕ= 4πε0 R
◇ 体电荷 ρdτ 、面电荷 σdS 、线电荷 ρldl 产生的电位 分别为
三维电磁场矢量有限元分析
单元刚度矩阵的计算
F i j d
e ij
f xx Fije 0 0
0 f yy 0
0 0 f zz
4 2 xyz f xx 36 2 1
2 4 1 2
2 1 4 2
1 2 2 4
总刚度矩阵的计算
里兹(Ritz)变分方法 LФ=f
伽辽金(Galerkin)方法
是边值问题的近似解，得到非零 假设 的残数： f 0 r L
残数加权方法类型，正如其名称所指，它通过对微 分方程的残数求加权方法来得到方程的解。
在伽辽金方法中，加权函数与近似解展开中所用的函数 相同。
最佳近似应能使残数R在Ω内所有点上 R有最小值。残数加权方法要求 Ri wi rd 0
exy eyy ezy
exz e yz ezz
T T exx x x x x d z z y y T T y y y y eyy d z z x x T T ezz z z z z d y y x x
s
齐次狄利克雷(Dirichlet)条件
ˆE 0 n
2 E kt d Jd i j i j j
正六面体单元插值形函数的表达式
1 x
x2
3 x
1 y z y y z z c c yz 2 2 1 y z y y z z c c yz 2 2
磁场有限元分析步骤
1 绪论1.1磁流变液概述1947年，美国学者WinsfowW.M.发明电流变材料。
第二年，美国学者Rabinow，J.[1]发明了磁流变材料应用装置(离合器)就出现了。
然而在随后的研究中，人们首先对电流变材料及其应用的研究给予了极大关注，各种不同母液和悬浮微粒的电流变材料和一些电流变器件[2,3]相继研究成功。
在这个时期，磁流变技术一直处于停滞不前的状态，很少引起人们对此领域的关注，直到上世纪九十年代，研究人员发现电流变材料屈服应力较低，且存在高压安全性问题等一系列难以解决技术问题，所以磁流变液又重新引起了研究者们的兴趣。
磁流变液(Magnetorheological fluid，即MRF)作为一种新型智能材料，主要由高磁导率、低磁滞性的微小铁磁性颗粒、非导磁性载液及稳定剂组成。
载液的作用是将固体粒子均匀地分散开，这种分散作用能保证在零磁场时，使磁流变液仍保持有牛顿流体的特性，而在有磁场作用时，则使粒子形成链化的结构，产生剪切屈服应力，并使磁流变液呈现Bingham弹塑性流体的特性。
好的载液应具有低零场粘度，大的温度稳定性，不污染环境等特性。
铁磁性微粒大多数由铁、钻、镍等磁性材料组成，在磁场的作用下粒子形成链状，产生磁流变效应。
因此，固体粒子材料的化学性质和物理性质，对磁流变液的性能优劣起着决定性的作用。
固体粒子一般用球形金属及铁氧体磁性材料，其尺寸范围为5-50μm，所占体积比为15%-30%。
稳定剂又称表面活性剂，其作用是促进磁流变效应，防止铁磁性微粒凝聚及沉淀。
按作用来说，有沉降稳定剂和悬浮稳定剂，最常见的稳定剂有油酸盐、山梨醇、甘油、二乙胺、低分子量的聚酞胺、琥珀酸酰亚胺、聚二甲基硅氧烷、苯酸盐、磺酸盐和磷酸盐等，其在磁流变液中的含量较少，一般低于5%。
优质磁流变液具有的性能特征是;沉淀稳定性好、易于再分散、动屈服应力高、零场粘度低、响应时间快、工作温度范围宽。
Phule[4]利用中等尺寸炭基铁或镍锌铁氧体的颗粒、极性有机载液与纳米级胶体和其他添加剂制备了相对稳定、易于再分散且力学性能优异的作了适当改进。
ANSYS教程：ANSYS电磁场分析
ANSYS教程：ANSYS电磁场分析静态磁场分析：用于分析不随时间变化的磁场，主要包括三类情况：用磁场的磁场，稳恒电流产生的磁场，匀速运动的导体所产生的磁场。
对于三位静态磁场分析，ansys程序采用了两种方法：标量势法（scalar method）和单元边法（edge-based-method），其中标量势法根据其标量势方程的不同又可分为三种不同的标量势分析方法：简化标量势法（RSP）、微分标量势法（DSP）和广义标量势法（GSP）。
使用单元边法时，电流源是作为整个系统的一部分一起进行网格划分的，由此使用该方法不仅能计算常规物流量（如磁场、磁动势等），还能计算诸如焦耳热损、洛伦兹力等。
根据以下原则选择不同的分析方法：当所分析的问题中不含铁芯区域或虽含铁芯区域但不含电流源时，采用RSP法，在含有铁芯和电流源的模型分析中通常不使用RSP 法。
对于“单连通”铁芯区域模型，使用DSP法，对于“多连通”铁芯区域模型，使用GSP法。
单连通区域指的是带有空气隙的磁路不封闭的铁芯系统，没有空气隙的则为磁路封闭多连通铁芯区域系统。
对于非连续介质模型一般采用单元边法进行求解。
提示：单元边法中使用的单元的节点自由度矢量磁势是沿单元边切向积分的结果，其求解精度高于标量势法的求解精度。
单元边法不仅适用于三维静态磁场分析中，也适用于三维谐性和瞬态磁场分析中。
1 电磁场分析中的默认单位制为MKS单位制，即米、安培和秒。
可以定义其他的单位制：main menu/preprocessor/material props/electromag units2 电磁场分析中大多材料的磁性能可以从ansys程序的材料库中读入，用于也可以自己定义材料性能，方法如下：2.1 定义路径main menu/preprocessor/material props/material library/library path2.2 读入材料参数main menu/preprocessor/material props/material library/import librarymain menu/preprocessor/loads/load step opts/change mat props2.3 修正材料参数main menu/preprocessor/material props/material library/export library2.4 定义材料B-H曲线main menu/preprocessor/material props/material models/electomagnetics/BH curve2.5 在模型上施加电流密度载荷main menu/preprocessor/loads/define loads/apply/magnetic/excitation/current density/on elements2.6 施加电压载荷main menu/preprocessor/loads/define loads/apply/magnetic/excitation/volt drop/on elements2.7 进行求解main menu/solution/solve/electromagnet/static analysis/opt&solv2.8 退出求解器main menu/finish谐性磁场分析：用于分析激励源按正弦或余弦规律变化的磁场问题，如变压器、感应式电机，感应加热炉等电磁装置引发的磁场均属于谐性磁场问题。
ANSYS三维磁场分析
δ = (π μ σ f) -1/2 (m)
式中
μ = 磁导率 = μr μ0
=1.256E-6
σ = 电导率= 1 /ρ = 5.83E-8
f = 场频率 (Hz)=1000
δ=2.1 mm
第4页/共23页
4.4-3
• 在内导体对称平面加VOLT 约束 • 施加 VOLT约束
Solu>apply>-electric-boundary>on areas
• 在模型末端( Z 最大值位置）加通 量平行边界条件作为一个远程边 界条件
• 面上加边界条件 • 使整个模型激活 • 选择外表面 • 不选Z=0的平面
第10页/共23页
4.4-9
• 加通量平行条件 Preproc>loads>apply>magnetic-boundary>-flux par’l-on areas
δ = (π μ σ f) -1/2 (m) 式中 μ = 磁导体= μr μ0 =50*1.256E-6 σ = 电导率 = 1 / ρ ρ = 电阻率 (Ohm-m) =1E-7 f = 场频率(Hz)=1000 δ = .77 mm
第3页/共23页
BAR的边界
4.4-2
外导体的内半径
• 由于集肤效应，外导体内半径也要网 格细化
• 在命令窗口输入ac3d建模 • 模型信息
• 单元类型: 边单元solid117 • 导体- AZ,VOLT自由度 • 空气-AZ 自由度
• 励磁导体：导体要建模 • 总电流加到截面的一个节点上
• 单元组件: • 内导体: BAR • 外导体: ELEMENT
第2页/共23页
4.4-1
• 靠近杆BAR外半径处要加密网格，以 满足集肤效应的计算要求（集肤厚度 (δ)上要有两个单元）
第五章 时变电磁场
解：1、 I J dS 2 10r 1.5 r 2 sin d d
S
00
40 r0.5
3.9738A
r 1mm
2、因为
J
1 r2
d dr
r 2 10r 1.5
dS
H dS
S
上式右边应用散度定理可以写为
S H dS V H dV 0
左边为
D
S
J
c
t
dS
Ic
Id
I
0
证毕
例5-3 坐标原点附近区域内传导电流为 J er 10r 1.5( A / m2 ) 试求：1、通过半径 r = 1mm的球面的电流值；
B
E
l
dl
S
t
dS
B
S
dS
0
D
S
dS
q
微分形式 H J D
t E B
t B 0
D
可见，时变电场是有旋有散的，时变磁场是有旋无散的。但是， 时变电磁场中的电场与磁场是不可分割的，因此，时变电磁场是有旋 有散场。
四、麦克斯韦方程组的辅助方程—本构关系 》一般媒质本构关系 》各向同性线性媒质本构关系
D B
0E 0 ( H
P M
)
J
E
D E
(整理)ANSYS电磁场分析指南第五章3-D静态磁场分析标量法.
第五章３-Ｄ静态磁场分析（标量法）5.1 在3-D静态磁场分析（标量法）中要用到的单元表1三维实体单元：表2三维界面单元表3三维连接单元表4三维远场单元SOLID96和SOLID97是磁场分析专用单元，SOLID62、SOLID5和SOLID98更适合于耦合场求解。
5.2 磁标量位（MSP）法介绍在磁标量位方法中，可使用三种不同的分析方法：简化标势法(RSP)、差分标势法(DSP)和通用标势(GSP)法。
·若模型中不包含铁区，或有铁区但无电流源时，用RSP法。
若模型中既有铁区又有电流源时，就不能用这种方法。
·若不适用RSP法，就选择DSP法或GSP法。
DSP法适用于单连通铁区，GSP 法适用于多连通铁区。
5.2.1 单连通区与多连通区单连通铁区是指不能为电流源所产生的磁通量提供闭合回路的铁区，而多连通铁区则可以构成闭合回路。
参见图1(a)、(b)“连通域”。
数学上，通过安培定律来判断单连通区或是多连通区，即磁场强度沿闭合回路的积分等于包围的电流（或是电动势降MMF）。
因为铁的磁导率非常大，所以在单连通区域中的MMF降接近于零，几乎全部的MMF降都发生在空气隙中。
但在多连通区域中，无论铁的磁导率如何，所有的MMF降都发生在铁芯中。
5.3 3-D静态磁标势分析的步骤该分析类型与2－D静态分析的步骤基本一样：1.建立物理环境2.建模、给模型区域赋属性和分网格3.加边界条件和载荷（激励）4.用RSP、DSP或GSP方法求解5.观察结果5.3.1创建物理环境首先设置分析参数为“Magnetic-Nodal”，并给出分析题目。
然后用ANSYS 前处理器定义物理环境包含的项目。
即单元类型、KEYOPT选项、材料特性等。
3D分析的大部分过程与2D分析一致，本章下面部分介绍3D分析中要特殊注意的事项。
· SOLID96单元可为模型所有的内部区域建模，包括：饱和区、永磁区和空气区（自由空间）。
ansys三维静态分析
[分享]ansys的help文件中“5.3. 三维静态标量分析”一段翻译5.3. 三维静态标量分析5.3.1创造物理环境对于一个3维静态分析,要记住以下几点:1、可以用SOLID96单元来代表你的模型所有的内部区域:可饱和区、永磁区、空气（自由空间）。
对于电流导通区域，用SOURC36电流源单元，可以参考Modeling Current Conduction Regions获取更多的详细情况。
2、在空气单元外部，推荐用4节点边界单元的INFIN47或者8节点和20节点边界单元。
这两种单元可用来建立磁场远场区衰减模型，通常比平行通量或垂直通量条件下得到更好的结果。
而INFIN111单元建立空气外部空间更精确。
3、默认的单位制是MKS即米、千克、秒。
你可以通过下面的方法换算成其他的单位制。
一旦你选定了一种单位制，所有输入的数据都以这种单位制为准。
为了方便起见，你可以用其他单位制创建几何模型，例如毫米、英寸，然后按比例放大到适合的单位制里。
命令：EMUNITGUI: Main Menu> Preprocessor> Material Props> Electromag Unitsl 根据你输入的单位，空气的磁导率就被自动确定了：μo = 4 π* 10-7 H/m 是按国际单位制确定的。
5.3.3.指定材料属性模型可能会有以下的材料区：空气，磁导材料，载流区，或永磁体。
每种材料区域都需要材料属性。
ANSYS材料库包含了一些磁性材料的定义，你可以把这些材料属性写入ANSYS数据库，必要时进行修改以更接近地匹配分析中的材料。
材料包含材料定义的属性文件铜 emagCopper.SI_MPLM3钢 emagM3.SI_MPLM54钢 emagM54.SI_MPLSA1010 钢 emagSa1010.SI_MPL硅钢 emagSilicon.SI_MPL铁钴钒钢 emagVanad.SI_MPL铜的属性由随温度变化的电阻率和相对磁导率表示。
第五章 三维恒定磁场分析 Ansys工程电磁场有限元分析 华科电气
第五章 三维恒定磁场分析（3D Static Magnetic Field Analysis ）5.1 基本理论5.1.1 恒定磁场基本方程：∇⋅=∇⨯=B H J5.1.2 矢量磁位公式：=∇⨯B A1μ∇⨯∇⨯=A J(1) 上述双旋度方程的解是不唯一的。
引入库仑规范0∇⋅=A(2) 如果μ为常数，利用2()∇⨯∇⨯=∇∇⋅-∇A A A 得：2μ∇=-A J(3) 上述三个方程的关系：（1）+（2）=（3）在合适的边界条件下，方程(3)的解是唯一的，而方程(1)的解是不唯一的。
但吊诡的是，方程（1）居然也可以解；不但可以解，而且结果比（3）还好！A的矢量泊松方程和双旋度方程两种表述，从数理方程的角度看，矢量泊松方程表述比较严格，因为divA规定为0，满足库仑规范，因此A的解答唯一。
双旋度方程未加库仑规范，所以存在A不唯一的问题。
但是在适当的边界条件下反而是双旋度方程不但有确定的数值解，而且B与实测值也吻合。
……这里的问题是，为何不加divA=0的条件，双旋度方程也能得到确定的数值解？作者认为，原来电磁场方程的解答是一系列调和函数和超越函数的叠加，解函数具有较高的光滑性。
在有限元离散的过程中，方程的维数先从无限维降为有限维，接着在一阶四面体单元中，又用线性函数作为单元内的插值函数，降低了对解光滑性的要求，结果使原来的微分方程转化为N维的线性代数方程组，相应地，微分方程解答的唯一性问题，也转化为联立线性代数方程组的可解性问题，于是在特定的边界条件下，双旋度方程将有确定的数值解。
——汤蕴璆，梁艳萍. 电机电磁场的分析与计算[M]. 北京：机械工业出版社，2010结论：在节点有限元中，用A 计算三维磁场，计算量大，效果不佳，缺乏吸引力。
类似的问题夜存在于三维涡流场的计算中。
这个问题的解决是棱边有限元的建立。
——但是，棱边有限元也有他的问题。
5.1.3 标量磁位公式0, ∇⋅=∇⨯=B H J令 m s =+H H H其中：m , 0s ∇⨯=∇⨯=H J Hs H —单纯由电流产生（即：不存在磁介质时）的磁场，由毕奥沙伐定律计算；可视作已知量。
5.3 三维静磁场的有限元分析
5.3 三维静磁场的有限元分析5.3.1 边值问题以标量磁位m ϕ表示的无源区磁场的边值问题与电位的拉普拉斯边值问题的数学表达形式完全一样，可以如前节所述的有限元分析。
在此，考虑有电流存在以矢量磁位A 作为待求变量的有限元分析。
设在线性媒质中，磁场满足的边界条件：边界1S 面上有0A A =，在边界S 2面上取某种形式的对称面作为第二类齐次边界，在该面上磁场强度H 的切向分量为零：()0=⨯⨯∇=⨯n m n e A e H γ，有如下边值问题：()()⎪⎩⎪⎨⎧∈=⨯⨯∇∈=∈=⨯∇⨯∇210s s V n mm e A A A J A γγ5.3.2 场域剖分与插值对于求解场域V ，根据其形状和场的定性分布，选择合适的单元(例如四面体单元)，进行场域剖分，得到0Z 个单元、0N 个节点。
在单元e 内，对位函数A 进行插值。
若采用四面体单元：∑==41j je j N A A ~式中ej N 是单元形状函数，分量形式z j zje j y j yj ej x j xj ej zz y y x x A N A N A N A A A e e e e ~e ~e ~A ~⎪⎪⎭⎫⎝⎛+⎪⎪⎭⎫⎝⎛+⎪⎪⎭⎫ ⎝⎛=++=∑∑∑===414141 以矩阵表示磁矢量位A 在单元节点上的各分量[][]()z y x l A A A A A Tl l l l e l ,,==4321),,(][][~41z y x l A N A N A el T e j lj e j l ===∑=于是，磁矢量位z z y y x x A A A e ~e ~e ~A ~++=的矩阵表达式∴[][][][][][][][][][][][]e Te z x Te T eT ee z T ee y T ee x Tez y x A A A N N N A N A N A N A A A A 041000000N ~~~~=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=式中：Tz z y y x x x e A A A A A A A A ][][4141421 =[][][][]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=e ee e N N N 000000N 5.3.3 单元分析将近似函数A ~代入方程有余量()J A ~R -⨯∇⨯∇=m γ 取矢量权函数为W i ，有加权余量方程0>=<R ,W i单元剖分后余量加权的和式为零()[]{}∑∑⎰===∑=-⨯∇⨯∇⋅>=<0110Z e Z e i VemieidV εγJ A W R,W由矢量恒等式()()()B A A B B A ⨯∇⋅-⨯∇⋅=⨯⋅∇上式分解为：d()()()()()()()⎰⎰⎰⎰⎰⎰⎰⎰⋅-⋅⨯∇⨯-⨯∇⋅⨯∇=⋅-⨯∇⨯⋅∇-⨯∇⋅⨯∇=⋅-⨯∇⨯∇⋅eeeeeeeeV V i S m i m i V V i V m i m i V i m V i VV V V V VV d d d d d d d d J W S A W A W J W A W A W J W A W γγγγγ分析上式第二项，若有一部分处在外边界上，则由边界条件2s 上：()⇒=⨯⨯∇0n m e A γ()[]()[]0=⨯⨯∇⋅=⋅⨯∇⨯n m i n m i e A W e A W γγ被积函数为零。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第五章３-Ｄ静态磁场分析（标量法）5.1 在3-D静态磁场分析（标量法）中要用到的单元表1三维实体单元：单元维数形状或特性自由度SOLID5 3-D 六面体，8个节点每节点6个：位移、电势、磁标量位或温度SOLID96 3-D 六面体，8个节点磁标量位SOLID98 3-D 四面体，10个节点位移、电势、磁标量位、温度表2三维界面单元单元维数形状或特性自由度INTER115 3-D 四边形，4个节点磁标量位，磁矢量位表3三维连接单元单元维数形状或特性自由度3D杆状( Bar)、弧状(Arc)、线圈SOURC363个节点无(Coil)基元表4三维远场单元单元维数形状或特性自由度四边形，4个节点；INFIN47 3-D磁标量位、温度或三边形，3个节点INFIN111 3-D 六面体，8个或20个节点磁矢量位、磁标量位、电势、温度SOLID96和SOLID97是磁场分析专用单元，SOLID62、SOLID5和SOLID98更适合于耦合场求解。
5.2 磁标量位（MSP）法介绍在磁标量位方法中，可使用三种不同的分析方法：简化标势法(RSP)、差分标势法(DSP)和通用标势(GSP)法。
·若模型中不包含铁区，或有铁区但无电流源时，用RSP法。
若模型中既有铁区又有电流源时，就不能用这种方法。
·若不适用RSP法，就选择DSP法或GSP法。
DSP法适用于单连通铁区，GSP法适用于多连通铁区。
5.2.1单连通区与多连通区单连通铁区是指不能为电流源所产生的磁通量提供闭合回路的铁区，而多连通铁区则可以构成闭合回路。
参见图1(a)、(b)“连通域”。
数学上，通过安培定律来判断单连通区或是多连通区，即磁场强度沿闭合回路的积分等于包围的电流（或是电动势降MMF）。
因为铁的磁导率非常大，所以在单连通区域中的MMF降接近于零，几乎全部的MMF降都发生在空气隙中。
但在多连通区域中，无论铁的磁导率如何，所有的MMF降都发生在铁芯中。
5.3 3-D静态磁标势分析的步骤该分析类型与2－D静态分析的步骤基本一样：1.建立物理环境2.建模、给模型区域赋属性和分网格3.加边界条件和载荷（激励）4.用RSP、DSP或GSP方法求解5.观察结果5.3.1创建物理环境首先设置分析参数为“Magnetic-Nodal”，并给出分析题目。
然后用ANSYS前处理器定义物理环境包含的项目。
即单元类型、KEYOPT选项、材料特性等。
3D分析的大部分过程与2D 分析一致，本章下面部分介绍3D分析中要特殊注意的事项。
· SOLID96单元可为模型所有的内部区域建模，包括：饱和区、永磁区和空气区（自由空间）。
对于电流传导区，需用SOURC36单元来表示，关于电流传导区建模，后面有详细讲述。
·对于空气单元的外层区域，推荐使用INFIN47单元（4节点边界单元）或INFIN111单元（8节点或20节点边界单元）。
INFIN47单元和INFIN111单元可很好地描述磁场的远场衰减，通常比使用磁力线垂直或磁力线平行条件得到的结果更准确。
二种单元中，INFIN111更精确一些。
·缺省单位制使用MKS单位制(米-千克-秒国际单位制)，可用下列方式改变成其他单位制。
一旦选定，所有输入数据都应该使用该单位制。
为了方便建模，可以先在其他单位制系统下面建模（如毫米或英寸），然后进行缩放。
用下列方式定义单位制：命令：EMUNITGUI: Main Menu>Preprocessor>Material Props>Electromag Units ·根据用户设定的单位制，自由空间的相对导磁率将自动设定：在MKS单位制中，，或者根据用命令EMUNIT来设定一个值。
5.3.2设置GUI菜单过滤如果你是通过GUI路径来运行ANSYS，当ANSYS被激活后第一件要做的事情是选择菜单路径：Main Menu>Preferences，在对话框出现后，选择Magnetic-Nodal。
因为ANSYS会根据你选择的参数来对单元进行过滤，选择Magnetic-Nodal以确保能够使用用于3D静态磁场分析的单元。
5.3.3定义材料属性分析模型可有一种或多种材料区域：空气、导磁材料、导电区和永磁体。
每种类型的材料区具有所要求的材料性质。
ANSYS材料库自身带有几种磁性材料，可以直接把这些材料性质读入数据库，不用再手动的逐点输入。
如有必要，可对它们进行修改，以便与所分析的课题相匹配。
在ANSYS材料库中定义的磁性材料如下：材料材料性质文件Copper（铜）emag Copper. SI_MPLM3 steel（钢）emag M3. SI_MPLM54 steel（钢）emag M54. SI_MPLSA1010 steel（钢）emag Sa1010. SI_MPLCarpenter steel（硅钢）emag Silicon. SI_MPLIron Cobalt Vanadium steel（铁－钴－钒－钢）emag Vanad. SI_MPL该表中铜的材料性质定义有与温度有关的电阻率和相对导磁率，所有其他材料的性质均定义为B－H曲线。
对于列表中的材料，在ANSYS材料库内定义的都是典型性质，而且已外推到整个高饱和区。
你所需的实际材料值可能与ANSYS材料库提供值有所不同，因此，必要时可修正所用ANSYS材料库文件以满足用户所需。
5.3.3.1访问材料库文件：下面介绍读写材料库文件的基本过程。
详细参见《ANSYS入门指南》和《ANSYS基本过程手册》。
读材料库文件，进行以下操作：1. 如果你还没有定义好单位制，用/UNITS命令定义。
注意：缺省单位制为MKS,GUI列表只列出当前被激活单位制的材料库文件。
2. 定义材料库文件所在的路径。
（你需要知道系统管理员放置材料库文件的路径）命令：/MPLIB,read,pathdataGUI: Main Menu>Preprocessor>Material Props>Material Library>Library Path3. 将材料库文件读入到数据库中。
命令：MPREAD,filename,,,LIBGUI:Main Menu>Preprocessor>Material Props>Material Library>Import Library Main Menu>Preprocessor>Loads>-Load Step Opts-Other>Change MatProps>Material Library>Import Library写材料库文件，进行以下操作：1. 用MP命令或菜单Main Menu>Preprocessor>Material Props>Isotropic编辑材料性质定义，然后将改后的材料特性写回到材料库文件当中去。
2.在前处理器中执行下列命令：命令：MPWRITE,filename,,,LIB,MATGUI:Main Menu>Preprocessor>Material Props>Material Library>Export Library5.3.4定义材料属性和实常数的一般原则下面讲述关于设置物理模型区域的一般原则。
在“2-D谐波（AC）分析”一章中详细描述了2-D模型中需要设定的一些特殊区域。
5.3.4.1空气：说明相对磁导率为1.0。
命令：MP,murxGUI: Main Menu>Preprocessor>Material Props>Material Models > Electromagnetics > Relative Permeability >Constant5.3.4.2自由空间导磁材料区：说明B-H曲线，可以从库中读出，也可以输入自己定义的B-H曲线：命令：MPREAD,filename,…GUI:Main Menu>Preprocessor>Material Props>Material Library>Import Library 命令：TB,TBPTGUI:Main Menu>Preprocessor>Material Props>Material Models>Electromagnetics>BH Curve如果该材料是线性的，那么只需说明相对磁导率（对于均匀或者各向异性介质）。
如果定义B-H曲线，应满足下列准则，以保证准确地模拟材料属性：*输入B-H曲线必须要遵守的规则：1. B与H要一一对应，且应B随H是单调递增，如第2章中图1所示。
B-H曲线缺省通过原点，即0.0点不输入。
用下面的命令验证B-H曲线：命令：TBPLOTGUI: Main Menu>Preprocessor>Material Props>Material Models>Electromagnetics>BH Curve2. ANSYS计算的µ-H曲线应该是光滑连续的，可用TBPLOT命令画出µ-H曲线来进行验证（参见图2（b））。
B-H曲线应覆盖材料的全部工作范围。
如果需要超出B-H曲线的点，程序按固定斜率进行外推处理，固定斜率应等于或者略大于µr，你可以按照如下方式改变X-轴的范围，并用TBPLOT命令画图来观察其外推情况。
命令：/XRANGEGUI: Utility Menu>PlotCtrls>Style>Graphs对于永磁体和非线性各向异性材料，材料特性的输入与2D情况一样。
5.3.5建立模型《ANSYS建模和分网指南》对3-D标标量位方法进行磁场分析的建模有详细描述，但是对于电流源的处理有一些特殊考虑。
5.3.5.1建立电流传导区可以用基元模拟电流传导区域，不需要材料性质。
在3D标量位方法分析中，电流源不是有限元模型的一个组成的部分（在2D矢量位方法分析中是一个组成部分）。
只需用一个有限元哑元单元，SOURC36，来指明电流源的形状和位置。
可以在模型中的任意位置定义线圈、杆状、弧状电流源，电流源的大小和其他电流源数据可以通过哑元单元的实常数定义给出。
图3为用SOURC36原始元表示的一个电流源。
注意：即使采用半对称或四分之一对称模型，也要为整个电流源建模。
线圈和弧单元的内半径不能设置成为0。
因为SOURC36单元并不是一个真正的有限元，只能通过直接生成来定义它们，而不能通过实体建模的方式。