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

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

内容 提要

本章首先概述了水利工程 ANSYS应用，其次介绍了ANSYS重力坝抗震性能分析步骤，最后用实例详细介绍了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.2 ANSYS重力坝抗震性能分析步骤

重力坝是一种古老而重要的坝型，主要依靠坝体自身重力来维持坝身的稳定。岩基上重力坝的基本剖面呈三角形，上游面通常是垂直的或者稍倾向下游的三角形断面。

重力坝的具有很多优点：

(1)安全可靠。但剖面尺寸较大，抵抗水的渗漏，洪水漫顶，地震或战争破坏的能力都比较强，因而失事率较低。

(2)对地形、地质条件适应性强，坝体作用于地基面上的压应力不高，所以对地质条件的要求也较低，低坝甚至可修建在土基上。

(3)枢纽泄洪容易解决，便于枢纽布置。

(4)施工方便，便于机械化施工。

(5)结构作用明确，应力计算和稳定计算比较简单。 鉴于重力坝有如此多优点，所以它得到了广泛应用。

但是许多大坝都是建在地震多发和高烈度地区，并且坝还要承受重力、水压力等长期载荷的作用，为确保工程和人民生命财产在偶发地震载荷作用下的安全，需对大坝做抗震安全分析。

重力坝抗震性能分析一般分以下五个步骤：

1、创建物理环境

2、建立模型和划分网格

3、施加边界条件和载荷（地震荷载） 4、求解

5、后处理（查看计算结果）

5.2.1 创建物理环境

在定义坝体抗震性能分析问题的物理环境时，进入ANSYS前处理器，建立这个坝体抗震性能分析的数学仿真模型。按照以下几个步骤来建立物理环境：

1、 设置GUT菜单过滤

如果你希望通过GUI路径来运行ANSYS，当ANSYS被激活后第一件要做的事情就是选择菜单路径：Main Menu>Preferences，执行上述命令后，弹出一个如图5-1所示的对话框出现后，选择Structural。这样ANSYS会根据你所选择的参数来对GUI图形界面进行过滤，选择Structural以便在进行坝体抗震性能分析时过滤掉一些不必要的菜单及相应图形界面。

2、 定义分析标题（／TITLE）

在进行分析前，可以给你所要进行的分析起一个能够代表所分析内容的标题，比如“Dam stability Analysis”，以便能够从标题上与其他相似物理几何模型区别。用下列方法定义分析标题。

命令：／TITLE

第5章 ANSYS水利工程应用实例分析 GUI：Utility Menu>File>Change Title

3、 说明单元类型及其选项（KEYOPT选项）

与ANSYS的其他分析一样，也要进行相应的单元选择。ANSYS软件提供了100种以上的单元类型，可以用来模拟工程中的各种结构和材料，各种不同的单元组合在一起，成为具体的物理问题的抽象模型。坝体用PLANE42单元来模拟。

大多数单元类型都有关键选项（KEYOPTS），这些选项用以修正单元特性。例如，PLANE42有如下KEYOPTS：

KEYOPT(2) 包含或抑制过大位移设置

KEYOPT(3) 平面应力、轴对称、平面应变或考虑厚度的平面应力设置 KEYOPT(5) 解输出控制

设置单元以及其关键选项的方式如下： 命令：ET

KEYOPT

GUI：Main Menu> Preprocessor> Element Type> Add/Edit/Delete

图5-1 GUI图形界面过滤

4定义单位

结构分析只有时间单位、长度单位和质量单位三个基本单位，则所有输入的数据都应当是这三个单位组成的表达方式。如标准国际单位制下，时间是秒（s），长度是米（m），质量是千克（kg），则导出力的单位是kg?m/s2（相当于牛顿N），材料的弹性模量单位是kg/m?s2（相当于帕Pa）。

命令：/UNITS

5、定义材料属性

ANSYS10.0 土木工程应用实例分析 大多数单元类型在进行程序分析时都需要指定材料特性，ANSYS程序可方便地定义各种材料的特性，如结构材料属性参数、热性能参数、流体性能参数和电磁性能参数等。

ANSYS程序可定义的材料特性有以下三种： （1）线性或非线性。

（2）各向同性、正交异性或非弹性。 （3）随温度变化或不随温度变化。

因为进行坝体抗震性能分析时，ANSYS默认谱分析将忽略材料非线性，因此，坝体抗震性能分析模型采用弹性模型，只需要定义坝体材料属性中：容重、弹性模量、泊松比。

命令：MP

GUI：Main Menu>Preprocessor>Material Props>Material Models

或 Main Menu>Solution>Load Step Opts>Other>Change Mat Props>Material Models ? 坝体静力分析时可考虑材料的非线性，但进行抗震性能分析时，需要将非线性参数：内摩

擦角和内聚力删除。

5.2.2 建立模型和划分网格

创建好物理环境，就可以建立模型。在进行坝体抗震性能分析时，需要建立模拟坝体的PLANE82单元。在建立好的模型指定特性（单元类型、选项和材料性质等）以后，就可以划分有限元网格了。

通过GUI为模型中的各区赋予特性：

1、选择Main Menu> Preprocessor> Meshing> Mesh Attributes> Picked Areas 2、点击模型中要选定的区域。

3、在对话框中为所选定的区域说明材料号、实常数号、单元类型号和单元坐标系号。 通过命令为模型中的各区赋予特性： ASEL（选择模型区域） MAT （说明材料号） TYPE（指定单元类型号）

? 建立大坝模型时，坝体和地基是赋予不同材料属性，本文只进行坝体抗震性能分析。 ? 进行大坝3-D模拟分析时，用SOLID65模拟混凝土单元和SOLID45模拟岩石单元的，2-D

模拟分析时，只需用一个PLANE42单元就可以。

5.2.3 施加约束和荷载

在施加边界条件和荷载时，既可以给实体模型（关键点、线、面）也可以给有限元模型（节点和单元）施加边界条件和荷载。在求解时，ANSYS程序会自动将加到实体模型上的边界条件和载荷转递到有限元模型上。

