隧道衬砌ANSYS强度检算
隧道强度检验ANSYS命令流
TITLE,Mechanical analysis on three spans continuum beam ! 确定分析标题/NOPR !菜单过滤设置/PMETH,OFF,0KEYW,PR_SET,1KEYW,PR_STRUC,1/COM,/COM,Preferences for GUI filtering have been set to display:/COM, Structural !只要结构分析部分菜单!*/PREP7 !* 进入前处理器ET,1,BEAM3 !设置梁单元类型R,1,0.18,0.0054,0.6, , , , !设置几何常数MPTEMP,,,,,,,, !设置材料模型MPTEMP,1,0MPDA TA,EX,1,,30e9 !输入弹性模量MPDA TA,PRXY,1,,0.2 !输入泊松比MPTEMP,,,,,,,,MPTEMP,1,0MPDA TA,DENS,1,,2500 !输入密度!*SA VE !保存数据库K,1,0,0,0, !创建关键点K,2,8,0,0,K,3,20,0,0,K,4,28,0,0,/REP,FAST/PNUM,KP,1 !显示关键点号LSTR, 1, 2 !创建直线LSTR, 2, 3LSTR, 3, 4SA VELESIZE,ALL,1, , , ,1, , ,1, !设置单元大小FLST,2,3,4,ORDE,2FITEM,2,1FITEM,2,-3LMESH,P51X !划分网格/PNUM,ELEM,1!显示单元内容/REPLOTSA VEFINISH/SOLFLST,2,1,3,ORDE,1FITEM,2,1DK,P51X, ,0, ,0,UX,UY, , , , , !施加关键点位移约束FLST,2,3,3,ORDE,2FITEM,2,2FITEM,2,-4/GODK,P51X, ,0, ,0,UY, , , , , ,FLST,2,1,1,ORDE,1FITEM,2,16FLST,2,7,1,ORDE,2FITEM,2,3FITEM,2,-9/GOF,P51X,FY,-20000 !施加节点力FLST,2,2,1,ORDE,2FITEM,2,1FITEM,2,-2/GOF,P51X,FY,-10000FLST,2,1,1,ORDE,1FITEM,2,16/GOF,P51X,FY,-100000/STA TUS,SOLUSOLVE !求解FINISH/POST1PLDISP,1!绘制变形图A VPRIN,0, ,ETABLE, ,SMISC, 6 !单元表数据设置A VPRIN,0, ,ETABLE, ,SMISC, 12A VPRIN,0, ,ETABLE, ,SMISC,2A VPRIN,0, ,ETABLE, ,SMISC, 8A VPRIN,0, ,ETABLE, ,SMISC, 1A VPRIN,0, ,ETABLE, ,SMISC, 7PLLS,SMIS6,SMIS12,-1,0 !绘制内力图PLLS,SMIS2,SMIS8,1,0PLLS,SMIS1,SMIS7,1,0Finish/SOLU !再一次进入求解器ACEL,0,10,0, !施加重力加速度/STA TUS,SOLUSOLVE !求解FINISH/POST1 !进入后处理器PLDISP,1 !绘制变形图ETABLE,REFL !更新PLLS,SMIS6,SMIS12,-1,0 !绘制弯矩图PLLS,SMIS2,SMIS8,1,0 !绘制剪力图SA VEFINISH。
隧道衬砌支护结构的ANSYS数值模拟
隧道衬砌支护结构的ANSYS数值模拟摘要：为了确保隧道施工及运行的安全性，必须对其支护结构进行受力分析。
本文以城市长大隧道为例，基于ANSYS有限元分析软件平台建立隧道支护的荷载—结构模型，并从结构变形、弯矩、轴力和剪力等方面实现对隧道支护结构的数值模拟，从分析结论及安全性的角度出发，为隧道结构的优化设计和现场施工提供依据和指导。
关键词：隧道；支护结构；ANSYS；数值模拟目前，伴随岩土力学的发展和计算机的普遍使用及其性能的不断提高，有限元数值分析已成为隧道结构分析中发展最迅速的方法。
在参数选取合理的情况下，通过对隧道开挖过程进行仿真分析，可判定隧道围岩应力大小以及应力区和塑性区的范围，能够预测隧道施工中的险情，保证隧道施工安全和稳定性。
一、有限元数值模拟方法有限元法的基本思想是将连续的结构离散成有限个单元并设定节点，将连续体看作是只在节点处相连接的一组单元的集合体；同时选定场函数的节点值作为基本未知量，在每一个单元中假设一近似差值函数以表示单元场中场函数的分布规律，利用力学中的某些变分原理去建立用以求解节点未知量的有限元方程，从而将一个连续域中的无限自由度问题化为离散域中的自由度问题，一经求解就可以利用解得的节点值和设定的插值函数确定单元上以至整个集合体上的场函数[1]。
在实际工程应用中，有限元法可以考虑岩土介质的非均匀性、各向异性、非连续性和几何非线性等，适用于各种边界条件，结合大型通用有限元软件ANSYS能较好实现隧道结构的数值计算。
基本建模流程包括选择分析模型类别、创建物理环境、建立模型和划分网格、施加约束和荷载、建立有限元模型、求解和后处理等。
当前，对隧道支护结构体系一般按照荷载—结构模型进行演算，分析过程中将围岩视为隧道结构上的荷载，且为结构本身的一部分，两者间的相互作用通过围岩的弹性支撑对结构施加约束来实现。
二、隧道结构受力分析实例2.1 设计概况目标隧道为双向六车道设计，含多种断面衬砌类型，围岩级别Ⅲ～Ⅵ级。
隧道内力及可靠度分析操作流程手册
隧道结构内力及可靠性分析操作流程手册目录第一章隧道结构内力分析1.1 Ansys菜单操作本例以标准图1202普货-200双线Ⅳ级围岩无砟隧道为例。
本隧道为二衬为c35混凝土，混凝土弹性模量为32E10，泊松比为0.2，Ⅳ级围岩弹性反力系数为350e6，计算按照深埋求得竖向及水平荷载。
第一步：定义工作文件名和工作标题⑴进入ANSYS/Multiphysics的程序界面后，选择菜单Utility Menu：File→Change Jobname，出现Change Jobname对话框。
在【/FILNAM】Enter new Jobname输入框中输入工作名称Support，单击OK按钮关闭该对话框。
⑵选择菜单Utility Menu：File→Change Title命令，出现Change Title对话框，在输入栏中输入Tunnel Support Structural Analysis，单击OK按钮关闭该对话框。
第二步：定义单元类型选择菜单Main Menu：Preprocessor→Element Type→Add/Edit/Delete命令，出现Element Types 对话框，单击Add按钮，出现Library of Element Types对话框。
在左侧滚动栏中选择Structural Beam，在右侧滚动栏中选择2D elestic 3，单击Apply按钮，定义Beam3单元，如图1-1所示。
最后单击Close按钮关闭对话框。
图1-1 单元类型库对话框第三步：定义单元实常数选择菜单Main Menu：Preprocessor→Real Constants→Add/Edit/Delete命令，出现Real Contants对话框，单击Add按钮，出现Element Type for Real Contants对话框，单击OK按钮，选择TYPE 1 BEAM3，单击OK按钮，最后在弹出的Real Constant for BEAM3对话框中分别输入隧道腰部和顶部衬砌支护结构BEAM3梁单元的横截面积AREA：0.45、惯性矩IZZ：0.00759375、高度HEIGHT：0.45，如图1-2所示。
隧道的结构计算公式及强度验算
在十九世纪末，混凝土已经是广泛使用的建筑材料，它具 有整体性好，可以在现场根据需要进行模注等特点。这时，隧 道衬砌结构是作为超静定弹性拱计算的，但仅考虑作用在衬砌 上的围岩压力，而未将围岩的弹性抗力计算在内，忽视了围岩 对衬砌的约束作用。由于把衬砌视为自由变形的弹性结构，因 而，通过计算得到的衬砌结构厚度很大，过于安全。大量的隧 道工程实践表明，衬砌厚度可以减小，所以，后来上述两种计 算方法不在使用了。
图4-3
脚反力的作用下围岩表面将发生弹性变形，使拱脚产生角位移 和线位移。拱脚位移将使拱圈内力发生改变，因而计算中除按 固端无铰拱考虑外，还必须考虑拱脚位移的影响。对于拱脚位 移，还可以作些具体分析，使计算图式得到简化。通常，拱脚 截面剪力很小，它与围岩之间的摩擦力很大，可以认为拱脚没 有径向位移只有切向位移，所以在计算图式中，在固端支座上 用一根径向刚性支承链杆加以约束，见图4-4(a)。切向位移可 以分解为垂直方向和水平方向两个分位移。在结构对称、荷载 对称条件下，两拱脚的位移也是对称的。对称的垂直分位移对 拱圈内力不产生影响。拱脚的转角 和切向位移的水平分位移
图4-4
式中： ——单位变位，即在基本结构上，因
作用时，在Xi方向
上所产生的变位；
——荷载变位，即基本结构因外荷载作用，在Xi方向的变位；
f ——拱圈的变位都能求出，则可用结构力学的力法
知识解算出多余未知力X1和X2，那时，拱圈内力即可求出。
三、拱脚位移计算
1. 单位力矩作用时 单位力矩作用在拱脚围岩上时，拱脚截面绕中心点a转过
一个角度 ，见图4-6，拱脚截面仍保持平面，其内(外)缘处 围岩的最大应力 为：
式中： ——拱脚截面厚度；
b ——拱脚截面纵向单位宽度，
隧道衬砌设计检算
小 导 管 加 强 支 护 。 格 栅 钢 架 纵 向 间 距 计算 结果分 析 为 ０ ８ ，超 前 小 导 管 纵 向 问距 １６ ， ．ｍ ＿ｍ
图 １和 图 ２ 别 为 二 衬 弯 矩 图 和 轴 分
每 环 ２ 根 ，每 根 长 ３ ５ ０ ． ｍ。 隧 道 施 工 均 力 图 ．可 知 ， 弯 矩 较 大 值 出 现 在 拱 顶 按 新 奥 法 采 用 光 面 爆 破 ．锚 喷 支 护 及 湿 喷 技 术 ，并 按 要 求 对 结 构 支 护 体
土 ３５ 。 左 右 方 向 受 水 平 固 定 约 束 ． ｍ
弹 性反力系数 Ｋ（ ＭＰ
ａｍ ） ｌ
弹 性模 量
Ｅ （ ａ） Ｇｐ
泊 松 比
容重
ＫＮ， ｍ
计 算 摩 擦 角
ｃ（ ） 。
土 柱 摩 擦 角 ０
（ ） 。
１０ ５ ２ ８
．
２ ２
５ ５
４ Ｏ
（ 射混凝 土 ） 喷 二 衬 Ｃ３ ０ ｆ 筋混凝 土 ） 钢
Ｋ抗 拉 安 全 系 数 ｌ
１９ ．２
按 单 筋 配 筋Ａ９
构造 配筋
控 制 状 态
受 拉 控 制
Ｂ Ｃ
一 ５５ ６ ． ６ ４
—６ １ ４ ７． —７ ５ １６
Ｏ１ ４ ０ １ ．１
０４ ．３ ０４ ．３
７４ ． ８７ ．
３．４ ４ ４． ７７
；桥 梁 隧 道
｝ ｌ 鹦ｌ － 《 鬻
隧道衬砌设计检算
工程 概况
某 公 路 隧 道 地 址 处 属 丘 陵 地 貌 ， 地 形 起 伏 较 大 相 对 高 差 约
隧道开挖ansys模拟分析
隧道台阶法开挖的有限元模拟分析1.力学模型的建立岩体的性质是十分复杂的，在地下岩体的力学分析中，要全面考虑岩体的所有性质几乎是不可能的。
建立岩体力学模型，是将一些影响岩石性质的次要因素略去，抓住问题的主要矛盾，即着眼于岩体的最主要的性质。
在模型中，简化的岩体性质有强度、变形、还有岩体的连续性、各项同性及均匀性等。
考虑岩石的性质和变形特性，以及外界因素的影响，采用的模型有弹性、塑性、弹塑性、粘弹性、粘弹塑性等。
根据对隧道的现场调查及试验结果分析，围岩具有明显的弹塑性性质。
因此，根据隧道的实际情况，考虑岩体的弹塑性性质，在符合真实施工工序和支护措施的基础上，在数值模拟过程中将计算模型简化成弹塑性平面应变问题，采用Drucker—Prager屈服准则来模拟围岩的非线性并且不考虑其体积膨胀，混凝土材料为线弹性且不计其非线性变形。
对地下工程开挖进行分析，一般有两种计算模型：（1）“先开洞，后加载”在加入初始地应力场前，首先将开挖掉的单元从整体刚度矩阵中删除，然后对剩余的单元加入初始地应力场进行有限元计算。
（2）“先加载，后开洞”这种方法是首先在整个计算区域内作用地应力场，然后在开挖边界上施加反转力，经过有限元计算得到所需要的应力、位移等物理量。
两种方法对线弹性分析而言，所得到的应力场是相同的，而位移场是不同的，模型（2）（即：“先加载，后开洞”）更接近实际情况。
在实际地下工程开挖中部分岩体已进入塑性状态，必须用弹塑性有限元进行计算分析，而塑性变形与加载的路径有关，所以模拟计算必须按真实的施工过程进行，即在对地下工程开挖进行弹塑性数值模拟过程中，必须遵循“先加载，后开洞”的原则。
在有限元法中，求解非线性问题最常采用的方法是常刚度初应力法。
对于弹塑性问题，由于塑性变形不可恢复，应力和应变不再是一一对应的关系，即应力状态与加载路径有关，因此应该用增量法求解。
弹塑性应力增量与应变增量之间的关系可近似地表示为}{}]{[}]){[]([}{][}{0σεεεσd d D d D D d D d p ep +=-== （1） 式中，][D —弹性矩阵，][p D —塑性矩阵。
隧道衬砌结构强度检算
隧道衬砌结构强度检算
课程设计任务书
一、课程设计的主要目的
通过本次课程设计，掌握隧道衬砌结构强度检算的过程，熟悉围岩压力的计算方法，了解荷载－结构模式的计算原理和计算方法，复习巩固钢筋混凝土结构的配筋计算方法，学习了解大型有限元程序ANSYS的建模、求解以及结果处理的过程，为毕业设计以及今后工作奠定基础。
二、主要任务内容
本次课程设计对象为交通隧道，具体设计内容包括：
1、荷载结构模式中围岩压力（荷载）的计算
2、利用ANSYS有限元程序建立模型，并进行求解计算
3、利用有限元计算结果对衬砌结构的强度进行检算，必要时，应进行配筋计算。
三、具体要求
1、课程设计分组进行，每人单独完成相关隧道衬砌结构的计算，隧道图纸由指导教师提供，要求每人至少完成3种断面型式的检算工作。
2、计算过程中所涉及的计算参数应参照现行相关隧道设计规范执行，但要求同一隧道小组成员之间计算参数应不同。
3、设计结束后应提交完整的检算报告，报告封面按网上学校要求的统一格式制作，一律采用A4纸书写或打印，要求页边距为2cm。
4、报告要求明确的计算分析过程和计算结果，所有工作必须独立完成。
5、报告中应附有工程概况、计算对象的横断面图、计算参数取值、计算模型图示、内力图和关键部位内力计算结果汇总表。
6、报告上交时间为2013年3月25日。
指导教师：李文江、孙星亮、韩现民。
隧道衬砌计算
第五章隧道衬砌结构检算5、1结构检算一般规定为了保证隧道衬砌结构的安全,需对衬砌进行检算。
隧道结构应按破损阶段法对构件截面强度进行验算。
结构抗裂有要求时,对混凝土应进行抗裂验算。
5、2 隧道结构计算方法本隧道结构计算采用荷载结构法。
其基本原理为:隧道开挖后地层的作用主要就是对衬砌结构产生荷载,衬砌结构应能安全可靠地承受地层压力等荷载的作用。
计算时先按地层分类法或由实用公式确定地层压力,然后按照弹性地基上结构物的计算方法计算衬砌结构的内力,并进行结构截面设计。
5、3 隧道结构计算模型本隧道衬砌结构验算采用荷载—结构法进行验算,计算软件为ANSYS10、0。
取单位长度(1m)的隧道结构进行分析,建模时进行了如下简化处理或假定:①衬砌结构简化为二维弹性梁单元(beam3),梁的轴线为二次衬砌厚度中线位置。
②围岩的约束采用弹簧单元(COMBIN14),弹簧单元以铰接的方式支撑在衬砌梁单元之间的节点上,该单元不能承受弯矩,只有在受压时承受轴力,受拉时失效。
计算时通过多次迭代,逐步杀死受拉的COMBIN14单元,只保留受压的COMBIN14单元。
图5-1 受拉弹簧单元的迭代处理过程③衬砌结构上的荷载通过等效换算,以竖直与水平集中力的模式直接施加到梁单元节点上。
④衬砌结构自重通过施加加速度来实现,不再单独施加节点力。
⑤衬砌结构材料采用理想线弹性材料。
⑥衬砌结构单元划分长度小于0、5m。
隧道结构计算模型及荷载施加后如图5-2所示。
5、4 结构检算及配筋本隧道主要验算明洞段、Ⅴ级围岩段与Ⅳ级围岩段衬砌结构。
根据隧道规范深、浅埋判定方法可知,Ⅴ级围岩段分为超浅埋段、浅埋段与深埋段。
Ⅳ级围岩段为深埋段。
根据所给的材料基本参数与修改后的程序,得出各工况下的结构变形图、轴力图、建立图与弯矩图。
从得出的结果可知,Ⅴ级围岩深埋段,所受内力均较大,故对此工况进行结构检算。
5、4、1 材料基本参数 (1)Ⅴ级围岩围岩重度318.5/kN m γ=,弹性抗力系数300/k MPa m =,计算摩擦角045ϕ=o ,泊松比u=0、4。
隧道衬砌设计检算
隧道衬砌设计检算
隧道衬砌设计检算是一种专门用于计算隧道衬砌的工作。
它包含对衬砌结构的检查、分析和计算，以确保衬砌的安全性能，确保衬砌的机械特性，以及确保衬砌的维护和维修方面的要求。
隧道衬砌设计检算的基本方法是通过给定衬砌材料的力学参数，根据衬砌结构的外形和尺寸，采用适当的数学方法，来计算和分析衬砌结构的承载能力。
衬砌设计检算的目的一般有2个：第一，确保衬砌结构能够抵御外界荷载；第二，确保衬砌结构的维护和维修要求。
衬砌设计检算的基本步骤主要包括：
1. 设计荷载的确定：衬砌设计的第一步就是确定设计荷载，即需要考虑的最大外部荷载，如：重力荷载、水位荷载、地震荷载等。
2. 基础模型和计算模型的确定：根据设计荷载的确定，再确定基础模型和计算模型，如线性、非线性、稳定性分析等。
3. 衬砌材料的力学参数的确定：根据衬砌材料的性能，确定衬砌材料的力学参数，如抗弯刚度、抗压强度、抗拉强度、屈服强度等。
4. 衬砌结构受力的分析：根据设计荷载和衬砌材料的力学参数，采用适当的计算模型，对衬砌结构受力情况进行分析，以及衬砌结构受力的最大值。
5. 衬砌结构的安全性能分析：根据衬砌结构受力的分析，采用适当的安全系数，确定衬砌结构的安全性能，并确定衬砌结构能够抵御外部荷载的能力。
6. 维护和维修要求的分析：根据衬砌结构的安全性能分析结果，确定衬砌的维护和维修要求，如衬砌的检测要求、衬砌的修理要求等。
隧道衬砌设计检算是隧道施工中必不可少的一个工作，其目的是确保衬砌的安全性能和维护维修要求。
在进行衬砌设计检算时，必须正确确定设计荷载、衬砌材料的力学参数及衬砌结构受力的分析，以便确保衬砌能够抵御外部荷载，保证衬砌的安全性能和维护维修要求。
浅谈ANSYS系统在隧道结构计算中的应用条件
浅谈ANSYS系统在隧道结构计算中的应用条件摘要：在公路隧道设计与施工中，为了提前判断在开挖和支护工程中隧道的结构安全性，隧道结构计算的数值研究方法就成为了一种重要的设计依据和施工控制措施。
本文提供了一种方法，即利用ANSYS软件模拟隧道结构在开挖个步骤中的计算模式与应用条件。
关键词：隧道结构ANSYS模拟隧道的结构分析是利用工程力学原理，选取合理的介质，通过相似模型体系对其结构进行计算，具体过程一般通过两个途径来进行，其一是利用相似性理论，采取合理的相似系数，在室内通过模型试验来模拟实际的工程问题。
其二是数值计算，这种方法伴随着计算机的发展有了长足的进步。
目前，伴随着岩土力学的发展，再加上计算机的普遍使用及其性能的不断提高，有限元法成为发展最迅速的用于隧道结构分析的数值计算方法。
有限元法先将结构分解为有限的小单元，在每一个单元上，利用弹性力学、弹塑性力学等力学理论建立力学性能参数之间的关系，然后根据位移或者应力协调条件把这些小单元组合起来，求出整体结构的力学特征。
因为有限元法是利用矩阵代数方法求解方程组，而矩阵代数建立的方程组非常方便与计算机的存储与求解，所以，有限元法非常适用于分析复杂的地下结构。
1模型的建立利用ANSYS来模拟隧道开挖过程，有两种建模方法，一个是建立真三维的模型，三维模型不仅可考虑围岩的流变特性，还能考虑开挖和支护的空间效应，能保证较好的计算精度。
但是建模复杂，计算时间长，且费用较高。
另一种建模方法是建立二维模型，即按平面应变问题来处理，隧道在长度方向的尺寸比横截面的尺寸大得多,在忽略掘进的空间效应及岩石流变效应的影响时,计算模型取为平面应变是可行的。
另外，可以通过各阶段相应的初始应力释放系数来考虑开挖过程和支护时间早晚对围岩及支护受力的影响。
本文采用后者建立有限元模型。
相对于整个岩体而言，开挖所引起的应力重分布的区域是有限的，因而要确定计算模型的范围。
实践和理论分析表明，对于地下洞室开挖后的应力应变，仅在洞室周围距洞室中心点3～5倍洞室开挖宽度（或高度）的范围内存在实际影响。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、衬砌结构的计算模型
隧道工程建筑物是埋置于地层中的结构物，它的受力和变形与围岩密切相关，支护结构与围岩作为一个统一的受力体系相互约束，共同工作。
这种共同作用正是地下结构与地面结构的主要区别。
根据本工程浅埋及松散地层的特点，使用阶段结构安全性检算采用“荷载—结构”模式，即将支护和围岩分开考虑，支护结构是承载主体，围岩作为荷载的来源和支护结构的弹性支承。
支护结构与围岩的相互作用是通过弹性支承对支护结构施加约束来实现的。
计算模型中，二衬结构采用弹性平面梁单元模拟，弹性抗力以及隧底地基均采用弹簧单元模拟。
组合荷载根据不同作用方向分别转换成等效节点力施加在相应的单元结点上。
具体计算模型见图1。
图1 计算模型
二、荷载计算
围岩压力计算参照课本中有关我国铁路隧道推荐的方法进行确定（双线隧道）或参照《铁路隧道设计规范》，深浅埋分别计算。
三、配筋计算
结构强度检算和配筋计算应按照现行《铁路隧道设计规范》的方法进行。
四、ANSYS操作过程
1、更改路径和工作名
2、进入前处理模块（preprocessor）
(1)定义单元类型element type
(2)定义实常数real constant
(3)定义材料参数material props
(4)定义梁的截面特性sections
(5)进入modeling进行建模，生成几何模型
(6)进行网格划分meshing
a)给几何模型赋属性
meshing>mesh attribute>picked lines
(7)施加弹性约束
Model>Creat>piping models>spring support
定义弹性抗力系数K和距离所选结构节点的相对距离DX, DY, DZ。
3、进入求解器solution
(1)定义分析类型analysis type>new analysis>static
(2)定义荷载define loads
(3)设置荷载添加形式setting>replace vs add>force，按如下图示设置
(4)施加等效节点力
define loads>apply>force(编程实现)
(5)施加重力
define loads>apply>inertia>gravity>global
(6)求解计算
Solve>current LS
4、后处理（general postproc）
(1)读入结果
Read results>last set
(2)查看变形图，（plot results>deformed shape）检查弹簧约束范围是否正确（所有弹簧均应受压，即处于抗力区）否则添加或删除弹簧单元，重新计算。
(3)定义单元表
Element table>define table，出现对话框后点击add按钮，出现下列对话框：
选择by sequence num 并在右下角SMISC,后分别键入数字“6”后点击apply按钮。
重复上述操作并分别键入数字“12”和“1”和，点击OK，完成弯矩和轴力的定义。
(4)查看内力图
plot results>contour plot>line elem res，查看弯矩云图
查看轴力图
(5)存储显示图形
Plotctrls>capture image,将文件另存即可。
(6)列表显示计算结果
List results>elem table data
出现对话框，选择单元表定义的内容，点击ok。
然后将文件另存即可。
附注：加载程序示意
1.均布荷载
(1)y方向荷载施加
*do, i, 1, 10
a=q*abs(nx(i)-nx(i+1))/2
f, i, fy, a
f, i+1, fy, a
*enddo
(2)x方向荷载施加
*do, i, 1, 10
a=q*abs(ny(i)-ny(i+1))/2
f, i, fx, a
f, i+1, fx, a
*enddo
2.梯形荷载
*set,shuxiang,292608.9
*set,shangce,78149.5
*set,xiace,127057.6
!竖向荷载
*do,i,5,18
a=-shuxiang*0.5*abs(nx(i)-nx(i+1))
f,i+1,fy,a
f,i,fy,a
*enddo
!左水平荷载(11.82为梯形的高度)
！3指梯形上底对应的节点号
！12指梯形下底对应的节点号
*do,i,3,11
a=abs(ny(i)-ny(12))*(xiace-shangce)/11.82+shangce
b=abs(ny(i+1)-ny(12))*(xiace-shangce)/11.82+shangce f,i+1,fx,(a/6+b/3)*abs(ny(i)-ny(i+1))
f,i,fx,(a/3+b/6)*abs(ny(i)-ny(i+1))
*enddo
*do,i,22,26
a=abs(ny(i)-ny(12))*(xiace-shangce)/11.82+shangce
b=abs(ny(i+1)-ny(12))*(xiace-shangce)/11.82+shangce f,i+1,fx,(a/6+b/3)*abs(ny(i)-ny(i+1))
f,i,fx,(a/3+b/6)*abs(ny(i)-ny(i+1))
*enddo
a=abs(ny(1)-ny(12))*(xiace-shangce)/11.82+shangce b=abs(ny(3)-ny(12))*(xiace-shangce)/11.82+shangce f,3,fx,(a/6+b/3)*abs(ny(i)-ny(i+1))
f,1,fx,(a/3+b/6)*abs(ny(i)-ny(i+1))
a=abs(ny(1)-ny(12))*(xiace-shangce)/11.82+shangce b=abs(ny(22)-ny(12))*(xiace-shangce)/11.82+shangce f,22,fx,(a/6+b/3)*abs(ny(i)-ny(i+1))
f,1,fx,(a/3+b/6)*abs(ny(i)-ny(i+1))
!右水平荷载
*do,i,12,20
a=abs(ny(i)-ny(12))*(xiace-shangce)/11.82+shangce
b=abs(ny(i+1)-ny(12))*(xiace-shangce)/11.82+shangce f,i+1,fx,-(a/6+b/3)*abs(ny(i)-ny(i+1))
f,i,fx,-(a/3+b/6)*abs(ny(i)-ny(i+1))
*enddo
*do,i,27,31
a=abs(ny(i)-ny(12))*(xiace-shangce)/11.82+shangce
b=abs(ny(i+1)-ny(12))*(xiace-shangce)/11.82+shangce f,i+1,fx,-(a/6+b/3)*abs(ny(i)-ny(i+1))
f,i,fx,-(a/3+b/6)*abs(ny(i)-ny(i+1))
*enddo
a=abs(ny(2)-ny(12))*(xiace-shangce)/11.82+shangce b=abs(ny(21)-ny(12))*(xiace-shangce)/11.82+shangce f,21,fx,(a/6+b/3)*abs(ny(i)-ny(i+1))
f,2,fx,(a/3+b/6)*abs(ny(i)-ny(i+1))
a=abs(ny(2)-ny(12))*(xiace-shangce)/11.82+shangce b=abs(ny(32)-ny(12))*(xiace-shangce)/11.82+shangce
f,32,fx,(a/6+b/3)*abs(ny(i)-ny(i+1))
f,2,fx,(a/3+b/6)*abs(ny(i)-ny(i+1))
上述程序中，i为节点号，q为竖向围岩压力，a为等效节点力计算中间值，nx(i)为i节点的水平坐标值，竖直坐标用ny(i)表示。