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

1、第5章ANSYS水利工程应用实例分析水利工程ANSYS重力坝抗震性能分析步骤ANSYS重力坝抗震性能用实例分析(fefJ本章典型效果图水利工程概述虽然我国水利资源非常丰富，但河流在地区和时间分配上很不均衡，许多地区在枯水季节容易出现干早，而在洪水季节又往往由于水量过多而形成洪涝灾害。为了解决这一矛盾，人们 修建了许多水利工程来达到防洪、灌溉、发电、供水、航运等目的，促进国民经济建设的发展。水利工程中各种建筑物按其在水利枢纽中所起的作用，可以分为以下几类 :(1) 挡水建筑物用以拦截河流，形成水库，如各种坝和水闸以及抵御洪水所用的堤防等。(2) 泄水建筑物用以宣泄水库 或渠道 )在洪水期间或其它

2、情况下的多余水量， 以保证坝 (或 渠道)的安全，如各种溢流坝、溢流道、泄洪隧道和泄洪涵管等。(3) 输水建筑物为灌溉、 发电或供水， 从水库 ( 或河道) 向库外 (或下游 ) 输水用的建筑物， 如 引水隧道、引水涵管、渠道和渡槽等。(4) 取水建筑物是输水建筑物的首部建筑， 如为灌溉、 发电、供水而建的进水闸、 扬水站等。(5) 整治建筑物用以调整水流与河床、 河岸的相互作用以及防护水库、 湖泊中的波浪和水流 对岸坡的冲刷，如丁坝、顺坝、导流堤、护底和护岸等。由于破坏后果的灾难性，大型水利工程建设的首要目标是安全可靠，其次才是经济合理。 所以说研究大坝等水工建筑物的安全分析、评价和监控，是

3、工程技术人员需要解决的课题，正 确分析大坝性态已经成为当务之急。当前对各种水利工程评价主要采用有限元分析方法，借助各种有限元软件对这些水利工程建筑物进行安全评价，其中应用比较广泛的是 ANSY软件。目前，ANSY软件在水利工程中主要应用以下几个方面：(1) 应用各种坝体工程的设计和施工利用ANSY软件，模拟各种坝体施工过程以及坝体在使用阶段受到各种载荷(如水位变化对坝体的压力、地震荷载等)下结构的安全性能进行评价，模拟坝体的温度场和应力场，借助模 拟结果修改设计或对坝体采取加固措施。(2) 应用于各种引水隧道、引水涵管等设计和施工利用ANSY软件，模拟这些工程开挖、支护、浇注、回填过程，分析结

4、构在载荷作用下的变形情况、结构的安全可靠度，以及衬砌支护结构在水压、温度发生变化后产生的变形情况和结 构内力，依靠ANSY模拟结果对结构安全性进行评价。(3) 应用于各种水库闸门的设计和施工水库闸门在上游水作用下将发生弯曲、扭转、剪切和拉压等组合变形，利用 ANSY中的 SHELL6单元来模拟闸门，禾U用大型结构有限元分析程序 ANSYS对闸门结构进行三维有限元分 析，根据分析结果进行强度校核。ANSYS 重力坝抗震性能分析步骤重力坝是一种古老而重要的坝型，主要依靠坝体自身重力来维持坝身的稳定。岩基上重力坝的基本剖面呈三角形，上游面通常是垂直的或者稍倾向下游的三角形断面。 重力坝的具有很多优点

5、：(1)安全可靠。但剖面尺寸较大，抵抗水的渗漏，洪水漫顶，地震或战争破坏的能力都比较 强，因而失事率较低。(2)对地形、地质条件适应性强，坝体作用于地基面上的压应力不高，所以对地质条件的要 求也较低，低坝甚至可修建在土基上。(3)枢纽泄洪容易解决，便于枢纽布置。(4)施工方便，便于机械化施工。(5)结构作用明确，应力计算和稳定计算比较简单。 鉴于重力坝有如此多优点，所以它得到了广泛应用。 但是许多大坝都是建在地震多发和高烈度地区，并且坝还要承受重力、水压力等长期载荷 的作用， 为确保工程和人民生命财产在偶发地震载荷作用下的安全，需对大坝做抗震安全分析。重力坝抗震性能分析一般分以下五个步骤：1、

6、创建物理环境2、建立模型和划分网格3、施加边界条件和载荷(地震荷载)4、求解5、后处理(查看计算结果)创建物理环境在定义坝体抗震性能分析问题的物理环境时，进入 ANSYS 前处理器，建立这个坝体抗震 性能分析的数学仿真模型。按照以下几个步骤来建立物理环境：1、 设置 GUT 菜单过滤如果你希望通过 GUI 路径来运行 ANSYS ，当 ANSYS 被激活后第一件要做的事情就是选 择菜单路径：Main MenuPreferences,执行上述命令后，弹出一个如图5-1所示的对话框出现后，选择Structuralo这样ANSYS会根据你所选择的参数来对GUI图形界面进行过滤，选择Structura

