ANSYS基坑弹性地基梁全程序即详解
ansys教程——梁解析
梁
5. 梁
• 梁单元 是线单元，用来创建3-D结构的一维理想体。 • 梁单元比实体和壳单元更有效，经常用于工业领域中:
– 建筑结构 – 桥梁和道路 – 公共交通 (有轨电车, 火车, 公共汽车) – 等
Training Manual
January 30, 2001 Inventory #001443 5-2
梁
...梁网格划分
步骤1: 线属性 • 梁网格划分的线属性包括:
– 材料号 – 横截面号 – 定位关键点
Training Manual
INTRODUCTION TO ANSYS 5.7 - Part 2
• 相对于梁轴线，横截面是怎样定位。
• 必须指定所有横截面类型. • 单个关键点可以分配给多条线 ( 即，不需要为每条线指定单个关键 点 )。
横截面 • 对 BEAM188 和 189单元的完整定义包括对横 截面属性的定义。 • BeamTool提供了方便的操作.
– Preprocessor > Sections > Common Sectns...
Training Manual
INTRODUCTION TO ANSYS 5.7 - Part 2
梁
...梁网格划分
步骤2 :线分隔 • 对 BEAM188 和 189 单元,不主张把整个梁当作一个单元 • 使用 Mesh Tool的 “Size Controls”指定想要的线分隔数 (或用 LESIZE 命令).
Training Manual
January 30, 2001 Inventory #001443 5-13
梁
...梁网格划分
Training Manual
ANSYS 入门教程 - 结构的弹性稳定性分析
ANSYS 入门教程- 结构的弹性稳定性分析2011-01-09 15:06:42| 分类：默认分类| 标签：|字号大中小订阅第7 章结构弹性稳定分析7.1 特征值屈曲分析的步骤7.2 构件的特征值屈曲分析7.3 结构的特征值屈曲分析一、结构失稳或结构屈曲：当结构所受载荷达到某一值时，若增加一微小的增量，则结构的平衡位形将发生很大的改变，这种现象叫做结构失稳或结构屈曲。
结构稳定问题一般分为两类:★第一类失稳：又称平衡分岔失稳、分枝点失稳、特征值屈曲分析。
结构失稳时相应的载荷可称为屈曲载荷、临界载荷、压屈载荷或平衡分枝载荷。
★第二类失稳：结构失稳时，平衡状态不发生质变，也称极值点失稳。
结构失稳时相应的载荷称为极限载荷或压溃载荷。
●跳跃失稳：当载荷达到某值时，结构平衡状态发生一明显的跳跃，突然过渡到非邻近的另一具有较大位移的平衡状态。
可归入第二类失稳。
★结构弹性稳定分析= 第一类稳定问题ANSYS 特征值屈曲分析（Buckling Analysis）。
★第二类稳定问题ANSYS 结构静力非线性分析，无论前屈曲平衡状态或后屈曲平衡状态均可一次求得，即“全过程分析”。
这里介绍ANSYS 特征值屈曲分析的相关技术。
在本章中如无特殊说明，单独使用的“屈曲分析”均指“特征值屈曲分析”。
7.1 特征值屈曲分析的步骤①创建模型②获得静力解③获得特征值屈曲解④查看结果一、创建模型注意三点：⑴仅考虑线性行为。
若定义了非线性单元将按线性单元处理。
刚度计算基于初始状态（静力分析后的刚度），并在后续计算中保持不变。
⑵必须定义材料的弹性模量或某种形式的刚度。
非线性性质即便定义了也将被忽略。
⑶单元网格密度对屈曲载荷系数影响很大。
例如采用结构自然节点划分时（一个构件仅划分一个单元）可能产生100% 的误差甚至出现错误结果，尤其对高阶屈曲模态的误差可能更大，其原因与形成单元应力刚度矩阵有关。
经验表明，仅关注第1 阶屈曲模态及其屈曲载荷系数时，每个自然杆应不少于 3 个单元。
弹性地基梁法(全面版)资料
弹性地基梁法（全面版）资料弹性地基梁法整体式平底板的平面尺寸远较厚度为大，可视为地基上的受力复杂的一块板。
目前工程实际仍用近似简化计算方法进行强度分析。
一般认为闸墩刚度较大，底板顺水流方向弯曲变形远较垂直水流方向小，假定顺水流方向地基反力呈直线分布，故常在垂直水流方向截取单宽板条进行内力计算。
按照不同的地基情况采用不同的底板应力计算方法。
相对密度Dr ＞0.5的砂土地基或粘性土地基,可采用弹性地基梁法。
相对密度Dr 0.5的砂土地基，因地基松软，底板刚度相对较大，变形容易得到调整，可以采用地基反力沿水流流向呈直线分布、垂直水流流向为均匀分布的反力直线分布法。
对小型水闸，则常采用倒置梁法。
（一）弹性地基梁法该法认为底板和地基都是弹性体，底板变形和地基沉降协调一致，垂直水流方向地基反力不呈均匀分布（图1），据此计算地基反力和底板内力。
此法考虑了底板变形和地基沉降相协调，又计入边荷载的影响，比较合理，但计算比较复杂。
当采用弹性地基梁法分析水闸闸底板应力时，应考虑可压缩土层厚度T 与弹性地基梁半长L ／2之比值的影响。
当L T 2小于0.25时，可按基床系数法（文克尔假定）计算；当L T 2大于2.0时，可按半无限深的弹性地基梁法计算；当2T ／L 为0.25-2.0时，可按有限深的弹性地基梁计算。
弹性地基梁法计算地基反力和底板内力的具体步骤如下：（1）用偏心受压公式计算闸底纵向(顺水流方向)地基反力。
（2）在垂直水流方向截取单宽板条及墩条，计算板条及墩条上的不平衡剪力。
以闸门槽上游边缘为界，将底板分为上、下游两段，分别在两段的中央截取单宽板条及墩条进行分析，如图1（a ）所示。
作用在板条及墩条上的力有：底板自重（q 1）、水重（q 2）、中墩重（G 1/b i ）及缝墩重（G 2/b i ），中墩及缝墩重中（包括其上部结构及设备自重在内），在底板的底面有扬压力（q 3）及地基反力（q 4），见图1（b ）所示。
ANSYS 入门教程 - 结构的弹性稳定性分析
ANSYS 入门教程- 结构的弹性稳定性分析2011-01-09 15:06:42| 分类：默认分类| 标签：|字号大中小订阅第7 章结构弹性稳定分析7.1 特征值屈曲分析的步骤7.2 构件的特征值屈曲分析7.3 结构的特征值屈曲分析一、结构失稳或结构屈曲：当结构所受载荷达到某一值时，若增加一微小的增量，则结构的平衡位形将发生很大的改变，这种现象叫做结构失稳或结构屈曲。
结构稳定问题一般分为两类:★第一类失稳：又称平衡分岔失稳、分枝点失稳、特征值屈曲分析。
结构失稳时相应的载荷可称为屈曲载荷、临界载荷、压屈载荷或平衡分枝载荷。
★第二类失稳：结构失稳时，平衡状态不发生质变，也称极值点失稳。
结构失稳时相应的载荷称为极限载荷或压溃载荷。
●跳跃失稳：当载荷达到某值时，结构平衡状态发生一明显的跳跃，突然过渡到非邻近的另一具有较大位移的平衡状态。
可归入第二类失稳。
★结构弹性稳定分析= 第一类稳定问题ANSYS 特征值屈曲分析（Buckling Analysis）。
★第二类稳定问题ANSYS 结构静力非线性分析，无论前屈曲平衡状态或后屈曲平衡状态均可一次求得，即“全过程分析”。
这里介绍ANSYS 特征值屈曲分析的相关技术。
在本章中如无特殊说明，单独使用的“屈曲分析”均指“特征值屈曲分析”。
7.1 特征值屈曲分析的步骤①创建模型②获得静力解③获得特征值屈曲解④查看结果一、创建模型注意三点：⑴仅考虑线性行为。
若定义了非线性单元将按线性单元处理。
刚度计算基于初始状态（静力分析后的刚度），并在后续计算中保持不变。
⑵必须定义材料的弹性模量或某种形式的刚度。
非线性性质即便定义了也将被忽略。
⑶单元网格密度对屈曲载荷系数影响很大。
例如采用结构自然节点划分时（一个构件仅划分一个单元）可能产生100% 的误差甚至出现错误结果，尤其对高阶屈曲模态的误差可能更大，其原因与形成单元应力刚度矩阵有关。
经验表明，仅关注第1 阶屈曲模态及其屈曲载荷系数时，每个自然杆应不少于 3 个单元。
梁ansys分析实例讲解课件
ansys软件操作流程
定义弹性模量、泊松比、密度等材料属性 建立几何体
在Model模块中选择“Model”选项卡
ansys软件操作流程
• 创建梁的几何体，输入梁的截面尺寸等信息
ansys软件操作流程
网格划分 在Model模块中选择“Mesh”选项卡 设置网格大小、网格类型等参数
定义约束和载荷
在简支梁的两个端点上定义约束和 载荷。
边界条件与载荷施加
固定约束
在简支梁的两个端点施加固定约 束，以模拟简支边界条件。
均布载荷
在简支梁的跨中施加均布载荷， 以模拟简支梁受到的集中力。
网格划分与求解
网格划分
对简支梁进行网格划分，可以选 择合适的网格密度以提高求解精
度。
求解设置
在“Solution”菜单下进行求解 设置，包括迭代次数、收敛准则
减少模型规模：减小 模型的规模，降低内 存需求；
如何设置ANSYS软 件中的单位制？
升级硬件：增加物理 内存或使用更好的计 算机配置。
a型时， 选择合适的单位制，如米制或英
制；
在模型树中选择模型名称，进入 Model模块，在Model模块中选 择“Model”选项卡，在弹出的 对话框中选择“Units”选项卡
求解设置
进行求解设置，包括迭代次数、收敛标准等。
求解过程
进行求解，得到桥梁结构的应力分布、位移分布等结果。
06
ansys软件操作流程及常见问题 解答
ansys软件操作流程
建立模型 启动ANSYS软件，选择Workbench模式
创建新的模型文件，命名并保存
ansys软件操作流程
Ansys Workbench详解教程全解
2018/10/21
16
向导
作用： 帮助用户设置分析过程中的基本步骤，如选择分析类型、定义材 料属性等基本分析步骤。 显示： 可以通过菜单View中的Windows选项或常用工具条中的图标 控制其显示。
2018/10/21
17
基本操作
创建、打开、保存文档 复制、剪切、粘贴
32
网格控制
整体网格 ： Relevance （-100~100 ） 、 Relevance Center （coarse~ fine ）
局部细化： 支撑处、载荷施加位置、应力变化较大的地方。
2018/10/21
33
网格控制
具体操作：选中结构树的Mesh项，点击鼠标右键，选择Insert，弹出 对网格进行控制的各分项，一般只需设置网格的形式（Method）和单元的 大小（Sizing）。
? y ?
? i? i
i
k e? e ? F e
c. 计算等效节点力
3. 进行单元集成； 4. 得到节点位移；
K? ? F
5. 根据弹性力学公式计算单元应变、应力。
ANSYS Workbench 软件介绍
运行软件 操作界面简介 基本操作 分析流程的各项操作
2018/10/21
9
运行软件
方法一：从CAD软件中进入
2018/10/21
14
属性窗口
属性窗口提供了输入数据的列表， 会根据选取分支的不同自动改变。
白色区域 : 显示当前输入的数据。 灰色区域 : 显示信息数据，不能
被编辑。 黄色区域 : 未完成的信息输入。
2018/10/21
15
图形窗口
模型和结果都将显 示在这个区域中， 包括：
手把手教你学ansys--钢筋混凝土梁
⼿把⼿教你学ansys--钢筋混凝⼟梁⼤家好，我是⽔哥。
⽔哥ansys使⽤经验三年多，既做过重⼤科研项⽬，也做过许多实际项⽬，对ansys的使⽤有⼀定的⼼得体会，本着分享经验的精神，今⽇以⼀个钢筋混凝⼟梁的建模求解过程来简单说明ansys的基本操作步骤。
（我的ansys14.0）总的说来，⽆论⼩项⽬还是⼤项⽬，总体过程⽆⾮定义单元、定义材料、物理模型、有限元模型、加载、约束、求解、查看这⼏个过程，和我们⼯程类的设计软件例如PKPM、SAP2000等基本过程都差不多，只不过最⼤的区别在于ansys的建模实在是蛋疼了⼀些。
废话不多说，以下⾯的⼀根悬臂钢筋混凝⼟梁来教新⼿如何快速进⼊ansys 的⼤堂，每⼀步都有GUI操作，完了之后会有相应部分的命令流，这⾥多说⼀句，⼀个ansys的使⽤⾼⼿必然是⼀个精通apdl编程的能⼿，所以我建议新⼿在学习的时候最好以apdl⼊⼿，GUI操作辅助，这样在学习的时候能节省⼤量时间，⽽不会浪费在GUI毫⽆意义的重复操作上。
此题如下：悬臂梁如下，梁宽200mm，梁顶有两根直径为16的钢筋，钢筋中⼼距梁边的距离为40mm，在梁端附近受集中⼒P=100KN的作⽤.要求对此悬臂梁进⾏完全线弹性分析，结果要显⽰主应⼒迹线。
材料参数：混凝⼟弹性模量为3000MPa,钢筋的弹性模量取200GPa，不考虑材料⾃重。
（建模时注意单位的协调性）⼀、题⽬解读与材料单元定义注意此题要求进⾏完全线弹性分析，此话的意义在于我们可以⽤除solid65 以外的其他实体单元。
在ansys单元中，solid65是专门⽤于模拟钢筋混凝⼟构件的单元，但⽤此单元模拟时，⼀般是需要考虑材料的⾮线性，也即是多⽤于构件的⾮线性分析，并且需要材料的本构模型。
此题要求我们做弹性分析，我们可以⽤link8单元模拟钢筋，⽤solid45单元模拟混凝⼟，注意⾼版本的ansys已经将许多低阶单元合并掉了，以ansys14为例，在link单元中只有180，⽽低阶的link8、link10等已被合并。
10-1 弹性地基梁的解析方法
2. 弹性地基梁法弹性地基梁内力计算：基床系数法和半无限弹性体法。
基床系数法：采用文克勒（Winkler）地基模型，地基由许多互不联系的弹簧所组成，某点的地基沉降仅由该点上作用的压力所产生。
通过求解弹性地基梁的挠曲微分方程，可求出基础梁的内力。
半无限弹性体法：假定地基为半无限弹性体，将柱下条形基础看作放在半无限弹性体表面上的梁，而基础梁在荷载作用下，满足一般的挠曲微分方程。
应用弹性理论求解基本挠曲微分方程，并引入基础与半无限弹性体满足变形协调的条件及基础的边界条件，求出基础的位移和基底压力，进而求出基础的内力。
半无限弹性体法的求解一般采用有限单元法等数值方法。
，根据微分梁单元力的平衡，则：∑Y＝M x w EI -=22d d 由材料力学知，梁的挠曲微分方程为：或2244d d d d xM x w EI -=根据截面剪力与弯矩的相互关系，即则：x x M d dQ d d 22=q bp x w EI +-=44d d q bkw x w EI =+44d d 引入文克勒地基模型及地基沉降s 与基础梁的挠曲变形协调条件，可得：。
w s =kw ks p ==代入上式，可得文克勒地基上梁的挠曲微分方程为：当梁上的分布荷载q =0时，梁的挠曲微分方程变为齐次方程：0d d 44=+bkw x w EI令，称为梁的柔度指标，其单位为(长度)-1。
的倒数值称为特征长度，值愈大，梁对地基的相对刚度愈大。
44EI kb =λλλλ1λ104d d 444=+w x w λ该微分方程的通解为)sin cos ()sin cos (4321x C x C e x C x C e w x x λλλλλλ+++=-于是，梁的挠曲微分方程可进一步写成如下形式：式中C 1、C 2、C 3、C 4为待定参数，根据荷载及边界条件定；为无量纲量，当x ＝L (L 为基础长度)，称为柔性指数，它反映了相对刚度对内力分布的影响。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
/prep7L1=30 !设置变量L2=30h=-25K, 1, 0, 0, 0,K, 2, L1, 0, 0,K, 3, L1, L2, 0,K, 4, 0, L2, 0,KWPA VE, 1 !将工作平面原点定义在1号点RECTNG, 0, L1, 0, L2,wpro, , -90, !将工作平面绕X轴Z到Y方向90度RECTNG, 0, L1, 0, -h,KWPA VE, 4 !将工作平面原点定义在4号点RECTNG, 0, L1, 0, -h,wpro, , ,90 !将工作平面绕y轴x到z方向90度RECTNG, 0, L2, 0, -h,KWPA VE, 3 !将工作平面原点定义在3号点RECTNG, 0, L2, 0, -h,AGLUE, all !粘结所有面ET, 1, SHELL43 !ET，ITYPE，Ename，KOPT1，~，KOPT6，INOPR（定义单元）!KOPT1～KOPT6为元素特性编码!shell43 4 节点塑性大应变单元ET, 2, COMBIN14 !COMBIN14弹簧-阻尼器Spring-DamperMPTEMP,,,,,,,, !删除系统中已存在的温度表MPTEMP, 1, 0 !定义一个温度表MPDATA, EX, 1, , 2.4E10 !指定与温度相应的材料性能数据弹性模量MPDATA, PRXY, 1, , 0.15 !主泊松比ESIZE, 1, 0 !指定单元边长AMESH, ALL !划分面生成面单元NSEL, S, LOC, Z, 0 !选择一组节点子集创建新集ESLN, S !选择已选节点上的单元NSEL, S, LOC, Z, -1 !选择z坐标值为-1的---ESLN, U !从已选集中删除此时剩下只支撑板CM, STRUT, ELEM !将选择集命名STRUT生成元件alls !all sel 全选CMSEL, U, STRUT !去除STRUT元件CM, W ALL, ELEM !将选择集命名wall生成元件NSEL, S, LOC, X, 0.1, L1-0.1 !选择一组节点子集创建新集NPLOT !显示节点NSEL, R, LOC, Y, 0 !从当前集选择一组节点子集ESLN, S !从已选集中选择NSEL, S, LOC, Y, 1 !从当前集选择一组节点子集ESLN, U !从已选集中删除ENSYM, , , , ALL !反转壳单元法线方向NSEL, S, LOC, Y, 0.1, L2-0.1 !选择一组节点子集创建新集NPLOT !显示节点NSEL, R, LOC, X, 0 !从当前集选择一组节点子集ESLN, S !从已选集中选择NSEL, S, LOC, X, 1 !从当前集选择一组节点子集ESLN, U !从已选集中删除ENSYM , , , , ALL !反转壳单元法线方向ALLSNUMCMP, ALL !所有实体进行重新编号!直接生成节点*DO, i, 1, L1-1 ! 从1到29进行循环CSYS, 0 !激活默认笛卡尔坐标系N, 100000+2*i-1, i, 0, -1 ! 节点编号后面为坐标N, 100000+2*i, i, 2.5, -1 ! 节点编号后面为坐标*enddo*DO, i, 1, L1-1 ! 从1到29进行循环!Modeling>Creat>Elements>Elem AttributesTYPE, 2 !设置单元类型属性指示器MAT , 1 !MP命令中的MA T即材料性能REAL, 0 !材料实常数ESYS, 0 !材料坐标系统属性指示器EN, 100000+i, 100000+2*i-1, 100000+2*i !根据给定的单元号和节点号生成单元*enddoNUMCMP , ALL !所有实体进行重新编号*DO, i, 1, L1-1CSYS, 0 !激活默认笛卡尔坐标系N, 100000+2*i-1, i, L2, -1N, 100000+2*i, i, L2-2.5, -1*enddo*DO, i, 1, L1-1TYPE, 2MAT , 1REAL, 0ESYS, 0EN, 100000+i, 100000+2*i-1, 100000+2*i*enddoNUMCMP , ALL !所有实体进行重新编号*DO, i, 1, L2-1CSYS, 0N, 100000+2*i-1, 0, i, -1N, 100000+2*i, 2.5, i, -1*enddo*DO, i, 1, L2-1TYPE, 2MAT , 1REAL, 0ESYS, 0EN, 100000+i, 100000+2*i-1, 100000+2*i*enddoNUMCMP , ALL !所有实体进行重新编号*DO, i, 1, L2-1CSYS, 0N,100000+2*i-1, L1, i, -1N,100000+2*i, L1-2.5, i, -1*enddo*DO, i, 1, L2-1TYPE, 2 $ MAT , 1 $ REAL , 0ESYS, 0EN, 100000+i, 100000+2*i-1, 100000+2*i*enddoNUMCMP , ALL !所有实体进行重新编号ESEL, S, TYPE, , 2 !选择单元类型号为2的单元EGEN, 25, 100000, ALL, , , , , , , , , , -1,ESEL, S, TYPE, , 2 !选择单元类型号为2的单元CM, SPRING, ELEM !生成一个叫SPRING由单元组成的元件ALLSALLSEL, ALL !选择所有实体NUMMRG, NODE, 0.01, 0.01, LOW !节点合并距离小于0.01则同保留编号底的点NUMCMP, ALL !所有实体进行重新编号*DO, i, 1, 25NSEL, S, LOC, Z, -1*i !从当前集选择一组节点子集深度-1以下ESLN, S !从已选集中选择NSEL, R, TYPE, , 2 !选择单元类型号为2的单元R, i, m1*i*b*h, , , !单位面积内受的力随深度增加而增加EMODIF, ALL, REAL, i, !对已存在单元进行修改*enddoR, 101, 0.6, 0.6, 0.6, 0.6, , !识别号+实常数R, 102, 0.1, 0.1, 0.1, 0.1, ,CMSEL, S, WALLEMODIF, ALL, REAL, 101, !将实常数101组赋给墙CMSEL, S, STRUTEMODIF, ALL, REAL, 102, !赋值给支撑NSEL, S, LOC, Z-25 !约束墙底竖向位移D, ALL, , , , , , UZ, , , , ,CMSEL, S, SPRING !选择土弹簧单元NSLE, S !以下命令从已选弹簧集合中选CMSEL, S, WALL !在已选集中选墙单元NSLE, U !从集合中删除刚选择的单元，即与墙有关的单元D, ALL, , , , , , ALL, , , , , !约束土弹簧单元端点的所有位移ALLSSA VE/SOLUALLSANTYPE, STATIC, NEW !分析模式静力NROPT, FULL !指定计算模式*AFUN, DEG !指定角度单位为度Q=2.0E4 !Q为超载GAMA1=0.9E4 !浮重度FAI1=20 !内摩擦角C1=10e3 !粘聚力m1=1500e3 !比例系数b=1 !单元宽度h=1 !墙体单元高度!无支撑开挖1m 坑内水位-1.0 坑外0NSEL, S, LOC, Z, 0, -0.9 !0到0.9的位置即开挖面以上ESLN, S !以下命令在当前集里选择CMSEL, R, wall !在当前集选wll 单元NSLE, S !以下命令在当前集里选择*GET, ZMIN, NODE, , MNLOC, Z, , , , !墙最浅节点处的位置*GET, ZMAX, NODE, , MXLOC, Z, , , , !墙最深节点处的位置LOCZ1=abs (ZMAX) !取绝对值LOCZ2=abs (ZMIN)KA= (TAN(45.0 - FAI1/2))**2 !主动土压力系数!采用水土分算PA1=(Q + GAMA1*LOCZ1)*KA-2.0*C1*SQRT(KA) !最深主动土压力公式* IF , PA1, LT, 0, THEN !去除小于零值的可能PA1=0*ENDIFPA2=10*1E3*LOCZ1 !水压力计算PA=PA1+PA2 !总应力SA1=(Q+GAMA1*LOCZ2)*KA-2*C1*SQRT(KA)SA2=10*1E3*LOCZ2WA=SA2SA=SA1+SA2pressure=SASLZER=ZMAXSLOPE=(SA-PA)/(LOCZ2-LOCZ1) !增长率SFGRAD, PRES, 0, z, SLZER, SLOPE !沿z方向从-1m开始的面力减少SFE, all, 1, PRES, , -PA, , , !alls!开挖面以下加载NSEL, S, LOC, Z, -1.1, -24.9ESLN, SGMSEL, R, wallNSLE, S*GET, ZMIN, NODE, , MNLOC, Z, , , ,*GET, ZMAX, NODE, , MXLOC, Z, , , ,SLZER=ZMAXSLOPE=0SFGRAD, PRES, O, z, SLZER, SLOPESFE, all, 1, PRES, , -pressure, , , !alls !将集合扩大到全集cmsel, S, strut !选择内支撑全部杀死EKILL, ALLallsNSEL, S, LOC, Z, 0.1, -1.1ESLN, SESEL, R, TYPE, , 2 !选择单元类型号为2的单元EKILL, ALL*do, i, 2, 25z=i-1 !开挖了1m 要在被动区减去相应的土压力K=m1*z*b*hR, i, K, , ,*enddoTIME, 1 !载荷步1allssolveNSEL, S, LOC, Z, 0, -4.9ESLN, SCMSEL, R, wallNSLE, S*GET, ZMIN, NODE, , MNLOC, Z, , , ,*GET, ZMAX, NODE, , MXLOC, Z, , , ,LOCZ1=abs (ZMAX)LOCZ2=abs (ZMIN)KA= (TAN(45.0 - FAI1/2))**2PA1=(Q + GAMA1*LOCZ1)*KA-2.0*C1*SQRT(KA)* IF , PA1, LT, 0, THENPA1=0*ENDIFPA2=10*1E3*LOCZ1PA=PA1+PA2SA1=(Q+GAMA1*LOCZ2)*KA-2*C1*SQRT(KA) SA2=10*1E3*LOCZ2WA=SA2SA=SA1+SA2pressure=SASLZER=ZMAXSLOPE=(SA-PA)/(LOCZ2-LOCZ1)SFGRAD, PRES, 0, z, SLZER, SLOPESFE, all, 1, PRES, , -PA, , ,allsNSEL, S, LOC, Z, -5.1, -24.9ESLN, SGMSEL, R, wallNSLE, S*GET, ZMIN, NODE, , MNLOC, Z, , , ,*GET, ZMAX, NODE, , MXLOC, Z, , , ,SLZER=ZMAXSLOPE=0SFGRAD, PRES, O, z, SLZER, SLOPESFE, all, 1, PRES, , -pressure, , ,allscmsel, s, strutealive, ALLallsNSEL, S, LOC, Z, 0.1, -5.1ESLN, SESEL, R, TYPE, , 2 !选择单元类型号为2的单元EKILL, ALL*do, i, 6, 25z=i-5K=m1*z*b*hR, i, K, , ,*enddoTIME, 2allssolve。