众所周知,在ANSYS中,预应力混凝土分析(有粘结)可采用等效荷载法和实体力筋法。所谓等效荷载法,就是将力筋的作用以荷载的形式作用于混凝土结构;所谓实体力筋法就是用solid模拟混凝土,而link模拟力筋。 1 等效荷载法的优缺点 优点是建模简单,不必考虑力筋的具体位置而可直接建模,网格划分简单;对结构的在预应力作用下的整体效应比较容易求得。 其主要缺点是: ①等效荷载法没有考虑力筋对混凝土的作用分布和方向,力筋对混凝土作用显然在各处是不同的,等效荷载法则无法考虑;水平均布分量没有考虑。 ②对某些线形的力筋模拟困难,例如通常采用的是直线(较短)+曲线+直线(很长)+曲线+直线(较短),这种形式的布筋等效起来麻烦,且可能不合理。 ③难以求得结构细部受力反映,否则荷载必须施加在力筋的位置上,这又失去建模的方便性。 ④在外荷载作用下的共同作用难以考虑,不能确定力筋在外荷载作用下的应力增量。 ⑤对张拉过程无法模拟。 ⑥无法模拟应力损失引起的力筋各处应力不等的因素。 其最大的一个缺点是:较粗!得到的结果与实际情况误差较大!最近做了点实际计算,经过比较发现,结果与实际的误差相差较多(可能是特例),所以采用该方法需要谨慎和校验一下。 2 实体力筋法的优缺点 拜年帖2-----预应力混凝土分析中实体力筋法的ansys处理过程 1 体分割法 2 独立建模耦合法 ①建立实体几何模型(不考虑力筋); ②建立力筋线的几何模型(不考虑体的存在); ③将几何模型按一定的要求划分单元(这时也是各自独立的); ④选择所有力筋线; ⑤选择与上述力筋相关的节点(nsll命令),并定义选择集; ⑥将上述力筋节点存入数组; ⑦选择所有节点,并去掉⑤中的节点集(这时是除力筋节点外的所有节点); ⑧按力筋节点数组搜寻所有最近的实体节点号,并存入数组中; ⑨耦合力筋节点与最近的节点,一一耦合(cp命令)(不能使用cpintf命令,这样可能耦合其它节点,且容易不耦合) ⑩选择所有,并施加边界条件和荷载,可以求解了。 预应力简支梁弹性分析--体线独立耦合法示例 !----------------------------------------- /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 模拟预应力主要分为两大块:一.模拟预应力的张拉过程;二.模拟预应力在结构中的作用. 对于一比较复杂,需要考虑摩擦,主要问题是参数的选取,而不是如何建立模型的问题,具体参数需要大量的试验才能确定. 对于二是结构分析比较关心的问题,不论采用耦合、共用节点,约束方程,主要问题是力的传递路径是否正确,以及预应力的分布是否正确。对于体外索,整根预应力束应力相同,可以采用统一的初始应变,或着降温,中间节点在和转向器连接的地方放松纵向自由度,径向同转向器节点耦合。而对于体内束,不管先张、后张,结构形成以后预应力和混凝土已经固节(除非你分析极限承载力,考虑滑移),预应力节点和混凝土节点应该完全耦合。其主要问题是整个预应力束的预应力分布如何模拟,施加分段初应变是可取的,但是比较繁琐(但完全可以很好的控制)。 还有一个问题,就是初应变的大小如何确定,实际上,确定初应变的大小就是模拟张拉的一个过程,我们最终要的是存留值,要想精确模拟所以必须通过反复调整初应变来模拟混凝土梁的弹性压缩损失。 还有“有限元梁杆组合结构分析方法” 个人愚见: 所以很多分析都是针对跨中数据的检验,好像正确(跨中刚好对称),若对离开跨中一段距离进行分析,就会发现不正确。所以若要对梁体进行局部应力分布分析,要考虑滑移和摩察的切向影响后,就会发现问题。 所以实际用ANSYS计算的预应力效果不如用结构设计原理和规范计算结果准确。桥梁博士是采用设计原理和规范计算方法,所以比较准确。 正在做连续梁桥的悬臂施工仿真分析,对于三向预应力钢筋的模拟是一个很关键的问题,总结了一些自己的体会. 1;用面面切割体生成线是很方便也是比较理想的,能准确的模拟力筋的位置,对线型把握得很好,但要考虑一个问题就是在力筋线很多,而模型本身又不是很规整的情况下会造成切割成太多的小实体(比如横向,竖向的钢筋靠得很近而且又多),在分网的时候就会出现问题以致没法做下去. 2:等效荷载法是一种很传统的方法,很多桥梁分析软件在考虑预应力的时候是在平面杆系结构中用等效荷载来做的..我觉得在做整体分析是用它来分析应该是可以的.要做局部分析的话,不是很好;比如在做箱梁分析的时候,纵向的预应力钢筋长而且线型复杂.用这种方法是不行的.而在考虑竖向钢筋的时候,直筋而且不是很长,(LMAX=4M)采取等效荷载法来做. 3: 当预应力筋的线型比较复杂,实体和预应力筋单独建模,分网,用约束方程法来做不错,我在考虑腹板的纵向预应力筋的时候就是这样来做的 4;在同一个模型中结合上面的三种方法来做,说来也是没办法啊,在用方法一实在是做不下去了,就用综合后面的方法来做,不过最后还算可以. 5;预应力的损失如何算,就按规范来做吗,可不可以用ANSYS来做这方面的分析,希望有更多的人来探讨 6.一点拙见,希望各位指教. 在ansys中施加预应力好像是个蛮重要的问题,希望大伙能仁者见仁智者见智。 我先说两句: 1.在构件表面可以加等效荷载,现在好像一般都采用这种方法。 2.把预应力筋的位置从混凝土体中挖去,然后在在洞口的表面上加上荷载,适于精确分析。这种方法我正在尝试,但比较繁,而且网格划分会很密,不适宜做整体的分析,但可用于局部分析。 3.加温度荷载。在预应力筋处降温,产生收缩,施加当量的预应力,这种方法目前正在尝试。请高手多多指教! 在一篇论文中有这样的描述:墙板,环梁,底板均采用ANSYS提供的三维实体SOLID65进行分析,预应力钢绞线采用连接元LINK8,此单元有拉伸压缩刚度,无抗弯刚度。划分单元时候,使LINK8元的节点和SOLID65元的节点置空间同一点。空间同点上的二节点在线单元LINK8垂直平面内的位移,利用约束方程处理为相同,而二节点沿线元纵向的位移各自独立。LINK8元和SOLID65元空间同点节点的约束方程处理,正确反映了后张预应力施工时候,钢绞线只能沿墙纵向相对滑移而不会走离墙中心位置饿真实过程。 It depends what do you want. If you want to determine the stress distribution at EACH stage you apply the loads, I advice you use link element with initial strain. If you want to decide the stress distributon ON If you want to design the structure at ultimate limit stage, maybe you can even igonre the effect of prestress, because all the reinforcement have reached their strength. (Of course you can include the seconde order moments and shears). 用温度应力代替预应力是可行的. 可能也是最好的模拟预应力方法. 可以做几个习题试一下 在ANSYS中通过变温施加预应力是一种简便实用的手段。先将预应力筋的形状用LINK单元表示好(转折处用约束方程规定),不必担心单元的重叠,然后根据需要定义LINK单元的线胀系数及其温变(其它单元不要规定温变)。 尝试计算一次,核对一下LINK单元的应力是否满足精度(因为相连的其它单元会一起变形),调整一两次就可以了 还有一个办法,就是直接加预应力荷载. 在solution中,有一个apply initial stress 的命令.不过好像只能用于部分单元:link180,beam188,beam189,etc... 预应力混凝土的分析方法可分为两大类:其一是将力筋的作用以荷载的形式作用于结构即所谓的等效荷载法,其二是力筋和混凝土分别用相应的单元模拟,预应力通过不同的模拟方法施加称之为实体力筋法,这两种方法都可根据不同的分析目的或需要而采用不同的单元进行模拟。 实体力筋 |
|