第5章 ANSYS水利工程应用实例分析 重力坝抗震性能分析中，主要是给坝体底部施加自由度约束。

命令：D

作用在重力坝上的荷载包含水压力、冰压力、泥沙压力、地震力及坝体自重荷载等。 (1)自重荷载。由坝体体积和材料的容重算出。

(2)静水压力。作用在坝面上的静水压力可根据静水力学原理计算，分为水平力及垂直力。 水平力： Psp?

垂直力： Pcz?1212?H1 (5-1)

2?nH1 (5-2)

2式中：?为水容重，H1代表上游水深，n代表上游坝面坡度系数。

同理可求得下游坝面的总静水压力的水平向分力及垂直分力。

(3)扬压力。重力坝坝体混凝土或浆砌石砌体不是绝对不透水的，它们的表面及内部存在着无数微小的孔隙，坝基岩石本身孔隙虽然很少，但是也存在着节理、裂隙。重力坝建成挡水后，在上下游水位差的长期作用下，上游的水将通过这些孔隙及坝体和坝基的接触面、坝基的节

理裂隙等向下游渗透，从而使得坝体内和坝底面产生渗透水压力。

(4)动水压力。在溢流面上作用有动水压力，坝顶曲线和下游面直线段上的动水压力很小， 可忽略不计。只计算反弧段上的动水压力。

(5)冰压力。在寒冷地区水库表面冬季结成冰盖，当气温回升时，冰盖发生膨胀，因而对 挡水建筑物上游面产生冰压力。

(6)泥沙压力。水库蓄水后，入库水流挟带泥沙，逐年淤积在坝前，对坝面产生泥沙压力。

(7)浪压力。浪压力与风速和水库吹程有关，但在荷载中所占比重较小，通常忽略。 (8 )地震荷载。主要是由建筑物质量引起的地震惯性力、地震动水压力和动土压力。

5.2.4 求解

1. 静力求解

首先对重力坝进行静力求解：在ANSYS程序根据现有选项的设置，从数据库获取模型和载荷信息并进行计算求解，将结果数据写入到结果文件和数据库中。得到坝体在静力荷载作用下的位移场与应力，了解坝体在设计条件下的工作形态，对混凝土重力坝方案的可靠性进行评价评价大坝

命令：SOLVE

GUI：Main Menu>Solution>Solve> Current LS 2.动力分析求解

由于地震时的地面运动以水平方向为主，在地震力作用下结构的振动也以水平振动为主，

ANSYS10.0 土木工程应用实例分析 故本次分析只考虑了水平方向的地震载荷的作用。对重力坝抗震性能计算分析可以采用以下几种方法：

1)拟静力法

拟静力法是一种把地震的影响用一种折算的静载荷来表示，求出这种地震荷载后，按照常规的静力法进行坝体的各项应力、位移的抗震分析方法。它是假定地震时与地面加速度相同的加速度作用在坝体各部位，求出地震时的惯性力，然后根据惯性力来评价大坝的安全性。

根据拟静力分析方法，大坝的水平地震惯性力可简化为：

QH?KHCZFW ( 5-3 )

式中：

KH——水平向地震系数，为地面最大水平加速度代表值与重力加速度的比值； CZ——综合影响系数，重力坝取为1/4； F——地震惯性力系数；

W——产生惯性力的建筑物的总重量。

采用拟静力法计算重力坝的地震作用效应时，水深h处的地震动水压力的代表值的计算： PW(h)?ah??(h)?WH0 ( 5-4 ) 式中：

PW(h)——作用在直立迎水坝面水深h处的动水压力代表值；

ah——水平向地震加速度的代表值，地震烈度为8时对应的值是0.2g；

?——地震作用效应折减系数，除另有规定外，取0.25； ?(h)——水深h处的地震动水压力分布系数； ?W——水体质量密度的标准值；

H0——水总深度。

与水平面夹角为?的倾斜迎水坝面，按公式(5-4)计算的动水压力代表值乘以折减系数： ?C??/90 ( 5-5 ) 2)反应谱分析法

反应谱分析法是以单质点弹性体系在实际地震过程中的反应为基础，来进行结构反应的分析，它通过反应谱巧妙地将动力问题静力化，使得复杂的结构地震反应计算变得简单易行。按照这一理论，应用地震谱曲线，就可以按照实际地面运动来计算建筑物的反应。反应谱是单点弹性体系对于实际地面运动的最大反应和体系自振周期的函数关系。对于复杂的结构可以简化为若干振型的叠加，每个振型又可转化为一个单质点来考虑。使用已经确定的设计反应谱计算重力坝在地震作用下的反应，就归结为寻求坝体的自振特性。

地震产生的破坏，与受力大小和受频谱的最大振动的持续时间的都有关系。在进行谱分析计算前，首先要计算大坝的自振特性。模态分析用于确定结构的振动特性，即结构的固有频率和振型，它们是结构承受动态荷载设计中的重要参数，也是更详细的动力分析的基础。

第5章 ANSYS水利工程应用实例分析 模态分析计算中采用了子空间迭代法提取模态。水深h处的地震动水压力的作用按公式(5-6)转化为相应的坝面附加质量。

PW(h)?78ah?WH0h ( 5-6 )

根据如图5-2所示的大坝设计的反应谱曲线图，可得大坝反应谱曲线方程：

??1?10T(??1)，........0.?..T?0.1max?? ????max，.........................0..1..?T?Tg ( 5-7)

?T?(0)0.9*?，..........T...g.?T.max??T

图5-2 大坝设计反应谱

本次重力坝抗震性能分析中，?max取值为2，Tg取值为0.2，特征周期T0取值为0.2S。 3)时程分析法

时程分析方法是将地震动记录或人工波作用在结构上.直接对结构运动方程进行积分，求得结构任意时刻地震反应的分析方法，所以动态时程分析方法也称为直接积分法。 ? 本次大坝抗震性能实例分析采用反应谱分析方法。

