14.1 什么是静电场分析
静电场分析用以确定由电荷分布或外加电势所产生的电场和电场标量位(电压)分布。该分析能加二种形式的载荷:电压和电荷密度。 静电场分析是假定为线性的,电场正比于所加电压。
静电场分析可以使用两种方法:h方法和p方法。本章讨论传统的h方法。下一章讨论p方法。
h方法静电分析使用如下ANSYS单元:
表1. 二维实体单元
单元 | 维数 | 形状或特征 | 自由度 |
PLANE121 | 2-D | 四边形,8节点 | 每个节点上的电压 |
| | | |
表2. 三维实体单元
单元 | 维数 | 形状或特征 | 自由度 |
SOLID122 | 3-D | 砖形(六面体),20节点 | 每个节点上的电压 |
SOLID123 | 3-D福建金融职业技术学院 | 砖形(六面体),20节点 | 每个节点上的电压 |
| | | |
表3. 特殊单元
单元 | 维数 | 形状或特征 | 自由度 |
MATRIX50 | 无(超单元) | 取决于构成本单元的单元 | 取决于构成本单元的单元类型 |
INFIN110 | 2-D | 4或8节点 | 每个节点1个;磁矢量位,温度,或电位 |
INFIN111 | 3-D | 六面体,8或20节点 | AX、AY、AZ磁矢势,温度,电势,或磁标量势 |
INFIN9 | 2-D | 平面,无界,2节点 | AZ磁矢势,温度 |
INFIN47 | 3-D | 四边形4节点或三角形3节点 | AZ磁矢势,温度 |
| | | |
14.3h方法静电场分析的步骤
静电场分析过程由三个主要步骤组成:
1.建模
3.观察结果
14.3.1 建模
命令:/FILNAME,/TITLE
GUI:Utility Menu>File>Change Jobname
Utility Menu>File>Change Title
如果是GUI方式,设置分析参考框:
GUI:Main Menu>Preferences>Electromagnetics:航空工业中心医院Electric
设置为Electric,以确保电场分析所需的单元能显示出来。之后就可以使用ANSYS前处理器来建立模型,其过程与其它分析类似,详见《ANSYS建模和分网指南》。
对于静电分析,必须定义材料的介电常数(PERX),它可能与温度有关,可能是各向同性,也可能是各向异性。
对于微机电系统(MEMS),最好能更方便地设置单位制,因为一些部件只有几微米大小。详见下面MKS制到µMKSV制电参数换算系数和MKS制到µMSVfA制电参数换算系数表
表4. MKS制到µMKSV制电参数换算系数表
电参数 | MKS制 | 量纲 | 乘数 | µMKSV制 | 量纲 |
电压 | V | (kg)(m)2/(A)(s)3 | 1 | V | (kg)(µm)2/(pA)(s)3 |
电流 | A | A | 1012 | pA | pA |
电荷 | C | (A)(s) | 1012 | pC | (pA)(s) |
导电率 | S/m | (A)2(s)3/(kg)(m)3 | 106 | pS/µm | (pA)2(s)3/(kg)(µm)3 |
电阻率 | Ωm | (kg)(m)3/(A)2(s)3 | 10-6 | TΩµm | (kg)(µm)3/(pA)2(s)3 |
介电常数1 | F/m | (A)2(s)4/(kg)(m)3 | 106 | pF/µm | (pA)2(s)2/(kg)(µm)3 |
能量 | J | (kg)(m)2/(s)2 | 1012 | pJ | (kg)(µm)2/(s)2 |
电容 | F | (A)2(s)4/(kg)(m)2 | 1012 | pF | (pA)2(s)4/(kg)(µm)2 |
电场 | V/m | (kg)(m)/(s)3(A) | 10-6 | V/µm | (kg)(µm)/(s)3(pA) |
通量密度 | C/(m)2 | (A)(s)/(m)2 | 1 | pC/(µm)2 | (pA)(s)/(µm)2 |
| | | | | |
自由空间介电常数等于8.0854E-6pF/µm
表5. MKS制到µMSVfA制电参数换算系数表
电参数 | MKS 制 | 量纲 | 乘数 | µMSVfA制 | 量纲 |
电压 | V | (kg)(m)2/(A)(s)3 | 1 | V | (g)(µm)2/(fA)(s)3 |
电流 | A | A | 1015 | fA | fA |
电荷 | C | (A)(s) | 1015 | fC | (fA)(s) |
导电率 | S/m | (A)2(s)3/(kg)(m)3 | 109 | fS/µm | (fA)2(s)3/(g)(µm)3 |
电阻率 | Ωm | (Kg)(m)3/(A)2(s)3 | 10-9 | -- | (g)(µm)3/(fA)2asmuchas(s)3 |
介电常数 | F/m | (A)2(s)4/(kg)(m)3 | 109 | fF/µm | (fA)2(s)2/(g)(µm)3 |
能量 | J | (kg)(m)2/(s)2 | 1015 | fJ | (g)(µm)2/(s)2 剩女进化论 |
电容 | F | (A)2(s)4/(kg)(m)2 | 1015 | fF | (fA)2(s)4/(g)(µm)2 |
电场 | V/m | (kg)(m)/(s)3(A) | 10-6 | V/µm | (g)(µm)/(s)3(fA) |
通量密度 | C/(m)2 | (A)(s)/(m)2 | 103 | fC/(µm)2 | (fA)(s)/(µm)2 |
| | | | | |
自由空间介电常数等于8.0854E-3fF/µm
14.3.2 加载荷和求解
本步定义分析类型和选项、给模型加载、定义载荷步选项和开始求解。
14.3.2.1 进入求解处理器
命令:/SOLU
GUI:Main Menu>Solution
14.3.2.2定义分析类型
选择下列方式之一:
·GUI:选菜单路径Main Menu>Solution>New Analysis并选择静态分析
·命令:ANTYPE,STATIC,NEW
·
如果你要重新开始一个以前做过的分析(例如,分析附加载荷步),执行命令ANTYPE,STATIC,REST。重启动分析的前提条件是:预先完成了一个静电分析,且该预分析的Jobname. EMAT,Jobname. ESAV和Jobname.DB文件都存在。
14.3.2.3定义分析选项
可以选择波前求解器(缺省)、预条件共轭梯度求解器(PCG)、雅可比共轭梯度求解器(JCG)和不完全乔列斯基共轭梯度求解器(ICCG)之一进行求解:
命令:EQSLV
GUI:Main Menu>Solution>Analysis Options
如果选择JCG求解器或者PCG求解器,还可以定义一个求解器误差值,缺省为1.0-8。
14.3.2.4 加载
静电分析中的典型载荷类型有:
14.3.2.4.1 电压(VOLT)
该载荷是自由度约束,用以定义在模型边界上的已知电压:
命令:D
GUI:Main Menu>Solution>Loads>-Loads-Apply>-Electric-Boundary> -Voltage-
14.3.2.4.2电荷密度(CHRG)
命令:F
GUI:Main Menu>Solution>Loads>-Loads-Apply>-Electric-Excitation>-Charge-On Nodes
14.3.2.4.3面电荷密度(CHRGS)
命令:SF
GUI:Main Menu>Solution>Loads>-Loads-Apply>-Electric-Excitation-Surf Chrg Den-
14.3.2.4.4 Maxwell 力标志(MXWF)
这并不是真实载荷,只是表示在该表面将计算静电力分布,MXWF只是一个标志。通常,MXWF定义在靠近“空气-电介质”交界面的空气单元面上,ANSYS使用Maxwell应力张量法计算力并存储在空气单元中,在通用后处理器中可以进行处理。
命令:FMAGBC
GUI:Main Menu>Solution>-Loads-Apply>-Electric-Flag>-Maxwell Surf-option
14.3.2.4.5 无限面标志(INF)
这并不是真实载荷,只是表示无限单元的存在,INF仅仅是一个标志。
命令:SF
GUI:Main Menu>Solution>-Loads-Apply>-Electric-Flag>-Infinite Surf-option
14.3.2.4.6 体电荷密度(CHRGD)
命令:BF,BFE
GUI:Main Menu>Solution>-Loads-Apply>-Electric-Excitation>-Charge Density-option
另外,还可以用命令BFL、BFL、BFV等命令分别把体电荷密度加到实体模型的线、面和体上。
14.3.2.4.7定义载荷步选项
对于静电分析,可以用其它命令将载荷加到电流传导分析模型中,也能控制输出选项和载荷步选项,详细信息可参见第16章“分析选项和求解方法”
14.3.2.4.8保存数据库备份
使用ANSYS工具条的SAVE_DB按钮来保存一个数据库备份。在需要的时候可以恢复模型数据:
命令:RESUME
GUI:Utility Menu>File>Resume Jobname.db
14.3.2.4.9 开始求解
命令:SOLVE
GUI:Main Menu>Solution>Current LS
14.3.2.4.10结束求解
命令:FINISH
GUI:Main Menu>Finish
14.3.3 观察结果
ANSYS和ANSYS/Emag程序把静电分析结果写到结果文件Jobname.RST中,结果中包括如下数据:
主数据:节点电压(VOLT)
导出数据:
·节点和单元电场(EFX,EFY,EFZ,EFSUM)
·节点电通量密度(DX,DY,DZ,DSUM)
·节点静电力(FMAG:分量X,Y,Z,SUM)
·节点感生电流段(CSGX,CSGY,CSGZ)
通常在POST1通用后处理器中观察分析结果:
命令:/POST1
GUI:Main Menu>General Postproc
对于整个后处理功能的完整描述,见ANSYS基本分析过程指南。
将所需结果读入数据库:
命令:SET,,,,,TIME
GUI:Utility Menu>List>Results>Load Step Summary
如果所定义的时间值处并没有计算好的结果,ANSYS将在该时刻进行线性插值计算。
对于线单元(LINK68),只能用以下方式得到导出结果:
命令:ETABLE
GUI:Main Menu>General Postproc>Element Table>Define Table
命令:PLETAB
GUI:Main Menu>General Postproc>Plot Results>Elem Table
Main Menu>General Postproc>Element Table>Plot Elem Table
命令:PRETAB
GUI:Main Menu>General Postproc>List Results>List Elem Table
Main Menu>General Postproc>Element Table>Elem Table Data
绘制等值线图:
命令:PLESOL,PLNSOL
私有化
GUI:Main Menu>General Postproc>Plot Results>Element Solution