电磁场数值分析作业

合集下载

工程电磁场数值分析试题

工程电磁场数值分析试题

工程电磁场数值分析试题一、一同轴电缆,内导体(铜)外半径为0.01m,外导体(铜)内半径为0.03m,外导体厚度为0.003m,内外导体间有两层电介质,一层电介质为聚乙烯(13rε=),另外一层电介质为聚氯乙烯(26rε=),两层电介质厚度均为0.01m,内导体电位为5kV,外导体电位为0V。

(1)试用有限元法求内外导体间的电位和电场强度分布,(2)求此电缆中最大场强的位置和最大值,能否击穿电介质或发生局部击穿,(3)在不击穿的前提下,此电缆能承载的最大电压为多少?分析:参数设置:铜相对介电常数ε=1,电阻率ρ=1e-7Ω/m聚乙烯相对介电常数ε=3,电阻率ρ=1e+13Ω/m聚氯乙烯相对介电常数ε=6,电阻率ρ=1e+14Ω/m(1)其中电位分布及场强分布如下:通过定义路径(-0.033,0)到(0.033,0)分析其场强分布如下图:电压分布如下图:(2)从图中可以看出其场强最大值位于内导体外半径附近处取得距离内导体圆心0.0132m处取最大值422728V/m。

聚乙烯击穿场强为35-50MV/m,聚氯乙烯击穿场强为20-35MV/m,计算可知无法击穿电介质。

(3)聚乙烯击穿场强为35-50MV/m,聚氯乙烯击穿场强为20-35MV/m,按最小值计算理论上0.01m距离上其耐压分别为350kv和200kv,所以电介质不会被击穿。

如不击穿理论上应能够承压200kv。

二、一同轴电缆,内导体(铜)外半径为0.01m,外导体(铜)内半径为0.03m,外导体厚度为0.003m,内外导体间有一层电介质,电介质为聚氯乙烯(26ε=),电介质厚度均为0.02m,内导体电位为5kV,外导体电位为0V。

(1)试用有限元法求内外导体间的电位和电场强度分布,(2)求此电缆中最大场强的位置和最大值,能否击穿电介质或发生局部击穿,(3)在不击穿的前提下,此电缆能承载的最大电压为多少?(4)通过一题和二题的对比,说明同轴电缆的内外导体间用一层还是二层电介质比较好,为什么?分析:参数设置:铜相对介电常数ε=1,电阻率ρ=1e-7Ω/m聚氯乙烯相对介电常数ε=6,电阻率ρ=1e+14Ω/m (1)其中电位分布及场强分布如下:通过定义路径(-0.033,0)到(0.033,0)分析其场强分布如下图:电压分布如下图:。

西安电子科技大学 研究生 电磁场数值分析期末考试题

西安电子科技大学 研究生 电磁场数值分析期末考试题

西安电子科技大学何超电磁场数值分析考点1:矩量法的一般过程(算子方程、离散化过程、选配过程、矩阵方程求解)。

给定算子方程和基函数,采用伽略金法,计算阻抗矩阵和激励电压矩阵,从而求得电流系数矩阵,即得到方程的近似解。

(矩阵维数一般为2×2,或3×3,便于计算)。

1/link?url=oRwkn_6gajdEKC3YUFvvipOKLuZJXnVk43odUwyDWYRao nT1SlZLKEq9PCQba5xPYg_7mXpK8pZW0R-_RfT5EOXLvj0BKqKmQ6cfXMuW8P7有3个矩量法例题考点2:ScaLAPACK 的矩阵分布方式。

给定进程网格,矩阵分块大小,要求能写出按ScaLAPACK矩阵分布方式,每个进程对应的矩阵元素。

?1 并行矩阵填充在PC集群系统中MPI并行矩量法研究36 37考点3:temporary block column 对active block column 分解产生的影响.对于当前活动列块(即正在进行LU分解的列块),要能够分析其左侧临时列块对其LU分解所产生的影响。

?英文书写得很详细了啊45--55有lu分解将系数矩阵A转变成等价两个矩阵L和U的乘积,其中L和U分别是下三角和上三角矩阵。

当A 的所有顺序主子式都不为0时,矩阵A可以分解为A=LU,且当L的对角元全为1时分解唯一。

其中L是下三角矩阵,U是上三角矩阵。

4阶矩阵的LU分解[1]高斯消元法见数值分析教材考点4:积分方程的建立要求掌握EFIE 、MFIF 、PMCHW(电场、磁场、表面积分方程)根据等效原理建立的过程,即对于给定的问题(PEC (理想导体)或介质)能根据等效原理建立积分方程(不要求写出场的位函数表达式,主要考察方程建立的思想)。

看矩量法的书那个英文书只有EFIE 等效原理EFIE考点 5:RWG 基函数考察 RWG 基函数的 表达式,以及其 特点,对于给定的一个三角形网格图要能够标出哪些地方( 公共边上) 存在基函数。

《电磁场问题数值算法》经典习题解答

《电磁场问题数值算法》经典习题解答
《电磁场问题数值算法》学习经典习题
吕波
矩量法
问题:求解金属圆柱在水平入射均匀平面波照射下,面电流分布及由此产生的远 区散射场分布。 分析:可将面电流密度用脉冲函数展开,即:
N J sz i i (r r i ) ……………………………………………………(1-1) i 1
它表示把导体横截面的周界分成 N 个小段, 且设每一个小段中的面电流密度 i 作 均匀分布。这样算子方程为:
jk ( x j cos i y j sin i )
(j=1,2,…,N)……………………(1-8)
Matlab 源程序如下:
%矩量法计算圆柱散射 %copyright by lvbo % clear; yita=120*pi%自由空间波阻抗 freq=5*10^8%频率为 500MHz k=(2*pi*freq)/(3*10^8)%该频率的传播常数 N=100%将圆柱导体横截面的周界 l 分成 N 个小段 r=0.1%圆柱导体的半径为 0.1 detaL=2*pi*r/N%每个小段的长度 2003,6,24
Hx1=zeros(N+1,N); Hy1=zeros(N,N+1); %迭代后新时刻场量值 Ez2=zeros(N+1,N+1); Hx2=zeros(N+1,N); Hy2=zeros(N,N+1); %设置脉冲源 tao=2*10^-10;%高斯脉冲的宽度 t0=tao/2%高斯脉冲峰值出现的时刻 deltaL=L/N;%空间间隔要求 deltaT=tao/(4*N);%时间间隔要求 t=0;%可以修改 t 的值一观察不同时刻的电磁场 %迭代过程 for dc=1:1000%迭代的次数 ca=1;%迭代系数 cb=deltaT/yps; cp=1; cq=deltaT/u; %激励源加入 It=exp(-4*pi*(t-t0).^2/(tao^2)); %非边界场值迭代 clear i j; for i=1:(N+1) for j=1:N Hx2(i,j)=cp*Hx1(i,j)-cq*(Ez1(i,j+1)-Ez1(i,j))/deltaL; end end clear i j; for i=1:N for j=1:(N+1) Hy2(i,j)=cp*Hy1(i,j)+cq*(Ez1(i+1,j)-Ez1(i,j))/deltaL; end end clear i j; for i=2:N for j=2:N