5.2.5 后处理

后处理的目的是以图和表的形式描述计算结果。对于大坝抗震性能分析中，进入后处理器

后，查看大坝变形图和节点的位移和应力。通过研究大坝的变形、位移和应力情况，来综合判

ANSYS10.0 土木工程应用实例分析 断大坝的抗震性能及安全性能。

命令：／POST1

GUI： Main Menu> General Postproc

? 首先查看大坝静力分析求解结果，再查看大坝动力分析求解结果。

5.3 ANSYS重力坝抗震性能实例分析

5.3.1 实例介绍

图5-3 重力坝断面结构

实例选取应用非常广泛的重力坝，断面结构如图5-3所示。坝高120米，坝底宽为76米，坝顶为10米，上游坝面坡度和下游坝面坡度如图中所示。

因为重力坝结构比较简单，垂直于长度方向的断面结构受力分布情况也基本相同，并且大坝的纵向长度远大于其横断面，因此大坝抗震性能分析选用单位断面进行平面应变分析是可行

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

大坝抗震性能分析的计算条件如下：

1） 假设大坝的基础是嵌入到基岩中，地基是刚性的。

2） 大坝采用的材料参数为：弹性模量E=35GPa，泊松比?=0.2，容重?=25KN/m3。 3） 计算分析大坝水位为120米。 4） 水的质量密度1000kg/m3。

5） 大坝设防地震烈度为8，水平方向地震加速度值为0.2g。

5.3.2 GUI操作方法

5.3.2.1 创建物理环境

1) 在【开始】菜单中依次选取【所有程序】/【ANSYS10.0】/【ANSYS Product Launcher】，得到“10.0ANSYS Product Launcher”对话框。 2）选中【File Management】，在“Working Directory”栏输入工作目录“D:\\ansys\\example5-1”，在“Job Name”栏输入文件名“Dam”。

3）单击“RUN”按钮，进入ANSYS10.0的GUI操作界面。

4）过滤图形界面：Main Menu> Preferences，弹出“Preferences for GUI Filtering”对话框，选中“Structural”来对后面的分析进行菜单及相应的图形界面过滤。

5）定义工作标题：Utility Menu> File> Change Title，在弹出的对话框中输入“Dam seismic Analysis”，单击“OK”，如图5-4。

图5-4 定义工作标题

6）定义单元类型：

a.定义PLANE42单元：Main Menu> Preprocessor> Element Type> Add/Edit/Delete，弹出一个单元类型对话框，单击“Add”按钮。弹出如图5-5所示对话框。在该对话框左面滚动栏中选择“Solid”，在右边的滚动栏中选择“Quad 4node 42”，单击“Apply”，就定义了“PLANE42”单元。

ANSYS10.0 土木工程应用实例分析

图5-5 定义PLANE42单元对话框

b.设定PLANE42单元选项：Main Menu> Preprocessor> Element Type> Add/Edit/Delete，弹出一个单元类型对话框，选中“Type 2 PLANE42”，单击“Options”按钮，弹出一个“PLANE42 element Type options”对话框，如图5-6所示。在“Element behavior K3”栏后面的下拉菜单中选取“Plane strain”，其它栏后面的下拉菜单采用ANSYS默认设置就可以，单击“OK”按钮。

图5-6 PLANE42单元库类型选项对话框

? 通过设置PLANE42单元选项“K3”为“Plane strain”来设定本实例分析采取平面应变模

型进行分析。因为大坝是纵向很长的实体，故计算模型可以简化为平面应变问题。

7）定义材料属性

执行Main Menu> Preprocessor> Material Props> Material Models，弹出“Define Material Model Behavior”对话框，如图5-7所示。

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

图 5-7 定义材料本构模型对话框

在图5-7中右边栏中连续双击“Structural> Linear> Elastic>Isotropic”后，又弹出如图5-8所示“Linear Isotropic Properties for Material Number 1”对话框，在该对话框中“EX”后面的输入栏输入“3.5E10”，在“PRXY” 后面的输入栏输入“0.2”，单击“OK”。再在选中“Density”并双击，弹出如图5-9所示“Density for Material Number 1”对话框，在“DENS”后面的栏中输入边坡土体材料的密度“2500”，单击“OK”按钮。

图 5-8 线弹性材料模型对话框 图 5-9材料密度输入对话框 5.3.2.2 建立模型和划分网格

1）创建大坝线模型

a.输入关键点：Main Menu> Preprocessor> Modeling> Create>Keypoints>In Active CS，弹出“Creae Keypoints in Active Cooedinate System”对话框，如图5-10所示。在“NPT keypoint number”

ANSYS10.0 土木工程应用实例分析 栏后面输入“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”栏后面输入“（76，0，0）、（15.6，104.1，0）、（15.6，120，0）、（5.6，120，0），最后单击“OK”按钮。

图5-10 在当前坐标系创建关键点对话框

b.创建坝体线模型：Main Menu> Preprocessor> Modeling> Create>Lines>Straight line，弹出“Creae straight lines”对话框，用鼠标依次点击关键点1、2，单击“Apply”按钮，这样就创建了直接L1，同样分别连接关键点“2、3”，“3、4”，“4、5”，“5、1”，最后单击“OK”按钮，就得到坝体线模型，如图5-11所示。

图5-11 坝体线模型

2）创建坝体面模型

a.打开关键点编号显示：Utility Menu> PlotCtrls> Numbering，弹出“Plot Numbering Controls”对话框，如图5-12所示。选中“Keypoint Numbers” 选项，后面的文字由“off”变为“on”,单击“OK”关闭窗口。

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

图5-12 打开关键点编号对话框

b.创建坝体面模型：Main Menu>Preprocessor>Modeling>Create>Areas>Arbitrary>Through KPs，弹出一个“Create Area by Keypoints”对话框，在图形中选取关键点1、2、3、4和5，点击“Apply ”按钮，就得到坝体模型的面模型，如图5-13所示。

