电磁场有限元分析2解读ppt课件.ppt

《电磁场有限元分析2解读ppt课件.ppt》由会员分享，可在线阅读，更多相关《电磁场有限元分析2解读ppt课件.ppt（47页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、西安交通大学西安交通大学No.1 8/2/2022M&ISISchool of MEXian JiaotongUniversity低频电磁场有限元分析（低频电磁场有限元分析（ANSYS）孙岩桦孙岩桦 副教授副教授M&ISI, School of MEXian Jiaotong Univ.Xian, Shaanxi, P.R. China, 西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.23 2D 静态分析静态分析3.1 所用单元所用单元3.2 分析区域类型分析区域类型3.3 常用材料属性常用材料属性3.4 分析步

2、骤分析步骤 3.5 实例电磁作动器实例电磁作动器西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.3 3.1 所用单元所用单元单元类型形状和特点自由度PLANE134节点四边形3节点三角形每节点最多4自由度：AZ，位移，温度，VOLTPLANE538节点四边形6节点三角形每节点最多4自由度：AZ，电流，EMF，VOLT(1) 2D实体单元单元类型形状和特点自由度INFIN92节点线元AZ，位移，温度，VOLTINFIN1104节点或8节点4边形AZ，温度，电势(2) 远场单元西安交通大学西安交通大学 8/2/2022

3、M&ISISchool of MEXian JiaotongUniversityNo.4单元类型形状和特点自由度CIRCU124最多6节点每节点最多3自由度：电势，电流，EMF(3) 通用电路单元单元类型形状和特点自由度TARGE169目标线段无接触分析中对目标区域建模CONTA171面-面接触元，2节点AZ接触分析中对接触区域建模CONTA172面-面接触元，3节点AZ，接触分析中对接触区域建模CONTA175点-面接触元，1节点AZ，接触分析中对接触区域建模(4) 接触单元西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversi

4、tyNo.51 1）空气）空气 自由度：自由度：AZ AZ 材料属性：材料属性：r r ( (MURX),MURX),(RSVX)(RSVX)2 2）铁磁性物质）铁磁性物质 自由度：自由度：AZAZ 材料属性：材料属性：r r ( (MURX)MURX)或或B-HB-H曲线曲线, ,3 3）永磁体）永磁体 自由度：自由度：AZAZ 材料属性：材料属性：r r ( (MURX)MURX)或或B-HB-H曲线曲线, Hc, Hc (MGXX,MGYY) (MGXX,MGYY) 3.2 分析区域类型分析区域类型西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian Ji

5、aotongUniversityNo.6/PREP7HC=3000 ! Coercive forceBR=4000 ! Residual inductionTHETA=30 ! Permanent magnet orientation*AFUN,DEG ! Angular parametric functions in degreesMP,MGXX,2,HC ! X component of coercive force西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.7! B-H curve:TB,BH,2 ! B

6、-H curve for material 2TBPT,DEFI,-3000+HC,0 ! Shifted B-H curveTBPT,-2800+HC,500 ! First field defaults to DEFITBPT,-2550+HC,1000TBPT,-2250+HC,1500TBPT,-2000+HC,1800TBPT,-1800+HC,2000TBPT,-1350+HC,2500TBPT,-900+HC,3000TBPT,-425+HC,3500TBPT,0+HC,4000TBPLOT,BH,2 ! Plot of B vs. H西安交通大学西安交通大学 8/2/2022M

7、&ISISchool of MEXian JiaotongUniversityNo.8西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.94 4）电流激励多股线圈）电流激励多股线圈 自由度：自由度：AZAZ 材料属性：材料属性：r r 用用 BFEBFE命令施加命令施加单元单元电流密度电流密度JSJS5 5）电压激励多股线圈）电压激励多股线圈(PLANE53)(PLANE53) 自由度：自由度：AZAZ，CURRCURR 材料属性：材料属性：r r ( (MURX),MURX), (RSVX)(RSVX) 实常数：实

8、常数：CARE,TURN,LENG,DIRZ,FILLCARE,TURN,LENG,DIRZ,FILL 用用 BFEBFE命令施加电压命令施加电压VLTGVLTG 耦合区域内所有单元的耦合区域内所有单元的CURRCURR自由度，自由度，因为因为CURRCURR是流过每匝线圈的电流是流过每匝线圈的电流西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.10(RSVX) : (RSVX) : 各向同性各向同性 实常数：实常数：CARECARE：线圈横截面积：线圈横截面积TURNTURN：线圈匝数：线圈匝数LENGLENG：

