第4章-ANSYS边坡工程计划应用实例分析.doc

《第4章-ANSYS边坡工程计划应用实例分析.doc》由会员分享，可在线阅读，更多相关《第4章-ANSYS边坡工程计划应用实例分析.doc（65页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、-_第第 4 4 章章 ANSYSANSYS 边坡工程应用实例分析边坡工程应用实例分析本章重点本章重点边坡工程概述 ANSYS 边坡稳定性分析步骤ANSYS 边坡稳定性实例分析本章典型效果图本章典型效果图第 4 章 ANSYS 边坡工程应用实例分析4.1 边坡工程概述4.1.1 边坡工程边坡指地壳表部一切具有侧向临空面的地质体，是坡面、坡顶及其下部一定深度坡体的 总称。坡面与坡顶面下部至坡脚高程的岩体称为坡体。 倾斜的地面称为斜坡，铁路、公路建筑施工中，所形成的路堤斜坡称为路堤边坡；开挖 路堑所形成的斜坡称为路堑边坡；水利、市政或露天煤矿等工程开挖施工所形成的斜坡也称 为边坡；这些对应工程就称

2、为边坡工程 对边坡工程进行地质分类时，考虑了下述各点。首先，按其物质组成，即按组成边坡的 地层和岩性，可以分为岩质边坡和土质边坡(后者包括黄土边坡、砂土边坡、土石混合边坡)。 地层和岩性是决定边坡工程地质特征的基本因素之一，也是研究区域性边坡稳定问题的主要 依据.其次，再按边坡的结构状况进行分类。因为在岩性相同的条件下，坡体结构是决定边 坡稳定状况的主要因素，它直接关系到边坡稳定性的评价和处理方法。最后，如果边坡已经 变形，再按其主要变形形式进行划分。即边坡类属的称谓顺序是:岩性 结构变形。 边坡工程对国民经济建设有重要的影响：在铁路、公路与水利建设中，边坡修建是不可 避免的，边坡的稳定性严重

3、影响到铁路、公路与水利工程的施工安全、运营安全以及建设成 本。在路堤施工中，在路堤高度一定条件下，坡角越大，路基所占面积就越小，反之越大。 在山区，坡角越大，则路堤所需填方量越少。因此，很有必要对边坡稳定性进行分析，4.1.2 边坡变形破坏基本原理4.1.2.1 应力分布状态边坡从其形成开始，就处于各种应力作用(自重应力、构造应力、热应力等)之下。在边 坡的发展变化过程中，由于边坡形态和结构的不断改变以及自然和人为营力的作用，边坡的 应力状态也随之调整改变。根据资料及有限元法计算，应力主要发生以下变化： （1)岩体中的主应力迹线发生明显偏转，边坡坡面附近最大主应力方向和坡而平行，而 最小主应力

4、方向则与坡面近于垂直，并开始出现水平方向的剪应力，其总趋势是由 内向外增多，愈近坡脚愈高，向坡内逐渐恢复到原始应力状态。 (2) 在坡脚逐渐形成明显的应力集中带。边坡愈陡，应力集中愈严重，最大最小主应力ANSYS10.0 土木工程应用实例分析的差值也愈大。此外，在边坡下边分别形成切向应力减弱带和水平应力紧缩带，而 在靠近边坡的表部所测得的应力值均大于按上覆岩体重量计算的数值。 (3) 边坡坡面岩体由于侧向应力近于零，实际上变为两向受力。在较陡边坡的坡面和顶 面，出现拉应力，形成拉应力带.拉应力带的分布位置与边坡的形状和坡面的角度 有关。边坡 应力的调整和拉应力带的出现，是边坡变形破坏最初始的征

5、兆。例如， 由于坡脚应力的集中，常是坡脚出现挤压破碎带的原因；由于坡面及坡顶出现拉应 力带，常是表层岩体松动变形的原因。 4.1.2.2 边坡岩体变形破坏基本形式边坡在复杂的内外地质营力作用下形成，又在各种因素作用下变化发展。所有边坡都在 不断变形过程中，通过变形逐步发展至破坏。其基本变形破坏形式主要有:松弛张裂、滑动、 崩塌、倾倒、蠕动和流动。4.1.3 影响边坡稳定性的因案影响边坡稳定性的主要因素有:(1)边坡材料力学特性参数： 包括弹性模量、泊松比、摩擦角、粘结力、容重、抗剪强度等参数。 (2)边坡的几何尺寸参数：包括边坡高度、坡面角和边坡边界尺寸以及坡面后方坡体的几何形状，即坡体的不