图5-13 坝体面模型

3）划分坝体单元网格

a.设置网格份数：Main Menu>Preprocessor>Meshing>Size Cntrls>ManualSize>Layers>Picked

ANSYS10.0 土木工程应用实例分析 Lines，弹出一个“Set Layer Controls”对话框，如图5-14所示，用鼠标选取线L1，单击“OK”按钮。弹出一个“Area Layer Mesh Control on Picked lines”对话框，如图5-15所示，在“No of line division”栏后面输入“20”，单击“OK”按钮。

图5-14 选取线对话框 图5-15 设置网格分数对话框

相同方法设置线L2分割份数为32；设置线L3、L4和L5线的分割份数分别为6、4、40。

b.划分单元网格：Main Menu>Preprocessor>Meshing>Mesh>Area>Free，弹出一个拾取面积对话框，拾取图形中面，单击拾取框上的“OK”按钮，得到坝体模型单元网格，如图5-16所示。

4)保存坝体单元网格

Utility Menu> File> Save as，弹出一个“Save Database”对话框，在“Save Database to”下面输入栏中输入文件名“dam-grid.db”，单击“OK”。

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

图5-16 坝体单元网格

图5-17 给坝体底部施加位移约束对话框

5.3.2.3 施加约束和荷载

1）给坝体模型底部施加位移约束

执行Main Menu>Solution>Define Loads>Apply>Structural>Displacement>on Nodes，弹出在节点上施加位移约束对话框，用鼠标选取隧道模型两侧边界上所有节点，单击“OK”按钮。弹出“Apply U,ROT on Nodes”对话框，如图5-17所示，在“DOFS to be constrained”栏后面中选取“ALL DOF”，在“Apply as”栏后面的下拉菜单中选取“Constant value”，在“Displacement

ANSYS10.0 土木工程应用实例分析 value”栏后面输入“0”值，然后单击“OK”按钮。

2）施加重力加速度：Main Menu>Solution>Define Loads>Apply>Structural>Inertia>Gravity，

弹出“Apply(Gravitational)Acceleration”对话框，如图5-18所示。在“Global Cartesian Y-comp”栏后面输入重力加速度值“9.8”就可以，单击“OK”按钮，就完成了重力加速度的施加。

图5-18 施加重力加速度对话框

3)施加水压力载荷：Main Menu>Solution>Define Loads>Apply>Structural>Structure>on Lines，

弹出一个对话框，用鼠标选中线L5，单击“OK”。弹出“Apply PRES on lines”对话框，如图5-19所示。分别输入数据“0”和“1101370”，单击“OK”按钮，就完成了水压力载荷的施加。

图5-19 施加水压力载荷对话框

? 本次加的荷载是水深为120米时作用在坝上的水压力，迎水面波度是87度。 5.3.2.4 求解

1.静力分析求解 1)求解设置

第5章 ANSYS水利工程应用实例分析 a.指定求解类型：Main Menu>Solution>Analysis Type>New Analysis，弹出一个如图5-20所示对话框，在“Type of analysis”栏后面选中“Static”，单击“OK”按钮。

图5-20指定求解类型对话框

b.设置载荷步： Main Menu>Preprocessor>Loads>Analysis Type>Sol'n 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 Menu>Solution>Analysis Type>Sol'n Controls，弹出一个“Solution Controls”对话框，用鼠标单击“Nonlinear”选项，如图5-22所示，在“Line search”栏后面下拉菜单选中“ON”，单击“OK”。

2）静力求解

a.求解：Main Menu>Solution>Solve>Current LS，弹出一个求解选项信息和一个当前求解载荷步对话框，检查信息无错误后，单击“OK”，开始求解运算，直到出现一个“Solution is done”

的提示栏，表示求解结束。

ANSYS10.0 土木工程应用实例分析 b. 保存求解结果；Utility Menu> File> Save as，弹出一个“Save Database”对话框，在“Save Database to”下面输入栏中输入文件名“Dam-static.db”，单击“OK”。

图5-22设置线性搜索对话框

图5-23 设置模态分析选项对话框

2.抗震性能分析求解 1）模态分析求解

a.设置分析类型：Main Menu>Solution>Analysis Type>New Analysis，弹出一个如图5-20所示对话框，在“Type of analysis”栏后面选中“Modal”，单击“OK”按钮。

第5章 ANSYS水利工程应用实例分析 b.设置模态分析选项：Main Menu>Solution>Analysis Type>Analysis Options，弹出一个如图5-23所示对话框，在“Mode extraction method”栏后面选中“Sunspace”，在“No. of modes to be extract”栏后面输入“18”，在“Expand mode shapes”后面小方框用鼠标选中，单击“OK”按钮，又弹出一个“Sunspace Modal Analysis”对话框，如图5-24所示，按图中设置后，单击“OK”按钮。

图5-24 模态分析求解时子空间设置对话框

图5-25 模态求解选项信息

c.模态分析求解：Main Menu>Solution>Solve>Current LS，弹出一个模态求解选项信息（如图5-25所示）和一个当前求解载荷步对话框，检查信息无错误后，单击“OK”，开始求解运算，直到出现一个“Solution is done”的提示栏，表示求解结束。

d. 保存求解结果；Utility Menu> File> Save as，弹出一个“Save Database”对话框，在“Save

ANSYS10.0 土木工程应用实例分析 Database to”下面输入栏中输入文件名“Dam-modal.db”，单击“OK”。

图5-26 模态分析各阶频率

e.调出模态分析各阶频率：Main Menu> General Postproc >Read Summary，弹出如图5-26所示对话框。

? 动力求解和静力求解的模型相同，约束条件也相同。 ? 谱分析时，ANSYS忽略材料非线性。

? 调出模态分析各阶频率是为后面求解反应谱值。

2）反应谱分析求解

a.求出反应谱值：由图5-26中前18阶频率值f，可以算出对应的周期T，再根据大坝反应谱曲线方程（式5-7），可以计算出前10阶的反应谱值，见表5-1。

