沙牌碾压混凝土拱坝温度徐变应的力仿真计算_水利工程论文.docVIP

沙牌碾压混凝土拱坝温度徐变应力仿真计算_水利工程论文 沙牌碾压混凝土拱坝温度徐变应力仿真计算_水利工程论文 摘要：根据沙牌工程混凝土徐变试验资料，按混凝土固化徐变理论，分解了沙牌碾压混凝土徐变度函数，得到了沙牌混凝土粘弹性相变形、粘性相变形的数学表达式，提出了混凝土的非线性徐变应力计算方法；根据沙牌碾压混凝土拱坝的材料参数与环境参数，模拟了混凝土的施工过程，得到了沙牌碾压混凝土拱坝的三维温度场与三维应力场的仿真计算成果；比较了混凝土线性徐变应力理论与非线性徐变应力理论下拱冠剖面不同高程、不同部位大坝混凝土应力随时间的变化过程，得出了一些有意义的结论，可供大坝温控设计参考。关键词：大坝仿真分析 温度应力 混凝土徐变 不可恢复徐变 　　对不设横缝或横缝间距很大的碾压混凝土拱坝，无论是在施工期，还是在运行期，温度荷载所占的比例都相当高，且具有准周期荷载的特性。在计算混凝土温度徐变应力时，应该考虑混凝土不可恢复徐变对坝体应力状态的影响。但由于混凝土不可恢复徐变的试验有一定的难度，一般的工程也不做，因此，从混凝土的已有徐变实验资料中，分离出其中的不可恢复部分，就具有重要的工程意义。Bazant固化徐变理论公式[1]是从混凝土组成的微观机制出发，根据各组成材料的物理性质推导出来的。具有概念明确、参数较少、方程线性等优良性质。文献[2]通过对沙牌工程碾压混凝土徐变资料的拟合计算表明：该公式拟合效果良好，拟合参数唯一，各参数的重要性处于同一水平。不同龄期、不同持荷时间下，老化粘弹性相徐变Ca(t,τ)、非老化粘弹性相徐变Cna(t,τ)、粘性流动相徐变Cf(t,τ)(不可复徐变)在混凝土总徐变C(t,τ)中所占的比例，与工程试验资料基本吻合，可以用于建立混凝土非线性徐变理论模型。这种考虑了不可复徐变在不同应力水平下的非线性性质的理论公式，对研究大坝混凝土温度徐变应力具有一定的优势。因为，分缝很少的大体积混凝土在温升过程中的预压应力被混凝土后期温降拉应力逐渐消解直至反超的过程，呈现出一个典型的加载又卸载的徐变应力问题，需要相应的非线性徐变理论来计算。 1　沙牌碾压混凝土徐变试验资料及其分解 　　按照Bazant固化徐变理论公式[1]，混凝土徐变度函数C(t,τ)可以分解为：C(t,τ)=Ca(t,τ)+Cna(t,τ)+Cf(t,τ)(1) 其中：Ca(t,τ)=q2Γ(t,τ)(2) Cna(t,τ)=q3ln[1+（t-τ/λ0）](3) Cf(t,τ)=q4ln（t/τ）(4) 　表1 “沙牌工程”碾压混凝土徐变度计算值与试验值 单位：10-6MPa-1 加荷龄期τ/d t-τ/d372890180 试 验 值3 7 14 30 60 90 180 36067 80 89 99 105 109 115 120 40 49 55 62 67 69 74 7824 29 33 37 42 45 48 5313 17 20 24 26 28 31 359 12 15 18 20 22 24 27 计 算 值 3 7 14 30 60 90 180 36066(71) 81(79) 86(85) 90(92) 97(99) 105(103) 118(111) 123(117)43(47) 53(53) 56(58) 58(64) 64(71) 67(74) 76(81) 85(88)23(25) 30(28) 31(31) 33(35) 36(39) 38(42) 44(48) 50(54)16(16) 20(18) 21(19) 23(21) 26(24) 28(26) 33(30) 38(35)13(13) 16(14) 18(15) 19(16) 22(18) 24(20) 29(22) 34(26) (5) 　　τ为混凝土的加载龄期，t-τ为混凝土的持荷时间；λ0、m、n为经验系数；q2、q3、q4为对具体工程试验数据进行拟合时的拟合系数。对于沙牌工程，其拟合结果为[2]： 　　q2=133.23, q3=5.44, q4=7.98，变异系数ωopt＝0.065. 　　沙牌碾压混凝土徐变度试验值与按式(1)得到的计算值列于表1中。为了和现行规范比较，表1的括号中还给出了按朱伯芳公式[3]得到的沙牌碾压混凝土徐变度计算值。由表1可见：二者的拟合效果都相当好。按照公式(2)、(3)、(4)分解式(1)得到的老化粘弹性项徐变Ca(t,τ)、非老化粘弹性项徐变Cna(t,τ)Q、混凝土的不可复徐变，即粘性流动项徐变Cf(t,τ)见表2.从表2中可以得出如下结论：表2 沙牌碾压混凝土各种徐变百分率随龄期与持荷时间变化规律 混凝土持荷 时间t-τ/d 加荷龄期τ/dCa(t,τ) C(t,τ) 3 7 14 30 60 9