6、连续面与开挖面的坡度及方向之间的几何关系，它将确定坡体的各个部分是否滑动 或塌落。 (3)边坡外部荷载：包括地震力、重力场、渗流场、地质构造地应力等。4.1.4 边坡稳定性的分析方法分析边坡稳定问题，基本上可以分为两种方法：极限平衡方法和数值分析方法。 4.1.4.1 极限平衡方法极限平衡方法的基本思想是:以摩尔一库仑抗剪强度理论为基础，将滑坡体划分成若干 垂直条块，建立作用在垂直条块上的力的平衡方程式，求解安全系数。 这种计算分析方法遵循下列基本假定: (1)遵循库仑定律或由此引伸的准则。第 4 章 ANSYS 边坡工程应用实例分析(2)将滑体作为均质刚性体考虑，认为滑体本身不变形，且可以传

7、递应力。因此只研究 滑动面上的受力大小，不研究滑体及滑床内部的应力状态。 (3)将滑体的边界条件大大简化。如将复杂的滑体型态简化为简单的几何型态；将滑面 简化为圆弧面、平面或折面；一般将立体问题简化为平面问题，取沿滑动方向的代 表性剖面，以表征滑体的基本型态；将均布力简化为集中力，有时还将力的作用点 简化为通过滑体重心。 极限平衡方法包括以下几种方法： (1)瑞典圆弧滑动法 (2)简化逼肖普法 (3)简布普通条分法(4)摩根斯坦-普赖斯法 (5)不平衡推力传递法以上各种方法都是假定土体是理想塑性材料，把土条作为一个刚体，按照极限平衡的原 则进行力的分析，最大的不同之处在于对相邻上条之间的内力作

8、何种假定，也就是如何增加 已知条件使超静定问题变成静定问题。这些假定的物理意义不一样，所能满足的平衡条件也 不相同，计算步骤有繁有简，使用时必须注意他们的适用场合。 极限平衡方法关键是对滑体的休型和滑面的形态进行分析、正确选用滑面的计算参数以 及正确引用滑体的荷载条件等。因为极限平衡方法完全不考虑土体本身的应力-应变关系， 不能真实地反映边坡失稳时的应力场和位移场，因此而受到质疑。 4.1.4.2 数值分析方法数值数值分析方法考虑土体应力应变关系，克服了极限平衡方法完全不考虑土体本身的 应力-应变关系缺点，为边坡稳定分析提供了较为正确和深入的概念。 边坡稳定性数值分析方法主要包含以下几种方法：

9、 (1)有限元法有限单元法是数值模拟方法在边坡稳定评价中应用最早的方法，也是目前最广泛使用的 一种数值方法，可以用来求解弹性、弹塑性、粘弹塑性、粘塑性等问题。目前 用有限元法 求解边坡稳定主要有两种方法。 a.有限元滑面搜索法:将边坡体离散为有限单元格，按照施加的荷载及边界条件进行有 限元计算可得到每个结点的应力张量。然后假定一个滑动面，用有限元数据给出滑动面任一 点的向正应力和剪应力，根据摩尔一库仑准则可得该点的抗滑力，由此即能求得滑动面上每 个结点的下滑力与抗滑力，再对滑动面上下滑力与抗滑力进行积分，就可以求得每一个滑动 面的安全系数。ANSYS10.0 土木工程应用实例分析b.有限元强度

10、折减法:首先选取初始折减系数，将岩土体强度参数进行折减，将折减后 的参数作为输入，进行有限元计算，若程序收敛，则岩土体仍处于稳定状态，然后再增加折 减系数，直到程序恰好不收敛，此时的折减系数即为稳定或安全系数。 (2)自适应有限元法自20世纪70年代开始自适应理论被引入有限元计算，主导思想是减少前处理工作量和实 现网格离散的客观控制。现已基本建立了一般弹性力学、流体动力学、渗流分析等领域的平 面自适应分析系统，能使计算较为快速和准确。 (3)离散单元法离散单元法的突出功能是它在反映岩块之间接触面的滑移、分离与倾翻等大位移的同时， 又能计算岩块内部的变形与应力分布。因此，任何一种岩体材料都可引入

