预应力建模(曲线预应力索及直线预应力索).docxVIP

2 独立建模耦合法 该法的基本思想是实体和力筋独立建几何模型，分别划分单元，然后采用耦合方程将力筋单元和实体单元联系起来，这种方法是基于有限元模型的处理。其基本步骤如下： ①建立实体几何模型（不考虑力筋）； ②建立力筋线的几何模型（不考虑体的存在）； ③将几何模型按一定的要求划分单元（这时也是各自独立的）； ④选择所有力筋线； ⑤选择与上述力筋相关的节点(nsll命令)，并定义选择集； ⑥将上述力筋节点存入数组； ⑦选择所有节点，并去掉⑤中的节点集(这时是除力筋节点外的所有节点)； ⑧按力筋节点数组搜寻所有最近的实体节点号，并存入数组中； ⑨耦合力筋节点与最近的节点，一一耦合(cp命令)（不能使用cpintf命令，这样可能耦合其它节点，且容易不耦合） ⑩选择所有，并施加边界条件和荷载，可以求解了。 这种方法建模特别简单，耦合处理也比较简单(APDL要熟悉些)，缺点是当实体单元划分不够密时，力筋节点位置可能有些走动，但误差在可接受范围之内！这种方法是解决力筋线形复杂且力筋数量很多时的较佳方法。 预应力简支梁弹性分析--体线独立耦合法示例 !----------------------------------------- /prep7 eg=2e5 ag=140 eh=4e4 r0=9345 yyl=200000 et,1,link8 et,2,solid95 r,1,ag,yyl/eg/ag*1.036258 r,2 mp,ex,1,eg mp,prxy,1,0.3 mp,ex,2,eh mp,prxy,2,0.2 blc4,,,100,200,3000 /view,1,1,1,1 /ang,1 vplot !------------定义力筋线 ksel,all *get,kp0,kp,0,num,max lsel,none k,kp0+1,50,160 k,kp0+2,50,160,3000 k,kp0+3,50,800,1500 larc,kp0+1,kp0+2,kp0+3,r0 kdele,kp0+3 *get,line1,line,0,num,min !-------------定义约束 lsel,s,loc,z,0 lsel,r,loc,y,0 dl,all,,uy lsel,s,loc,z,3000 lsel,r,loc,y,0 dl,all,,all allsel,all !-----------单元划分 lsel,s,,,line1 latt,1,1,1 lesize,all,,,50 lmesh,all vsel,all vatt,2,2,2 lsel,s,loc,z,0 lsel,r,loc,y,10,140 lesize,all,,,8 lsel,s,loc,z,0 lsel,u,loc,y,10,140 lesize,all,,,4 lsel,s,loc,y,0 lsel,r,loc,x,0 lesize,all,,,50 vsweep,all allsel,all !耦合自由度 lsel,s,,,line1 nsll,s,1 cm,cmljnod,node *get,max1,node,0,count *dim,ojd,,max1 *dim,jd,,max1 *get,nod1,node,0,num,min ojd(1)=nod1 *do,i,2,max1 ojd(i)=ndnext(ojd(i-1)) *enddo allsel,all nsel,all cmsel,u,cmljnod *do,i,1,max1 nod1=ojd(i) j=nnear(nod1) jd(i)=j *enddo nsel,all ji=1 *do,i,1,max1 cp,ji,ux,ojd(i),jd(i) cp,ji+1,uy,ojd(i),jd(i) cp,ji+2,uz,ojd(i),jd(i) ji=ji+3 *enddo allsel,all ji= i= max1= nod1= ojd= jd= j= ag= eg= eh= kp0= r0= yyl= line1= finish /solu solve finish /post1 pldisp,1 etable,sigi,ls,1 plls,sigi,sigi,1 /prep7??*SET,eg,2.95e5? ? !钢束弹模*SET,ag,139? ? !钢束面积*SET,eh,3.45e4? ? !混凝土的弹模*SET,r0,9345? ?!曲线钢束的半径? ?*SET,yyl,180000??!钢束的预应力大小et,1,link8? ?? ? !定于预应力钢束单元类型为link8et,2,solid95? ???!定义混凝土单