7、l 以便在进行坝体抗震性能分析时过滤掉一些不必要的菜单及相应图形界面。2、 定义分析标题( TITLE )在进行分析前，可以给你所要进行的分析起一个能够代表所分析内容的标题，比如“Damstability Analysis ”,以便能够从标题上与其他相似物理几何模型区别。用下列方法定义分析标题。命令： TITLEGUI ： Utility MenuFileChange Title3、 说明单元类型及其选项( KEYOPT 选项)与 ANSYS 的其他分析一样， 也要进行相应的单元选择。 ANSYS 软件提供了 100种以上的 单元类型，可以用来模拟工程中的各种结构和材料，各种不同的单元组合在一

8、起，成为具体的 物理问题的抽象模型。坝体用 PLANE42 单元来模拟。大多数单元类型都有关键选项 (KEYOPTS ),这些选项用以修正单元特性。 例如，PLANE42 有如下 KEYOPTS ：KEY0PT（2）包含或抑制过大位移设置KEY0PT（3）平面应力、轴对称、平面应变或考虑厚度的平面应力设置KEYOPT （5） 解输出控制设置单元以及其关键选项的方式如下：命令：ETKEYOPTGUI : Main MenuPreprocessorElement TypeAdd/Edit/Delete图5-1 GUI图形界面过滤4定义单位结构分析只有时间单位、长度单位和质量单位三个基本单位，则所有

9、输入的数据都应当是这三个单位组成的表达方式。如标准国际单位制下，时间是秒（s）,长度是米（m），质量是千克（kg），则导出力的单位是kg?m/s2 （相当于牛顿N），材料的弹性模量单位是 kg/m?s2 （相当 于帕Pa）。命令：/UNITS5、定义材料属性大多数单元类型在进行程序分析时都需要指定材料特性，ANSYS程序可方便地定义各种材料的特性，如结构材料属性参数、热性能参数、流体性能参数和电磁性能参数等。ANSYS程序可定义的材料特性有以下三种：（1）线性或非线性。（2）各向同性、正交异性或非弹性。（3）随温度变化或不随温度变化。因为进行坝体抗震性能分析时，ANSYS默认谱分析将忽略材料非

10、线性，因此，坝体抗震性能分析模型采用弹性模型，只需要定义坝体材料属性中：容重、弹性模量、泊松比。命令：MPGUI: Main MenuPreprocessorMaterial PropsMaterial Models或 Mai n Men uSolutio nLoad Step OptsOtherCha nge Mat PropsMaterial Models坝体静力分析时可考虑材料的非线性，但进行抗震性能分析时，需要将非线性参数：内摩 擦角和内聚力删除。建立模型和划分网格创建好物理环境，就可以建立模型。在进行坝体抗震性能分析时，需要建立模拟坝体的 PLANE82单元。在建立好的模型指定特性（

11、单元类型、选项和材料性质等）以后，就可以划分 有限元网格了。通过GUI为模型中的各区赋予特性：1、选择 Main?Menu PreprocessorMeshing Mesh Attributes Picked Areas2、点击模型中要选定的区域。3、在对话框中为所选定的区域说明材料号、实常数号、单元类型号和单元坐标系号。通过命令为模型中的各区赋予特性：ASEL （选择模型区域）MAT （说明材料号）TYPE （指定单元类型号）建立大坝模型时，坝体和地基是赋予不同材料属性，本文只进行坝体抗震性能分析。进行大坝3-D模拟分析时，用SOLID65模拟混凝土单元和SOLID45模拟岩石单元的，2-D

12、 模拟分析时，只需用一个 PLANE42单元就可以。施加约束和荷载在施加边界条件和荷载时， 既可以给实体模型（关键点、线、面）也可以给有限元模型（节 点和单元）施加边界条件和荷载。在求解时，ANSYS程序会自动将加到实体模型上的边界条件 和载荷转递到有限元模型上。重力坝抗震性能分析中，主要是给坝体底部施加自由度约束。命令：D作用在重力坝上的荷载包含水压力、冰压力、泥沙压力、地震力及坝体自重荷载等。（1）自重荷载。由坝体体积和材料的容重算出。（2）静水压力。作用在坝面上的静水压力可根据静水力学原理计算，分为水平力及垂直力。1 2水平力：Psp - H!（5-1）垂直力：Pcz 1 nH；（5-2

13、）2式中： 为水容重，H1代表上游水深，n代表上游坝面坡度系数。同理可求得下游坝面的总静水压力的水平向分力及垂直分力。（3）扬压力。重力坝坝体混凝土或浆砌石砌体不是绝对不透水的，它们的表面及内部存在 着无数微小的孔隙，坝基岩石本身孔隙虽然很少，但是也存在着节理、裂隙。重力坝建成挡水PW (h) ah (h) W H 0( 5-4 )后，在上下游水位差的长期作用下，上游的水将通过这些孔隙及坝体和坝基的接触面、坝基的 节理裂隙等向下游渗透，从而使得坝体内和坝底面产生渗透水压力。(4)动水压力。在溢流面上作用有动水压力，坝顶曲线和下游面直线段上的动水压力很小， 可忽略不计。只计算反弧段上的动水压力。

