《预应力建模(曲线预应力索和直线预应力索)》由会员分享，可在线阅读，更多相关《预应力建模(曲线预应力索和直线预应力索)（10页珍藏版）》请在金锄头文库上搜索。
1、2 独立建模耦合法该法的基本思想是实体和力筋独立建几何模型，分别划分单元，然后采用耦合方程将力筋单 元和实体单元联系起来，这种方法是基于有限元模型的处理。其基本步骤如下： 建立实体几何模型（不考虑力筋）； 建立力筋线的几何模型（不考虑体的存在）； 将几何模型按一定的要求划分单元（这时也是各自独立的）； 选择所有力筋线； 选择与上述力筋相关的节点（nsll命令），并定义选择集； 将上述力筋节点存入数组； 选择所有节点，并去掉中的节点集（这时是除力筋节点外的所有节点）； 按力筋节点数组搜寻所有最近的实体节点号，并存入数组中； 耦合力筋节点与最近的节点，一一耦合（cp命令）（不能使用cpintf命令
2、，这样可能耦合其 它节点，且容易不耦合） 选择所有，并施加边界条件和荷载，可以求解了。这种方法建模特别简单，耦合处理也比较简单（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.
3、2 blc4,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,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,uylsel,s,loc,z,3000 lsel,r,loc,y,0 dl,all,all allsel,a
4、ll !单元划分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,
5、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,allcmsel,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=
6、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!定于预应力钢束单元类型为 link8e
7、t,2,solid95!定义混凝土单元类型为 solid95r,1,ag,yyl/eg/ag*1.1194364!定义单元实常数，对于Iink8单元为钢束面积、以及钢束的初始应 变（除模拟预应力张拉还会用于应力刚度矩阵的计算）r,2!solid95的单元实常数仅为x轴的调整，并且该实常数仅仅在单元坐标系统选择为 1 的时候填写。一般不填。mp,ex,1,eg!预应力钢束单元性质，包括之前定义的弹模，泊松比mp,prxy,1,0.3mp,ex,2,eh!混凝土单元性质，包括弹模、泊松比mp,prxy,2,0.2MPTEMP,!钢束和混凝土的密度MPTEMP,1,0MPDATA,DENS,1,7.
8、921e-6MPTEMP,MPTEMP,1,0MPDATA,DENS,2,2.3e-6blc4,100,200,3000 !定义长方体的长宽高（混凝土块体）/view,1,1,1,1!调整观察的角度/ang,1vplot!显示选择的体ksel,all!选择所有的关键点*get,kpO,kp,O,num,max !将最大的kp编号存于kpO中，因为事先不知道最大的关键点号是多 少，所以先要存在kpO中lsel,none!不选择所有的线k,kp0+1,50,160!定义关键点，kpO是已有的最大关键点编号。k,kpO+2,5O,16O,3OOOk,kpO+3,5O,8OO,15OOlarc,kpO
9、+1,kpO+2,kpO+3,rO!将定义的以上三点连成一条圆弧，半径为 rO。 kpO+3 用于定义圆弧所在面及曲线边中点kdeIe,kpO+3!删除kpO+3这个关键点*get,Iine1,Iine,O,num,min !将最小的线编号存于linel中，因为linel是最先定义的预应力 钢束lsel,s,loc,z,0选择z、y坐标均为0的线，用于施加约束（该线为右侧截面的梁底线）lsel,r,loc,y,Odl,all,uy lsel,s,loc,z,3000 lsel,r,loc,y,0 dl,all,allallsel,all lsel,s,lineilatt,i,i,i lesiz
10、e,all,50lmesh,all vsel,all!对线施加y方向的约束!选择左侧截面梁底线（即z=3000，而y=0的线）! 对选择的线施加约束（ UX 、UY 、 UZ ） !选择所有对象!选择线，最小值 linei!将单元性质与选择的、尚未网格化的线联系起来 !线单元划分的数目!对linei进行分网! 选择所有的体vatt,2,2,2 lsel,s,loc,z,0 lsel,r,loc,y,i0,i40 lesize,all,8 lsel,s,loc,z,0 lsel,u,loc,y,i0,i40 lesize,all,4 lsel,s,loc,y,0 分网的扫掠线。 lsel,r,l
11、oc,x,0 lesize,all,50 vsweep,all allsel,all lsel,s,linei nsll,s,i!给体赋予相应的属性（混凝土）选择z=0,而y坐标从10到140范围内的所有的线!分网数目 8!再次选择z=0,而不选y坐标从10到140范围内的所有线! 分网数目为 4!选择x=0、y=0坐标处的线（仅有一根线被选中）注意这根线是一根!分网数目为 50!通过扫掠网格的方法来用单元填充体!选择所有的对象!选择 line1 最小值（预应力钢绞线）!选择与所选择的线对应的所有的节点，即钢绞线上面的节点，因为网格数为 50，所以节点总数为51cm,cmljnod,node!
12、将钢绞线节点编组，其组名为 cmljnod*get,max1,node,0,count!选择钢绞线中的节点单元数目，将其赋给变量maxi,此变量可以作*do循环的上界来确定循环的次数*dim,ojd,max1*dim,jd,max1!定义数组序列。数组名为 ojd， max1 为行的范围!定义数组序列。数组名为 jd， max1 为行的范围*get,nod1,node,0,num,min!选择集中最小的节点号，赋给变量nodi,作为*do循环的上界*SET,ojd(1),nod1*do,i,2,max1!将nod1的值赋给ojd（1）!do循环，从2到max1*SET,ojd(i),ndnex
13、t(ojd(i-1) !将 ndnext(ojd(i-1)赋给 ojd(i)*enddo allsel,all nsel,all cmsel,u,cmljnod *do,i,1,max1 *SET,nod1,ojd(i) *SET,j,nnear(nod1)!结束循环!选择所有!选择所有节点!不选择cmljnod，即钢绞线对应的节点!i从1循环开始到max1!将 ojd（i）赋给 nod1!将nnear（nod1）赋给j,是指ojdi邻近的节点吗？ nnear（）是一个*SET,jd(i),j *enddo nsel,all *SET,ji,1 *do,i,1,max1 cp,ji,ux,oj
14、d(i),jd(i) 方向的位移。 cp,ji+1,uy,ojd(i),jd(i) cp,ji+2,uz,ojd(i),jd(i) *SET,ji,ji+3 *enddo allsel,all *SET,ji, *SET,i, *SET,max1, *SET,nod1, *SET,ojd, *SET,jd, *SET,j, *SET,ag, *SET,eg, *SET,eh, *SET,kp0, *SET,r0, *SET,yyl, *SET,line1,!将 j 赋给 jd(i)!结束循环! 选择所有节点!将1赋给ji!i 从1 到 max1!定义耦合自由度，ji只是随机设置的耦合编号。ux
15、为耦合节点的x!ji+3 赋给 ji! 结束循环! 选择所有单元!将所有的参数删除(第三项value为空格时标识删除)finish/solu! 建模结束，进入求解allsel,all acel,9.8 solve! 考虑自重函数，表示选择最近的节点。!进入后处理FINISH/POST1pldisp,1 etable,sigi,ls,1 plls,sigi,sigi,1 plnsol,s,1建立直线型tender的预应力梁/prep7*SET,eg,2.95e5 !钢束弹模*SET,ag,139 !钢束面积*SET,eh,3.45e4 !混凝土的弹模*SET,r0,9345 !曲线钢束的半径*SET,yyl,140000 !钢束的预应力大小et,1,link8!定于预应力钢束单元类型为 link8et,2,solid95 !定义混
