电场分析的步骤和ANSYShelp中的例子

1、关于电场分析的步骤和例子电场分析要计算的典型物理量有：电场，电流密度，电荷密度，传导焦耳热纯电场分析：包括，稳态电流传导分析，静电场分析，电路分析静电场的分析基础是泊松方程：主要的未知量（节点自由度）是标量电位（电压）可用于电场分析的单元：表1 传导杆单元单元维数形状或特性自由度LINK683-D单轴，2节点温度和电压表 2 2-D实体单元单元维数形状或特性自由度PLANE672-D四边形，4节点温度和电压PLANE1212-D四边形，8节点电压表 3 3-D实体单元单元维数形状或特性自由度使用注意SOLID53-D六面体,8节点每个节点6个自由度; 可以是位移、温度、电压、磁标量位可用作热电

2、耦合单元或作为电磁耦合场单元SOLID693-D六面体,8节点节点温度、电压可用作热电耦合单元SOLID983-D四面体,10节点每个节点6个自由度; 可以是位移、温度、电压、磁标量位可用作热电耦合单元或作为电磁耦合场单元SOLID1223-D六面体,20节点电压SOLID1233-D四面体,10节点电压表 4 壳单元单元维数形状或特性自由度SHELL1573-D四边形壳，4节点温度和电压表 5 特殊单元单元维数形状或特性自由度MATRIX50无(超单元)根据结构中包括的单元确定根据包含的单元类型决定INFIN1102-D4或8节点每节点一个，可以是磁矢量位、温度、电压INFIN1113-D六

3、面体, 8 或20节点AX, AY, AZ磁矢量位、温度、标量电位或标量磁位表 6 通用电路单元单元维数形状或特性自由度CIRCU124无通用电路单元，最多可6节点每节点三个；可以是电势、电流或电动势降表7 带电压自由度单元的反作用力ElementKEYOPT (1)DOFsMaterial Property Input for VOLT DOFReaction ForcePLANE67N/ATEMP, VOLTRSVX, RSVYElectric Current Flabel=AMPS LINK68N/ATEMP, VOLTRSVXElectric Current Flabel=AMPS S

4、OLID69N/ATEMP, VOLTRSVX, RSVY, RSVZElectric Current Flabel=AMPS SHELL157N/ATEMP, VOLTRSVX, RSVYElectric Current Flabel=AMPS PLANE531VOLT, AZRSVX, RSVYElectric Current Flabel=AMPSSOLID971AX ,AY, AX, VOLTRSVX, RSVY, RSVZElectric Current Flabel=AMPS4AX, AY, AZ, VOLT, CURRSOLID1171AZ,VOLTRSVX, RSVY, RSV

5、ZElectric Current Flabel=AMPS PLANE121N/AVOLTPERX, PERYElectric Charge Flabel=CHRG SOLID122N/AVOLTPERX, PERY, PERZElectric Charge Flabel=CHRG SOLID123N/AVOLTPERX, PERY, PERZElectric Charge Flabel=CHRG SOLID127N/AVOLTPERX, PERY, PERZElectric Charge Flabel=CHRG SOLID128N/AVOLTPERX, PERY, PERZElectric

6、Charge Flabel=CHRG CIRCU124012VOLT, CURR, EMFN/AElectric Current Flabel=AMPS CIRCU1250 or 1VOLTN/AElectric Current Flabel=AMPS TRANS126N/AUX-VOLT, UY-VOLT, UZ-VOLTN/AElectric Current Flabel=AMPS PLANE136VOLT, AZRSVX, RSVYElectric Current Flabel=AMPS 7UX, UY, UZ, VOLTPERX, PERYNegative Electric Charg

7、e Flabel=AMPS SOLID50UX,UY, UZ, TEMP, VOLT, MAGRSVX, RSVY, RSVZElectric Current Flabel=AMPS PERX, PERY, PERZNegative Electric Charge Flabel=AMPS 1TEMP, VOLT, MAGRSVX, RSVY, RSVZElectric Current Flabel=AMPS 3UX, UY, UZ, VOLTPERX, PERY, PERZNegative Electric Charge Flabel=AMPS 9VOLTRSVX, RSVY, RSVZEle

