附录四坝体混凝土温度和温度应力计算

1、附录四 坝体混凝土温度和温度应力计算 、温度计算 1.坝体混凝土的初期温度计算(有内热源的温度场计算) (1)计算目的： 坝体混凝土的初期温度计算目的，主要是确定基础块混凝土(或靠近老混凝土块的混凝土)中的最高温度T，以便控制基础温差，最高温度T可按下式计算： (附79)式中 Tj混凝土的浇筑温度(，以下均同)，或称入仓温度； Tr混凝土因水化热和其他原因产生的最高温升。 (2)混凝土的浇筑温度计算： (附80)式中 To混凝土的拌和(即出机口)温度(忽略拌和中的热量损失或热量流入影响)； Tq混凝土浇筑时的平均气温； 考虑混凝土在拌和、装卸、运输、转运和浇筑过程中热量损失或倒罐的系数。在一般

2、的现场条件下，=0.20.3，当运距较长，转运手续较多以及采用人工方法浇筑时，h =0.40.5。混凝土的拌和温度按下述公式计算： (附81)式中 Wi每立米混凝土中各种原材料的重量，kg/m3； Ci混凝土各种原材料的比热，kcal/(kg)； Ti混凝土各种原材料的温度。注：在公式(附81)中未考虑骨料含水率的影响，当骨料含水率较大，不宜忽略时，应在公式中加以考虑。 当在混凝土拌和中加入冰屑时，应考虑冰的潜热(80kcal/kg)和有效利用系数0.70.8。 应考虑混凝土拌和时，拌和机发出的机械热，在没有实测资料情况下，可用350kcal。 在缺乏具体资料时，各种原材料的比热C可按附表16

3、采用。附表16材料比热C20304050水石英砂普通水泥1.00.1660.1081.00.1760.1231.00.1860.1481.00.2000.180粗骨料玄武岩白云岩花岗岩石灰岩石英岩粗面岩0.1830.1920.1720.1780.1640.1830.1830.1950.1720.1810.1720.1850.1850.2000.1730.1850.1790.1890.1940.2080.1800.1920.1860.196 (3)混凝土的温升计算： 混凝土入仓后的温升Tr，主要由水化热引起，此外混凝土入仓温度Tj和气温Tq的温差；浇筑块顶面(有时顶面加侧面)和冷却水管的散热以及

4、基岩的吸热作用也对Tr有一定的影响。 对于重要的混凝土重力坝，较长的浇筑块以及在较复杂的施工条件下(预冷骨料，坝块并不均匀上升或在浇筑期中气温有较大变化者)，宜用差分法计算浇筑块的温度变化过程，直到求出浇筑块上的最高温度值，并将最高温度沿浇筑块高程画成曲线。 用差分法计算混凝土浇筑块内温度分布的基本公式为： 1)单向差分法： (附82)式中 计算点(点n)在第时段的温度(参见附图14)； 计算点在下一时段(相距)时的温度； ，与计算点上下相邻的两点(点n-1，点n+1)在第时段的温度，这些点子与计算点相距各为； a混凝土的导温系数，m2/d； 计算中所取的时间分段，d； 计算中所取的距离分段，

5、m； 在第时段()中由于水化热所产生的温升。水泥的水化热过程可以用公式或表示时，可按下式计算： (附83) (附84)以上式中 在龄期时的累积水化热，kcal/kg； 水泥最终发热量，kcal/kg； m水泥发热速率参数，d-1； n和m性质类似的常数； Th水化热绝热温升()； W每立米混凝土中的水泥用量，kg/m3； Ch混凝土的比热，kcal/(kg)。 2)双向差分法： (附85)式中 计算点在第时段的温度； 计算点在下一时段(相隔)时的温度；与计算点左右相邻两点(点m-1，点m+1)在第时段的温度，这些点子与计算点相距各为。 进行差分法计算时应注意以下各点： 参变数。在单向差分法公式