14、(5)冰压力。在寒冷地区水库表面冬季结成冰盖，当气温回升时，冰盖发生膨胀，因而对 挡水建筑物上游面产生冰压力。(6)泥沙压力。水库蓄水后，入库水流挟带泥沙，逐年淤积在坝前，对坝面产生泥沙压力。(7)浪压力。浪压力与风速和水库吹程有关，但在荷载中所占比重较小，通常忽略。(8 )地震荷载。主要是由建筑物质量引起的地震惯性力、地震动水压力和动土压力。求解1. 静力求解首先对重力坝进行静力求解：在 ANSYS 程序根据现有选项的设置，从数据库获取模型和 载荷信息并进行计算求解，将结果数据写入到结果文件和数据库中。得到坝体在静力荷载作用 下的位移场与应力，了解坝体在设计条件下的工作形态，对混凝土重力坝方

15、案的可靠性进行评 价评价大坝命令： SOLVEGUI ：Main MenuSolutionSolve Current LS2.动力分析求解 由于地震时的地面运动以水平方向为主，在地震力作用下结构的振动也以水平振动为主， 故本次分析只考虑了水平方向的地震载荷的作用。对重力坝抗震性能计算分析可以采用以下几 种方法：1)拟静力法 拟静力法是一种把地震的影响用一种折算的静载荷来表示，求出这种地震荷载后，按照常 规的静力法进行坝体的各项应力、位移的抗震分析方法。它是假定地震时与地面加速度相同的 加速度作用在坝体各部位，求出地震时的惯性力，然后根据惯性力来评价大坝的安全性。根据拟静力分析方法，大坝的水平地

16、震惯性力可简化为：QH KHCZFW( 5-3 )式中：KH 水平向地震系数，为地面最大水平加速度代表值与重力加速度的比值；CZ 综合影响系数，重力坝取为 1/4；F 地震惯性力系数；W 产生惯性力的建筑物的总重量。采用拟静力法计算重力坝的地震作用效应时，水深h处的地震动水压力的代表值的计算：、Hoh110T( max 1)，.0T 0.1max，(5-5 )(5-6)(5-6 )(5-7)Tg T.式中：FW(h)作用在直立迎水坝面水深 h处的动水压力代表值；ah 水平向地震加速度的代表值，地震烈度为8时对应的值是；地震作用效应折减系数，除另有规定外，取；(h)水深h处的地震动水压力分布系数

17、；W水体质量密度的标准值；H 0 水总深度。与水平面夹角为 的倾斜迎水坝面，按公式(5-4)计算的动水压力代表值乘以折减系数:C /902)反应谱分析法反应谱分析法是以单质点弹性体系在实际地震过程中的反应为基础，来进行结构反应的分 析，它通过反应谱巧妙地将动力问题静力化，使得复杂的结构地震反应计算变得简单易行。按 照这一理论，应用地震谱曲线，就可以按照实际地面运动来计算建筑物的反应。反应谱是单点 弹性体系对于实际地面运动的最大反应和体系自振周期的函数关系。对于复杂的结构可以简化 为若干振型的叠加，每个振型又可转化为一个单质点来考虑。使用已经确定的设计反应谱计算 重力坝在地震作用下的反应，就归结

18、为寻求坝体的自振特性。地震产生的破坏，与受力大小和受频谱的最大振动的持续时间的都有关系。在进行谱分析 计算前，首先要计算大坝的自振特性。模态分析用于确定结构的振动特性，即结构的固有频率 和振型，它们是结构承受动态荷载设计中的重要参数，也是更详细的动力分析的基础。模态分析计算中采用了子空间迭代法提取模态。水深h处的地震动水压力的作用按公式转化为相应的坝面附加质量。Pw(h)根据如图5-2所示的大坝设计的反应谱曲线图，可得大坝反应谱曲线方程:0.1 T TgzT0 0.9()求得图5-2大坝设计反应谱本次重力坝抗震性能分析中，max取值为2, Tg取值为，特征周期To取值为。3）时程分析法时程分析

19、方法是将地震动记录或人工波作用在结构上直接对结构运动方程进行积分,结构任意时刻地震反应的分析方法，所以动态时程分析方法也称为直接积分法。本次大坝抗震性能实例分析采用反应谱分析方法。后处理后处理的目的是以图和表的形式描述计算结果。对于大坝抗震性能分析中，进入后处理器 后，查看大坝变形图和节点的位移和应力。通过研究大坝的变形、位移和应力情况，来综合判 断大坝的抗震性能及安全性能。命令：/ POST1GUI: Main MenuGeneral Postproc首先查看大坝静力分析求解结果，再查看大坝动力分析求解结果。ANSYS重力坝抗震性能实例分析实例介绍图5-3重力坝断面结构实例选取应用非常广泛的

