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 为基础长度)，称为柔性指数，它反映了相对刚度对内力分布的影响。
ANSYS的基本步骤讲解
ANSYS的基本步骤讲解1.问题定义：在使用ANSYS之前，您需要明确要解决的问题。
定义问题包括确定您要分析的物理现象，所使用的材料，边界条件和所需的结果等。
2.创建几何模型：根据问题定义，您需要创建一个几何模型来表示分析的对象。
ANSYS 提供了各种建模工具，可以用于创建二维和三维的几何形状。
您可以使用ANSYS自带的CAD工具或导入其他CAD软件创建的模型。
3.划分网格：对几何模型进行网格划分是进行仿真分析的关键步骤。
ANSYS提供了各种网格划分工具，可以根据需要选择不同的划分技术和网格密度。
良好的网格划分可以提高仿真的准确性和效率。
4.材料属性定义：根据问题定义，您需要为模型中的不同部分定义适当的材料属性。
ANSYS提供了一个材料库，可以选择多种不同的材料，并根据需要定义其属性参数，如弹性模量，热导率等。
5.添加约束条件：在仿真中，机械结构通常受到约束条件的限制。
您需要添加适当的约束条件来代表物理世界中的限制。
ANSYS提供了各种约束选项，包括固定支撑、自由支撑、弹簧等。
通过添加这些约束条件，您可以更准确地模拟实际场景。
6.应用载荷：在仿真中，您需要明确地定义作用在模型上的载荷。
这可能是一个力，一个压力，一个温度或者其他物理效应。
您需要在模型的相关位置上添加适当的载荷。
ANSYS提供了各种载荷选项，可以精确描述应用的载荷。
7.设置仿真参数：在进行仿真之前，您需要设置一些模拟参数。
这些参数可以控制仿真过程和结果的精度。
例如，您可以设置时间步长、迭代次数、收敛标准等。
8.运行仿真：配置完所有参数后，可以开始运行仿真。
ANSYS将根据所设置的仿真参数对模型进行计算。
这可能需要一段时间，具体取决于模型的大小和复杂程度。
9.结果分析：仿真运行结束后，可以进行结果分析。
ANSYS提供了丰富的结果可视化和后处理工具，可以帮助您更好地理解和解释结果。
您可以查看模型的位移、应力、应变、温度分布等结果。
10.结果验证：最后，您需要对仿真结果进行验证。
ANSYS桥梁工程应用实例分析(详细)(图文)
本章介绍桥梁结构的模拟分析。
桥梁是一种重要的工程结构，精确分析桥梁结构在各种受力方式下的响应有较大的工程价值。
模拟不同类型的桥梁需要不同的建模方法，分析内容包括静力分析、动荷载响应分析、施工过程分析等等。
在本章中着重介绍桁架桥、刚架桥和斜拉桥三种类型桥梁。
内容 提要 第6章 ANSYS 桥梁工程应用实例分析本章重点结构分析具体步骤结构静力分析 桁架结构建模方法 结构模态分析本章典型效果图6.1 引言ANSYS通用有限元软件在土木工程应用分析中可发挥巨大的作用。
我们用它来分析桥梁工程结构,可以很好的模拟各种类型桥梁的受力、施工工况、动荷载的耦合等。
ANSYS程序有丰富的单元库和材料库,几乎可以仿真模拟出任何形式的桥梁。
静力分析中,可以较精确的反应出结构的变形、应力分布、内力情况等；动力分析中,也可精确的表达结构的自振频率、振型、荷载耦合、时程响应等特性。
利用有限元软件对桥梁结构进行全桥模拟分析,可以得出较准确的分析结果。
本章介绍桥梁结构的模拟分析。
作为一种重要的工程结构,桥梁的精确分析具有较大的工程价值。
桥梁的种类繁多,如梁桥、拱桥、钢构桥、悬索桥、斜拉桥等等,不同类型的桥梁可以采用不同的建模方法。
桥梁的分析内容又包括静力分析、施工过程模拟、动荷载响应分析等。
可以看出桥梁的整体分析过程比较复杂。
总体上来说,主要的模拟分析过程如下:(1) 根据计算数据,选择合适的单元和材料,建立准确的桥梁有限元模型。
(2) 施加静力或者动力荷载,选择适当的边界条件。
(3) 根据分析问题的不同,选择合适的求解器进行求解。
(4) 在后处理器中观察计算结果。
(5) 如有需要,调整模型或者荷载条件,重新分析计算。
桥梁的种类和分析内容众多,不同类型桥梁的的分析过程有所不同,分析侧重点也不一样。
在这里仅仅给出大致的分析过程,具体内容还要看具体实例的情况。
6.2 典型桥梁分析模拟过程6.2.1 创建物理环境建立桥梁模型之前必须对工作环境进行一系列的设置。
最新ANSYS基坑弹性地基梁全程序即详解汇总
A N S Y S基坑弹性地基梁全程序即详解/prep7L1=30 !设置变量L2=30h=-25K, 1, 0, 0, 0,K, 2, L1, 0, 0,K, 3, L1, L2, 0,K, 4, 0, L2, 0,KWPAVE, 1 !将工作平面原点定义在 1号点RECTNG, 0, L1, 0, L2,wpro, , -90, !将工作平面绕X轴Z到Y方向90度RECTNG, 0, L1, 0, -h,KWPAVE, 4 !将工作平面原点定义在 4号点RECTNG, 0, L1, 0, -h,wpro, , ,90 !将工作平面绕y轴x到z方向90度RECTNG, 0, L2, 0, -h,KWPAVE, 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-Damper MPTEMP,,,,,,,, !删除系统中已存在的温度表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, WALL, 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命令中的MAT即材料性能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, , , , , !约束土弹簧单元端点的所有位移ALLSSAVE/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, , , *enddo TIME, 2 alls solve。
弹性地基梁.
图3.3 弹性地基梁的微元分析
1.弹性地基梁的挠度曲线微分方程式
为建立 y x 应满足的挠曲微分方程，在梁中截取一微段 d x ，考察该段 的平衡有： Y 0, 得：
Q (Q dQ) kydx q( x) d x 0
dQ 化简得： ky q( x) dx
（3-2）
由微分方程理论知，上述方程的通解由四个线性无关的特解组合而成。为 寻找四个线性无关的特解，令 y e rx 并代入上式有：
K EI
4
或
K cos i sin EI
4
由复数开方根公式得：
rk 4
令
4
K 2k 2k i sin COS k 0,1,2,3 EI 4 4 Kb K ， 若地基梁宽度为b，则有 4 EI EI
二阶
(8)
y py qy 0
(9)
设想(9)有形式解 y = erx (为什么?)
代入得 (r2 + pr + q ) erx = 0
故有
r2 + pr + q = 0
(10)
(10)式称为(9)的特征方程, 分三种情形讨论 (i) = p2– 4q > 0, (10)有两个不等实根 r1, r2.
优点：
1、地基的连续整体性；2、几何物理上简化模型
缺点：
1、地基土非连续；2、地基土非均质；
图3.2 弹性地基梁的受力和变形
3.弹性地基梁的挠度曲线微分方程 式及其初参数解
基本假设：
除局部弹性地基模型假设外，还需作假设： （1）地基梁在外荷载作用下产生变形的过程中，梁底面与地基 表面始终紧密相贴，即地基的沉陷或隆起与梁的挠度处处相等；
ansys模仿地基开挖的竖向弹性地基梁的方法
运用ANSYS模拟多撑式深基坑开挖的研究王一鸣张小平（南京航空航天大学土木工程学院，江苏南京，210016）摘要：以竖向弹性地基梁理论为基础，以ANSYS为平台，采用APDL语言建立了多撑式深基坑开挖的计算模型。
将其应用于上海丽晶苑大厦基坑支护结构的计算，取得了满意的结果。
同时以此工程为原型，对该模型的主要影响因素进行了分析,提出要对墙背土压力进行修正，以进一步完善该模型。
关键词：弹性地基梁ANSYS深基坑土压力中国分类号：TU398 文献标识码：A 文章编号：The research of the application of ANSYS in emulation of the excavation of deep foundation pits with multiple bracesWang Yiming Zhang Xiaopin(School of Civil Engineering, Nan Jing University of Aeronautics and Astronautics, 210016, China ) Abstract:With the theory of beams on elastic foundation for foundation, using ANSYS system, an approach of using finite element method for simulation the behavior of deep excavation pit with multiple knightheads is presented in this paper. This technique is applied to design the deep excavation work of the LiJingyuan mansion of Shanghai, and the perfect reliability and the high efficiency are proved by the design results. What’s more ,through the analysis of this work , the factors that have important influence to the calculation results are studied, and Put forward that the earth pressure behind the pile should be modified to perfect that model. Key words:beams on elastic foundation,ANSYS, deep foundation pit, Earth pressure前言对于多撑式地下连续墙的内力计算分析，在计算机普及之前，以日本工程界提出的“山肩邦男法”、“弹性法”、“弹塑性法”等解析法为主。
ANSYS基坑弹性地基梁全程序即详解6页word文档
/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-Damper MPTEMP,,,,,,,, !删除系统中已存在的温度表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。
ANSYS程序简介及基本使用方法
b）确定粱单元的截面参数 main menu >preprocessor >real constants >add/edit/delete
c）确定材料参数 ANSYS main menu >preprocessor >material props >material models 在弹出的“Define Material Model Behavior”界面中双击 structural>linear>elastic>isotropic
谢谢大家！
网络教育专业课
有 限 元 原理
9 程序简介及基本使用方法
本部分包括如下内容： 1. 有限元分析的基本方法 2. ANSYS简介 3. 梁的受力与变形分析－简支粱的变形分析 4. 弹性力学三维问题的分析－支架的受力分析
5 温度场分析
6 讨论（计算结果分析）
1. 有限元分析的基本方法 1）建立实际工程问题的计算模型 利用几何、载荷的对称性简化模型 建立等效模型 2）选择适当的分析工具 侧重考虑以下几个方面： 多物理场耦合问题 大变形 网格重划分 3）前处理（Preprocessing） 建立几何模型（Geometric Modeling，自下而上，或 基本单元组合） 有限单元划分（Meshing）与网格控制
在“Plot Numbering Controls”界面中，可以选择显示几 何实体和网格的编号。
在“Size and Shape”对话框中选择单元的显示方式。在 [/SHRINK]域的下拉列表中选择将单元缩小显示的百分比； 在[/ESHAPE]域中选择，是否用实参数来显示单元。
3）施加边界条件及加载 a) 在支撑位置施加位移约束
结点支反力如下图所示：
弹性地基梁程序说明书
弹性地基梁程序研究设计人员：郑楠指导老师：刘川顺老师武汉大学水利水电学院2003年6月9日1弹性地基梁基本原理假定地基为半无限的连续弹性体，应用弹性理论进行计算。
本程序中，弹性地基梁采用的是链杆法原理，并考虑边荷载的作用。
其具体表述如下：链杆法的基本作法是：将地基梁分成若干段，并且在每一梁段的中心设置一根不可压缩的刚性链杆将梁和地基联系起来，这样，就可以把一个原有无限多个支承的地基梁的计算问题，转化为一个支承在有限个可沉陷支座上的连续梁的计算问题。
通过梁的平衡条件和梁与地基的变形协调条件，就可以建立求解超静定结构的典型方程，从而计算出各个链杆的内力，进而计算地基梁各计算截面的内力。
采用混合法建立典型方程。
为了使建立的联立方程的工作简化，混合法采用悬臂梁作为基本结构，也就是除了把n个竖向力链杆切断外，还在梁的一端附加一个竖向链杆和一个刚臂，控制梁在这端的移动和转动，使其成为固定端。
因而，基本未知量不但包括n个未知力（链杆内力），而且还包括两个位移（梁端的挠度和转角）。
如图所示：基本结构和原结构相比，可以知道基本结构中每一个竖向链杆的切口处的X链杆切口处的相对位移应该等于零的条件，即相对线位移应该为零，根据k=0（设使切口张开的位移为正）可以成立下列方程：k0001=∆+--∑=kp k ni kiia y X ϕσ式中： ki σ ：对单位力分别作用在悬臂梁和地基上时, k X 链杆切口处引起的相对位移kp ∆ ：外荷载使悬臂梁在k X 链杆切口处方向引起的位移0y ：悬臂梁固定端处的竖向位移，以0y 向下为正，所以沿k X 方向的位移为-0y0ϕ：悬臂梁固定端处的转角，以顺时针方向为正，所以由0ϕ引起的沿方向的位移为-0ϕk a显然，对于n 个竖向链杆的切口，可以列出n 个方程。
另外，通过否定梁端附加 链杆和刚臂的存在（即0,000==M R ），还可以列出如下的两平衡方程，即由∑=0y ，有 021=-+++++∑P X X X X n i由00=∑M ，有∑∑=-=01M X a ni i i式中： ∑P ：外荷载在竖直方向投影代数和，以向下为正；∑M ：外荷载对梁左端弯矩的代数和，以顺时针为正。
ANSYS基坑弹性地基梁全程序即详解
(/ =20kPa/prep7L1=30 !设置变量L2=30h=-25K, 1,0, 0, 0,K, 2, L1, 0, 0,K, 3, L1, L2, 0,K, 4, 0, L2, 0,KWPAVE, 1 !将工作平面原点定义在1号点RECTNG, 0, L1, 0, L2,wpro, , -90, !将工作平面绕X轴Z到Y方向90度RECTNG, 0, L1, 0, -h,KWPAVE, 4 !将工作平面原点定义在4号点RECTNG, 0, L1, 0, -h,wpro, , ,90 !将工作平面绕y轴x到z方向90度RECTNG, 0, L2, 0, -h,KWPAVE, 3 !将工作平面原点定义在3号点RECTNG, 0, L2, 0, -h,AGLUE, all !粘结所有面ALLS NUMCMP , ALL! 直接生成节点 *D0, i, 1, L1-1 CSYS, 0!激活默认笛卡尔坐标系ET, 1, SHELL43元） !ET ，ITYPE ，Ename ，KOPT1，~，KOPT6 ，INOPR （定义单ET, 2, COMBIN14!K0PT1〜K0PT6为元素特性编码!shell43 4 节点塑性大应变单元 !C0MBIN14弹簧 -阻尼器 Spring-Damper MPTEMP,,,,,,,,! 删除系统中已存在的温度表 MPTEMP, 1, 0! 定义一个温度表 MPDATA, EX, 1, , 2.4E10! 指定与温度相应的材料性能数据 弹性模量MPDATA, PRXY , 1, , 0.15!主泊松比 ESIZE, 1, 0!指定单元边长 AMESH, ALL! 划分面生成面单元 NSEL, S, L0C, Z, 0ESLN, SNSEL, S, L0C, Z, -1ESLN, UCM, STRUT, ELEMallsCMSEL, U, STRUTCM, WALL, ELEM! 选择一组节点子集创建新集 ! 选择已选节点上的单元 选择 z 坐标值为 -1 的 --- ! 从已选集中删除此时剩下只支撑板 ! 将选择集命名 STRUT 生成元件 !all sel 全选 ! 去除 STRUT 元件 ! 将选择集命名 wall 生成元件 NSEL, S, L0C, X, 0.1, L1-0.1NPL0TNSEL, R, L0C, Y, 0ESLN, SNSEL, S, L0C, Y, 1ESLN, UENSYM, , , , ALLNSEL, S, L0C, Y, 0.1, L2-0.1NPL0TNSEL, R, L0C, X, 0ESLN, SNSEL, S, L0C, X, 1ESLN, UENSYM , , , , ALL ! 选择一组节点子集创建新集 ! 显示节点 ! 从当前集选择一组节点子集 ! 从已选集中选择 ! 从当前集选择一组节点子集 ! 从已选集中删除 ! 反转壳单元法线方向 ! 选择一组节点子集创建新集 ! 显示节点 ! 从当前集选择一组节点子集 ! 从已选集中选择 ! 从当前集选择一组节点子集! 从已选集中删除! 反转壳单元法线方向 ! 所有实体进行重新编号从1 到 29进行循环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, 2MAT , 1REAL, 0 ESYS, 0!设置单元类型属性指示器!MP 命令中的MAT 即材料性能! 材料实常数!材料坐标系统属性指示器EN, 100000+i, 100000+2*i-1, 100000+2*i ! 根据给定的单元号和节点号生成单元*enddoNUMCMP , ALL !所有实体进行重新编号*DO, i, 1, L1-1CSYS, 0N, 100000+2*i-1, i, L2, -1 N, 100000+2*i, i, L2-2.5, *enddo ! 激活默认笛卡尔坐标系1*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 $ MA T , 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 由单元组成的元件! 从当前集选择一组节点子集 深度-1 以下 !从已选集中选择 !选择单元类型号为 2 的单元 !单位面积内受的力随深度增加而增加 !对已存在单元进行修改R, 101, 0.6, 0.6, 0.6, 0.6, ,R, 102, 0.1, 0.1, 0.1, 0.1, ,CMSEL, S, W ALLALLSALLSEL, ALLNUMMRG , NODE, 0.01, 0.01, LOWNUMCMP , ALL!选择所有实体 ! 节点合并 距离小于 0.01 则同保留编 !所有实体进行重新编号 底的点 *DO, i, 1, 25NSEL, S, LOC, Z, -1*iESLN, SNSEL, R, TYPE, , 2R, i, m1*i*b*h, , , EMODIF, ALL, REAL, i, *enddo !识别号 +实常数EMODIF, ALL, REAL, 101, CMSEL, S, STRUT EMODIF, ALL, REAL, 102, !将实常数101 组赋给墙! 赋值给支撑NSEL, S, LOC, Z-25D, ALL, , , , , , UZ, , , , ,!约束墙底竖向位移CMSEL, S, SPRING NSLE, SCMSEL, S, W ALL NSLE, UD, ALL, , , , , , ALL, , , , ,!选择土弹簧单元!以下命令从已选弹簧集合中选! 在已选集中选墙单元! 从集合中删除刚选择的单元，即与墙有关的单元!约束土弹簧单元端点的所有位移ALLS SAVE/SOLUALLSANTYPE, STA TIC, NEW NROPT, FULL*AFUN, DEGQ=2.0E4GAMA1=0.9E4FAI1=20C1=10e3m1=1500e3b=1h=1! 分析模式静力!指定计算模式! 指定角度单位为度!Q 为超载!浮重度! 内摩擦角!粘聚力!比例系数!单元宽度!墙体单元高度! 无支撑开挖1m 坑内水位-1.0 坑外0NSEL, S, LOC, Z, 0, -0.9 ESLN, SCMSEL, R, wall NSLE, S!0 到0.9 的位置即开挖面以上!以下命令在当前集里选择! 在当前集选wll 单元!以下命令在当前集里选择*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, 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,乙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, 25 z=i-1 !开挖了 1m 要在被动区减去相应的土压力 !水压力计算 ! 总应力R, i, K, , , *enddoTIME, 1 alls solve! 载荷步1NSEL, 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+SA2 pressure=SA SLZER=ZMAXSLOPE=(SA-PA)/(LOCZ2-LOCZ1) SFGRAD, PRES, 0, z, SLZER, SLOPE SFE, all, 1, PRES, , -PA, , , allsNSEL, S, LOC, Z, -5.1, -24.9ESLN, SGMSEL, R, wall*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, S2 的单元ESEL, R, TYPE, , 2 !选择单元类型号为EKILL, ALL*do, i, 6, 25z=i-5K=m1*z*b*hR, i, K, , ,*enddoTIME, 2allssolve。
- 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。