11、到模型中，例如弹 性、粘弹性或断裂等均可考虑，故该法对块状结构、层状破裂或一般破裂结构岩体边坡比较 合适。并且，它利用显式时间差分法(动态差分法)求解动力平衡方程，求解非线性大位移与 动力问题比较容易。 离散元法在模拟过程中考虑了边坡失稳破坏的动态过程，允许岩土体存在滑动、平移、 转动和岩体的断裂及松散等复杂过程，具有宏观上的不连续性和单个岩块休运动的随机性， 可以较真实、动态地模拟边坡在形成和开挖过程中应力、位移和状态的变化，预测边坡的稳 定性，因此在岩质高边坡稳定性的研究中得到广泛的应用。 (4)拉格朗日元法为了克服有限元等方法不能求解大变形问题的缺陷，人们根据有限差分法的原理，提出 了F

12、LAC数值分析方法。该方法较有限元法能更好地考虑岩土体的不连续和大变形特性，求解 速度较快。缺点是计算边界、单元网格的划分带有很大的随意性。 (5)界面元法界面元法是一种基于累积单元变形于界面的界面应力元法模型，建立适用于分析不连续、 非均匀、各向异性和各类非线性问题、场问题，以及能够完全模拟各类锚杆复杂空间布局和 开挖扰动的方法。 4.1.4.3 有限元法用于边坡稳定性分析优点 有限元法考虑了介质的变形特征，真实地反应了边坡的受力状态。它可以模拟连续介质， 也可以模拟不连续介质；能考虑边坡沿软弱结构面的破坏，也能分析边坡的整体稳定破坏。 有限元法可以模拟边坡的圆弧滑动破坏和非圆弧滑动破坏。同

13、时它还能适应各种边界条件和 不规则几何形状，具有很广泛的适用性。 有限元法应用于边坡工程，有其独特的优越性。与一般解析方法相比，有限元法有以下 优点：第 4 章 ANSYS 边坡工程应用实例分析(1) 它考虑了岩体的应力-应变关系，求出每一单元的应力与变形，反映了岩体真实工 作状态。 (2)与极限平衡法相比，不需要进行条间力的简化，岩体自始至终处于平衡状态。 (3)不需要像极限平衡法一样事先假定边坡的滑动面，边坡的变形特性、塑性区形成都 根据实际应力应变状态“自然”形成。 (4)若岩体的初始应力己知，可以模拟有构造应力边坡的受力状态。 (5)不但能像极限平衡法一样模拟边坡的整体破坏，还能模拟边

14、坡的局部破坏，把边坡 的整体破坏和局部破坏纳入统一的体系。 (6)可以模拟边坡的开挖过程，描述和反应岩体中存在的节理裂隙、断层等构造面。鉴于有限元法具有如此多优点，本章借助通用有限元软件ANSYS来实现对边坡稳定性分 析，用具体的边坡工程实例详细介绍应用ANSYS软件分析边坡稳定性问题。4.2 ANSYS 边坡稳定性分析步骤ANSYS 边坡稳定性分析一般分以下几个步骤： 1、创建物理环境 2、建立模型，划分网格，对模型的不同区域赋予特性 3、加边界条件和载荷 4、求解 5、后处理（查看计算结果）4.2.1 创建物理环境在定义边坡稳定性分析问题的物理环境时，进入 ANSYS 前处理器，建立这个边

15、坡稳定 性分析的数学仿真模型。按照以下几个步骤来建立物理环境：1、 设置 GUT 菜单过滤 如果你希望通过 GUI 路径来运行 ANSYS，当 ANSYS 被激活后第一件要做的事情就是 选择菜单路径：Main MenuPreferences，执行上述命令后，弹出一个如图 4-1 所示的对话框 出现后，选择 Structural。这样 ANSYS 会根据你所选择的参数来对 GUI 图形界面进行过滤， 选择 Structural 以便在进行边坡稳定性分析时过滤掉一些不必要的菜单及相应图形界面。 2、 定义分析标题（TITLE） 在进行分析前，可以给你所要进行的分析起一个能够代表所分析内容的标题，比

