用物理环境法实现静电结构耦合，怎样用Ldread读取静电力啊？

zhanfei

我根据ansys耦合场分析指南中一个静电-结构耦合的例子，这个例子使用ESSOLV宏实现物理环境的迭代的，我想自己是实现两个物理环境的迭代。但是不知道该怎么读取静电场分析得到的静电力。我想在梁上加一个麦克斯韦面标志，但是得到的静电力为0，不知道怎么回事？
**************************************分析指南中例子的源代码*******************************
/batch,list
finish
/clear
/prep7
et,1,121  !air
et,2,121  !Beam
emunit,epzro,8.854e-6
mp,perx,2,1  !空气介电常数
bl=150
bh=2.0
glc=bl/2  
gl=90
gh=1.5
gap=4.5
vltg=120
/pnum,line,1
/pnum,area,1
rectng,0,bl,gap,bh+gap
rectng,glc-gl/2,glc+gl/2,-gh,0
rectng,-10,170,-20,30
aovlap,all
asel,s,area,,1
aatt,1,,1
asel,s,area,,4  !Air
cm,air,area
aatt,2,,2
allsel,all
smrtsiz,2
amesh,1 !mesh beam
mshape,1
amesh,4
asel,s,area,,1
lsla,s
dl,all,,volt,vltg  !Voltage
asel,s,area,,2
lsla,s
dl,all,,volt,0 !Ground
allsel,all
et,1,0
physics,write,electros,electros
physics,clear
et,1,82,,,2
et,2,0
mp,ex,1,170e3
mp,nuxy,1,0.34
dl,4,,ux,0
dl,4,,uy,0
allsel
finish
physics,write,structure,structure
ESSOLV,'electros','structure',2,0,'air',,,,10
finish
physics,read,electros
/post1
set,last
esel,s,mat,,2
etable,fx,fmag,x
etable,fy,fmag,y
ssum
finish


**********************************我自己写的代码*****************************************
/batch,list
finish
/clear
/prep7
et,1,121  !beam
et,2,121  !air
emunit,epzro,8.854e-6
mp,perx,2,1  !空气介电常数
!!!!(一)、建立模型区
bl=150
bh=2.0
glc=bl/2  
gl=90
gh=1.5
gap=4.5
vltg=120
!显示面和边的编号
/pnum,line,1
/pnum,area,1
rectng,0,bl,gap,bh+gap
rectng,glc-gl/2,glc+gl/2,-gh,0
rectng,-10,170,-20,30
aovlap,all
!划分网格
asel,s,area,,1
aatt,1,,1
asel,s,area,,4  !Air
cm,air,area    !空气区设定为一个部件
aatt,2,,2
allsel,all
smrtsiz,2
amesh,1 !mesh beam
mshape,1
amesh,4

!（二）、建立静电物理环境
et,1,0
asel,s,area,,1
lsla,s
dl,all,,volt,vltg  !Voltage
asel,s,area,,2
lsla,s
dl,all,,volt,0 !Ground
allsel,all
!将梁的底边上的节点作为一个部件
lsel,s,,,1
nsll,,1
cm,beambutton,node
ESEL,S,type,,1                        ! 选择单元类型号为2的Beam
CM,beam,ELEM
FMAGBC,'beam'
/com,+++++++++++++++++STARTING beam coordinates------------
alls
/title,Elec Analysis
physics,write,electros,electros
!!!!(三)、创建结构物理环境
physics,clear
et,1,82,,,2
et,2,0
mp,ex,1,170e3
mp,nuxy,1,0.34
dl,4,,ux,0
dl,4,,uy,0
allsel
/title,structural analysis
finish
/solu
antype,static
nlgeom,on
cnvtol,f,,,,-1
physics,write,structure,structure
physics,clear
save

!!!!(四)、静电-结构耦合
loop=25  !最大迭代次数
toler=0.005
*dim,dismax,array,loop    !最大位移数组
*dim,strcri,array,loop   !收敛值数组
*dim,index,array,loop
*do,i,1,loop
/solu
!读入静电物理环境
physics,read,electros
solve
fini
!end of elec portion
!读入结构物理环境
physics,read,structure
/assign,esave,structure,esav
/assign,emat,structure,emat
!如果不是第一次迭代，将重启结构分析
*if,i,gt,1,then
parsave,all   !save parameters for convergence check
resume  !Resume DB-to return original node positions
parresume !Resume parameters needed for convergence check
/prep7
antype,stat,rest
fini
*endif
/solu
solc,off
lsel,s,,,1,3
nsll,,1
esel,s,type,,1
ldread,forc,last,,,,,rst  !读取静电力
alls
rescontrol,,none
solve

还没写完，不过后面的都是后处理的了。到这里不知道该怎么读取静电力。各位大虾帮忙看一下吧！！！！