張業(yè)民,周傳璐,宋 濤
(1.沈陽大學 建筑工程學院,遼寧 沈陽 110044;2.大連交通大學 土木與安全工程學院,遼寧 大連 116028;3.聊城市公路局 養(yǎng)護中心,山東 聊城 252000)
京滬高鐵隧道長度約16km,占正線長度的1.2%.全線鋪設無砟正線約1 268km,占線路長度的96.2%.有砟軌道正線約50km,占線路長度的3.8%.京滬高鐵全線鋪設無砟軌道,無砟軌道是一種少維護的軌道結構,它利用成型的組合材料代替道砟,將輪軌力分布并傳遞到路基基礎上.與有砟軌道相比,無砟軌道具有以下優(yōu)點:良好的軌道穩(wěn)定性、連續(xù)性和平順性,使用壽命長;結構耐久性好,維修工作量少;可避免高速行駛下有砟軌道的道砟飛濺;有利于適應地形選線,減少線路的工程投資;可減輕橋梁二期恒載,降低隧道凈空;載荷特性和車輛的平穩(wěn)性都較好;維護費用比有砟軌道低等.但是,同時也存在很多缺點:初期投資大;一旦基礎變形下沉,修復困難,改進的可能性受到限制;無砟軌道不能在黏土深路塹、松軟土路堤或地震區(qū)域鋪設;無砟軌道減震性較差;無砟軌道噪聲水平比有砟軌道高5dB;混凝土無砟軌道為剛性承載層,當達到承載強度極限時將產生斷裂,并引起軌道幾何尺寸的突然變化及難以預見的惡化[1].另外,高速鐵路上的隧道不同于一般的鐵路隧道,高速鐵路的空氣動力學效應是一個非常值得關注的問題.
由于高速鐵路的以上眾多特點,高速列車通過隧道時所產生的震動對路基的影響是研究的一個重要方面.高速鐵路中,隧道多以盾構法施工,其斷面為圓形,成形后隧道頂棚以上的空間和路基板以下的空間用作通風管道、電纜管道和其他設施等.成形后的盾構隧道斷面如圖1所示.
圖1 圓形隧道斷面示意圖Fig.1 The diagram of circular tunnel section
不可避免地,在路基板上因排水及通風等工程實際需要會設有圓形孔洞,因此,本文以研究路基板上的圓形空洞所產生的動應力集中為主要目的.依此,可進一步研究路基板上矩形、裂紋等不規(guī)則形狀開口所產生的動應力集中.
研究對象是在無限均勻的平板(即路基板),SV 波在路基板傳播的示意圖如圖2所示.
圖2 SV波在路基板傳播的示意圖Fig.2 The schematic diagram of the spread of SV wave in the board
穩(wěn)態(tài)平面SV 波沿x 軸正方向入射到路基板中.入射波在xOy 平面內與時間的依賴關系為exp(-iωt),受彎板橫向位移w的控制方程為[2]
在極坐標中,彎矩和剪力可寫成
SV 波沿路基板通過圓孔時發(fā)生散射,產生的散射場[3]可描述為
略去時間因子之后,總波場可寫成
式中,Jn(·)為第一類Bessel函數.
路基板與圓孔之間的邊界條件為[4]
根據各系數間的正交性[5]可得模式系數,由下式決定:
式中
在空洞表面上,用M0=w0Dk2無量綱化[6]后的非零彎矩為
通過給式中的參數定義不同的值,利用MATLAB將式(13)進行數值計算,研究路基板SV 波[7]的散射與動應力集中.
當k 取不同值時,計算結果如圖3 和圖4所示.
圖3 不同波數下的動應力分布Fig.3 Dynamic stress distribution under different wave number
圖4 不同波數下動應力隨泊松比的變化關系Fig.4 Relationship of dynamic stress that changes with the Poisson's ratio under different wave number
從圖3和圖4可以看出,孔與孔洞邊界處的動應力隨著無量綱波數[8]k的變化而變化,但是,同時最大動應力的位置有所偏移.總體而言,在k=0.05附近時達到最大值.從圖4 可以看出,當泊松比[9]ν 取不同值時,動應力會有不同的波動,除了在k=0.05附近時有最大值外,在其他位置會有波動而出現(xiàn)峰值,如在k=1.8附近時.
本文研究了隧道中板型路基上的圓孔對SV波的散射和動應力集中,得到了問題的解析解,給出了路基中動應力集中因子的數值解,分析了路基的剪切模量的大小、波數等對動應力集中的影響.通過數值計算分析表明,板型路基的剪切模量的大小、波數等是影響動應力集中的重要因素.
數值計算結果表明,在波數為0.05時動應力集中達到最大值,并且波數在1.8附近時會出現(xiàn)峰值,這些位置在進行結構設計時應作動應力集中方面的考慮.
通過對板型路基上的圓孔對SV 波的散射和動應力集中的研究,為隧道中板式路基的高速列車震動評價和地震評價提供了理論依據.
[1]張業(yè)民,蘇居儒,王芳,等.土木工程新理論新技術[M].北京:高等教育出版社,2003:7-21.(Zhang Yemin,Su Juru,Wang Fang,et al.The New Theory of Civil Engineering New Technologies[M].Beijing:Higher Education Press,2003:7-21.)
[2]鮑亦興,毛昭宙.彈性波的衍射與動應力集中[M].北京:科學出版社,1993:160-178.(Bao Yixing,Mao Zhaozhou.Diffraction of Elastic Waves and Dynamic Stress Concentration[M].Beijing:Science Press,1993:160-178.)
[3]Zhang Yemin.Building Resoname Device[J].Recent Development of Geoteachnical and Geo-Environmental Engineering In Asia,2006,12:585-588.
[4]AndersenK H,Jacobus H P,Stephen F B,et al.Cyclic and Static Laboratory Tests on Drammen Clay[J].Journal of Geotechnical Engineering Division,1980,106(5):499-529.
[5]Zhang Yemin,Zhang Xinhua.Study of Cement-based Composite Materials on the Force Senstived Effect[C].The Third International Conference on Transportation Engineering(ICTE 2011),2011(27):97-99.
[6]陸全康,趙慧芬.數學物理方法[M].北京:高等教育出版社,2003:398-430.(Lu Quankang,Zhao Huifen.Methods of Mathematical Physics[M].Beijing:Higher Education Press,2003:398-430.)
[7]Waterman P C.Matrix Theory of Elastic Waves Cattering[J].Journal of the Acoustical Society of America,1976,60(3):567-578.
[8]Rice J M,Sadd M H.Propagation and Scattering of SH Wavesin Semi-infinite Domains using a Time-dependent Boundary Element Method[J].ASME Journal of Applied Mechanics,1984,51(5):641-645.
[9]楊賀娟,馬驥.基于MATLAB的直流調速系統(tǒng)參數研究[J].沈陽大學學報,2011,23(6):115-117.(Yang Hejuan,Ma Ji.Parameters Research of Direct Current Governor System based on MATLAB[J].Journal of Shenyang University,2011,23(6):115-117.)