表5-1 大坝动力计算前10阶振动频率及反应谱值 振型 1 2 3 4 5 6 7 振动频率1/s 3.5138 8.0730 11.118 14.30 21.776 24.980 30.734 振动周期T 0.2846 0.1239 0.0899 0.0699 0.0459 0.0400 0.0325 反应谱值 1.456 2.0 1.899 1.699 1.459 1.400 1.325 第5章 ANSYS水利工程应用实例分析 8 9 10 11 12 13 14 15 16 17 34.119 36.263 39.043 40.852 43.237 47.907 49.639 52.449 55.521 56.661 0.0293 0.0276 0.0256 0.0245 0.0231 0.0209 0.0201 0.0191 0.0180 0.0176 0.0171 1.293 1.276 1.256 1.245 1.231 1.209 1.201 1.191 1.180 1.176 1.171 18 58.613 b.设置反应谱分析求解选项

? 设置分析类型：Main Menu>Solution>Analysis Type>New Analysis，弹出一个如图5-27所示对话框，在“Type of analysis”栏后面选中“Spectrum”，单击“OK”按钮。

图5-27 定义反应谱分析对话框

? 设置反应谱分析选项：Main Menu>Solution>Analysis Type>Analysis Options，弹出一个“Spectrum Analysis”对话框，如图5-28所示，在“Type of spectrum”后面栏中选取“Single-pt resp”，在“No. of modes for solu”后面输入“10”，在“Calculate elem stress?”后面选中“Yes”，单击“OK”按钮。

? 设置反应谱单点分析选项：Main Menu>Solution>Load Step Opts>Spectrum>SinglePt> Settings，弹出一个“Setting for single-point Response Spectrum”对话框，如图5-29所示。在“Type of response spectrum”栏后面下拉菜单选中“Seismic accel”，在“SEDX，SEDY，SEDZ”栏后面依次输入“0、1、0”，单击“OK”按，钮。

ANSYS10.0 土木工程应用实例分析

图5-28 设置反应谱分析对话框

图5-29 设置单点反应谱分析对话框

?定义反应谱分析频率表：Main Menu>Solution>Load Step Opts>Spectrum>SinglePt>Freq Table，弹出一个“Frequency Table”对话框，如图5-30所示。根据表5-1依次输入大坝的前18阶振动频率。

?定义反应谱值：Main Menu>Solution>Load Step Opts>Spectrum>SinglePt>Spectr Values，弹出一个“Spectrum Values Damping Ration”对话框，单击“OK”按钮。弹出一个“Spectrum Values”对话框，如图5-30所示。根据表5-1依次输入大坝的前18阶反应谱值。 ? 注意，输入的振动频率必须按升序排列。 ? FREQ1必须大于零。

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

图5-30 定义频率表对话框

图5-31 定义反应谱值对话框

c.反应谱分析求解：Main Menu>Solution>Solve>Current LS，弹出一个模态求解选项信息和一个当前求解载荷步对话框，检查信息无错误后，单击“OK”，开始求解运算，直到出现一个“Solution is done”的提示栏，表示求解结束。

d. 保存求解结果；Utility Menu> File> Save as，弹出一个“Save Database”对话框，在“Save Database to”下面输入栏中输入文件名“Dam-spectrum.db”，单击“OK”。

3) 模态扩展分析求解

a.设置分析类型：Main Menu>Solution>Analysis Type>New Analysis，弹出一个“New Analysis”对话框，在“Type of analysis”栏后面选中“Modal”，单击“OK”按钮。

ANSYS10.0 土木工程应用实例分析 b.设置模态扩展分析求解选项

?定义模态扩展分析：Main Menu>Solution>Analysis Type>ExpansionPass，弹出一个“Expansion Pass”对话框，如图5-32所示，选中“Expansion Pass” 选项，后面的文字由“off”变为“on”,单击“OK”关闭窗口。

图5-32 模态扩展对话框

?设置模态扩展分析：Main Menu>Preprocessor>Loads>Load Step Opts>ExpansionPass>Single Expand>Expand Modes，弹出一个“Expand Modes”对话框，如图5-33所示，在“No. of modes to expand”栏后面输入“18”，其它如图中设置，单击“OK”按钮。

图5-33 模态扩展分析设置对话框

c.模态扩展分析求解：Main Menu>Solution>Solve>Current LS，弹出一个模态求解选项信息和一个当前求解载荷步对话框，检查信息无错误后，单击“OK”，开始求解运算，直到出现一个“Solution is done”的提示栏，表示求解结束。

d. 保存求解结果；Utility Menu> File> Save as，弹出一个“Save Database”对话框，在“Save Database to”下面输入栏中输入文件名“Dam-expand.db”，单击“OK”。

4)合并模态分析求解

a.设置分析类型：Main Menu>Solution>Analysis Type>New Analysis，弹出一个“New Analysis”对话框，在“Type of analysis”栏后面选中“spectrum”，单击“OK”按钮。

b.按平方和方根法进行组合：Main Menu>Solution>Load Step Opts> Spectrum> SinglePt> Mode Combine，弹出一个“Mode Combination Methods”对话框，如图5-34所示，在“Mode Combination Method”栏后面下拉菜单中选取“SRSS”，在“Signification threshold”栏后面输入

第5章 ANSYS水利工程应用实例分析 “0.1”，在“Type of output”栏后面下拉菜单中选取“Displacement”，单击“OK”按钮。

图5-34 合并模态求解设置对话框

c.合并模态分析求解：Main Menu>Solution>Solve>Current LS，弹出一个模态求解选项信息和一个当前求解载荷步对话框，检查信息无错误后，单击“OK”，开始求解运算，直到出现一个“Solution is done”的提示栏，表示求解结束。

d. 保存求解结果；Utility Menu> File> Save as，弹出一个“Save Database”对话框，在“Save Database to”下面输入栏中输入文件名“Dam- combination.db”，单击“OK”。 5.3.2.5 后处理

