在不规则底层下三维隧道开挖阶段的APDL分析Word文件下载.docx

1、软岩E_R_2=50000/TONF_N !U_R_2=0.3 !DENS_R_2=2.3/TON_KG !密度，KG/M-3C_R_2=2/TONF_N !FI_R_2=33 !FTK_2=20/TONF_N !硬岩E_R_3=200000/TONF_N !U_R_3=0.25 !DENS_R_3=2.5/TON_KG !C_R_3=20/TONF_N !FI_R_3=35 !FTK_3=300/TONF_N !抗拉强度 喷射混凝土E_1=1.5E5/TONF_N !U_1=0.2 !DENS_1=2.4/TON_KG !ALPX_1=1E-6 !热膨胀系数ALPY_1=1E-6 !ALPZ

2、_1=1E-6 !PH_H=0.16 !锚喷的高度M25锚杆，ANCHOR RODE_2=2E7/TONF_N !U_2=0.3 !DENS_2=7.85/TON_KG !ALPX_2=1.2E-5 !AM=0.0005067 !锚杆截面面积M2LM=4 !锚杆长度M=前处理=/PREP7=定义单元类型=ET,1,95 !定义SOLID95单元，岩土单元ET,2,93 !定义SHELL93单元，喷射混凝土单ET,3,8 !定义LINK8单元，锚杆单元ET,4,82 !定义辅助平面82单元 =定义材料属性=岩层断面材料编号：土为1，软岩为2，硬岩为3*DO,I,1,3,1 MP,EX,I,E_R

3、_%I% !围岩弹性模量 MP,PRXY,I,U_R_%I% !围岩泊松比 MP,DENS,I,DENS_R _%I% !围岩密度 TB,DP,I TBDATA,1,C_R_%I%,FI_R_%I%, !内聚力和摩擦角*ENDDOTB,BKIN,I !等向随动强化TBDATA,1,FTK_%I%喷射混凝土和锚杆混凝土喷射为4，锚杆为5MP,EX,4,E_1 !混凝土弹性模量MP,PRXY,4,U_1 !混凝土泊松比MP,DENS,4,DENS_1 !混凝土密度MP,ALPX,4,ALPX_1 !热传递系数MP,ALPY,4,ALPY_1 !MP,ALPZ,4,ALPZ_1 !MP,EX,5,E

4、_2 !MP,PRXY,5,U_2 !MP,DENS,5,DENS_2 !MP,ALPX,5,ALPX_2 !=定义实常数=R,1, !岩层深度R,2,PH_H,PH_H,PH_H ,PH_H !喷射混凝土厚度的实常数PI=ACOS（-1）R,3,AM !锚杆=几何模型创建=右下角为坐标系原点=三心园的方式创建创建隧道截面D0=12 !岩层长度L1=11.7 !隧道底部宽度R1=6.5 !R2=6 !A1=60A2=55D_MG=1.5 !锚杆间距NM=13 !奇数土层尺寸H1=60 !硬岩层高度W=110 !土层宽度X1=40$Y1=30 !隧道位置 施工阶段的参数D_SD=2 !每个施工阶

5、段的深度DN=D0/D_SD土层的参数NX=11 !X方向点的个数NZ=DN+1 !Z方向点的个数DDX=W/（NX-1）DDZ=D0/（NZ-1）创建隧道截面采用高阶几何创建这样就可以不用考虑关键点编号的问题*AFUN,DEG/PNUM,KP,ONNM1=NINT（A1*R1*PI/180/D_MG-MOD（A1*R1*PI/180,D_MG）NMT=NM1*2+1J1=D_MG/R1*180/PIJ2=A1-NM1*J1PL1=（A1-J1*NM1）*R1*PI/180PL2=D_MG-PL1J3=D_MG/R2*180/PIJ4=PL2/R2*180/PINM2=（NM-1-NM1*2）

6、/2NMB=NM2*2PCIRC,R1,R1+LM,90-A1,90+A1CSYS,1K,100,R1-R2,90-A1$K,101,R1-R2,90+A1CSYS,0LOCAL,11,1,KX（100）,KY（100）,WPCSYS,11PCIRC,R2,R2+LM,90-A1,90-A1-J4-（NM2-1）*J3LOCAL,12,1,KX（101）,KY（101）,WPCSYS,12PCIRC,R2,R2+LM,90+A1,90+A1+J4+（NM2-1）*J3AADD,ALLNUMMRG,KPALLSELLSEL,U,LENGTH,LMLSEL,R,LOC,X,0,R1LCOMB,AL

