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

2 独立建模耦合法该法的基本思想是实体和力筋独立建几何模型，分别划分单元，然后采用耦合方程将力筋单元和实体单元联系起来，这种方法是基于有限元模型的处理。其基本步骤如下：①建立实体几何模型（不考虑力筋）；②建立力筋线的几何模型（不考虑体的存在）；③将几何模型按一定的要求划分单元（这时也是各自独立的）；④选择所有力筋线；⑤选择与上述力筋相关的节点(nsll命令)，并定义选择集；⑥将上述力筋节点存入数组；⑦选择所有节点，并去掉⑤中的节点集(这时是除力筋节点外的所有节点)；⑧按力筋节点数组搜寻所有最近的实体节点号，并存入数组中；⑨耦合力筋节点与最近的节点，一一耦合(cp命令)（不能使用cpintf命令，这样可能耦合其它节点，且容易不耦合）⑩选择所有，并施加边界条件和荷载，可以求解了。这种方法建模特别简单，耦合处理也比较简单(APDL要熟悉些)，缺点是当实体单元划分不够密时，力筋节点位置可能有些走动，但误差在可接受范围之内！这种方法是解决力筋线形复杂且力筋数量很多时的较佳方法。预应力简支梁弹性分析--体线独立耦合法示例!-----------------------------------------/prep7eg=2e5ag=140eh=4e4r0=9345yyl=200000et,1,link8et,2,solid95r,1,ag,yyl/eg/ag*1.036258r,2mp,ex,1,egmp,prxy,1,0.3mp,ex,2,ehmp,prxy,2,0.2blc4,,,100,200,3000/view,1,1,1,1/ang,1vplot!------------定义力筋线ksel,all*get,kp0,kp,0,num,maxlsel,nonek,kp0+1,50,160k,kp0+2,50,160,3000k,kp0+3,50,800,1500larc,kp0+1,kp0+2,kp0+3,r0kdele,kp0+3*get,line1,line,0,num,min!-------------定义约束lsel,s,loc,z,0lsel,r,loc,y,0dl,all,,uylsel,s,loc,z,3000lsel,r,loc,y,0dl,all,,allallsel,all!-----------单元划分lsel,s,,,line1latt,1,1,1lesize,all,,,50lmesh,allvsel,allvatt,2,2,2lsel,s,loc,z,0lsel,r,loc,y,10,140lesize,all,,,8lsel,s,loc,z,0lsel,u,loc,y,10,140lesize,all,,,4lsel,s,loc,y,0lsel,r,loc,x,0lesize,all,,,50vsweep,allallsel,all!耦合自由度lsel,s,,,line1nsll,s,1cm,cmljnod,node*get,max1,node,0,count*dim,ojd,,max1*dim,jd,,max1*get,nod1,node,0,num,minojd(1)=nod1*do,i,2,max1ojd(i)=ndnext(ojd(i-1))*enddoallsel,allnsel,allcmsel,u,cmljnod*do,i,1,max1nod1=ojd(i)j=nnear(nod1)jd(i)=j*enddonsel,allji=1*do,i,1,max1cp,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*enddoallsel,allji=i=max1=nod1=ojd=jd=j=ag=eg=eh=kp0=r0=yyl=line1=finish/solusolvefinish/post1pldisp,1etable,sigi,ls,1plls,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? ???!定义混凝土单元类型为solid95r,1,ag,yyl/eg/ag*1.1194364!定义单元实常数，对于link8单元为钢束面积、以及钢束的初始应变（除模拟预应力张拉还会用于应力刚度矩阵的计算）? ?r,2? ?? ?? ?? ???!solid95的单元实常数仅为x轴的