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