对坝体进行了静力分析和动力分析，通过对重力坝静力有限元分析，可以知道坝体在静力载荷作用下的位移场和应力场，从而可以了解坝体的安全性能。通过对重力坝抗震性能有限元分析，可以了解大坝的地震载荷作用下的动力响应特性，从而可以评价大坝在地震载荷作用下的安全性能。

1.静力分析求解结果

a.读入静力分析求解结果数据：Utility Menu> resume…，弹出一个“Resume Database”对话框，选中刚才保存的文件“Dam-static.db”，单击“OK”按钮。

图5-35 绘制变形图对话框

b.绘制坝体变形图：Main Menu>General Postproc>Plot Results>Deformed Shape，弹出一个“Plot Deform Shape”对话框，如图5-35所示。选中“Def +undeformed”，单击“OK”按钮，得到坝体变形图，如图5-36所示。

ANSYS10.0 土木工程应用实例分析

图5-36 坝体变形图

5-37 节点解云图绘制对话框

c. 显示坝体x方向位移云图：Main Menu>General Postproc>Plot Results>Contour Plot>Nodal Solu，弹出一个“Contour Nodal Solution Data”对话框，如图5-37所示，用鼠标依次点击“Nodal Solution/DOF Solution/ X-Compoment of displacement”，再单击“OK”按钮，就得到坝体X方向位移云图，如图5-38所示。此时，坝体水平方向最大位移为12.111mm，位置发生在坝顶。

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

图5-38 坝体x方向位移云图

d.显示坝体y方向位移云图：Main Menu>General Postproc>Plot Results>Contour Plot>Nodal Solu，弹出一个“Contour Nodal Solution Data”对话框，如图5-37所示，用鼠标依次点击“Nodal Solution/Stress/ Y-Compoment of displacement”，再单击“OK”按钮，就得到坝体Y方向位移云图，如图5-39所示。此时，坝体垂直方向最大位移为0.44mm，发生在坝腹处。

图5-39 坝体y方向位移云图

e.显示坝体x方向应力云图：Main Menu>General Postproc>Plot Results>Contour Plot>Nodal Solu，弹出一个“Contour Nodal Solution Data”对话框，如图5-37所示，用鼠标依次点击“Nodal Solution/stress/ x-Compoment of stress”，再单击“OK”按钮，就得到坝体x方向应力云图，如图5-40所示。

f.显示坝体y方向应力云图：Main Menu>General Postproc>Plot Results>Contour Plot>Nodal Solu，弹出一个“Contour Nodal Solution Data”对话框，如图5-37所示，用鼠标依次点击“Nodal Solution/stress/ y-Compoment of stress”，再单击“OK”按钮，就得到坝体x方向应力云图，如图5-41所示。

ANSYS10.0 土木工程应用实例分析

图5-40 坝体x方向应力云图

图5-41 坝体y方向应力云图

图5-42 第一主应力云图

g.显示坝体第一主应力云图：Main Menu>General Postproc>Plot Results>Contour Plot>Nodal

第5章 ANSYS水利工程应用实例分析 Solu，弹出一个“Contour Nodal Solution Data”对话框，如图5-37所示，用鼠标依次点击“Nodal Solution/stress/ 1st principal stress”，再单击“OK”按钮，就得到坝体第一主应力云图，如图5-42所示。

2. 抗震性能分析求解结果

1）绘制坝体振型图

a.读入模态分析分析求解结果数据：Utility Menu> resume…，弹出一个“Resume Database”对话框，选中刚才保存的文件“Dam-modal.db”，单击“OK”按钮。

b.绘制坝体第1阶震型

? 读入第1阶数据：Main Menu>General Postproc>Read Results>First set.

? 绘制坝体第1阶震型：Main Menu>General Postproc>Plot Results>Deformed Shape，弹出一个“Plot Deform Shape”对话框，选中“Def +undeformed”，单击“OK”按钮，得到坝体第1阶震型，如图5-43所示。

c. 绘制坝体第2阶震型

? 读入第2阶数据：Main Menu>General Postproc>Read Results>Next set.

? 绘制坝体第2阶震型：Main Menu>General Postproc>Plot Results>Deformed Shape，弹出一个“Plot Deform Shape”对话框，选中“Def +undeformed”，单击“OK”按钮，得到坝体第2阶震型，如图5-44所示。

图5-43 坝体第1阶振型 图5-44 坝体第2阶振型

依次重复执行c步骤，就可以依次画出坝体第3到第18阶振型，如图5-45~图5-60所示。 d.调出模态分析结果的固有频率和振型：Main Menu> General Postproc >Read Summary，弹出如图5-26所示对话框。

ANSYS10.0 土木工程应用实例分析 图5-45 坝体第3阶振型 图5-47 坝体第5阶振型 图5-49 坝体第7阶振型

图5-46 坝体第4阶振型

图5-48 坝体第6阶振型

图5-50 坝体第8阶振型

第5章 ANSYS水利工程应用实例分析 图5-51 坝体第9阶振型 图5-53 坝体第11阶振型 图5-55 坝体第13阶振型

图5-52 坝体第10阶振型

图5-54 坝体第12阶振型

图5-56 坝体第14阶振型

ANSYS10.0 土木工程应用实例分析

图5-57 坝体第15阶振型 图5-58 坝体第16阶振型

图5-59 坝体第17阶振型 图5-60 坝体第18阶振型

2)绘制坝体位移云图

a.读入抗震性能分析最终求解结果数据：Utility Menu> resume…，弹出一个“Resume Database”对话框，选中刚才保存的文件“Dam- comination.db”，单击“OK”按钮。 b.绘制合并模态求解后坝体位移云图：

因为响应谱分析是在频域内进行的，对于结构动力特性依赖于频率而变化，因此在模态分析后要进行模态合并求解，才能得到坝体结构真实的总体效应。合并模态求解后，得到坝体在各阶频率的真实位移云图。

