FINISH?

/clear?

/filename,tunnel_tx

/PREP7?

!總體參數(shù)

ONAME='OUT'? !節(jié)點(diǎn)內(nèi)力輸出文件名稱"***"

PI=ACOS(-1)? !返回π值

?

!二襯承受圍巖壓力比例

?

hzblq=1!二襯承受圍巖壓力比例(明洞)

?

?

!計(jì)算需輸入的參數(shù)

?

cqhd=0.8!二襯厚度

Bkd=10.5 !坑道寬度(包括初支)

hkd=6.5!坑道高度(包括初支)

cqtm=30e9 !襯砌彈性模量

cqzd=25e3 !襯砌重度

dycd=0.5 !襯砌單元計(jì)算長(zhǎng)度(每個(gè)單元)

?

?

!圍巖有關(guān)參數(shù)

?

FAIc=35*PI/180 明洞周邊土體內(nèi)摩擦角

Wyzd=22e3 !明洞周邊土體容重

FAIc1=35*PI/180 回填土內(nèi)摩擦角

Wyzd1=12e3 ! 回填土容重

Hp=4.0 !洞頂回填土厚度

pstg=0 ! M10回填土底部距離邊墻底部

?

?

!明洞隧道側(cè)壓力系數(shù)λ

rmtq=0.487325 !明洞隧道周邊土體側(cè)壓力系數(shù)λ

rmtq1=0.487325 !明洞隧道回填土側(cè)壓力系數(shù)λ

?

!高速公路車輛荷載

?qq1=50e3?? !此值是均布車輛荷載和人均均布荷載組成

jzhz=0

jBhz=0E3 !此值是由集中荷載按30度擴(kuò)散角計(jì)算而得

QQQ1=qq1+jBhz

?

?!明洞隧道荷載計(jì)算(車道均布荷載換算成相應(yīng)的土體高度計(jì)算側(cè)壓力）

?qq=Wyzd1*Hp !明洞隧道回填土豎直荷載(已考慮承載比例)

?eq1=rmtq1*wyzd1*(hp+QQQ1/wyzd1)*hzblq !明洞隧道回填土水平荷載最小值(已考慮承載比例)

?eq2=rmtq1*wyzd1*(hp+hkd-pstg+QQQ1/wyzd1)*hzblq !明洞隧道回填土水平荷載最大值(已考慮承載比例)

?eq3=rmtq*wyzd*((wyzd1/wyzd)*(hp+hkd+QQQ1/wyzd1-pstg))*hzblq !明洞隧道水平荷載最小值(已考慮承載比例)

?eq4=rmtq*wyzd*(pstg+(wyzd1/wyzd)*(hp+hkd+QQQ1/wyzd1-pstg))*hzblq !明洞隧道水平荷載最大值(已考慮承載比例)

?

?

?

?

!單元,材料,實(shí)常數(shù)定義

?!單元定義

?

?ET,1,BEAM3 !1號(hào)單元為梁?jiǎn)卧?模擬襯砌

?ET,2,LINK10 !2號(hào)單元為彈簧連桿,模擬彈性抗力

?KEYOPT,2,2,0 !2號(hào)單元參數(shù)設(shè)定

?KEYOPT,2,3,1 !2號(hào)單元參數(shù)設(shè)定

?ET,3,SURF153 !3號(hào)單元為表面荷載單元,加載用

?

?!材料定義

?MP,EX,1, cqtm !第1種材料彈性模量

?MP,NUXY,1,0.2 !第1種材料泊松比

?MP,DENS,1, cqzd !第1種材料重度

?

?!第2種材料彈性模量

? MP,EX,2,dycd*100E6

?

!第3種材料彈性模量

?MP,EX,3,0 !第3種材料彈性模量

?

?!實(shí)常數(shù)定義

?

?TYPE,1

?MAT,1

?R,1 ,cqhd,cqhd**3/12,cqhd !第1種材料實(shí)常數(shù)

?TYPE,2

?r,2,dycd !第2種材料實(shí)常數(shù)

?type,3

?r,3 !第3種材料實(shí)常數(shù)

?

!***********讀入節(jié)點(diǎn)信息***********

NREAD,'jjl','txt',' '

!***********讀入節(jié)點(diǎn)信息***********

?

?

?

?

!**********************建立襯砌梁?jiǎn)卧?********

TYPE, 1??

MAT, 1

Real,1??

eSYS, 0????

*DO,I,1,63

E,I,I+1

*ENDDO

E,64,1

?

!**********************建立連桿單元*********

?

TYPE,?? 2?

MAT,?????? 2

REAL,????? 2

ESYS,?????? 0?????

?

*DO,II,1,64

