ANSYS混凝土收缩徐变.docVIP

WORD格式可编辑 专业知识整理分享 ANSYS和MIDAS混凝土徐变模拟比较 简述：本文主要对比ANSYS和MIDAS这两种有限元软件在模拟混凝土收缩徐变上的差异，包括计算精度、计算方式、计算时间等方面。计算模型为10m长的C50方形柱顶施加1kN的集中力，柱截面为1m×1m。 混凝土徐变 混凝土徐变是混凝土结构在长期荷载作用下随着时间的增长混凝土中产生的应变变化 目前尚未对混凝土徐变有比较统一的说法，在此不去讨论具体有何说法，关键在于理解混凝土徐变与应力是有关系的。而通常我们计算结构时大部分是按照线性徐变处理的。 混凝土徐变本构关系 老化理论本构关系 根据迪辛格尔法可知徐变函数可定义为在t0时刻作用于混凝土的单位应力(即σt0=1)至时刻ｔ所产生的总应变。如采用徐变系数φ( 如采用第二种定义，则可表示为： ANSYS立柱计算模型 由于ANSYS并没有专门板块来混凝土徐变模拟，故而需要借助金属蠕变的计算机理来 等效模拟混凝土徐变效应。ANSYS提供两种方法计算徐变：显式计算和隐式计算。显式计算需要细分较多的时间步长，计算时间长；隐式计算计算精度高，计算时间短。但是在实践中也发现，涉及到单元生死情况时，隐式计算可能出现异常现象。下面将会对这两种方法进行详细的比较。 ANSYS显式计算 显式计算对时间步长是有要求的，尤其是在徐变系数曲线变化剧烈的时间段需要细分子步以减小误差和帮助收敛。因而，时间步长的划分方式、时间点的数目对计算结果都会有较大的影响。 等间距时间步长和对数时间步长 假设混凝土的龄期是7天，徐变变化速率为0.005，考虑收缩徐变10年（3650天），若3650天时刻的徐变系数为1，那么按照等间距时间步长划分，则时间步长间距，（3650-7）/500=7.286。按照对数时间步长划分，若采用30个数据点，具体数据如下所示。 表 SEQ 表 \* ARABIC 1 对数时间步长数据表 编号 时间 编号 时间 编号 时间 编号 时间 编号 时间 1 7 11 24.46454 21 85.50194 31 298.8236 41 1044.369 2 7.933082 12 27.7256 22 96.89913 32 338.656 42 1183.58 3 8.990541 13 31.42135 23 109.8155 33 383.798 43 1341.349 4 10.18896 14 35.60973 24 124.4537 34 434.9573 44 1520.147 5 11.54712 15 40.35642 25 141.043 35 492.936 45 1722.779 6 13.08632 16 45.73582 26 159.8437 36 558.6431 46 1952.421 7 14.83069 17 51.83229 27 181.1504 37 633.1087 47 2212.673 8 16.80758 18 58.7414 28 205.2973 38 717.5005 48 2507.617 9 19.04799 19 66.57147 29 232.6629 39 813.1415 49 2841.876 10 21.58704 20 75.44528 30 263.6763 40 921.5311 50 3650 徐变变化速率 徐变变化速率影响着徐变变化曲线的陡缓，将会对不同徐变变化速率值进行比较：0.001、0.002、0.003、0.004、0.005、0.006、0.007、0.008、0.009、0.010。 计算结果对比 两种时间步长划分方式和不同徐变变化速率柱顶徐变10年位移计算结果如下表所示： 表 SEQ 表 \* ARABIC 2: 柱顶徐变10年计算结果对比（单位:m） 徐变变化速率 理论值 等间距时间步长 对数时间步长 计算值 相对误差 计算值 相对误差 0.001 5.721e-7 5.71e-7 0.19% 5.66e-7 1.07% 0.002 5.795e-7 5.78e-7 0.26% 5.74e-7 0.95% 0.003 5.797e-7 5.77e-7 0.47% 5.75e-7 0.81% 0.004 5.797e-7 5.77e-7 0.47% 5.75e-7 0.81% 0.005 5.797e-7 5.76e-7 0.64% 5.74e-7 0.98% 0.006 5.797e-7 5.75e-7 0.81% 5.74e-7 0.98% 0.007 5.797e-7 5.74e-7 0.98% 5.74e-7 0.98% 0.0