作业三 规则平面域的自动剖分.docVIP

作业三：规则平面域的自动剖分 设长直接地金属槽的横截面如图5-3所示，其侧壁与底面电位均为零，顶盖电位为10（相对值） 分析： 对于此槽中间区段的电场分析，可理想化为二维平行平面场问题。选定直角坐标系如图5-3所示。为提高计算效率，注意到槽内电场分布的对称性，可将计算的场域归结为左半（右半）空间，因而槽内电位函数满足拉普拉斯方程,构成如下的混合型边值问题： =0 （0x0.5a, 0y0.5a） |(x=0,0y0.5a)=|(0x0.5a,y=0)=0 |(0x0.5a,y=0.5a)=10 |(x=0.5a,0y0.5a)=0 根据变分原理，可知与上述边值问题等价的条件变分问题是： J[]=]dxdy=min |(x=0,0y0.5a)=|(0x0.5a,y=0)=0 |(0x0.5a,y=0.5a)=10 按有限元法的手算步骤是： 场域剖分，并对所有的节点和三角元逐个按一定的顺序编号。剖分后整理得出下列原始数据信息： ANSYS有限元分析部分 第一部分： 首先建立平面模型如图1： 图1 实体模型 该长方形的长为10，宽为5。在本模型中视金属槽壁的厚度为零，槽内为空气，其介电常数为1，在材料的参数设置中将PERX的值设置为1。 在定义了单元类型和材料参数（无实常数）后即可对该模型进行网格划分。对该模型进行网格划分后的图形如图2： 在网格划分的控制中，将网格该模型划分成边长为1的正方形，这样计算得到的接点的电位值便如和其他方式计算的结果进行比较： 图2 网格划分 在网格划分完成后，ANSYS的前处理工作基本完成。 第二部分： ANSYS的加载数据和求解。在模型的上边加载10V的电位值，在两个侧边和底边加载0V的电压值。这样就可进行求解了。加载求解后的实体模型如图3： 图3 加载求解后 第三部分： 结果分析，完成计算以后，可以通过ANSYS的后处理模块来查看计算得到的结果，在本题中最主要的目的之一就是查看电位值在模型中的分布以及接点的电位值，电位等值线分布图见最后彩图。 图5 接点图 图6 电位矢量图 以上就是本题的ANSYS分析的全过程，还有一些相关物理量的分析不是题目考察重点，由于篇幅有限，这里不一一列出来！ 以上通过控制网格划分的计算是比较粗糙的，可见软件不是万无一失的，以下给出控制网格划分（X轴200等份，Y轴400等份）电位等值线分布图！ （X轴5等份，Y轴10等份） （X轴200等份，Y轴400等份） 自动剖分技术 对应于边界为分段直线构成，本法是一种简单、实用的自动剖分方法。设这类规则平面 域的基本子域为矩形域，现构造如图1所示采用三角元的自动剖分方案，其单元可归并为图1所示的Ⅰ、Ⅱ两类单元生成所需的按单元局部编码序形成的各单元节点编码：NODEI（E）、NODEJ（E）、NODEM（E），同时，亦可直观地生成所需的节点坐标：X（N）、Y（N）。 自动剖分C程序 C主程序 //*************************** //** 有限元计算程序 ** //*************************** #includeiostream #includefstream using namespace std; struct node { double x; double y; };//节点结构，包括节点的x、y坐标 struct element { double nodei; double nodej; double nodem; };//单元结构，包括三角元e的三顶点编号i,j,k编号 struct boundarynodetype { int bnode; double bptal; };//强制边界条件结构，包括节点编号bnode及其电位值bptal class EMF { int NSX1,NSX2,NSY1,NSY2; int NSX,NSY,LASTND,LASTEL; double RX1,RX2,RY1,RY2; int *Address; double *KMatric,*PMatric; node *Node; element *Element; boundarynodetype *Boundarynode; public: EMF(int,int