!这是该书上的一个例子,程序运行了很多遍出现了两个大问题(出现出有注释),程序如下:
/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
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,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 !1、在加载荷后求解时,会出现327个警告,在错误报告中说 Finer mesh required near center line for element即在线单元划分时有问题
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
*cfoen,u10,txt,d:\ansys\result,
*vwrite,u1(1)
(f10.8)
*cfclos
finish
/fillname,logwell_K
resum
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
*cfoen,u20,txt,d:\ansys\result,
*vwrite,u2(1)
(f10.8)
*cfclos
finish
/fillname,logwell_K
resum
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
*cfoen,u30,txt,d:\ansys\result,
*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) !2、矩阵求解时显示Matrix is singular for SOLV command |