数值分析方法在电磁场计算中的应用

数值分析方法在电磁场计算中的应用

数值分析方法在电磁场计算中的应用电磁场是物理学中最重要的一部分之一,它广泛应用于现代工业、交通、通信、能源和医疗设备等领域。

因此,研究电磁场的行为对于建立新技术和改进现有技术非常重要。

不过由于电磁场是一个非线性的动态系统,因此分析它的行为非常困难。

为了解决这个问题,我们需要数值分析方法来帮助我们更好地理解电磁场的行为。

电磁场的计算方法有很多种,常见的有有限元法、有限差分法等等。

本文将着重介绍有限差分法在电磁场计算中的应用。

有限差分法是经典的数值计算方法,它是一种数值求解偏微分方程的方法。

它的基本原理是将要求解的偏微分方程转化为差分方程,然后利用计算机来求解这个差分方程。

有限差分法的求解过程是离散化的,因此它更便于计算机的处理,同时它的数值误差也比较小。

有限差分法在电磁场计算中的应用非常广泛。

我们可以利用有限差分法来计算电磁场的强度、分布、辐射等参数。

下面我们将介绍一些在电磁场计算中使用有限差分法的实例。

首先,我们来看一个简单的电磁场问题:平面电容器之间的电场强度。

在这个问题中,我们需要求解电场的分布情况。

我们可以利用有限差分法来求解这个问题。

将计算区域离散化成若干个网格点,然后利用电场的高斯定理,将它的积分式子转化为差分式子,最后用差分方程来求解电场值。

在电磁场计算中,还有一些需要注意的问题。

首先是边界条件的处理。

由于有限差分法是一种离散的方法,因此我们需要在计算区域的外部放置边界条件。

这些边界条件包括电场的值、电势的值、电荷密度等等。

其次是计算精度的问题。

由于有限差分法是一种数值方法,因此它的计算精度有时会受到误差的影响。

我们可以通过适当地选择网格点的数量和大小来提高计算精度。

总体来说,有限差分法在电磁场计算中的应用非常广泛,并且具有很好的计算效果。

在实际应用中,我们需要根据具体问题选择合适的数值计算方法,并且在计算时注意处理边界条件和计算精度的问题。

工程电磁场数值分析(有限差分法)

工程电磁场数值分析(有限差分法)

进一步的参考书:
胡之光. 电机电磁场的分析与计算. 北京:机械工业出版 社,1989
从有限差分法看数值解的基本思想
离散解(数值解)的概念
方程的离散——化无限维问题为有限维问题 (化微分方程为代数方程组,借助计算机求解)
解的离散——离散点上的数值解
数值法的一般步骤
求解区域的离散(前处理) 方程的离散 代数方程组的求解 离散数据的分析(后处理)
h
2
F0
对所有的节点都建立一个方程,N 个 齐次第二类边界条件 节点有 N 个未知数,建立 N 个方程。
③ 求解线性代数方程组
N个方程联立成为线性代数方程组
[ A][ ] [ F ]
求解得到节点上的电位值。
④ 后处理
计算电场强度E,
Ex 0
1 3
2h
Ey
0
2 4
2h
绘制场图;计算电容、受力、能量分布等。
一些实际需要考虑的问题
精度问题与计算量
媒质交界面条件问题
边界条件处理问题
不规则网格问题
三维问题
有限差分法是最古老、最直观的一种数值方法,直至现
在仍有强大的生命力,在许多学科领域广为应用。在电磁场 领域,目前最受关注的是时域有限差分法(Finite Difference Time-Domain Method, FDTD)和有限体积法(Finite Volume Method, FVM)。
各种数值方法的不同之处
在于离散方程所依据的原
理不同,从而导致方程求
解技术、求解效率、适用 对象等的不
有限差分法(Finite Differential Method,FDM) 是基于差分原理的一种数值计算法。其基本思想是: 将场域离散为许多小网格,用差分代替微分,用差商 代替求导,将求解连续函数 的泊松方程的问题转换 为求解网格节点上 的差分方程组的问题。

数值分析求解电磁场问题

数值分析求解电磁场问题

《数值分析》课程设计课程设计报告一、实验目的通过设计,使得学生熟悉和掌握所学习的数值分析理论知识和运算方法,并掌握相关的程序设计方法,培养学生基于程序设计运用理论能力。

1. 巩固和掌握数值分析理论课程所学习的理论知识。

2. 学习和掌握基于MATLAB解决问题的基本方法和程序设计方法。

二、实验原理MATLAB 作为一款优秀的程序仿真开发工具,集科学计算、图像处理和声音处理于一体,具有许多突出优点:(1)功能强大。

主要包括数值计算和符号计算、计算结果和编程可视化、数字和文字统一处理、离线和在线计算,针对不同领域提供了丰富的工具箱。

还具有可扩展的特征,用户可根据自己的需要任意扩充函数工具库。

(2)界面友好,编程简洁、效率高。

MATLAB 是一种以矩阵为基本单元的可视化程序设计语言,语法结构简单,数据类型单一,指令表达和标准教材的数学表达式相近。

(3)强大的绘图功能。

具有强大的二维、三维绘图功能,可以直观的从图形衡量程序的效果,这是任何其它语言无法比拟的。

(4)简单好学。

特别适合于初学者和需要快速创建应用程序的用户,用户可以在短时间内掌握其主要内容和基本操作。

正是由于其强大的功能和广泛的适用性,MATLAB 才得到了用户的普遍认可,在自动控制、信号处理、神经网络等诸多方面,都有广泛的应用。