9、线圈长度，用来计算电特性：线圈长度，用来计算电特性DIRZDIRZ：电流方向：电流方向FILLFILL：填充系数，影响线圈电阻：填充系数，影响线圈电阻西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.116 6）运动导体）运动导体 自由度：自由度：AZAZ 材料属性：材料属性：r r ( (MURX)MURX)或或B-HB-H曲线曲线, , (RSVX) (RSVX) （各向同性各向同性） 实常数：实常数：VELOX, VELOY,OMEGAZ,XLOC,YLOCVELOX, VELOY,OMEGAZ,XLOC,YL

10、OC 导体运动时区域的空间位置和属性不发生变化，典型应用：导体运动时区域的空间位置和属性不发生变化，典型应用：u实心转子感应电机实心转子感应电机u直线感应电机直线感应电机u涡流刹车系统涡流刹车系统u磁轴承磁轴承西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.12运动导体分析时所用实常数运动导体分析时所用实常数a)a)VELOX, VELOYVELOX, VELOY：全局笛卡尔坐标系中速度分量：全局笛卡尔坐标系中速度分量Z Zb)b)OMEGAZOMEGAZ：旋转角速度：旋转角速度, ,单位：单位：HzHzc)c)X

11、LOC,YLOCXLOC,YLOC：旋转轴在全局笛卡尔坐标系中位置：旋转轴在全局笛卡尔坐标系中位置分析时计算精度受集肤效应影响，与下面因素有关：分析时计算精度受集肤效应影响，与下面因素有关：a)a)网格剖分的精细程度网格剖分的精细程度b)b)相对磁导率相对磁导率c)c)电导率电导率d)d)速度速度可以用磁雷诺数来衡量（可以用磁雷诺数来衡量（数量级1.0之内有较高精度）：）：磁导率运动速度电阻率沿运动方向特征长度reVdMVd西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.137 7）不同网格边界）不同网格边界 自由

12、度：自由度：AZAZ 接触单元：接触单元： TARGE169, CONTA171, CONTA172, CONTA175TARGE169, CONTA171, CONTA172, CONTA175西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.14ANSYS ANSYS 材料库中预定义了如下材料的属性：材料库中预定义了如下材料的属性： 如果同一材料同时给定如果同一材料同时给定B-HB-H曲线和曲线和urur，那么使用，那么使用urur 可以指定各向异性相对磁导率，可以指定各向异性相对磁导率，MURX,MURY,MU

13、RZMURX,MURY,MURZ 可以同时用可以同时用urur和和B-HB-H曲线指定各方向磁导率，若曲线指定各方向磁导率，若urur=0,=0,则使用则使用B-HB-H曲线曲线mp,murx,2,1000 mp,mury,2,0! read B-H curve for material 2mp,murz,2,1000 3.3 常用材料属性常用材料属性材料特性材料属性文件名铜电阻率/温度曲线，uremagCopper.SI_MPLM3钢B-H曲线emagM3.SI_MPLM54钢B-H曲线emagM54.SI_MPLSA1010钢B-H曲线emagSa1010.SI_MPL硅钢B-H曲线ema

14、gSilicon.SI_MPL钴钢B-H曲线emagVanad.SI_MPL西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.15 设置分析环境设置分析环境(/prep7)(/prep7)： 建立分析模型、剖分网格并给每个区域设定建立分析模型、剖分网格并给每个区域设定属性参数属性参数 施加边界条件和载荷（或激励）施加边界条件和载荷（或激励） 求解求解 后处理（后处理（/post1/post1） 3.4 分析步骤分析步骤西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian Jiaoton

15、gUniversityNo.161）边界条件和载荷）边界条件和载荷a a）边界条件）边界条件 磁通垂直（磁通垂直（缺省缺省）：）： 磁通平行磁通平行: AZ=0: AZ=0 远场远场:INFIN9,INFIN110:INFIN9,INFIN110 周期性周期性: : 外部场外部场: :直接给定节点直接给定节点AZAZ值值b b）激励）激励 源电流密度（源电流密度（JSJS）：）：BFE, BFABFE, BFA 电压降（电压降（VLTGVLTG）：）： BFE, BFABFE, BFAc c）标记（）标记（FLAGFLAG） 力标记：力标记：FMAGBCFMAGBC宏宏可以直接给物体上施加麦克