6、中要凑成1/2或1/4，在双向差分法公式中要凑成1/4或1/8，值须能将各浇筑块均匀分段。 起始温度。浇筑在基岩上的第一层混凝土，其基岩内各点的起始温度可采用浇筑月份的原始地温，在浇筑块内可采用浇筑温度Tj。 表面温度。初期由于有水泥水化热的影响，混凝土表常较气温为高，即 T为混凝土表面温度高于气温的差值，其值可采用：顶面不盖草袋时，T=35；顶面盖一层草袋时，T10；顶面流水养护时，(气温+流水温度)。 水化热绝热温升。混凝土水化热绝热温升，在混凝土内部各计算点按公式(附83)计算，在基岩内各计算点=0；在混凝土与基岩接触点取；在新老混凝土接触点取。 初期水管冷却。混凝土重力坝坝内埋设蛇形冷

7、却水管时，在浇筑收仓数小时开始通水进行初期冷却，以削减水泥水化热温升，这时计算中要考虑水管和浇筑块顶面同时散热，可用双向差分法公式进行计算。 在浇筑块中取一长条，其宽度和冷却水管的水平间距相同，在选择时应使水管正处在中间的计算点上，但冷却水管所在的那个计算点，通水过程中始终保持常温通水水温，其值可取1/2(进口水温+出口水温)，而两侧的边界作对称处理见附图15。附图15附图16 3)简捷计算法： 热传导的微分方程和辅助条件均为线性的，可以用迭加的原理先把复杂的条件简化成各种不同单元，然后迭加而成(见附图16)。对于初期有顶面散热和内部通水冷却的浇筑块，在计算浇筑块的平均最高温度Tm值时，可采用

8、以下公式计算。 无水管冷却 (附86) 有水管冷却 (附87)式中 Tm混凝土坝块平均最高温度； Txi下层混凝土的温度； Tb混凝土表面温度； Ts冷却水管初期通水的水温； X冷却水管散热残留比，曲线见附图17； E1底部不绝缘，上层新混凝土接受下层混凝土传热并向顶面散热的残留比，曲线见附图18； E2底部不绝缘，上层新混凝土向下层混凝土及顶面散热的残留比，曲线见附图19； Tr通过顶面及冷却水管散热之后的水化热温升，其值可用时差法求得。附图17中 混凝土的导热系数，kcal/(mh)； D、R冷却圆柱体的直径，m，半径，m； r冷却水管的半径，m； Cs水的比热，kcal/(kg)； qs

9、冷却水的流量，m3/h； l每根连续的冷却水管的长度，m。附图17 埋置水管冷却从始点到l长的全部圆柱体平均温度计算图不同值的 在短间歇均匀上升的情况下，为简化计算可令Txi=Tm则公式(附86)和(附87)可简化为： 无水管冷却 (附88) 有水管冷却 (附89) 在具体计算中往往只需求最高温度，不必求各天的平均温度，由于发生水化热最高温升的龄期即为混凝土发生最高温度的龄期，因此只需要把逐日发生的水化热温升列表算出，从中确定发生最高温升的龄期，根据这个龄期从附图17附图19查得E1、E2、X值代入公式(附86)、(附87)或(附88)、(附89)即可一次求得坝块平均最高温度。举例说明如下：

10、某坝块浇筑层厚度3.0m，间歇期7天，华新矿渣大坝500#水泥214kg/m3，浇筑温度22，浇筑后旬平均气温28，水管间距水平2.0m，垂直3.0m，初期回水水温20，=2.77kca1/(mh)，a=0.126m2/d，水管直径d=2.54cm，Cs=1.0kcal/(kg)，=1000kg/m3，qs=1.06m2/h，l=200m，求坝块平均最高温度。附图18 新混凝土受老混凝土固定热源作用的残留比曲线附图19 新混凝土固定热源向空气和老混凝土传热的残留比曲线 冷却圆柱体的直径按等周长计算见公式(附93)。(为化引的混凝土的导温系数，m2/d) 查附图17得X见附表17。附表17(d)