16、如ANSYS10.0 土木工程应用实例分析“Slope stability Analysis” ，以便能够从标题上与其他相似物理几何模型区别。用下列方法定 义分析标题。 命令：TITLE GUI：Utility MenuFileChange Title 3、 说明单元类型及其选项（KEYOPT 选项） 与 ANSYS 的其他分析一样，也要进行相应的单元选择。ANSYS 软件提供了 100 种以 上的单元类型，可以用来模拟工程中的各种结构和材料，各种不同的单元组合在一起，成为 具体的物理问题的抽象模型。例如，不同材料属性的边坡土体用 PLANE82 单元来模拟。 大多数单元类型都有关键选项（KE

17、YOPTS） ，这些选项用以修正单元特性。例如， PLANE82 有如下 KEYOPTS： KEYOPT(2) 包含或抑制过大位移设置 KEYOPT(3) 平面应力、轴对称、平面应变或考虑厚度的平面应力设置 KEYOPT(9) 用户子程序初始应力设置 设置单元以及其关键选项的方式如下： 命令：ETKEYOPT GUI：Main Menu Preprocessor Element Type Add/Edit/Delete图 4-1 GUI 图形界面过滤 4 定义单位第 4 章 ANSYS 边坡工程应用实例分析结构分析只有时间单位、长度单位和质量单位三个基本单位，则所有输入的数据都应当 是这三个单

18、位组成的表达方式。如标准国际单位制下，时间是秒（s） ，长度是米（m） ，质量 是千克（kg） ，则导出力的单位是 kgm/s2（相当于牛顿 N） ，材料的弹性模量单位是kg/ms2（相当于帕 Pa） 。 命令：/UNITS 5、定义材料属性 大多数单元类型在进行程序分析时都需要指定材料特性，ANSYS 程序可方便地定义各 种材料的特性，如结构材料属性参数、热性能参数、流体性能参数和电磁性能参数等。 ANSYS 程序可定义的材料特性有以下三种： （1）线性或非线性。 （2）各向同性、正交异性或非弹性。 （3）随温度变化或不随温度变化。 因为分析的边坡模型采用理想弹塑性模型（D-P 模型） ，因

19、此边坡稳定性分析中需要定 义边坡中不同土体的材料属性：容重、弹性模量、泊松比、凝聚力以及摩擦角。 命令：MP GUI：Main MenuPreprocessorMaterial PropsMaterial Models或 Main MenuSolutionLoad Step OptsOtherChange Mat PropsMaterial Models 进行边坡稳定性分析计算时，采用强度折减法来实现。首先选取初始折减系数 F，然后 对边坡土体材料强度系数进行折减，折减后凝聚力以及摩擦角分别见式 4-1 和式 4-2。(4-FCC 1)(4-Ftantan2)和为边坡土体的初始凝聚力和摩擦角。

20、C 对和进行折减，输入边坡模型计算，若收敛，则此时边坡是稳定的；继续增C大折减系数 F，直到程序恰好不收敛，此时的折减系数即为稳定或安全系数。 4.2.2 建立模型和划分网格 创建好物理环境，就可以建立模型。在进行边坡稳定性分析时，需要建立模拟边坡土体ANSYS10.0 土木工程应用实例分析的 PLANE82 单元。在建立好的模型各个区域内指定特性（单元类型、选项、实常数和材料 性质等）以后，就可以划分有限元网格了。 通过 GUI 为模型中的各区赋予特性： 1、选择 Main Menu Preprocessor Meshing Mesh Attributes Picked Areas 2、点击

21、模型中要选定的区域。 3、在对话框中为所选定的区域说明材料号、实常数号、单元类型号和单元坐标系号。 4、重复以上三个步骤，直至处理完所有区域。 通过命令为模型中的各区赋予特性： ASEL（选择模型区域） MAT （说明材料号） REAL（说明实常数组号） TYPE（指定单元类型号） ESYS（说明单元坐标系号）4.2.3 施加约束和荷载在施加边界条件和荷载时，既可以给实体模型（关键点、线、面）也可以给有限元模型 （节点和单元）施加边界条件和荷载。在求解时，ANSYS 程序会自动将加到实体模型上的 边界条件和载荷转递到有限元模型上。 边坡稳定性分析中，主要是给边坡两侧和底部施加自由度约束。 命令