7、LLSEL,R,LOC,X,R1+1,R1+LMWPROTA,90WPROTA,90-A1-J4-（NM2-1）*J3*DO,I,1,NM2-1,1WPROTA,J3ASBW,ALLWPCSYS,0WPROTA,90-A1+J2$ASBW,ALL*DO,I,1,2*NM1,1WPROTA,J1WPROTA,90+A1WPROTA,J4改变线的方向CSYS,11$LSEL,S,LOC,X,R2$LSEL,R,LOC,Y,90-A1-J4-（NM2-1）*J3,90-A1-J4-（NM2-1）*J3+J3LREVERSE,ALLCSYS,12$LSEL,S,LOC,X,R2$LSEL,R,LOC,

8、Y,90+A1+J4+（NM2-1）*J3,90+A1+J4+（NM2-1）*J3-J3CSYS,11K,200,R2,90-A1-A2L,200,KP（R2,90-A1-J4-（NM2-1）*J3,0）CSYS,12K,201,R2,90+A1+A2L,KP（R2,90+A1+J4+（NM2-1）*J3,0）,201L,201,200ADELE,ALLWPCSYS,-1=定义硬岩层=LOCAL,13,0,KX（200）+X1-W,KY（200）-Y1,-D0 !定义13坐标系 LOCAL,14,0,KX（200）+X1-W,KY（200）-Y1,0 !定义14坐标系WPCSYS,14 !RE

9、CT,0,W,0,H1ASBL,ALL,ALLASEL,S,14CM,AA1,AREAASEL,INVEASEL,U,15CM,AA2,AREA切割 水平切割CSYS,0$WPCSYS,-1WPROTA,90$WPOFFS,6$ASBW,ALLWPOFFS,-（11+R1+LM）$ASBW,ALL竖直切割WPROTA,90CSWPLA,19,0ASEL,U,LOC,X,0,20WPOFFS,-（R1+LM+5）$ASBW,ALLWPOFFS,2*（R1+LM+5）$ASBW,ALLCSYS,14LSEL,S,LOC,Y,H1$LDIV,ALL,NX-1ASEL,S,25$CM,AA3,AREA

10、ASEL,INVE$CMSEL,U,AA1$CMSEL,U,AA2CM,AA4,AREA=拉伸成直线隧道=K,1E3,$K,1E3+1,-D0$L,1E3,1E3+1LSEL,S,LENGTH,D0$LSEL,R,LOC,Y,0$LDIV,ALL,DNCM,LL1,LINEK,2E3,$K,2E3+1,-D0$L,2E3,2E3+1CM,LL2,LINEK,3E3,$K,3E3+1,-D0$L,3E3,3E3+1CM,LL3,LINECMSEL,S,AA1VDRAG,ALL,LL1CM,SD_V,VOLU !为隧道部分编组VSEL,NONECMSEL,S,AA2CM,MG_V,VOLU !为锚

11、杆部分编组CMSEL,S,AA3VDRAG,ALL,LL2CM,YYC_V0,VOLU !为周围硬岩层编组CMSEL,S,AA4VDRAG,ALL,LL3CM,YYC_V1,VOLU !这样可以实现一些不重要关键点的编号进行重新编号KSEL,U,200,201NUMCMP,KP/PNUM,KP,OFF/VIEW,1,1,1,1CMSEL,S,SD_VCMSEL,S,MG_VCMSEL,S,YYC_V0CMSEL,S,YYC_V1定义锚杆LSEL,S,LENGTH,LM$LSEL,R,TAN1,ZCM,MG_L,LINE*DO,I,1,DN,1CMSEL,S,MG_LLSEL,R,LOC,Z,-

12、D_SD-（I-1）*D_SDCM,MG_L%I%,LINE定义施工开挖施工隧道节段VSEL,R,LOC,Z,-（I-1）*D_SD,-I*D_SDCM,SD_V%I%,VOLU定义喷射混凝土ASLV,S$LSLA,SLSEL,R,TAN1,ZASLL,U,1CM,PH_A,AREACM,PH_L,LINECMSEL,S,PH_AASEL,R,LOC,Z,-（I-1）*D_SD,-I*D_SDCM,PH_A%I%,AREACM,YYC_V,VOLU !定义隧道及其周围部分的体积组=生成其他岩层=软岩层=先生成导线和引线D1=60NX1=11 !NZ1=11 !DDX1=W/（NX1-1）DDZ

13、1=D1/（NZ1-1）CSYS,13$WPCSYS,13ALLSEL$ASEL,NONE$VSEL,NONE$LSEL,NONEK,1E4+1,0,80K,1E4+2,30,78K,1E4+3,60,84K,1E4+4,90,79K,1E4+5,110,80K,1E4+6,115,80.1BSPLIN,1E4+1,1E4+2,1E4+3,1E4+4,1E4+5,1E4+6*GET,DX1,LINE,0,NUM,MINK,1E4+11,0,80,0K,1E4+12,0,78,20K,1E4+13,0,83,40K,1E4+14,0,80,60K,1E4+15,0,80.1,65BSPLIN,1

