双线隧道施工过程ANSYS数值模拟分析

/NOPR双线铁路隧道施工过程仿真分析KEYW,PR_STRUC,1!保留结构分析部分菜单/PREP7!进入前处理器ET,1,PLANE42ET,2,BEAM3KEYOPT,1,3,2!设置成平面应变模型!定义模型中的相关参数MP,EX,1,6E6!地表上MP,PRXY,1,0.26MP,DENS,1,1760MP,EX,2,10E6!土层2MP,PRXY,2,0.35MP,DENS,2,1910MP,EX,3,15E6!土层3MP,PRXY,3,0.28MP,DENS,3,2040MP,EX,4,14E6!土层4MP,PRXY,4,0.29MP,DENS,4,2075MP,EX,5,30E6!土层5MP,PRXY,5,0.3MP,DENS,5,2150MP,EX,6,50E6!超前小导管加固区MP,PRXY,6,0.2MP,DENS,6,3000MP,EX,7,20E6!锚杆加固区MP,PRXY,7,0.2MP,DENS,7,1910MP,EX,8,210E9!钢支撑MP,PRXY,8,0.17MP,DENS,8,7800MP,EX,9,21E9!初期支护MP,PRXY,9,0.2MP,DENS,9,2200!定义实常数。模型中采用梁单元来模拟初期支护和临时支护，其中初期支护为C20喷射混凝土，临时支护包括中隔墙的立柱和临时仰拱，为25a号工字钢。定义其实常数的命令流如下：R,1,0.009546,0.0001005,0.25,0,0,0!工字钢实常数，按每米2根算R,2,0.35,3.573E-3,0.35!初期支护实常数!建立几何模型。由于CRD法开挖隧道的过程相当复杂，建立几何模型时通常采用自底向上的建模方法，即先根据施工设计图纸找出一些关键点，再由这些关键点生成圆弧等段线，再由这些线段生成面，最后再对这些面进行单元划分。K,1,-5.65,4.35,0!定义关键点K,2,-2.84,6.17,0K,3,0.75,6.65,0K,4,3.47,5.90,0K,5,5.65,4.35,0K,6,6.24,3.66,0K,7,6.83,2.69,0K,8,-0.45,4.54,0K,9,-1.11,1.77,0K,10,-6.21,3.71,0K,11,-6.83,2.69,0K,12,-7.41,0.73,0K,13,-7.31,-1.3,0K,14,-7.16,-1.86,0K,15,-6.98,-2.37,0K,16,-6.46,-3.19,0K,17,-5.63,-3.84,0K,18,0,-5.23,0K,19,1.10,-5.18,0K,20,4.02,-4.54,0K,21,5.63,-3.84,0K,22,6.44,-3.21,0K,23,6.98,-2.37,0K,24,7.17,-1.83,0K,25,7.31,-1.3,0K,26,7.41,0.72,0K,27,-1.14,0.65,0K,28,-0.6,-2.23,0K,29,-0.3,-2.98,0K,30,-5.79,4.49,0K,31,-2.92,6.35,0K,32,0.75,6.85,0K,33,3.66,6.02,0K,34,5.79,4.49,0K,35,-7.6,0.79,0K,36,-7.16,-2.44,0K,37,-6.6,-3.34,0K,38,-5.72,-4.02,0K,39,-1.84,-5.29,0K,40,1.12,-5.38,0K,41,4.09,-4.73,0K,42,5.72,-4.02,0K,43,6.59,-3.34,0K,44,7.16,-2.44,0K,45,7.60,0.76,0K,501,-3.72,2.02,0K,502,2.73,1.89,0K,503,-2.9,-2.09,0K,504,3.30,-2.04,0LARC,11,1,10!由关键点生成各段圆弧LARC,1,3,2!以关键点1和3为端点，关键点2为中间点生成圆弧LARC,3,5,4LARC,5,7,6LARC,9,3,8LARC,13,11,12LARC,15,13,14LARC,17,15,16LARC,19,17,18LARC,21,19,20LARC,23,21,22LARC,25,23,24LARC,7,25,26LARC,28,9,27LARC,19,28,29LARC,7,9,502LARC,9,11,501LARC,28,13,503LARC,25,28,504LARC,30,32,31LARC,32,34,33LARC,36,30,35LARC,38,36,37LARC,40,38,39LARC,42,40,41LARC,44,42,43LARC,34,44,45SAVEL,3,32!由关键点生成各段直线L,1,30L,5,34L,15,36L,19,40L,23,44LCOMB,23,24,0!将线23和24合并LCOMB,25,26,0!将线25和26合并NUMMRG,ALL,,,,LOW!合并所有元素NUMCMP,ALLK,48,-8.12,6.82,0!再一次定义关键点K,49,-10.61,2.55,0K,50,-10.21,-3.71,0K,51,-5.51,-8.03,0K,52,0,-9.36,0K,53,6.90,-7.19,0K,54,10.21,-3.71,0K,55,10.61,2.55,0K,56,8.12,6.82,0K,57,-8.12,9.27,0K,58,0.75,9.27,0K,59,8.12,9.27,0K,60,-14,9.27,0K,61,-14,-11.33,0K,62,14,-11.33,0K,63,14,9.27,0K,64,-14,-3.71,0K,65,0,-11.33,0K,66,14,-3.71,0LARC,48,50,49!再一次通过3点生成圆弧LARC,50,52,51LARC,52,54,53LARC,54,56,55L,48,57!再一次由关键点生成各段直线L,57,58L,58,59L,59,56L,34,56L,42,54L,39,52L,36,50L,30,48L,32,58L,57,60L,60,64L,64,61L,61,65L,65,62L,62,66L,66,63L,63,59L,50,64L,52,65L,54,66AL,1,2,5,17!由线生成面AL,5,3,4,16AL,14,16,13,19AL,15,19,12,11,10AL,7,8,9,15,18AL,6,17,14,18AL,27,20,26,2AL,26,21,28,3AL,28,25,31,12,13,4AL,31,24,30,10,11AL,30,23,29,8,9AL,29,22,27,1,6,7AL,44,36,37,45,20AL,45,38,39,40,21AL,40,35,41,25AL,41,34,42,24AL,33,42,23,43AL,32,43,22,44AL,54,47,46,36,32AL,54,48,49,55,33AL,55,34,56,51,50AL,56,35,39,53,52SAVE!划分隧道部分单元TYPE,1!选择单元类型MAT,1!选择材料类型LESIZE,1,,,4!设定线划分的单元数LESIZE,2,,,12!将2号线划分成12个单元LESIZE,5,,,8LESIZE,17,,,8LESIZE,3,,,10LESIZE,4,,,4LESIZE,16,,,10LESIZE,18,,,8LESIZE,19,,,10LESIZE,7,,,2LESIZE,12,,,2LESIZE,8,,,4LESIZE,11,,,4LESIZE,9,,,10LESIZE,10,,,8LESIZE,15,,,8LESIZE,6,,,6LESIZE,13,,,6LESIZE,14,,,6MSHAPE,1,2D!设置单元划分的形状，对面采用三角形单元划分MSHKEY,0!采用FREE方式划分单元AMESH,1,2!采用上述设置划分面1和面2AMESH,4,5!采用上述设置划分面4及面5MSHAPE,0,2D!采用四边形单元划分MSHKEY,1!采用MAP映射网格方式划分单元AMESH,3!采用上述设置划分面3及6AMESH,6LSEL,S,,,26,31,1!选择线LESIZE,ALL,,,1!设置单元数ALLSEL,ALL!选择所有元素LESIZE,21,,,10!设置单元数LESIZE,20,,,12LESIZE,22,,,12LESIZE,23,,,14LESIZE,24,,,12LESIZE,25,,,12MSHAPE,0,2D!采用四边形单元划分MSHKEY,0!采用FREE方式划分单元AMESH,7,12,1!采用上述设置划分面LSEL,S,,,40,45,1!选择线LESIZE,ALL,,,5!设置单元数ALLSEL,ALL!选择所有元素LESIZE,32,,,12!设置单元数LESIZE,33,,,14LESIZE,34,,,12LESIZE,35,,,12LESIZE,32,,,12LESIZE,36,,,4LESIZE,37,,,8LESIZE,38,,,6LESIZE,39,,,4MSHAPE,0,2D!采用四边形单元划分MSHKEY,1!采用MAP方式划分单元AMESH,15,18,1如果此处使用映射网格划分对划分的结果会产生怎样的结果SAVE!AMESH,13,14,1!此命令不可执行!？MSHAPE,0,2D$MSHKEY,0$AMESH,13,14,1次命令与AMAP,13,30,32,58,48AMAP,14,32,34,56,58的本质区别是什么？AMAP,13,30,32,58,48!采用MAP方式划分面13AMAP,14,32,34,56,58LSEL,S,,,54,56,1!选择线LESIZE,all,,,4!设置单元数ALLSEL,ALL!选择所有元素LESIZE,46,,,4!设置单元数LESIZE,53,,,4LESIZE,48,,,5LESIZE,51,,,5LESIZE,47,,,8LESIZE,52,,,8LESIZE,49,,,9LESIZE,50,,,7AMAP,19,64,50,48,57!采用MAP方式划分面19到面22AMAP,20,64,50,52,65AMAP,21,65,52,54,66AMAP,22,66,54,56,59RECTNG,-14,-60,-3.71,9.27!绘制矩形面RECTNG,-14,-60,-11.33,-3.71RECTNG,-14,-60,-11.33,-40RECTNG,-14,0,-11.33,-40RECTNG,14,0,-11.33,-40ALLSEL,ALL!选择所有元素NUMMRG,ALL,,,,LOW!合并所有元素NUMCMP,ALL!压缩所有元素编号LESIZE,59,,,8LESIZE,61,,,5!设置单元数LESIZE,65,,,9LESIZE,67,,,7LESIZE,58,,,12LESIZE,57,,,12LESIZE,60,,,12LESIZE,62,,,12LESIZE,64,,,8LESIZE,63,,,8LESIZE,66,,,8LESIZE,68,,,8MSHAPE,0,2D!采用四边形单元划分MSHKEY,1!采用MAP方式划分单元AMESH,23,27,1!对面划分单元ASEL,S,AREA,,23,25,1ARSYM,X,ALL!对称所有元素ALLSEL,ALL!选择所有元素NUMMRG,ALL,,,,LOW!合并所有元素NUMCMP,ALL!压缩所有元素编号ALLSEL,ALL!选择所有元素NSEL,S,LOC,Y,-40,-22.1!为第5层土赋材料参数ESLN,SMPCHG,5,ALLNSEL,S,LOC,Y,-22,-11.5!为第4层土赋材料参数ESLN,SMPCHG,4,ALLNS