ANSYS电磁场分析指南
第十一章磁宏
11.1 什么是电磁宏
电磁宏是ANSYS宏命令，其主要功能是帮助用户方便地建立分析模型、方便地获取想要观察的分析结果。目前，ANSYS提供了下列宏命令，可用于电磁场分析：
·CMATRIX：计算导体间自有和共有电容系数
·CURR2D：计算二维导电体内电流
·EMAGERR：计算在静电或电磁场分析中的相对误差
·EMF：沿预定路径计算电动力（emf）或电压降
·FLUXV：计算通过闭合回路的通量
·FMAGBC：对一个单元组件加力边界条件
·FMAGSUM：对单元组件进行电磁力求和计算
·FOR2D：计算一个体上的磁力
·HFSWEEP：在一个频率范围内对高频电磁波导进行时谐响应分析，并进行
相应的后处理计算
·HMAGSOLV：定义2-D谐波电磁求解选项并进行谐波求解
·IMPD：计算同轴电磁设备在一个特定参考面上的阻抗
·LMATRIX：计算任意一组导体间的电感矩阵
·MAGSOLV：对静态分析定义磁分析选项并开始求解
·MMF：沿一条路径计算磁动力
·PERBC2D：对2—D平面分析施加周期性约束
·PLF2D：生成等势的等值线图
·PMGTRAN：对瞬态分析的电磁结果求和
·POWERH：在导体内计算均方根（RMS）能量损失
·QFACT：根据高频模态分析结果计算高频电磁谐振器件的品质因子
·RACE：定义一个“跑道形”电流源
·REFLCOEF：计算同轴电磁设备的电压反射系数、驻波比、和回波损失
·SENERGY：计算单元中储存的磁能或共能
·SPARM：计算同轴波导或TE10模式矩形波导两个端口间的反射参数
·TORQ2D：计算在磁场中物体上的力矩
·TORQC2D：基于一个圆形环路计算在磁场中物体上的力矩
·TORQSUM：对2-D平面问题中单元部件上的Maxwell力矩和虚功力矩求和
本章对这些宏有详细描述。在ANSYS命令手册和理论手册对这些宏有更详细的描述。
下面的表格列出了这些电磁宏的使用范畴。
1这些宏也应用于静电场问题
2能用在通过界面单元INTER115连接的MVP区域11.2 使用电磁宏
电磁宏根据其实现的功能，可以分为如下四类：
·建模类·求解类
·后处理类·高频分析类磁宏11.2.1 建模类
有三个宏可用作帮助建模：RACE、PERBC2D和FMAGBC
1）RACE产生一个由条形和弧形基元（SOURCE36单元）组成的“跑道”形电流源
命令：RACE
GUI：Main Menu>Preprocessor>-Modeling-Create>Racetrack Coil
Main
Menu>Preprocessor>Loads>-Loads-Apply>-Magnetic-Excitation>Racetrack Coil
RACE宏要求的参数如前面图1所示。“跑道”由二个参数XC和YC定位，这些值是在工作平面内分别沿X和Y轴到线圈厚度中点的距离。执行该宏时，可以把构成线圈的这些SOURCE36单元定义为一个部件，将部件名作为该宏的一个输入参数即可。
2）PERBC2D宏通过生成两个周期性对称面所必须的约束方程或节点耦合来施加周期对称边界条件，调用该宏的方式如下：
命令：PERBC2D
GUI：Main
Menu>Preprocessor>Loads>-Loads-Apply>-Magnetic-Boundary>-Vector Poten- Periodic BCs
下面的图形描述了该宏的三种选项型式：
奇对称选项表示一个半周期对称条件，偶对称条件表示全周期对称条件（重复结
构）。
3）FMAGBC用于对单元组件施加Maxwell面标志和虚功边界条件：
命令：FMAGBC
GUI：Main Menu>Preprocessor>Loads>-Loads-Apply>-Magnetic-Flag> CompForce/Torq
Main Menu>Solution>Apply>Comp. Force/Torq
11.2.2 求解类
有四个宏可帮助求解：MAGSOLV, HMAGSOLV, CMATRIX和LMATRIX
11.2.2.1MAGSOLV
MAGSOLV宏对大多数静磁分析问题能很快地定义求解选项并开始求解。它可应用于2D和3D模型，标量势法、矢量势法和棱边单元求解方法，以及线性和非线性分析。该宏不需要用户使用MAGOPT命令和二步或三步求解顺序（一定情况所要求），它也允许定义非线性收敛标准，并提供选项来控制电流源Biot-Savart 积分的重新计算。
命令：MAGSOLV
GUI：Main Menu>Solution>-Solve-Electromagnet>-Static
Analysis-Opt&Solv
11.2.2.2HMAGSOLV
HMAGSOLV宏对谐波分析能很快定义求解选项和开始求解。对2-D模型，它使用磁矢势（MVP）方法求解。它可用于线性和非线性分析。对于非线性分析，此宏不需要定义二步求解过程，并允许用户自己定义收敛标准。
用下列方法之一，调用HMAGSOLV宏：
命令：HMAGSOLV
GUI：Main Menu>Solution>-Solve-Electromagnet>-Harmonic
Analys-Opt&Solv
11.2.2.3 CMATRIX
CMATRIX可计算“对地”和“集总”电容矩阵。“对地”电容值表示一个导体的电荷与导体对地电压之比。“集总”电容值表示二个导体之间的电容值。实例详见多导体系统求取电容和本手册电容计算实例（命令方法）。详细情况见ANSYS理论手册第五章。
用下面方法之一，调用CMATRIX宏命令：
命令：CMATRIX
GUI：Main Menu>Solution>Solve-Electromagnet>Capac Matrix
11.2.2.4 LMATRIX
LMATRIX宏可以计算任意线圈组中每个线圈的微分电感矩阵和总磁链。参见《ANSYS理论手册》第5章。
LMATRIX宏用于在静磁场分析的一个“工作点”上计算任意一组导体间的微分电感矩阵和磁链。“工作点”被定义为在系统上加工作（名义）电流所得到的解，该宏命令既可用于线性求解也可用于非线性求解。
必须用波前求解器来计算“工作点”的解。
LMATRIX宏的计算依赖于对工作点进行求解的过程中建立的多个文件。该宏在执行求解之前在这些文件前面加一个前缀OPER来重命名文件，并在完成求解后自动保存这些文件。用户自己也可以保存这些文件的拷贝以进行备份。
该宏命令返回一个N×N＋1矩阵参数，N×N部分表示N-绕组系统的微分电感值，此处N表示系统中的线圈数。N＋1列表示总磁链。第I行表示第I个线圈。另外，电感矩阵的值还以文本文件的格式输出，以供外部使用。文件中第一个列表表示每个线圈的磁链。第二个列表表示微分电感矩阵的上三角部分。
命令：LMATRIX
GUI：Main Menu>Solution>-Solve-Electromagnet>-Static
Analysis-Induct Matrix
在调用LMATRIX宏之前，还需要给线圈单元赋一个名义电流值。对于使用磁矢势（MVP）法或基于棱边元方法进行求解的静磁分析，可以使用BFV、BFA或BFE命令来给线圈单元赋名义电流（以电流密度的方式）。对于使用简化标势法（RSP）、差分标势法（DSP）和通用标势法（GSP）的静磁分析，可以使用SOURCE36单元的实常数来给线圈单元赋名义电流。
为了使用LMATRIX宏，必须事先用*DIM命令定义一个N阶数组，N为线圈数，数组的每行都表示一个线圈。数组的值等于线圈在工作点时每匝的名义电流值，且电流值不能为零，当确实有零电流时，可以用一个很小的电流值来近似。另外，还需用CM命令把每个线圈的单元组合成一个部件。每组独立线圈单元的部件名必须是用一个前缀后面再加线圈号来定义。一个线圈部件可由标量
（RSP/DSP/GSP）或矢量单元（MVP）混合组成，最重要的一点是这些单元的激励电流与前面数组中所描述的电流相同。
在LMATRIX宏中需定义一个用于保存电感矩阵的数组名，用LMATRIX宏的对称系数（symfac）来定义对称性。如果由于对称性而只建了n分之一部分模型，则计算出的电感乘以n就得到总的电感值。
当工作点位于BH曲线的弯点处时，切向磁导率变化最快，会导致计算的感应系数随收敛标准而变化。为了获得更加准确的解，收敛标准要定义得更加严格一些，不仅仅是缺省值1.0×10－3。一般在执行MAGSOLV命令时，选择1.0×10－4或1.0×10－5。
在使用LMATRIX命令前，不要施加（或删除）非均匀加载，非均匀加载由以下原因生成：
·自由度命令(D, DA,等)在节点或者实体模型上定义非0值
·带有非0约束的CE命令
不要在不包含在单元组件中的单元上施加任何载荷(如current)
下面的例子是一个3线圈系统，每个线圈的名义电流分别为1.2、1.5和1.7安/匝，其分析的命令流如下。在这个例子中，数组名为“curr”，线圈部件名前缀为“wind”，电感矩阵的计算值存贮在名为“ind”数组中。值得注意的是，在LMATRIX命令行中，这些名字必须用单引号引起来。
*dim，cur，，3！3个线圈系统数组
cur（1）=1.2！线圈1的名义电流为1.2安培/匝
cur（2）=1.5！线圈2的名义电流为1.5安培/匝
cur（3）=1.7！线圈3的名义电流为1.7安培/匝
esel，s……！选择线圈1的单元
cm，wind1，elem！给选出的单元赋予部件名wind1
esel，s……！选择线圈2的单元
cm，wind2，elem！给选出的单元赋予部件名wind2
esel，s……！选择线圈3的单元
cm，wind3，elem！给选出的单元赋予部件名wind3
symfac=2！对称系数
Imaxtrix，symfac，’wind’，’curr’，’ind’ ！计算微分电感矩阵和总磁链
*stat，ind！列出ind电感矩阵
11.2.2.5 下面是以命令流方式进行的一个计算电感矩阵的例子
该例计算一个二线圈系统（永磁电感器件）在非线性工作点下的微分电感矩阵和
总磁链，其示意图如下：
几何性质：x1=0.1, x2=0.1, x=0.1, y=0.1
材料性质：μ
r =1.0（空气），H
c
=25（永磁体），B-H曲线（永磁体，见输
入参数）
线圈1：名义电流＝0.25安/匝，匝数＝10
线圈2：名义电流＝0.125安/匝，匝数＝20
目标值：L11=4, L22=16, L12=8
命令流如下：
/batch,list
/title, Two-coil inductor with a permanent magnet /nopr
!geometry data
!
n=1! meshing parameter
x=0.1! width (x size) of core
y=0.1! hight of core, y size of window
z=1! thickness of iron in z direction
x1=0.1! width (x size) of coil 1
x2=0.1! width (x size) of coil 2
Hcy=25! coercive magnetic field in y direction
n1=10! number of turns in coil1
n2=20! number of turns in coil2
!
!excitation data used by LMATRIX.MAC
!
symfac=1! symmetric factor for inductance computation nc=2! number of coils
*dim,cur,array,nc! nominal currents of coils
*dim,coils,char,nc! names of coil components
!
cur(1)=0.25! nominal current of 1st coil
coils(1)="wind1"! name of coil 1 component
!
cur(2)=-0.125! nominal current of 2nd coil
coils(2)="wind2"! name of coil 2 component
!
!auxiliary parameters
!
mu0=3.1415926*4.0e-7
x3=x1+x2! x coordinate right to coil2 left
x4=x3+2*x! x coordinate right to core
x5=x4+x2! x coordinate right to coil2 right
x6=x5+x1! x coordinate right to coil1 right
js1=cur(1)*n1/(x1*y)! nominal current density of coil1 js2=cur(2)*n2/(x2*y)! nominal current density of coil2 !
/prep7
et,1,53
!
mp,murx,1,1! air/coil
mp,mgyy,2,Hcy! coercive term
Bs=2! saturation flux density
Hs=100! saturation magnetic field
TB,BH,2! core: H = Hs (B/Bs)^2; BS=2T;HS=100A/m
*do,qqq,1,20
B=qqq/10*Bs
tbpt,,Hs*(B/Bs)**2,B
*enddo
!
rect, 0,x1,0,y! coil1 left
rect,x1,x3,0,y! coil2 left
rect,x3,x4,0,y! core
rect,x4,x5,0,y! coil2 right
rect,x5,x6,0,y! coil1 right
!
aglue,all
!
asel,s,loc,x,x1/2! coil 1 volume attribute
aatt,1,1,1
asel,s,loc,x,x5+x1/2
aatt,1,2,1
asel,s,loc,x,x1+x2/2! coil 2 volume attribute
aatt,1,3,1
asel,s,loc,x,x4+x2/2
aatt,1,4,1
asel,s,loc,x,x3+x! iron volume attribute
aatt,2,5,1
asel,all
!
esize,,n
amesh,all
!
nsel,s,loc,x,x6! flux parallel Dirichlet at symmetry plain, x=x6 !! homogeneous Neumann flux normal at yoke, x=0
d,all,az,0
nsel,all
!
esel,s,real,,1! coil 1 left component
bfe,all,JS,,,,js1! unite current density in coil 1
!
esel,s,real,,2! coil 1 right component
bfe,all,JS,,,,-js1! return unite current density in coil 1 !
esel,s,real,,1,2
cm,coils(1),elem
!
esel,s,real,,3! coil 2 left component
bfe,all,JS,,,,js2! unite current density in coil 2
!
esel,s,real,,4! coil 2 right component
bfe,all,JS,,,,-js2! return unite current density in coil 2 !
esel,s,real,,3,4
cm,coils(2),elem
!
!
allsel
!
fini
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! /com
/com obtain operating solution
/com
!
/solu
cnvtol,csg,,1.0e-4
/out,scratch
solve
fini
!
/post1
!
/out
!
/com,
/com,
senergy,! Stored electromagnetic energy savelen=S_ENG
senergy,1! Co-energy
savelce=C_ENG
!
fini
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! compute inductance
lmatrix,symfac,"wind","cur","ind", ! compute inductance matrix and flux
!
/com
finish
你将得到如下结果：
SUMMARY OF STORED ENERGY CALCULATION
Load Step Number:1.
Substep Number:1.
Time:0.1000E+01
Material Number ofStored EnergyMaterial Description
NumberElements(J/m)
1.4.0.52360E-05LinearIsotrp...
2.1.-0.33314E+00Nonlin.MagnetIsotrp.
_____________________________________________________________________
T O T A L5.-0.33313E+00
Note:
The energy density for the active elements used in the energy
calculation is stored in the element item "MG_ENG" for display
and printing. The total stored energy is saved as parameter (S_ENG)
_____________________________________________________________________
SUMMARY OF COENERGY CALCULATION
Load Step Number:1.
Substep Number:1.
Time:0.1000E+01
Material Number ofCoenergyMaterial Description
NumberElements(J/m)
1.4.0.52360E-05LinearIsotrp...
2.1.0.33314E+00Nonlin.MagnetIsotrp.
_____________________________________________________________________ T O T A L5.0.33314E+00
Note:
The co-energy density for the active elements used in the co-energy calculation is stored in the element item "MG_COENG" for display and printing. The total coenergy is saved as parameter (C_ENG)
_____________________________________________________________________ ________________ LMATRIX SOLUTION SUMMARY ___________________
Flux linkage of coil1. =0.19989E+01
Flux linkage of coil2. =0.39978E+01
Self inductance of coil1. =0.39976E+01
Self inductance of coil2. =0.15989E+02
Mutual inductance between coils1. and2. =0.79948E+01
Inductance matrix is stored in array parameter ind( 2., 3.)
Inductance matrix is stored in file ind.txt
11.2.3 后处理类
1）TORQ2D沿一条预先定义好的路径用面积分计算在磁场内一个物体上的力矩。该宏需要定义通过围绕在感兴趣物体周围的空气单元的一条路径，如图6
所示。
定义路径用PPATH命令（Main Menu>General Postproc>Elec&Mag
Calc>Define Path），定义路径时按反时针方向选择节点，这样就可以得到力矩的正确方向。路径必须完全在物体周围的空气域内。该宏只适用于2D分析。
命令：TORQ2D
GUI：Main Menu>General Postproc>Elec & Mag Calc>Torque
2）TORQC2D沿以总座标原点为中心的圆周线（路径）用面积分计算在磁场内一个物体上的力矩，它常用于计算作用在旋转机械上的力矩：
命令：TORQC2D
GUI：Main Menu>General Postproc>Elec & Mag Calc>Circular Torq
TORQC2D宏根据用户定义的半径（RAD）自动生成所需的路径。为了确保得到正确的解，在使用该宏前，只能选择感兴趣圆形体周围的空气单元，如图7所示。
3）FOR2D沿一条预先定义好的路径用面积分计算在磁场内一个物体上的力，该宏要求路径定义在感兴趣物体四周的空气中，使用下面的命令或GUI路径实现
FOR2D：
命令：FOR2D
GUI：Main Menu>General Postproc>Elec & Mag Calc>Mag Forces
图8显示二个平面分析实例的路径定义。定义路径用PPATH命令（Main Menu>General Postproc>Elec&Mag Calc>Define Path），定义路径时按反时针方向选择节点，这样就可以得到力的正确方向。
4）MMF宏计算磁动力，磁动力就是沿一条预先定义好的路径（用PATH、PPATH 命令或其等效GUI路径定义）对磁场强度H进行线积分。定义路径时按反时针方向选择节点，这样就可以得到MMF的正确方向。MMF宏自动设置PMAP命令的“ACCURATE”映射和“MAT”不连续项。执行宏命令后，ANSYS程序保持PMAP命令的这些设置。如果路径跨越多种材料，则每种材料应至少有1个路径点（如图9（b））。对于高频分析，MMF宏仅对TEM或TM模式有效。
命令：MMF
GUI：Main Menu>General Postproc>Elec & Mag Calc>2D and 3D-MMF
5）EMF宏计算电动力（emf），电动力就是沿一条预先定义好的路径（用PATH、PPATH命令或其等效GUI路径定义）对电场强度E或电压降进行线积分。它能用在2-D和3-D电流传导分析、静电场分析以及高频电磁场分析中。对于高频分析，EMF宏仅对TEM或TE模式有效。计算出的电动力保存在EMF参数中。在调用EMF 宏之前，必须先定义一条路径，该宏采用计算出的电场和路径操作进行电动力计算，当宏执行完毕后，所有路径项都将被清除。EMF宏自动设置PMAP命令的“ACCURATE”映射和“MAT”不连续项。执行宏命令后，ANSYS程序保持PMAP命令的这些设置。
命令：EMF
GUI：Main Menu>General Postproc>Elec & Mag Calc>-2D and 3D-EMF
6）POWERH宏在谐波分析中计算一个导体内的时间平均能量损失。在调用宏之前，必须先选择要进行计算的导体区域的单元。当导体区的单元足够细密时，该宏的计算结果最准确。
命令：POWERH
GUI：Main Menu>General Postproc>Elec & Mag Calc>-2D and 3D-Power Loss
7）FLUXV宏计算通过一个预定义回路的通量。在2-D分析中，路径只少要2个点定义。在3-D MVP分析中，路径必须为一条封闭围线，即第1点和最后一点必须是同一点。定义路径时按反时针方向选择点，这样就可以得到通量的正确方向。下图描述了在2-D和3-D分析中调用FLUXV宏时的路径选择。该宏只能用于磁矢量势法（MVP）的分析中。
命令：FLUXV
GUI：Main Menu>General Postproc>Elec & Mag Calc>Path Flux
8）PLF2D宏显示自由度AZ的等值线，在轴对称分析中，显示的等值线为：
半径×AZ=常数。此宏仅适用于2-D分析，等值线是平行于通量线的，它很好地描述了磁通图形。
命令：PLF2D
GUI：Main Menu>General Postproc>Plot Results>2D Flux Lines
Utility Menu>Plot >Results>Flux Lines
9）SENERGY宏计算模型中贮存的磁能和共能。能量密度贮存在单元表中，供图形显示和列表显示用。
图11说明如何确定非永磁材料的磁能和共能。
对于永磁体，磁能和共能按照如下计算：
·能量是曲线右边部分 (参见图12 (a))，注意图示中能量是负值。
·共能是曲线下部分(参见
图12 (b)).
第六章３-Ｄ静态磁场分析（棱边单元方法） 6.1何时使用棱边元方法 在理论上，当存在非均匀介质时，用基于节点的连续矢量位A来进行有限元计算会产生不精确的解，这种理论上的缺陷可通过使用棱边元方法予以消除。这种方法不但适用于静态分析，还适用于谐波和瞬态磁场分析。在大多数实际3-D 分析中，推荐使用这种方法。在棱边元方法中，电流源是整个网格的一个部分，虽然建模比较困难，但对导体的形状没有控制，更少约束。另外也正因为对电流源也要划分网格，所以可以计算焦耳热和洛伦兹力。 用棱边元方法分析的典型使用情况有： ·电机 ·变压器 ·感应加热 ·螺线管电磁铁 ·强场磁体 ·非破坏性试验 ·磁搅动 ·电解装置 ·粒子加速器 ·医疗和地球物理仪器 《ANSYS理论手册》不同章节中讨论了棱边单元的公式。这些章节包括棱边分析方法的概述、矩阵列式的讨论、棱边方法型函数的信息。 对于ANSYS的SOLID117棱边单元，自由度是矢量位A沿单元边切向分量的积分。物理解释为：沿闭合环路对边自由度（通量）求和，得到通过封闭环路的磁通量。正的通量值表示单元边矢量是由较低节点号指向较高节点号（由单元边连接）。磁通量方向由封闭环路的方向根据右手法则来判定。 在ANSYS中，AZ表示边通量自由度，它在MKS单位制中的单位是韦伯（Volt·Secs），SOLID117是20节点六面体单元，它的12个边节点（每条边
的中间节点）上持有边通量自由度AZ。单元边矢量是由较低节点号指向较高节点号。在动态问题中，8个角节点上持有时间积分电势自由度VOLT。 ANSYS程序可用棱边元方法分析3-D静态、谐波和瞬态磁场问题。（实体模型与其它分析类型一样，只是边界条件不同），具体参见第7章，第8章。 6.2单元边方法中用到的单元 表 1三维实体单元 6.3物理模型区域的特性与设置 对于包括空气、铁、永磁体、源电流的静态磁场分析模型，可以通过设置不同区域不同材料特性来完成。参见下表，详情在后面部分叙述。
基于ANSYS软件的电机电磁场有限元分析 发表时间：2007-9-11 作者: 黄劭刚夏永洪张景明来源: 万方数据 关键字: APDL语言同步发电机电磁场有限元 介绍了应用ANSYS自带的APDL编程语言进行软件开发，将该软件应用于同步发电机空载磁场分析中，在电机的电磁场计算中实现了电机的自动旋转、自动施加载荷的功能，使用、修改方便，并且计算速度快。通过对电磁场计算结果的后处理，得出了同步发电机的旋转磁场波形和电压波形。样机测试结果验证了分析结果的正确。 1 前言 ANSYS软件是一个功能强大、灵活的，融结构、热、流体、电磁、声学于一体的大型通用有限元分析软件。广泛用于核工业、石油化工、航空航天、机械制造、土木工程等一般工业及科学研究领域的设计分析。 在实际的电机电磁场分析中，电机的转子磁极形状、定子齿槽形状、气隙大小以及铁磁材料均已确定，但是当转子相对十定子齿槽的位置不同时一，其计算结果也不相同。为了分析电机电磁场问题，若把定、转子相对位置固定不变进行求解，再对电磁场计算结果进行傅立叶级数分解来计算电机绕组的电势则误差太大。为此，需要对定、转子不同位置时一分别进行计算，然后通过电磁场的计算结果求出电机何个定子齿部磁通随转角变化的关系，然后根据磁通的变化率求出电机基波绕组的电势。ANSYS软件是目前应用最为广泛、使用最方便的通用有限元分析软件之一，应用ANSYS软件来分析电机电磁场是非常有效的。但是当采用ANSYS软件的图形用户界面( GUI)操作方式时，每次定、转子之间的旋转、网格剖分、施加载荷进行求解、查看计算结果等都需要人工进行重复操作，使用起来非常繁琐，并且效率低。为此，木文采用ANSYS软件的APDL语言编写的软件对同步发电机的空载磁场进行研究，实现了电机定、转子之间的自动旋转，自动网格剖分，自动施加载荷以及自动求解的功能。整个电磁场分析过程无需人工进行干预，使用方便，便于修改，并且大大提高了计算速度。通过对同步发电机电磁场计算结果进行后处理，得出了同步发电机的旋转磁场波形和电压波形。 2 软件实现 ANSYS软件提供了图形用户界面与命令流两种方式来分析电机电磁场问题。在电机电磁场计算中，命令流方式和图形用户界面方式相比，具有以下优点:(1)通用性好，对于同系列、同型号的电机电磁场计算只要对电机的尺寸参数进行修改即可，而采用ANSYS的图形用户界面方式进行电机电磁场计算，每次计算都要重新输入图形，没有通用性;(2)通过合理应用ANSYS的APDL语言编写一个两重循环程序就可实现转子自动旋转和自动施加励磁电流的功能，与ANSYS 的图形用户界面方式相比，减少了人机交互的次数，缩短了计算时间。 2.1软件编写
第九章3-D静态、谐波和瞬态分析（节点法） 9.1节点法（MPV）进行3D静态磁场分析 3-D节点法磁场分析的具体步骤与2-D静态分析类似，选择GUI参数路径Main Menu>Preferences> Magnetic-Nodal，便于使用相应的单元和加载。与2-D静态分析同样的方式定义物理环境，但要注意下面讨论的存在区别的地方。 9.1.1 选择单元类型和定义实常数 对于节点法3 –D静磁分析，可选的单元为3D 矢量位SOLID97单元，与2D单元不同。自由度为：AX,AY,AZ。3D矢量位方程中，用INFIN111远场单元（AX、AY、AZ三个自由度）来为无限边界建模。 对于载压和载流绞线圈（只有SOLID97单元），必须定义如下实常数：
速度效应 可求解运动物体在特定情况下的电磁场，2-D静磁分析讨论了运动体的应用和限制，在3-D中，只有SOLID97单元类型能通过设置单元KEYOPT选项来考虑速度效应。 9.1.2 定义分析类型 用与2D静态磁场分析相同的方式定义3D静态磁场分析，即，可以通过菜单路径Main Menu>Solution>New Analysis、或者用命令ANTYPE,STATIC,NEW来定义一个新的静态磁场分析；或者用ANTYPE,STATIC,REST 命令来重启动一个3-D分析。 如果使用了速度效应，不能在3D静态分析(ANTYPE,STATIC)中直接求解具有速度效应的静态直流激励场，而要用具有很低频率的时谐分析(ANTYPE,HARMIC)来完成。 9.1.3 选择方程求解器 命令：EQSLV GUI：Main Menu>Solution>Analysis Options 3D模型建议使用JCG或PCG法进行求解。而对于载压模型、载流模型、或有速度效应的具有非对称矩阵的模型，只能使用波前法、JCG法、或ICCG法求解。电路激励模型只能用稀疏矩阵法或波前法求解器。 9.1.4 加载和求解 进入求解器： 命令：/SOLU GUI：Main Menu>Solution 3-D 静态MVP分析的载荷与2-D静态分析稍有些不同，但其菜单路径是一样的。下面是关于3-D静态磁场分析的一些加载： 该载荷用以定义磁力线垂直、磁力线平行、远场、周期边界条件、以及强加外磁场等载荷和边界条件，下表描述了各种边界条件相应的磁矢量势值：
(1) 如果你模拟结构体中裂缝扩展过程的模拟，在Ansys中可以用全解耦损伤分析方法来近似模拟裂缝扩展，我曾用Ansys软件中提供的可以定义10,000个材料参数和单元ekill/alive 功能完成了层状路面体中表面裂缝和反射裂缝在变温作用下的扩展过程的模拟。我模拟的过程相对来说比较简单，模拟过程中我们首先要知道裂缝的可能扩展方向，这样在裂缝可能扩展的带内进行网格加密处理，加密到什么程度依据计算的问题来确定。 (2) 如果采用断裂力学理论计算含裂缝结构体的应力强度因子，建模时只需在裂尖通过命令kscon生成奇异单元即可。Ansys模块中存在的断裂力学模块可以计算I、II、III型应力强度因子(线弹性断裂力学)和J积分(弹塑性断裂力学)，在Ansys中verification里面有一个计算I型应力强度因子的例子vm143，参见该例子就可以了。 (3) 如果通过断裂力学模拟裂缝的扩展过程，需要采用动态网格划分，这方面我没有做，通过Ansys的宏命令流应该可以实现。技术参考可参阅文献：杨庆生、杨卫.断裂过程的有限元模拟.计算力学学报，1997,14(4). (4) 我现在做动荷载作用下路面结构体中应力强度因子的分布规律，我是通过位移插值得到不同时间点处的应力强度因子。如果想这样做，可参阅理论参考中关于应力强度因子计算说明。 1. 讨论两种Ansys求极限荷载的方法 （1）力加载 可以通过对应的方法（比如说特征值屈曲）估计结构的极限荷载的大致范围，然后给结构施加一个稍大的荷载，打开自动荷载步二分法进行非线性静力分析，最后计算会因不收敛终止，则倒数第二个子步对应的就是结构的极限荷载；另外，也可以选择弧长法，采用足够的子步（弧长法可以一直分析到极限承载力之后的过程）同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。 （2）位移加载 给结构施加一个比较大的位移，打开自动荷载步二分法进行非线性分析，保证足够的子步数，这样也可以分析到极限荷载以后，通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。 希望众高手讨论一下 （1）弧长法求极限荷载的收敛性问题，如何画到荷载位移曲线的下降段？ （2）位移法求极限荷载的具体步骤？ 2. 需要注意的问题 1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据，因此在很多情况下会因为应力集中而使混凝土提前破坏，从而和试验结果不相吻合，因此，在实际应用过程中应该对单元分划进行有效控制，根据作者经验，当最小单元尺寸大于5cm 时，就可以有效避免应力集中带来的问题； 2. 支座是另一个需要注意的问题。在有限元分析中，很多时候约束是直接加在混凝土节点上，这样很可能在支座位置产生很大的应力集中，从而使支座附近的混凝土突然破坏，造成求解失败。因此，在实际应用过程中，应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块，避免支座的应力集中；
第十六章电路分析 16.1 什么是电路分析 电路分析可以计算源电压和源电流在电路中引起的电压和电流分布。分析方法由源的类型来决定： 源的类型分析方法 交流（AC）谐波分析 直流（DC）静态分析 随时间变化瞬态分析 要在电磁学分析中用有限元来模拟全部电势，就必须提供足够的灵活性来模拟载流电磁设备。ANSYS程序对于电路分析有如下性能： ·用经过改进的基于节点的分析方法来模拟电路分析 ·可以将电路与绕线圈和块状导体直接耦合 ·2-D和3-D模型都可以进行耦合分析 ·支持直流、交流和时间瞬态模拟 ANSYS程序中先进的电路耦合模拟功能精确地模拟多种电子设备，： ·螺线管线圈 ·变压器 ·交流机械 16.2 使用CIRCU124单元 ANSYS提供一种通用电路单元CIRCU124对线性电路进行模拟，该单元求解未知的节点电压（在有些情况下为电流）。电路由各种部件组成，如电阻、电感、互感、电容、独立电压源和电流源、受控电压源和电流源等，这些元件都可以用CIRCU124单元来模拟。 注：本章只描述CIRCU124单元的某些最重要的特性，对该单元的详细描述参见《ANSYS 单元手册》。 16.2.1 可用CIRCU124单元模拟的电路元件
对CIRCU124单元通过设置KEYOPT（1）来确定该单元模拟的电路元件，如下表所示。例如，把KEYOPT（1）设置为2，就可用CIRCU124来模拟电容。对所有的电路元件，正向电流都是从节点I流向节点J。 表1CIRCU124单元能模拟的电路元件 注意：全部的电路选项如上表和下图图1所示，ANSYS的电路建模程序自动生成下列实常数：R15（图形偏置，GOFFST）和R16（单元识别号，ID）。本章下一节将详细讨论电路建模程序。 下图显示了利用不同的KEYOPT（1）设置建立的不同电路元件，那些靠近元件标志的节点是“浮动”节点（即它们并不直接连接到电路中）。 16.2.2 CIRCU124单元的载荷类型 对于独立电流源和独立电压源可用CIRCU124单元KEYOPT（2）选项来设置激励形式，可以定义电流或电压的正弦、脉冲、指数或分段线性激励。详细的载荷函数图和相应的实常数请参见《ANSYS单元手册》。 16.2.3 将FEA（有限元）区耦合到电路区 可将电路分析的三种元件耦合到FEA区，图2所示的这三种元件直接连接到有限元模型的导体上（耦合是在矩阵中进行耦合的，因此只能为线性的）：
回旋加速器 在一般电磁场分析中关心的典型的物理量为： -磁通密度?能量损耗 -磁场强度?磁漏 ?磁力及磁矩? S-参数 ?阻抗?品质因子Q ?电感?回波损耗 ?涡流?本征频率 存在电流、永磁体和外加场都会激励起需要分析的磁场 1.2ANSYS 如何完成电磁场分析计算 ANSYSU Maxwell 方程组作为电磁场分析的出发点。有限元方法计算的未知 量（自由度）主要是磁位或通量，其他关心的物理量可以由这些自由度导出。根 ANSY 电磁场分析指南第一章 发表时间： 2007-9-20 作者 : 安世亚太 来源 : e-works 关键字 : ANSYS 电磁场分析 CAE 教程 第一章磁场分析概述 1.1 磁场分析对象 利用ANSYS/Ema 或ANSYS/Multiphysics 模块中的电磁场分析功能，ANSYS 可分析计算下列的设备中的电磁场，如： 电力发电机 磁带及磁盘驱动器 变压器 波导 螺线管传动器 谐振腔 电动机 连接器 磁成像系统 天线辐射 图像显示设备传感器 滤波器
据用户所选择的单元类型和单元选项的不同， ANSYS+算的自由度可以是标量磁 位、矢量磁位或边界通量。 1.3 静态、谐波、瞬态磁场分析 利用ANSY 测以完成下列磁场分析： ?2-D 静态磁场分析，分析直流电（DC ）或永磁体所产生的磁场，用矢量位方 程。参见本书“二维静态磁场分析” ?2-D 谐波磁场分析，分析低频交流电流（AC ）或交流电压所产生的磁场，用 矢量位方程。参见本书“二维谐波磁场分析” ?2 -D 瞬态磁场分析，分析随时间任意变化的电流或外场所产生的磁场，包 含永磁体的效应，用矢量位方程。参见本书“二维瞬态磁场分析” ?3-D 静态磁场分析，分析直流电或永磁体所产生的磁场，用标量位方法。 参见本书“三维静态磁场分析（标量位方法）” ?3-D 静态磁场分析，分析直流电或永磁体所产生的磁场，用棱边单元法。 参见本书“三维静态磁场分析（棱边元方法）” ?3-D 谐波磁场分析，分析低频交流电所产生的磁场，用棱边单元法。建议 尽量用这种方法求解谐波磁场分析。参见本书“三维谐波磁场分析（棱边元方法） ?3-D 瞬态磁场分析，分析随时间任意变化的电流或外场所产生的磁场，用 棱边单元法。建议尽量用这种方法求解谐波磁场分析。参见本书“三维瞬态磁场 分析（棱边元方法）” 3-D 静态磁场分析，用矢量位方法。参见“基于节点方法 ?基于节点方法的3-D 谐波磁场分析，用矢量位方法。参见“基于节点方法 的 3-D 谐波磁场分析” 1.4 关于棱边单元、标量位、矢量位方法的比较 什么时候选择2-D 模型，什么时候选择3-D 模型？标量位方法和矢量位方 法有何不同？棱边元方法和基于节点的方法求解 3-D 问题又有什么区别？在下面 将进行详细比较。 1.4.12-D 分析和 3-D 分析比较 3-D 分析就是用 3-D 模型模拟被分析的结构。现实生活中大多数结构需要 3- D 模型来进行模拟。然而3-D 模型对建模的复杂度和计算的时间都有较高要求。 所以，若 ?基于节点方法的 的 3-D 静态磁场分析” ?基于节点方法的 的 3-D 瞬态磁场分析” 3-D 瞬态磁场分析， 用矢量位方法。参见“基于节点方法
一学习ANSYS需要认识到的几点 相对于其他应用型软件而言，ANSYS作为大型权威性的有限元分析软件，对提高解决问题的能力是一个全面的锻炼过程，是一门相当难学的软件，因而，要学好ANSYS，对学习者就提出了很高的要求，一方面，需要学习者有比较扎实的力学理论基础，对ANSYS分析结果能有个比较准确的预测和判断，可以说，理论水平的高低在很大程度上决定了ANSYS使用水平；另一方面，需要学习者不断摸索出软件的使用经验不断总结以提高解决问题的效率。在学习ANSYS的方法上，为了让初学者有一个比较好的把握，特提出以下五点建议：（1）将ANSYS的学习紧密与工程力学专业结合起来 毫无疑问，刚开始接触ANSYS时，如果对有限元，单元，节点，形函数等《有限元单元法及程序设计》中的基本概念没有清楚的了解话，那么学ANSYS很长一段时间都会感觉还没入门，只是在僵硬的模仿，即使已经了解了，在学ANSYS之前，也非常有必要先反复看几遍书，加深对有限元单元法及其基本概念的理解。 作为工程力学专业的学生，虽然力学理论知识学了很多，但对许多基本概念的理解许多人基本上是只停留于一个符号的认识上，理论认识不够，更没有太多的感性认识，比如一开始学ANSYS时可能很多人都不知道钢材应输入一个多大的弹性模量是合适的。而在进行有限元数值计算时，需要对相关参数的数值有很清楚的了解，比如材料常数，直接关系到结果的正确性，一定要准确。实际上在学ANSYS时，以前学的很多基本概念和力学理论知识都忘得差不多了，因而遇到有一
定理论难度的问题可能很难下手，特别是对结果的分析，需要用到《材料力学》，《弹性力学》和《塑性力学》里面的知识进行理论上的判断，所以在这种情况下，复习一下《材料力学》，《弹性力学》和《塑性力学》是非常有必要的，加深对基本概念的理解，实际上，适当的复习并不要花很多时间，效果却很明显，不仅能勾起遥远的回忆，加深理解，又能使遇到的问题得到顺利的解决。 在涉及到复杂的非线性问题时（比如接触问题），一方面，不同的问题对应着不同的数值计算方法，求解器的选择直接关系到程序的计算代价和问题是否能顺利解决；另一方面，需要对非线性的求解过程有比较清楚的了解，知道程序的求解是如何实现的。只有这样，才能在程序的求解过程中，对计算的情况做出正确的判断。因此，要能对具体的问题选择什么计算方法做出正确判断以及对计算过程进行适当控制，对《计算方法》里面的知识必须要相当熟悉，将其理解运用到ANSYS的计算过程中来，彼此相互加强理解。要知道ANSYS是基于有限元单元法与现代数值计算方法的发展而逐步发展起来的。因此，在解决非线性问题时，千万别忘了复习一下《计算方法》。此外，对《计算固体力学》也要有所了解（一门非常难学的课），ANSYS对非线性问题处理的理论基础就是基于《计算固体力学》里面所讲到的复杂理论。 作为学工程力学的学生，提高建模能力是非常急需加强的一个方面。在做偏向于理论的分析时，可能对建模能力要求不是很高，但对于实际的工程问题，有限元模型的建立可以说是一个最重要的问题，而后
期末大作业 题目：简单直流致动器 ANSYS电磁场分析及与ansoft仿真分析结果比较作者姓名：柴飞龙 学科（专业）：机械工程 学号：21225169 所在院系：机械工程学系 提交日期2013 年 1 月
1、 背景简述： ANSYS 软件是融结构、流体、电场、磁场、声场分析于一体的大型通用软件有限元分析软件，是现代产品设计中的高级CAE 工具之一。而ansoft Maxwell 软件是一款专门分析电磁场的分析软件，如传感器、调节器、电动机、变压器等。 本人在实验室做的课题涉及到电机仿真，用的较多的是ansoft 软件，因为其对电机仿真的功能更强大，电机功能模块更多，界面友好。 现就对一电磁场应用实例，用ANSYS 进行仿真分析，得到的结果与ansoft 得到的结果进行简单核对比较。 2、 问题描述： 简单直流致动器由2个实体圆柱铁芯，中间被空气隙分开的部件组成，线圈中心点处于空气隙中心。衔铁是导磁材料，导磁率为常数（即线性材料，r μ=1000），线圈是可视为均匀材料，空气区为自由空间（1=r μ），匝数为2000，线圈励磁为直流电流：2A 。模型为轴对称。 3、 ANSYS 仿真操作步骤： 第一步：Main menu>preferences
第二步：定义所有物理区的单元类型为PLANE53 Preprocessor>Element type>Add/Edit/Delete 第三步：设置单元行为 模拟模型的轴对称形状，选择Options（选项） 第四步：定义材料 Preprocessor>Material Props> ?定义空气为1号材料(MURX = 1) ?定义衔铁为2号材料(MURX = 1000) ?定义线圈为3号材料(自由空间导磁率，MURX=1)
三维螺线管静态磁场分析 要求计算螺线管，如下图所示，衔铁所受磁力，线圈为直流激励，产生力驱动衔铁。线圈电流为6A，500匝。由于对称性，只分析1/4的模型，如图1所示： 图1螺线管制动器 在仿真分析时，空气相对磁导系数为1.0；使用智能网格划分（LVL=8）；设定全部面为通量平行，这是自然边界条件，自动得到满足。因为是采用的1/4对称模型，所以磁力的计算结果要乘以4。
施加边界条件： ! /SOLU D,2,MAG,0 ! !SOLVE ! ALLSEL,ALL MAGSOLV,3,,,,,1 FINISH ! 建立的模型如下图所示：
对模型进行智能网格划分，如下图所示： 仿真分析所得磁场强度分布图为：
衔铁所受磁力分布图为： 衔铁所受磁力分布图为：
计算所得衔铁所受磁力为： SUMMARY OF FORCES BY VIRTUAL WORK Load Step Number: 2. Substep Number: 1. Time: 0.2000E+01 Units of Force: ( N ) Component Force-X Force-Y Force-Z ARM 0.14339E+02 0.11359E+02 -0.12846E+02 ___________________________________________________ SUMMARY OF FORCES BY MAXWELL STRESS TENSOR Units of Force: ( N ) Component Force-X Force-Y Force-Z ARM 0.79007E+01 0.55769E+01 -0.11511E+02 _____________________________________________________ Note: Maxwell forces are in the Global Cartesian coordinate system. Virtual work forces are in the element ESYS coordinate system.
1. 讨论两种Ansys求极限荷载的方法 （1）力加载 可以通过对应的方法（比如说特征值屈曲）估计结构的极限荷载的大致范围，然后给结构施加一个稍大的荷载，打开自动荷载步二分法进行非线性静力分析，最后计算会因不收敛终止，则倒数第二个子步对应的就是结构的极限荷载；另外，也可以选择弧长法，采用足够的子步（弧长法可以一直分析到极限承载力之后的过程）同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。 （2）位移加载 给结构施加一个比较大的位移，打开自动荷载步二分法进行非线性分析，保证足够的子步数，这样也可以分析到极限荷载以后，通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。 希望众高手讨论一下 （1）弧长法求极限荷载的收敛性问题，如何画到荷载位移曲线的下降段？ （2）位移法求极限荷载的具体步骤？ 2. 需要注意的问题 1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据，因此在很多情况下会因为应力集中而使混凝土提前破坏，从而和试验结果不相吻合，因此，在实际应用过程中应该对单元分划进行有效控制，根据作者经验，当最小单元尺寸大于5cm 时，就可以有效避免应力集中带来的问题； 2. 支座是另一个需要注意的问题。在有限元分析中，很多时候约束是直接加在混凝土节点上，这样很可能在支座位置产生很大的应力集中，从而使支座附近的混凝土突然破坏，造成求解失败。因此，在实际应用过程中，应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块，避免支座的应力集中； 3. 六面体的SOLID 65 单元一般比四面体的单元计算要稳定且收敛性好，因此，只要条件允许，应该尽量使用六面体单元； 4. 正确选择收敛标准，一般位移控制加载最好用位移的无穷范数控制收敛，而用力控制加载时可以用残余力的二范数控制收敛。在裂缝刚刚出现和接近破坏的阶段，可以适当放松收敛标准，保证计算的连续性； 3. 关于下降段的问题 1）在实际混凝土中都有下降段，但是在计算的时候要特别小心下降段的问题。 2）下降段很容易导致计算不收敛，有时为了计算的收敛要避免设置下降段，采用rush模型。 3）利用最大压应变准则来判断混凝土是否破坏。 4. Solid65单元中的破坏准则 1）采用Willam&Warnke五参数破坏准则 2）需要参数： 单轴抗拉强度，单轴，双轴抗压强度，围压压力，在围压作用下双轴，单轴抗压强度 5. 近来我对混凝土单元进行了一点思考，有一些想法，贴在下面，共同探讨： 1）分析混凝土结构，选择合理的材料特性是建立模型的关键，所以有必要弄清混凝土的材料特性。混凝土是脆性材料,并具有不同的拉伸和压缩特性。典型混凝土的抗拉强度只有抗压强度的8％-15%。 在ANSYS中，对于混凝土单元，材料特性ANSYS要求输入以下数据（为了清楚起见，我将几个系数均译为了中文）：弹性模量、泊松比、张开与闭合滑移面的剪切强度缩减系数、抗拉与抗压强度、极限双轴抗压强度、周围静水应力状态、静水应力状态下单轴与双轴压缩的
《大型结构分析软件的应用及开发》 学习报告 学院：建筑工程学院 专业班级：工程力学141 姓名：付贤凯 指导老师：姚激 学号：201411012111
1．模型介绍 如下图所示的一桁架结构，受一集中力大小为800N的作用，杆件的弹性模量为200GPa，泊松比为0.3。杆件的截面为正方形达长为1m，横截面面积为1m2。现求它的变形图与轴力图。 图1 桁架模型与受力简图（单位：mm） 2．建模与划分网格 利用大型有限元软件ANSYS，采用Link，2Dspar 1的单元进行模拟，通过网格的划分得到如图2所示的有限元模型。 图2 有限元模型
结合有限元模型中的约束条件为左侧在X与Y方向铰支固定，荷载条件为最右侧处施加向下的集中力P=800N。施加约束与荷载后的几何模型如图4所示。 图3 施加荷载与约束的几何模型 3．位移与轴力图 因在Y方向受力，所以主要做Y方向的位移图，又因为杆件在轴线方向有变形，故在X 方向仍有一定的位移。则图5为变形前后的板件形状。图6为模型沿Y方向的位移图，图7为模型沿X方向的位移图，图8为模型的总位移图。 图4 桁架变形前后形状图
图5 Y方向位移图 图6 X方向位移图
图7总位移图 分析所有的位移图可以看出从以看出左端变形最小，为零，右端变形最大。从总位移图可以看出最大的位移在左下点处，大小为0.164×10?5m。从X方向位移图可以看出，左下点处在X方向位移最大为0.36×10?6。从Y方向位移图可以看出最大位移在左下点处为0.164×10?5。都符合实际情况，图9为模型的轴力图。 图8 轴力图
有限元与CAE分析报告 专业： 班级： 学号： 姓名： 指导教师： 2016年 1 月 2 日
简支梁的静力分析 一、问题提出 长3m的工字型梁两端铰接中间1.5m位置处受到6KN的载荷作用，材料弹性模量E=200e9，泊松比0.28，密度7850kg/㎡ 二、建立模型 1.定义单元类型 依次单击Main Menu→Preprocessor→Elementtype→Add/Edit/Delete，出现对话框如图，单击“Add”，出现一个“Library of Element Type”对话框，在“Library of Element Type”左面的列表栏中选择“Structural Beam”，在右面的列表栏中选择3 node 189,单击“OK”。
2设置材料属性 依次单击Main Menu→Preprocessor→MaterialProps>Material Modes，出现“Define Material ModelBehavior”对话框，在“Material Model Available”下面的对话框中，双击打开“Structural→Linear→Elastic→Isotropic”,出现对话框，输入弹性模量EX=2E+011,PRXY=0.28，单击“OK”。 依次单击Main Menu→Preprocessor→MaterialProps>Material Modes，出现“Define Material ModelBehavior”对话框，在“Material Model Available”下面的对话框中，双击打开“Structural→Density”弹出对话框，输入DENS为7850 3.创建几何模型 1)设定梁的截面尺寸
!!!!!~~~~~!!!!!~~~~~!!!!!~~~~~!!!!!~~~~~!!!!!~~~~~!!!!! !!!!!~~~~~~~~~ansys数据处理的相关命令流~~~~~~~~~~~!!!!! !(1)数据输入的相关命令 !利用*TREAD命令读取数据文件并填充TABLE表格 *TREAD, Par, Fname, Ext, --, NSKIP !以下利用*TREAD命令读取1维数据表格 !tdata.txt文本文件含有如下内容 STRAIN STRESS 00 0.0025 0.0046 0.0067 *DIM,Ttxy,table,4,1,,TIME,ACEL *TREAD,Ttxy,tdata,txt,,1 !以下利用*TREAD命令读取2维数据表格 !要特别注意2维数据的行数 !tdata.txt文本文件含有如下内容 TIME X Y Z 0000 0.020.10.20.3 0.040.20.40.6 0.060.30.60.9 !希望输入地震波激励，X、Y、Z三个方向 *DIM,Ttxy,table,3,3,,TIME,ACEL *TREAD,Ttxy,tdata,txt,,1 !以下利用*TREAD命令读取3维数据表格 !tdata.txt文本文件含有如下内容 TEMP X Y Z 0000 0.020.10.20.3 0.040.20.40.6 0.060.30.60.9 5000 0.030.20.30.4 0.050.40.60.8 0.070.60.90.9 !希望读取不同温度下，不同时刻的泊松比 *DIM,Ttxy,table,3,3,2,TIME,NUXP,TEMP *TREAD,Ttxy,tdata,txt,,1 !利用*SREAD命令读取字符文件 *SREAD, StrArray, Fname, Ext, --, nChar, nSkip, nRead 页: 1
第十四章 静电场分析（h方法） 14.1 什么是静电场分析 静电场分析用以确定由电荷分布或外加电势所产生的电场和电场标量位（电压）分布。该分析能加二种形式的载荷：电压和电荷密度。 静电场分析是假定为线性的，电场正比于所加电压。 静电场分析可以使用两种方法：h方法和p方法。本章讨论传统的h方法。下一章讨论p方法。 14.2 h方法静电场分析中所用单元 h方法静电分析使用如下ANSYS单元： 表1. 二维实体单元 单元维数形状或特征自由度 PLANE1212-D四边形，8节点每个节点上的电压 表2. 三维实体单元 单元维数形状或特征自由度 SOLID1223-D砖形（六面体），20节点每个节点上的电压 SOLID1233-D砖形（六面体），20节点每个节点上的电压 表3. 特殊单元 单元维数形状或特征自由度 MATRIX50无(超单元)取决于构成本单元的单元取决于构成本单元的单元类型 INFIN1102-D4或8节点每个节点1个；磁矢量位，温度，或电位 INFIN1113-D六面体，8或20节点AX、AY、AZ磁矢势，温度，电势，或磁标量势 INFIN92-D平面，无界，2节点AZ磁矢势，温度INFIN473-D四边形4节点或三角形3节点AZ磁矢势，温度 14.3 h方法静电场分析的步骤 静电场分析过程由三个主要步骤组成： 1.建模 2.加载和求解 3.观察结果 14.3.1 建模 定义工作名和标题： 命令：/FILNAME，/TITLE GUI：Utility Menu>File>Change Jobname Utility Menu>File>Change Title 如果是GUI方式，设置分析参考框：
学习ansys的一些心得 学习ansys的一些心得（送给初学者和没有盟币的兄弟） 1 做了布尔运算后要重画图形（删除实体）时：需拾取Utility Menu>Plot>Replot 2 标点的输入是在英文状态下，―,‖。 3 线段中点的建立：Modling>Creat>Keypoints>Fill between kps 4 还不会环形阵列。 5 所谓杆系结构指的是长度远远大于其他方向尺寸（10：1）的构件组成的结构，如连续梁，桁架，钢架等。 6 静力学分析的结果包括结构的位移，应变，应力和反作用力等，一般是使用POST1处理（普通后处理器）和查看这些结果。 7 干系结构的静力学分析—平面桁架的建模，用NODE（节点），ELEMENT(元素)创建。复杂体积的建模一般用KPS(关键点)，LINE（Straight line—直线），再生成面，再生成体。 8 如果输入的数据单位是国际单位制单位，则输出的数据单位也是国际制单位。 9 创建正六边形：Creat>Areas>Polygon>Hexagon.指定中心和半径。 10 由面沿线挤出体：Modling>Operate>Extrude>Areas>Along Lines. 11 Ansys中没有Undo命令.需及时保存数据库文件. Def Shape Only:只显示变形图.Def + Undeformed:显示未变形的图.Def + Udef egde:显示未变形的图形的边界. 13 用等高线显示：Plot Results>Contour Plot>Nodal Solu.
14 模态分析用于分析结构的振动特性，即确定结构的固有频率和振型，它也是谐响应分析，瞬态动力学分析以及谱分析等其他动力学分析的基础。 15 Ansys的模态分析是线型分析。任何非线型分析，例如，塑性，接触单元等，即使被定义了也将被忽略。 16 平面桁架：Beam(2D elastic 3) 厚壁圆筒：Solid(8 node 13)>Options(K3—Plane strain) 17 一般材料的弹性模量（EX）：2e11.泊松比(PRXY)：0.3.密度：7800 18 做完静力学分析后，再做模态分析时，要再次求解，同时预应力效果也应该打开（PSTRES,on）.可以在命令行中输入：pstres,on 也可以用菜单路径：Solution>Analysis Type>Analysis Options. 19 弹簧阻尼器单元：Combination-Spring damper 14. 20 接触问题属于状态非线性问题，是一种高度非线性行为，需要较多的计算资源。接触问题有两个基本类型：刚体-柔体的接触，柔体-柔体的接触（许多金属成型的接触问题）。在刚体-柔体的接触问题中，有的接触面与它接触的变形体相比，有较大的刚度而被当做刚体。而柔体-柔体的接触，是一种更普遍的类型，此时两个接触体具有近似的刚度，都为变形体。 21 1 点-点接触：过盈装配问题是用点点接触单元模拟面面接触的典型例子。 2 点-面接触：不必预先知道准确的接触位置，接触面之间也不需要保持一致的网格，并且允许有较大的变形和相对滑动。典型实例：模拟插头插入插座里。 3 面-面接触：刚性面作为目标面，柔性面作为接触面。 22 打开自动时间步长：Solution>Load Step Opts>Time Frequenc>Time And Substps.
第十一章磁宏 11.1 什么是电磁宏 电磁宏是ANSYS宏命令，其主要功能是帮助用户方便地建立分析模型、方便地获取想要观察的分析结果。目前，ANSYS提供了下列宏命令，可用于电磁场分析： ·CMATRIX：计算导体间自有和共有电容系数 ·CURR2D：计算二维导电体内电流 ·EMAGERR：计算在静电或电磁场分析中的相对误差 ·EMF：沿预定路径计算电动力（emf）或电压降 ·FLUXV：计算通过闭合回路的通量 ·FMAGBC：对一个单元组件加力边界条件 ·FMAGSUM：对单元组件进行电磁力求和计算 ·FOR2D：计算一个体上的磁力 ·HFSWEEP：在一个频率范围内对高频电磁波导进行时谐响应分析，并进行 相应的后处理计算 ·HMAGSOLV：定义2-D谐波电磁求解选项并进行谐波求解 ·IMPD：计算同轴电磁设备在一个特定参考面上的阻抗 ·LMATRIX：计算任意一组导体间的电感矩阵 ·MAGSOLV：对静态分析定义磁分析选项并开始求解 ·MMF：沿一条路径计算磁动力 ·PERBC2D：对2—D平面分析施加周期性约束 ·PLF2D：生成等势的等值线图 ·PMGTRAN：对瞬态分析的电磁结果求和 ·POWERH：在导体内计算均方根（RMS）能量损失
·QFACT：根据高频模态分析结果计算高频电磁谐振器件的品质因子 ·RACE：定义一个“跑道形”电流源 ·REFLCOEF：计算同轴电磁设备的电压反射系数、驻波比、和回波损失 ·SENERGY：计算单元中储存的磁能或共能 ·SPARM：计算同轴波导或TE10模式矩形波导两个端口间的反射参数 ·TORQ2D：计算在磁场中物体上的力矩 ·TORQC2D：基于一个圆形环路计算在磁场中物体上的力矩 ·TORQSUM：对2-D平面问题中单元部件上的Maxwell力矩和虚功力矩求和 本章对这些宏有详细描述。在ANSYS命令手册和理论手册对这些宏有更详细的描述。 下面的表格列出了这些电磁宏的使用范畴。
ANSYS使用心得体会 本次结构力学课程设计是学习使用ANSYS软件对框架结构内力进行计算，在未学习该软件前，对于此类问题，通常会采用力矩分配法来进行计算，计算过程繁复，计算量大。导致过程缓慢。 通过对ANSYS软件的学习和了解，知道了它的一些明显的优点。 相对于其他应用型软件而言，ANSYS作为大型权威性的有限元分析软件，对提高解决问题的能力是一个全面的锻炼过程，是一门相当难学的软件，因而，要学好ANSYS，对我们提出了很高的要求，一方面，需要我们有比较扎实的力学理论基础，对ANSYS分析结果能有个比较准确的预测和判断，可以说，理论水平的高低在很大程度上决定了ANSYS使用水平；另一方面，需要我们不断摸索出软件的使用经验不断总结以提高解决问题的效率。 刚开始接触ANSYS时，没有限元，单元，节点，形函数等的基本概念没有清楚的了解话，会感觉还没入门，只是在僵硬的模仿，即使已经了解了，必要先反复看几遍书，加深对有限元单元法及其基本概念的理解。 ANSYS在对结构力学的静力学分析非常方便，用来求解外载荷引起的位移、应力和力。静力分析很适合求解惯性和对结构的影响并不显著的问题。ANSYS 程序中的静力分析不仅可以进行线性分析，而且也可以进行非线性分析，如塑性、膨胀、大变形、大应变及接触分析。 但是学习的过程是充满烦恼和惊喜的，因为总是会碰到许多的新问题，需要较好的耐心去解决这些问题，这是在学习过程中遇到的最大的难题。然而，在解决问题之后，就会有恍然大悟的喜悦，可以说是痛苦和快乐并存的。所以对于初学者，缺乏经验是非常难的。必须保持良好的心态，对于不断出现的ERROR提示要坚定自己的信心，坚信自己可以解决这些问题。所有困难都会迎刃而解。 本次的学习让我认识到了提高建模能力是非常急需加强的一个方面。在做偏向于理论的分析时，可能对建模能力要求不是很高，但对于实际的工程问题，有限元模型的建立可以说是一个最重要的问题，而后面的工作变得相对简单。建模能力的提高，需要掌握好的建模思想和技巧。 ANSYS软件是一款在建模等方面非常实用的软件，本次的学习我其实并没有完全熟练地掌握它的应用，以后还要加强对它的学习，相信在以后的学习和工作中会带来巨大的便利。
ANSYS分析实例详解 姓名：XXX 学号：XXX 专业：XXX 内容：空调支架的有限元分析 本次作业为对一空调支架的有限元分析，其主要内容包括空调支架的建模、有限元分析、强度校核以及结构优化等。下图为空调支架一侧的实物图片： 1、空调支架的特点分析 由于空调支架为一个完全对称结构，空调的重量均匀分部在两侧对称支架上，因此只要对空调支架的一侧进行分析即可达到对整体空调支架的分析，同时也达到了简化空调支架分析的目的。本次作业可以分三部分来完成：一，空调支架一侧的建模；二，利用商业化有限元分析软件对建好的空调支架模型进行有限元分析；三，根据空调支架模型有限元分析的结果对支架进行强度校核以及结构优化。 2、空调支架的建模 空调支架的具体尺寸图如下图所示：
考虑到空调支架模型结构简单，故在此没有利用三维软件建模而是直接在有限元分析软件中进行建模，本次作业采用的有限元分析软件为美国ANSYS公司研制的大型通用有限元分析（FEA）软件ANSYS10.0。建立模型包括设定分析作业名和标题，定义单元类型、定义材料属性、建立三维模型、划分有限元网格。 2.1设定分析作业名和标题 打开ANSYS软件进入ANSYS操作界面，首先从主菜单中选择【Preferences】命令，勾选Structural。然后从实用菜单中选择【Change Jobname】命令，将文件名修改为Ktiao2，从实用菜单中选择【Change Title】命令，将标题修改为Ktiao2。如下图所示： 2.2定义单元类型 在进行有限元分析时，首先应根据分析问题的几何结构、分析类型和所分析的问题精度要求等，选定适合具体分析的单元类型。本文中选用8节点六面体单元Solid185。如下图所示：
ANSYS电磁场分析指南第十六章 发表时间：2007-9-20 作者: 安世亚太来源: e-works 关键字: ANSYS 电磁场分析 CAE教程 第十六章电路分析 16.1 什么是电路分析 电路分析可以计算源电压和源电流在电路中引起的电压和电流分布。分析方法由源的类型来决定：源的类型分析方法 交流（AC）谐波分析 直流（DC）静态分析 随时间变化瞬态分析 要在电磁学分析中用有限元来模拟全部电势，就必须提供足够的灵活性来模拟载流电磁设备。ANSYS 程序对于电路分析有如下性能： ·用经过改进的基于节点的分析方法来模拟电路分析 ·可以将电路与绕线圈和块状导体直接耦合 ·2-D和3-D模型都可以进行耦合分析 ·支持直流、交流和时间瞬态模拟 ANSYS程序中先进的电路耦合模拟功能精确地模拟多种电子设备，： ·螺线管线圈 ·变压器 ·交流机械 16.2 使用CIRCU124单元 ANSYS提供一种通用电路单元CIRCU124对线性电路进行模拟，该单元求解未知的节点电压（在有些情况下为电流）。电路由各种部件组成，如电阻、电感、互感、电容、独立电压源和电流源、受控电压源和电流源等，这些元件都可以用CIRCU124单元来模拟。 注：本章只描述CIRCU124单元的某些最重要的特性，对该单元的详细描述参见《ANSYS单元手册》。 16.2.1 可用CIRCU124单元模拟的电路元件 对CIRCU124单元通过设置KEYOPT（1）来确定该单元模拟的电路元件，如下表所示。例如，把KEYOPT（1）设置为2，就可用CIRCU124来模拟电容。对所有的电路元件，正向电流都是从节点I流向节点J。