电磁场数据获取系统中,往往涉及大量的数值计算,这些数据可从电磁场数值计算模拟和实验测试获得。

但是手工计算结果往往会很繁琐、甚至出现错误,徒手绘制图形,则更加复杂。

然而,利用MATLAB复数运算功能则可以方便的将计算结果显示出来。

另外,采用MATLAB 语言,能够直观地模拟和演示各种电磁现象,给出直观的图像,更加形象、精确、可视性很强。

三、数学模型静电场场量的计算和场图的绘制是求解电磁场问题的基础,在涉及到矢量的积分时求解极其复杂,但如果应用MATLAB 进行辅助求解绘图时则很容易方便。

本设计建立数学模型如下:已知空间电位分布为:V (x,y,z)=log(x2+y2),计算空间各点电场,并画出电位线和电力线。

电磁场数值分析期末

电磁场数值分析期末

《电磁场数值分析》(期末作业)--- 2019学年 ---学院:学号:姓名:联系方式:任课教师:2019年5月作业1模拟真空中二维TM 电磁波的传播,边界设置为一阶Mur 吸收边界,观察电磁波的传播过程。

波源为正弦函数:sin()sin(2)25z t cE t n t ωπ==∆代码: clc clear close allxmesh =150; ymesh =150;mu0=4*pi*1.0E-7; eps0=8.85E-12;C= 3.0E8; dx=1.0; dt=0.7*dx/C; timestep=150; ez( 1:xmesh+1,1:ymesh+1 ) = 0.0; hx( 1:xmesh+1,1:ymesh ) = 0.0; hy( 1:xmesh,1:ymesh+1 ) = 0.0;coef1 = dt/( mu0 * dx ); coef2 =dt/( eps0 * dx );coef3=(C*dt-dx)/(C*dt+dx); ez1=ez;for now = 1 : timestephx = hx - coef1 * ( ez( :, 2 : ymesh+1 ) - ez( :, 1 : ymesh ) ); hy = hy + coef1 * ( ez(2 : xmesh+1, : ) - ez(1 : xmesh, : )); ez( 2 : xmesh , 2 : ymesh ) = ez( 2 : xmesh , 2 : ymesh ) - ... coef2 * ( hx( 2 : xmesh, 2 : ymesh ) - hx( 2 : xmesh , 1 :ymesh - 1) ) + ...coef2 * ( hy( 2 : xmesh ,2 : ymesh ) - hy( 1 : xmesh - 1,2 : ymesh) );ez(1,:)=ez1(2,:)+coef3*(ez(2,:)-ez1(1,:));ez(xmesh+1,:)=ez1(xmesh,:)+coef3*(ez(xmesh,:)-ez1(xmesh+1 ,:));ez(:,1)=ez1(:,2)+coef3*(ez(:,2)-ez1(:,1));ez(:,ymesh+1)=ez1(:,ymesh)+coef3*(ez(:,ymesh)-ez1(:,ymesh +1));ez( xmesh/2+1, ymesh/2+1) = sin( now * dt * 2 * pi * C / 25.0 ); mesh(ez);pause(0.05)ez1=ez;end结果与分析:第10时间步第100时间步第150时间步作业2基于Pocklington方程用MoM分析半波对称振子天线:观察天线线径和分段数目分别取不同值对天线阻抗和辐射特性的影响(半径分别取0.001λ,0.0001λ,0.00001λ,分段数取11,21,31,可列表说明)代码:clear all; close all; clc;% 初始化参数c=3e8; % 光速r=1 % 波长f=c/r; % 频率w=2*pi*f; % 角频率e0=8.85e-12; % 介电常数u0=4*pi*1e-7; % 磁导率a=0.00001*r; % 半径L=0.5*r; % 振子长度k=2*pi/r; % 波数N=11; % 分段数(奇数段)dl=L/(N+1); % 每段长度(分母中+1 为两头半段之和)l=L/2-dl/2; % 两头空出半段,满足电流为0的边界条件lz=-l:dl:l;lzs=lz(1:N); % 每一小段的起点坐标lzm=lz(1:N)+dl/2; % 每一小段的中点坐标lze=lz(2:N+1); % 每一小段的终点坐标%阻抗矩阵元素求解fi=log(dl/a)/(2*pi*dl)-k/(4*pi)*1i;fi_1=exp(-k*dl*1i)/(4*pi*dl);fi_2=exp(-k*2*dl*1i)/(8*pi*dl);z=ones(N,N);for m=1:Nfor n=1:Nif m==nfi1=fi;fi2=fi_1;fi3=fi_1;z(m,n)=((k^2*dl^2-2)*fi1+fi2+fi3);elseif abs(m-n)==1fi1=fi_1;fi2=fi;fi3=fi_2;z(m,n)=((k^2*dl^2-2)*fi1+fi2+fi3);elsefi1=exp(-k*abs(m-n)*dl*1i)/(4*pi*abs(m-n)*dl);fi2=exp(-k*abs(m+1-n)*dl*1i)/(4*pi*abs(m+1-n)*dl); fi3=exp(-k*abs(n+1-m)*dl*1i)/(4*pi*abs(n+1-m)*dl); z(m,n)=((k^2*dl^2-2)*fi1+fi2+fi3);endendend%电压矩阵求解V=zeros(N,1);V((N+1)/2)=-1*(1i*w*e0);% 计算电流系数矩阵I=z\V;% 计算输入阻抗Z_in=1/I((N+1)/2);disp(['输入阻抗 = ',num2str(Z_in)]);% 计算振子上归一化电流分布I_amp=abs(I); Max=max(I_amp);Iunit2=[0;I_amp/Max(1);0]; % 两端零电流figure(1);h=0:dl/r:L/r;Ithe=sin(pi*h*r/L); % 半波振子电流解析值plot(h,Iunit2,'b',h,Ithe,'r','linewidth',2);legend('pocklinton','解析值');grid on;xlabel('电长度L/\lambda');ylabel('归一化电流');% 画方向图theta=0:0.01:2*pi;abs_f=zeros(1,length(theta));for n=1:1:Nabs_f=abs_f+I(n)*exp(k*(n*dl-L/2)*cos(theta)*1i);endabs_f=abs(sin(theta)*dl.*abs_f);Max_f=abs(sum(I)*dl);Far_patten2=abs_f/Max_f(1);theta_2=0:0.1:2*pi;Far_theory=abs((cos(k*(L/2)*cos(theta_2))-cos(k*L/2))./si n(theta_2));figure(2);polar(theta,Far_patten2,'-b');hold on;polar(theta_2,Far_theory,'or');hold off;legend('pocklinton','解析值');title('半波振子天线E面方向图');figure(3);polar(theta,ones(1,length(theta)),'-b');title('半波振子天线H面方向图');% 半波振子增益I_in=I((N+1)/2);A=(w*u0)^2/(4*pi*sqrt(u0/e0)*real(Z_in)*(abs(I_in))^2); G_theta=A*abs_f.^2;Max_gain=max(G_theta)Max_gain_dB=10*log10(Max_gain);disp(['半波振子增益 = ',sprintf('%.4fdBi', Max_gain_dB)]); 结果与分析:作业3基于电场积分方程用MoM分析对称振子天线:计算振子总长度分别为0.25λ ,0.5λ,λ,1.5λ时,振子的输入阻抗和E面方向图。