E,II,II+100

*ENDDO

allsel

eplot

?

?

?

?

?

?

!創(chuàng)建單元組

?

?!創(chuàng)建襯砌單元組

?

?esel,s,type,,1

?cm,e-cq,elem

?

?!創(chuàng)建連桿單元組

?esel,s,type,,2

?cm,e-lgx,elem

?

?

!創(chuàng)建荷載單元組

? !創(chuàng)建右側(cè)荷載單元

!**********************建立連表面荷載單元*********

TYPE, 3

MAT, 3

Real,3

ESYS, 0

*DO,I,21,36

E,I,I+1

*ENDDO?

ESEL,S,TYPE,,3

eplot

cm,e-hz-y,elem

?

!創(chuàng)建左側(cè)荷載單元

!**********************建立連表面荷載單元*********

TYPE, 3

MAT, 3

Real,3

ESYS, 0

*DO,I,5,20

E,I,I+1

*ENDDO?

ESEL,S,TYPE,,3

eplot

cmsel,u,e-hz-y

cm,e-hz-z,elem

?

?!創(chuàng)建右側(cè)下方荷載單元

TYPE, 3

MAT, 3

Real,3

ESYS, 0

*DO,I,37,43

E,I,I+1

*ENDDO?

ESEL,S,TYPE,,3

eplot

cmsel,u,e-hz-y

cmsel,u,e-hz-z

cm,e-hz-yx,elem

?

?!創(chuàng)建左側(cè)下方荷載單元

TYPE,3

MAT, 3

Real,3

ESYS, 0

*DO,I,1,4

E,I,I+1

*ENDDO

E,1,64

E,64,63

E,63,62?

ESEL,S,TYPE,,3

eplot

cmsel,u,e-hz-y

cmsel,u,e-hz-z

cmsel,u,e-hz-yx

cm,e-hz-zx,elem

?

?

?

?!創(chuàng)建上部豎向荷載單元

TYPE, 3

MAT, 3

Real,3

ESYS, 0

*DO,I,5,20

E,I,I+1

*ENDDO?

*DO,I,21,36

E,I,I+1

*ENDDO?

ESEL,S,TYPE,,3

eplot

cmsel,u,e-hz-y

cmsel,u,e-hz-z

cmsel,u,e-hz-yx

cmsel,u,e-hz-zx

cm,e-hz-sc,elem

?

?

!創(chuàng)建上部豎向集中車輛擴(kuò)散影響荷載單元

TYPE, 3

MAT, 3

Real,3

ESYS, 0

*DO,I,5,20

E,I,I+1

*ENDDO?

*DO,I,21,36

E,I,I+1

*ENDDO?

ESEL,S,TYPE,,3

eplot

cmsel,u,e-hz-y

cmsel,u,e-hz-z

cmsel,u,e-hz-yx

cmsel,u,e-hz-zx

cmsel,u,e-hz-sc

cm,e-hz-sJ,elem

?

?

?

!求襯砌單元節(jié)點(diǎn)X、Y的極值

?cmsel,s,e-cq

?ALLSEL,BELOW,ELEM

?cm,n-cq,node

?*get,xmax,node,,mxloc,x

?*get,xmin,node,,mnloc,x

?*get,ymax,node,,mxloc,y

?*get,ymin,node,,mnloc,y

xmax1=nx(5)

xmax2=nx(63)

ymax1=ny(5)

ymax2=ny(63)

?

?

!將連桿外側(cè)節(jié)點(diǎn)設(shè)為固定端

?

?NSEL,ALL

?cmsel,u,n-cq

?d,all,all

allsel

eplot

?

*DIM,Xy,ARRAY,64,2?????????????? !定義2個(gè)數(shù)組

*DIM,Nodes,ARRAY,100

*GET,Nd,NODE,0,NUM,MIN??????????? !得到最小的節(jié)點(diǎn)編號(hào)??

*DO,I,1,64,1??

Nodes(I)=Nd?????????? !將節(jié)點(diǎn)列表放到數(shù)組NODES中??

Xy(I,1)=NX(Nd)????????? !取出節(jié)點(diǎn)的X坐標(biāo)??

Xy(I,2)=NY(Nd)???????? !取出節(jié)點(diǎn)的Y坐標(biāo)??

Nd=NDNEXT(Nd)???????? !讀出下一個(gè)節(jié)點(diǎn)編號(hào)??

*ENDDO??

cmsel,s,e-cq

*GET,NELM,ELEM,0,COUNT?????????? !得到所選擇的單元總數(shù)??

*DIM,Elems,ARRAY,NELM????????? !定義2個(gè)數(shù)組??

*DIM,Ndlst,ARRAY,NELM,2

