水工结构计算力学-第三章-周伟-温度渗流原理.pptVIP

水工结构计算力学—热传导和渗流问题的有限单元法 周 伟 　 副教授、博士 武汉大学水利水电学院 主要讲授内容 变分原理 力学问题的有限元法一般原理和表达格式 热传导问题的有限单元法 渗流问题的有限单元法 数值分析方法的分类 1 有限差分法－用差分方程代替微分方程 2 建立等效（弱）积分提法，求近似解——加权余量法（配点法、最小二乘法、伽辽金法）——可以用来构造泛函。解的收敛性没有严格的理论证明，近似解不具备明确的上下界性质 数值分析方法的分类 3 里兹法 要求解的问题的微分方程具有特定的性质（线性、自伴随、偶数阶微分），则等效积分提法可归结为某个泛函的变分。相应的近似解就是求泛函的驻值 缺点：在这个求解区域上假设近似函数，几何形状复杂的求解域，找不到合乎要求的近似函数 数值分析方法的分类 4 有限单元法 有限元法是以变分原理为基础，将要求解的微分方程型数学模型，首先转化为相应的变分问题，即泛函求极值问题；然后，利用分片插值将变分问题离散化为普通多元函数的极值问题，最终归结为一组多元的代数方程组，求解该方程组，即可获得问题的数值解 等效（弱）积分提法 加权余量法 权函数选择的不同：配点法、最小二乘法、伽辽金法（利用试探函数作为权函数） 变分原理 变分原理 原问题的微分方程和边界条件等效于泛函的变分等于零，即泛函取驻值 如果泛函取驻值，等效于满足原问题的微分方程和边界条件 变分原理中泛函的构造 采用伽辽金法提出原问题微分方程和边界条件的等效积分方程 对积分方程实行分部积分，得到等效弱积分方程，降低微分算子的阶数，以降低试探函数的连续性要求 基于变分原理的有限元提法 有限元法是以变分原理为基础 找出问题的微分方程（数学模型）——构造泛函——泛函变分——单元平衡方程——组装整体平衡方程——求解方程组——获得问题的数值解 力学问题变分原理 力学问题变分原理 热传导问题变分原理 基本微分方程 单位时间内沿着x方向流入的净热量 热量平衡 单位时间内固体温度升高所吸收的热量 单位时间内单位体积放出的热量 绝热状态 热传导问题 热传导方程建立了温度与时间、空间的关系，但满足热传导方程的解有无限个。 为了确定我们所需要的温度场，还必须知道初始条件和边界条件。 热传导问题初始条件 在初始瞬时固体内部的温度分布规律 当t=0时 热传导问题边界条件 固体表面与周围介质相互作用的规律 第一类边界条件：固体表面温度是时间t的已知函数 热传导问题边界条件 第二类边界条件：固体表面温度与流体（如空气）接触时，通过固体表面的热流密度，与固体表面温度之差成正比 热传导问题泛函的构造 用分部积分证明上述构造的泛函是正确的？ 泛函变分求热传导问题的平衡方程 渗流问题的基本方程 渗流问题的边界条件 第一类边界条件——水头已知 第二类边界条件——法向流速已知 渗流问题的边界条件 自由面边界 溢出面的边界条件 渗流问题泛函的构造 第三章结束 ！！ 水工结构计算力学—重力坝温控分析 周 伟 武汉大学水利水电学院 主要讲授内容 重力坝温控计算的目的 温控标准和温控措施 混凝土热学力学参数输入 重力坝温控分析的边界和初始条件 重力坝有限元温控分析的网格剖分 重力坝有限元温控分析注意事项 重力坝有限元温控成果分析整理 重力坝温控计算的目的 提出现有混凝土材料参数条件下的温控标准、温控措施、混凝土浇筑方案 温控标准 温度方面 基础温差、内外温差、上下层温差，最高允许温度 温度应力方面 最大允许拉应力 最高温度 最大温度应力 混凝土热学力学参数输入 水化热、绝热温升 导热系数 、导温系数、比热、表面放热系数 弹模、泊松比、徐变、自身体积变形、线胀系数 抗拉强度（劈拉、轴拉）、抗压强度、极限拉伸值 水库水温 稳定温度场 浇筑方案的初步拟定 开浇日期 浇筑温度 浇筑层厚 间歇时间 冷却水管的布置 通水水温、通水时间 表面保护措施 重力坝温控分析的边界和初始条件 重力坝有限元温控分析的网格剖分 计算范围 重力坝有限元温控分析的网格剖分 重力坝有限元温控分析注意事项 考虑坝体混凝土浇筑过程 温度初始条件 边界条件 ：水温、气温为时间的函数 水化热 、绝热温升：时间的函数 水管冷却：等效热传导方程 重力坝有限元温控分析注意事项 材料力学性质 温度及温度应力仿真计算时，混凝土的水化热、弹模和自身体积变形是随时间变化的函数，可根据实际工程的实验资料拟合为相应的曲线形式。徐变度是持荷时间和加载龄期的函数 重力坝有限元温控成果