8、ctric Current Flabel=AMPS SOLID980UX,UY, UZ, TEMP, VOLT, MAGRSVX, RSVY, RSVZ Electric Current Flabel=AMPS PERX, PERY, PERZNegative Electric Charge Flabel=AMPS 1TEMP, VOLT, MAGRSVX, RSVY, RSVZElectric Current Flabel=AMPS 3UX, UY, UZ, VOLTPERX, PERY, PERZNegative Electric Charge Flabel=AMPS 9VOLTRSVX,

9、 RSVY, RSVZElectric Current Flabel=AMPS SOLID62N/AUX, UY,UZ, AX, AY, AZ, VOLTRSVX, RSVY, RSVZElectic Current Flabel=AMPS INFIN1101VOLTPERX, PERYElectric Charge Flabel=CHRG INFIN1112VOLTPERX, PERY, PERZElectric Charge Flabel=CHRG 稳态电流传导分析简介稳态电流传导分析可以分析计算直流电流和电压降产生的电流密度和电位分布。可以进行两种加载：电压和电流。稳态电流传导分析认为电

10、压和电流成线性关系，即电流与所加电压成正比。稳态电流传导分析的步骤稳态电流传导分析有三个主要的步骤：1.建立模型2.加载并求解3.观察结1. 建立模型建立模型，定义工作文件名和标题：命令：/FILNAME, /TITLEGUI：Utility Menu>File>Change JobnameUtility Menu>File>Change Title在GUI参数选择框中选择Electric选项。以便能够选择需要的单元。GUI:Main Menu>Preferences>Electromagnetics>Electric然后按照ANSYS建模与分网指南中

11、的描述定义单元类型、定义材料特性并建立几何模型。在电流传导分析中，可以使用下列单元：·LINK68：三维二节点热/电线单元·PLANE67：二维四节点热/电四边形单元·SOLID5：三维八节点结构/热/磁/电六面体单元·SOLID69：三维八节点热/电六面体单元·SOLID98：三维十节点结构/热/磁/电四面体单元·SHELL157：三维四节点热/电壳单元·MATRIX50：三维超单元单元的详细介绍可参看前面单元表。必须只定义一种材料特性：电阻(RSVX)，它可以是和温度有关的。2. 加载并求解此步骤定义分析类型及其选项、给

12、模型加载、定义载荷步选项并求解：进入SOLUTION处理器命令：/SOLUGUI:Main Menu>Solution定义分析类型作下列任何一个操作： 在GUI方式下，选择路径：Main Menu>Solution>New Analysis并选择Static分析。 如果是一个新的分析，执行下列命令：ANTYPE,static,new如果是需要重启动一个前面做过的分析（如施加了另外一种激励），使用命令ANTYPE,STATIC,REST。如果先前分析的结果文件Jobname.EMAT, Jobname.ESAV, 和Jobname.DB还可用，就可以重启动分析。定义分析选项选择