20、重力坝，断面结构如图5-3所示。坝高120米，坝底宽为76米，坝顶为10米，上游坝面坡度和下游坝面坡度如图中所示。因为重力坝结构比较简单，垂直于长度方向的断面结构受力分布情况也基本相同，并且大 坝的纵向长度远大于其横断面，因此大坝抗震性能分析选用单位断面进行平面应变分析是可行 的。大坝抗震性能分析的计算条件如下：1） 假设大坝的基础是嵌入到基岩中，地基是刚性的。32） 大坝采用的材料参数为：弹性模量 E=35GPa,泊松比 =，容重 =25KN/m 。3） 计算分析大坝水位为120米。34） 水的质量密度1000kg/ m。5）大坝设防地震烈度为8，水平方向地震加速度值为。GUI操作方法创建物

21、理环境1）在【开始】菜单中依次选取【所有程序】/【】/【ANSYS Product Launcher】，得到“ Product Laun cher ”对话框。2）选中【File Management】，在 “Working Directory ” 栏输入工作目录 “D:ansysexample5-1”， 在Job Name栏输入文件名Dam ”。3） 单击“ RUN ”按钮，进入的GUI操作界面。4）过滤图形界面： Main MenuPreferences,弹出“ Preferences for GUI Filtering”对话框， 选中“ Structural来对后面的分析进行菜单及相应的图

22、形界面过滤。5） 定义工作标题：Utility Menu File Change Title，在弹出的对话框中输入“ Dam seismic Analysis”，单击“ OK”，如图 5-4。6）定义单元类型:a.定义 PLANE42 单元：Ma in Me nu Preprocessor Eleme nt Type Add/Edit/Delete，弹出一 个单元类型对话框，单击“ Add”按钮。弹出如图5-5所示对话框。在该对话框左面滚动栏中选 择“ Solid ”，在右边的滚动栏中选择“ Quad 4node 42”单击“ Apply ”，就定义了“ PLANE42 ” 单元。IMH Ib

23、ra评 of 匚loment Type雪Cnly才 fib-w t 可e shzi jmLbrary 口 f E eirnent Type=Etomerit 対 pc rwltrurce nu 恤 or图5-5定义PLANE42单元对话框b.设定 PLANE42 单元选项：Main Menu Preprocessor Element Type Add/Edit/Delete，弹 出一个单元类型对话框， 选中“Type 2 PLANE42 ”，单击“Options”按钮，弹出一个“PLANE42 element Type options”对话框，如图 5-6所示。在“ Element behav

24、ior K3 ”栏后面的下拉菜单中图5-4定义工作标题2L选取“ Plane strain”其它栏后面的下拉菜单采用ANSYS默认设置就可以，单击“ 0K”按钮。图5-6 PLANE42单元库类型选项对话框通过设置PLANE42单元选项“ K3 ”为“ Plane strain”来设定本实例分析采取平面应变模 型进行分析。因为大坝是纵向很长的实体，故计算模型可以简化为平面应变问题。7）定义材料属性执行 Main MenuPreprocessorMaterial PropsMaterial Models，弹出 “ Define Material Model Behavior ”对话框，如图 5-

25、7所示。图5-7定义材料本构模型对话框在图5-7中右边栏中连续双击“ Structural Linear Elasticlsotropic”后，又弹出如图 5-8 所示“ Lin ear Isotropic Properties for Material Number 1 ” 对话框，在该对话框中“ EX ”后面的输 入栏输入“”，在“ PRXY ”后面的输入栏输入“”，单击“0K”。再在选中“ Density ”并双击， 弹出如图5-9所示“ Density for Material Number 1 对话框，在“ DENS”后面的栏中输入边坡土 体材料的密度“ 2500”单击“ OK”按钮

26、。风| igj前叩PropeTies foivtohAba鬪 MffMEeaPIgEBIL izar ：50h卩 t later tai Pr卯白 Des fcr 甘出孑 iall Nu nhs-1Density for Maoeral Mumber 1Tlidc -eirprratLrs IteeteTEmerstirenr图5-8线弹性材料模型对话框TerrperatLresCEhG2500Add Tenperatjn& Delete enperaLirelOKCancel图5-9材料密度输入对话框建立模型和划分网格1）创建大坝线模型a.输入关键点：Main Me nu Preproces

27、sorModel in g CreateKeypoi ntsln Active CS，弹出 “ Creae Keypoints in Active Cooedinate System”对话框，如图 5-10所示。在NPT keypoint number ” 栏后面输入 “1”，在“X，Y，Z Location in active CS ” 栏后面输入“（0, 0，0）”，单击“ Apply ” 按钮，这样就创建了关键点1。再依次重复在“ NPT keypoint number”栏后面输入“ 2、3、4、5”，在对应“ X，Y，Z Location in active CS ” 栏后面输入“（7

