工程电磁场.doc

1
求解内外导体间的电位和电场分布，其中外导体接地，内导体的电压为5000V，外球体半径为12CM，，正方形的边长为10CM。
剖分图：
ECTNG,-,,-,,
FLST,3,2,5,ORDE,2
5
FITEM,3,1
FITEM,3,-2
ASBA, 1,P51X
SMRT,6
SMRT,1
MSHAPE,0,2D
MSHKEY,0
!*
CM,_Y,AREA
ASEL, , , , 3
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
FLST,2,4,4,ORDE,2
FITEM,2,1
FITEM,2,-4
/GO
!*
DL,P51X, ,VOLT,0
FLST,2,16,4,ORDE,2
FITEM,2,13
FITEM,2,-28
/GO
!*
DL,P51X, ,VOLT,5000
FINISH
/SOL
/STATUS,SOLU
SOLVE
FINISH
/POST1
!*
/VSCALE,1,1,0
!
!*
PLVECT,EF, , , ,VECT,ELEM,ON,0
!*
6
/EFACET,1
PLNSOL, VOLT,, 0
SAVE
SAVE
SAVE
SAVE
SAVE
SAVE
FINISH
! /EXIT,NOSAV
：
7
程序流程图：
8
程序：
定义函数m1(a,b):
function y=m1(a,b)
% m1(a,b) is a function;
% y=m1(a,b)
% x1 is variable;
% x2 is variable;
% a is constant;
% b is constant;
x1=a;x2=b;
y=[126743*x2-127627*x2*exp(-x1)+127627*126743*exp(-x1)-127627*126743;126743*x2-129988*x2*exp(-4*x1)+129988*126743*exp(-4*x1)-126743*129988];
定义函数m11(a,b):
function y=m11(a,b)
x1=a;
x2=b;
y=-127627*exp(-x1)*(126743-x2);
定义函数m12(a,b):
function y=m12(a,b)
x1=a;
x2=b;
y=126743-127627*exp(-x1);
定义函数m21(a,b):
function y=m21(a,b)
x1=a;
x2=b;
y=(-4)*129988*exp(-4*x1)*(126743-x2);
定义函数m22(a,b):
function y=m22(a,b)
x1=a;
x2=b;
y=126743-129988*exp(-4*x1);
定义函数mJ(a,b):
function y=mJ(a,b)
% mJ(x1,x2) is a function
% y=mJ(a,b)
% x1 is variable
% x2 is variable
% a is constant
% b is constant
% A is matrix
x1=a;
9
x2=b;
y=[m11(x1,x2) m12(x1,x2);m21(x1,x2) m22(x1,x2)];
定义函数mJd(a,b):
function y=mJd(a,b)
% mJd(x1,x2) is a function
% y=mJd(a,b)
% x1 is variable
% x2 is variable
% a is constant
% b is constant
% A1 is matrix
% A2 is matrix
x1=a;
x2=b;
error=;
while erro