13、方程求解器(系统缺省使用Frontal求解器)。命令：EQSLVGUI:Main Menu>Solution>Analysis Options加载电流电流(AMPS)是经常加在模型边界上的集中节点载荷(AMPS仅仅是一个载荷标志，和单位制无关），正值代表电流流入节点，负值代表流出节点。如果是均匀电流密度分布，应该耦合节点上的VOLT自由度，再将总电流加到某一个节点上去。命令：FGUI: Main Menu>Solution>-Loads-Apply>-Electric-Excitation>Current电压电压(VOLT)是经常加在模型边界上的DOF约束，

14、一个典型的应用是说明导体的一端电压值为零（接地端），另一端为一给定电压。命令：DGUI: Main Menu>Solution>-Loads-Apply>-Electric-Excitation>Charge备份数据用工具条中的SAVE_DB按钮来备份数据库，如果计算机出错，可以方便的恢复需要的模型数据。恢复模型时，用下面的命令：命令：RESUMEGUI:Utility Menu>File>Resume Jobname.db开始求解命令：SOLVEGUI:Main Menu>Solution>Current LS施加其它载荷条件如果希望进行其他加

15、载情况的计算，可以从这里再按照上述步骤操作即可。完成求解命令：FINISHGUI:Main Menu>Finish观看结果ANSYS把结果文件写入Jobename.RST中，数据有：主数据：节点电压(VOLT)导出数据：节点和单元电场(EFX,EFY,EFZ,EFSUM)单元电流密度(JSX,JSY,JSZ,JSSUM)单元焦耳热(JHEAT)节点感生电流进入后处理器：命令：/POST1GUI：Main Menu>General Postproc在POST1中读结果在POST1中后处理数据时，数据库中的模型数据一定要与结果数据相统一，且存在Jobname.RST文件。用RESUME

16、命令读入模型数据，用SET命令读入结果数据。用下列命令把希望的时间点的结果读入数据库：命令：SET,TIME GUI: Utility Menu>List>Results>Load Step Summary如果没有数据和指定的时间点对应，程序自动进行线性插值以得到在指定的时间点处的数据。处理杆单元(LINK68)时，为了得到导出数据，必须使用下列命令读结果到数据库中：命令：ETABLEGUI:Main menu>General Postproc>Element Table>Define Table命令：PLETABGUI:Main menu>Gener

17、al Postproc>Plot Results>Elem TableMain menu>General Postproc>Element Table>Plot Elem Table命令：PRETABGUI:Main menu>General Postproc>List Results>List Elem TableMain menu>General Postproc>Element Table>Elem Table Data等值线显示：命令：PLESOL, PLNSOLGUI:Main menu>General Post

18、proc>Plot Results>Element SolutionMain menu>General Postproc>Plot Results>Nodal Solu矢量（箭头）显示：命令：PLVECTGUI:Main menu>General Postproc>Plot Results>PredifinedGUI:Main menu>General Postproc>Plot Results>User Defined列表显示：命令：PRESOL, PRNSOL, PRRSOLGUI:Main menu>General

19、Postproc>List Results>Element SolutionMain menu>General Postproc>List Results>Nodal SoluMain menu>General Postproc>List Results>Reaction Solu其他后处理，请参见ANSYS基本过程手册。电流传导分析的其他例题VM117 网路中的电流VM170 正方形电流环中的磁场VM173 电线中心线温度静电场分析（h方法）静电场分析用以确定由电荷分布或外加电势所产生的电场和电场标量位（电压）分布。该分析能加二种形式的载荷：电

20、压和电荷密度。静电场分析是假定为线性的，电场正比于所加电压。静电场分析可以使用两种方法：h方法和p方法。h方法静电场分析中所用单元表1. 二维实体单元单元维数形状或特征自由度PLANE1212-D四边形，8节点每个节点上的电压表2. 三维实体单元单元维数形状或特征自由度SOLID1223-D砖形（六面体），20节点每个节点上的电压SOLID1233-D砖形（六面体），20节点每个节点上的电压表3. 特殊单元单元维数形状或特征自由度MATRIX50无(超单元)取决于构成本单元的单元取决于构成本单元的单元类型INFIN1102-D4或8节点每个节点1个；磁矢量位，温度，或电位INFIN1113-D

21、六面体，8或20节点AX、AY、AZ磁矢势，温度，电势，或磁标量势INFIN92-D平面，无界，2节点AZ磁矢势，温度INFIN473-D四边形4节点或三角形3节点AZ磁矢势，温度h方法静电场分析的步骤静电场分析过程由三个主要步骤组成：1.建模2.加载和求解3.观察结果建模定义工作名和标题：命令：/FILNAME，/TITLEGUI：Utility Menu>File>Change JobnameUtility Menu>File>Change Title如果是GUI方式，设置分析参考框：GUI：Main Menu>Preferences>Electroma

22、gnetics：Electric对于静电分析，必须定义材料的介电常数（PERX），它可能与温度有关，可能是各向同性，也可能是各向异性。对于微机电系统（MEMS），最好能更方便地设置单位制，因为一些部件只有几微米大小。详见下面MKS制到µMKSV制电参数换算系数和MKS制到µMSVfA制电参数换算系数表表4 MKS制到µMKSV制电参数换算系数表电参数MKS制量纲乘数µMKSV制量纲电压V(kg)(m)2/(A)(s)31V(kg)(µm)2/(pA)(s)3电流AA1012pApA电荷C(A)(s)1012pC(pA)(s)导电率S/m(A)2

23、(s)3/(kg)(m)3106pS/µm(pA)2(s)3/(kg)(µm)3电阻率m(kg)(m)3/(A)2(s)310-6Tµm(kg)(µm)3/(pA)2(s)3介电常数1F/m(A)2(s)4/(kg)(m)3106pF/µm(pA)2(s)2/(kg)(µm)3能量J(kg)(m)2/(s)21012pJ(kg)(µm)2/(s)2电容F(A)2(s)4/(kg)(m)21012pF(pA)2(s)4/(kg)(µm)2电场V/m(kg)(m)/(s)3(A)10-6V/µm(kg)(&#

24、181;m)/(s)3(pA)通量密度C/(m)2(A)(s)/(m)21pC/(µm)2(pA)(s)/(µm)2自由空间介电常数等于8.0854E-6pF/µm 表5 MKS制到µMSVfA制电参数换算系数表电参数MKS 制量纲乘数µMSVfA制量纲电压V(kg)(m)2/(A)(s)31V(g)(µm)2/(fA)(s)3电流AA1015fAfA电荷C(A)(s)1015fC(fA)(s)导电率S/m(A)2(s)3/(kg)(m)3109fS/µm(fA)2(s)3/(g)(µm)3电阻率m(Kg)(m)3

25、/(A)2(s)310-9-(g)(µm)3/(fA)2(s)3介电常数F/m(A)2(s)4/(kg)(m)3109fF/µm(fA)2(s)2/(g)(µm)3能量J(kg)(m)2/(s)21015fJ(g)(µm)2/(s)2电容F(A)2(s)4/(kg)(m)21015fF(fA)2(s)4/(g)(µm)2电场V/m(kg)(m)/(s)3(A)10-6V/µm(g)(µm)/(s)3(fA)通量密度C/(m)2(A)(s)/(m)2103fC/(µm)2(fA)(s)/(µm)2自由空间介

26、电常数等于8.0854E-3fF/µm加载荷和求解本步定义分析类型和选项、给模型加载、定义载荷步选项和开始求解。进入求解处理器命令：/SOLUGUI：Main Menu>Solution定义分析类型选择下列方式之一：·GUI：选菜单路径Main Menu>Solution>New Analysis并选择静态分析·命令：ANTYPE，STATIC，NEW·如果你要重新开始一个以前做过的分析（例如，分析附加载荷步），执行命令ANTYPE，STATIC，REST。重启动分析的前提条件是：预先完成了一个静电分析，且该预分析的Jobname. E

27、MAT，Jobname. ESAV和Jobname.DB文件都存在。定义分析选项可以选择波前求解器（缺省）、预条件共轭梯度求解器（PCG）、雅可比共轭梯度求解器（JCG）和不完全乔列斯基共轭梯度求解器（ICCG）之一进行求解：命令：EQSLVGUI：Main Menu>Solution>Analysis Options如果选择JCG求解器或者PCG求解器，还可以定义一个求解器误差值，缺省为1.0-8。加载静电分析中的典型载荷类型有：14.3.2.4.1 电压（VOLT）该载荷是自由度约束，用以定义在模型边界上的已知电压：命令：DGUI：Main Menu>Solution&g

28、t;Loads>-Loads-Apply>-Electric-Boundary> -Voltage-电荷密度（CHRG）命令：FGUI：Main Menu>Solution>Loads>-Loads-Apply>-Electric-Excitation>-Charge-On Nodes面电荷密度（CHRGS）命令：SFGUI：Main Menu>Solution>Loads>-Loads-Apply>-Electric-Excitation-Surf Chrg Den-Maxwell 力标志（MXWF）这并不是真实载荷，只

29、是表示在该表面将计算静电力分布，MXWF只是一个标志。通常，MXWF定义在靠近“空气-电介质”交界面的空气单元面上，ANSYS使用Maxwell应力张量法计算力并存储在空气单元中，在通用后处理器中可以进行处理。命令：FMAGBCGUI：Main Menu>Solution>-Loads-Apply>-Electric-Flag>-Maxwell Surf-option无限面标志（INF）这并不是真实载荷，只是表示无限单元的存在，INF仅仅是一个标志。命令：SFGUI：Main Menu>Solution>-Loads-Apply>-Electric-F

30、lag>-Infinite Surf-option体电荷密度（CHRGD）命令：BF，BFEGUI：Main Menu>Solution>-Loads-Apply>-Electric-Excitation>-Charge Density-option另外，还可以用命令BFL、BFL、BFV等命令分别把体电荷密度加到实体模型的线、面和体上。定义载荷步选项对于静电分析，可以用其它命令将载荷加到电流传导分析模型中，也能控制输出选项和载荷步选项，详细信息可参见第16章“分析选项和求解方法”保存数据库备份使用ANSYS工具条的SAVE_DB按钮来保存一个数据库备份。在需要的

31、时候可以恢复模型数据：命令：RESUMEGUI：Utility Menu>File>Resume Jobname.db开始求解命令：SOLVEGUI：Main Menu>Solution>Current LS结束求解命令：FINISHGUI：Main Menu>Finish观察结果ANSYS和ANSYS/Emag程序把静电分析结果写到结果文件Jobname.RST中，结果中包括如下数据：主数据：节点电压（VOLT）导出数据：·节点和单元电场（EFX，EFY，EFZ，EFSUM）·节点电通量密度（DX，DY，DZ，DSUM）·节点静电力

32、（FMAG：分量X，Y，Z，SUM）·节点感生电流段（CSGX，CSGY，CSGZ）通常在POST1通用后处理器中观察分析结果：命令：/POST1GUI：Main Menu>General Postproc对于整个后处理功能的完整描述，见ANSYS基本分析过程指南。将所需结果读入数据库：命令：SET,TIMEGUI：Utility Menu>List>Results>Load Step Summary如果所定义的时间值处并没有计算好的结果，ANSYS将在该时刻进行线性插值计算。对于线单元（LINK68），只能用以下方式得到导出结果：命令：ETABLEGUI：M

33、ain Menu>General Postproc>Element Table>Define Table命令：PLETABGUI：Main Menu>General Postproc>Plot Results>Elem TableMain Menu>General Postproc>Element Table>Plot Elem Table命令：PRETABGUI：Main Menu>General Postproc>List Results>List Elem TableMain Menu>General Post

34、proc>Element Table>Elem Table Data绘制等值线图：命令：PLESOL，PLNSOLGUI：Main Menu>General Postproc>Plot Results>Element SolutionMain Menu>General Postproc>Plot Results>Nodal Solu绘制矢量图：命令：PLVECTGUI：Main Menu>General Postproc>Plot Results>PredefinedMain Menu>General Postproc&g

35、t;Plot Results>User Defined以表格的方式显示数据：命令：PRESOL，PRNSOL，PRRSOLGUI：Main Menu>General Postproc>List Results>Element SolutionMain Menu>General Postproc>List Results>Nodal SolutionMain Menu>General Postproc>List Results>Reaction SoluPOST1执行许多其他后处理功能，包括按路径和载荷条件的组合绘制结果图。更详细信息见

36、ANSYS基本分析过程手册。多导体系统提取电容（很像翻译的HELP， 故我把它全部copy过来）14.4静电场分析求解的一个主要参数就是电容。在多导体系统中，包括求解自电容和互电容，以便在电路模拟中能定义等效集总电容。CMATRIX宏命令能求得多导体系统自电容和互电容。详见ANSYS理论手册5.10节。14.4.1 对地电容和集总电容有限元仿真计算，可以提取带（对地）电压降导体由于电荷堆积形成的“对地”电容矩阵。下面叙述一个三导体系统（一个导体为地）。方程式中Q1和Q2为电极1和2上的电荷，U1和U2分别为电压降。Q1= (Cg)11(U1)+(Cg)12(U2)Q2= (Cg)12(U1)+

37、(Cg)22(U2)式中Cg称作为“对地电容”矩阵。这些对地电容并不表示集总电容（常用于电路分析），因为它们不涉及到二个导体之间的电容。使用CMATRIX宏命令能把对地电容矩阵变换成集总电容矩阵，以便用于电路仿真。图2描述了三导体系统的等效集总电容。下面二个方程描述了感应电荷与电压降之间形成的集总电容：Q1=(C1)11(U1)+(C1)12(U1U2)Q2=(C1)12(U1U2)+(C1)22(U2)式中C1称为集总电容的电容矩阵。14.4.2 步骤CMATRIX宏命令将进行多元模拟，可求得对地电容矩阵和集总电容矩阵值。为了便于CMATRIX宏命令使用，必须把导体节点组成节点部件，而且不要

38、加任何载荷到模型上（电压、电荷、电荷密度等等）。导体节点的部件名必须包括同样的前缀名，后缀为数字，数字按照到系统中所含导体数目进行编号。最高编号必须为地导体（零电压）。应用CMATRIX宏命令步骤如下：1.建模和分网格。导体假定为完全导电体，故导电体区域内部不需要进行网格划分，只需对周围的电介质区和空气区进行网格划分，节点部件用导体表面的节点表示。2.选择每个导体面上的节点，组成节点部件。命令：CMGUI:Utility Menu >Select >Comp/Assembly >Create Component导体节点的部件名必须包括同样的前缀名，后缀为数字，数字按照到系统中

39、所含导体数目进行编号。例如图2中，用前缀“Cond”为三导体系统中的节点部件命名，分别命名为为“Condl”、“Cond2”和“Cond3”，最后一个部件“Cond3”应该为表示地的节点集。3.用下列方法之一，进入求解过程：命令：SoluGUI：Main Menu>Solution4.选择方程求解器（建议用JCG）：命令：EQSLVGUI：Main Menu>Solution>Analysis Options5.执行CMATRIX宏：命令：CMATRIXGUI：Main Menu >Solution >Electromagnet >Capac MatrixC

40、MATRIX宏要求下列输入：·对称系数（SYMFAC）：如果模型不对称，对称系数为1（缺省）。如果你利用对称只建一部分模型，乘以对称系数得到正确电容值。·节点部件前缀名（Condname）。定义导体节点部件名。上例中，前缀名为“Cond”。宏命令要求字符串前缀名用单引号。因此，本例输入为Cond，在GUI菜单中，程序会自动处理单引号。·导体系统中总共的节点部件数（NUMCON），上例中，导体节点部件总数为“3”。·地基准选项（GRNDKEY）。如果模型不包含开放边界，那么最高节点部件号表示“地”。在这种情况下，不需特殊处理，直接将“地”作为基准设置为零（

41、缺省状态值）。如果模型中包含开放边界（使用远场单元或Trefttz区域），而模型中无限远处又不能作为导体，那么可以将“地”选项设置为零（缺省）。在某些情况下，必须把远场看作导体“地”（例如，在空气中单个带电荷球体，为了保持电荷平衡，要求无限远处作为“地”）。用INFIN111单元或Trefftz区域表示远场地时，把“地”选项设置为“1”·输入贮存电容值矩阵的文件名（Capname）。宏命令贮存所计算的三维数组对地电容和集总电容矩阵值。其中“i”和“j”列代表导体编号，“k”列表示对地（k=1）或集总（k=2）项。缺省名为CMATRIX。例如，CMATRIX（i,j,1）为对地项，CM

42、ATRIX（i,j,2）为集总项。宏命令也建立包含矩阵的文本文件，其扩展名为.TXT。注意：在使用CMATRIX命令前，不要施加非均匀加载。以下操作会造成非均匀加载：·在节点或者实体模型上施加非0自由度值的命令(D, DA, 等)·在节点、单元或者实体模型定义非0值的命令(F, BF, BFE, BFA, 等)·带非0项的CE命令CMATRIX执行一系列求解，计算二个导体之间自电容和互电容，求解结果贮存在结果文件中，可以便于后处理器中使用。执行后，给出一个信息表。如果远场单元（INFIN110和INFIN111）共享一个导体边界（例如地平面），可以把地面和无限远边

43、界作为一个导体（只需要把地平面节点组成一个节点部件）。下图图3描述了具有合理的NUMCOND和GRNDKEY选项设置值的各种开放和闭合区域模型。后面有例题详细介绍如何利用CMATRIX做电容计算。14.5 开放边界的Trefftz方法模拟开放区域的一种方法是利用远场单元(INFIN110和INFIN111)，另一种方法为混合有限元Trefftz方法（称作Trefftz方法）。Trefftz方法以边界元方法的创立者名字命名。 Trefftz方法使用与有限元类似的正定刚度矩阵高效处理开放区域的边界问题。它可处理大纵横比的复杂面几何体，它很易生成Trefftz完整函数系统。对于处理静电问题中的开放边

44、界条件是一种易用而精确的方法。Trefftz方法的理论分析参见ANSYS理论手册。本手册有“用Trefftz方法进行静电场分析”的例题。14.5.1概述使用Trefftz方法需要建立一个Trefftz区域，Trefftz区域由下列部分组成：·在有限元区域内的一个Trefftz源节点部件，但与有限元模型无关；·带有标记的有限元区域的外表面；·由Trefftz源节点部件和带有标记的有限元外表面共同创建的子结构矩阵；·由子结构定义的超单元；·连同子结构产生的一组约束方程；与远场单元法相比，Trefftz方法有许多优点，也有一些缺点。Trefftz方法

45、有如下正面特征：·本方法形成对称矩阵；·处理开放边界时，不存在理论上的限制；·不存在奇异积分；·未知数最少（20100个未知量就可得到可靠结果）；·可用于大纵横比边界；·允许灵活的生成格林（Greens）函数；·利用Trefftz区域，可以在两个无关联的有限元区之间建立联系；Tefftz方法与远场单元比较有如下优点；·通常具有更高的精确度；·远场区不要求建模和划分单元；·可用于大纵横比有限元区域，并且具有很好精度；·远场单元区不必按一般有限元要求的那样，把有限元区扩展到超出装置模型区很

46、多；Trefftz方法与远场单元比较有如下缺点：·只能用于全对称模型；·只对三维分析有效；·模型外表面单元只能是四面体单元；·要求定义有限元区内Trefftz源节点部件，并生成子结构和约束方程（当然，这一过程是程序自动完成）。Trefftz方法有如下限制：·Trefftz节点最大数为1000；·最高容许的节点号为1,000,000；·最高容许的外表面节点数为100,000；·外表面容许的最大单元面数（小平面）为100,000；Trefftz方法假设无限远处是0电位。因此，在处理具有不同电位的多电极系统时，使用本方法

47、要注意建立不同的节点部件。当然，对于使用CMATRIX命令宏来提取电容，程序已经完全考虑，已经把无限远处设成了0电位或者接近0电位。14.5.2 步骤在3-D静电分析中建立一个Trefftz区域，定义Trefftz区域按下列过程进行：1）建立一个静电区域的有限元模型（包括导体、介质和四周空气）。对有限模型加上全部必需的边界条件（电压、电荷、电荷密度等）2）对有限元区域的外表面加上标志，作一个无限面来处理。加无限面标志（INF Label）,使用如下方法：命令：SF，SFA，SFEGUI：Main Menu>Preprocessor>Loads>-Loads- Apply>

48、;-Electric- Flag>-Infinite Surf-On Nodes Main Menu>Preprocessor>Loads>-Loads- Apply>-Electric- Flag>-Infinite Surf-On Areas Main Menu>Preprocessor>Trefftz-Domain>Infinite Surf-On Areas3）建立Trefftz源节点，源节点作为Trefftz 区域的未知量。这些未知量表示Trefftz方法的源电荷，用CURR自由度计算且储存Trefftz节点上的这些源电荷。如图4

49、“定义 Trefftz的区域”中步骤3所示，应在模型装置与有限元区外表面之间设置Trefftz源节点。Trefftz源节点离模型装置的距离应该小于到有限元模型外表面的距离，这样Trefftz 方法计算所得的结果会更精确。Trefftz源节点离有限元模型外表面表面越远，得到的结果越精确。例如，X方向上，Trefftz节点正好包围模型装置（b/c>1）,有限元外边界设置到较远距离处(a/b>2)。对Y和Z方向应用大致相同的规则。若Trefftz源节点不接近于装置或在有限元区域表面上，会导致一个近似奇异解而产生不正确的结果。利用定义一个简单实模型体(如六面体、球、园柱体或它们的布尔运算组

50、合体等)，很容易地建立包围模型装置并在有限元区域内的Trefftz节点。但是它应该在有限元外表面的内部，如图4所示。一旦定义了简单模型体，可以采用下列方法之一把简单实模型划分网格并建立Trefftz节点：命令：TZAMESHGUI：Main Menu>Preprocessor>Trefftz Domain>Mesh TZ Geometry用TZAMESH命令对体表面进行网格划分，然后删除非求解单元，只留下Trefftz节点。它把 Trefftz节点组成命名为TZ-NOD的节点部件，以备在Trefftz子结构生成中调用。Trefftz方法只要求很少的源节点。缺省时，TZAMES

51、H命令把简单实模型体各边分成二段。对大纵横比几何体，可按规定的长度划分实体。这二种选项在TZAMESH命令中都有效。它会提供很多Trefftz节点，但是并不是节点越多精度越高。精度也受外表面单元数和Trefftz 源项近似的影响。一般例题将不超过20到100Trefftz节点。利用下列方法，可删除Trefftz节点：命令：CMSEL,TZ_NOD NDELE,ALLCMDELE,TZ_NODGUI：Main Menu>Preprocessor>Trefftz Domain>Delete TZ Nodes4）建立Trefftz 子结构、超单元、和约束方程。Trefftz方法使用

52、有限元模型的外表面和Trefftz节点建立子结构矩阵。用MATRIX50超单元将该矩阵组合到模型中。另外，需要一组约束方程来完善Trefftz 区域。利用TZEGEN宏命令，可自动完成建立子结构、用超单元组合到模型、定义约束方程等过程。建立子结构并使其以超单元的方式组合到模型中，用下列方式：命令：TZEGENGUI：Main Menu>Preprocessor>Trefftz Domain>-Superelement-Generate TZTZEGEN命令也自动定约束方程。一旦建立了Trefftz区域，就可利用标准求解步骤来解题。如果分网面上的单元发生了变动或要建立一个新的T

53、refftz 区域，则已定义的Trefftz区域应被删除。在求解模型内只能同时存在一个 Trefftz区域。采用下列方法可删除Trefftz 超单元、相应的约束方程和全部Trefftz 文件：命令：TZDELEGUI：Main Menu>Preprocessor>Trefftz Domain>-Superelement-Delete TZTZDELE命令删除在生成超单元过程中产生的全部Trefftz文件，包括如下文件：·Jobname.TZN Trefftz源节点· Jobname.TZE 在有限元边界上的Trefftz 表面·Jobname.

54、TZX 在有限元边界上表面节点·Jobname.TZM Trefftz 材料文件详见本手册例题“用Trefftz方法进行静电分析(命令方法)”CALL，下面竟然是例子，HELP中的，怎么办，也copy吧14.6 用h方法进行静电场分析的实例（GUI方式）14.6.1 问题描述本节描述如何做一个屏蔽微带传输线的静电分析，该传输线是由基片、微带和屏蔽组成。微带电势为V1，屏蔽的电势为V0，确定传输线的电容。该算例的描述见下图。材料和几何参数14.6.2 分析方法与建模提示通过能量和电位差的关系可以求得电容：We = 1/2C(V1-V0)2，We是静电场能量，C为电容。在后处理器中对所有

55、单元能量求和可以获得静电场的能量。后处理器中还可以画等位线和电场矢量图等。14.6.3 目标结果目标电容, pF/m 178.1步骤 1: 开始1.进入ANSYS程序. 2.选择菜单路径Utility Menu>File>Change Title. 3.输入"Microstrip transmission line analysis."4.点击 OK.5.选择Main Menu>Preferences. 6.点击Magnetic-Nodal和Electric.7.点击 OK.步骤 2: 定义参数1.选择Utility Menu>Parameters&

56、gt;Scalar Parameters.2.输入下列参数，若发生输入错误，重新输入即可V1 = 1.5V0 = 0.53.点击 Close步骤 3: 定义单元类型1.选择Main Menu>Preprocessor>Element Type>Add/Edit/Delete. 2.点击 Add.3.点击高亮度的"Electrostatic"和"2D Quad 121."4.点击 OK.5.点击 Close.步骤 4: 定义材料属性1.选择Main Menu>Preprocessor>Material Props>Mat

57、erial Models.2.在材料窗口,依次双击以下选项: Electromagnetics, Relative Permittivity, Constant3.MURX (Relative permeability)输入 1,点击OK.在定义材料的窗口的左边区域显示的材料号为1.4.选择菜单路径 Edit>Copy. 点击OK。把材料1拷贝到材料2.5.在材料框中，双击2号材料和Permittivity (constant). 6在PERX区域输入10, 点击OK.7选择菜单路径 Material>Exit 8点击 SAVE_DB on the ANSYS Toolbar.步骤

58、 5: 建立几何模型和压缩编号1.选择Main Menu>Preprocessor>-Modeling-Create>-Areas-Rectangle> By Dimensions.2.输入下列值(用TAB键，在输入区域间切换)X1 域0X2域.5Y1域0Y2域13.点击 Apply. 4.创建第2个矩形，输入下列值：X1域.5X2域5Y1域0Y2域15.点击 Apply. 6.创建第3个矩形，输入下列值：X1域0X2域.5Y1域17.点击 Apply.8.创建第4个矩形，输入下列值：X1域.5X2域5Y1域1Y2域109.点击 OK.10.粘接所有面, 选择Main Menu>Preprocessor>-Modeling-Operate>-Booleans-Glue>Areas.11.点击 Pick All.12.选择Main Menu>Preprocessor>Numbering Ctrls>Compress Numbers