工程电磁场数值计算教程方案.docx

杨旭东 3120103501 工程电磁场数值计算 大作业报告 一、大作业要求 运用FEM法求解算题5—8，删去要求（2），设其具有平行平面磁场分布的特征。 作业题目如下所示： 二、问题分析及建立模型 根据P149对平行平面场的静电场和磁场统一的数学模型的描述 我们可以得到此问题对应的偏微分方程及相应的定解问题为： 进而可以求得此题对应的泛函及等价的变分问题为： 根据以上条件，我们可以把此题与例5-2作比较，他们的边界条件形式已经基本一致了，所以我们可以利用EMF2D的程序对此题进行计算。 下面所以下我们的主要解题思路。 1、由于是一个圆形区域，且是对称的，所以我们只需求1/4圆周即可。我们运用圆域剖分程序CAMG对整个区域进行剖分。这里我们需要注意的是最外层的边界条件，我们选用选定10倍半径，即1米，进行三段剖分。 2、运用程序EMF2D，把圆域剖分出来的结果当作此程序的输入。需要注意的是需要对剖分出来的最外层的点，进行“手动输入”。我们需要注意两个程序的输入输出的格式进行统一，修改EMF2D的强制边界条件程序FB。  三、程序及结果 1、圆域剖分 我们并没有改变什么CAMG程序，程序如下 PARAMETER(NSEG=5,LND=1000,LEL=600) PARAMETER(NINPT=10,NOUPT=20) DIMENSION RADIUS(NSEG+1),Q(NSEG) INTEGER SECTION_R,SUB_SEG(NSEG),SUMELEM,SUMNODE COMMON/NE/X(LND),Y(LND),II(LEL,3) OPEN(NINPT,FILE=CAMGIN1.DAT,FORM=FORMATTED) READ(NINPT,70)SECTION_R 70 FORMAT(I3) READ(NINPT,90)RADIUS(1) 90 FORMAT(F6.2) DO 5 I=1,SECTION_R READ(NINPT,80)RADIUS(I+1),Q(I),SUB_SEG(I) 80 FORMAT(F6.2,F6.2,I3) 5 CONTINUE PI=3.14159 R=0. X(1)=0. Y(1)=0. NR=0 IJ=0 N=1 DO 30 I=1,SECTION_R NR=NR+SUB_SEG(I) SUM_F=0 DO 10 K=1,SUB_SEG(I) SUM_F=SUM_F+Q(I)**(K-1) 10 CONTINUE DO 20 J=1,SUB_SEG(I) IJ=IJ+1 DR=(RADIUS(I+1)-RADIUS(I))*Q(I)**(J-1)/SUM_F R=R+DR DF=0.5*PI/IJ AF=0 IJ1=IJ+1 DO 40 K=1,IJ1 N=N+1 X(N)=R*COS(AF) Y(N)=R*SIN(AF) AF=DF+AF 40 CONTINUE 20 CONTINUE 30 CONTINUE SUMNODE=N CALL IIFILL(1,1,2,3) IL=2 IR=4 IE=1 DO 50 K=2,NR K1=K-1 DO 60 KK=1,K1 IE=IE+1 CALL IIFILL(IE,IL,IR,IR+1) IE=IE+1 CALL IIFILL(IE,IL,IR+1,IL+1) IL=IL+1 IR=IR+1 60 CONTINUE IE=IE+1 CALL IIFILL(IE,IL,IR,IR+1) IL=IL+1 IR=IR+2 50 CONTINUE SUMELEM=IE OPEN(NOUPT,FILE=CAMGOUT.DAT) WRITE(NOUPT,(2I6))SUMNODE,SUMELEM WRITE(NOUPT,100)(N,X(N),Y(N),N=1,SUMNODE) 100 FORMAT(3(1X,I5,2F10.5)) WRITE(NOUPT,101)(I,(II(I,J),J=1,3),I=1,SUMELEM) 101 FORMAT(3(4X,4I5)) END SUBROUTINE IIFILL(IE,I1,I2,I3) PARAMETER(LND=1000,LEL=600) COMMON/NE/X(LND),Y(LND),II(LEL,3) II(IE,1)=I1 II(IE,2)=I2 II(IE,3)=I3 RETURN END 我们的输入数据如下： 3 0. 0.1 1.0 7 0.2 1.0 8 1.0 1.0 6 由输入可以知道我们内环分7段，中环分8段，外环分6段。 得到的输出结果CAMGOUT结果如下： 190 324 1 0.00000