? 读入第1阶数据：Main Menu>General Postproc>Read Results>First set.

? 绘制坝体第1阶位移云图：Main Menu>General Postproc>Plot Results>Contour Plot>Nodal Solu，弹出一个“Contour Nodal Solution Data”对话框，用鼠标依次点击“Nodal Solution/DOF Solution/ x-Compoment of displacement得到坝体第1阶x方向位移云图，如图5-61所示.同法选择y-Compoment of displacement, 得到坝体第1阶y方向位移云图，如图5-62所示。 ? 读入下一阶数据：Main Menu>General Postproc>Read Results>Next set.

? 绘制坝体下一阶位移云图：Main Menu>General Postproc>Plot Results>Contour Plot>Nodal Solu，弹出一个“Contour Nodal Solution Data”对话框，用鼠标依次点击“Nodal Solution/DOF

第5章 ANSYS水利工程应用实例分析 Solution/ x-Compoment of displacement得到坝体下一阶x方向位移云图，如图5-63所示.同法选择y-Compoment of displacement, 得到坝体下一阶阶y方向位移云图，如图5-64所示。重复读入下一个数据，可绘制出合并模态求解后坝体真实位移云图，如图5-65~图5-72所示。

图5-61 坝体第1阶x方向位移云图 图5-63 坝体第2阶x方向位移云图 图5-65 坝体第3阶x方向位移云图

图5-62 坝体第1阶y方向位移云图

图5-64 坝体第2阶y方向位移云图

图5-66 坝体第3阶y方向位移云图

ANSYS10.0 土木工程应用实例分析 图5-67 坝体第4阶x方向位移云图 图5-69 坝体第5阶x方向位移云图 图5-71 坝体第7阶x方向位移云图

图5-68 坝体第4阶y方向位移云图

图5-70 坝体第5阶y方向位移云图

图5-72 坝体第7阶y方向位移云图

第5章 ANSYS水利工程应用实例分析 3）绘制坝体应力/应变云图

? 读入第1阶数据：Main Menu>General Postproc>Read Results>First set.

? 绘制坝体第1阶应力：Main Menu>General Postproc>Plot Results>Contour Plot>Nodal Solu，弹出一个“Contour Nodal Solution Data”对话框，用鼠标依次点击“Nodal Solution/ stress/ x-Compoment of stress”得到坝体第1阶x方向应力云图，如图5-73所示。同法选择y-Compoment of stress/1st principal stress/, 得到坝体第1阶y方向应力/第1主应力云图，如图5-74和图5-75所示。

? 绘制坝体第1阶主应变云图：Main Menu>General Postproc>Plot Results>Contour Plot>Nodal Solu，用鼠标依次点击“Nodal Solution/ strain/1st principal strain”得到坝体第1阶第1主应变云图，如图5-76所示。

? 读入下一阶数据：Main Menu>General Postproc>Read Results>Next set.

? 绘制坝体下一阶应力/应变云图：与绘制第1阶一样，就可以绘制出下一阶坝体应力/应变云图。重复读入下一个数据，可绘制出合并模态求解后坝体应力云图，如图5-77~图5-7

图5-73 坝体第1阶x方向应力云图 图5-74 坝体第1阶y方向应力云图

图5-75 坝体第1阶第1主应力云图 图5-76 坝体第1阶第1主应变云图

ANSYS10.0 土木工程应用实例分析 图5-77 坝体第2阶x方向应力云图 图5-79 坝体第2阶第1主应力云图 图5-81 坝体第3阶第1主应力云图

图5-78 坝体第2阶y方向应力云图

图5-80 坝体第2阶第1主应变云图

图5-82 坝体第3阶第1主应变云

第5章 ANSYS水利工程应用实例分析 图5-84 坝体第3阶x方向应力云图 图5-86 坝体第4阶x方向应力云图 图5-88 坝体第4阶第1主应力云图

图5-85 坝体第3阶y方向应力云图

图5-87 坝体第4阶y方向应力云图

图5-89 坝体第4阶第1主应变云

ANSYS10.0 土木工程应用实例分析 图5-90 坝体第5阶第1主应力云图 图5-92 坝体第5阶x方向应力云图 图5-94 坝体第7阶x方向应力云图

图5-91 坝体第5阶第1主应变云

图5-93 坝体第5阶y方向应力云图

图5-95 坝体第7阶y方向应力云图

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

图5-96 坝体第7阶第1主应力云图 图5-97 坝体第7阶第1主应变云图 进行模态求解时并没有同时进行模态扩展，而是在获得谱解后又单独进行了扩展模态求解。这样，有明显意义的模态为1, 2, 3, 4, 5和7阶模态。 大坝水平方向最大位移发生在坝顶。

5.3.2 命令流实现

(1)创建物理环境

/COM, Structural ! 指定结构分析 /TITLE, Dam seismic Analysis ! 定义工作标题 /FILNAM,dam,1 ! 定义工作文件名 (2)建立模型 !进入前处理器 /PREP7

!定义单元类型

ET,1,PLANE42 !定义坝体单元 KEYOPT,1,3,2 !定义平面应变 !定义坝体材料属性

MP,EX,1,3.5E10 !弹性模量 MP,PRXY,1,0.2 !泊松比 MP,DENS,1,2500 !密度 !创建坝体线模型

K,1,,,, !创建关键点 K,2,76,,, K,3,15.6,104.1,,

K,4,15.6,120,,

ANSYS10.0 土木工程应用实例分析 K,5,5.6,120,,

l,1,2 !连接关键点1、2生成直线 l,2,3 l,3,4 l,4,5

l,5,1

!创建边坡面模型

/PNUM, kp, on !划分网格生成有限元模型

a,1,2,3,4,5 ! 划分坝体单元网格 !设置网格份数

lsel,s,,,1 lesize,all,,,20 lsel,s,,,2 lesize,all,,,32 lsel,s,,,3 lesize,all,,,6 lsel,s,,,4 lesize,all,,,4 lsel,s,,,5 lesize,all,,,40 amesh,all save,dam-grid.db (3)施加约束和荷载

