高坝结构分析有限单元法常晓林重力坝.pptVIP

有限单元法在水工结构分析中的应用——高坝结构分析的有限元法 常晓林 武汉大学水利水电学院 主要讲授内容 1 重力坝应力与稳定分析 2 拱坝的应力与稳定分析 3 土石坝的应力与稳定分析 计算目的 校核坝体的强度 校核坝基稳定 为孔口配筋提供依据 论证坝基加固措施的可行性 高坝结构有限元分析过程（重力坝） 合理地概化几何模型 有限元网格的剖分 边界条件和荷载的施加 材料本构模型的选择 施工工程的模拟 计算成果的整理、评价 几何模型的概化 坝体的孔口闸墩 坝基——岩块、断层、节理、夹层 计算范围的选取 平面计算/三维计算 有限元网格的剖分 原则——关心的部位网格应该密，疏密均匀过渡 单元选择——尽量采用四节点的接近于正方形单元（平面），八节点的立方体单元（三维），以保证单元插值精度 在同一个连同域内单元应该共节点、共边、共面 边界条件和荷载的施加 计算域的边界－施加法向位移约束 荷载：自重、水荷载（面力或者体力）、泥沙压力、温度荷载 材料本构模型的选择 拱坝、重力坝坝体——弹性、弹塑性（混凝土材料模型） 混凝土 单轴本构 模型 混凝土 双轴模型 混凝土三维破坏面 混凝土材料模型 破坏准则 莫尔库仑准则 Drucker－Prager准则（带拉断） Willam－Warnke三参数、五参数准则 混凝土材料本构模型 坝基材料模型 Drunker－Prager准则 莫尔库仑准则 莫尔库仑准则 为避免在屈服面出现尖点，采用修圆措施 断层、节理、夹层材料模拟 采用各向异性材料：沿断层、夹层和节理的顺层方向：采用M-C、D-P（包括软化特性），法向不抗拉/低拉 单元类型选择 坝体——实体高精度等参单元（平面4边形、三维六边形） 坝基岩体——实体高精度等参单元 坝基断层、夹层——有厚度节理单元 坝基节理——无厚度节理单元 施工工程的模拟 坝基地应力的模拟 构造应力强烈：采用实测资料进行地应力反演 无构造应力：近似施加自重荷载模拟 坝体施工过程 对于重力坝，可以采取一次激活坝体的方式（自重全部由坝体本身承担） 坝体蓄水过程的模拟 分若干荷载步施加上游表面水压力和淤沙应力 计算成果的整理、评价 铅直向位移分布（单位：米） 计算成果的整理、评价 顺河向位移分布（单位：米） 计算成果的整理、评价 第一主应力分布（单位：Pa） 计算成果的整理、评价 第三主应力分布（单位：Pa） 计算成果的整理、评价 主应力矢量 计算成果的整理、评价 主应力矢量 计算成果的整理、评价 建基面第一主应力沿顺河向的变化曲线 计算成果的整理、评价 中孔进水口段横河向、铅直正应力分布图 计算成果的整理、评价 铅直向位移分布（单位：米） 重力坝有限元应力计算控制标准 《混凝土重力坝设计规范》（DL5108-1999）规定用有限元计算时，“当计扬压力时，拉应力区宽度应小于坝底宽度的0.07倍（垂直拉应力分布宽度／坝底面宽度），或不超过帷幕中心线” 重力坝渐进破坏过程模拟方法 超载法 超载法主要考虑作用荷载的不确定性，以此研究结构承受超载作用的能力 该方法较直观，便于在结构物理模型试验中采用，从而使数值模拟与物理模拟结果相互印证，且积累了较多的工程经验 但要使结构达到最终整体失稳的极限状态，其相应的超载系数是很大的，而实际上结构的这种荷载状态几乎是不可能出现的，故这种方法求得的超载系数只是结构安全度的一个表征指标 重力坝渐进破坏过程模拟方法 超载法 超载法认为，作用在坝上的外荷载由于某些特殊原因有可能超过设计荷载，超过的总荷载与设计总荷载之比称为超载系数 用逐渐增加超载系数研究坝从局部到整体破坏的渐进破坏过程的方法，称为超载法。超载法分为超水容重Kγ(三角形超载)和超水位KH(矩形超载)两种方法 重力坝渐进破坏过程模拟方法 强度储备系数法 主要考虑材料强度的不确定性和可能的弱化效应，以此研究结构在设计上的强度储备程度 天然岩体由于成因和结构构造运动其不均匀性非常明显，节理、裂隙和断层发育且分布规律复杂，很难准确地把握其工程尺度范围内的物理力学性能，各局部材料参数相差数倍是完全可能的，因此强度储备系数法从这种意义上能较真实地反映结构破坏的实质和可能的失稳模式 但这种方法目前工程应用经验积累还显不足 重力坝渐进破坏过程模拟方法 强度储备系数法 用降低强度参数的方法，研究大坝失稳的渐进破坏过程，如令K表示强度储备系数，K为大于1.0的值，f、c为实际的抗剪强度参数，降强度就是用f／K、c／K代替f、c值进行计算，随着K值的逐渐增大，可以求出大坝从局部破坏到全部破坏的破坏全过程，所算得的整体破坏时的K值的大小也能反映大坝安全的程度 强度储备系数法包括对f、c值采用等比例降强度和不等比例降强度（等