11、0.51.52.53.54.5X0.00590.990.01770.970.02950.9550.04130.9350.05310.922 华新矿渣大坝500#水泥的水化热试验值和混凝土的绝热温升(单位水泥用量214kg/m3)见附表18。附表18龄期(d)1234567水泥水化热试验值(cal/g)混凝土绝热温升()29.8712.1038.4315.5644.4217.9848.2019.5251.2920.7853.2321.5554.4722.05查附图19得E2(见附表19)。附表19(d)0.51.52.53.54.5F0E20.0070.860.0210.760.0350.690

12、.0490.630.0630.58 通过顶面和冷却水管同时散热后的水化热温升用附表20计算。 从附表20可见最高温升发生在第三天，其混凝土最高温度也同样发生在第三天。则=3天时残留比如下：E1=0.115 E2=0.655 X=0.945E1X=0.1150.945=0.109(1E1X)=0.891E2X=0.6550.945=0.618附表20Tr(E2)0.5X0.5(E2X)0.5d1d2d3d4d5d6d1234512.1015.5617.9819.5220.7812.103.462.421.541.270.860.760.690.630.580.990.970.9550.9350.

13、9220.850.740.660.590.5310.308.952.948.002.562.067.152.281.791.316.422.041.601.141.68Tr 10.3011.8912.6212.5312.28 注：附表20中(E2X)0.5d的说明：在用附表20计算水化热温升的残留比E2和X必须采用计算时段的中值，即第一个时段(以一天为一个时段)采用t=0.5d时的E2和X值，第二个时段采用t=1.5d时的E2和X值，余此类推。在E2和X的右下角注以0.5的注脚，即(E2)0.5、X0.5和(E2X)0.5d以示区别。=0.0550.655=0.036=28+5=33 代入公式

14、(附89)得： 2.坝体混凝土的后期温度计算(无热源的温度场计算) (1)计算目的： 坝体混凝土的天然冷却计算目的，主要在确定坝体长期暴露部位的内外温差；长间歇继续上升时的上下层温差以及坝内埋设冷却水管确定接缝灌浆时间或在指定时间内下降至指定温度所需水管间距和流量。 (2)坝体混凝土的天然冷却计算： 设坝体具有初始的均匀初温T0，两侧暴露在稳定的气温Tq中，其平均温度的变化可以用下式计算： (附90)式中 E温度残留比，为参变数的函数(l为坝块宽度)E值可取自附图20。附图20 平板厚度为l初温为均匀分布的平均温度计算图不同值的残留比E值 实际上坝体冷却时，边界气温Tq不是常数，也有变化，此时

15、可按照数值法进行计算，即将整个冷却过程划分为若干时段，对于每一个时段，假定气温为常值，而应用上述公式(附90)计算，然后用迭加法求出最终成果。 用数值法计算时，可以采用附表21的形式进行计算。 当坝体具有两个散热面时(如宽缝重力坝)，其计算原理和方法同上，只需以代替E即可，式中lx和ly为坝块两个方向的宽度。 上述浇筑块的中心温度Tc可以采用相同的步骤计算，只须把公式(附90)中的E值改为Ec值即可，Ec值可自附图21中查得。附图21 求出坝块的平均温度Tm和中心温度Tc后，内外温差即为T=TcT，T为边界温度，作为近似估计，可取T=Tq(气温)，其成果稍偏安全，如欲估计边界温度和气温的差别，

