您的位置:网站首页 > Ansys教程

Ansys排水管道计算

时间:2008-09-15 09:49:49 来源:
埋设在地下的排水管道  
rpc管道混凝土管道用solid45单元模拟,link8模拟纵向钢筋和箍筋,直接建立单元,采用复制的方式直接生成实体。让link单元和solid单元共用节点,(nummrg和numcmp命令)这样就省去了耦合和添加约束方程,预应力采用施加初始应变的方式。约束采用简直方式。计算以后在支座处有比较名显的应力集中现象。这是我做的思路,有命令流,有问题大家可以交流。
附:ansys log file
/prep7
et,1,solid45
mp,ex,1,4.5e10   
mp,prxy,1,0.25   
et,2,link8   
mp,ex,2,2.0e11   
mp,prxy,2,0.3   
R,1,137.4e-6, ,  
csys,1
n,1,0.52
n,2,0.56
n,3,0.6  
NGEN,2,3,P51X, , , ,18, ,1,
NGEN,2,6,P51X, , , , ,0.05,1,  
/VIEW, 1 ,1,1,1  
e,1,2,5,4,7,8,11,10  
e,2,3,6,5,8,9,12,11  
eplot   
ENGEN,40,40,12,P51X, , , , , , , , ,18, ,   
nummrg,all   
numcmp,all   
ENGEN,4000,100,120,P51X, , , , , , , , , ,0.05,  
nummrg,all   
numcmp,all   
nsel,s,loc,x,0.56   
nsel,r,loc,z,0.25   
nplot,1  
l,302,305   
*repeat,19,3,3   
l,359,302   
nsel,s,loc,x,0.56   
nsel,r,loc,z,0.55   
nplot,1  
e,662,665   
*repeat,19,3,3   
e,719,662   
nsel,s,loc,x,0.56   
nsel,r,loc,z,0.85   
nplot,1  
e,1022,1025  
*repeat,19,3,3   
e,1079,1022  
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,z,1.15   
nplot,1  
e,1382,1385  
*repeat,19,3,3   
e,1439,1382  
nsel,s,loc,x,0.56   
nsel,r,loc,z,1.45   
nplot,1  
e,1742,1745  
*repeat,19,3,3   
e,1799,1742  
nsel,s,loc,x,0.56   
nsel,r,loc,z,1.75   
nplot,1  
e,2102,2105  
*repeat,19,3,3   
e,2159,2102  
nsel,s,loc,x,0.56   
nsel,r,loc,z,2.05   
nplot,1  
e,2462,2465  
*repeat,19,3,3   
e,2519,2462  
nsel,s,loc,x,0.56   
nsel,r,loc,z,2.35   
nplot,1  
e,2822,2825  
*repeat,19,3,3   
e,2879,2822  
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,z,2.65   
nplot,1  
e,3182,3185  
*repeat,19,3,3   
e,3239,3182  
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,z,2.95   
nplot,1  
e,3542,3545  
*repeat,19,3,3   
e,3599,3542  
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,z,3.25   
nplot,1  
e,3902,3905  
*repeat,19,3,3   
e,3959,3902  
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,z,3.55   
nplot,1  
e,4262,4265  
*repeat,19,3,3   
e,4319,4262  
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,z,3.85    #p#分页标题#e#
nplot,1  
e,4622,4625  
*repeat,19,3,3   
e,4679,4622  
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,z,4.15   
nplot,1  
e,4982,4985  
*repeat,19,3,3   
e,5039,4982  
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,z,4.45   
nplot,1  
e,5342,5345  
*repeat,19,3,3   
e,5399,5342  
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,z,4.75   
nplot,1  
e,5702,5705  
*repeat,19,3,3   
e,5759,5702  
SAVE
ALLSEL,ALL   
ESEL,S,TYPE,,2   
EPLOT   
ALLSEL,ALL   
nsel,s,loc,x,0.56   
nsel,r,loc,y,0
e,2,8   
e,8,122  
e,122,182   
*repeat,98,60,60
   
nsel,s,loc,x,0.56   
nsel,r,loc,y,36  
nplot,1  
e,14,17  
e,17,128
e,128,188   
*repeat,98,60,60
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,y,72  
nplot,1  
e,26,29  
e,29,134
e,134,194   
*repeat,98,60,60
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,y,108
nplot,1  
e,38,41  
e,41,140
e,140,200   
*repeat,98,60,60
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,y,144
nplot,1  
e,50,53  
e,53,146
e,146,206   
*repeat,98,60,60
nsel,s,loc,x,0.56   
nsel,r,loc,y,180
nplot,1  
e,62,65  
e,65,152
e,152,212   
*repeat,98,60,60
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,y,216
nplot,1  
e,74,77  
e,77,158
e,158,218   
*repeat,98,60,60
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,y,252
nplot,1  
e,86,89  
e,89,164
e,164,224   
*repeat,98,60,60
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,y,288
nplot,1  
e,98,101
e,101,170   
e,170,230   
*repeat,98,60,60
SAVE
nsel,s,loc,x,0.56   
nsel,r,loc,y,324
nplot,1  
e,110,113   
e,113,176   
e,176,236   
*repeat,98,60,60
SAVE
allsel,all   
ESEL,S,TYPE,,2   
EPLOT   
nummrg,all   
numcmp,all   
nlist     
ESEL,S,TYPE,,2  
BFE,P51X,TEMP,1,-375, , ,!施加预应力,已改为用初应变initial strain施加   
  
施加其它荷载和约束以及后处理的命令省略,在后处理中注意转换结果坐标系到cylinder coordinate system     
  
结果1:每隔20cm配一道箍筋,施加1250 a的预应力。
  
   
   
结果2:每隔30cm配一道箍筋,施加1250 a的预应力。