500kV变电站工频电磁场数值分析

500kV变电站工频电磁场数值分析
错 。该变 电 站有 两 回母 线 四 回出线 , 线 离 地 高 度 母
1 变 电站 电场 强 度 分 析 原 理
交流高压变电站内的带电导体 比线路复杂, 地
1m, 间距离 65 出线 的最低点离 地高度是 2 m。 8 相 .m, 3 为 了与文 献 . 的实 测 结 果 对 比 , 文 选 取 4 中 本 了该 变 电站 内 50 V 开关 场 I组母 线 下 方 作 为 研 0k
利用前述方法计算 了该 区域 的工 频磁场 , 计算结
果 如 图 8所 示 。
从 图 8可 以看 出 , 考 虑 设备 对 磁 场 强度 的影 不
响时 , 由于载 流导 线 的离 地距 离 比较 高 , 电流 较 小 ,
度 的分 布规 律相 似 高场 强点 主要分 布在 通流 量最
大 , 离地 面最 近 的导线 附近 , 些地 方如 变压 器和 距 有
L分为 两段 , 别 对 应 两个 角 分
原 理得 到此 处 的磁 感应 强度 为 :
数 十年来 , 国持 续对 10~70 V不 同 电压等 我 1 5k
和 :, 用 叠 加 应
级 的变电站 的工 频电场做 了大量 的现场 测量 工作 , 测 得 了不 同电压 等级 的变 电站 电场分布实例 。 J 测量的 1 k 1 V变 电站 中 电场 强度 的最大值 为 1. 0 6 6k/ 20V变 电站中电场强度的最大值为 1.9V 1V m, k 2 33k/
该站中大部分 区域的磁场强度都 比较小。连接 C 、 T P T的管母 线 由于离 地高度 较小 , 以其 下方 的磁 场 所
电抗器 , 因为有漏 磁 现象也 可能 产生较 大 的场强 。
强度最高 , l T左右 , 在 5 远远低于 INT C IP导则规

电磁场数值计算上机题报告

电磁场数值计算上机题报告

电磁场数值计算上机题报告第一题计算长直接地金属槽中的电场分布。

金属槽横截面如图1所示,其侧壁与底面电位均为零,顶盖电位相对值为10。

槽内电位函数满足拉普拉斯方程。

计算槽内电位分布。

要求:(1)先用正方形网格粗分,每边取4个网格计算,取不同的松弛因子,比较其收敛速度。

取计算精度为千分之一。

(2)划分网格加倍,计算电位分布,并与上面计算结果比较。

(一)建立问题的数学物理模型首先列出方程及其边界条件 槽内的电位满足二维的拉普拉斯方程:222220x y ϕϕϕ∂∂∇=+=∂∂其中,x y 的范围是:0,0x a y a <<<<边界条件是:0000;10x x ay y aϕϕϕϕ========图1 (二)算法设计及其实现在本题中,因为区域为正方形区域,网格采用正方形网格,每边四个网格,因此,每边要有5个结点, 网格数m=n=4,比较少,不能用课本中的公式计算收敛因子,取收敛因子为 1.4α=,计算的程序的流程为:①选取计算的场域,并划分网格,网格划分如右图2所示:右图表示网格的划分,共16个网格,一共划分了25个结点,每个节点用相应的下标(i,j )来表示,对应的电位为(,)u i j 。

根据题意,边界条件的 处理如下: 图 2(1,)(,1)(,1)0(,)10(21)u j u i u n u k j k n ===⎧⎨=≤≤-⎩其中n 为一行对应的节点数,根据题意,这里n=5② 用(,)u i j 表示节点的电位,设经过第n 次迭代之后的结点电位用(),n i j u 来表示,则超松弛迭代法的差分格式(记,(,)i j u u i j =)为:(1)()()()(1)(1)(),,1,,11,,1,(4)4n n n n n n n i j i j i j i j i j i j i j u u u u u u u α+++++--=++++- 边界条件前面已经给出 ③给各个节点赋初值,对于非边界上的点(对于边界点的值前面已经赋过初值),如下10(,)(1)1u i j j n =--④迭代计算,直到已经满足精度条件为止,这里精度为0.0001,最后输出计算的结果,输出结果见生成的数据文件shuju.txt 中 ⑤计算框图如图3所示图 3 计算框图⑥用fortran90语言编写计算机程序,计算各点的电位,程序清单如下:(三)计算结果及数据分析当网格数目为4(节点数目为5),加速收敛因子为1.4,计算精度为0.0001时,计算出各结点电位如下所示(对不同的收敛因子的比较后面进行)不可能完全对称,只能是近似对称的),这点不难从理论上进行分析得到。

工程电磁场数值计算

工程电磁场数值计算

工程电磁场数值计算大作业报告一、大作业要求运用FEM法求解算题5—8,删去要求(2),设其具有平行平面磁场分布的特征。