16、可用下式估算。 内外温差=TcT=(TcTq) 式中为折减系数，是参变数/和温度奕化周期的函数，可取附表22，为混凝土表面对空气的热交换系数。附表21时段气温()温差()时间(d)Ei平 均 温 度 计 算第1时段第2时段第3时段第4时段1Tq1T0-Tq1E1E1(T0-Tq1)E2(T0-Tq1)E3(T0-Tq1)E4(T0-Tq1)2Tq2Tq1- Tq22E2E1(Tq1-Tq2)E2(Tq1-Tq2)E3(Tq1-Tq2)3Tq3Tq2- Tq33E3E1(Tq2-Tq3)E2(Tq2-Tq3)4Tq4Tq3- Tq44E4E1(Tq3-Tq4)MMMMMMMMMMMM合 计E1(

17、T0-Tq1)+Tq1附表22/温度变化瞬时降温或日变化长期气温变化15(d)1(a)0.10(m)0.20(m)0.610.420.870.770.970.94 (3)后期水管冷却计算： 冷却水管间距一般为1.5m1.5m3.0m3.0m，管径一般为2.54cm或1.9cm，通水流量为1420L/min。冷却直径D的确定，当水管呈方形布置时： (附91) 当呈长方形布置时： (附92) 等周长换算公式为： (附93) 当梅花形布置S1=1.1547S2时，按等面积换算公式为： (附94)式中 S1水管水平间距； S2水管垂直间距。 已知开始通水冷却时混凝土初温为T0，蛇形水管进口水温为Ts，

18、则通水冷却时的混凝土平均温度为： (附95) X曲线见附图17。 当水管冷却的同时有侧面自然散热时，可以采用时差法计算，具体计算方法和步骤举例说明如下： 某宽缝重力坝第坝块后期水管冷却，坝块宽l=16m，水管间距2.0m2.0m，12月1日开始通水，混凝土初温24，通水水温4，坝块的稳定温度8，需通水冷却多少天始能进行接缝灌浆? =2.77kcal/(mh)，a=0.126m2/d，水管直径d=2.54cm，Cs=1.0kcal/(kg)，=1000kg/m3，qs=1.08m3/h，水管长度l=200m。 从附图17查得水管散热残留比X填于附表23中的第11栏中。 从附图20中查得自然散热残

19、留比E填于附表23中的第7栏中。 从附表23的计算结果可以看出要达到稳定温度8的通水时间为40天，即元月10日可以进行接缝灌浆。附表23 初温T0=24(1)时段(2)月份(3)平均气温()(4)气温差()(5)时间(d)(6)当量时间(7)残留比E(8)(9)平均温度计算(10)当量时间(11)残留比X(12)12345112月上旬6.217.8100.004920.840.8414.9310.659.806.746.284.193.952.542.400.20160.750.752中旬4.71.5200.009840.7750.9221.261.161.081.020.40320.560.

20、743下旬4.00.7300.014760.7230.9320.590.540.500.60480.4150.7441月上旬3.10.9400.019680.6850.9440.760.700.80640.3080.745中旬3.10500.024600.6460.94401.00800.230.746下旬3.2-0.1600.029520.6100.9441.20960.1680.74(a)旬平均气温6.24.74.03.13.1(b)不考虑水管冷却达到的温度21.1315.7612.039.437.72(c)水温4.04.04.04.04.0(d)与水温之差=(b)(c)17.1311.7

21、68.035.433.72(e)水管冷却之温度=(d)(1-)4.283.062.091.410.97(f)最终平均温度(b)-(e)16.8512.709.948.026.75 附表23中的计算步骤和注意事项如下： 表中第8栏、第12栏中的，为下一时段的残留比和本时段残留比之比，如。 计算从第一时段开始，混凝土初温为24，12月上旬平均气温为6.2，相差17.8，此值乘以平均第一时段之残留比0.84，得14.93，故若不计水管影响，在第一时段末的平均温度为14.93+6.2=21.13，此值与水管水温4相差17.13，再根据水管散热的第一时段残留比0.75，可知道在第一时段中水管能散去温度为

22、17.13(1-0.75)=4.28，故第一时段末之最后平均温度是21.13-4.28=16.85。这些计算都依次列在附表23下部(a)(f)等六行中。 既然第一时段中从水管散去4.28，则原来不考虑水管影响而得出的14.93中应减去4.28，改正为10.65，这改正就记在表上，并以表示。 在第二时段，将上述10.65乘以平板第二时段残留比(注意需用第8栏值，不可用第7栏值，因为现在是从第一时段末的温度去计算第二时段温度)0.922得9.80，此外第一、二时段的气温差6.2-4.7=1.5乘以平板第一时段残留比0.84，得1.26，故不考虑水管影响的第二时段末之平均温度为9.80+1.26+4

