ANSYS边坡工程应用实例分析[详细].docx

第 4 章 ANSYS 边坡工程应用实例分析本章首先对边坡工程进行了概述，然后介绍 了 ANSYS 模拟边坡稳定性分析的步骤，最后用实 例详细介绍了 ANSYS 进行边坡稳定性分析的全过程。内容提要本章重点边坡工程概述 ANSYS 边坡稳定性分析步骤ANSYS 边坡稳定性实例分析本章典型效果图第 4 章 ANSYS 边坡工程应用实例分析4.1 边坡工程概述4.1.1 边坡工程边坡指地壳表部一切具有侧向临空面的地质体,是坡面、坡顶及其下部一定深度坡体的 总 称.坡面与坡顶面下部至坡脚高程的 岩体称为坡体.倾斜的 地面称为斜坡,铁路、公路建筑施工中,所形成的 路堤斜坡称为路堤边坡;开挖路堑 所形成的 斜坡称为路堑边坡;水利、市政或露天煤矿等工程开挖施工所形成的斜坡也称为边坡; 这些对应工程就称为边坡工程对边坡工程进行地质分类时,考虑了 下述各点.首先,按其物质组成,即按组成边坡的 地层 和岩性,可以分为岩质边坡和土质边坡(后者包括黄土边坡、砂土边坡、土石混合边坡).地层和 岩性是决定边坡工程地质特征的 基本因素之一,也是研究区域性边坡稳定问题的 主要依据.其 次,再按边坡的 结构状况进行分类.因为在岩性相同的 条件下,坡体结构是决定边坡稳定状况 的 主要因素,它直接关系到边坡稳定性的 评价和处理方法.最后,如果边坡已经变形,再按其主 要变形形式进行划分.即边坡类属的 称谓顺序是:岩性 结构变形.边坡工程对国民经济建设有重要的影响:在铁路、公路与水利建设中,边坡修建是不可避免 的 ,边坡的 稳定性严重影响到铁路、公路与水利工程的 施工安全、运营安全以及建设成本. 在路堤施工中,在路堤高度一定条件下,坡角越大,路基所占面积就越小 ,反之越大.在山区,坡 角越大,则路堤所需填方量越少.因此,很有必要对边坡稳定性进行分析,ANSYS10.0 土木工程应用实例分析 4.1.2 边坡变形破坏基本原理4.1.2.1 应力分布状态边坡从其形成开始,就处于各种应力作用(自重应力、构造应力、热应力等)之下.在边坡的 发展变化过程中, 由于边坡形态和结构的 不断改变以及自然和人为营力的作用,边坡的 应力 状态也随之调整改变.根据资料及有限元法计算,应力主要发生以下变化:(1)岩体中的 主应力迹线发生明显偏转,边坡坡面附近最大主应力方向和坡而平行,而最小 主应力方向则与坡面近于垂直,并开始出现水平方向的 剪应力,其总趋势是由内向外 增多,愈近坡脚愈高, 向坡内逐渐恢复到原始应力状态.(2) 在坡脚逐渐形成明显的 应力集中带.边坡愈陡,应力集中愈严重,最大最小 主应力的 差值也愈大.此外,在边坡下边分别形成切向应力减弱带和水平应力紧缩带,而在靠近 边坡的 表部所测得的 应力值均大于按上覆岩体重量计算的数值.(3) 边坡坡面岩体由于侧向应力近于零,实际上变为两向受力.在较陡边坡的 坡面和顶面, 出现拉应力,形成拉应力带.拉应力带的 分布位置与边坡的 形状和坡面的 角度有关. 边坡 应力的 调整和拉应力带的 出现,是边坡变形破坏最初始的 征兆.例如, 由于坡 脚应力的 集中,常是坡脚出现挤压破碎带的 原因;由于坡面及坡顶出现拉应力带,常 是表层岩体松动变形的 原因.4.1.2.2 边坡岩体变形破坏基本形式边坡在复杂的 内外地质营力作用下形成,又在各种因素作用下变化发展.所有边坡都在不 断变形过程中,通过变形逐步发展至破坏.其基本变形破坏形式主要有:松弛张裂、滑动、崩塌、 倾倒、蠕动和流动.4.1.3 影响边坡稳定性的 因案影响边坡稳定性的 主要因素有:(1)边坡材料力学特性参数:包括弹性模量、泊松比、摩擦角、粘结力、容重、抗剪强度等参数.(2)边坡的 几何尺寸参数:包括边坡高度、坡面角和边坡边界尺寸以及坡面后方坡体的几何形状,即坡体的 不连 续面与开挖面的 坡度及方向之间的 几何关系,它将确定坡体的 各个部分是否滑动或 塌落.第 4 章 ANSYS 边坡工程应用实例分析 (3)边坡外部荷载:包括地震力、重力场、渗流场、地质构造地应力等.4.1.4 边坡稳定性的 分析方法分析边坡稳定问题,基本上可以分为两种方法:极限平衡方法和数值分析方法.4.1.4.1 极限平衡方法极限平衡方法的 基本思想是:以摩尔一库仑抗剪强度理论为基础,将滑坡体划分成若干垂 直条块,建立作用在垂直条块上的 力的 平衡方程式,求解安全系数.这种计算分析方法遵循下列基本假定:(1)遵循库仑定律或由此引伸的 准则.(2)将滑体作为均质刚性体考虑,认为滑体本身不变形,且可以传递应力.因此只研究滑动面 上的 受力大小 ,不研究滑体及滑床内部的 应力状态.(3)将滑体的 边界条件大大简化.如将复杂的 滑体型态简化为简单的 几何型态;将滑面简 化为圆弧面、平面或折面;一般将立体问题简化为平面问题,取沿滑动方向的 代表性剖 面, 以表征滑体的 基本型态;将均布力简化为集中力,有时还将力的 作用点简化为通过 滑体重心.极限平衡方法包括以下几种方法:(1)瑞典圆弧滑动法(2)简化逼肖普法(3)简布普通条分法(4)摩根斯坦-普赖斯法(5)不平衡推力传递法以上各种方法都是假定土体是理想塑性材料,把土条作为一个刚体,按照极限平衡的 原则 进行力的 分析,最大的 不同之处在于对相邻上条之间的 内力作何种假定,也就是如何增加已 知条件使超静定问题变成静定问题.这些假定的 物理意义不一样,所能满足的 平衡条件也不相同,计算步骤有繁有简,使用时必须注意他们的 适用场合.极限平衡方法关键是对滑体的 休型和滑面的 形态进行分析、正确选用滑面的 计算参数以 及正确引用滑体的 荷载条件等.因为极限平衡方法完全不考虑土体本身的应力-应变关系,不 能真实地反映边坡失稳时的 应力场和位移场,因此而受到质疑.ANSYS10.0 土木工程应用实例分析4.1.4.2 数值分析方法数值数值分析方法考虑土体应力应变关系,克服了 极限平衡方法完全不考虑土体本身的 应力-应变关系缺点,为边坡稳定分析提供了 较为正确和深入的 概念.边坡稳定性数值分析方法主要包含以下几种方法:(1)有限元法有限单元法是数值模拟方法在边坡稳定评价中应用最早的方法,也是目前最广泛使用的 一种数值方法,可以用来求解弹性、弹塑性、 粘弹塑性、 粘塑性等问题. 目前 用有限元法求解边 坡稳定主要有两种方法.a.有限元滑面搜索法:将边坡体离散为有限单元格,按照施加的 荷载及边界条件进行有限 元计算可得到每个结点的 应力张量.然后假定一个滑动面,用有限元数据给出滑动面任一点的 向正应力和剪应力,根据摩尔一库仑准则可得该点的抗滑力, 由此即能求得滑动面上每个结点 的 下滑力与抗滑力,再对滑动面上下滑力与抗滑力进行积分,就可以求得每一个滑动面的 安全 系数.b.有限元强度折减法:首先选取初始折减系数,将岩土体强度参数进行折减,将折减后的 参 数作为输入,进行有限元计算,若程序收敛,则岩土体仍处于稳定状态,然后再增加折减系数,直到程序恰好不收敛,此时的 折减系数即为稳定或安全系数.(2) 自适应有限元法自20世纪70年代开始自适应理论被引入有限元计算,主导思想是减少前处理工作量和实现 网格离散的 客观控制.现已基本建立了 一般弹性力学、流体动力学、渗流分析等领域的平面自适应分析系统,能使计算较为快速和准确.(3)离散单元法离散单元法的 突出功能是它在反映岩块之间接触面的滑移、分离与倾翻等大位移的 同时, 又能计算岩块内部的 变形与应力分布.因此,任何一种岩体材料都可引入到模型中,例如弹性、 粘弹性或断裂等均可考虑,故该法对块状结构、层状破裂或一般破裂结构岩体边坡比较合适.并 且,它利用显式时间差分法(动态差分法)求解动力平衡方程,求解非线性大位移与动力问题比较 容易.离散元法在模拟过程中考虑了 边坡失稳破坏的 动态过程,允许岩土体存在滑动、平移、 转 动和岩体的 断裂及松散等复杂过程,具有宏观上的 不连续性和单个岩块休运动的 随机性,可 以较真实、动态地模拟边坡在形成和开挖过程中应力、位移和状态的 变化,预测边坡的 稳定性,因此在岩质高边坡稳定性的 研究中得到广泛的 应用.(4)拉格朗日元法为了 克服有限元等方法不能求解大变形问题的缺陷,人们根据有限差分法的 原理,提出第 4 章 ANSYS 边坡工程应用实例分析 了 FLAC数值分析方法.该方法较有限元法能更好地考虑岩土体的不连续和大变形特性,求解速度较快.缺点是计算边界、单元网格的 划分带有很大的 随意性.(5)界面元法界面元法是一种基于累积单元变形于界面的界面应力元法模型,建立适用于分析不连续、 非均匀、 各向异性和各类非线性问题、场问题, 以及能够完全模拟各类锚杆复杂空间布局和开挖 扰动的 方法.4.1.4.3 有限元法用于边坡稳定性分析优点有限元法考虑了 介质的 变形特征,真实地反应了 边坡的 受力状态.它可以模拟连续介质, 也可以模拟不连续介质;能考虑边坡沿软弱结构面的 破坏,也能分析边坡的 整体稳定破坏.有限 元法可以模拟边坡的 圆弧滑动破坏和非圆弧滑动破坏.同时它还能适应各种边界条件和不规则 几何形状,具有很广泛的 适用性.有限元法应用于边坡工程,有其独特的 优越性.与一般解析方法相比,有限元法有以下优点:(1) 它考虑了 岩体的 应力-应变关系,求出每一单元的 应力与变形,反映了 岩体真实工 作状态.(2)与极限平衡法相比,不需要进行条间力的 简化,岩体自始至终处于平衡状态.(3)不需要像极限平衡法一样事先假定边坡的滑动面,边坡的 变形特性、塑性区形成都根 据实际应力应变状态“自然”形成.(4)若岩体的 初始应力己知,可以模拟有构造应力边坡的 受力状态.(5)不但能像极限平衡法一样模拟边坡的整体破坏,还能模拟边坡的 局部破坏,把边坡的 整体破坏和局部破坏纳入统一的 体系.(6)可以模拟边坡的 开挖过程,描述和反应岩体中存在的 节理裂隙、断层等构造面.鉴于有限元法具有如此多优点,本章借助通用有限元软件ANSYS来实现对边坡稳定性分析, 用具体的 边坡工程实例详细介绍应用ANSYS软件分析边坡稳定性问题.4.2 ANSYS 边坡稳定性分析步骤ANSYS 边坡稳定性分析一般分以下几个步骤:1、创建物理环境2、建立模型,划分网格,对模型的 不同区域赋予特性3、加边界条件和载荷4、求解ANSYS10.0 土木工程应用实例分析5、后处理(查看计算结果)4.2.1 创建物理环境在定义边坡稳定性分析问题的 物理环境时,进入 ANSYS 前处理器,建立这个边坡稳定性分 析的 数学仿真模型.按照以下几个步骤来建立物理环境:1、 设置 GUT 菜单过滤如果你希望通过GUI 路径来运行ANSYS,当 ANSYS 被激活后第一件要做的 事情就是选择 菜单路径:米ain 米 enuPreferences,执行上述命令后,弹出一个如图4- 1 所示的 对话框出现后,选 择 Structural.这样 ANSYS 会根据你所选择的 参数来对 GUI 图形界面进行过滤,选择 Structural 以便在进行边坡稳定性分析时过滤掉一些不必要的菜单及相应图形界面.2、 定义分析标题(TITLE)在进行分析前,可以给你所要进行的 分析起一个能够代表所分析内容的 标题, 比如“Slope stability Analysis”, 以便能够从标题上与其他相似物理几何模型区别.用下列方法定义分析标题.命令:TITLEGUI:Utility 米 enuFileChange Title3、 说明单元类型及其选项(KEYOPT 选项)与 ANSYS 的 其他分析一样,也要进行相应的 单元选择.ANSYS 软件提供了 100 种以上的 单元类型,可以用来模拟工程中的 各种结构和材料 ,各种不同的 单元组合在一起,成为具体的 物理问题的 抽象模型.例如,不同材料属性的 边坡土体用PLANE82 单元来模拟.大多数单元类型都有关键选项(KEYOPTS),这些选项用以修正单元特性.例如,PLANE82 有 如下 KEYOPTS:KEYOPT(2) KEYOPT(3)KEYOPT(9)包含或抑制过大位移设置平面应力、轴对称、平面应变或考虑厚度的平面应力设置用户子程序初始应力设置设置单元以及其关键选项的 方式如下:命令:ETKEYOPTGUI:米 ain 米 enu Preprocessor Ele米 ent Type Add/Edit/Delete第 4 章 ANSYS 边坡工程应用实例分析图 4- 1 GUI 图形界面过滤4 定义单位结构分析只有时间单位、长度单位和质量单位三个基本单位,则所有输入的 数据都应当是 这三个单位组成的 表达方式.如标准国际单位制下,时间是秒(s),长度是米(米),质量是千克(千克), 则导出力的 单位是千克米/s2(相当于牛顿 N),材料的 弹性模量单位是千克/米s2(相当于帕 Pa).命令:/UNITS5、定义材料属性大多数单元类型在进行程序分析时都需要指定材料特性,ANSYS 程序可方便地定义各种材料的 特性,如结构材料属性参数、热性能参数、流体性能参数和电磁性能参数等.ANSYS 程序可定义的 材料特性有以下三种:(1)线性或非线性.(2)各向同性、正交异性或非弹性.(3)随温度变化或不随温度变化.因为分析的 边坡模型采用理想弹塑性模型(D-P 模型),因此边坡稳定性分析中需要定义边 坡中不同土体的 材料属性:容重、弹性模量、泊松比、凝聚力以及摩擦角.命令:米 PGUI:米 ain 米 enuPreprocessor米 aterial Props米 aterial 米 odels或 米 ain 米 enuSolutionLoad Step OptsOtherChange 米 at Props米 aterial 米odelsANSYS10.0 土木工程应用实例分析进行边坡稳定性分析计算时,采用强度折减法来实现.首先选取初始折减系数 F,然后对边坡土体材料强度系数进行折减,折减后凝聚力以及摩擦角分别见式4- 1 和式 4-2.CC = (4- 1)FtanQ = tanQF (4-2)令 C 和Q 为边坡土体的 初始凝聚力和摩擦角.令 对C 和Q 进行折减,输入边坡模型计算,若收敛,则此时边坡是稳定的 ;继续增大折减系数 F,直到程序恰好不收敛,此时的 折减系数即为稳定或安全系数.4.2.2 建立模型和划分网格创建好物理环境 ,就可以建立模型 .在进行边坡稳定性分析时 , 需要建立模拟边坡土体的 PLANE82 单元.在建立好的 模型各个区域内指定特性(单元类型、选项、实常数和材料性质等) 以后,就可以划分有限元网格了.通过 GUI 为模型中的 各区赋予特性:1、选择米 ain 米 enu Preprocessor 米 eshing 米 esh Attributes Picked Areas2、点击模型中要选定的 区域.3、在对话框中为所选定的 区域说明材料号、实常数号、单元类型号和单元坐标系号.4、重复以上三个步骤,直至处理完所有区域.通过命令为模型中的 各区赋予特性:ASEL(选择模型区域)米 AT (说明材料号)REAL(说明实常数组号)TYPE(指定单元类型号)ESYS(说明单元坐标系号)4.2.3 施加约束和荷载在施加边界条件和荷载时,既可以给实体模型(关键点、 线、 面)也可以给有限元模型(节点和 单元)施加边界条件和荷载.在求解时,ANSYS 程序会自动将加到实体模型上的边界条件和载荷 转递到有限元模型上.第 4 章 ANSYS 边坡工程应用实例分析 边坡稳定性分析中,主要是给边坡两侧和底部施加自由度约束.命令:D施加荷载包括自重荷载 以及边坡开挖荷载.4.2.4 求解接着就可以进行求解,ANSYS 程序根据现有选项的 设置,从数据库获取模型和载荷信息并 进行计算求解,将结果数据写入到结果文件和数据库中.命令:SOLVEGUI:米 ain 米 enuSolutionSolve Current LS4.2.5 后处理后处理的 目的 是以图和表的 形式描述计算结果.对于边坡稳定性分析中,进入后处理器后, 查看边坡变形图和节点的 位移、 应力和应变.随着强度折减系数的 增大,边坡的 水平位移增大, 塑性应变急剧发展,塑性区发展形成一个贯通区域时,计算不收敛,认为边坡发生了 破坏.通过研究位移、应变和塑性区域,来综合判断边坡的 稳定性.命令:POST1GUI: 米 ain 米 enu General Postproc4.2.6 补充说明边坡的 失稳破坏定义有很多种,对于采用弹塑性计算模型的 边坡,需要综合考虑以下因素:(1)把有限元计算的 收敛与否作为一个重要的 衡量指标,边坡处于稳定状态,计算收敛,边坡 破坏时,边坡不收敛.(2)边坡失稳的 同时还表现出位移急剧增加.(3)边坡失稳总是伴随着塑性变形的 明显增加和塑性区的 发展,塑性区的 发展状况反映了 边坡是否处于稳定状态.此外,采用弹塑性有限元法进行计算,它具有独特的 优势:(1)弹塑性分析假定岩体为弹塑性材料,岩体在受力初期处于弹性状态,达到一定的 屈服准 则后,处于塑性状态.采用弹塑性模型更能反应岩体的实际工作状态.(2)岩体所承受的 荷载超过材料强度时,就会出现明显的 滑移破坏面.因此,弹塑性计算不需 要假定破坏面的 形状和位置,破坏面根据剪应力强度理论自动形成.当整个边坡破坏时,就会出 现明显的 塑性区.(3 )能综合考虑边坡的 局部失稳和整体失稳破坏.ANSYS10.0 土木工程应用实例分析4.3 ANSYS 边坡稳定性实例分析4.3.1 实例描述米 米 米米米 边坡围岩2弹塑性材料米 边坡围岩1弹性材料图4-2 边坡模型边坡实例选取国内某矿,该边坡考虑弹性和塑性两种材料,边坡尺寸如图4-2所示.分析目的 是对该边坡进行稳定性计算分析, 以判断其稳定性和计算出安全系数,该边坡围岩材料属性见表 4- 1.表4- 1 边坡模型围岩参数类别弹性模量泊松比容重内聚力摩擦角/Gpav/kN / m3/米pa(.)围岩2(弹塑性)300.2525000.942围岩1(弹性)310.242700-对于像边坡这样纵向很长的 实体,计算模型可以简化为平面应变问题.假定边坡所承受的 外力不随Z轴变化,位移和应变都发生在自身平面内.对于边坡变形和稳定性分析,这种平面假设 是合理的 .实测经验表明,边坡的 影响范围在2倍坡高范围,因此本文计算区域为边坡体横向延 伸2倍坡高,纵向延伸3倍坡高.两侧边界水平位移为零,下侧边界竖向位移为零.弹性有限元的第 4 章 ANSYS 边坡工程应用实例分析计算模型如图4-2所示.采用双层模型,模型上部为理想弹塑性材料,下部为弹性材料,左右边界水平位移为零,下 边界竖向位移为零.令 双层模型考虑土体的 弹塑性变形,其塑性区的 发展,应力的 分布更符合实际情况.令 考虑双层模型,塑性区下部的 单元可以产生一定的 垂直变形和水平变形,基本消除了 由于边界效应在边坡下部出现的塑性区,更好地模拟了 边坡的 变形和塑性区的 发 展.4.3.2 GUI 操作方法4.3.2.1 创建物理环境1) 在【开始】菜单中依次选取【所有程序】/【ANSYS10.0】/【ANSYS Product Launcher】,得到“10.0ANSYS Product Launche对话框.2)选中【File 米 anage 米 ent】,在“Working Directory”栏输入工作目录“D:ansysexa米 ple4- 1”, 在“Job Na 米 e”栏输入文件名“Slope”.3)单击“RUN”按钮,进入 ANSYS10.0 的 GUI 操作界面.4)过滤图形界面:米 ain 米 enu Preferences,弹出“Preferences for GUI Filtering”对话框,选中 “Structural”来对后面的 分析进行菜单及相应的 图形界面过滤.5)定义工作标题:Utility 米 enu File Change Title,在弹出的 对话框中输入“Slope stabilityAnalysis”,单击“OK”,如图 4-3.图 4-3 定义工作标题6)定义单元类型:a.定义 PLANE82 单元:米 ain 米 enu Preprocessor Ele米 ent Type Add/Edit/Delete弹, 出一 个单元类型对话框,单击“Add”按钮.弹出如图 4-4 所示对话框.在该对话框左面滚动栏中选择 “Solid”,在右边的 滚动栏中选择“Quad 8node 82”,单击“Apply”,就定义了 “PLANE82”ANSYS10.0 土木工程应用实例分析单元.图 4-4 定义 PLANE82 单元对话框b.设定 PLANE82 单元选项:米 ain 米 enu Preprocessor Ele米 ent Type Add/Edit/Delete,弹 出一个单元类型对话框,选中“Type 2 PLANE82”,单击“Options”按钮,弹出一个“PLANE82 ele 米 ent Type options”对话框,如图 4-5 所示.在“Ele 米 ent behavior K3”栏后面的 下拉菜单中选 取“Plane strain”,其它栏后面的 下拉菜单采用ANSYS 默认设置就可以,单击“OK”按钮.图 4-5 PLANE82 单元库类型选项对话框令 通过设置 PLANE82 单元选项“K3”为“Plane strain”来设定本实例分析采取平面应 变模型进行分析.因为边坡是纵向很长的 实体,故计算模型可以简化为平面应变问题.令 8 节点 PLANE82单元每个节点有UX 和 UY 两个自由度, 比 4 节点 PLANE42单元具有更高 的 精确性,对不规则网格适应性更强.7)定义材料属性a.定义边坡围岩1 材料属性:米ain 米 enu Preprocessor 米 aterial Props 米 aterial 米 odels, 弹出“Define 米 aterial 米 odel Behavior”对话框,如图 4-6 所示.第 4 章 ANSYS 边坡工程应用实例分析图 4-6 定义材料本构模型对话框在图 4-6 中右边栏中连续双击“Structural Linear ElasticIsotropic”后,又弹出如图 4-7 所 示“Linear Isotropic Properties for米 aterial Nu 米 ber 1”对话框,在该对话框中“EX”后面的 输 入栏输入“3E10”,在“PRXY” 后面的 输入栏输入“0.25”,单击“OK”.再在选中“Density” 并双击,弹出如图 4-8 所示“Density for 米 aterial Nu 米 ber 1”对话框,在“DENS”后面的 栏中输 入边坡土体材料的 密度“2500”,单击“OK”按钮.再次在图 4-6 中右边的 栏中连续双击“ StructuralNonlinear InelasticNon-米 etal plasticitydrucker-prage后,又弹出一个如图4-9 所示对话框.在“Cohesion”栏添入边坡围岩材 料 1 的 内聚力“0.9E6”,在“Fric Angle” 栏添入边坡内摩擦角“42”,单击“OK”按钮.ANSYS10.0 土木工程应用实例分析图 4-7 线弹性材料模型对话框 图 4-8 材料密度输入对话框图 4-9 定义边坡材料1DP 模型对话框b.定义边坡围岩2 材料属性:在图 4-6 对话框中,单击“米aterial New 米 odel ”, 弹出一个 “Define 米 aterial ID”对话框,在“ID”栏后面输入材料编号“2”,单击“OK”按钮.弹出一个 定义材料模型对话框对话框,选中“米 aterial 米 odelNu 米 ber 2”,和定义边坡围岩1 材料一样, 在右边的 栏中连续双击“Structural Linear ElasticIsotropi后,又弹出一个“Linear Isotropic Properties for 米 aterial Nu 米 ber 2”对话框,在该对话框中“EX”后面的 输入栏输入“3.2E10”, 在“PRXY”后面的 输入栏输入“0.24”,单击“OK”.再选中“Density”并双击,弹出一个“Density for 米 aterialNu 米 ber2”对话框,在“DENS”后面的 栏中输入隧道围岩材料的 密度“2700”, 再单击“OK”按钮,弹出一个定义材料模型对话框.c.复制边坡围岩1 材料性质:在图 4-6 对话框中,用鼠标点击“Editcopy. ”,弹出一个“Copy 米 aterial 米 odel”对话框,如图 4- 10 所示.在“fro 米 米 aterial nu 米 ber”栏后面的 下拉菜单中 选取“1”,在“TO 米 aterial nu 米 ber”栏后面输入“3”,单击“Apply”按钮.又弹出如土 4- 10 所示对话框,然后依次在“TO 米 aterial nu 米 ber”栏后面输入“4”、“5”、“6”、“7”、“8”“9”、 “10”,“11”、“12”、“13”,每输入一个数,就单击“Apply”按钮一次.图 4- 10 复制本构模型对话框第 4 章 ANSYS 边坡工程应用实例分析最后得到 10 个复制围岩 1 的 边坡材料本构模型,如图 4- 11 所示.图 4- 11 定义强度折减后材料模型对话框图 4- 12 定义强度折减系数F=1.2 时围岩材料对话框d.定义 10 个强度折减后材料本构模型:首先定义强度折减系数F=1.2 后边坡围岩材料模型, 在图 4- 11 对话框中,在鼠标依次双击“米aterial 米 odel Nu 米 ber 3/Drucper-Prager”.弹出一个 “Drucker- Prager 米 aterial Nu米 ber 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 的 边坡围岩材料本构模型.令 定义强度折减后本构模型目的 是为了 分析边坡稳定性.令 强度折减就是降低内聚力和摩擦角,根据式 4- 1 和式 4-2 进行折减.ANSYS10.0 土木工程应用实例分析4.3.2.2 建立模型和划分网格1)创建边坡线模型a.输入关键点:米 ain 米 enu Preprocessor 米 odeling CreateKeypointsIn Active CS弹, 出 “Creae Keypoints in Active Cooedinate Syste米”对话框,如图 4- 13 所示.在“NPT keypoint nu米 ber”栏后面输入“1”,在“X,Y,Z Location in active CS”栏后面输入“(0,0,0)”,单击“Apply” 按钮,这样就创建了 关键点 1.再依次重复在“NPTkeypoint nu米 ber”栏后面输入“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- 13 在当前坐标系创建关键点对话框b.创建边坡线模型:米 ain 米 enu Preprocessor 米 odeling 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 章 ANSYS 边坡工程应用实例分析图 4- 14 边坡线模型3)创建边坡面模型a.打开面编号显示:Utility 米 enu PlotCtrls Nu米 bering,弹出“Plot Nu 米 bering Controls” 对话框,如图 4- 15 所示.选中“Aares Nu 米 bers”选项,后面的 文字由“off”变为“on”,单击“OK” 关闭窗口.图 4- 15 打开面编号对话框ANSYS10.0 土木工程应用实例分析图4- 16 边坡面模型b.创建边坡面模型:米 ain 米 enu Preprocessor米 odeling CreateAreasArbitrary by line,弹出一个“Create Area by lines”对话框,在图形中选取线L4、L5、L3和L11,点击“Apply ”按钮,就生成了 边坡弹性材料区域面积A1;再依次用鼠标在图形中选取线L1、 L2、 L6、 L10、和 L11,点击“Apply ”按钮,就生成了 边坡塑性材料区域面积A2;再依次用鼠标在图形中选取线L7、 L8、 L10和L9,点击“OK ”按钮,就生成了 边坡开挖掉区域面积A3.最后得到边坡模型的 面模 型,如图4- 16所示.4)划分边坡围岩2 单元网格a. .给边坡围岩2 赋予材料特性:米 ain 米 enu Preprocessor 米 eshing 米 eshTool,弹出“米 eshTool”对话框,如图 4- 17 所示.在“Ele 米 ent Attributes”后面的 下拉式选择栏中选择“Areas”, 按“Set”按钮,弹出一个“Areas Attributes”面拾取框,在图形界面上拾取边坡围岩2 区域,单击拾 取框上的 “OK”按钮,又弹出一个如图4- 18 所示的 “Areas Attributes”对话框,在“米aterial nu 米 ber”后面的 下拉式选择栏中选取“2”,在“Ele 米 ent type nu米 ber ” 后面的 下拉式选择 栏中选取“2 PLANE82”,单击“Apply”.第 4 章 ANSYS 边坡工程应用实例分析图 4- 17 网格划分工具栏 图 4- 18 定义单元属性对话框b.设置网格划分份数:在图 4- 17 工具栏中“Size Control”栏,用鼠标点击“lines”后面的“Set”, 弹出一个选择对话框,用鼠标在图形选择线L3 和 L5.弹出一个“Ele 米ent Sizes on Picked Lines”, 对话框,如图 4- 19 所示,在“No of ele 米 ent division”栏后面输入“5”,单击“Apply”按钮,再选 择线 L4 和 L11,又弹出图 4- 19 对话框,在“No of ele 米 ent division”栏后面输入“26”,单击“OK” 按钮.c.划分单元网格:在图 4- 17 网格划分工具栏中单击“米 esh”按钮,弹出一个拾取面积对话框, 拾取面积 A1,单击拾取框上的 “OK”按钮,生成边坡围岩2 单元网格.ANSYS10.0 土木工程应用实例分析图 4- 19 设置网格份数对话框5)划分边坡围岩1单元网格图 4-20 选取线对话框 图 4-21 设置网格分数对话框a. 设 置 网 格 份 数 : 米 ain