16、斯韦应力和虚位可以直接给物体上施加麦克斯韦应力和虚位移表面标记，以识别需要计算电磁力和力矩的物体。移表面标记，以识别需要计算电磁力和力矩的物体。 无限大表面：用来表示指向无限大开区域的单元无限大表面：用来表示指向无限大开区域的单元d d）其它）其它 线电流（线电流（CSGXCSGX） 麦克斯韦表面：在包围铁磁体的空气单元上施加麦克斯韦表面：在包围铁磁体的空气单元上施加 磁虚位移：在和空气接触的铁磁体边界节点上施加磁虚位移：在和空气接触的铁磁体边界节点上施加西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.172） 分

17、析选项分析选项a a）分析类型：）分析类型：ANTYPE,STATICANTYPE,STATICb b）求解算法）求解算法 SparseSparse Frontal solver(Frontal solver(缺省缺省) ) JacobiJacobi Conjugate Gradient(JCG Conjugate Gradient(JCG) solver) solver JCG out-of-memory solverJCG out-of-memory solver Incomplete CholeskyIncomplete Cholesky Conjugate Gradient(ICCG

18、Conjugate Gradient(ICCG) solver) solver Preconditioned Conjugate Gradient solver (PCG)Preconditioned Conjugate Gradient solver (PCG) PCG out-of-memory solverPCG out-of-memory solverc c）注意事项）注意事项 一般采用一般采用SparseSparse或或Frontal solverFrontal solver 模型很大是，可以采用模型很大是，可以采用JCGJCG或或PCG solverPCG solver 如果是电压

19、激励或包含速度效应时，只能用如果是电压激励或包含速度效应时，只能用sparse,frontalsparse,frontal, , JCG,ICCG solverJCG,ICCG solver 电路激励时，只能用电路激励时，只能用SparseSparse或或Frontal solver Frontal solver 西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.183） 后处理后处理基本结果：基本结果： 节点电磁自由度节点电磁自由度 (AZ, CURR)间接结果：间接结果： 节点磁密节点磁密 (BX, BY, BS

20、UM) 节点磁场强度节点磁场强度 (HX, HY, HSUM) 节点电磁力节点电磁力 (FMAG: X, Y分量分量, SUM) 节点感生线电流节点感生线电流 (CSGZ) 单元源电流密度单元源电流密度 (JSZ) 单位体积焦耳热单位体积焦耳热 (JHEAT) . 等等.西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.19结果显示：结果显示：磁力线：等磁力线：等AZ线线等高线等高线向量图向量图电磁力、力矩：电磁力、力矩：电参数：电参数：线圈电阻线圈电阻电感电感西安交通大学西安交通大学 8/2/2022M&ISISc

21、hool of MEXian JiaotongUniversityNo.20其它宏：其它宏：CURR2D ：calculates current flow in a 2-D conductor.EMAGERR ：calculates the relative error in an electrostatic or electromagnetic field analysis.FLUXV ：calculates the flux passing through a closed contour.FMAGSUM ：summarizes electromagnetic force calculat

22、ions on element components.FOR2D ：calculates magnetic forces on a body.MMF：calculates magnetomotive force along a path.PLF2D ：generates a contour line plot of equipotentials.SENERGY ：determines the stored magnetic energy or co-energy.TORQ2D：calculates torque on a body in a magnetic field.TORQC2D：cal

23、culates torque on a body in a magnetic field based on a circular path.TORQSUM ：summarizes electromagnetic Maxwell and Virtual work torque calculations on element components for 2-D planar problems.西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.21假定：假定：n没有磁饱和没有磁饱和n周围没有漏磁周围没有漏磁 3.5 实例

24、电磁作动器实例电磁作动器西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.22n=650Number of turns in the coil; used in postprocessingI=1.0Current per turnta=.75Thickness of inner leg of magnetic circuittb=.75Thickness of lower leg of magnetic circuittc=.50Thickness of outer leg of magnetic circuitt

25、d=.75Armature thicknesswc=1Width of coilhc=2Height of coilgap=.25Gapspace=.25Space around coilws=wc+2*space hs=hc+.75 w=ta+ws+tcTotal width of modelhb=tb+hs h=hb+gap+tdTotal height of modelacoil=wc*hcCoil areajdens=n*i/acoilCurrent density of coil西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUn

26、iversityNo.23！设置分析环境！设置分析环境/PREP7 /TITLE, 2-D Solenoid Actuator Static Analysis ！定义单元类型和材料属性！定义单元类型和材料属性ET,1,PLANE53 ! Define PLANE 53 as element type KEYOPT,1,3,1 ! Use axisymmetric analysis option MP,MURX,1,1 ! Define material properties (permeability) MP,MURX,2,1000 ! Permeability of backiron MP,

