盾构隧道结构ansys计算方法
ANSYS计算算例(结构-荷载法)
（3）直边墙的输入
第二段圆弧的终点即是直边墙的顶点，查出该点坐标为（-3.520017，-3.077515）， 并由此计算出直边墙底端点的坐标为（-3.520017，-7.507515）
GUI：Main Menu> Preprocessor> Create>Keypoints>In Active CS
计算中考虑结构自重对结构内力的影响，还应在此输入材料的密度。 Nonlinear-Density
4、模型建立 （1）第一段圆弧的输入（圆心：0，－3.41，起点： 0，0，圆心角：45º ）
Menu menu：Preprocessor>Material Props>Modeling Create>Lines Arcs>By Cent & Radius
三、加载及约束
1、荷载计算及转化 (1) 荷载计算 对深埋隧道用《隧规》推荐的统计法公式计算得： 竖向围岩压力：q=54084.25Pa； 水平围岩压力：e＝0.15 q=8112.64Pa (2) 匀布荷载转化为等效节点荷载 将匀布荷载按如下步骤转化成等效节点力。 A、列出结构上所有节点坐标 GUI：Utility Menu>List>Nodes
拾取基底两节点
选取UX,UY并输入限制位移值0
3、惯性荷载加入 惯性荷载是与质量相关的荷载，该处是加入重力加速度。为了利用惯 性效果来模拟重力，应当在重力的相反方向施加惯性项，即在Y方向 施加一个加速度荷载相当于模拟-Y方向的重力： GUI：Menu Main> Solution>Apply>Gravity
ANSYS计算算例
• 拱形结构算例 – 铁路隧道 – 公路隧道 – 地铁区间隧道(矿山法施工) 矩形框架结构 – 地下街 – 地下停车场 – 地铁车站 圆形结构 – 地铁区间隧道(盾构法施工)
盾构法隧道等代层参数反演的ANSYS方法
张宏洲 ， 张金伟 ，翟建华
（． １ 中国地质大学工程技术学院， 北 ０ 河南 洛 阳 ４ １０ ） ７ ０９
摘要 ： 在盾 构法隧道 中将衬砌周 围 由于盾 尾空隙的闭合 、 注浆充填作用和盾构推进对周 围地层 的扰动 等形成 的特殊 圈层 概化 为具
中图分类 号 ： ４ １ Ｕ ５ 文献标识码 ： Ａ
ＡＮＳ ｅ ｈ ｄ ｆ ｒ Ｄｅ ｕ ｔｏ ｆＰａ ａ ｅ ｅ ｓ ＹＳ Ｍ ｔ ｏ ｏ ｄ ｃ ｉ ｎ ｏ ｒ ｍ ｔ ｒ
ｏ ｕ ｖ ｌｎｔＣｉｃ ｌ ｒ Ｚｏ ｆＳｈ ｅｄ Ｔｕ ｎｅ ｆＥｑ ｉａ ｅ ｒ ｕ ａ ｎｅ ｏ ｉｌ ｎ ｌ
在地层位移分析中 ， 对上述过渡层 的处理都是 问 题的关键 ， 也是问题 的难点。前人对这一问题 的处理 包括 ：） 为注浆材料 完全充填盾尾 空隙 ， １认 隧道壁面 土体完全没有 向盾尾 空隙的移动。２ 完 全不考 虑注 ）
间、 为衬砌的拼装操作所 留的建筑空隙和盾构推进时 部分土体被粘附于盾构外壳上 ， 在衬砌环背面与实际
ｈ ｏ ｃ ｐ ｑ ｉａｅ ｔｃｒ ｌｒｚ ｎ ｕ ｐ ｓ ｄ ｔ ｅ ｃ ｎ ｅ ｔ ｆｅ ｕｖ ｌｎ ｉ ｕ ａ ｏ ｅ ｉ ｓ ｐ ｏ ｅ ．Ｗ ｉ ｈ ｅｐ ｏ ｌ ａ ｕ ｅ ｉｐ ａ ｅ ｎ ｓ ｈ ａ ａ ｔｒ ｆｔ ｅ ｏ ｃ ｓ ｔ ｔ ｅ ｈ ｌ ｆｆ ｅ ｍｅ ｓ ｒ ｄ ｄ ｓ ｌ ｃ ｍｅ ｔ ，ｔ ｅ ｐ ｒ ｍｅｅ ｓ ｏ ｈ ｉ ｈ ｅ ｕ ｖ ｌｎ ｉ ｕ ａ ｏ ｅ Ｃ ｅｏ ｔｉ ｅ ｘ ｃｌ ．Ｔ ｉ ｐ ｐ ｒｕ ｅ ｑ ｉ ａｅ ｔ ｒ ｌｒｚ ｎ ａ ｂ ｂ ａｎ ｄ ｅ ａ ｔ ｃ ｃ ｎ ｙ ｈｓ ａ ｅ ｓ ｓＡＮＳ ｎ ｔ ｌ ｍｅ ｔ ｔ ｏ ｏ ｄ ｄ ｃ ｅ ｐ ｒ ｍｅ ｅｓ ＹＳ ｆ ｉ ｅｅ ｎ ｈ ｄ ｔ ｅ ｕ ｔ ａ ａ ｔｒ ｉ ｅ ｍｅ ｈ ｔ ｏ ｅ ｅ ｕ ｖ ｌｎ ｉ ｕ ａ ｏ ｅ ｆ ｑ ｉ ａｅ ｔ ｒ ｌ ｚ ｎ ，ｗｈ ｃ ａ ｒ ｖ ｄ ｅ ｅ ｅ ｃ ｏ ｅ ｄ ｓｇ ｉ ｇ ａ ｄ ｃ ｎ ｔ ｃｉ ｎ ｏ ｈ ｅｄ ｔ ｎ ｅ ｓ ｈ ｔ ｃ ｃ ｒ ｉｈ ｃ ｎ ｐ ｏ ｉ ｅ ｒｆ ｒ ｎ ｅ ｆｒｔ ｅ ｉ ｎ ｎ ｎ ｏ ｓｒ ｔ ｆｓ ｉｌ ｕ ｎ ｌ． ｈ ｕ ｏ Ｋｅ ｒ ｓ ｈ ｅｄ ｔ ｎ ｅ ：ｅ ｕｖ ｌｎ ｉ ｕ ａ ｏ ｅ ｎ ｅ ｓ ｎ ｙ ｉ ；ＡＮＳ ｙ ｗｏ ｄ ：ｓ ｉｌ ｕ ｎ ｌ ｑ ｉａｅ ｔｃｒ ｌｒｚ ｎ ：ｉ ｖ ｒ ｅ ａ ａ ｓ ｃ ｌ ｓ ＹＳ
地铁盾构隧道掘进过程ANSYS数值模拟分析
地铁盾构隧道掘进过程数值模拟分析具体做法如下：⑴采用在掘削面施加顶进压力的方法来模拟开挖面土体的移动；⑵采用施加注浆压力的方法来模拟盾尾注浆过程；⑶采用更换注浆层材料参数的方法来模拟盾尾注浆对盾尾空隙的填充效果。
由于目前计算软件的限制，难以模拟盾构机推进过程中对土体的扰动，这里简化处理。
即假定盾构隧道开挖后，随机进行注浆。
计算时，只需将开挖不断地向前推进，同时在后面进行注浆、换材料参数等操作，即可实现盾构隧道的动态开挖过程，详细的计算操作见后面的求解过程。
工程问题的描述地铁盾构隧道管片衬砌内径为 5.4m，外径为D=6m，埋深为12m 自上至下，根据土层的物理性质参数不同将其分为 3 层，各层的材料参数和厚度如下：第一层：厚8m，E=3.94Mpa，v=0.35，ρ=18.28KN/m3第二层：厚18m，E=20.6Mpa，v=0.3，ρ=20.62KN/m3第三层：厚15m，E=500Mpa，v=0.33，ρ=21.6KN/m3 施工中掘削面顶进压力为0.3Mpa，盾尾注浆压力为0.15Mpa 模型的建立!进入前处理器FINISH/CLE/PREP7 !进入前处理器ET,1,SOLID45 !定义实体单元ET,2,MESH200,6 !定义非求解单元，辅助面网格的划分! 定义模型中的材料参数。
模型中共有 5 种材料，其中土体有 3 种，即地表浅层覆土、盾构隧道所在土层和基岩及管片衬砌和注浆层。
其中，管片衬砌为管片式的拼装结构，为了计算方便，将其等效为一均质体，等效时对原有刚度进行折减。
定义材料参数的命令流如下：!土体材料参数MP,EX,1,3.94E6 !第一层土层材料参数MP,PRXY,1,0.35MP,DENS,1,1828MP,EX,2,20.6E6 !第二层土层材料参数MP,PRXY,2,0.30MP,DENS,2,2160MP,EX,3,500E6 !第三层土层材料参数MP,PRXY,3,0.33MP,DENS,3,2160!管片材料参数，管片衬砌按各向同性计算MP,EX,4,27.6E9MP,PRXY,4,0.2MP,DENS,4,2500!注浆层，参数按水泥土取值MP,EX,5,1E9MP,PRXY,5,0.2MP,DENS,5,2100!建立平面内模型并划分单元!在隧道中心线定义局部坐标，便于后来的实体选取LOCAL,11,0,0,0,0 ! 设置局部直角坐标系原点坐标(0,0,0)LOCAL,12,1,0,0,0 ! 设置局部柱坐标系原点坐标(0,0,0) CSYS,11WPCSYS,-1CYL4,,,,,2.7,90CYL4,0,0,2.7,0,3,9 0CYL4,0,0,3,0,3.2,9 0RECTNG,0,4.5,0,4.5 AOVLAP,ALL NUMMRG ,ALL !创建开挖土体所在的面!创建管片!压缩编号!将当前坐标系转化为局部直角RECTNG,4.5,31.5,0,4.5NUMMRG ,ALLNUMCMP,ALL!划分单元MSHAPE,0,2D !采用四边形单元划分网格MESHKEY ,1 !映射网格划分TYPE,2LESIZE,1,,,6LESIZE,2,,,6LESIZE,3,,,6AMESH,1LESIZE,4,,,6LESIZE,8,,,2LESIZE,9,,,2AMESH,2LESIZE,5,,,6LESIZE,10,,,1LESIZE,11,,,1AMESH,3LESIZE,12,,,3LESIZE,13,,,3LESIZE,6,,,3LESIZE,7,,,3LESIZE,14,,,8,2LESIZE,16,,,8,0.5AMAP,4,7,6,8,10AMAP,5,9,8,11,12!利用对称性，得到下半部分模型ARSYM,Y ,ALL !通过坐标轴对称建立面NUMMRG ,ALLNUMCMP,ALLALLSEL,ALL!建立隧道下方土层模型RECTNG,0,4.5,-4.5,-26RECTNG,4.5,31.5,-4.5,-26NUMMRG ,ALLNUMCMP,ALLLESIZE,28,,,3LESIZE,29,,,5,0.5LESIZE,30,,,5,2LESIZE,32,,,5,0.5LESIZE,31,,,8,2ASEL,S,AREA,,11,12,1!建立隧道土方土层模型AMESH,ALLRECTNG,0,4.5,4.5,15RECTNG,4.5,31.5,4.5,15NUMMRG ,ALLNUMCMP,ALLLESIZE,34,,,3LSEL,S,LINE,,35,36,1LSEL,A,LINE,,33LESIZE,ALL,,,4LSEL,S,LINE,,37LESIZE,37,,,8,0.5AMESH,13AMESH,14NUMMRG ,ALLNUMCMP,ALL!利用对称性得到平面内的全部模型ALLSEL,ALLARSYM,X,ALLNUMMRG ,ALLNUMCMP,ALL!建立体模型。
隧道ansys计算程序算例_荷载结构模式
选取新建铁路宜昌（宜）－万州（万）铁路线上的别岩槽隧道某断面.该断面设计单位采用的支护结构如图3-3所示。
为保证结构的安全性.采用了荷载—结构模型.利用ANSYS 对其进行计算分析。
主要参数如下：●隧道腰部和顶部衬砌厚度是65cm.隧道仰拱衬砌厚度为85cm。
●采用C30钢筋混凝土为衬砌材料。
●隧道围岩是Ⅳ级,洞跨是5.36米.深埋隧道。
●隧道仰拱下承受水压.水压0.2MPa。
图 3-3 隧道支护结构断面图隧道围岩级别是Ⅳ级.其物理力学指标及衬砌材料C30钢筋混凝土的物理力学指标见表3-3所示。
根据《铁路隧道设计规范》.可计算出深埋隧道围岩的垂直匀布力和水平匀布力。
对于竖向和水平的分布荷载.其等效节点力分别近似的取节点两相临单元水平或垂直投影长度的一般衬砌计算宽度这一面积范围内的分布荷载的总和。
自重荷载通过ANSYS程序直接添加密度施加。
隧道仰拱部受到的水压0.2MPa按照径向方向载置换为等效节点力,分解为水平竖直方向加载。
3.3.3 GUI操作方法3.3.3.1 创建物理环境1) 在【开始】菜单中依次选取【所有程序】/【ANSYS10.0】/【ANSYS Product Launcher】.得到“10.0ANSYS Product Launcher”对话框。
2）选中【File Management】.在“Working Directory”栏输入工作目录“D:\ansys\example301”.在“Job Name”栏输入文件名“Support”。
3）单击“RUN”按钮.进入ANSYS10.0的GUI操作界面。
4）过滤图形界面：Main Menu> Preferences.弹出“Preferences for GUI Filtering”对话框.选中“Structural”来对后面的分析进行菜单及相应的图形界面过滤。
5）定义工作标题：Utility Menu> File> Change Title.在弹出的对话框中输入“Tunnel Support Structural Analysis”.单击“OK”.如图3-4所示。
盾构隧道衬砌圆环内力计算实例(含命令流)
截 面 惯 性 矩 ： IZZ=1*0.5*0.5*0.5/12 截 面 高 度 ： HEIGHT=0.5， 单 击 OK， 则 1 号 实 常 数 定 义 完 成 。 11 定 义 衬 砌 材 料 属 性 ○ 执 行 Main Menu>Preprocessor>Material Props>Material Models 命 令 ， 弹 出 定 义 材 料 模 型 行 为 对 话 框 。选 择 Structural>Linear>Elastic>Isotropic, 在 弹 出 的 对 话 框 中 输 入 弹 性 模 量 EX=36e9,泊 松 比 PRXY=0.33，单 击 OK。再 选 择 Structural>Density,在 弹 出 的 对 话 框 中 输 入 密 度 DENS=2600。 12 ○ 单元网格划分
4 ○ 5 ○
进 入 前 处 理 器 。 单 击 Main Menu>Preprocessor。 画 点 。 执 行 Main Menu>Preprocessor>Modeling>Create>Kdypoint>In
Active CS 命 令 ， 设 置 关 键 点 1 作 为 圆 心 ， 输 入 坐 标 （ 0， 0） 。
图 5.2.6 节 点 编 号 图
14 ○
加载
本例中，衬砌所承受的荷载分为主动压力与被动抗力。盾构法施工中，开 挖 与 架 设 衬 砌 几 乎 同 时 进 行 ，所 以 荷 载 来 不 及 释 放 ，全 部 作 用 在 衬 砌 结 构 上 形 成 主 动 压 力 ，这 在 荷 载 计 算 小 节 里 面 已 经 算 出 。在 主 动 荷 载 作 用 下 ，衬 砌 结 构 的一部分将发生向着围岩方向的变形， 具有一定刚度的围岩就会对衬砌结构 产 生 反 作 用 力 来 抵 制 它 的 变 形 ，即 为 被 动 抗 力 。这 种 被 动 抗 力 的 施 加 通 过 弹 簧 支撑来体现的。 a． 计 算 主 动 压 力 产 生 的 等 效 节 点 荷 载 。
盾构隧道荷载计算方法
盾构隧道荷载计算方法我折腾了好久盾构隧道荷载计算方法，总算找到点门道。
咱先说说土压力这一块的荷载计算吧。
一开始我真是瞎摸索，就知道土压力肯定和土的重量啊，隧道的深度之类的有关系。
我做过这么个尝试，就按照土力学里最基础的公式来算静止土压力，也就是认为土层是静止不动的这么个状态去计算压力。
我就把土的重度，隧道埋深那些数值往里套，感觉挺简单的是不？但实际上，这可差点让我吃了大亏。
后来我才发现这种计算方法太理想化了，像在软土地层里就非常不准确。
再来说说盾构推进过程中的附加荷载。
这个就很复杂啦。
你可以想象盾构机就像个在土里慢慢往前挤的大虫子。
我试过很多方法去估算这个附加荷载。
我有一次按照平均扩散角的方法去计算附加荷载，就是假设盾构机给土层的力量是沿着一个角度扩散开的，就像是把一个东西放在沙堆上，周围的沙子呈现一定角度散开那样。
可是我在考虑土层性质的时候就搞混了，不同的土层刚度不一样，这个扩散角肯定也不一样啊，我当时就简单地用了同一个数值，结果算得那叫一个离谱。
后来我就慢慢总结出来，在计算盾构隧道荷载的时候，一定要搞清楚地层的情况。
比如地层是软土还是岩石土，这个区别可大了。
如果是软土，很可能需要考虑土体受到盾构机挤压后的变形，这个变形又会反过来影响荷载。
如果是岩石土，那可能更多要考虑岩石的节理裂隙那些状况。
关于计算水压力的荷载，这也是我走过弯路的地方。
有水的情况下可不能简单的就用静水压力的公式。
像在有渗透的情况下，水压力就不是单纯的只和水深有关系了。
我之前没考虑到这个，结果就完全不对。
就是说在实际工程中，地下水的流动啊，排水效果啊这些都会影响水压力荷载。
还有就是建筑物附加荷载。
要是盾构隧道周边有建筑物，那就得考虑建筑物的自重给地层增加的压力，以及这个压力会怎么传播到盾构隧道上。
我只记得大概的方法，就是要根据建筑物基础的类型，比如是桩基础还是筏板基础，再根据它们和隧道的相对位置啊之类的去计算。
但是这里面具体的系数和计算方法在不同规范或者实际工程里还可能有变化，这个我还不是特别确定，不过我觉得多看实际工程案例的计算书肯定能搞明白。
隧道及地下工程ANSYS实例分析_第5章_双线铁路隧道衬砌结构设计力学分析[1]
![隧道及地下工程ANSYS实例分析_第5章_双线铁路隧道衬砌结构设计力学分析[1]](https://img.taocdn.com/s3/m/85e80cf87c1cfad6195fa70a.png)
5.1二次衬砌结构力学分析/TITLE，Mechanical analysis on railway tunnel 2nd lining ! 确定分析标题/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 !设置梁单元类型ET，2，COMBIN14 !设置弹簧单元类型R，1，0.4，0.0053333，0.4，，，，!设置梁单元几何常数R，2，400e6，，，!设置弹簧单元几何常数MPTEMP，，，，，，，，!设置材料模型MPTEMP，1，0MPDATA，EX，1，，30.0e9 !输入弹性模量MPDATA，PRXY，1，，0.2 !输入泊松比MPTEMP，，，，，，，，!设置材料模型MPTEMP，1，0MPDATA，DENS，1，，2500 !输入密度SAVE !保存数据库1．建立几何模型K, 1, , , , !创建关键点（隧道二次衬砌）K, 2, 4.71, -1.82, ,K, 3, -4.71, -1.82, ,K, 4, 0, -3.75, ,K, 5, 0, 5.05, ,K, 20, 5.6, -2.162, , !创建关键点（地层弹簧）K, 30, -5.6, -2.16, ,K, 40, 0, -4.75, ,K, 50, 0, 6, ,!创建隧道衬砌线LARC, 2, 5, 1, 5.05, !创建圆弧线（拱顶部）LARC, 5, 3, 1, 5.05,LARC, 2, 3, 4 !创建圆弧线（仰拱部）第1章大型有限元软件ANSYS简介2 ! 创建地层弹簧线LARC, 20, 50, 1, 6, !创建圆弧线（拱顶部）LARC,50, 30, 1, 6,LARC, 20, 30, 40 !创建圆弧线（仰拱部）SAVE !保存数据2．单元网格划分设置单元大小并将所有直线划分单元，其单元图如图5-8所示。
同济大学Ansys隧道荷载结构法算例
隧道荷载结构法算例 1． 问题描述：
主要截面尺寸如图 1 所示， R 1 5.2m ， R2 8.8m 。内轮廓采用曲墙拱三心圆拱。弹性抗力 系数取为 k 1.5 10 kN / m 。衬砌厚度为 0.6m，截面的物理参数如表所示。
建立局部柱坐标系，施加径向弹簧（用 PSPRNG 命令） 。根据计算结果中弹簧是否收拉，对结 果进行迭代计算。第 1 次计算结果如下图所示。
图5
结构变形图
删除受拉弹簧后的计算结果如下图。
图6 命令流： /prep7 et,1,beam3 mp,ex,1,2.85e10 mp,prxy,1,0.25 r,1,0.6,0.018,0.6 k,1,0,5.2 k,2,-5.2,0 k,3,-5.0797,-1.45 k,4,0,0 k,5,3.6,0 larc,1,2,4,5.2 larc,2,3,5,8.8 lesize,1,,,8 lesize,2,,,2 lmesh,all /pnum,node,1 nplot nlist !compute nodal force in excel f,1,fx,1.249456357 f,3,fx,4.949809522 f,4,fx,9.709400641
finish /solu solve finish /post1 pldisp,1 !check springs with traction (+ means traction) etable,sforce,smisc,1 pretab,sforce !delete springs with tractions and recalculate here ! etable,m1,smisc,6 etable,m2,smisc,12 plls,m1,m2,-1
盾构法隧道等代层参数反演的ansys方法
盾构法隧道等代层参数反演的ansys方法
隧道工程中，代层参数是对地质地形的描述，是隧道地质施工的重要依据。
在盾构施工中，代层参数可以帮助工程师确定掘进工作面的工作条件和防止隧道坍塌等安全问题。
因此，正确地测量和分析代层参数非常重要。
ANSYS是一款强大的工程仿真软件，可以帮助工程师模拟隧道代层参数的反演过程。
通过ANSYS，工程师可以建立一个三维模型，模拟隧道的设计和施工过程，并掌握每个代层参数的变化。
首先，工程师需要利用现有数据建立一个三维的地质模型。
这个模型包括地层、构造、岩性、孔隙率等信息，能够模拟隧道施工地理环境的所有要素。
在这个模型基础上，可进行力学分析、有限元分析等操作，定量测量只能测量特定点位的所有参数。
然后，结合测量数据，利用ANSYS进行数值分析和优化。
ANSYS 能够根据已有的数据，反演出每个代层参数的数值，并根据实际的建构过程，对数值进行反复地调整和优化，以适应现场环境的变化。
最后，工程师可基于ANSYS的分析结果，制定相应的代层参数处理方案，以便更好地指导实际的隧道施工。
ANSYS方法为隧道施工的精细化管理提供了重要支持，可提高施工质量、降低风险，为隧道工程的安全高效完成奠定了坚实的基础。
浅谈ANSYS系统在隧道结构计算中的应用条件
浅谈ANSYS系统在隧道结构计算中的应用条件摘要：在公路隧道设计与施工中，为了提前判断在开挖和支护工程中隧道的结构安全性，隧道结构计算的数值研究方法就成为了一种重要的设计依据和施工控制措施。
本文提供了一种方法，即利用ANSYS软件模拟隧道结构在开挖个步骤中的计算模式与应用条件。
关键词：隧道结构ANSYS模拟隧道的结构分析是利用工程力学原理，选取合理的介质，通过相似模型体系对其结构进行计算，具体过程一般通过两个途径来进行，其一是利用相似性理论，采取合理的相似系数，在室内通过模型试验来模拟实际的工程问题。
其二是数值计算，这种方法伴随着计算机的发展有了长足的进步。
目前，伴随着岩土力学的发展，再加上计算机的普遍使用及其性能的不断提高，有限元法成为发展最迅速的用于隧道结构分析的数值计算方法。
有限元法先将结构分解为有限的小单元，在每一个单元上，利用弹性力学、弹塑性力学等力学理论建立力学性能参数之间的关系，然后根据位移或者应力协调条件把这些小单元组合起来，求出整体结构的力学特征。
因为有限元法是利用矩阵代数方法求解方程组，而矩阵代数建立的方程组非常方便与计算机的存储与求解，所以，有限元法非常适用于分析复杂的地下结构。
1模型的建立利用ANSYS来模拟隧道开挖过程，有两种建模方法，一个是建立真三维的模型，三维模型不仅可考虑围岩的流变特性，还能考虑开挖和支护的空间效应，能保证较好的计算精度。
但是建模复杂，计算时间长，且费用较高。
另一种建模方法是建立二维模型，即按平面应变问题来处理，隧道在长度方向的尺寸比横截面的尺寸大得多,在忽略掘进的空间效应及岩石流变效应的影响时,计算模型取为平面应变是可行的。
另外，可以通过各阶段相应的初始应力释放系数来考虑开挖过程和支护时间早晚对围岩及支护受力的影响。
本文采用后者建立有限元模型。
相对于整个岩体而言，开挖所引起的应力重分布的区域是有限的，因而要确定计算模型的范围。
实践和理论分析表明，对于地下洞室开挖后的应力应变，仅在洞室周围距洞室中心点3～5倍洞室开挖宽度（或高度）的范围内存在实际影响。
盾构管片修正惯用法内力计算 ansys命令流
!匀质圆环,在原程序上加了静水压力，故本命令流适合水土分算fini/cle*dim,wxn,array,2000 !!定义名为wxn的数组，2000行1列，下同*dim,wyn,array,2000*dim,xn,array,2000*dim,yn,array,2000*dim,fxn,array,2000*dim,fyn,array,2000*dim,aa1,array,2000*dim,bb1,array,2000*dim,aa2,array,2000*dim,bb2,array,2000*dim,jx1,array,2000*dim,jy1,array,2000*dim,jx2,array,2000*dim,jy2,array,2000定义参数并赋值nodesum=120 单元数pi=3.1415927height=0.35 环厚0.35米width=1.0 环宽area=height*width 面积inertia=width*height*height*height/12 惯性矩emod=34500*1000000 弹性模量，按C50混凝土计radius=2.925 计算半径density=25*1000/10 材料密度wradius=radius+0.5 弹簧节点所处半径tankang=12*1000000*width 弹簧刚度,MPa/m赋荷载值ptop=200*1000*width 顶板压力，单位按KN计pltop=80*1000*width 侧向压力plbot=50*1000*width 侧向沿竖向增加的侧压，即墙底水平压力为pltop+plbot pbot=140*1000*width 底板压力water=8.0*1000*width 顶板处水位埋深/prep7et,1,beam3 定义单元类型,弹模,密度，实常数mp,ex,1,emodmp,dens,1,densitymp,prxy,1,0.2r,1,area,inertia,heightet,2,combin39 定义弹簧单元keyopt,2,1,0 0为沿与加载相同的路径卸载keyopt,2,2,0 按照所定义的受压曲线加压力荷载keyopt,2,3,0 沿节点X轴方向平动keyopt,2,4,3 二维轴向单元keyopt,2,6,0 基本单元输出*do,i,1,nodesum,1xn(i)=radius*cos(pi/2+(i-1)*(2*pi/nodesum)) 定义管片节点I的x坐标yn(i)=radius*sin(pi/2+(i-1)*(2*pi/nodesum)) 定义管片节点I的y坐标n,i,xn(i),yn(i) 生成管片节点I*enddo*do,i,1,nodesum,1wxn(i)=wradius*cos(pi/2+(i-1)*(2*pi/nodesum)) 定义弹簧节点I的x坐标wyn(i)=wradius*sin(pi/2+(i-1)*(2*pi/nodesum)) 定义弹簧节点I的y坐标n,i+nodesum,wxn(i),wyn(i) 生成弹簧节点I*enddo生成单元，类型为1，实常数为1type,1real,1*do,i,1,nodesum-1,1e,i,i+1*enddoe,nodesum,1生成弹簧单元len=radius*(2*pi/nodesum)wytankang=len*tankangr,2,-1,-wytankang,0,0,1,wytankang/1000000type,2real,2*do,i,1,nodesum,1e,i,i+nodesum*enddo*do,i,2,nodesum,1aa1(i)=abs(xn(i)-xn(i-1)) 取两节点X坐标差的绝对值，并赋给数组aa1bb1(i)=abs(yn(i)-yn(i-1))jx1(i)=xn(i)*xn(i-1) 取两节点X坐标的乘积，并赋给数组jx1 jy1(i)=yn(i)*yn(i-1)*enddoaa1(1)=abs(xn(1)-xn(nodesum))bb1(1)=abs(yn(1)-yn(nodesum))jx1(1)=xn(1)*xn(nodesum)jy1(1)=yn(1)*yn(nodesum)*do,i,1,nodesum-1,1aa2(i)=abs(xn(i)-xn(i+1)) 相邻两节点坐标差的绝对值bb2(i)=abs(yn(i)-yn(i+1))jx2(i)=xn(i)*xn(i+1) 相邻两节点x坐标乘积jy2(i)=yn(i)*yn(i+1)*enddoaa2(nodesum)=abs(xn(nodesum)-xn(1))bb2(nodesum)=abs(yn(nodesum)-yn(1))jx2(nodesum)=xn(nodesum)*xn(1)jy2(nodesum)=yn(nodesum)*yn(1)*do,i,1,nodesum,1*if,yn(i),lt,0,then 如yn(i)小于0（x轴以下的节点），则a1=aa1(i)a2=aa2(i)fyn(i)=pbot*(a1+a2)/2*elseif,yn(i),gt,0 如yn(i)大于0（x轴以上的节点），则a1=aa1(i)a2=aa2(i)fyn(i)=-ptop*(a1+a2)/2*elsefyn(i)=0*endif*enddo*do,i,1,nodesum,1*if,xn(i),lt,0,then 如xn(i)小于0(y轴左侧节点)，则b1=bb1(i)b2=bb2(i)px=pltop+plbot*(radius-yn(i))/(2*radius)fxn(i)=px*(b1+b2)/2*elseif,xn(i),gt,0 如xn(i)大于0(y轴右侧节点)，则b1=bb1(i)b2=bb2(i)px=-(pltop+plbot*(radius-yn(i))/(2*radius))fxn(i)=px*(b1+b2)/2*elsefxn(i)=0*endif*enddo进入求解层/solud,nodesum+1,ux,0,,nodesum+nodesum,1 弹簧节点x,y轴进行约束d,nodesum+1,uy,0,,nodesum+nodesum,1d,1+nodesum/2,ux,0 将圆环最下侧节点进行水平约束ACEL,0,10,0*do,i,1,nodesum,1 对管片节点施加节点力f,i,fx,fxn(i)f,i,fy,fyn(i)*enddo施加静水压力*do,i,1,nodesum,1ni=nelem(i,1) 得到单元I的起始节点号nj=nelem(i,2) 得到单元I的末节点号water1=(water-ny(ni)+radius)*10 单元I处i节点静水压力water2=(water-ny(nj)+radius)*10sfbeam,i,1,pres,-water1,-water2*enddoNROPT,FULL,, !采用全牛顿－拉普森法进行求解Allsel !选择所有内容Outres,all,all !输出所有内容solvefinish进入后处理/post1allsel,allesel,s,type,,1/AUTO,1 相当于fit/post1/pnum,node,0/RGB,INDEX,100,100,100, 0/RGB,INDEX, 80, 80, 80,13/RGB,INDEX, 60, 60, 60,14/RGB,INDEX, 0, 0, 0,15/REPLOT/COLOR,PBAK,OFF ! 关闭背景颜色! 画变形结果图/TITLE,deform diagrampldisp,1! 保存为图片/image,save,deform_1,jpg建立单元表etable,imoment,smisc,6etable,jmoment,smisc,12etable,ishear,smisc,2etable,jshear,smisc,8etable,iforce,smisc,1etable,jforce,smisc,7!更新单元表ETABLE,REFL!画弯矩分布图/TITLE,Bending moment diagramplls,imoment,jmoment,-1/image,save,Bending_moment_1,jpg! 画剪力分布图/TITLE,Shearing force diagram plls,ishear,jshear/image,save,Shearing_force_1,jpg ! 画轴力分布图/TITLE,Axial force diagram plls,iforce,jforce/image,save,Axial_force_1,jpg! 更新单元表!allsel,all!esel,s,type,,2!ETABLE,REFL!etable,f,smisc,1!pretab,f。
ANSYS隧道的衬砌计算
fk
tg
c
式中，τ—岩石抗剪强度，φ—岩石内摩擦角，
σ—正应力， c—粘结力
一、荷载及荷载组合
实际工程中：
fk
A
Rc 100
Rc—岩石单轴抗压强度 A—小于1的修正系数
近似计算为：
fk
A
Rc 100
普氏观测试验资料，制成各种岩石 fk分类表。 最坚硬的岩石 fk达20；最差(松、软) f只k 有0.3。
时的内水压力和外水压力。
一、荷载及荷载组合
特殊荷载：校核水位时内、外水压力、灌浆压力、温度荷 载、地震力、施工荷载等。
衬砌计算时常采用下列荷载组合： 基本组合：1、正常运行情况：山岩压力+衬砌自重+宣泄设
计洪水时内水压力+外水压力。 不同地段岩 计弹性抗力 北方 考虑温度荷载 石情况不同 不考虑弹性抗力 非寒冷地区 不考虑温度荷载 特殊组合：2、施工、检修情况：山岩压力+衬砌自重+可能出
（二）围岩的弹性抗力 当衬砌承受荷载向围岩方向变形时，将受到围岩的抵抗，这个抵抗
力叫弹性抗力。 弹性抗力的大小和性质与工程地质条件有密切的关系，坚固完整
的岩石，弹性抗力大；围岩软弱破碎，弹性抗力小，甚至不能利用。 为了减少山岩压力，有效地利用弹性抗力，常对围岩进行灌浆加
固，并填实衬砌与围岩间的空隙，以保证衬砌与围岩紧密相接。 弹性抗力的计算： 通常假定岩石为理想弹性体，按文克尔假定，认为岩石
缺点：与实际情况不太吻合。 2、采用衬砌常微分方程边值问题数值解法。 ★二类：有限元法 将围岩与衬砌视为整体，其准确性取决于计 算模型和原始参数E、φ、C等。 岩体性态复杂多变，上两者难以准确确定。
一、荷载及荷载组合
隧道及地下工程ANSYS分析地铁盾构隧道掘进施工过程三维仿真分析
第十章/TITLE, 3D analysis on shield tunnel in Metro ! 拟定分析标题/NOPR !菜单过滤设立/PMETH, OFF, 0KEYW, PR_SET, 1KEYW, PR_STRUC, 1 !保存结构分析部分菜单/COM,/COM, Preferences for GUI filtering have been set to display: 1．/COM, Structural2．材料、实常数和单元类型定义/clear !更新数据库/prep7 !进入前解决器et,1,solid45 !设立单元类型et,2,mesh200,6save !保持数据（2）定义模型中的材料参数。
!土体材料参数mp,ex,1,3.94e6 !地表层土弹性模量mp,prxy,1,0.35 !地表层土泊松比mp,dens,1,1828 !地表层土密度mp,ex,2,20.6e6 !盾构隧道所在地层参数mp,prxy,2,0.30mp,dens,2,2160mp,ex,3,500e6 !基岩地层参数mp,prxy,3,0.33mp,dens,3,2160!管片材料参数, 管片衬砌按各向同性计算mp,ex,4,27.6e9 !管片衬砌弹性模量mp,prxy,4,0.2 !管片衬砌泊松比mp,dens,4,2500 !管片衬砌密度!注浆层, 参数按水泥土取值mp,ex,5,1e9 !注浆层弹性模量mp,prxy,5,0.2 !注浆层泊松比mp,dens,5,2100 !注浆层密度save !保持数据3．建立平面内模型并划分单元（1）在隧道中心线定义局部坐标, 便于后来的实体选取。
local,11,0,0,0,0 !局部笛卡儿坐标local,12,1,0,0,0 !局部极坐标csys,11 !将当前坐标转换为局部坐标wpcsys,-1 !同时将工作平面转换到局坐标cyl4,,,,,2.7,90 !画部分圆半径为2.7cyl4,0,0,2.7,0,3,90 !画管片层部分圆cyl4,0,0,3,0,3.2,90 !画注浆层部分圆rectng,0,4.5,0,4.5 !画外边界矩形aovlap,all !做面递加nummrg,all !合并所有元素numcmp,all !压缩所有元素编号rectng,4.5,31.5,0,4.5 !画矩形面nummrg,all !合并所有元素numcmp,all !压缩所有元素编号save !保持数据（2）划分单元, 如图10-1所示。
ansys隧道荷载结构法命令流及设计师对隧道设计的经验总结
/prep7*set,shuxiang,276170 !竖向荷载*set,zuoce,125500 !水平荷载*set,youce,-125500*set,jizhongli,-100000 !地震集中力*set,uxishu,136.966e6*0.5 !水平弹性抗力系数*set,vxishu,87.4348e6*0.5 !竖向弹性抗力系数!定义单元类型、实常数及材料性质et,1,beam3r,1,0.3,0.00225,0.3mp,ex,1,3.1e10 !弹性模量mp,dens,1,2500 !密度mp,prxy,1,0.167 !泊松比!创建几何模型画点k,1,0,2.55k,2,-2.04,1.53k,3,-2.6,-0.15k,4,-2.319721,-2.023237k,5,-1.648852,-2.71152k,6,0,-2.98k,7,1.648852,-2.71152k,8,2.319721,-2.023237k,9,2.6,-0.15k,10,2.04,1.53k,1000,0,0 !中心点!画隧道轮廓线larc,1,2,1000,2.55larc,2,3,1000,2.8larc,3,4,1000,6.4larc,4,5,1000,1.04969larc,5,6,1000,5.2larc,6,7,1000,5.2larc,7,8,1000,1.04969larc,8,9,1000,6.4larc,9,10,1000,2.8larc,10,1,1000,2.55!选择左边的所有线allsel !选择所有的实体lsel,s,loc,x,0,-5 !选择左边的所有线:x从0到－5 lcomb,all !合并所选的线lesize,all,,,30 !把上面合并的线等分30段!选择右边的所有线allsellsel,s,loc,x,0,5lcomb,alllesize,all,,,30allsel !全选!划分单元lmesh,allnplot !显示节点!加弹簧!y方向*do,i,23,40PSPRNG,i,TRAN,vxishu,,-0.3,, , !tran-直线的*enddoPSPRNG,2,TRAN,vxishu,,-0.3,, , !固定结构用的!x方向*do,i,15,26PSPRNG,i,TRAN,uxishu,-0.3,,, ,*enddo*do,i,37,48PSPRNG,i,TRAN,uxishu,0.3,,, ,*enddoallselfinish!进入求解层，施加荷载，定义荷载步等/solu!施加约束nsel,s,,,2 !选择约束的节点nsel,s,loc,x,0dnsel,r,loc,y,-2.98d,all,ux !施加水平方向的约束allselfcum,add,, !一定要，使荷载能叠加!竖向荷载*do,i,3,15a=-shuxiang*0.5*abs(nx(i)-nx(i+1)) f,i+1,fy,af,i,fy,a*enddo*do,i,47,59a=-shuxiang*0.5*abs(nx(i)-nx(i+1)) f,i+1,fy,af,i,fy,a*enddo!对未循环的节点施加竖向荷载a=-shuxiang*0.5*abs(nx(1)-nx(3)) f,3,fy,af,1,fy,2*af,60,fy,a!水平荷载*do,i,3,30a=zuoce*0.5*abs(ny(i)-ny(i+1))f,i+1,fx,af,i,fx,a*enddo*do,i,32,59a=youce*0.5*abs(ny(i)-ny(i+1))f,i+1,fx,af,i,fx,a*enddof,1,fx,jizhongli !加集中力!对未循环的节点施加水平荷载a=youce*0.5*abs(ny(31)-ny(2))f,31,fx,-af,32,fx,a!施加重力acel,,9.8solve !求解finish!进入后处理，定义荷载工况并组合，输出图片和文本文件/post1!显示弯矩图etable,mi,smisc,6etable,mj,smisc,12plls,mi,mj,-1!显示轴力图etable,fi,smisc,1etable,fj,smisc,7plls,fi,fj,1!显示变形图pldisp,1综合该讨论的各家建议：隧道设计流程各家观点：1、我们做设计的时候也是这样做的,先用同济曙光采用地层结构法模拟开挖过程,对喷锚支护及开挖方法进行分析,选择一个合适的开挖顺序,然后用荷载结构法计算初衬的变形, 和二衬的承载能力和裂缝宽度。
基于ANSYS的隧道二衬内力计算[详细]
![基于ANSYS的隧道二衬内力计算[详细]](https://img.taocdn.com/s3/m/78a2f25b6294dd88d1d26b24.png)
单击Size Controls栏下Lines右边的Set，为线 的网格划分设置单元尺寸控制参数。选中L1 单击OK，在NDIV栏输入40，即一条线划分 为40个单元。单击OK，回到MeshTool对话 框。单击Mesh，选中L1，单击OK,网格划分 完毕。
15、查看单元
执行PlotCtrls>Numbering命令，如图
输入参数时，轴力两端为1和7，剪力为2和8， 弯矩为6和12。
• 执行Main Menu>General Postproc>Plot Results>Contour Plot>Line Element Res命令， 显示图形。
完
• 运用统计法计算，竖向均布作用为
•
• 式中 S ——围岩级别
• γ——围岩重度
•
ω——宽度影响系数，其值为：ω=1+i（B-5）
•
B ——坑道宽度
•
i ——此处取i=0.1。
• 注：此公式规定应用条件限制在采用在钻 爆法施工的深埋隧道，且H/B＜1.7的范围内 （H为坑道的高度，B为坑道宽度）。
—围岩的摩擦角。
ii) 求侧压力系数λ：
•
• iii)求围岩垂直均布松动压力q： •
• 式中 B —坑道跨度；
•
h —洞顶岩体覆盖层厚度；
•
γ —围岩重度。
• iv) 求围岩水平均布松动压力：
• 水平压力按梯形分布，则作用在隧道顶部和底 部的水平压力可直接写为：
e1 =γhλ
e2=γHλ
• 式中 H —地面到洞底岩体高度。
执行Main Menu>Preprocessor>Modeling>Delete
盾构隧道结构ansys计算方法
一、盾构隧道结构计算模型1、惯用法(自由圆环变形法)惯用法的想法早在1960年就提出了，在日本国内得到了广泛的应用。
惯用法假设管片环是弯曲刚度均匀的环，不考虑管片接头部分的柔性特征和弯曲刚度下降，管片截面具有同样刚度，并且弯曲刚度均匀的方法。
这种方法计算出的管片环变形量偏小，导致在软弱地基中计算出的管片截面内力过小，而在良好地基条件下计算出的内力又过大。
地层反力假设仅在水平方向上下45°范围内按三角形规律分布，这种模型可以计算出解析解。
P 0k δ2、修正惯用法在采用惯用法的60年代，怎样评价错缝拼装效应是一个问题。
如果错缝拼装管片，可弥补管片接头存在造成的刚度下降。
于是，在对带有螺栓接头的管片环进行多次核对研究时，首次引入了η-ξ对错缝拼装的衬砌进行内力计算，即为修正惯用法。
该法将衬砌视为具有刚度ηEI的均质圆环，将计算出的弯矩增大即（1+ξ）M，得到管片处的弯矩；将求出的弯矩减少即（1-ξ）M，得到接头处的弯矩。
其中η称为弯曲刚度有效率，ξ称为弯矩增加率，它为传递给邻环的弯矩与计算弯矩之比。
管片接头由于存在一些铰的作用，所以可以认为弯矩并不是全部经由管片接头传递，其一部分是利用环接头的剪切阻力传递给错缝拼装起来的邻接管片。
隧道纵向接头传递弯矩示意图二、管片计算荷载的确定1、荷载的分类衬砌设计所考虑的各种荷载，应根据不同的地质条件和设计方法进行假定并根据隧道的用途加以考虑。
衬砌设计所考虑的各种荷载见表所示。
衬砌设计荷载分类表2、计算断面选择●埋深最大断面●埋深最小断面●埋深一般断面●水位3、水土压力计算对于粘性土层，如西安地铁黄土地层、成都地铁二号线膨胀土地层等，地下水位以上地层荷载用湿容重计算，地下水位以下用饱和容重计算。
对于透水性较好的砂性地层，如西安地铁粗砂、中砂地层，成都地铁卵石土地层等，此时地下水位以上地层荷载用湿容重计算，地下水位以下用浮容重计算。
水土压力合算与分算，主要影响管片结构侧向荷载。
ANSYS地铁盾构隧道管片结构设计修改版
9.1管片衬砌结构力学分析过程9.1.1 前处理确定分析标题。
/TITLE, Mechanical analysis on segment lining of shield tunnel in Metro !分析标题/NOPR !菜单过滤设置/PMETH,OFF,0KEYW,PR_SET,1KEYW,PR_STRUC,1 !保留结构分析部分菜单/COM,/COM,Preferences for GUI filtering have been set to display:/COM, Structural1．材料、实常数和单元类型定义/PREP7 !进入前处理器ET,1,BEAM3 !设置梁单元类型，模拟管片衬砌ET,2, COMBIN14 !设置弹簧单元类型，模拟衬砌与地层相互作用R,1,0.3, 0.00225,0.3, , , , !设置梁单元几何常数，单位为mR,2,30e6, , , !设置弹簧单元几何常数，单位为Pa MPTEMP,,,,,,,, !设置材料模型MPTEMP,1,0MPDATA,EX,1,,34.5e9 !输入弹性模量，单位为PaMPDATA,PRXY,1,,0.2 !输入泊松比MPTEMP,,,,,,,, !设置材料模型MPTEMP,1,0MPDATA,DENS,1,,2500 !输入密度，单位为kgSAVE !保存数据库2．建立几何模型（1）创建关键点。
K,100,0,0,0, !通过编号和坐标创建关键点，圆心位置K,1,0.5925,2.7877,0,K,2,-0.5925,2.7877,0,K,3,2.7105,0.8807,0,K,4,-2.7105,0.8807,0,K,5,1.6752,-2.3057,0,K,6,-1.6752,-2.3057,0, !以上6个关键点为管片环向接头位置SAVE !保存数据（2）创建盾构隧道轮廓线。
第1章大型有限元软件ANSYS简介2 Larc,1,2,100,2.85 !通过两个端点和圆弧内侧任意一点以及半径画圆弧Larc,2,4,100,2.85Larc,4,6,100,2.85Larc,6,5,100,2.85Larc,5,3,100,2.85Larc,3,1,100,2.85 !以上为盾构隧道轮廓线的绘制SAVE !保存数据3．4. 创建弹簧节点所在的轮廓线（1）创建关键点。
盾构隧道结构计算模型简述
盾构隧道结构计算模型简述发布时间：2021-06-24T08:22:11.008Z 来源：《防护工程》2021年6期作者：武鹏[导读] 传统的隧道于地下工程结构计算方式主要有荷载-结构模型，地层结构模型。
近些年来，随着大量盾构隧道工程的出现，对于隧道结构的计算提出了新的要求。
虽然各种计算模型已经百花齐放，但各计算模型的优缺点，适用条件，在实际工程设计中仍然存在一定的混淆，本文从荷载计算、结构模拟的角度分析不同计算模型的特点、分类、适用条件，指出了其在实际工程设计中的适用性。
武鹏中国公路工程咨询集团有限公司北京市 100089摘要：传统的隧道于地下工程结构计算方式主要有荷载-结构模型，地层结构模型。
近些年来，随着大量盾构隧道工程的出现，对于隧道结构的计算提出了新的要求。
虽然各种计算模型已经百花齐放，但各计算模型的优缺点，适用条件，在实际工程设计中仍然存在一定的混淆，本文从荷载计算、结构模拟的角度分析不同计算模型的特点、分类、适用条件，指出了其在实际工程设计中的适用性。
1、盾构隧道荷载的计算理论地下工程结构的荷载计算，目前主要分为两类：荷载-结构模型和地层-结构模型。
1.1 荷载-结构模型荷载-结构模型默认围岩是一种松散体，是荷载的来源，而结构的作用只是被动承受荷载的荷载—结构模型；而地层-结构模型则认为围岩虽然是荷载的来源，但本身具有一定的承载能力，而结构的作用是对围岩的保护与补强，两者协同作用，共同承担荷载。
荷载-结构模型的前提是围岩因为工程的开挖而发生了较大的松弛或者崩塌，其已失去了承载能力，简言之，围岩是一种松散体，为支护结构“松动”压力的来源。
隧道结构设计的关键，即为确定围岩作用在支撑结构的主动荷载，长久以来，各国工程师，科研人员根据埋深不同，提出了太沙基理论、普氏理论等计算主动荷载，这些理论具有取值简单，适用性强的特点，在工程领域取得了广泛的应用。
确定了荷载后，即可运用结构力学、弹性力学等知识求解超静定结构的内力与变形，并由此确定安全系数。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、盾构隧道结构计算模型
1、惯用法(自由圆环变形法)
惯用法的想法早在1960年就提出了，在日本国内得到了广泛的应用。
惯用法假设管片环是弯曲刚度均匀的环，不考虑管片接头部分的柔性特征和弯曲刚度下降，管片截面具有同样刚度，并且弯曲刚度均匀的方法。
这种方法计算出的管片环变形量偏小，导致在软弱地基中计算出的管片截面内力过小，而在良好地基条件下计算出的内力又过大。
地层反力假设仅在水平方向上下45°范围内按三角形规律分布，这种模型可以计算出解析解。
P 0
k δ
2、修正惯用法
在采用惯用法的60年代，怎样评价错缝拼装效应是一个问题。
如果错缝拼装管片，可弥补管片接头存在造成的刚度下降。
于是，在对带有螺栓接头的管片环进行多次核对研究时，首次引入了η-ξ对错缝拼装的衬砌进行内力计算，即为修正惯用法。
该法将衬砌视为具有刚度ηEI的均质圆环，将计算出的弯矩增大即（1+ξ）M，得到管片处的弯矩；将求出的弯矩减少即（1-ξ）M，得到接头处的弯矩。
其中η称为弯曲刚度有效率，ξ称为弯矩增加率，它为传递给邻环的弯矩与计算弯矩之比。
管片接头由于存在一些铰的作用，所以可以认为弯矩并不是全部经由管片接头传递，其一部分是利用环接头的剪切阻力传递给错缝拼装起来的邻接管片。
隧
道
纵
向
接头传递弯矩示意图
二、管片计算荷载的确定
1、荷载的分类
衬砌设计所考虑的各种荷载，应根据不同的地质条件和设计方法进行假定并根据隧道的用途加以考虑。
衬砌设计所考虑的各种荷载见表所示。
衬砌设计荷载分类表
2、计算断面选择
埋深最大断面
埋深最小断面
埋深一般断面
水位
3、水土压力计算
对于粘性土层，如西安地铁黄土地层、成都地铁二号线膨胀土地层等，应采用水土压力合算的方式进行荷载计算。
此时，地下水位以上地层荷载用湿容重计算，地下水位以下用饱和容重计算。
对于透水性较好的砂性地层，如西安地铁粗砂、中砂地层，成都地铁卵石土地层等，应采用水土压力分算的方式进行荷载计算。
此时地下水位以上地层荷载用湿容重计算，地下水位以下用浮容重计算。
水土压力合算与分算，主要影响管片结构侧向荷载。
一般水土分算时侧向压力更大。
4、松弛土压力
将垂直土压力作为作用于衬砌顶部的均布荷载来考虑。
其大小宜根据隧道的覆土厚度、隧道的断面形式、外径和围岩条件等来决定。
考虑长期作用于隧道上的土压力时，如果覆土厚度小于隧道外径，一般不考虑地基的拱效应而采用总覆土压力。
但当覆土厚度大于隧道外径时，地基中产生拱效应的可能性比较大，可以考虑在计算时采用松弛土压力，一般采用泰沙基公式计算。
P 0
101
0/.tan 0/.tan 0110)1(tan )/1(B H K B H K e
P e K B c B h ϕϕγ
ϕγ--+--=
1
010/.tan 0/.tan 011)1(tan )/(B H K B H K v e
P e K B c B ϕϕϕ
γσ--+--=
)2
2
/4/cot(
01ϕπ+=R B
式中：—土的松动高度—0h
力—太沙基竖向松弛土压—v σ
ϕ
—内摩擦角
—
—
P
—上覆荷载
γ
—
—土的容重
c
—土的粘聚力
—
三、管片内力有限元计算单元选择
1、beam3单元
输入参数：
（1）实常数（Real Constants）
AREA –横截面积
IZZ –横截面惯性矩
HEIGHT –梁高
（2）材料参数（Material Properties）
EX——弹模
DENS——密度
PRXY——泊松比
其它选项采用默认即可
2、Link10单元
（1）Real Constants
AREA –横截面积，取单位面积1。
（2）Material Properties
EX ——弹模，需满足E*A=地层抗力刚度 其它选项采用默认即可 3、link1单元
（1）Real Constants
AREA – 横截面积，取单位面积1。
（2）Material Properties
EX ——弹模，为了不影响计算结果取足够小，如1。
四、结点力计算
1、均布荷载时
中间结点：
2/)(11-+-=n n x x p Fn
两端结点：
2/)(121x x p F -= 2/)(1--=N N N x x p F
2、梯形荷载时
中间结点：
[]))(2())(2(6
1
1111n n n n n n n n y y p p y y p p Fn -++-+=++--
两端结点：
[]))(2(61
12211y y p p F -+=
[]))(2(6
1
11---+=N N N N N y y p p F
五、圆形结构内力计算
1、计算参数
R=4.5m；B=1.5m；H=0.5m(厚)；
Pw1=;Pe1=;
Pw2=;Pe2=;
Q1=
Q2=
2、在自重情况下的结构内力分布
计算命令流为：
/prep7
et,1,beam3
et,2,link10
keyopt,2,3,1 !*****只受压******
R,1,,,, , , ,
R,2,1, ,
!定义混凝土衬砌材料
mp,ex,1,
mp,prxy,1,
mp,dens,1,2600
mp,ex,2,5e7
mp,prxy,2,
!建立几何模型
k,1,0,0,0
circle,1,
!设置划分单元大小为0.5m/个
Lsel,all
LESIZE,all,, , , , , , ,1
Type,1
Mat,1
Real,1
Lmesh,all
!复制建立地层弹簧
CSYS,1
wpcsys,-1
nsel,s,loc,x,,
NGEN,2,500,all, , ,1 , , ,1,
allsel,all
*do,i,1,60
type,2
real,2
mat,2
e,i,i+500
*enddo
nsel,s,loc,x,,
d,all,all
nsel,s,,,32
D,all, , , , , ,UX, , , , ,
allsel,all
csys,0
wpcsys,-1
/sol
ACEL,0,10,0,!********施加重力命令********* solve
/post1
ETABLE, ,SMISC, 1
ETABLE, ,SMISC, 7
ETABLE, ,SMISC, 2
ETABLE, ,SMISC, 8 ETABLE, ,SMISC, 6 ETABLE, ,SMISC, 12
工况一：πg＋地层反力作为主动荷载施加
工况二：πg＋地层反力作为被动荷载施加
结果比较：。