22、：D 施加荷载包括自重荷载以及边坡开挖荷载。 4.2.4 求解接着就可以进行求解，ANSYS 程序根据现有选项的设置，从数据库获取模型和载荷信 息并进行计算求解，将结果数据写入到结果文件和数据库中。 命令：SOLVE GUI：Main MenuSolutionSolve Current LS 4.2.5 后处理后处理的目的是以图和表的形式描述计算结果。对于边坡稳定性分析中，进入后处理器 后，查看边坡变形图和节点的位移、应力和应变。随着强度折减系数的增大，边坡的水平位 移增大，塑性应变急剧发展，塑性区发展形成一个贯通区域时，计算不收敛，认为边坡发生第 4 章 ANSYS 边坡工程应用实例分析了破

23、坏。通过研究位移、应变和塑性区域，来综合判断边坡的稳定性。 命令：POST1 GUI： Main Menu General Postproc 4.2.6 补充说明边坡的失稳破坏定义有很多种，对于采用弹塑性计算模型的边坡，需要综合考虑以下因 素: （1）把有限元计算的收敛与否作为一个重要的衡量指标，边坡处于稳定状态，计算收 敛，边坡破坏时，边坡不收敛。（2）边坡失稳的同时还表现出位移急剧增加。（3）边坡失稳总是伴随着塑性变形的明显增加和塑性区的发展，塑性区的发展状况反 映了边坡是否处于稳定状态。 此外，采用弹塑性有限元法进行计算，它具有独特的优势： (1)弹塑性分析假定岩体为弹塑性材料，岩体在受

24、力初期处于弹性状态，达到一定的屈服 准则后，处于塑性状态。采用弹塑性模型更能反应岩体的实际工作状态。 (2)岩体所承受的荷载超过材料强度时，就会出现明显的滑移破坏面。因此，弹塑性计算 不需要假定破坏面的形状和位置，破坏面根据剪应力强度理论自动形成。当整个边坡破坏时， 就会出现明显的塑性区。 (3 )能综合考虑边坡的局部失稳和整体失稳破坏。ANSYS10.0 土木工程应用实例分析4.3 ANSYS 边坡稳定性实例分析4.3.1 实例描述米米米米米米边坡围岩1弹性材料边坡围岩2弹塑性材料图4-2 边坡模型 边坡实例选取国内某矿，该边坡考虑弹性和塑性两种材料，边坡尺寸如图4-2所示。分 析目的是对该

25、边坡进行稳定性计算分析，以判断其稳定性和计算出安全系数，该边坡围岩材 料属性见表4-1。 表4-1 边坡模型围岩参数类别弹性模量 /Gpa泊松比 v容重/ 3/mkN内聚力 /Mpa摩擦角（。）围岩2（弹塑性）300.2525000.942围岩1（弹性）310.242700-对于像边坡这样纵向很长的实体，计算模型可以简化为平面应变问题。假定边坡所承受 的外力不随Z轴变化，位移和应变都发生在自身平面内。对于边坡变形和稳定性分析，这种 平面假设是合理的。实测经验表明，边坡的影响范围在2倍坡高范围，因此本文计算区域为第 4 章 ANSYS 边坡工程应用实例分析边坡体横向延伸2倍坡高，纵向延伸3倍坡高

26、。两侧边界水平位移为零，下侧边界竖向位移为 零。弹性有限元的计算模型如图4-2所示。 采用双层模型，模型上部为理想弹塑性材料，下部为弹性材料，左右边界水平位移 为零，下边界竖向位移为零。 双层模型考虑土体的弹塑性变形，其塑性区的发展，应力的分布更符合实际情况。 考虑双层模型，塑性区下部的单元可以产生一定的垂直变形和水平变形，基本消除 了由于边界效应在边坡下部出现的塑性区，更好地模拟了边坡的变形和塑性区的发 展。4.3.2 GUI 操作方法4.3.2.1 创建物理环境1) 在【开始】菜单中依次选取【所有程序】/【ANSYS10.0】/【ANSYS Product Launcher】 ，得到“10

27、.0ANSYS Product Launcher”对话框。 2）选中【File Management】 ，在“Working Directory”栏输入工作目录“D:ansysexample4-1” ，在“Job Name”栏输入文件名“Slope” 。 3）单击“RUN”按钮，进入 ANSYS10.0 的 GUI 操作界面。 4）过滤图形界面：Main Menu Preferences，弹出“Preferences for GUI Filtering”对话 框，选中“Structural”来对后面的分析进行菜单及相应的图形界面过滤。 5）定义工作标题：Utility Menu File Ch

