求助：关于鞍形索网建模

jia35

m=30
n=30
s1=2
s2=2
f1=3
f2=3
a=m*s1/2
b=n*s2/2
/prep7
et,1,link10
mp,ex,1,1.7e11
r,1,10.56e-4,0.001818738
r,2,10.56e-4,0.001818738
*dim,nod,,(m+1)*(n+1)
*do,i,1,(m+1)*(n+1)
nod(i)=0
*enddo
！形成节点
*do,j,1,n+1
*do,i,1,m+1
k=(j-1)*(m+1)+i
x=(i-1)*s1-a
y=(j-1)*s2-b
z=f1*(x/a)**2-f2*(y/b)**2
real,1
*do,j,1,n+1
*do,i,1,m
e1=(j-1)*(m+1)+i
e2=e1+1
*if,nod(e1),eq,1,and,nod(e2),eq,1,then
e,e1,e2
*endif
*enddo
*enddo
real,2
*do,i,1,m+1
*do,j,1,n
e1=(j-1)*(m+1)+i
e2=j*(m+1)+i
*if,nod(e1),eq,1,and,nod(e2),eq,1,then
e,e1,e2
*endif
*enddo
*enddo
numcmp,node
*if,x,ge,0,and,y,ge,0,then
*if,y,le,(-b/a*x+b),then
nod(k)=1
n,k,x,y,z
*endif
*elseif,x,ge,0,and,y,le,0,then
*if,y,ge,(b/a*x-b),then
nod(k)=1
n,k,x,y,z
*endif
*elseif,x,le,0,and,y,ge,0,then
*if,y,le,b/a*x+b,then
nod(k)=1
n,k,x,y,z
*endif
*elseif,x,le,0,and,y,le,0,then
*if,y,ge,(-b/a*x-b),then
nod(k)=1
n,k,x,y,z
*endif
*endif
*enddo
*enddo
这是我在书上看到的建立鞍形索网节点的程序
可是输进去，什么都没有:'(
哪位高人能指点一下，真的万分感谢

shilw

。。。不应该的哦。

jia35

大家快帮帮忙吧