28、6，0，0）、（，0）、（，120，0）、（，120，0），最后单击“ OK”按钮。图5-10在当前坐标系创建关键点对话框b.创建坝体线模型：Main Menu Preprocessor Modeling CreateLinesStraight line，弹出“ Creae straight lines对话框，用鼠标依次点击关键点 1、2，单击“ Apply ”按钮，这样就创建 了直接L1，同样分别连接关键点“ 2、3”，“3、4”，“4、5”，“5、1”，最后单击“ OK ”按钮， 就得到坝体线模型，如图 5-11所示。nompcrEXpfycrGraprL1UL3L iuL rrjdiD砂

29、 sex图5-11 坝体线模型2）创建坝体面模型a.打开关键点编号显示：Utility MenuPlotCtrls Numbering ,弹出 “ Plot Numbering Controls 对话框，如图5-12所示。选中Keypoint Numbers”选项，后面的文字由off”变为on” 单击“ OK”关闭窗口。jjjplotNumbering ControlsIrPNlJM Pin t Ncimzer nq CcnlrDl=KT DycohTntmtcrsP iDnLIFE L he rumbenr ZIP=AREAriurrtasrftr errvoiu vou 叫rreer*r

30、offNODE hJode numbersr offEbm / Attrib mfrterlrg|MO nurbiirirgTABTJ Tankf Ma -issr沖SVALEFWlF 回iMr DfF| kumoEr rq Ehoiwi wrth| 厂nlnrs.也RnhprsIPLSTJ Pellet gen 俳询出|RsqbtCK Afpk心紳 |help图5-12打开关键点编号对话框b.创建坝体面模型：Main MenuPreprocessorModelingCreateAreasArbitraryThroughKPs，弹出一个“ Create Area by Keypoints”对话

31、框，在图形中选取关键点 1、2、3、4和5,点击 “Apply ”按钮，就得到坝体模型的面模型，如图5-13所示。L1KE.3TaIz ETO.AIM+YW 3 juaae Ji domic5图5-13坝体面模型3）划分坝体单元网格a.设置网格份数：Mai n Men uPreprocessorMeshi ngSize Cn trlsMa nualSizeLayersPicked Lines,弹出一个Set Layer Controls”对话框，如图5-14所示，用鼠标选取线 L1，单击0K” 按钮。弹出一个Area Layer Mesh Control on Picked lines ” 对

32、话框，如图 5-15 所示，在No ofline divisio n ”栏后面输入“ 20”，单击“OK”按钮。Sot Layer ControIs住广ITiijjlckf*S云阳営1.龟厂BoxLTlmyz厂Ciiicle广 LeapC QUiTltr0MSTC-S w iwi=01出竝3&1L liiw-仔a f JSKIEL TIzLiz* Ha玉* 工口匚图5-14选取线对话框相同方法设置线L2分割份数为图5-15设置网格分数对话框32;设置线L3、L4和L5线的分割份数分别为6、 4、 40。b.划分单元网格：Main MenuPreprocessorMeshingMeshAreaF

33、ree，弹出一个拾取面积 对话框，拾取图形中面，单击拾取框上的“0K ”按钮，得到坝体模型单元网格，如图5-16所示。4）保存坝体单元网格Utility MenuFile Save as，弹出一个“ Save Database 对话框，在“ Save Database to 下 面输入栏中输入文件名“”，单击“ 0K”。图5-16坝体单元网格图5-17给坝体底部施加位移约束对话框施加约束和荷载1）给坝体模型底部施加位移约束执行 Main MenuSolutionDefine LoadsApplyStructuralDisplacementon Nodes ，弹出在 节点上施加位移约束对话框，用

34、鼠标选取隧道模型两侧边界上所有节点，单击“0K ”按钮。弹出Apply U,ROT on Nodes ” 对话框，如图 5-17 所示，在DOFS to be constrained”栏后面中 选取 “ ALL DOF ”，在 “ Apply as ” 栏后面的下拉菜单中选取 “ Con sta nt value,在 “ Displaceme nt value”栏后面输入“ 0”值，然后单击“ OK ”按钮。2）施加重力加速度： Main Me nuSolutio nDefi ne LoadsApplyStructurall nertiaGravity ， 弹出“ Apply（Gravitat

35、ional）Acceleration ” 对话框，如图 5-18 所示。在“ Global Cartesian Y-comp” 栏后面输入重力加速度值“”就可以，单击“OK ”按钮，就完成了重力加速度的施加。hlAppI#Aoz巳le吕tioriACEL Apply (Gravitational) AccelerationAC ELX Qlobal Cartes ian X-compACELY Global Cartesian Y-compACELZ Glob日I Cartesian Z-cornp图5-18施加重力加速度对话框3）施加水压力载荷： Main Me nuSolutio nDef