*GET,El,ELEM,0,NUM,MIN??????????? !得到最小的單元編號(hào)??

*DO,I,1,NELM,1??

Elems(I)=El??????????????? !將單元編號(hào)存入到數(shù)組ELEMS中??

*DO,J,1,2??

?? Ndlst(I,J)=NELEM(El,J)???????????? !讀出單元的2個(gè)節(jié)點(diǎn)編號(hào)??

*ENDDO??

El=ELNEXT(El)?????????????????????? !讀出下一個(gè)單元編號(hào)??

*ENDDO

ALLSEL

?

?

finish

save,tunnel_model

?

/solu

ANTYPE,0? !分析類型為靜態(tài)

NLGEOM,0

ACEL,0,1,0,

?

FCUM,ADD,1,

SFCUM,PRES,ADD,1,1,

?

!加拱部圍巖垂直荷載

F,1,FY,-jzhz

SFGRAD,PRES, ,X,0,22*0.9075/5.25*1000,

cmsel,s,e-hz-sc? !選擇上部荷載單元

SFE,all,4,PRES,,qq+QQ1+JBHZ,1,-1E10

!cmsel,s,e-hz-sJ? !選擇上部荷載單元

!SFE,all,4,PRES,,JBHZ,1,-1E10

?

?

?

?

eqtd=(eq2-eq1)/(ymax-ymax1) !求淺埋隧道水平力梯度

SFGRAD,PRES, ,Y,ymax,-eqtd,

cmsel,s,e-hz-z? !選擇左邊荷載單元

sfe,all,4,pres,,eq1,1e10,1

cmsel,s,e-hz-y? !選擇右邊荷載單元

sfe,all,4,pres,,eq1,-1e10,1

?

eqtd1=(eq4-eq3)/(ymax1-ymax2) !求淺埋隧道水平力梯度

SFGRAD,PRES, ,Y,ymax1,-eqtd1,

cmsel,s,e-hz-zx? !選擇左邊下荷載單元

sfe,all,4,pres,,eq3,1e10,1

cmsel,s,e-hz-yx? !選擇右邊下荷載單元

sfe,all,4,pres,,eq3,-1e10,1

?

?

?

?

TIME,1

/STATUS,SOLU

NSUBST,10,0,0?

ALLSEL

SOLVE

finish

?

?

!后處理

/POST1

?

!創(chuàng)建結(jié)果文件

cmsel,s,e-cq

ETAB,WJ11,SMISC,6

ETAB,WJ12,SMISC,12??

ETAB,ZL11,SMISC,1

ETAB,ZL12,SMISC,7

ETAB,JL11,SMISC,2

ETAB,JL12,SMISC,8

?

!顯示結(jié)果圖形并保存圖形

plls,wj11,wj12,-1!顯示并保存彎距圖

/ui,copy,save,jpeg,full,color,reverse,portrait,yes,100

plls,zl11,zl12,-1 !顯示并保存軸力圖

/ui,copy,save,jpeg,full,color,reverse,portrait,yes,100

plls,jl11,jl12,-1 !顯示并保存剪力圖

/ui,copy,save,jpeg,full,color,reverse,portrait,yes,100

*DIM,jdnl,ARRAY,100,7

*DO,I,1,100,1??

jdnl(i,1)=ndlst(i,1)

*GET,jdnl(i,2),elem,%i%,SMISC,6 !i端軸力

*GET,jdnl(i,3),elem,%i%,SMISC,12 !j端軸力

*GET,jdnl(i,4),elem,%i%,SMISC,1 !i端彎矩

*GET,jdnl(i,5),elem,%i%,SMISC,7? !J端彎矩

*GET,jdnl(i,6),elem,%i%,SMISC,2 !i端剪力

*GET,jdnl(i,7),elem,%i%,SMISC,8 !J端剪力

*ENDDO

*cfopen,oname,sqq

*vwrite,

(//)

*vwrite

('*********************************襯砌節(jié)點(diǎn)彎矩、軸力，剪力結(jié)果輸出***********************')

*vwrite,

(//)

*vwrite,

('**************************************************************************************')

*vwrite,

('單元編號(hào)',4x'I端彎矩',4x,'J端彎矩',6x,'I端軸力',6x,'J端軸力',6x,'I端剪力',6x,'J端剪力')

*vwrite,jdnl(1,1),jdnl(1,2),jdnl(1,3),jdnl(1,4),jdnl(1,5),jdnl(1,6),jdnl(1,7)

(f4.0,2x,f12.4,1x,f14.4,1x,f14.4,1x,f14.4,1x,f14.4,1x,f14.4)

*vwrite,

('**************************************************************************************')

*cfclos