本帖最后由 livshich 于 2017-11-26 14:30 编辑
*dim,canshu,,3,4
*dim,weight,,3,1
*dim,vtot,,3,1
*vread,canshu(1,1),canshu,txt,,jik,4,3
(f4.0,f4.0,f4.0,f4.0)
*creat,define_length,mac
/prep7
et,1,link1
r,1,a1
r,2,a2
r,3,a3
mp,ex,1,2.1e6
mp,prxy,1,0.3
n,1,-b,0,0
n,2,0,0,0
n,3,b,0,0
n,4,0,-1000,0
/pnum,node,1
/number,0
nplot
real,1
e,1,4
real,2
e,2,4
real,3
e,3,4
finish
/solu
d,1,all,0,,3
f,4,fx,200000
f,4,fy,-200000
solve
finish
*end
*do,i,1,3
a1=canshu(i,1)
a2=canshu(i,2)
a3=canshu(i,3)
b=canshu(i,4)
define_length
/post1
set,last
etable,evol,volu
ssum
*get,vtot(i,1),ssum,,item,evol
rho=2.85e-4
weight(i,1)=rho*vtot(i,1)
*enddo |