36、i ne LoadsApplyStructuralStructureo n Lin es ， 弹出一个对话框，用鼠标选中线L5，单击“OK”。弹出“Apply PRES on lines”对话框，如图5-19 所示。分别输入数据“ 0”和“ 1101370”，单击“ OK”按钮，就完成了水压力载荷的施加。求解1静力分析求解1）求解设置a.指定求解类型：Main MenuSolutionAnalysis TypeNew Analysis，弹出一个如图 5-20 所图5-19施加水压力载荷对话框本次加的荷载是水深为120米时作用在坝上的水压力，迎水面波度是87度。示对话框，在Type of ana

37、lysis”栏后面选中Static,单击OK”按钮。图5-20指定求解类型对话框b.设置载荷步： Main MenuPreprocessorLoadsAnalysis TypeSoln Controls，弹出一个Solution Controls对话框,用鼠标单击 “Basic选项，如图 5-21 所示，在Number of Substeps” 栏后面输入5,在Max no. of substeps” 栏后面输入100”，在Min no. of substeps”栏后 面输入“ 1 ”，单击“ OK ”按钮。图5-21设置载荷步对话框c. 设置线性搜索： Main MenuSolutionAn

38、alysis TypeSoln Controls，弹出一个“ Solution Controls”对话框，用鼠标单击“ Nonlinear”选项，如图5-22所示，在“ Line search栏后面下 拉菜单选中“ ON ”，单击“ OK ”。2）静力求解a. 求解：Main MenuSolutionSolveCurrent LS，弹出一个求解选项信息和一个当前求解载 荷步对话框，检查信息无错误后，单击“OK ”，开始求解运算，直到出现一个“ Solution is done” 的提示栏，表示求解结束。b. 保存求解结果；Utility Menu File Save as ,弹出一个 “Sav

39、e Database 对话框，在“ Save Database td下面输入栏中输入文件名“”，单击“ OK”。图5-22设置线性搜索对话框图5-23设置模态分析选项对话框2抗震性能分析求解1）模态分析求解a. 设置分析类型：Main MenuSolutionAnalysis TypeNew Analysis，弹出一个如图 5-20所示 对话框，在Type of analysis”栏后面选中Modal”，单击OK”按钮。b. 设置模态分析选项： Main MenuSolutionAnalysis TypeAnalysis Options，弹出一个如图5-23 所示对话框，在Mode extra

40、ction method ” 栏后面选中Sun space，在No. of modes to be extract”栏后面输入18”，在Expand mode shapes后面小方框用鼠标选中，单击 OK”按 钮，又弹出一个Sun space Modal An alysis”对话框，如图5-24所示，按图中设置后，单击OK ” 按钮。图5-24模态分析求解时子空间设置对话框图5-25模态求解选项信息c. 模态分析求解：Main MenuSolutionSolveCurrent LS，弹出一个模态求解选项信息（如图5-25所示）和一个当前求解载荷步对话框，检查信息无错误后，单击“0K”开始求解运

41、算，直到出现一个Solution is done 的提示栏，表示求解结束。d. 保存求解结果；Utility Menu File Save as，弹出一个 “Save Database 对话框，在“Save Database td下面输入栏中输入文件名“”，单击“ OK”。11反应谱值e调出模态分析各阶频率：Main Me nu Ge neral Postproc Read Summary，弹出如图5-26所示对话框。动力求解和静力求解的模型相同，约束条件也相同。谱分析时，ANSYS忽略材料非线性。调出模态分析各阶频率是为后面求解反应谱值。2）反应谱分析求解a.求出反应谱值：由图5-26中前1

42、8阶频率值f,可以算出对应的周期T，再根据大坝反应谱曲线方程（式5-7），可以计算出前10阶的反应谱值，见表5-1。表5-1大坝动力计算前10阶振动频率及反应谱值振型 振动频率1/S 振动周期T12345678910图5-26模态分析各阶频率12131415161718b.设置反应谱分析求解选项设置分析类型：Main MenuSolutionAnalysis TypeNew Analysis，弹出一个如图 5-27所示对话框，在Type of analysis”栏后面选中Spectrum”，单击OK”按钮。图5-27定义反应谱分析对话框设置反应谱分析选项：Main MenuSolutionAn

43、alysis TypeAnalysis Options，弹出一个Spectrum Analysis对话框，如图5-28所示，在Type of spectrum后面栏中选取 “Single-pt resp”， 在No. of modes for solu ”后面输入10”，在Calculate elem stress?后面选中Yes,单 击“OK”按钮。设置反应谱单点分析选项：Main Me nuSolutio nLoad Step OptsSpectrumS in glePtSett in gs,弹出一个Setti ng for si ngle-poi nt Respo nse Spectru

44、m” 对话框，如图 5-29 所示。在“ Type of response spectrum栏后面下拉菜单选中“ Seismic accel,在“ SEDX, SEDY, SEDZ”栏后 面依次输入“ 0、1、0,单击“ OK ”按，钮。2J.SFCPT SpQJrurTl AraH腑 CfltfG Bpt|pe Tiipe nf Epeclr u -Tt 卜# of rrcdPi forfclcz CaliciblE tlsn 曲w貂炉(for P .E.D. only)OK ICancel ajB U1 IUB EB.EBAB -kilfi-ir i 心耳 tt vTn vaennBrar