23、.7=15.76，比水管水温高11.76，利用水管散热第二时段的残留比0.74(注意需用第12栏而非第11栏)，知第二时段中水管可散去温度11.76(1-0.74)=3.06，因此第二时段最后平均温度为15.76-3.06=12.70，将水管散去的3.06全部放在9.80中改正，即9.80-3.06=6.74，并以表示。 第三时段计算时将6.74和1.26分别乘0.932和0.922得6.28和1.16，将气温差0.7乘以0.84得0.59，于是不计水管影响的平均温度为6.28+1.16+0.59+4.0=12.03，较水温高8.03，因而从水管散去8.03(1-0.74)=2.09，则第三时

24、段末的最终平均温度即为12.03-2.09=9.94。 再将水管散去之温度2.09放在6.28改正得4.19记在表中以表示，以供第四时段计算用，以下时段的计算仿上述的步骤，不再详述。 二、温度应力计算 1.基础浇筑块温度应力计算 基础浇筑块的温度应力主要指基础温差Tj+Tr-Td所引起的浇筑块中央断面上的水平正应力，可按下式计算： (附96)式中 由均布温差Tj-Td所产生的温度应力； 由非均布的水化热温升Tr所产生的温度应力； Td坝体稳定温度。 应力可按下式计算： (附97)式中 kp由于混凝土徐变引起的应力松弛系数，当缺乏专门试验资料时，kp可取为0.50； 混凝土的泊松比； a混凝土的

25、线膨胀系数，1/； R基础约束系数，当混凝土的弹性模量Eh等于基岩的弹性模量E时，R可按附表24采用。附表24y/l00.10.20.30.40.5R0.610.440.270.160.100注：l浇筑长度，m； y计算点离基础面高度。 当混凝土弹性模量Eh与基岩弹性模量E不等时，R值可按附表25所列数值增减。附表25Eh/E00.51.01.52.03.04.0R(y=0)1.00.720.610.510.440.360.32 应力可按下式计算(参阅附图22)： (附98)式中 T(y)应力计算点y处的温度值； y坐标y的增量； l浇筑块的长度，m； 在y=处加一对单位荷载P=1对计算点y所

26、产生的正应力影响系数，可由附图23及附图24查取； T()在y=处的温度； kq考虑早期升温阶段的压应力的折减系数，可采用0.750.85。附图22 基础浇筑块附图23 浇筑块温度应力影响线Eh=E附图24 浇筑块温度应力影响线Eh=E/2 应力也可按下式估算： (附99)式中系数A可从附图25中查得。附图25 水化温升应力系数A 2.上下层温差引起的温度应力计算： 首先，用差分法计算上下层中的温度分布，然后按下式计算任一点y处的温度应力(参阅附图26)： (附100)式中 在点处作用一对单位水平力对y点正应力的影响系数，可自附表26查得。附图26 上下层温差引起的应力计算附表26 上下层应力

27、影响系数By()i00.10.20.30.40.50.60.70.80.91.000.10.20.30.40.50.60.70.80.91.01.11.21.36.684.552.681.330.4230-0.394-0.448-0.421-0.327-0.2324.543.562.501.570.8260.3400-0.20-0.25-0.20-0.122.682.492.231.771.260.7800.4280.190.100-0.101.331.551.761.781.611.230.8570.5410.2610.1000.4280.801.101.421.481.581.240.91

28、0.4310.2200.1000.3220.7951.241.501.681.591.290.8740.5550.290-0.37400.4310.8761.301.591.781.601.230.8740.539-0.44-0.200.190.540.911.291.601.701.601.200.850.44-0.414-0.250.100.2470.5310.8751.231.561.741.531.220.850.44-0.33-0.2000.100.220.560.811.201.601.681.591.200.85-0.234-0.122-0.1000.110.2960.5380.