28、ange Title，在弹出的对话框中输入“Slope stability Analysis” ，单击“OK” ，如图 4-3。图 4-3 定义工作标题6）定义单元类型：a.定义 PLANE82 单元：Main Menu Preprocessor Element Type Add/Edit/Delete，弹出 一个单元类型对话框，单击“Add”按钮。弹出如图 4-4 所示对话框。在该对话框左面滚动ANSYS10.0 土木工程应用实例分析栏中选择“Solid” ，在右边的滚动栏中选择“Quad 8node 82” ，单击“Apply” ，就定义了“PLANE82”单元。图 4-4 定义 PLAN

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

30、元选项“K3”为“Plane strain”来设定本实例分析采取平面 应变模型进行分析。因为边坡是纵向很长的实体，故计算模型可以简化为平面应变 问题。 8 节点 PLANE82 单元每个节点有 UX 和 UY 两个自由度，比 4 节点 PLANE42 单元具有 更高的精确性，对不规则网格适应性更强。第 4 章 ANSYS 边坡工程应用实例分析7）定义材料属性 a.定义边坡围岩 1 材料属性：Main Menu Preprocessor Material Props Material Models，弹出“Define Material Model Behavior”对话框，如图 4-6 所示。图

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

32、。 再次在图 4-6 中右边的栏中连续双击“StructuralNonlinear InelasticNon-metal plasticitydrucker-prager”后，又弹出一个如图 4-9 所示对话框。在“Cohesion”栏添入边坡 围岩材料 1 的内聚力“0.9E6” ，在“Fric Angle” 栏添入边坡内摩擦角“42” ，单击“OK” 按钮。ANSYS10.0 土木工程应用实例分析图 4-7 线弹性材料模型对话框 图 4-8 材料密度输入对话框图 4-9 定义边坡材料 1DP 模型对话框b.定义边坡围岩 2 材料属性：在图 4-6 对话框中，单击“Material New

33、Model” ， 弹出 一个“Define Material ID”对话框，在“ID”栏后面输入材料编号“2” ，单击“OK”按钮。 弹出一个定义材料模型对话框对话框，选中“Material Model Number 2” ，和定义边坡围岩 1 材料一样，在右边的栏中连续双击“Structural Linear ElasticIsotropic”后，又弹出一个“Linear Isotropic Properties for Material Number 2”对话框，在该对话框中“EX”后面的输 入栏输入“3.2E10” ，在“PRXY” 后面的输入栏输入“0.24” ，单击“OK” 。再选中

34、“Density”并双击，弹出一个“Density for Material Number 2”对话框，在“DENS”后面的 栏中输入隧道围岩材料的密度“2700” ，再单击“OK”按钮，弹出一个定义材料模型对话框。c.复制边坡围岩 1 材料性质：在图 4-6 对话框中，用鼠标点击“Editcopy.” ，弹出一第 4 章 ANSYS 边坡工程应用实例分析个“Copy Material Model”对话框，如图 4-10 所示。在“from Material number”栏后面的下 拉菜单中选取“1” ，在“TO Material number”栏后面输入“3” ，单击“Apply”按钮。又

35、弹 出如土 4-10 所示对话框，然后依次在“TO Material number”栏后面输入“4” 、 “5” 、 “6” 、“7” 、 “8” “9” 、 “10” ， “11” 、 “12” 、 “13” ，每输入一个数，就单击“Apply”按钮一次。图 4-10 复制本构模型对话框最后得到 10 个复制围岩 1 的边坡材料本构模型，如图 4-11 所示。图 4-11 定义强度折减后材料模型对话框ANSYS10.0 土木工程应用实例分析图 4-12 定义强度折减系数 F=1.2 时围岩材料对话框 d.定义 10 个强度折减后材料本构模型：首先定义强度折减系数 F=1.2 后边坡围岩材料模

36、 型，在图 4-11 对话框中，在鼠标依次双击“Material Model Number 3/Drucper-Prager” 。弹出 一个“Drucker- Prager Material Number 3” ，如图 4-12 所示，在“Cohesion”栏添入强度折减 系数 F=1.2 后边坡围岩材料 1 的内聚力“0.75E6” ，在“Fric Angle” 栏添入折减后边坡内摩 擦角“37.7” ，单击“OK”按钮。 用相同方法定义强度折减系数分别为： F=1.4、F=1.6、F=1.8、F=2.0、F=2.2、F=2.4、F=2.6、F=2.8、F=3.0 的边坡围岩材料本构模 型。

