/batch,list
/title, DC-current distribution in a thin conducting disk
a=20.e-2 ! disk radius, m
b=10.e-2 ! electrode distance from the center, m
I=1.e-3 ! current, A
rho=100 ! resistivity, Ohm*m
/nopr
/PREP7
! Model and meshing
et,1,PLANE230 ! electric element type
mp,rsvx,1,rho ! resistivity
cyl4,0,0,0,0,a,360 ! circular area
esize,,64
msha,1,2-D ! mesh with triangles
amesh,1
csys,1
d,node(0,0,0),volt,0
f,node(b,0,0),amps,-I
f,node(b,180,0),amps,I
fini
/solu
antype,static ! steady-state current conduction
solve
fini
/post1
plnsol,volt ! plot electric potential
plnsol,jc,sum ! plot current density vector magnitude
*dim,value,,10,2
*dim,coord,,10,2
*do,i,1,10
r=a $ phi=(i-1)*10
coord(i,1)=r
coord(i,2)=phi
value(i,1)=volt(node(a,(i-1)*10,0))
*enddo
上面这个是ANSYS中help中的例子,但是结果计算却不太正确 |