*vwrite语句疑问

      小可根据《ansys工程结构数值分析》一书中9.2.2 节中的例题，加了几道*vwrite语句，按理说应该输出四行结果，但却输出了七行，小可分析是那句简单的循环语句的问题，可是不知错在哪？请各位指教，谢谢！
命令流：
finish$/clear$/prep7
et,1,beam3$mp,ex,1,2.1e11$mp,prxy,1,0.3
mp,dens,1,7800$r,1,0.06,0.00045,0.3
k,1$k,2,6$l,1,2
lesize,all,,,20$lmesh,all
dk,1,all
finish
/solu
antype,2
modopt,lanb,3,,,,1
mxpand,3,,,yes
lumpm,off
solve
*dim,fi,,3$*dim,pfi,,3$*dim,mci,,3$*dim,dai,,3
*do,i,1,3
*get,fi(i),mode,i,freq
*get,pfi(i),mode,i,pfact
*get,mci(i),mode,i,mcoef
*get,dai(i),mode,i,damp
*enddo
finish
/post1
！以下语句为小可加的
*cfopen,ex9.1mode,txt
*vwrite,
('*********模态分析结果***********')
*do,i,1,3,1
*vwrite,i,fi(i),pfi(i),mci(i),dai(i)
(f6.0,2x,e15.6,2x,e15.6,2x,e15.6,2x,e15.6)
*enddo
*cfclose

输出的文本结果为：
*********模态分析结果***********
    1.     0.698154E+01     0.180898E+03     0.180898E+03     0.000000E+00
    1.     0.436267E+02    -0.289074E+02    -0.289074E+02     0.000000E+00
    1.     0.121595E+03     0.103480E+02     0.103480E+02     0.000000E+00
    2.     0.436267E+02    -0.289074E+02    -0.289074E+02     0.000000E+00
    2.     0.121595E+03     0.103480E+02     0.103480E+02     0.000000E+00
    3.     0.121595E+03     0.103480E+02     0.103480E+02     0.000000E+00

!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*do,i,1,3,1
*vwrite,i,fi(i),pfi(i),mci(i),dai(i)
(f6.0,2x,e15.6,2x,e15.6,2x,e15.6,2x,e15.6)
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
这一段改成
*vwrite,fi(i),pfi(1),mci(1),dai(1)
(f6.0,4x,f6.0,4x,f6.0,4x,f6.0,4x)
这样，基本就能看了。
因为原来那样写每循环一次i，接下来的数组fi(i)又自身进行了循环，所以会出现那样的结果。

你好！按你所说的更改后，输出结果为：
*********模态分析结果***********
  122.       10.       10.        0.
这好像也不太对的。

回复 mufan 的帖子

把输出的格式改一下，
*vwrite,sequ,fi(1),pfi(1),mci(1),dai(1)
(f6.0,2x,e15.6,2x,e15.6,2x,e15.6,2x,e15.6)

*********模态分析结果***********
    1.     0.698154E+01     0.180898E+03     0.180898E+03     0.000000E+00
    2.     0.436267E+02    -0.289074E+02    -0.289074E+02     0.000000E+00
    3.     0.121595E+03     0.103480E+02     0.103480E+02     0.000000E+00

果然如此，谢谢！