37、 定义强度折减后本构模型目的是为了分析边坡稳定性。 强度折减就是降低内聚力和摩擦角，根据式 4-1 和式 4-2 进行折减。 4.3.2.2 建立模型和划分网格1）创建边坡线模型 a.输入关键点：Main Menu Preprocessor Modeling CreateKeypointsIn Active CS，弹 出“Creae Keypoints in Active Cooedinate System”对话框，如图 4-13 所示。在“NPT keypoint number”栏后面输入“1” ，在“X，Y，Z Location in active CS”栏后面输入 “（0，0，0） ”

38、，单击“Apply”按钮，这样就创建了关键点 1。再依次重复在“NPT keypoint number”栏后面输入“2、3、4、5、6、7、8、9” ，在对应“X，Y，Z Location in active CS”栏后面输入“（-800，0，0） 、 （-800，-800，0） 、 （-800，-1200，0） 、 （1200，- 1200，0） 、 （1200，-800，0） 、 （1200，0，0） 、 （1200，378，0） 、 （430，378，0） ，最后单击“OK”按钮。第 4 章 ANSYS 边坡工程应用实例分析图 4-13 在当前坐标系创建关键点对话框 b.创建边坡线模型：

39、Main Menu Preprocessor Modeling CreateLinesStraight line，弹 出“Creae straight lines”对话框，用鼠标依次点击关键点 1、2，单击“Apply”按钮，这样 就创建了直接 L1，同样分别连接关键点“2、3” ， “3、4” ， “4、5” ， “5、6” ， “6、7” ，“7、4” ， “7、8” ， “8、9” ， “9、1” ， “9、2” ，最后单击“OK”按钮，就得到边坡线模型， 如图 4-14 所示。图 4-14 边坡线模型 3）创建边坡面模型 a.打开面编号显示：Utility Menu PlotCtrls

40、 Numbering，弹出“Plot Numbering Controls”对话框，如图 4-15 所示。选中“Aares Numbers” 选项，后面的文字由“off”变 为“on”,单击“OK”关闭窗口。ANSYS10.0 土木工程应用实例分析图 4-15 打开面编号对话框图4-16 边坡面模型b.创建边坡面模型：Main Menu Preprocessor Modeling CreateAreas Arbitraryby line，弹出一个“Create Area by lines”对话框，在图形中选取线L4、L5、L3和L11，点第 4 章 ANSYS 边坡工程应用实例分析击“Appl

41、y ”按钮，就生成了边坡弹性材料区域面积A1；再依次用鼠标在图形中选取线 L1、L2、L6、L10、和L11，点击“Apply ”按钮，就生成了边坡塑性材料区域面积A2；再 依次用鼠标在图形中选取线L7、L8、L10和L9，点击“OK ”按钮，就生成了边坡开挖掉区 域面积A3。最后得到边坡模型的面模型，如图4-16所示。4）划分边坡围岩 2 单元网格a. .给边坡围岩 2 赋予材料特性：Main Menu Preprocessor Meshing MeshTool，弹出“MeshTool”对话框，如图 4-17 所示。在“Element Attributes”后面的下拉式选择栏中选择 “Are

42、as” ，按“Set”按钮，弹出一个“Areas Attributes”面拾取框，在图形界面上拾取边坡 围岩 2 区域，单击拾取框上的“OK”按钮，又弹出一个如图 4-18 所示的“Areas Attributes”对话框，在“Material number”后面的下拉式选择栏中选取“2” ，在“Element type number ” 后面的下拉式选择栏中选取“2 PLANE82” ，单击“Apply” 。图 4-17 网格划分工具栏 图 4-18 定义单元属性对话框b.设置网格划分份数：在图 4-17 工具栏中“Size Control”栏，用鼠标点击“lines”后面 的“Set” ，