!坝体底部施加约束

nsel,s,loc,y,0 d,all,uy d,all,ux

acel,,9.8 sfl,5,pres,0,1101370 (4)求解 1)静力分析求解

/solu

antype,static !打开关键点号开关 !由关键点生成面 !选择线L1

!把所选择线分为20段 !选择线L2

!把所选择线分为32段 !选择线L1

!把所选择线分为6段 !选择线L2

!把所选择线分为4段 !选择线L2

!把所选择线分为40段 !划分坝体单元 !保存网格模型 !选择Y=0线上所有节点

!对所选择节点约束X、Y方向位移 !施加重力加速度 !施加水压力荷载 !设定为静力求解

第5章 ANSYS水利工程应用实例分析 nsubst,100 !设定最大子步数为100 lnsrch,on !打开线性搜索 allsel

solve !进行求解

save,dam-static,db !保存求解结果数据 finish

2)抗震性能分析求解 !模态分析

/solu

antype,modal !设定为模态求解

modopt,subsp,18 !采用子空间迭代法提取18阶模态 solve !进行模态求解

save,dam-modal,db !保存求解结果数据 set,list !调出模态分析各阶频率 finish

!反应谱分析

/solu

antype,spectrum !设定为反应谱求解

spopt,sprs,18,yes !设置单点反应谱分析选项 svtyp,2 !设定为地震加速度谱分析 sed,1,0,0 !定义反应谱激磁方向是x方向 freq,3.5138,8.073,11.118,14.3,21.776,24.98,30.734,34.119,36.263,39.043,40.852,43.237,47.907,49.679,52.521,56.661,58.613 ！输入振动频率

sv,1.456,2,1.899,1.699,1.459,1.4,1.325,1.293,1.276,1.256,1.245,1.231,1.209,1.201,1.191,1.18,1.176,1.171 ！输入对应反应谱值

solve !进行求解 save,Dam-spectrum,db !保存求解结果 !模态扩展求解

/solu

antype,modal !设定为模态求解 expass,on !定义为模态扩展分析 mxpand,18,,,yes,0.001 !设置模态扩展求解选项 solve !进行求解

save,Dam-expand,db !保存求解结果数据 finish

!合并模态求解

ANSYS10.0 土木工程应用实例分析 /solu

antype,spectrum !设定为反应谱求解

srss,0.15,disp !设置位移响应谱平方和平方根输出

solve

save,Dam-combination,db !保存求解结果数据 finish

（5）后处理 1)静力分析求解结果

/post1

Resume,'dam-static','db' pldisp,1 plnsol,u,x plnsol,u,y plnsol,s,x plnsol,s,y plnsol,s,1 finish

2)抗震性能分析求解结果 !绘制坝体振型图

/post1

Resume,'dam-modal','db' set,first pldisp,1 set,next pldisp,1 set,next pldisp,1 !......

set,next pldisp,1 finish

!绘制坝体位移/应力/应变云图

/post

Resume,'dam-conbination','db' set,first plnsol,u,x !读入静力分析求解结果数据 !绘制坝体变形图 !绘制坝体x方向位移云图 !绘制坝体y方向位移云图 !绘制坝体x方向应力云图 !绘制坝体y方向应力云图 !绘制坝体第一主应力云图 读入模态分析求解结果数据 读入第一阶数据 !绘制坝体第一阶振型图 !读入第2阶数据

!绘制坝体第2阶振型图 !读入第3阶数据 !绘制坝体第3阶振型图 !读入第18阶数据

!绘制坝体第18阶振型图 !读入合并模态分析求解结果数据 !读入第一阶数据

!绘制坝体第1阶x方向位移云图

! ! 第5章 ANSYS水利工程应用实例分析 plnsol,u,y !绘制坝体第1阶y方向位移云图 plnsol,s,x !绘制坝体第1阶x方向应力云图 plnsol,s,y !绘制坝体第1阶y方向应力云图 plnsol,s,1 !绘制坝体第1阶第一主应力云图 set,next !读入第2阶数据

plnsol,u,x plnsol,u,y plnsol,s,x plnsol,s,y plnsol,s,1 set,next plnsol,u,x plnsol,u,y plnsol,s,x plnsol,s,y plnsol,s,1 set,next

plnsol,u,x plnsol,u,y plnsol,s,x plnsol,s,y plnsol,s,1 set,next plnsol,u,x plnsol,u,y plnsol,s,x plnsol,s,y plnsol,s,1 set,next plnsol,u,x plnsol,u,y plnsol,s,x plnsol,s,y plnsol,s,1 !绘制坝体第2阶x方向位移云图 !绘制坝体第2阶y方向位移云图 !绘制坝体第2阶x方向应力云图 !绘制坝体第2阶y方向应力云图 !绘制坝体第2阶第一主应力云图 !读入第3阶数据

!绘制坝体第3阶x方向位移云图 !绘制坝体第3阶y方向位移云图 !绘制坝体第3阶x方向应力云图 !绘制坝体第3阶y方向应力云图 !绘制坝体第3阶第一主应力云图 !绘制坝体第4阶x方向位移云图 !绘制坝体第4阶y方向位移云图 !绘制坝体第4阶x方向应力云图 !绘制坝体第4阶y方向应力云图 !绘制坝体第4阶第一主应力云图 !读入第5阶数据

!绘制坝体第5阶x方向位移云图 !绘制坝体第5阶y方向位移云图 !绘制坝体第5阶x方向应力云图 !绘制坝体第5阶y方向应力云图 !绘制坝体第5阶第一主应力云图 !读入第7阶数据

!绘制坝体第7阶x方向位移云图 !绘制坝体第7阶y方向位移云图 !绘制坝体第7阶x方向应力云图 !绘制坝体第7阶y方向应力云图 !绘制坝体第7阶第一主应力云图

本文来源：https://www.bwwdw.com/article/sn5w.html