45、 tpt 广 W|L|ll3piT舉冃YTH r D.D A.M.r PSSJD lift-P YKO ICSEBI Ii* |p图5-29设置单点反应谱分析对话框定义反应谱分析频率表：Main MenuSolutionLoad Step OptsSpectrumSinglePtFreq Table ,弹出一个“ Frequency Table”对话框，如图5-30所示。根据表5-1依次输入大坝的前18阶振动频 率。定义反应谱值： Main MenuSolutionLoad Step OptsSpectrumSinglePtSpectr V alues，弹出 一个 “Spectrum Valu

46、es Damping Ration” 对话框，单击 “ OK” 按钮。弹出一个 “Spectrum Values 对话框，如图5-30所示。根据表5-1依次输入大坝的前18阶反应谱值。?注意，输入的振动频率必须按升序排列。? FREQ1必须大于零。图5-28设置反应谱分析对话框Xlrri h nprirlr im *ili *x 5印ESEftFDCH chwt刖(glctial CartesLi i)图5-30定义频率表对话框图5-31定义反应谱值对话框c. 反应谱分析求解：Main MenuSolutionSolveCurrent LS，弹出一个模态求解选项信息和一个当前求解载荷步对话框，

47、检查信息无错误后，单击“OK”，开始求解运算，直到出现一个a Solution is done”的提示栏，表示求解结束。d. 保存求解结果；Utility Menu File Save as，弹出一个 “Save Database 对话框，在“Save Database td下面输入栏中输入文件名“”，单击“ 0K”。3）模态扩展分析求解a. 设置分析类型： Main MenuSolutionAnalysis TypeNew Analysis，弹出一个“ NewAnalysis”对话框，在“ Type of analysis”栏后面选中“ Modal”，单击“ OK”按钮。b. 设置模态扩展分

48、析求解选项定义模态扩展分析： Main MenuSolutionAnalysis TypeExpansionPass ,弹出一个 “ ExpansionPass对话框，如图5-32所示，选中“ Expansion Pass选项，后面的文字由“ off”变为“ on” ,单击“ OK”关闭窗口。图5-32模态扩展对话框设置模态扩展分析： Main MenuPreprocessorLoadsLoad Step OptsExpansionPassSingle ExpandExpand Modes，弹出一个 “Expand Modes”对话框，如图5-33所示，在“No. of modes to ex

49、pand”栏后面输入“18”，其它如图中设置，单击“ OK”按钮。图5-33模态扩展分析设置对话框c. 模态扩展分析求解： Main MenuSolutionSolveCurrent LS，弹出一个模态求解选项信息和一个当前求解载荷步对话框，检查信息无错误后，单击“OK ”，开始求解运算，直到出现一个“ Solution is done”的提示栏，表示求解结束。d. 保存求解结果；Utility Menu File Save as，弹出一个 “Save Database 对话框，在“Save Database td下面输入栏中输入文件名“”，单击“ OK”。4）合并模态分析求解a. 设置分析类

50、型： Main MenuSolutionAnalysis TypeNew Analysis，弹出一个“ NewAnalysis” 对话框，在“ Type of analysis” 栏后面选中“ spectrum,单击“ OK” 按钮。b. 按平方和方根法进行组合： Main Men uSolutio nLoad Step Opts Spectrum Sin glePtMode Combine,弹出一个“ Mode Combination Methods” 对话框，如图 5-34所示，在“ Mode Combination Method ”栏后面下拉菜单中选取“ SRSS,在“ Signific

51、ation threshold ”栏后面输入“”，在“ Type of output”栏后面下拉菜单中选取“ Displacement”，单击“ OK”按钮。图5-34合并模态求解设置对话框c. 合并模态分析求解： Main MenuSolutionSolveCurrent LS，弹出一个模态求解选项信息和一个当前求解载荷步对话框，检查信息无错误后，单击“0K ”，开始求解运算，直到出现一个Solution is done”的提示栏，表示求解结束。d. 保存求解结果；Utility Menu File Save as，弹出一个 “Save Database 对话框，在“Save Databas

52、e tc”下面输入栏中输入文件名“ Dam- ”，单击“ 0K”。后处理对坝体进行了静力分析和动力分析，通过对重力坝静力有限元分析，可以知道坝体在静力载荷作用下的位移场和应力场，从而可以了解坝体的安全性能。通过对重力坝抗震性能有限元 分析，可以了解大坝的地震载荷作用下的动力响应特性，从而可以评价大坝在地震载荷作用下 的安全性能。1.静力分析求解结果a.读入静力分析求解结果数据：Utility Menuresume，弹出一个“ Resume Database对话框，选中刚才保存的文件“”，单击“ 0K”按钮。r tersh pj Jiy石 DeT * unLfcftrredL * UPtM Hc

