渠道衬砌结构冻胀破坏力学模型及冻胀数值模拟备课讲稿.ppt

在前面渠道冻胀数值模拟的基础上，以某U 形渠观测段为研究对象，利用有限元方法分析了一个冻融周期内渠道的温度场、应力场变化，以期寻求动态温度场引起渠道冻胀破坏的控制状态，分析和预测在冻融过程中基土温度场、应力场及其位移特别是衬砌板上冻胀力和冻胀变形的变化规律，为渠道抗冻胀设计提供更加科学的理论依据。 2.2.2 周期性温度变化下的混凝土衬砌渠道冻胀的数值模拟 图2-23 有限元模型 图2-22 衬砌大U型渠道断面 (cm) 边界条件 （1）瞬态边界条件 上边界温度表示为： 下边界深度足够达到稳定温度，取值为10℃，左右边界视为绝热条件，位移边界条件为渠基左右边界水平方向位移为0，下边界竖直方向位移为0。 （2）稳态边界条件 温度场计算近似取原型各部位月平均表面温度最大值作为有限元模型各部位边界温度，边界约束条件同瞬态分析。 2-24 监测点冻胀量变化曲线 2-25法向应力变化曲线 2-26切向应力变化曲线 注：时间为1月6日至4月6日，共90天 (1) 对渠道冻融过程渠基温度变化进行仿真分析，可知阴坡冻胀期明显长于阳坡，得到的渠道各边坡的冻融日与实测情况一致，阳坡滞后阴坡约15天冻结，正是这种温度场的不对称、不均匀、不同步现象引起衬砌不均匀的冻胀与融沉病害。 (2) 运用有限元方法对渠道冻融过程应力、变形场进行模拟可知：冻胀量、冻胀力与冬季日温度的变化规律比较一致，但最大冻胀量及最大冻胀力滞后日平均最低气温约4天；衬砌底面法向、切向冻胀力发展虽然不同，但是与冻胀量一样都显现出滞后现象，其时间大小分布与相关文献及实测资料一致，说明瞬态数值模拟能掌握渠道抗冻胀设计的控制状态。 因此，运用瞬态结果进行渠道抗冻胀设计更能掌握冻胀破坏的极限状态才是更加合理的。 主要结论： 2.2.3 考虑日照热辐射及热传导的渠道冻胀数值模拟 混凝土渠系结构通常是暴露在外面的，表面与空气或水接触，冬季昼夜太阳辐射对渠道衬砌冻胀影响较大。 本节考虑日照对衬砌渠道冻胀量及冻胀力的影响，并按瞬态温变模式加载温度,对衬砌渠道冻胀过程进行动态数值模拟。 本方法克服了渠道冻胀数值模拟分析依赖于实测温度场的缺陷。 考虑日辐射的温度变化 综合分析热传导方程、定解条件、太阳辐射热的规律等，日照引起的温度升值可以表示如下： 日变化温度为： 计算实例 某U形渠观测段断面形式（同图2-22） 正弦环境温度变化： 辐射的温度增值 渠底： 阴坡： 阳坡： 注：其中 int()表示取整操作。 U形渠道有限元模型（同图2-23） b. 不考虑热辐射的正弦温度变化 a.考虑热辐射的正弦温度变化 图2-27 温度、冻胀量和法向冻胀力比较 主要结论： 考虑到太阳辐射之后，基土和衬砌板的温度场、冻胀力、冻胀量与不考虑差别显著，阳坡差异最大、渠底其次、阴坡最小； 不考虑热辐射的模拟值比考虑之后结果都偏大； 考虑热辐射的渠道冻胀数值模拟方法，克服了通过实测混凝土板各点昼夜温度来确定边界条件的缺陷，使衬砌渠道冻胀预测及动态模拟分析更加方便科学。 2.2.4 渠道冻胀敏感性数值模拟分析 工程实践中已陆续提出了多种行之有效的防冻胀技术措施。如何对防冻胀技术措施进行优选，精确评价各种防冻胀措施的应用效果，探寻效果最显著的防冻胀措施， 本节利用正交实验设计理论，基于渠道冻胀有限元数值模拟，选取渠道冻胀不均匀系数作为评价指标，以渠道断面形式、混凝土衬砌板厚度、分缝位置、渠基土冻胀系数和混凝土衬砌板弹性模量E作为影响因素，对渠道冻胀进行敏感性分析。 　　因素 水平 渠 道 断面形式 衬砌板 厚度（cm） 分缝位置 冻胀系数 砼弹性模量 E（E10 Pa） 1 梯形渠道 8 Ⅰ ×0.5 2.4 2 弧底梯形 10 Ⅱ ×1.0 2.6 3 U形渠道 12 Ⅲ ×1.5 2.8 表2-1 各因素取值 图2-29 分缝位置（以梯形渠道为例） 图2-28 断面形式（cm） 表2-2 极差分析表 表2-3 方差分析表 图2-28影响因素趋势图 (1) 通过极差分析，得出影响渠道冻胀的五个因素的敏感性由大到小依次排列为：渠道断面形式＞基土冻胀系数＞分缝位置＞混凝土衬砌板厚度＞混凝土衬砌板弹性模量E，基本符合冻胀过程中的实际情况。 (2) 由趋势图可得：冻胀不均匀程度随渠基土冻胀系数增大而增大；随混凝土衬砌板厚度增大而减小，即衬砌板越厚冻胀越均匀；随着弹性模量增大而减小，这说明混凝土弹模越高，对冻胀基土约束越强，则冻胀变位越小；在1/3处设缝时，冻胀变形分布最均匀；采用U形渠道能有效地降低渠道冻胀变形。 (3) 根据各因