10-2_2.txt

ANSYS电磁分析的命令流文件

/fillname,logwell_K
keyw,magelc,1

/prep7
r0=0.0445                      ! 仪器半径
k,1,0,-35
k,2,0,-1.8
k,3,r0,-1.8
k,4,r0,-0.7
k,5,r0,-0.5
k,6,r0,-0.1
k,7,r0,-0.075
k,8,r0,0.075
k,9,r0,0.1
k,10,r0,0.5
k,11,r0,0.7
k,12,r0,1.8
k,13,0,1.8
k,14,0,35
k,15,0.1,35
k,16,0.1,-35
a,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16   ! 挖去仪器的井眼
rectng,0.1,35,1,35
rectng,0.1,35,-1,1
rectng,0.1,35,-1,-35          ! 右边地层
rectng,0.1,50,1,51
rectng,35,50,-1,1
rectng,0.1,50,-1,-51          ! 右边远场地层
rectng,0,0.1,35,51
rectng,0,0.1,-35,-51          ! 上下远场井
aovlap,all
aglue,all                     ! 布尔操作
nummrg,all
numcmp,all
/pnum,area,1
/replot
et,1,67,,,1                  ! 设置单元类型为plane67
mp,rsvx,1,1                  ! 设置两个相同的电阻率
mp,rsvx,2,1
asel,s,area,,1,3,2
asel,a,area,,5,9
aatt,1,,1                    ! 给围岩和井内泥浆附属性
allsel
asel,s,area,,2,4,2
aatt,2,,1                   ! 给地层附属性
allsel
lesize,1,,,80,0.05,,,,1     ! 设置线单元个数
lesize,2,,,4
lesize,3,,,80
lesize,4,,,15
lesize,5,,,35
lesize,6,,,2
lesize,7,,,12
lesize,8,,,2
lesize,9,,,35
lesize,10,,,15
lesize,11,,,80
lesize,12,,,4
lesize,13,,,80,20,,,,1
lesize,14,,,1
lesize,15,,,1
lesize,24,,,1
lesize,27,,,1
lesize,16,,,85,20,,,,0
lesize,28,,,85,0.05,,,,0
lesize,29,,,14,,,,,0
lesize,17,,,14,,,,,0
lesize,18,,,85,20,,,,0
lesize,32,,,85,0.05,,,,0
lesize,33,,,14,,,,,0
lesize,34,,,14,,,,,0
lesize,30,,,35,,,,,1
!lesize,31,,,35,,,,,1
lesize,20,,,14,,,,,0
lesize,19,,,14,,,,,0
lesize,21,,,14,,,,,0
lesize,38,,,14,,,,,0
lesize,23,,,7,,,,,0
lesize,35,,,7,,,,,0
lesize,36,,,7,,,,,0
lesize,26,,,7,,,,,0
lccat,19,20
lccat,17,29
lccat,33,34
lccat,21,38
mshape,1,2d              ! 网格划分
mshkey,2
amesh,3,5,2
mshape,1,2d
mshkey,0
amesh,4
mshape,1,2d
mshkey,1
amesh,7,8
mshape,1,2d
mshkey,1
amesh,6
mshape,1,2d
mshkey,1
amesh,1,2
smrtsize,1
mshape,1
mshkey,0
amesh,9
lsel,,,,7
nsll,s,1
cp,1,volt,all             ! 主电极上节点电压自由度耦合
allsel
lsel,,,,9
lsel,a,,,5
nsll,s,1
cp,2,volt,all            ! 屏蔽电极上节点电压自由度耦合
allsel 
lsel,,,,11
lsel,a,,,3
nsll,s,1
cp,3,volt,all            ! 屏蔽电极上节点电压自由度耦合
allsel
*dim,aa,,3               ! 取出三个电极上关键点同位置处节点号并赋予一个数组
*vget,aa(1),kp,7,attr,node
*vget,aa(2),kp,9,attr,node
*vget,aa(3),kp,11,attr,node
finish

/solu
lsel,s,loc,x,50
lsel,a,loc,y,-51
lsel,a,loc,y,51
dl,all,,volt,0            ! 施加无穷远边界条件
allsel,all
fk,7,amps,1               ! 在主电极关键点上施加电流激励
fk,9,amps,0
fk,11,amps,0
sbctran
solve                     ! 求解运算
save                      ! 保存数据库
finish

/post1
set,last
*dim,u1,,3               ! 取出三个电极上的电压值并赋予一个数组
*vget,u1(1),node,aa(1),volt
*vget,u1(2),node,aa(2),volt
*vget,u1(3),node,aa(3),volt
*cfopen,u10,txt,d:\ansys\result,   ! 按指定路径（该路径要存在）输出一分场结果
*vwrite,u1(1)
(f10.8)
*cfclos
finish

/fillname,logwell_K
resum                    ! 恢复数据库
/solution
fk,7,amps,0
fk,9,amps,1              ! 在屏蔽电极关键点上施加电流激励
fk,11,amps,0
solve                    ! 求解运算
save                     ! 保存数据库
finish

/post1
set,last
*dim,u2,,3                         ! 取出三个电极上的电压值并赋予一个数组
*vget,u2(1),node,aa(1),volt
*vget,u2(2),node,aa(2),volt
*vget,u2(3),node,aa(3),volt
*cfopen,u20,txt,d:\ansys\result,    ! 按指定路径输出二分场结果
*vwrite,u2(1)
(f10.8)
*cfclos
finish

/fillname,logwell_K
resum                     ! 恢复数据库
/solution
fk,7,amps,0
fk,9,amps,0
fk,11,amps,1              ! 在屏蔽电极关键点上施加电流激励
solve
save
finish

/post1
!set,last
*dim,u3,,3                        ! 取出三个电极上的电压值并赋予一个数组
*vget,u3(1),node,aa(1),volt
*vget,u3(2),node,aa(2),volt
*vget,u3(3),node,aa(3),volt
*cfopen,u30,txt,d:\ansys\result,append  ! 按指定路径输出二分场结果
*vwrite,u3(1)
(f10.8)
*cfclos
*dim,u1,,3
*dim,u2,,3
*vread,u1(1),u10,txt,d:\ansys\result, ! 三个分场解完，读入前面二个分场的值
(f10.8)
*vread,u2(1),u20,txt,d:\ansys\result,
(f10.8)
*dim,a,,2,2
*dim,b,,2
*dim,c,,2
a(1,1)=u2(1)-u2(2)      ! 将三个分场线性组合叠加，求解方程组得比例系数
a(1,2)=u3(1)-u3(2)
a(2,1)=u2(1)-u2(3)
a(2,2)=u3(1)-u3(3)
b(1)=u1(2)-u1(1)
b(2)=u1(3)-u1(1)
*moper,c(1),a(1,1),solve,b(1)
k=1/(u1(1)+c(1)*u2(1)+c(2)*u3(1))      ! 求仪器常数
*cfopen,c0,txt,d:\ansys\result,        ! 输出比例系数到一个文本文件
*vwrite,c(1)
(4f15.10)
*cfclos
*cfopen,k,txt,d:\ansys\result,         ! 输出仪器常数到一个文本文件
*vwrite,k
(f15.10)
*cfclos
finish