27、MURX,3,1 ! Permeability of coil MP,MURX,4,2000 ! Permeability of armature 西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.24！定义常量和变量参数！定义常量和变量参数/com, ! Set parameter values for analysis n=650 ! Number of coil turns i=1.0 ! Current per turn ta=.75 ! Model dimensions (centimeters) tb=.

28、75 tc=.50 td=.75 wc=1 hc=2 gap=.25 space=.25 ws=wc+2*space hs=hc+.75 w=ta+ws+tc hb=tb+hs h=hb+gap+td acoil=wc*hc ! Cross-section area of coil (cm*2) jdens=n*i/acoil ! Current density (A/cm*2)西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.25！建立实体模型！建立实体模型/PNUM,AREA,1 RECTNG,0,w,0,tb

29、! Create rectangular areas RECTNG,0,w,tb,hb RECTNG,ta,ta+ws,0,hRECTNG,ta+space,ta+space+wc,tb+space,tb+space+hc AOVLAP,ALL RECTNG,0,w,0,hb+gap RECTNG,0,w,0,h AOVLAP,ALL NUMCMP,AREA ! Compress out unused area numbers APLOT 西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.26！设定区域材料属性！设定

30、区域材料属性ASEL,S,AREA,2 ! Assign attributes to coil AATT,3,1,1,0 ASEL,S,AREA,1 ! Assign attributes to armature ASEL,A,AREA,12,13 AATT,4,1,1 ASEL,S,AREA,3,5 ! Assign attributes to backiron ASEL,A,AREA,7,8 AATT,2,1,1,0 /PNUM,MAT,1 ! Turn material numbers on ALLSEL,ALL APLOT ! Plot areas 西安交通大学西安交通大学 8/2/2

31、022M&ISISchool of MEXian JiaotongUniversityNo.27！剖分单元网格！剖分单元网格SMRTSIZE,4 ! Set smart size meshing level 4 (fine) AMESH,ALL ! Mesh all areas ！定义力边界！定义力边界ESEL,S,MAT,4 ! Select armature elements CM,ARM,ELEM ! Define armature as a component FMAGBC,ARM ! Apply force boundary conditions to armature ALLSEL

32、,ALL ARSCAL,ALL,.01,.01,1,1 ! Scale model to MKS (meters) FINISH 西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.28！施加边界条件和载荷！施加边界条件和载荷/SOLU ESEL,S,MAT,3 ! Select coil elements BFE,ALL,JS,1,jdens/.01*2 ! Apply current density (A/m*2) ESEL,ALL NSEL,EXT ! Select exterior nodes D,ALL,AZ

33、,0 ! Set potentials to zero (flux-parallel) ALLSEL,ALL FINISH！求解！求解 /SOLU MAGSOLV ! Solve magnetic field SAVE FINISH 西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.29！后处理，检查计算结果！后处理，检查计算结果/POST1 PLF2D ! Plot flux lines in the model FMAGSUM ! Summarize magnetic forces PLVECT,B,VECT,E

34、LEM,ON ! Plot flux density as vectors /GRAPHICS,POWER ! Turn PowerGraphics on AVRES,2 ! Dont average results across materials PLNSOL,B,SUM ! Plot flux density magnitude FINISH 西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.304 2D 时谐分析时谐分析4.1 导体模型导体模型4.2 振幅、初相位和频率振幅、初相位和频率西安交通大学西安交通大

35、学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.31 正旋变化的交变磁场在导体内会感应出涡流，在时谐分析正旋变化的交变磁场在导体内会感应出涡流，在时谐分析中提供了更多的导体模型中提供了更多的导体模型西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.321 1）短路实心导体）短路实心导体 自由度：自由度：AZAZ 材料属性：材料属性：r r (RSVX)(RSVX) 产生的涡流在回路内流动，无电压降产生的涡流在回路内流动，无电压降4.1 导体模型导体模型西安交通

36、大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.332 2）开路实心导体）开路实心导体 自由度：自由度：AZAZ，VOLTVOLT 材料属性：材料属性：r r (RSVX)(RSVX) 所有的所有的VOLTVOLT自由度需耦合自由度需耦合西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.343 3）电流激励实心导体）电流激励实心导体 自由度：自由度：AZAZ，VOLTVOLT 材料属性：材料属性：r r (RSVX)(RSVX) 所有的所有的VO

37、LTVOLT自由度需耦合自由度需耦合, ,在其中一个节点上给定电流在其中一个节点上给定电流西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.354 4）电流激励多股线圈）电流激励多股线圈 自由度：自由度：AZAZ 材料属性：材料属性：r r 用用 BFEBFE命令施加命令施加单元单元电流密度电流密度JSJS5 5）电压激励多股线圈）电压激励多股线圈(PLANE53)(PLANE53) 自由度：自由度：AZAZ，CURRCURR 材料属性：材料属性：r r ( (MURX),MURX), (RSVX)(RSVX) 实常