作业题目如下所示:二、问题分析及建立模型根据P149对平行平面场的静电场和磁场统一的数学模型的描述我们可以得到此问题对应的偏微分方程及相应的定解问题为:322220000300;;0;ρρμρϕ===⎧∂∂+=⎪∂∂⎪⎪==⎨⎪∂⎪=⎪∂⎩-y x H A A s y A A Ain x n进而可以求得此题对应的泛函及等价的变分问题为:2422221()221min(0;0)2S l l S A A A F A JA dxdy dl x y n A A A dxdy J x y n μ+⎡⎤⎛⎫∂∂∂⎛⎫=+--⎢⎥ ⎪ ⎪∂∂∂⎝⎭⎢⎥⎝⎭⎣⎦⎡⎤⎛⎫∂∂∂⎛⎫=+===⎢⎥ ⎪ ⎪∂∂∂⎝⎭⎢⎥⎝⎭⎣⎦⎰⎰⎰⎰⎰00;==y A 3003;ρρμρϕ==-H sin A根据以上条件,我们可以把此题与例5-2作比较,他们的边界条件形式已经基本一致了,所以我们可以利用EMF2D的程序对此题进行计算。

下面所以下我们的主要解题思路。

1、由于是一个圆形区域,且是对称的,所以我们只需求1/4圆周即可。

我们运用圆域剖分程序CAMG对整个区域进行剖分。

这里我们需要注意的是最外层的边界条件,我们选用选定10倍半径,即1米,进行三段剖分。

2、运用程序EMF2D,把圆域剖分出来的结果当作此程序的输入。

需要注意的是需要对剖分出来的最外层的点,进行“手动输入”。

我们需要注意两个程序的输入输出的格式进行统一,修改EMF2D 的强制边界条件程序FB。

三、程序及结果1、圆域剖分我们并没有改变什么CAMG程序,程序如下我们的输入数据如下:由输入可以知道我们内环分7段,中环分8段,外环分6段。

得到的输出结果CAMGOUT结果如下:前面表示节点坐标,后面表示每个三角元的顶点编号。

根据结果,我们得知了内环剖分了1~49个节点,中环剖分了49~169个节点,外环剖分了169~190的节点。

工程电磁场数值分析解读

工程电磁场数值分析解读

工程电磁场数值分析解读工程电磁场数值分析是一种应用有限元法来计算和解决电磁场问题的方法。

该方法通过将电磁场的连续性方程离散化为有限个小单元,再通过求解矩阵方程组来获取数值解。

这种分析方法能够定量计算电磁场的分布和特性,并为工程设计和优化提供重要的参考依据。

对于电磁场数值分析的解读,可以从以下几个方面进行讨论:首先,可以对电磁场的分布进行解读。

通过数值计算,可以得到电磁场在不同位置的数值结果,可以用来表示电磁场的强弱、方向和空间分布特性。

可以对电磁场的分布情况进行比较和分析,以评估电磁场的均匀性和一致性,为设计提供优化方案。

其次,可以对电磁场的特性进行解读。

通过数值计算,可以计算并分析电磁场的一些重要参数,如电场强度、磁场强度、电势、电感、电容等。

这些参数能够揭示电磁场的基本特性,并对电磁设备和系统的工作性能进行评估和优化。

另外,可以对电磁场的影响进行解读。

电磁场数值分析能够计算出电磁场对物体的作用效果,如力、热、电磁感应等。

通过对电磁场的影响进行解读,可以预测电磁场对设备、器件和系统的影响,并为电磁兼容性设计提供技术支持。

此外,还可以对电磁场数值分析方法和结果的准确性进行解读。

电磁场数值分析是一种近似求解的方法,所得数值结果可能与实际情况存在一定差异。

因此,在解读时需要对数值结果进行验证和确认,通过模型实验或其他可靠手段来验证模型的准确性和可靠性。

总之,工程电磁场数值分析是一种重要的工程设计方法,能够定量计算和解决电磁场问题。

通过对电磁场分布、特性和影响等方面进行解读,可以为工程设计和优化提供重要的参考依据。

同时也需要关注分析方法的准确性和结果的可靠性,以确保分析结果的准确性。

电磁场数值分析与仿真-无刷直流电机空载磁场分析

电磁场数值分析与仿真-无刷直流电机空载磁场分析

电磁场数值分析与计算仿真大作业题目:无刷直流电机空载瞬态磁场分析姓名:学号:完成时间: 2019年11月01日目录1问题描述 (1)1.1仿真模型参数 (1)2仿真软件简介 (1)3模型建立 (2)3.1 创建项目 (2)3.2 绘制电机几何模型 (3)3.3 绘制电机定子槽几何模型 (4)3.4 创建电机定子冲片模型 (5)3.5 创建导线模型 (6)3.6 创建永磁体模型 (7)3.7 创建转子轭模型 (8)3.8 旋转部分设定 (9)4材料属性 (11)4.1 材料指定 (11)4.2 新材料添加 (14)5网格剖分 (16)6激励加载 (17)7边界条件处理 (22)8求解器选择 (24)9仿真结果 (27)10仿真结果讨论 (38)1问题描述无刷直流电机由电机主体和驱动器构成,是一种典型的机电一体化产品。

无刷直流电机是采用半导体开关器件来实现电子换向的,即用电子开关器件代替传统的接触式换向器和电刷。

它具有可靠性高、无换向火花机械噪声低等优点,广泛应用于高档录音座、录像机、电子仪器及自动化办公设备中。

本例使用Maxwell 16.0版本的电磁仿真软件来对无刷电机进行空载瞬态磁场分析。

1.1仿真模型参数表1 两相无刷直流电机参数表模型参数额定功率0.55kw额定直流工作电压220v极数 4定子槽数24额定转速3000rpm轴向长度100mm定子外径110mm定子内径70mm定子槽口宽度2mm定子槽底宽度15mm定子槽深35mm绕组导线直径8mm转子轭半径55mm2仿真软件简介Maxwell 16.0是一款由Ansys公司推出的电磁场仿真设计软件,主要适用于各类机电和电气设备的仿真操作,软件拥有直观化的用户界面和强大、高效的仿真功能,可以将各类复杂的设计和分析变得更加简单化,Maxwell 16.0可以广泛地应用于汽车、致动器、变压器、传感器和线圈等设计领域。

图1.1 Maxwell 16.03模型建立3.1 创建项目启动Maxwell 16.0并建立新的项目文件。

工程电磁场数值分析(有限元法)