14、E4+11,1E4+12,1E4+13,1E4+14,1E4+15*GET,DX2,LINE,0,NUM,MAXADRAG,DX1,DX2CSYS,13$WPCSYS,1,13WPOFFS,D1ASEL,S,LOC,Z,D1,D1+200ASEL,R,LOC,Y,H1+1,H1+50ADELE,ALL,1ASEL,S,LOC,Y,H1+1,H1+50WPOFFS,W$WPROTA,90ASEL,S,LOC,X,W,W+200CM,RYC_A1,AREA !定义线LSLA,SKSLL,SVEXT,ALL,0,-50NUMCMP,LINECM,RYC_V1,VOLU删除不需要的部分WPOFFS,D

15、0$VSBW,ALLVSEL,R,LOC,Z,0,D0VSEL,S,LOC,Z,D0,D1VDELE,ALL,1=土层和软岩层的生成=CSYS,13WPCSYS,13WPOFFS,H1-4WPROTA,90,CMSEL,S,RYC_V1VSBW,ALLVSEL,S,LOC,Y,H1-4,200CMSEL,A,YYC_VVSEL,INVECSYS,19$WPCSYS,1,19 !WPOFFS,-（R1+LM+5）$VSBW,ALLWPOFFS,2*（R1+LM+5）$VSBW,ALLCM,RYC_V,VOLUVSEL,S,LOC,X,0,100VPTN,ALLVSEL,R,LOC,Y,0,H1C

16、M,YYC_V,VOLU*DO,I,1,NZ-1,1CMSEL,U,YYC_VWPOFFS,DDZCMSEL,U,SD_VCMSEL,U,MG_VCMSEL,U,YYC_V0CM,YYC_V1,VOLU=网格划分=NMG=2NPH=1E_SIZE1=2E_SIZE2=7E_EIZE3=5隧道及附近的划分TYPE,1$MAT,3$REAL,1CMSEL,A,SD_VCMSEL,A,MG_VASLV,S$LSLA,S$KSLL,SLSEL,R,TAN1,X$LSEL,R,TAN1,YLESIZE,ALL,1LSEL,R,LENGTH,LMLESIZE,ALL,NMG隧道周围CMSEL,S,YYC_

17、VCMSEL,U,YYC_V1ASEL,R,LOC,Z,0TYPE,4ESIZE,E_SIZE1MSHKEY,2MSHAPE,0,2DAMESH,ALLCMSEL,A,RYC_VESIZE,E_SIZE2隧道及其周围的划分TYPE,1VSWEEP,ALL软岩层的划分CMSEL,S,RYC_VESLV,SMPCHG,4,ALL清楚面单元ASEL,S,LOC,Z,0ESEL,S,TYPE,4ASLE,RACLEAR,ALLLATT,5,3,3$LMESH,ALLCMSEL,S,PH_LLESIZE,ALL,NPHAATT,4,2,2/AUTO,1NUMMRG,NODENUMCMP,NODENUMC

18、MP,ELEM=求解分析=/SOLUNSEL,S,LOC,Y,0$D,ALL,UYNSEL,S,LOC,X,0$NSEL,A,LOC,X,W$D,ALL,UXNSEL,S,LOC,Z,0$NSEL,A,LOC,Z,D0$D,ALL,UZANTYPE,STATIC !指定为静力学求解类型NLGEOM,ON !打开大变形求解NROPT,FULL !设定全N-R求解OUTRES,ALL,ALL !输出所有项，每一步都输出NSUB,60,100 !设定子步数为6，最大不超过10载荷步1，求解自重应力场ACEL,G !施加重力加速度9.8ESEL,S,TYPE,2,3,1 !选择EKILL,ALL !杀

19、死单元ALLSEL !选择所有元素SOLVESAVE,S1,DB载荷步2，采用循环进行施工阶段分析CMSEL,S,SD_V%I%,挖去左侧隧道下部的土 CMSEL,S,MG_L%I%$CMSEL,S,PH_A%I%ESLL,SESLA,AEALIVE,ALL !激活锚杆和软喷NSLE,S !选择当前单元中所有节点DDELE,ALL,ALL !删除当前选择集中所有节点的约束这样可以实现内部节点的选择和约束ESEL,ALL !选择所有单元ESEL,S,LIVE !选择所有活的单元NSEL,INVE !反向选择D,ALL,ALL !约束当前选择集中所有节点所有自由度NS=I+1SAVE,S%NS%,DB载荷步3，第二段施工段TIME,3CMSEL,S,SD_V2CMSEL,S,MG_L2$CMSEL,S,PH_A2N