1 概 述

某大坝采用混凝土重力坝型，坝顶高程257.50 m，坝轴线处坝顶全长205.00 m，最大坝高46.00 m，自左至右依次布置为左岸非溢流坝段、溢流坝段、排漂孔坝段、冲砂泄洪底孔坝段和右岸非溢流坝段。溢流坝闸室采用开敞式结构，共3孔，布置于河床中部，堰顶高程241.00 m，单孔净宽12.00 m，堰体采用WES曲线实用堰型。闸墩顺水流方向长48.00 m，闸墩顶部设有交通桥、门机轨道桥、人行桥、液压管道桥各一座。本分析针对典型坝段，在给定施工进度条件下进行施工期的温度场和徐变应力场仿真计算，推荐合理的施工混凝土温控方案。

2 温度场及温度应力计算原理

2.1 温度场计算原理

为全面反映温度对坝体结构特性的作用与影响，需要研究坝体施工期的温度场、初期蓄水过程中坝体随气温与水温等因素变化的变化温场、运行蓄水期的稳定(准稳定)温度场。根据热量平衡原理，可导出固体热传导基本方程[1]：

 

(1)

式中：为导温系数，λx、λy、λz为导热系数，c为材料比热，ρ为材料容重，τ，T代表叙述任意时间及温度。将求解区域R划为有限个单元Ωe，引入单元形函Ni，则单元内任意点的温度可由构成单元m个节点温度插值：

 

(2)

基于变分理论推导出有限元支配方程：

[H]{T}+{F}=0

(3)

式中：

日粮蛋白质水平对“京红1号”蛋种鸡育成期（9～13周龄）生长性能的影响见表2。由表2可知，日粮蛋白质水平对13周龄体重有极显著的影响（P＜0.01），随日粮蛋白质水平增加，体重呈线性增加趋势（P＜0.01）。P4组（CP 15.7%）13周龄体重最大且与其他各组差异极显著（P＜0.01），日粮蛋白质水平对蛋鸡平均日增重的影响也呈现相同趋势。

越南肉桂是中国肉桂的变种，味型、香气与中国肉桂类似。由于越南与防城地区接壤，相较而言，越南桂精油无论香气感官评价还是成分分析，都与防城桂精油更接近。

(4)

 

(5)

2.2 温度应力计算原理

溢流坝顺河向应力包络线、横河向应力包络线及其坝体沿高程分布图8～10。从图中可以看出在溢流坝体底部垫层出现较大的应力，主要是由于浇注的混凝土受到地基的强烈的约束所致，最大拉应力达到1.09 MPa，小于允许应力。在溢流坝的挑流圆弧部位出现了较大应力，应加强布设钢筋，以防开裂。

而：

水化热以体积力的形式施加在混凝土单元上，实际计算时取前后两个时间步的水化热之差，并视坝体混凝土材料分区不同而作相应的拟合处理。

(6)

式中：[K]为整体刚度矩阵；{ΔPn}L、{ΔPn}C、{ΔPn}T、{ΔPn}0、{ΔPn}S分别为外荷载、徐变、温度、自身体积变形和干缩引起的结点荷载增量。

卡帕系数(Kappa)是一种计算分类精度的方法[19].kappa计算结果为-1～1，但通常kappa是落在 0～1 间，可分为5组来表示不同级别的一致性：0.0～0.20极低的一致性(slight)、0.21～0.40一般的一致性(fair)、0.41～0.60 中等的一致性(moderate)、0.61～0.80 高度的一致性(substantial)和0.81～1几乎完全一致(almost perfect).

求出{Δδn}后，结构应力易于求得。

溢流坝体型随混凝土分层浇筑的同时亦在不断变化，在仿真分析时采用生死单元进行定义，根据整体建立的模型，依据溢流坝体混凝土实际进度而依次激活单元，相应的各种荷载同步施加，这样就可以仿真计算溢流坝体从施工期到运行期全过程的温度场和应力场。为使死单元不影响计算收敛及激活单元的计算结果，在温度场分析中，将比热矩阵和热传导矩阵乘以一个因子，且单元热通量设置为0。

 

(7)

逐时段累加后，得到各个单元τn时刻的应力为：

土壤环境及质量直接关系到农产品安全、生态安全和人体健康。随着信息技术的发展，信息化监测技术在土壤监测上发展极为迅速，具有快速、便捷、实时的特点。检测出土壤污染物类型及其含量，了解土壤基本情况及受污染状况[5]。遥感技术、地理信息系统、全球定位系统可快速调查土壤状况、采样布点的功能，可动态监测土壤侵蚀、沙化、次生盐碱化等问题。传感器网络技术可监测土壤温度、湿度、盐分等，该技术通过将监测的数据传到数据处理中心，监测人员可随时读取数据。

{σn}=∑{Δσn}

(8)

