西南交大工程计算实践作业

工程计算实践第一次作业

班级：专隧1班

姓名：XXXXX

学号：xxxxxxxxxx

2017年4月

1.工程背景

选取一相似实际工程，该隧道为单洞双线隧道。区域为中低山区，地形起伏较大，山坡自然坡度约30~60度，相对高差300~860m左右，主要岩层地层岩性较单一，地表零星分布有第四系坡残积层，主要为花岗岩、花岗闪长岩、全~弱风化，全风化层厚度约为20~30m，强风化层厚3~5m，下伏弱风化层。埋深约为30m。掌子面围岩主要为强风化花岗岩、花岗闪长岩。掌子面揭露呈黄褐色，岩体破碎，节理裂隙较发育，裂缝密集且部分张开；毛开挖面能够自稳，偶有掉块现象。初期支护喷射35cm厚，C25混凝土，钢架间距60cm。每三棍在拱部1400范围内施做一环4m长小导管，间距30cm。二衬厚度60cm, C30混凝土。

2.计算模型建立

根据隧道支护结构尺寸信息，采用ANSYS建模，仅考虑二次衬砌受力，二次衬砌承担50%的荷载。

2.1计算参数

二衬C30混凝土以及V级围岩的参数见表1-1。根据《铁路隧道设计规范》，可以计算出深埋隧道围岩的垂直和水平均布荷载，如表1-2所示。对于竖向和水平的均布荷载其等效节点力分别近似取节点两相邻单元水平或垂直投影长度的一半衬砌计算宽度这一面积范围内的分布荷载总和。结构自重直接通过ANSYS 添加密度施加。

2. 2模拟步骤

2.2.1建立模型

打开ANSYS，进入前处理器，输入如下命令可初步建立模型。

图1 模型示意图

2.2.2为模型赋予属性

模型材料选择为梁单元(beam3 )，截面宽度lm，高度0.4m，弹性模量32Gpa泊松比0.2，密度2500。同时选择弹簧单元（combin14）设置其相关物理和几何参数，如图2所示

2.2.3网格划分

分别对用梁单元（BEAM3）弹簧单元（combin14）对其进行网格划分。用杆单元来模拟弹簧单元。

2.2.4施加约束及荷载

施加节点约束，分别对弹簧单元施加X,Y方向的约束，约束值设置为0，同时在

隧道底部约束住隧道水平方向的位移，约束值设置为0。:如图3所示

图3 约束设置图2.2.5施加节点荷载

节点荷载计算方法

节点荷载计算示意图如下

图4 单元ij受力图

衬砌结构承受主动垂直线性荷载和水平线性荷载，计算时将其转换为节点荷载。假设单元ij受力如图1-10所示荷载作用，则节点i和节点j的等效节点荷载可分别由如下的荷载阵列求出。

命令流如下

*SET,_Q1, 99.404e3

*SET,_Q2, 99.404e3

*SET,_E1, 24.851e3

*SET,_E2, 24.851e3

*SET,_E3, 24.851e3

*SET,_E4, 24.851e3

*GET,_E_MAX,ELEM,,NUM,MAX

*GET,_E_MIN,ELEM,,NUM,MIN

*GET,_E_COUNT,ELEM,,COUNT

*GET,_N_MAX,NODE,,NUM,MAX

*GET,_N_MIN,NODE,,NUM,MIN

*GET,_N_COUNT,NODE,,COUNT

*GET,_LOC_X_MAX,NODE,,MXLOC,X

*GET,_LOC_X_MIN,NODE,,MNLOC,X

*GET,_LOC_Y_MAX,NODE,,MXLOC,Y

*GET,_LOC_Y_MIN,NODE,,MNLOC,Y

*SET,_WIDTH,_LOC_X_MAX-_LOC_X_MIN

*SET,_HEIGHT,_LOC_Y_MAX-_LOC_Y_MIN

*SET,_CENT_X,(_LOC_X_MAX+_LOC_X_MIN)/2

*SET,_CENT_Y,(_LOC_Y_MAX+_LOC_Y_MIN)/2

*SET,_I_ND,_N_MIN

*SET,_I_EL,_E_MIN

*DIM,_NF_INFO,ARRAY,_N_COUNT,4

*DIM,_ND_INDEX,ARRAY,_N_MAX

*DO,I,1,_N_COUNT,1

*IF,NX(_I_ND),EQ,_LOC_X_MAX,THEN

_CENT_Y=NY(_I_ND)

*ENDIF

*IF,NY(_I_ND),EQ,_LOC_Y_MAX,THEN

_CENT_X=NX(_I_ND)

*ENDIF

_NF_INFO(I,1)=_I_ND

_ND_INDEX(_I_ND)=I

_I_ND=NDNEXT(_I_ND)

*ENDDO

*DO,I,1,_E_COUNT,1

_NI=NELEM(_I_EL,1)

_NJ=NELEM(_I_EL,2)

_LOC_XI=NX(_NI)-_LOC_X_MIN

_LOC_YI=NY(_NI)-_LOC_Y_MIN

_LOC_XJ=NX(_NJ)-_LOC_X_MIN

_LOC_YJ=NY(_NJ)-_LOC_Y_MIN

*GET,_E_CENT_X,ELEM,_I_EL,CENT,X

*GET,_E_CENT_Y,ELEM,_I_EL,CENT,Y

*IF,_E_CENT_X,GT,_CENT_X,THEN

*IF,_E_CENT_Y,GT,_CENT_Y,THEN

_E_I=-_E3-(_E4-_E3)*(1-_LOC_YI/_HEIGHT)

_Q_I=-_Q1-(_Q2-_Q1)*_LOC_XI/_WIDTH

_E_J=-_E3-(_E4-_E3)*(1-_LOC_YJ/_HEIGHT)

_Q_J=-_Q1-(_Q2-_Q1)*_LOC_XJ/_WIDTH *ELSE

_E_I=-_E3-(_E4-_E3)*(1-_LOC_YI/_HEIGHT)

_Q_I=0

_E_J=-_E3-(_E4-_E3)*(1-_LOC_YJ/_HEIGHT)

_Q_J=0

*ENDIF

*ELSE

*IF,_E_CENT_Y,GT,_CENT_Y,THEN

_E_I= _E1+(_E2-_E1)*(1-_LOC_YI/_HEIGHT)

_Q_I=-_Q1-(_Q2-_Q1)*_LOC_XI/_WIDTH