请帮助看一下用梁单元模拟的管道的命令流是否正确啊？

shimly521

请帮助看一下用梁单元模拟的管道的命令流是否正确啊？管道两端的边界条件该如何确定才比较合适啊？
finish
/clear,start
/filname,pipeline-5,1
/com,structural
/unit,si
/config,nres,5000
*set,mat_nonlinear,0
*set,max_deffect,0.001
/PREP7
l=3.6
e=0.36     
ET,1,BEAM188      
MP,EX,1,2e11
MP,PRXY,1,0.3   
SECTYPE,1,BEAM, CTUBE,,0   
SECOFFSET, CENT  
SECDATA,0.295,0.305,,0,0,0,0,0,0,0  
k,1,
k,2,0,0,l
l,1,2
/view,1,1,1,1
TYPE,1            
MAT,1
REAL,1     
esize,e
lmesh,all
MP,EX,2,2e6
MP,PRXY,2,0.2
MP,EX,3,2e6
MP,PRXY,3,0.2
MP,EX,4,2e6
MP,PRXY,4,0.2
ET,2,COMBIN14,,,0  
ET,3,COMBIN14,,,0  
ET,4,COMBIN14,,,0  
Kx=10000        
Ky=20000
Kz=30000
r,2,Kx
r,3,Ky
r,4,Kz
ngen,2,11,3,11,1,0.1
ngen,2,22,3,11,1,,0.1
ngen,2,33,3,11,1,,,0.1
*do,i,3,11,1
type,2,
real,2
e,i,11+i
*enddo
*do,i,3,11,1
type,3,
real,3
e,i,22+i
*enddo
*do,i,3,11,1
type,4,
real,4
e,i,33+i
*enddo
gplot
save
finish
/solu
d,1,ux
d,1,uy
d,2,ux
d,2,uy
d,2,uz
nsel,s,loc,y,0.1
nsel,a,loc,x,0.1
d,all,ux
d,all,uy
d,all,uz
nsel,s,node,,36,44,1
d,all,ux
d,all,uy
d,all,uz
fk,1,fz,1
allsel
pstres,on
solve
save
finish