29、8541.201.591.681.651.220.85 如果温差系由水化热引起，所得应力须再乘以折减系数kq，以考虑早期升温阶段压应力的影响，kq可取0.750.85。 计算时将浇筑块沿高度分为若干层，每层高度为y，一般各层取为等高，公式中的符号表示各层的总和。 3.外界温度年变化引起的应力： 外界温度(水温、气温)的年变化可用正弦函数表示如下： (附101)式中 T0外界温度变幅，； (为温度变化周期，=365d)。 在深度为y处的温度为： (附102) 式中 (附103) 由(附102)式算出温度后，再由上下层温度应力计算公式(附100)可计算任一点在任意时间的应力，对于年变化温度应力，可

30、取应力松弛系数kp=0.65。 最大应力一般出现在坝块表面(y=0)，参阅附图27其值可按下式估算： (附104) 式中系数C可自附图27查得，应力松弛系数kp可采用0.65，附图27是按照导温系数a=0.10m2/d计算的，如果导温系数a相差较远，可按下式计算等效坝块长度：附图27 年变化温度应力系数C 4.气温骤降引起的混凝土表面应力： 如附图28所示，设日平均气温在期间急剧降低T0度，然后回升，气温骤降在混凝土表面产生的最大温度徐变应力可按下式估算： (附105)式中T0为气温降低幅度，D为系数可自附图28查得。图中为降温历时(天数)，为混凝土导热系数，一般可取=2.4kcal/(m2h

31、)，为表面放热系数，对于裸露表面，=1020kcal/(m2h)。附图28 气温骤降产生的表面应力系数D 当混凝土表面覆盖了保温材料时，应采用等效放热系数以代替，等效放热系数可按下式计算： (附106)式中 h1保温材料厚度，m； 1保温材料导热系数，kcal/(mh)。 例如，表面有2.5cm厚木板时3.7，表面覆盖5cm厚度草袋时，2.5，表面覆盖1cm厚度泡沫塑料时，22.3，单位均为kcal/(m2h)。 5.弹性基础梁温度应力： 浇筑块早期温度应力可按弹性基础梁计算。如附图29所示，设梁的长度为2l，高度为2h，宽度为1，混凝土弹性模量为Eh，梁内温度为T(y)，基础弹性模量为E，混

32、凝土和基础的泊松比假定均等于1/6。附图29 弹性基础梁 首先沿梁和基础的接触面切开，使梁脱离基础约束，按照平面应变和平截面假设，可得梁内部互相制约产生的自生应力如下： (附107) (附108) (附109)式中Tm为梁的平均温度，为平均温度梯度，它们可用数值积分公式(如梯形公式或辛浦生公式)由公式(附108)和(附109)求出。 基础约束作用在梁内x断面上产生的弯矩M(x)和轴向力P(x)分别如下： (附110) (附111) 梁中央断面上的弯矩M0和轴向力P0为 (附112) 求出弯矩M和轴向力P后，由下式可求出由于基础约束作用在梁内产生的温度应力如下： (附113) 约束应力在上缘和下

33、缘之间呈线性变化，自生应力一般说来是非线性变化的。约束应力和自生应力迭加后就可得到梁内温度应力。以上计算中未考虑混凝土徐变，因此按以上各式计算温度应力后还应再乘以松弛系数kp。 公式(附110)和(附111)中的系数B1和A2可由下列二式计算： (附114) (附115) (附116)式中 r=h/l；=Eh/E。 在计算早期温度应力时，为了考虑混凝土弹性模量随着龄期而变化，如附图30所示，将时间划分为许多时段，在每一时段内取一平均弹性模量。 (附117) 在每一时段内的温度增量为该时段首末各点的温度差，即 (附118) 根据平均弹性模量Ehi和温度增量Ti(y)，用以上各式计算第i时段内产生的应力增量，将各时段的应力增量累积起来，即得到梁内温度应力。 (附119) 6.温度应力控制： 为了防止裂缝，坝体温度应力必须满足下列条件： (附120)式中 由于各种温差所产生的温度应力，kg/cm2； 混凝土的极限拉伸值，对于重要的工程应通过试验测定，无试验数据时可