38、数：实常数：CARE,TURN,LENG,DIRZ,FILLCARE,TURN,LENG,DIRZ,FILL 用用 BFEBFE命令施加电压命令施加电压VLTGVLTG 耦合区域内所有单元的耦合区域内所有单元的CURRCURR自由度，自由度，因为因为CURRCURR是流过每匝线圈的电流是流过每匝线圈的电流西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.366 6）单端共地多根实心导体）单端共地多根实心导体 自由度：自由度：AZAZ，VOLTVOLT 材料属性：材料属性：r r (RSVX)(RSVX) 所有的所有的

39、VOLTVOLT自由度需耦合自由度需耦合用来仿真鼠笼转子等器件，忽略端部效应用来仿真鼠笼转子等器件，忽略端部效应西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.377 7）电路激励多股线圈）电路激励多股线圈 自由度：自由度：AZAZ，CURR,EMFCURR,EMF 材料属性：材料属性：r r (RSVX)(RSVX) 实常数：实常数：CARE, TURN, LENG, DIRZ, FILL CARE, TURN, LENG, DIRZ, FILL 所有的所有的CURRCURR自由度需耦合自由度需耦合,EMF,EM

40、F自由度需耦合自由度需耦合西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.388 8）电路激励实心导体）电路激励实心导体 自由度：自由度：AZAZ，CURR,EMFCURR,EMF 材料属性：材料属性：r r (RSVX)(RSVX) 实常数：实常数：CARE, LENG, CARE, LENG, 所有的所有的CURRCURR自由度需耦合自由度需耦合,EMF,EMF自由度需耦合自由度需耦合西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.

41、399 9）叠片结构铁磁体）叠片结构铁磁体 自由度：自由度：AZAZ 材料属性：材料属性：r r ( (MURX)MURX)或或B-HB-H曲线曲线, , 忽略区域内的涡流忽略区域内的涡流1010）空气（同静态分析）空气（同静态分析）1111）运动导体（同静态分析）运动导体（同静态分析）西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.404.2 振幅、初相位和频率振幅、初相位和频率西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.415

42、2D 瞬态分析瞬态分析5.1 边界条件和载荷（激励）边界条件和载荷（激励）5.2 后处理后处理5.3 实例实例西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.425.1 振幅、初相位和频率振幅、初相位和频率 边界条件和载荷都随时间变化，可以按照其变化确定载荷边界条件和载荷都随时间变化，可以按照其变化确定载荷步，建立载荷文件。每一载荷步需指定：步，建立载荷文件。每一载荷步需指定： 边界条件或载荷以及对应的时间边界条件或载荷以及对应的时间 载荷斜坡变化或阶跃变化载荷斜坡变化或阶跃变化 是否自动时间步长是否自动时间步长,

43、 , 西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.435.2 后处理后处理1 1）POST1 POST1 某一时间点的分析结果某一时间点的分析结果SET,TIME 2) POST262) POST26，变量的时间历程结果，变量的时间历程结果载入结果文件：载入结果文件：RESUME 定义变量定义变量 西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.44假定：假定：n没有磁饱和没有磁饱和n周围没有漏磁周围没有漏磁 5.3 实例电磁作动器

44、实例电磁作动器西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.45 线圈由上图所示的电压驱动，通过分析，得到线圈由上图所示的电压驱动，通过分析，得到如下结果：如下结果： 线圈电流变化线圈电流变化 作用在衔铁上电磁力的变化作用在衔铁上电磁力的变化 线圈电感的变化线圈电感的变化 0 0.01 0.03 0.06 t(s)v(V)12 0step1 step2 step3西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.46！step1TIME

45、,.01DELTIM,.002OUTRES,ALL,ALLSOLVE！step2TIME,.03NSUBST,1SOLVE！step3TIME,.03DELTIM,.005TIME,.06NSUBST,1SOLVEFINISH西安交通大学西安交通大学 8/2/2022M&ISISchool of MEXian JiaotongUniversityNo.47！求电磁力、线圈电流和电感！求电磁力、线圈电流和电感/POST26PMGTRAN,ARM,COIL,COILFINISH 00.010.020.030.040.050.06-30-25-20-15-10-5051015T(s)F(N)00.010.020.030.040.050.060200400600800100012001400T(s)I(A)