3 设计温控措施和温控标准

3.1 混凝土主要浇筑温控措施

3.1.2 坝体混凝土初步分析采用施工进度计划

对照品：磷酸肌酸钠（批号100875‐201302）、肌酸（批号100876‐200901）、三磷酸腺苷二钠（批号140674‐200401）、二磷酸腺苷二钠（批号 140809‐201501）、一磷酸腺苷钠（批号 14719‐200501）、肌酐（批号100877‐200901）均由中国食品药品检定研究院提供。

从加快混凝土施工进度，简化温控措施角度出发，考虑坝体埋水管冷却措施，优先考虑预冷骨料、加冰拌和等措施控制混凝土浇筑温度，实现控制坝体温度应力的目的。混凝土采用自然入仓，不专门采用保温和通水措施。

3.1.1 混凝土浇筑上升方式及浇筑温度

坝体混凝土浇筑施工按施工进度计划进行，每个坝段采取通仓浇筑，基础约束区(0.2～0.4 L、L为浇筑块长边尺寸)采用浇筑层厚1.50 m，间歇5 d；脱离约束区采用浇筑层厚2.80 m，间歇6 d。根据溢流坝段(排漂孔坝段、冲沙泄洪底孔坝段)的断面结构，223.50 m高程以下为基础约束区，约束区高10.50 m，分7层浇筑，每层厚1.50 m，平均每8 d浇筑一层，需56 d，223.50～237.50 m分5层浇筑，每层厚2.80 m，平均每9 d浇筑一层，需45 d，浇筑至237.50 m高程需101 d。

3.2 混凝土温控标准

为了研究初步可行的温控措施，需要依据具体材料热力学特性分析温度及温度应力，从而确定初步的温差控制标准。经分析，在安全系数取为1.8时，对溢流坝段通仓浇筑时基础允许温差强约束区0～0.2 L内部混凝土16℃，0.2～0.4 L为19℃。表1为根据各标号混凝土不同龄期的弹性模量以及相应的极限拉伸值，计算出的溢流坝混凝土的允许拉应力。

 

表1 溢流坝混凝土施工期容许温度应力

  

项目极限拉伸值/(εp×10-6)抗拉弹模/(104MPa)允许应力/MPa7d28d7d28d7d28dC9015W90665.277.11.822.460.721.21C9020W90670.582.32.322.611.001.41

4 仿真分析计算处理

4.1 溢流坝段混凝土浇筑模拟

徐变度计算：

第二类，商业聚集型。大多是河边的码头。如酉水流域的王村、里耶、洗车河、石堤、大溪、百福司、沙道沟；乌江流域的龚滩、淇滩、思渠等；清江流域的兹丘、景阳河、庆阳坝。

4.2 温度初始条件

溢流坝体新浇筑混凝土为浇筑温度，新老混凝土交界处采用新老混凝土层的温度平均值。

4.3 边界条件

地基模型上下游侧及底部为绝热边界条件处理。溢流坝上游面跟库水相接触边界为第一类边界条件，取河水温度(时间函数)，空气接触边界采用第三类边界。

4.4 水化热

[K]{Δδn}={ΔPn}L+{ΔPn}C+{ΔPn}T+{ΔPn}0+{ΔPn}S

王振辉指出，京东物流提出GSSC的时候也提出了一个目标，叫双48小时。第一个“48”是指从中国到全球或者全球到中国实现48小时内通达。第二个“48”是各国本地的交付，在中国次日达已经成为物流的标配产品，当日达即“211”（当天11点前下单，当日达）在七八年前就已经达成了。通过全网的覆盖、全通路的覆盖，让京东物流和合作伙伴能获得全球客户和消费者充分的信任。从Copyto China到Leadingthe Change，中国的物流模式已经开始在世界领先，京东物流有能力把硬件技术和软件技术覆盖到全世界。

4.5 混凝土热力学性能

在仿真分析计算应力时，导温、导热及线膨胀系数采用试验值，抗拉弹模根据实验数据拟合成指数(时间函数)，C9015及C9020徐变度参照类似工程经验，采用朱伯芳院士提出相关函数拟合。

4.6 初始地应力

考虑到初始地应力状态下，基岩基本上处于弹性状态，为了简化，可直接采用现有河谷形状的基岩在自重作用下产生的应力作为初始地应力。

5 溢流坝段施工期温度及应力场仿真分析

溢流坝体稳定温度场及应力场三维计算网格如图1、图2所示，其中建基面高程以下基岩厚度约1.5倍坝高，坝轴线上下游侧顺河向范围约1.5倍坝高[3]。共离散为9 290个节点和7 048个单元，采用空间8节点等参实体单元。

由施工期最高温度包络线图3可知，在溢流坝上部出现一个高温区，主要是该区域浇注温度相对较高，溢流坝的最高温度达到47.4℃。分析选取了典型部位的节点进行温度和应力分析，节点选取见图4。

  