43、弹出一个选择对话框，用鼠标在图形选择线 L3 和 L5。弹出一个“Element Sizes on Picked Lines” ，对话框，如图 4-19 所示，在“No of element division”栏后面输入“5” ，ANSYS10.0 土木工程应用实例分析单击“Apply”按钮，再选择线 L4 和 L11，又弹出图 4-19 对话框，在“No of element division”栏后面输入“26” ，单击“OK”按钮。 c.划分单元网格：在图 4-17 网格划分工具栏中单击“Mesh”按钮，弹出一个拾取面积 对话框，拾取面积 A1，单击拾取框上的“OK”按钮，生成边坡围岩 2

44、 单元网格。图 4-19 设置网格份数对话框 5）划分边坡围岩 1 单元网格第 4 章 ANSYS 边坡工程应用实例分析图 4-20 选取线对话框 图 4-21 设置网格分数对话框 a.设置网格份数：Main MenuPreprocessorMeshingSize CntrlsManualSizeLayersPicked Lines，弹出一个“Set Layer Controls”对话框，如图 4-20 所示，用鼠标选取线 L1、L1 和 L6 单击“OK”按钮。弹出一个“Area Layer Mesh Control on Picked lines”对话框，如图 4-21 所示，在“No o

45、f line division”栏后面输入“10” ，单击 “OK”按钮。 相同方法设置线 L8 和 L10 分割份数为 16；设置线 L7 和 L9 线的分割份数为 12。b. 给边坡围岩 1 赋予材料特性：Main Menu Preprocessor Meshing MeshTool，弹出“MeshTool”对话框，如图 4-17 所示。在“Element Attributes”后面的下拉式选择栏中选择 “Areas” ，按“Set”按钮，弹出一个“Areas Attributes”面拾取框，在图形界面上拾取面 A2 和 A3，单击拾取框上的“OK”按钮，又弹出一个 “Areas Attr

46、ibutes”对话框，在“Material number”后面的下拉式选择栏中选取“1” ，在“Element type number ” 后面的下 拉式选择栏中选取“2 PLANE82” ，单击“Apply” 。c.划分单元网格：在图 3-97 网格划分工具栏中单击“Mesh”按钮，弹出一个拾取面积 对话框，拾取围岩，单击拾取框上的“OK”按钮，生成边坡围岩 1 单元网格。 最后得到边坡模型单元网格，如图 4-22 所示。ANSYS10.0 土木工程应用实例分析图 4-22 边坡模型单元网格 4.3.2.3 施加约束和荷载1）给边坡模型施加约束 a.给边坡模型两边施加约束：执行 Main M

47、enuSolutionDefine LoadsApplyStructural Displacementon Nodes，弹出在节点上施加位移约束对话框，用鼠标选取隧道模型两侧边界 上所有节点，单击“OK”按钮。弹出“Apply U,ROT on Nodes”对话框，如图 4-23 所示， 在“DOFS to be constrained”栏后面中选取“UX” ，在“Apply as”栏后面的下拉菜单中选取“Constant value” ，在“Displacement value”栏后面输入“0”值，然后单击“OK”按钮。第 4 章 ANSYS 边坡工程应用实例分析图 4-23 给模型两侧施加

48、位移约束对话框 图 4-24 给模型底部施加位移约束对话框 b.给模型底部施加约束：执行 Main MenuSolutionDefine Loads ApplyStructural Displacementon Nodes，弹出在节点上施加位移约束对话框，用鼠标选取隧道模型底部边界 上所有节点，单击“OK”按钮。弹出图 4-24 所示对话框，在“DOFS to be constrained”栏 后面中选取“UX、UY” ，在“Apply as”栏后面的下拉菜单中选取“Constant value” ，在“Displacement value”栏后面输入“0”值，然后单击“OK”按钮。 节点选择

49、，可以先选择节点上线，再选择附在线上的节点。2）施加重力加速度：Main MenuSolutionDefine LoadsApplyStructuralInertiaGravity，弹出“Apply(Gravitational)Acceleration”对话框， 如图 4-25 所示。只需在“Global Cartesian Y-comp”栏后面输入重力加速度值“9.8”就可以， 单击“OK”按钮，就完成了重力加速度的施加。图 4-25 施加重力加速度对话框这时就可以得到施加约束和重力加速度后隧道有限元模型，如图 4-25 所示。ANSYS10.0 土木工程应用实例分析图 4-25 施加约束和重力荷载后边坡模型 4.3.2.4 求解1）求解设置 a.指定求解类型：Main MenuSolutionAnalysis TypeNew Analy