利用ANSYS软件进行电场边值问题的求解.docVIP

实验报告 课程名称： 工程电磁场与波 指导老师： 熊素铭 成绩：__________________ 实验名称： ANSYS软件仿真练习 实验类型：________________ 一、实验目的和要求（必填） 二、实验内容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 七、讨论、心得 一、实验目的 1. 初步接触Ansys软件,用其解决简单的初值边界问题. 2. 加深对初值边界问题的理解. 3. 熟悉Ansys软件,为后面较复杂的仿真实验打下基础. 二、实验内容和原理: 1. 分析二维电场电势的初值边界问题:其模型如图: 其中a = 0.1m; b=0.1m; =200V. 2. 分析其电势的分布,画出等势图; 3. 得到特殊电的电势值,与标准公式进行对比,标准公式为: 三、主要仪器设备 ANSYS10.0 Matlab6.5 四、操作方法和实验步骤: 1. 用ANSYS10.0进行仿真实验: (1)设置材料属性为相对介电常数为1,单元类型为plane121,单元参考号为1. (2) 建立几何模型,先建立关键点,考虑到200V点和0V点不能重合,于是设置一个点,将200V和0V隔开.并在关键点间连线,建立好模型具体几何模型图为: (3) 设置几何模型的属性,材料号为1,单元号为1,单元类型为plane121,单元坐标系为直角坐标系. (4) 剖分,建立网格,应用smart size精度为2和7时进行对比. (4) 对线进行赋初值处理:61线、12线、23线赋值电压为常量0V.45线赋值电压为常量200V. (5) 分别得到两种网格密度的等势图. (6) 分别得到两种网格密度下随机取下5个点的电势值,以便接下来的比对. 2. 用Matlab 6.5进行仿真. 编写.m文件.用标准公式计算产生各点的电势值,并画出和1000条等势线下的等势图.其中.m文件的代码为: function [z] = phi(phi0, a, b, x, y, N) z=0; for i=0:100 c=pi*(2*i+1)/0.1; z=z+(1/((2*i+1)*sinh(c*0.1))).*(sin(c.*x).*sinh(c.*y)); end z=z.*(4*200/pi); 得到Matlab仿真下对应5个点的电势,并于上面的结果进行比对. 所用.m文件代码为: phi0=200; a=0.1; b=0.1; rx=0:a/100:a; ry=0:b/100:b; n=100; [x, y]=meshgrid(rx, ry); z=phi(phi0, a, b, x, y, n); m=10; [C, h]=contour(x, y, z, m); 3. 分析误差产生的原因 五、实验数据记录和处理: 1. ANSYS仿真下”smart mesh”密度为7下的等势图为: ANSYS仿真下”smart mesh”密度为2下的等势图为: matlab仿真下1000条等势线的等势图: 比对Matlab和ANSYS的电势值数据(ANSYS网格精细程度为2) 点标号 ANSYSX值/m ANSYSY值/m MatlabX值/m MatlabY值/m ANSYS电势/V Matlab电势/V 1 0.089132 0.098119 0.089 0.098 177.58 176.6130 2 0.098786 0.097602 0.099 0.098 54.441 58.9897 3 0.098209 0.080892 0.098 0.081 11.131 12.5512 4 0.014036 0.091101 0.014 0.091 124.98 124.4652 5 0.08221 0.075746 0.082 0.076 71.309 72.7436 比对Matlab和ANSYS的电势值数据(ANSYS网格精细程度为7) 点标号 ANSYSX值/m ANSYSY值/m MatlabX值/m MatlabY值/m ANSYS电势/V Matlab电势/V 1 0.096019 0.097795 0.096 0.098