求助大神！ANSYS模态综合之困惑

写了个模态综合法的程序 与workbench的结果 相对误差太大，不知道是为什么，希望大家帮忙看看程序哪里出错了！！！
!generate pass
finish
/clear
/FILNAME,gen
/PREP7
/units,MPA
ET,1,PLANE82
KEYOPT,1,3,3
R,1,5
MP,EX,1,7e4
MP,PRXY,1,0.3
MP,DENS,1,2.64e-9
RECTNG,-200,200,-100,100
ESIZE,10
MSHAPE,0,2D
MSHKEY,0
amesh,all
finish

/SOLU
ANTYPE,subst
csys,0
seopt,gen,2
nsel,s,loc,x,-200
nsel,a,loc,x,200
M,all,all
allsel,all
save
solve
finish

!use pass
/clear
/FILNAME,use
/PREP7
/UNITS,MPA
ET,1,matrix50
r,1,5
MP,EX,1,7e4
MP,PRXY,1,0.3
MP,DENS,1,2.64e-9
ET,2,shell63
r,2,5
MP,EX,2,7e4
MP,PRXY,2,0.3
MP,DENS,2,2.64e-9
type,2
wpoffs,200
wprota,,,90
rectng,0,200,-100,100
wpoffs,,,-400
rectng,0,200,-100,100
wpoffs,,,400
wprota,,,-90
wpoffs,-200
type,2
ESIZE,10
MSHKEY,0
amesh,all
TYPE,1
SETRAN,gen,0,10000,'part1'
SE,part1
SETRAN,part1, ,1000,'part2',' ',' ',,,-100
SE,part2
eplot
nsel,s,loc,x,-200
nsel,a,loc,x,200
nsel,r,loc,z,0
CPINTF,ALL
allsel,all
nsel,s,loc,x,-200
nsel,a,loc,x,200
nsel,r,loc,z,-100
CPINTF,ALL
allsel,all
FINISH
/SOLU
ANTYPE,modal
modopt,lanb,10
mxpand,10,,,1
nsel,s,loc,z,-200
d,all,all
ALLSEL,all
save
solve
finish