工程电磁场数值分析(有限元法)
使用适当的数值方法求解离散方程组,得到场函数的近似解 。
04
有限元法在工程电磁场中的应用
静电场问题
总结词
有限元法在静电场问题中应用广泛,能够准确模拟和预测静电场 的分布和特性。
详细描述
静电场问题是指电荷在静止状态下产生的电场,有限元法通过将 连续的静电场离散化为有限个单元,对每个单元进行数学建模和 求解,能够得到精确的解。这种方法在电力设备设计、电磁兼容 性分析等领域具有重要应用。
单元分析
对每个单元进行数学建模,包 括建立单元的平衡方程、边界 条件和连接条件等。
整体分析
将所有单元的平衡方程和连接 条件组合起来,形成整体的代 数方程组。
求解代数方程组
通过求解代数方程组得到离散 点的场量值。
有限元法的优势和局限性
02
01
03
优势 可以处理复杂的几何形状和边界条件。 可以处理非线性问题和时变问题。
传统解析方法难以解决复杂电磁场问题,需要采用数值分析方法 进行求解。
有限元法的概述
有限元法是一种基于离散化的数值分 析方法,它将连续的求解域离散为有 限个小的单元,通过求解这些单元的 近似解来逼近原问题的解。
有限元法具有适应性强、精度高、计 算量小等优点,广泛应用于工程电磁 场问题的数值分析。
02
静磁场问题
总结词
有限元法在静磁场问题中同样适用,能够有效地解决磁场分布、磁力线走向等问题。
详细描述
静磁场问题是指恒定磁场,不随时间变化的磁场问题。有限元法通过将磁场离散化为有限个磁偶极子,对每个磁 偶极子进行数学建模和求解,能够得到静磁场的分布和特性。这种方法在电机设计、磁力泵设计等领域具有重要 应用。
有限元法的基本步骤
01

电磁场数值分析(西电)

电磁场数值分析(西电)
a heterogeneous space
volumetric equivalence principle
However, the equivalent sources are unknowns to be determined.
The derivatives in (1.19)and (1.20) must be interpreted in the context of generalized functions.
2
Maxwell’s equations
E j0r H H j0 r E
(1.1) (1.2) (1.3) (1.4)
0r H 0
0 r E 0




source-free region with relative permittivity r and permeability r Have specialized that the medium is linear and isotropic and j t the electromagnetic field has time dependence e
(1.17) (1.18) (1.19) (1.20)
10
Volumetric equivalence principle for penetrable scatterers (III)
Equations (1.13)-(1.16)
a homogeneous space
Equations (1.1)-(1.4)
7
Two-dimensional problems and the scalar Helmholtz equations (II)

电磁上数值分析作业

电磁上数值分析作业

有限元法解平行板电容器电位欧阳文辉(云南师范大学理 论物理)摘要 有限元法是以变分原理和部分插值为基础的方法,即将定解区域划分为许多小单元,然后按单元分别插值合并起来得到总插值,在以求泛函极值的方法来得出我们所需的近似解答。

这种方法不仅保有原来变分方法的优点,而且兼有差分方法的灵活性。

关键词 三角元 泛函 电位1引言 很多物理上的问题,虽然数学上已能用一些微分方程来表示。

但是有时直接求解这些微分方程式困难的,例如求解电磁波在任意截面波导中传播的问题就是这样,这时就需要采用近似的方法。

有一种方法就是将所研究的微分方程边值问题转化成相应的变分问题,也就是找出与该微分方程等价的泛函,然后用求极值的方法来定出该微分方程的解答。

可以发现用这种方法可以带来很大的方便。

变分方法就是研究求泛函极大值或极小值的方法。

在许多实际问题中,静电场是由带电导体决定的,例如电容器内部的电场是由作为电极的两个导体板上所带电荷决定;又如电子光学系统的静电透镜内部,电场是由分布在电极上的自由电荷决定的。

这些问题的特点是自由电荷只分布在导体的表面,在空间没有其他自由电荷分布。

因此我们选择这些表面作为区域V 的边界,则在V 内部自由电荷密度ρ=0,因而泊松方程化为比较简单的拉普拉斯方程20ϕ∇=,产生这电场的电荷分布在区域边界V 上,他们的作用通过边界条件反应出来。

2拉普拉斯方程和其等价的泛函接下来我们讨论电极结构如图1所示,上极板的电位为50V ,下极板接地点位为0V 则在区域S 内,电位满足拉普拉斯方程222220u uu x y∂∂∇=+=∂∂,(1) 它服从的边界条件为y=0,u=50V;y=2,u=50v,它的等价泛函为22()()s u u J dxdy x y ⎡⎤∂∂=+⎢⎥∂∂⎣⎦⎰ (2)设划分的三角形有k 个,则有2211()()kketet t t u u J J dxdy x y ==⎡⎤∂∂==+⎢⎥∂∂⎣⎦∑∑⎰, (3)来研究如图2所示的et 元,将它的三个点顶按逆时针方向编号为i 、j 、m 并设三个顶点的坐标分别为(x i ,y i )、(x j ,y j )、(x m ,y m )。

电磁场数值计算上机作业

电磁场数值计算上机作业

电磁场数值计算上机作业报告一、 有限差分法及原理有限差分法基本思想是把连续的定解区域用有限个离散点构成的网格来代替,这些离散点称作网格的节点;把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似;把原方程和定解条件中的微商用差商来近似,积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组 ,解此方程组就可以得到原问题在离散点上的近似解.然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解.在采用数值计算方法求解偏微分方程时,若将每一处导数由有限差分近似公式替代,从而把求解偏微分方程的问题转换成求解代数方程的问题。

求解拉普拉斯方程:22220x y ϕϕ∂∂+=∂∂为简单起见,将场域分成足够小的正方形网格,网格线之间的距离为h ,0h →。

节点0、1、2、3、4上的电位分别用ϕ、1ϕ、2ϕ、3ϕ和4ϕ表示。

点1、点3在x 0处可微,沿x 方向在x 0处的泰勒级数展开式为2323100002311()()()()().2!3!h h h x x x ϕϕϕϕϕ∂∂∂=-+-+-+∂∂∂2323300002311()()()()().2!3!h h h x x xϕϕϕϕϕ∂∂∂=-+-+-+∂∂∂点2、点4在y0处可微,沿y 方向在y0处的泰勒级数展开式为2323200002311()()().2!3!h h h y y yϕϕϕϕϕ∂∂∂=++++∂∂∂ 2323400002311()()()()().2!3!h h h y y yϕϕϕϕϕ∂∂∂=-+-+-+∂∂∂忽略高次项22212340000224()()()4h x y ϕϕϕϕϕϕϕϕ⎡⎤∂∂+++=++=⎢⎥∂∂⎣⎦稍作变化得到拉普拉斯方程的五点差分格式:123404ϕϕϕϕϕ+++=利用超松弛迭代法求解以上差分方程,二维场拉普拉斯方程等距剖分差分格式公式为:其中为超松驰因子,. 计算流程如图1所示:图1 超松驰迭代法计算流程利用上述方法求解两题。