图1 溢流坝段整体有限元计算网格

  

图2 溢流坝段坝体有限元网格

  

图3 施工期最高温度包络线图

  

图4 节点选取示意图

典型高程温度历时曲线见图5～7中的典型部位温度历时曲线。从图中可以看出，溢流坝体的表面在浇注后温度略有上升，迅速与气温达到一致。溢流坝体中心部位在混凝土浇注后，温度上升存在很长一段时间。溢流坝中心混凝土，在基础约束部位，在浇注50 d后达到最高温度，溢流坝体上部在20 d左右后达到最高温度，然后由于水化热速度降低，随着坝体表面的散热，温度逐渐下降。从图中可以看出，溢流坝体中心部位的温度受到外界环境的温度影响较小。

由于混凝土浇筑后弹性模量和徐变度都随时间而变化，可以采用增量法计算。对于此法的变应力作用下混凝土徐变应力的计算[2]，给出了隐式解法，与初应变法相比计算精度大大提高。根据虚功原理，不难导出计算混凝土结构在热、力作用下单元在时段n时的温度徐变应力矩阵方程：

溢流坝体中心部位顺河向应力、上游表面横河向应力历时曲线见图11～12，可知溢流坝新老混凝土因互相约束，局部层面产生较大拉应力。坝面新浇混凝土在最初迅速达到应力最大值，在水泥水化热及外界环境温度交换影响下，产生拉应力区，但小于允许应力，后随环境温度逐渐升高，受到膨胀约束，局部产生压应力区，在夏季高温时尤其明显。

  

图5 典型部位温度历时曲线(1)

  

图6 典型部位温度历时曲线(2)

  

图7 典型部位温度历时曲线(3)

  

图8 溢流坝段应力图

  

图9 施工期顺河向应力包络线图

  

图10 施工期横河向主应力包络线图

  

图11 典型部位横河向应力历时曲线

  

图12 典型部位顺河向应力历时曲线

距溢流坝面相近部位，因绝热温升较大且散热条件差，达到温度峰值，后随坝面热交换，温度逐渐降低，但同时溢流坝中心部位混凝土因散热慢，温度仍较高，因此极易出现较大温度梯度变化，在温差及周围混凝土约束条件下，距溢流坝面相近部位混凝土产生较大的拉应力区，受环境温度及溢流坝面混凝土热胀冷缩双重影响下，该拉应力一直维持较高值，导致在坝体中下部的距溢流坝面相近部位和坝体上部的中心部位出现了较大的应力。从上下相邻部位的温度历时曲线可以看出，上下两个浇注层的温差不是产生拉应力的主要原因。坝体的各区域混凝土的应力结果见表2。

 

表2 溢流坝段施工期应力结果 MPa

  

混凝土坝体区域C15C200～0.2L0.2～0.4L0～0.2L0.2～0.4L顺河向应力0.881.001.320.79横河向应力0.901.121.291.16允许应力1.201.201.401.40

注：允许应力取为龄期28 d的允许应力。

6 结 语

从仿真计算结果可以看出，溢流坝各个部位应力大部分小于允许应力，局部略大于允许应力。夏季在距离表面5 m左右的部位会产生应力，但应力小于允许应力，冬季表面产生的应力相对较小。坝体基础垫层由于受到地基的约束，产生了较大拉应力，但小于允许应力。坝体由于相对较小，基于仿真分析结果可不需专门采取通水和保温等温控设计措施，亦能满足工程要求。

参考文献：

[1] 朱伯芳.大体积混凝土温度应力与温度控制[M].北京：中国电力出版社，1999.

当其它变量不变时，物流服务水平敏感系数θ增大，产品销售量增大，产品零售价格提高，由此带来总利润的增大；同时，产品销量的增大导致物流服务规模的增大，规模经济效应显著，从而使得物流服务水平提高；另外，物流服务水平敏感系数θ增大，代表客户对物流服务水平的敏感性提高，即物流服务水平的提高会实现客户效用的最大化，因此为实现客户效用最大化的目标，物流服务水平会随之提高。

[2] 吴桐舟.混凝土徐变对坝体温度应力场的影响[J].水利与建筑工程学报，2016(4)：147～151.

[3] 吕文丽.碾压混凝土重力坝整体三维仿真分析[J].水利规划与设计，2011(3)：40～43.

我一帆风顺，一点也没想到会有厄运降临。同宿舍好友反映：李汝庆早晨唱东方红，不是唱毛主席是人民大救星，而是人民大舅子。组织上立即找我核实，其实真没这事，我相信党组织会公正处理，所以把自己的一切都向党坦白。半年过去了，没想到处分决定是遣送回家！更没想到这只是厄运的开始，“海外关系”与香港“特嫌”成了一直笼罩在我头上的阴影。