53、fcJiJml图5-35绘制变形图对话框b.绘制坝体变形图： Main MenuGeneral PostprocPlot ResultsDeformed Shape,弹出一个“Plot Deform Shape” 对话框，如图 5-35 所示。选中“ Def +undeformed”，单击“ OK” 按钮， 得到坝体变形图，如图 5-36所示。图5-36坝体变形图5-37节点解云图绘制对话框c.显示坝体 x 方向位移云图：Ma in Men uGe neral PostprocPlot ResultsCo ntour PlotNodal Solu，弹出一个“ Con tour Nodal So

54、lution Data ”对话框，如图5-37所示，用鼠标依次点击“ Nodal Solution/DOF Solution/ X-Compoment of displacement ”，再单击OK” 按钮，就得到坝体 X 方向 位移云图，如图5-38所示。此时，坝体水平方向最大位移为，位置发生在坝顶。nrprijSOE- TDIaHVTkfcVB)Pffl啣tax -.duzzjJIB A ZDQT-.0011J-d.显示坝体 y 方向位移云图： Main MenuGeneral PostprocPlot ResultsContour PlotNodal Solu，弹出一个“ Con tou

55、r Nodal Solution Data ”对话框，如图5-37所示，用鼠标依次点击“ Nodal Solution/Stress/ Y-Compoment of displacement”，再单击OK” 按钮，就得到坝体 Y 方向位移云 图，如图5-39所示。此时，坝体垂直方向最大位移为，发生在坝腹处。-.Xupii& deLSLAb： IfALTSiB图5-39坝体y方向位移云图e. 显示坝体 x 方向应力云图： Ma in Men uGe neral PostprocPlot ResultsCo ntour PlotNodal Solu，弹出一个“ Con tour Nodal Sol

56、ution Data ”对话框，如图5-37所示，用鼠标依次点击“ Nodal Solution/stress/ x-Compoment of stress,再单击OK”按钮，就得到坝体 x方向应力云图，女口 图5-40所示。f.显示坝体 y 方向应力云图：Ma in Men uGe neral PostprocPlot ResultsCo ntour PlotNodalSolu，弹出一个“ Con tour Nodal Solution Data ”对话框，如图5-37所示，用鼠标依次点击“ Nodal Solution/stress/ y-Compoment of stress,再单击OK

57、”按钮，就得到坝体 x方向应力云图，女口 图5-41所示。图5-38 坝体x方向位移云图Mht MkmiBn* iTina fT型訂ICTMl5 DU 岀 sq -_rzs*i7 t emunWP 1 :皿 4 ZHnTHFIT图5-41坝体y方向应力云图/WTSYijn 4 1M7i JIRD-w 禹如戸-2fcZ7图5-42第一主应力云图g.显示坝体第一主应力云图：Mai n Me nuGe neral PostprocPlot ResultsCo ntour PlotNodalEu Lmlc JLNIR.llL+m-2au+m.1771-Mi?.MH-HTZlVKtffl.IVEbflF

58、*4GLLJfibMZ-iLSolu，弹出一个“ Con tour Nodal Solution Data ”对话框，如图5-37所示，用鼠标依次点击“ Nodal Solution/stress/ 1st principal stress,再单击OK”按钮，就得到坝体第一主应力云图， 如图5-42 所示。图5-40坝体x方向应力云图HHCNkL MUmM mA3 IB -5-图5-43坝体第1阶振型依次重复执行c步骤，就可以依次画出坝体第图5-44坝体第2阶振型3到第18阶振型，如图5-45图5-60所示。2.抗震性能分析求解结果1）绘制坝体振型图a. 读入模态分析分析求解结果数据： Uti

59、lity Menu resume，弹出一个Resume Database对话框，选中刚才保存的文件“”，单击“ 0K ”按钮。b. 绘制坝体第1阶震型读入第 1 阶数据：Main Me nuGe neral PostprocRead ResultsFirst set.绘制坝体第 1 阶震型：Main MenuGeneral PostprocPlot ResultsDeformed Shape，弹出一 个“ Plot Deform Shape对话框，选中“ Def +undeformed”，单击“ OK ”按钮，得到坝体第 1 阶震型，如图5-43所示。c.绘制坝体第2阶震型读入第 2 阶数据：

60、Main MenuGeneral PostprocRead ResultsNext set.绘制坝体第 2 阶震型：Main MenuGeneral PostprocPlot ResultsDeformed Shape，弹出一 个“ Plot Deform Shape对话框，选中“ Def +undeformed”，单击“ OK ”按钮，得到坝体第 2 阶震型，如图5-44所示。d.调出模态分析结果的固有频率和振型：Main Me nu Ge neral Postproc Read Summary,弹出如图5-26所示对话框。图5-49坝体第7阶振型图5-50坝体第8阶振型so - i nih