电磁场数值分析(西电)

电磁场数值分析(西电)

Electromagnetic theory
Functional and matrix
Method of moment (MoM) (矩量法)
Finite element method (FEM) (有限元法)
Finite-difference time-domain methods (时域有限差分方法)
(1.12)
8
Volumetric equivalence principle for penetrable scatterers (I)
Differential equations Integral equations
To simplify the formulation of integral equations
a heterogeneous space
volumetric equivalence principle
However, the equivalent sources are unknowns to be determined.
The derivatives in (1.19)and (1.20) must be interpreted in the context of generalized functions.
1 2 E k rE 0 r 2 1 H k rH 0 r
(1.9)
(1.10)
where k 2 2 0 0 ,the parameter k is known as the wavenumber of the medium.
Convert the original problem into an equivalent problem

电磁场数值分析期中

电磁场数值分析期中

《电磁场数值分析》(期中作业)--- 2019学年 ---学院:学号:姓名:联系方式:任课教师:2019年5月作业1一个二维正方形(边长a=10mm)的静电场区域,电位边界条件如图所示(单位:V),求区域内的电位分布。

要求用超松弛迭代法求解差分方程组进行计算。

代码:clc;clear;close all;hx=11;hy=11;v1=ones(hy,hx);n=10;v1(1,:)=ones(1,hx)*50;v1(hy,:)=ones(1,hx)*100;for i=1:hy;v1(i,1)=0;v1(i,hx)=100;end%计算超松弛因子w=2/(1+sin(pi/n));v2=v1;maxt=1;t=0;%初始化k=0while(maxt>1e-6)k=k+1;maxt=0;for i=2:hy-1for j=2:hx-1v2(i,j)=(1-w)*v1(i,j)+(v1(i,j+1)+v1(i+1,j)+v2(i-1,j)+v2(i,j-1))*w/4;%差分式t=abs(v2(i,j)-v1(i,j));if(t>maxt) maxt=t; endendendv1=v2endsubplot(1,2,1),mesh(v2);axis([0,11,0,11,0,100]);subplot(1,2,2),contour(v2,20);hold on结果与分析:作业2如图微带线, , 。

试用有限差分法求其有效介电常数 及特性阻抗 。

(保角变换法结果 =6.5, =49.5)代码:clc clearclose allhx=11; %设置x 方向网格节点数,间隔0.5w hy=11; %设置y 方向网格节点数v1=zeros(hy,hx); %设置二维数组,赋初值 v1(hy,:)=zeros(1,hx); %y=a 边界条件 v1(1,:)=zeros(1,hx); %y=0边界条件 v1(:,hx)=zeros(hy,1);%x=hx 边界条件 p=2;q=3;e=9.5;%相对介电常数 for j=1:p;v1(q,j)=1; %微带线边界条件 endmaxt=1;t=0; %设置误差和最大误差参量 v2=v1;n=0%迭代次数n 赋初值while(maxt>1e-6) %由v1迭代算出v2,精度为10-6n=n+1 %计算迭代次数maxt=0;for i=2:q-1; %对称轴v2(i,1)=(2*v1(i,2)+v1(i+1,1)+v2(i-1,1))/4;endfor i=q+1:hy-1;v2(i,1)=(2*v1(i,2)+v1(i+1,1)+v2(i-1,1))/4;endfor i=2:q-1; %均匀区域for j=2:hx-1;v2(i,j)=(v1(i,j+1)+v1(i+1,j)+v2(i,j-1)+v2(i-1,j))/4;endendfor j=p+1:hx-1; %介质分界面v2(q,j)=(v1(q,j+1)+2/(1+e)*v1(q+1,j)+v2(q,j-1)+2/(1+e)*e* v2(q-1,j))/4;endfor i=q+1:hy-1; %均匀区域for j=2:hx-1;v2(i,j)=(v1(i,j+1)+v1(i+1,j)+v2(i,j-1)+v2(i-1,j))/4;endendfor i=1:hy;for j=1:hx;t=abs(v2(i,j)-v1(i,j)); %收敛精度判据if(t>maxt)maxt=t;endendendv1=v2end结果与分析:填充介质时=9.5填充为空气时=1由有限差分法求出电位分布,进而推出电场强度分布,由高斯定理得中心带上单位长度的总电荷,最后由电荷推出电容C=Q/φ=(2.2164εrε0+1.3448ε0)C0=3.4504ε0εe=C/C0=(2.2164εrε0+1.3448ε0)/3.4504ε0=6.492Z0=1/(vp*C)=√εe/(c*C)=42.84Ω。

电磁场数值分析实验报告

电磁场数值分析实验报告

电器数值分析实验报告
姓名:班级:学号:1、电磁体设计(请将计算结果填在相应位置)
磁通量等势线:
沿Z轴的磁感应强度分布曲线:
磁感应强度分布云图:
2、同轴线电感分析(请将计算结果填在相应位置)
电感值:0.00026485
画出相应的场量云图:
等势线图:
电场强度幅度云图:
4、基于Ansoft软件的电磁场分析基本步骤
1.建模
2.材料属性
3.载荷及边界条件
4.网格剖分
5.计算机解算
6.后处理
5、列举一到两个实验中遇到的问题及解决的方法
问题:在实验中,做好了磁通量等势线图,磁感应强度分布云图后,不知道怎么得到沿Z轴的磁感应强度分布曲线。

解决方案:自己在百度上查找相关资料,然后又找班里学的比较好的的同学问问,知道可以在第一步建立模型时,在Z轴上画一条线段,继而再次对工程进行求解,然后通过在field对话框中选择line1得到沿Z轴的磁密B的分布曲线。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

注:考生属哪种类别请划“√”
(博士、在校硕士、工程硕士、师资硕士、同等学力、研究生班)
辽宁工程技术大学
研究生考试试卷
考试科目:电磁场数值分析
考生班级:电控研
考生姓名:
学号:
考试分数:
注意事项
1、考前研究生将上述项目填写清楚
2、字迹要清楚,保持卷面清洁
3、试题、试卷一齐交监考老师
4、教师将试题、试卷、成绩单,一起送研究生学院;
专业课报所在院、系
直流无刷电机的内部电磁分析
1提出问题
在电磁学里,电磁场是一种由带电物体产生的一种物理场。

