第5章ANSYS水利工程应用实例分析

本章首先概述了水利工程ANSYS应用，其次介绍了ANSYS重力坝抗震性能分析步骤，最后用实例详细介绍了ANSYS抗震性能分析过程。内容提要第5章ANSYS水利工程应用实例分析本章重点水利工程ANSYS重力坝抗震性能分析步骤ANSYS重力坝抗震性能用实例分析本章典型效果图第5章ANSYS水利工程应用实例分析5.1水利工程概述虽然我国水利资源非常丰富，但河流在地区和时间分配上很不均衡，许多地区在枯水季节容易出现干早，而在洪水季节又往往由于水量过多而形成洪涝灾害。为了解决这一矛盾，人们修建了许多水利工程来达到防洪、灌溉、发电、供水、航运等目的，促进国民经济建设的发展。水利工程中各种建筑物按其在水利枢纽中所起的作用，可以分为以下几类:(1)挡水建筑物用以拦截河流，形成水库，如各种坝和水闸以及抵御洪水所用的堤防等。(2)泄水建筑物用以宣泄水库〔或渠道)在洪水期间或其它情况下的多余水量，以保证坝(或渠道)的安全，如各种溢流坝、溢流道、泄洪隧道和泄洪涵管等。(3)输水建筑物为灌溉、发电或供水，从水库(或河道)向库外(或下游)输水用的建筑物，如引水隧道、引水涵管、渠道和渡槽等。(4)取水建筑物是输水建筑物的首部建筑，如为灌溉、发电、供水而建的进水闸、扬水站等。(5)整治建筑物用以调整水流与河床、河岸的相互作用以及防护水库、湖泊中的波浪和水流对岸坡的冲刷，如丁坝、顺坝、导流堤、护底和护岸等。由于破坏后果的灾难性，大型水利工程建设的首要目标是安全可靠，其次才是经济合理。所以说研究大坝等水工建筑物的安全分析、评价和监控，是工程技术人员需要解决的课题，正确分析大坝性态已经成为当务之急。当前对各种水利工程评价主要采用有限元分析方法，借助各种有限元软件对这些水利工程建筑物进行安全评价，其中应用比较广泛的是ANSYS软件。目前，ANSYS软件在水利工程中主要应用以下几个方面：(1)应用各种坝体工程的设计和施工利用ANSYS软件，模拟各种坝体施工过程以及坝体在使用阶段受到各种载荷（如水位变化对坝体的压力、地震荷载等）下结构的安全性能进行评价，模拟坝体的温度场和应力场，借助模拟结果修改设计或对坝体采取加固措施。(2)应用于各种引水隧道、引水涵管等设计和施工利用ANSYS软件，模拟这些工程开挖、支护、浇注、回填过程，分析结构在载荷作用下的变形情况、结构的安全可靠度，以及衬砌支护结构在水压、温度发生变化后产生的变形情况和结构内力，依靠ANSYS模拟结果对结构安全性进行评价。(3)应用于各种水库闸门的设计和施工水库闸门在上游水作用下将发生弯曲、扭转、剪切和拉压等组合变形，利用ANSYS中的SHELL63单元来模拟闸门，利用大型结构有限元分析程序ANSYS，对闸门结构进行三维有限元分析，根据分析结果进行强度校核。ANSYS10.0土木工程应用实例分析5.2ANSYS重力坝抗震性能分析步骤重力坝是一种古老而重要的坝型，主要依靠坝体自身重力来维持坝身的稳定。岩基上重力坝的基本剖面呈三角形，上游面通常是垂直的或者稍倾向下游的三角形断面。重力坝的具有很多优点：(1)安全可靠。但剖面尺寸较大，抵抗水的渗漏，洪水漫顶，地震或战争破坏的能力都比较强，因而失事率较低。(2)对地形、地质条件适应性强，坝体作用于地基面上的压应力不高，所以对地质条件的要求也较低，低坝甚至可修建在土基上。(3)枢纽泄洪容易解决，便于枢纽布置。(4)施工方便，便于机械化施工。(5)结构作用明确，应力计算和稳定计算比较简单。鉴于重力坝有如此多优点，所以它得到了广泛应用。但是许多大坝都是建在地震多发和高烈度地区，并且坝还要承受重力、水压力等长期载荷的作用，为确保工程和人民生命财产在偶发地震载荷作用下的安全，需对大坝做抗震安全分析。重力坝抗震性能分析一般分以下五个步骤：1、创建物理环境2、建立模型和划分网格3、施加边界条件和载荷（地震荷载）4、求解5、后处理（查看计算结果）5.2.1创建物理环境在定义坝体抗震性能分析问题的物理环境时，进入ANSYS前处理器，建立这个坝体抗震性能分析的数学仿真模型。按照以下几个步骤来建立物理环境：1、设置GUT菜单过滤如果你希望通过GUI路径来运行ANSYS，当ANSYS被激活后第一件要做的事情就是选择菜单路径：MainMenuPreferences，执行上述命令后，弹出一个如图5-1所示的对话框出现后，选择Structural。这样ANSYS会根据你所选择的参数来对GUI图形界面进行过滤，选择Structural以便在进行坝体抗震性能分析时过滤掉一些不必要的菜单及相应图形界面。2、定义分析标题（／TITLE）在进行分析前，可以给你所要进行的分析起一个能够代表所分析内容的标题，比如“DamstabilityAnalysis”，以便能够从标题上与其他相似物理几何模型区别。用下列方法定义分析标题。命令：／TITLE第5章ANSYS水利工程应用实例分析GUI：UtilityMenuFileChangeTitle3、说明单元类型及其选项（KEYOPT选项）与ANSYS的其他分析一样，也要进行相应的单元选择。ANSYS软件提供了100种以上的单元类型，可以用来模拟工程中的各种结构和材料，各种不同的单元组合在一起，成为具体的物理问题的抽象模型。坝体用PLANE42单元来模拟。大多数单元类型都有关键选项（KEYOPTS），这些选项用以修正单元特性。例如，PLANE42有如下KEYOPTS：KEYOPT(2)包含或抑制过大位移设置KEYOPT(3)平面应力、轴对称、平面应变或考虑厚度的平面应力设置KEYOPT(5)解输出控制设置单元以及其关键选项的方式如下：命令：ETKEYOPTGUI：MainMenuPreprocessorElementTypeAdd/Edit/Delete图5-1GUI图形界面过滤4定义单位结构分析只有时间单位、长度单位和质量单位三个基本单位，则所有输入的数据都应当是这三个单位组成的表达方式。如标准国际单位制下，时间是秒（s），长度是米（m），质量是千克（kg），则导出力的单位是kg•m/s2（相当于牛顿N），材料的弹性模量单位是kg/m•s2（相当于帕Pa）。命令：/UNITS5、定义材料属性ANSYS10.0土木工程应用实例分析大多数单元类型在进行程序分析时都需要指定材料特性，ANSYS程序可方便地定义各种材料的特性，如结构材料属性参数、热性能参数、流体性能参数和电磁性能参数等。ANSYS程序可定义的材料特性有以下三种：（1）线性或非线性。（2）各向同性、正交异性或非弹性。（3）随温度变化或不随温度变化。因为进行坝体抗震性能分析时，ANSYS默认谱分析将忽略材料非线性，因此，坝体抗震性能分析模型采用弹性模型，只需要定义坝体材料属性中：容重、弹性模量、泊松比。命令：MPGUI：MainMenuPreprocessorMaterialPropsMaterialModels或MainMenuSolutionLoadStepOptsOtherChangeMatPropsMaterialModels坝体静力分析时可考虑材料的非线性，但进行抗震性能分析时，需要将非线性参数：内摩擦角和内聚力删除。5.2.2建立模型和划分网格创建好物理环境，就可以建立模型。在进行坝体抗震性能分析时，需要建立模拟坝体的PLANE82单元。在建立好的模型指定特性（单元类型、选项和材料性质等）以后，就可以划分有限元网格了。通过GUI为模型中的各区赋予特性：1、选择MainMenuPreprocessorMeshingMeshAttributesPickedAreas2、点击模型中要选定的区域。3、在对话框中为所选定的区域说明材料号、实常数号、单元类型号和单元坐标系号。通过命令为模型中的各区赋予特性：ASEL（选择模型区域）MAT（说明材料号）TYPE（指定单元类型号）建立大坝模型时，坝体和地基是赋予不同材料属性，本文只进行坝体抗震性能分析。进行大坝3-D模拟分析时，用SOLID65模拟混凝土单元和SOLID45模拟岩石单元的，2-D模拟分析时，只需用一个PLANE42单元就可以。5.2.3施加约束和荷载在施加边界条件和荷载时，既可以给实体模型（关键点、线、面）也可以给有限元模型（节点和单元）施加边界条件和荷载。在求解时，ANSYS程序会自动将加到实体模型上的边界条件和载荷转递到有限元模型上。第5章ANSYS水利工程应用实例分析重力坝抗震性能分析中，主要是给坝体底部施加自由度约束。命令：D作用在重力坝上的荷载包含水压力、冰压力、泥沙压力、地震力及坝体自重荷载等。(1)自重荷载。由坝体体积和材料的容重算出。(2)静水压力。作用在坝面上的静水压力可根据静水力学原理计算，分为水平力及垂直力。水平力：2121HPsp(5-1)垂直力：2121nHPcz(5-2)式中：为水容重，H1代表上游水深，n代表上游坝面坡度系数。同理可求得下游坝面的总静水压力的水平向分力及垂直分力。(3)扬压力。重力坝坝体混凝土或浆砌石砌体不是绝对不透水的，它们的表面及内部存在着无数微小的孔隙，坝基岩石本身孔隙虽然很少，但是也存在着节理、裂隙。重力坝建成挡水后，在上下游水位差的长期作用下，上游的水将通过这些孔隙及坝体和坝基的接触面、坝基的节理裂隙等向下游渗透，从而使得坝体内和坝底面产生渗透水压力。(4)动水压力。在溢流面上作用有动水压力，坝顶曲线和下游面直线段上的动水压力很小，可忽略不计。只计算反弧段上的动水压力。(5)冰压力。在寒冷地区水库表面冬季结成冰盖，当气温回升时，冰盖发生膨胀，因而对挡水建筑物上游面产生冰压力。(6)泥沙压力。水库蓄水后，入库水流挟带泥沙，逐年淤积在坝前，对坝面产生泥沙压力。(7)浪压力。浪压力与风速和水库吹程有关，但在荷载中所占比重较小，通常忽略。(8)地震荷载。主要是由建筑物质量引起的地震惯性力、地震动水压力和动土压力。5.2.4求解1.静力求解首先对重力坝进行静力求解：在ANSYS程序根据现有选项的设置，从数据库获取模型和载荷信息并进行计算求解，将结果数据写入到结果文件和数据库中。得到坝体在静力荷载作用下的位移场与应力，了解坝体在设计条件下的工作形态，对混凝土重力坝方案的可靠性进行评价评价大坝命令：SOLVEGUI：MainMenuSolutionSolveCurrentLS2.动力分析求解由于地震时的地面运动以水平方向为主，在地震力作用下结构的振动也以水平振动为主，ANSYS10.0土木工程应用实例分析故本次分析只考虑了水平方向的地震载荷的作用。对重力坝抗震性能计算分析可以采用以下几种方法：1)拟静力法拟静力法是一种把地震的影响用一种折算的静载荷来表示，求出这种地震荷载后，按照常规的静力法进行坝体的各项应力、位移的抗震分析方法。它是假定地震时与地面加速度相同的加速度作用在坝体各部位，求出地震时的惯性力，然后根据惯性力来评价大坝的安全性。根据拟静力分析方法，大坝的水平地震惯性力可简化为：FWCKQZHH(5-3)式中：HK——水平向地震系数，为地面最大水平加速度代表值与重力加速度的比值；ZC——综合影响系数，重力坝取为1/4；F——地震惯性力系数；W——产生惯性力的建筑物的总重量。采用拟静力法计算重力坝的地震作用效应时，水深h处的地震动水压力的代表值的计算：0)()(HhahPWhW(5-4)式中：)(hPW——作用在直立迎水坝面水深h处的动水压力代表值；ha——水平向地震加速度的代表值，地震烈度为8时对应的值是0.2g；——地震作用效应折减系数，除另有规定外，取0.25；)(h——水深h处的地震动水压力分布系数；W——水体质量密度的标准值；0H——水总深度。与水平面夹角为的倾斜迎水坝面，按公式(5-4)计算的动水压力代表值乘以折减系数：90/C(5-5)2)反应谱分析法反应谱分析法是以单质点弹性体系在实际地震过程中的反应为基础，来进行结构反应的分析，它通过反应谱巧妙地将动力问题静力化，使得复杂的结构地震反应计算变得简