第6章 高速公路分离式偏压隧道施工过程仿真分析

66..11施工过程仿真分析1．材料、实常数和单元类型定义/TITLE，Mechanicalanalysisonrailwaytunnel1ndlining!确定分析标题/NOPR!菜单过滤设置/PMETH，OFF，0KEYW，PR_SET，1KEYW，PR_STRUC，1!保留结构分析部分菜单/COM，/COM，PreferencesforGUIfilteringhavebeensettodisplay:/COM，Structural/PREP7!进入前处理器ET,1,PLANE42!设置实体单元类型KEYOPT,1,1,0KEYOPT,1,2,0KEYOPT,1,3,2!设置为平面应变模式KEYOPT,1,5,0KEYOPT,1,6,0ET,2,PLANE42!设置实体单元类型KEYOPT,1,1,0KEYOPT,1,2,0KEYOPT,1,3,2!设置为平面应变模式KEYOPT,1,5,0KEYOPT,1,6,0ET，3，BEAM3!设置梁单元类型R,1,0.25,1.5e-3,0.25,,,,!设置梁单元几何常数MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,0.8e9!输入弹性模量（围岩）MPDATA,PRXY,1,,0.38!输入泊松比（围岩）MPDATA,DENS,1,,2000!输入密度（围岩）TB,DP,1,,,!采用DP准则进行弹塑性分析TBMODIF,1,1,0.2e6!输入凝聚力（围岩）TBMODIF,1,2,30!输入摩擦角（围岩）TBMODIF,1,3,MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,2,,3.2e9!输入弹性模量（加固圈）MPDATA,PRXY,2,,0.32!输入泊松比（加固圈）第1章大型有限元软件ANSYS简介2MPDATA,DENS,2,,2200!输入密度（加固圈）TB,DP,2,,,!采用DP准则TBMODIF,1,1,0.5e6!输入凝聚力（加固圈）TBMODIF,1,2,35!输入摩擦角（加固圈）TBMODIF,1,3,MPDATA,EX,3,,27.5e9!输入弹性模量（喷射混凝土）MPDATA,PRXY,3,,0.2!输入泊松比（喷射混凝土）MPDATA,DENS,3,,2500!输入密度（喷射混凝土）SAVE!保存数据库2．建立几何模型创建关键点。K,1,0,0,,!创建计算区域轮廓线关键点K,2,65,-22.3,,K,3,-65,-22.3,,K,4,65,37.1,,K,5,19.6,37.1,,K,6,-65,-11.8,,K,7,-44.6,-22.3,,K,8,-44.6,0,,K,9,65,0,,K,10,15,-22.3,,K,11,-15,-22.3,,K,12,-15,17.1,,K,13,15,34.4,,K,14,-15,1.169,,!创建隧道轮廓线关键点，上半圆圆心K,15,-18.408,0.22,,!左下拱脚圆心K,16,-15.022,7.887,,!仰拱圆心K,17,0,25.742,,K,18,0,-22.3,,K,27,-11.636,0.22,,!右下拱脚圆心SAVE!保存数据创建线。Circle,14,5.4!画上半圆Circle,15,1.9!画左下拱脚圆Circle,27,1.9!画右下拱脚圆LDELE,5,6,1!删除左右下拱脚圆的部分线LDELE,8,11,1LOVLAP,7,3!作叠加布尔运算LOVLAP,12,4LDELE,6,10,2LDELE,7LDELE,13,16,3NUMMRG,ALL,,,,LOW!合并所有元素第1章大型有限元软件ANSYS简介3NUMCMP,ALL!压缩所有编号FLST,2,3,4,ORDE,3!作加运算FITEM,2,4FITEM,2,5FITEM,2,6LCOMB,P51X,,0FLST,2,3,4,ORDE,3!作加运算FITEM,2,3FITEM,2,7FITEM,2,9LCOMB,P51X,,0NUMMRG,ALL,,,,LOW!合并所有元素NUMCMP,ALL!压缩所有编号circle,16,10.178!画仰拱圆LOVLAP,3,7!作叠加布尔运算LOVLAP,4,8LDELE,19!删除部分线LDELE,5,7,1LDELE,10,11,1LDELE,14,16,2FLST,2,4,4,ORDE,4!做线加布尔运算FITEM,2,9FITEM,2,12FITEM,2,-13FITEM,2,15LCOMB,P51X,,0FLST,2,4,4,ORDE,4!做线加布尔运算FITEM,2,3FITEM,2,17FITEM,2,-18FITEM,2,20LCOMB,P51X,,0Circle,14,8.9!画上半圆(加固圈)Circle,15,5.4!画左下拱脚圆(加固圈)Circle,27,5.4!画右下拱脚圆(加固圈)ldele,8,12,2ldele,13,15,1LOVLAP,6,11!作叠加布尔运算LOVLAP,7,16ldele,10,14,2!删除部分线ldele,17,19,2ldele,11circle,16,13.678!画仰拱圆(加固圈)LOVLAP,11,13!作叠加布尔运算第1章大型有限元软件ANSYS简介4LOVLAP,12,18ldele,10,11,1!删除部分线ldele,23,26,3ldele,7,17,10ldele,14,22,8FLST,2,5,4,ORDE,5!做线加布尔运算FITEM,2,8FITEM,2,15FITEM,2,-16FITEM,2,19FITEM,2,21LCOMB,P51X,,0FLST,2,5,4,ORDE,5!做线加布尔运算FITEM,2,6FITEM,2,13FITEM,2,20FITEM,2,24FITEM,2,-25LCOMB,P51X,,0FLST,3,8,4,ORDE,4!对称操作，由左隧道镜像成右隧道。FITEM,3,1FITEM,3,-6FITEM,3,8FITEM,3,-9LSYMM,X,P51X,,,,0,0NUMMRG,ALL,,,,LOW!合并所有元素NUMCMP,ALL!压缩所有编号SAVE!所画出的隧道轮廓图见图6-7所示。k,56,0,1.169,0,!创建关键点，为画其它直线作准备k,57,30.044,1.169,0,k,58,30.044,-22.3,0,k,59,30.044,37.1,0,k,60,15,1.169,0,LSTR,3,7!画计算区域线和分块线LSTR,7,8LSTR,8,6LSTR,6,3LSTR,7,11LSTR,11,18LSTR,18,10LSTR,10,58LSTR,58,2LSTR,2,9LSTR,9,4LSTR,4,59第1章大型有限元软件ANSYS简介5LSTR,59,5LSTR,5,13LSTR,13,17LSTR,17,12LSTR,12,8LSTR,21,14LSTR,14,19LSTR,14,20LSTR,14,31LSTR,44,60LSTR,60,46LSTR,60,45LSTR,60,47LSTR,45,53LSTR,46,54LSTR,47,55LSTR,44,52LSTR,52,56LSTR,54,57LSTR,57,59LSTR,57,58LSTR,55,10LSTR,56,17LSTR,56,18LSTR,53,13LSTR,19,32LSTR,32,56LSTR,51,11LSTR,33,12LSTR,20,33LSTR,21,34LSTR,34,8LSTR,31,51LSTR,57,9/PNUM,KP,0!不显示关键点编号/PNUM,LINE,1!显示线编号/PNUM,AREA,0!不显示面编号/PNUM,VOLU,0!不显示体编号/PNUM,NODE,0!不显示节点编号!*/PNUM,ELEM,0!不显示单元编号/REPLOT!重新绘制!*LPLOT!绘制线图，带编号，如图6-8所示。SAVE!保存数据库第1章大型有限元软件ANSYS简介6创建面。本次的计算区域为横向130m，竖向为60m，即左右两侧计算边界为4倍左右双线隧道总跨度，下部边界为2倍隧道总高度。Al,35,1,36,!采用线创建面，依次创建27个面Al,36,2,34,Al,34,9,37,Al,37,3,35,!左隧道内核心土4个面绘制完成Al,54,4,58,1,Al,58,5,59,2,Al,59,8,61,9,Al,61,6,54,3,!左隧道加固圈土4个面绘制完成Al,39,10,40,Al,40,7,38,Al,38,11,41,Al,41,16,39,!右隧道内核心土4个面绘制完成Al,43,13,42,10,Al,42,12,45,7,Al,45,14,44,11,Al,44,15,43,16,!右隧道加固圈土4个面绘制完成Al,55,51,32,57,4,Al,57,33,60,5,Al,60,18,21,56,8,Al,56,22,52,55,6,!左隧道外围土体4个面绘制完成Al,47,48,29,30,53,13,Al,53,31,51,46,12,Al,46,52,23,50,14,Al,50,24,49,47,15,!右隧道外围土体4个面绘制完成Al,18,19,20,17,!左侧边境1个面绘制完成Al,49,25,26,62,Al,62,27,28,48,!右侧边境2个面绘制完成/PNUM,KP,0/PNUM,LINE,0/PNUM,AREA,1APLOT!绘制带编号的面,生成的面如图6-9所示。SAVE!保存数据,3．创建网格模型设置线被划分的单元数。LESIZE,_Y1,,,8,,,,,1lesize,1,,,8,,,,,1!设置单元大小,即L1线划分成8个单元TYPE,1!设置将要创建单元的类型MAT,1!设置将要创建单元的材料REAL,1!设置将要创建单元的几何常数/PNUM,LINE,1!显示线编号LPLOT!绘制线图，带编号和单元数。第1章大型有限元软件ANSYS简介7Save!按照“lesize,1,,,8,,,,,1”格式，将其他61条线也进行被划分单元数设置，其线的编号和被划分的单元数如表6-3所示，划分好后如图6-10所示。表6-1设置线单元数NL(NDIV)NL(NDIV)NL(NDIV)NL(NDIV)NL(NDIV)358398544554514368408584574324348388594604224378418614564524187843447418428108424534214981184544642043816844450419448128314301174581382642922446814825104854948815862102752810338234将所有面划分单元，其单元图如图6-11所示。AMAP,1,14,19,20,!划分面1，依次单击关键点14、19、20进行映射划分AMAP,2,14,20,21,AMAP,3,14,21,31,AMAP,4,14,31,19,!面1到4，为左隧道内开挖土体AMAP,9,60,46,45,AMAP,10,60,45,44,AMAP,11,60,44,47,AMAP,12,60,47,46,!面9到12，为右隧道内开挖土体AMAP,7,21,34,51,31AMAP,8,31,51,32,19!面7和8，为左隧道仰拱下土体AMAP,15,44,52,55,47AMAP,16,47,55,54,46!面15和16，为右隧道仰拱下土体TYPE,1!转换为加固圈土体单元类型和材料MAT,2AMAP,5,19,32,33,20!划分左隧道加固圈土体AMAP,6,20,33,34,21AMAP,13,46,54,53,45!划分右隧道加固圈土体AMAP,14,45,53,52,44TYPE,1!转换为普通土体单元类型和材料MAT,1AMAP,17,32,56,12,33!划