处于电磁场的带电物体会感受到电磁场的作用力。

电磁场与带电物体之间的相互作用可以用麦克斯韦方程和洛伦兹力定律来描述。

电磁场是有内在联系、相互依存的电场和磁场的统一体的总称。

随时间变化的电场产生磁场,随时间变化的磁场产生电场,两者互为因果,形成电磁场。

电磁场可由变速运动的带电粒子引起,也可由强弱变化的电流引起,不论原因如何,电磁场总以光速向四周传播,形成电磁波。

电磁场是电磁作用的媒介,具有能量和动量,是物质存在的一种形式,电磁场的性质、特征及其运动变化规律由麦克斯韦方程确定。

ANSYS软件提供了图形用户界面与命令流两种方式来分析电机电磁场问题。

在电机电磁场计算中,命令流方式和图形用户界面方式相比,具有以下优点:通用性好,对于同系列、同型号的电机电磁场计算只要对电机的尺寸参数进行修改即可,而采用ANSYS的图形用户界面方式进行电机电磁场计算,每次计算都要重新输入图形,没有通用性;
通过合理应用ANSYS的APDL语言编写一个两重循环程序就可实现转子自动旋转和自动施加励磁电流的功能,与ANSYS的图形用户界面方式相比,减少了人机交互的次数,缩短了计算时间。

电机的电磁分析,常用的软件是Maxwell,他是一个功能强大、灵活的,融结构、热、流体、电磁、声学于一体的大型通用有限元分析软件。

广泛用于核工业、石油化工、航空航天、国防军工、机械制造、土木工程等一般工业及科学研究领域的设计分析。

本次作业中,将对直流无刷电机的内部电磁进行分析,采用Maxwell3D来建模,并进行磁场分析。

2直流无刷电机
直流无刷电机被广泛的用于日常生活用具、汽车工业、航空、消费电子、医学电子、工业自动化等装置和仪表。

顾名思义,直流无刷电机不使用机械结构的换向电刷而直接使用电子换向器,在使用中直流无刷电机相比有刷电机有许多
的优点,比如:
·能获得更好的扭矩转速特性;
·高速动态响应;
·高效率;
·长寿命;
·低噪声;
·高转速。

直流无刷电机由定子转子以及霍尔传感器组成。

直流无刷电机定子是由许多硅钢片经过叠压和轴向冲压而成,每个冲槽内都有一定的线圈组成了绕组,定子是2至8对永磁体按照N 极和S 极交替排列在转子周围构成的(内转子型),如果是外转子型直流无刷电机那么就是贴在转子内壁,与有刷直流电机不同,无刷直流电机使用电子方式换向。

要使无刷直流电机转起来,必须要按照一定的顺序给定子通电,那么我们就需要知道转子的位置以便按照通电次序给相应的定子线圈通电。

定子的位置是由嵌入到定子的霍尔传感器感知的。

通常会安排3个霍尔传感器在转子的旋转路径周围。

无论何时,只要转子的磁极掠过霍尔元件时,根据转子当前磁极的极性霍尔元件会输出对应的高或低电平,这样只要根据3个霍尔元件产生的电平的时序就可以判断当前转子的位置,并相应的对定子绕组进行通电。

3矢量泊松方程
1)选取第一类边界条件
边界上的物理条件直接规定了物理量u 在边界上的值,即
1Γ:()()G u s u s =1()
s ∈Γ式中,()G u s 为边界上给定的函数分布。

2)矢量泊松方程
为了达到矢量磁位A 与磁感应强度B 之间关系的相互单值性,对于稳态磁场,引入库仑规范
A ∇⨯=利用矢量恒等式()()()A A A ∇⨯∇⨯=∇∇-∇∇,则双旋度方程可转化为矢量泊松,即
()A J
ν∇∇=-对于直角坐标系,可写作
A A A J x x y y z z ννν⎡⎤⎛⎫∂⎡∂⎤∂∂∂⎡∂⎤⎛⎫⎛⎫++=-⎢⎥ ⎪ ⎪ ⎢⎥⎢⎥∂∂∂∂∂∂⎝⎭⎝⎭⎣⎦⎣⎦
⎝⎭⎣⎦写成分量形式,可得到三个方程
x x x x A A A J x x y
y z z ννν⎡⎤⎛⎫⎡∂⎤∂⎡∂⎤∂∂∂⎛⎫⎛⎫++=-⎢⎥ ⎪ ⎪ ⎢⎥⎢⎥∂∂∂∂∂∂⎝⎭⎝⎭⎣⎦⎣⎦⎝⎭⎣⎦y y y y A A A J x x
y y z z ννν⎡∂⎤⎡∂⎤⎡∂⎤⎛⎫⎛⎫⎛⎫∂∂∂++=-⎢⎥⎢⎥⎢⎥ ⎪ ⎪ ⎪∂∂∂∂∂∂⎝⎭⎝⎭⎝⎭⎣⎦⎣⎦⎣⎦z z z z A A A J x x y y z z
ννν⎡⎤⎛⎫⎡∂⎤∂⎡∂⎤∂∂∂⎛⎫⎛⎫++=-⎢⎥ ⎪ ⎪ ⎢⎥⎢⎥∂∂∂∂∂∂⎝⎭⎝⎭⎣⎦⎣⎦⎝⎭⎣⎦对于均匀线性介质,可简化为
2A=-J
ν∇4建立导体模型
采用maxwell3D 建立无刷直流电机3D 模型,通过静态磁场分析模型,该模型可调用模型库中的现有模型,其参数如下图:
图1参数图
图2模型图
5仿真计算磁场
1)定义材料属性
设背景为静电场
图3静电场
设置电机的材料为镍,其数据如下图
图4材料属性
2)施加边界条件与激励
设置激励电流100A。

设置计算域为10。

设置边界条件。

设置自适应计算参数:最大迭代次数15,误差要求百分之一,每次迭代加密剖分单元比例百分之三十。

图5设置激励电流
图6施加计算域
图7自适应计算参数
3)求解计算
进行检查并仿真后得到模型的收敛信息数据图和划分网格信息。

图8收敛图
图9网格表
对电机XY轴进行电磁分析
图10分析设定
图11开始分析
图12模型磁力线分布图。

相关文档
最新文档