ANSYS电磁场分析指南 第八章 3

第八章3-D瞬态磁场分析（棱边单元法）8.1用棱边元方法进行3-D瞬态分析用棱边元方法进行3-D瞬态分析的步骤与下一章的用MVP方法进行瞬态分析的步骤大体一样，只不过它们使用不同的自由度和运算法则，棱边元方法使用Frontal、sparse、JCG和ICCG求解器。参见《ANSYS基本过程手册》中求解器的详细介绍。ANSYS支持3D静态、时谐、瞬态棱边单元分析。前面两章讨论的是静态和时谐分析，本章讨论瞬态分析。8.23-D瞬态磁场分析（棱边元方法）的步骤在GUI菜单过滤中选定Magnetic-Edge项定义任务名和题目命令：/FILNAME和/TITLEGUI:UtilityMenu>File>ChangeJobnameUtilityMenu>File>ChangeTitle进入ANSYS前处理器命令：/PREP7GUI:MainMenu>Preprocessor选择SOLID117单元命令：ET,solid117GUI:MainMenu>Preprocessor>ElementType>Add/Edit/Delete定义SOLID117单元选项对导电区用AZ-VOLT自由度，对不导电区用AZ自由度。命令：KEYOPTGUI:MainMenu>Preprocessor>ElementType>Add/Edit/Delete与第七章“谐波分析”描述的类似6．定义材料特性对涡流区必须说明电阻值RSVX，材料定义的过程详见第二章。7．建立模型用MainMenu>Preprocessor>-Modeling-界面。赋予特性命令：VATTGUI:Mainmenu>Preprocessor>-Attributes-Define划分网格（用Mapped网格）命令：VMESHGUI:MainMenu>Preprocessor>-Meshing-Mesh>-Volumes-Mapped进入求解器命令：/SOLUGUI:MainMenu>Solution11．给模型边界加磁力线平行和磁力线垂直边界条件命令：DAGUI:MainMenu>Solution>-Loads-Apply>-Magnetic-Boundary用AZ=0来模拟磁力线平行边界条件，磁力线垂直边界条件是自然边界条件无需说明。给绞线圈加电流密度载荷命令：BFE,jsGUI:MainMenu>Solution>-Loads-Apply>-Magnetic-Excitation除了加电流密度载荷外，还可以给一个块导体加总电流：命令：F,,ampsGUI:MainMenu>Solution>-Loads-Apply>-Electric-Excitation>-ImpressedCurrent-注意：在加总电流之前需耦合节点VOLT自由度。选择瞬态分析类型命令：ANTYPE,transient,newGUI:MainMenu>Solution>NewAnalysis>Transient注意：若要重启动一个先前做过的分析（如施加了另外一种载荷），可使用命令ANTYPE,TRANSIENT,REST。重启动分析要求先前分析的几个相关文件Jobname.EMAT,Jobname.ESAV,和Jobname.DB还可用。定义分析选项定义求解方法和求解器。选择求解方法：命令：TRNOPTGUI：Mainmenu>Solution>AnalysisOptions瞬态磁场分析要求选择全波方法（full）棱边元方法可以选用FRONT（缺省值），sparse,JCG,或者ICCG求解器。用下列命令选择求解器：命令：EQSLVGUI:Mainmenu>Solution>AnalysisOptions选择载荷步选项。如下对载荷步选项作了阐述：时间选项该选项说明载荷步终结的时间。命令：TIMEGUI:MainMenu>Preprocessor>Loads>-LoadStepOpts-Time/Frequenc>Time-TimeStepMainMenu>Solution>-LoadStepOpts-Time/Frequenc>Time-TimeStep子步数和时间步长积分时间步长是瞬态分析时间积分方程中的时间增量，可用DELTIM命令或其等效菜单路径直接定义，或用NSUBST或其等效菜单路径间接定义。时间步步长决定了求解精度，时间步长越小，精度就越高。当载荷出现较大的阶跃时，紧跟其后的第一个时间积分步长是尤为关键的。通过减小时间步步长，可以减少求解大阶跃变化(如温度过热加载)时的误差。注意：时间步步长也不能过小，尤其是在建立初始化条件时。太小的数值会使ANSYS在计算时产生数值误差。比如，如果使用小于1E-10的时间步步长会产生数值误差。如选择阶跃(Stepped)加载模式，在第一个子步上程序就加上全部载荷并一直保持常数，如选择斜坡(Ramped)加载模式，程序在每个子步上增加载荷值。命令：NSUBST,DELTIMGUI:MainMenu>Preprocessor>Loads>-LoadStepOpts-Time/Frequenc>TimeandSubstpsMainMenu>Solution>-LoadStepOpts-Time/Frequenc>TimeandSubstpsMainMenu>Solution>-LoadStepOpts-Tiem/Frequenc>Time-TimeStep自动计算时间步长在瞬态分析中也叫时间步长最优化，它使程序自动调整两个子步骤间的载荷增量，并在求解过程中根据模型的响应情况增加或减小时间步长。在大多数情况下，你需要打开这个选项，此外为了更好地控制时间步长的变化幅度，还要输入积分时间步长的上限和下限。命令：AUTOTSGUI:MainMenu>Preprocessor>Loads>-LoadStepOpts-Time/frequenc>TimeMainMenu>Solution>-LoadStepOpts-Time/Frequenc>TimeandSubstpsMainMenu>Solution>-LoadStepOpts-Tiem/Frequenc>Time-TimeStepNewton-Raphson选项这些选项定义在非线性求解过程中，切向刚度矩阵的更新频率，选项为：•程序自动选择(Program-chosen)(系统缺省)・全方法(Full)•修正法(Modified)•初刚度法(Initial-stiffness)在非线性分析中，推荐使用Full选项，其adaptivedescent选项可有助于瞬态分析收敛。命令：NROPTGUI:MainMenu>solution>AnalysisOptions平衡迭代数该选项使得在每一个子步都能得到一个收敛解，缺省为25次迭代，但是根据所处理问题的非线性程度，可适当增加这个数值。对线性瞬态分析，只需1次迭代。命令：NEQITGUI:MainMenu>Preprocessor>Loads-LoadStepOptions-Nonlinear>EquilibriumIterMainMenu>Solution-LoadStepOptions-Nonlinear>EquationIter收敛容差只要运算满足所说明的收敛判据，程序就认为它收敛，收敛判据可以基于磁势(A),也可以是磁电流段(CSG),或二者都有。在实际定义时，需要说明一个典型值(VALUE)和收敛容差(TOLER)，程序将VALUE*TOLER的值视为收敛判据。例如，如说明磁电流段的典型值为5000，容差为0.001，那么收敛判据则为5.0。ANSYS公司推荐VALUE值由缺省确定，TOLER的值为1.0e-3。命令：CNVTOLGUI:MainMenu>Preprocessor>Loads>Nonlinear>convergenceCritMainMenu>Solution>Nonlinear>ConvergenceCrit终止不收敛解当程序在指定的平衡迭代数内无法收敛，则程序根据你所指定的终止判据或终止求解、或直接进行第二个载荷步的求解。控制打印输出该选项控制你将哪些数据输出到打印输出文件中去。命令：OUTPRGUI:MainMenu>Preprocessor>Loads>OutputCtrls>SoluPrintoutMainMenu>Solution>OutputCtrls>SoluPrintout控制数据库和结果文件输出该选项控制你将哪些数据输出到结果文件中去。根据缺省，程序将每个载荷步中的最后一个子步的数据写入到结果文件中去。命令：OUTRESGUI:MainMenu>Preprocessor>OutputCtrls>DB/ResultFileMainMenu>Solution>OutputCtrls>DB/ResultFile备份数据用工具条中的SAVE_DB按钮来备份数据库，如果计算机出错，可以方便的恢复需要的模型数据。恢复模型时，用下面的命令：命令：RESUMEGUI:UtilityMenu>File>ResumeJobname.db开始求解如果对所有加载情况一并求解，用下列方式：命令:LSSOLVEGUI:MainMenu>Solution>FromLSFiles如果只有一个载荷步，也可以用下列方式：命令:SOLVEGUI:MainMenu>Solution>CurrentLS当使用棱边元方程时，在缺省情况下，ANSYS程序先估算待分析区域所有单元和节点。估算时，把不需要的自由度值设置为零，使计算更快进行：命令：GAUGEGUI:MainMenu>Solution>LoadStepOpts-magnetics>-OptionsOnly-Gauging使用棱边元方法做电磁分析必须要求估算，因此，在大多数情况下，不要关闭自动估算。如果是专家级用户，可以用GAUGE,OFF关闭估算，并用自己的经验程序进行相应估算。估算后ANSYS程序删除进行估算时设置的约束，因此，估算过程对用户是透明的。16．离开求解器命令：FINISHGUI:MainMenu>Finish17.后处理，后面详细阐述8.3观察结果ANSYS/Emag程序将棱边元方法瞬态磁场分析的数据结果写入到Jobname.RMG文件中，结果数据包括：主数据：节点自由度（AZ,VOLT）导出数据：•节点磁通量密度（BX,BY,BZ,BSUM）•节点磁场强度（HX,HY,HZ,HSUM）•节点洛仑兹磁力（FMAG:X,Y,Z分量和SUM）•单元总电流密度（JTX,JTY,JTZ）•单位体积生成的焦耳热（JHEAT）•单元磁能（仅对线性材料才有效）在后处理磁场分析的解时，在POST1通用后处理器中可观察整个模型在给定时间点处的解；在POST26时间历程后处理器中可观察模型中某一指定点在整个瞬态分析中的历程解，如下选择后处理器。命令：/POST1,/POST26GUI:MainMenu>GeneralPostprocMainMenu>TimeHistPostpro后处理命令的完整描述，参见《ANSYS基本过程手册》。8.3.1在P0ST26时间历程后处理器中观察结果数据可参看本指南第6章的详细阐述。在后处理数据时，数据库中的模型数据一定要与结果数据相统一，结果文件必须可用。如果模型不在数据库中，用下列恢复操作，再用SET命令读入想了解的结果集数据。命令：RESUMEGUI:UtilityMenu>File>ResumeJobname.dbP0ST26后处理中研究的是各种变量和时间的关系，每个变量有一个确定的参考号，时间变量为参考号1，其他变量需定义。定义主数据变量：命令：NS0LGUI：MainMenu>TimeHistPostpro>DefineVariables定义导出数据变量：命令：ES0LGUI：MainMenu>TimeHistPostpro>DefineVariables定义感生数据变量：命令：RF0RCEGUI：MainMenu>TimeHistPostpro>DefineVariables定义好一个变量后，可将它根据时间或其他任何变量进行图形显示：命令：PLVARGUI：MainMenu>TimeHistPostpro>GraphVariables对定义好的一个变量列表：命令：PRVARGUI：MainMenu>TimeHistPostpro>ListVariables仅列出最大变量值：命令：EXTREMGUI：MainMenu>TimeHistPostpro>ListExtremes可用PMGTRAN命令计算某个单元组元(Component)总的磁力、能量损耗、能量和电流值：命令：PMGTRANGUI：MainMenu>TimeHistPostpro>Elec&Mag>Magnetics在POST26后处理器中，还可以对各种变量进行数学运算，也可将变量的值赋给某数组参数中。参见《ANSYS基本过程手册》8.3.2在POST1时间历程后处理器中观察结果数据在POST1中后处理数据