ANSYS混凝土收缩徐变.docx

1、ANSYS混凝土收缩徐变ANSYS和MIDAS混凝土变模拟比较简述：本文主要对比ANSYS和MIDAS这两种有限元软件在模拟混凝土收缩变上的差异，包括计算精度、计算式、计算时间等面。计算模型为10m长的C50形柱顶施加1kN的集中力，柱截面为1m1m。1. 混凝土变混凝土变是混凝土结构在长期荷载作用下随着时间的增长混凝土中产生的应变变化目前尚未对混凝土变有比较统一的说法，在此不去讨论具体有说法，关键在于理解混凝土变与应力是有关系的。而通常我们计算结构时大部分是按照线性变处理的。2. 混凝土变本构关系2.1 老化理论本构关系根据迪辛格尔法可知变函数可定义为在t0时刻作用于混凝土的单位应力(即t0

2、=1)至时刻所产生的总应变。如采用变系数的第一种定义，则可表示为：如采用第二种定义，则可表示为：3. ANSYS立柱计算模型由于ANSYS并没有专门板块来混凝土变模拟，故而需要借助金属蠕变的计算机理来等效模拟混凝土变效应。ANSYS提供两种法计算变：显式计算和隐式计算。显式计算需要细分较多的时间步长，计算时间长；隐式计算计算精度高，计算时间短。但是在实践中也发现，涉及到单元生死情况时，隐式计算可能出现异常现象。下面将会对这两种法进行详细的比较。3.1 ANSYS显式计算显式计算对时间步长是有要求的，尤其是在变系数曲线变化剧烈的时间段需要细分子步以减小误差和帮助收敛。因而，时间步长的划分式、时间

3、点的数目对计算结果都会有较大的影响。（1） 等间距时间步长和对数时间步长假设混凝土的龄期是7天，变变化速率为0.005，考虑收缩变10年（3650天），若3650天时刻的变系数为1，那么按照等间距时间步长划分，则时间步长间距，（3650-7）/500=7.286。按照对数时间步长划分，若采用30个数据点，具体数据如下所示。表1 对数时间步长数据表编号时间编号时间编号时间编号时间编号时间171124.464542185.5019431298.8236411044.36927.9330821227.72562296.8991332338.656421183.5838.9905411331.4213

4、523109.815533383.798431341.349410.188961435.6097324124.453734434.9573441520.147511.547121540.3564225141.04335492.936451722.779613.086321645.7358226.843736558.6431461952.421714.830691751.8322927.150437633.1087472212.673816.807581858.741428205.297338717.5005482507.617919.047991966.5714729232.662939813

5、.1415492841.8761021.587042075.4452830263.676340921.5311503650（2） 变变化速率变变化速率影响着变变化曲线的陡缓，将会对不同变变化速率值进行比较：0.001、0.002、0.003、0.004、0.005、0.006、0.007、0.008、0.009、0.。（3） 计算结果对比两种时间步长划分式和不同变变化速率柱顶变10年位移计算结果如下表所示：表2: 柱顶变10年计算结果对比（单位:m）变变化速率理论值等间距时间步长对数时间步长计算值相对误差计算值相对误差0.0015.721e-75.71e-70.19%5.66e-71.07%0

6、.0025.795e-75.78e-70.26%5.74e-70.95%0.0035.797e-75.77e-70.47%5.75e-70.81%0.0045.797e-75.77e-70.47%5.75e-70.81%0.0055.797e-75.76e-70.64%5.74e-70.98%0.0065.797e-75.75e-70.81%5.74e-70.98%0.0075.797e-75.74e-70.98%5.74e-70.98%0.0085.797e-75.74e-70.98%5.74e-70.98%0.0095.797e-75.73e-71.16%5.74e-70.98%0.5.7

7、97e-75.72e-71.33%5.74e-70.98%从上表可以看出，不同的时间步划分式对结果产生较大的差异，等距时间步长随着变变化速率的增大精度不断降低，对数时间步长则随着变变化速率的增大精度先上升后趋于平稳。且对数时间步长的计算时间要短，精度也能满足工程要求，且在较大变变化速率区间，采用对数时间步长更容易获得较好结果。故而，建议采用对数时间步长进行显式计算。3.2 ANSYS隐式计算隐式计算也需要区分等间距时间步长和对数时间步长两种划分式。不过这里不对隐式计算进行详细的探讨，光是一个简单的立柱还不能够说清楚ANSYS模拟变问题。下面将会讨论考虑施工过程的两端固结梁的变问题。4. ANS

