高密度电法温纳装置模拟方法与结果分析

finish

/clear

/filename,yanrong

/title,analyzing the yanrong using ansys

/prep7

……

*set,currt,1

*set,pi,acos(-1)

*set,zn,19

*set,m,60

*dim,ab,,zn,1

*vread,ab(1),ab,dat,,jik,1,zn

(f5.0)

*dim,kk,,zn,1

*dim,r,,zn*(m-3)-3*zn*(zn-1)/2,1

*dim,v,, zn*(m-3)-3*zn*(zn-1)/2),1

*do,n,1,zn,1

*do,cd,3/2*n*a-(m-1)/2*a,(m -1)/2*3-3/2*n*a,a

/solu
!求解运算

lsclear,all

nsel,s,ext
!选择模拟无穷远边界上的节点

nsel,u,loc,y

*get,nnod,node,,count

*get,nmin,node,,num,min

*do,j,1,nnod,1
!逐点赋电位值

dist1=sqrt((nx(nmin)-(-1.5*ab(n)+cd))**2+ny(nmin)**2)

dist2=sqrt((nx(nmin)-(1.5*ab(n)+cd))**2+ny(nmin)**2)

v0=1*1000*(log(1/dist1)-log(1/dist2))/pi

d,nmin,volt,v0

nmin=ndnext(nmin)

*enddo

nsel,all

nsel,s,node,,node(-1.5*ab(n)+cd,0,0)

f,all,amps,1
!供电点赋电流值

nsel,all

nsel,s,node,,node(1.5*ab(n)+cd,0,0)

f,all,amps,-1

allsel,all

solve

finish

/post1

v((cd-1*(3/2*n*3-(m-1)/2*3))/3+(n-1)*(m-3)-3/2*(n-1)*(n-1-1)+1)=volt(node(-1*0.5*ab(n)+cd,0,0))-volt(node(0.5*ab(n)+cd,0,0))

kk(n)=pi/log(4/1)
！温纳a装置的装置系数

r((cd-1*(3/2*n*3-(m-1)/2*3))/3+(n-1)*(m-3)-3/2*(n-1)*(n-1-1)+1)=v((cd-1*(3/2*n*3-(m-1)/2*3))/3+(n-1)*(m-3)-3/2*(n-1)*(n-1-1)+1)*kk(n)

*enddo

*enddo

Finish