ANSYS混凝土收缩徐变.doc

(完整word版)ANSYS混凝土收缩徐变(完整word版)ANSYS混凝土收缩徐变(完整word版)ANSYS混凝土收缩徐变ANSYS和MIDAS混凝土徐变模拟比较简述:本文主要对比ANSYS和MIDAS这两种有限元软件在模拟混凝土收缩徐变上的差异，包括计算精度、计算方式、计算时间等方面。计算模型为10m长的C50方形柱顶施加1kN的集中力，柱截面为1m1m.混凝土徐变混凝土徐变是混凝土结构在长期荷载作用下随着时间的增长混凝土中产生的应变变化目前尚未对混凝土徐变有比较统一的说法,在此不去讨论具体有何说法，关键在于理解混凝土徐变与应力是有关系的.而通常我们计算结构时大部分是按照线性徐变处理的.混凝土徐变本构关系老化理论本构关系根据迪辛格尔法可知徐变函数可定义为在t0时刻作用于混凝土的单位应力（即t0=1）至时刻ｔ所产生的总应变。如采用徐变系数的第一种定义，则可表示为：如采用第二种定义，则可表示为：ANSYS立柱计算模型由于ANSYS并没有专门板块来混凝土徐变模拟，故而需要借助金属蠕变的计算机理来等效模拟混凝土徐变效应。ANSYS提供两种方法计算徐变:显式计算和隐式计算。显式计算需要细分较多的时间步长，计算时间长;隐式计算计算精度高,计算时间短。但是在实践中也发现，涉及到单元生死情况时，隐式计算可能出现异常现象。下面将会对这两种方法进行详细的比较。ANSYS显式计算显式计算对时间步长是有要求的，尤其是在徐变系数曲线变化剧烈的时间段需要细分子步以减小误差和帮助收敛。因而，时间步长的划分方式、时间点的数目对计算结果都会有较大的影响.等间距时间步长和对数时间步长假设混凝土的龄期是7天,徐变变化速率为0。005,考虑收缩徐变10年（3650天），若3650天时刻的徐变系数为1，那么按照等间距时间步长划分,则时间步长间距,（3650—7）/500=7.286。按照对数时间步长划分，若采用30个数据点，具体数据如下所示.表SEQ表\*ARABIC1对数时间步长数据表编号时间编号时间编号时间编号时间编号时间171124。464542185.5019431298。8236411044。36927。9330821227.72562296.8991332338.656421183。5838.9905411331.4213523109。815533383。798431341.349410.188961435。6097324124。453734434。9573441520.147511。547121540。3564225141。04335492.936451722。779613.086321645。7358226159。843736558.6431461952.421714.830691751.8322927181.150437633。1087472212.673816.807581858。741428205。297338717。5005482507.617919。047991966.5714729232。662939813。1415492841.8761021。587042075。4452830263.676340921.5311503650徐变变化速率徐变变化速率影响着徐变变化曲线的陡缓，将会对不同徐变变化速率值进行比较：0。001、0.002、0。003、0。004、0。005、0.006、0.007、0。008、0.009、0.010.计算结果对比两种时间步长划分方式和不同徐变变化速率柱顶徐变10年位移计算结果如下表所示：表SEQ表\*ARABIC2：柱顶徐变10年计算结果对比(单位：m）徐变变化速率理论值等间距时间步长对数时间步长计算值相对误差计算值相对误差0.0015.721e-75.71e-70.19%5.66e-71.07％0.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.0105.797e—75。72e-71。33%5。74e-70.98%从上表可以看出，不同的时间步划分方式对结果产生较大的差异,等距时间步长随着徐变变化速率的增大精度不断降低，对数时间步长则随着徐变变化速率的增大精度先上