8、YS两端固结梁计算模型只考虑显式计算。仍然考虑等间距划分和对数划分两种式，比较这两种划分式的计算结果、收敛情况、计算耗时等。此处附上显式计算命令流，不作详细说明，因为显式计算碰到收敛问题无法解决！! 计算悬臂梁转化为固定梁的变效应finish/clear/prep7Ec=3.45e10! time array*dim,tt,array,50,1*vread,tt(1),tt1,txt(1F8.3)! creep coefficient=0.005nn=50*dim,fi,array,nn,1*do,i,1,nnfi(i)=1-2.7182818*(-0.005*(tt(i)-7)*enddo*

9、dim,C1,array,nn,1f1=0k1=3dt=0*do,i,1,nnC1(i)=(fi(i)-f1)/(1+fi(i)/(tt(i)-dt)f1=fi(i)dt=tt(i)mp,ex,k1,Ecmp,prxy,k1,0.2tb,creep,k1tbdata,1,C1(i),0,1,0k1=k1+2*enddo*dim,C2,array,nn,1f1=0k2=4dt=0*do,i,1,nnC2(i)=2*(fi(i)-f1)/(1+2*fi(i)/(tt(i)-dt)f1=fi(i)dt=tt(i)mp,ex,k2,Ecmp,prxy,k2,0.2tb,creep,k2tbdata,1

10、,C2(i),0,1,0k2=k2+2*enddoet,1,23keyopt,1,6,0r,1,1,1/12,1! 左边材料属性mp,ex,1,Ecmp,prxy,1,0.2tb,creep,1tbdata,1,0,0,0,0! 右边材料属性mp,ex,2,Ecmp,prxy,2,0.2tb,creep,2tbdata,1,0,0,0,0! 建立有限元模型*do,i,1,19n,i,(i-1)/2*enddo*do,i,1,19n,i+19,(i-1)/2+9*enddotype,1mat,1*do,i,1,18e,i,i+1*enddotype,1mat,2*do,i,18,35e,i+2,

11、i+2+1*enddod,1,alld,38,allesel,s,1,18sfbeam,all,pres,1000esel,s,19,36sfbeam,all,pres,1000/soluoutres,all,allallsel,allcrplim,0.25bfunif,temp,100time,1e-6solvek1=3k2=4*do,i,1,nnesel,s,1,18mpchg,k1,allk1=k1+2allsel,allesel,s,19,36mpchg,k2,allk2=k2+2allsel,alltime,tt(i)nsubst,10solve*enddo隐式计算命令流如下：! 计

12、算悬臂梁转化为固定梁的变效应finish/clear/prep7Ec=3.45e10! time array*dim,tt,array,50,1*vread,tt(1),tt1,txt(1F8.3)! creep coefficient=0.005nn=50*dim,fi,array,nn,1*do,i,1,nnfi(i)=1-2.7182818*(-0.005*(tt(i)-7)*enddo*dim,C1,array,nn,1f1=0k1=3dt=0*do,i,1,nnC1(i)=(fi(i)-f1)/Ec/(tt(i)-dt)f1=fi(i)dt=tt(i)mp,ex,k1,Ecmp,pr

13、xy,k1,0.2tb,creep,k1,11tbdata,1,C1(i),1,0,0,0k1=k1+2*enddo*dim,C2,array,nn,1f1=0k2=4dt=0*do,i,1,nnC2(i)=2*(fi(i)-f1)/Ec/(tt(i)-dt)f1=fi(i)dt=tt(i)mp,ex,k2,Ecmp,prxy,k2,0.2tb,creep,k2,11tbdata,1,C1(i),1,0,0,0k2=k2+2*enddoet,1,sectype,1,beam,rect,a1secdata,1,1keyopt,1,1,0keyopt,1,3,3! 左边材料属性mp,ex,1,Ec

14、mp,prxy,1,0.2tb,creep,1,11tbdata,1,0,0,0,0,0,0! 右边材料属性mp,ex,2,Ecmp,prxy,2,0.2tb,creep,2,11tbdata,1,0,0,0,0,0,0! 建立有限元模型*do,i,1,19n,i,(i-1)/2*enddo*do,i,1,19n,i+19,(i-1)/2+9*enddotype,1mat,1*do,i,1,18e,i,i+1*enddotype,1mat,2*do,i,18,35e,i+2,i+2+1*enddod,1,alld,38,allesel,s,1,18sfbeam,all,pres,1000esel,s,19,36sfbeam,all,pres,1000/soluoutres,all,allrate,offallsel,allcrplim,0.25bfunif,temp,100time,1e-6solvecpintf,allrate,onk1=3k2=4*do,i,1,nnesel,s,1,18mpchg,k1,allk1=k1+2allsel,allesel,s,19,36mpchg,k2,allk2=k2+2allsel,alltime,tt(i)nsubst,40solve*enddo5. MIDAS计算模型