寶蘭客運專線路堤段影響車致地面振動的土體參數敏感度分析及反演

王東

（中鐵第一勘察設計院集團有限公司，西安 710043）

列車引起的地面振動是復雜非線性問題，而土體參數的敏感度對地面振動數值分析具有重要作用。參數敏感度是指參數發生微小變化時對系統響應的擾動程度，對實際工程有指導意義。車輛-軌道-路基-地基耦合系統響應與參數之間的關系無法通過解析的方式呈現，因此利用智能方法逼近這種復雜的非線性關系。

神經網絡的特性是選取唯一的輸入源對應唯一的輸出變量。利用神經網絡求解列車引起的地面振動問題。當土體參數一定時，輸出的響應也是唯一的。同時，神經網絡能夠逼近任何非線性函數，準確地描述多輸入與多輸出的問題。

以往的研究多運用控制變量的思想，采用數值分析或現場試驗分析單一參數的變化對地面振動響應的影響［1-3］。如王福彤等［4］通過現場試驗研究乘客滿載率及衰減距離對地面振動水平的影響，羅晶［5］通過數值模型分析了蘭州地區黃土含水率、路堤高度、車速等因素對高速鐵路引起地面振動的影響。

這些研究僅分析了某些單一參數的變化對地面振動響應水平的影響，還缺乏各項參數之間對地面振動響應的影響程度。因此，有必要對列車致地面振動過程中土體各參數的敏感度進行分析。本文結合有限元方法和BP神經網絡，依據現場試驗工況建立數值模型，結合敏感度理論對土體各參數敏感度進行分析。在參數反演時就能夠選擇出較敏感參數，也能對參數較敏感的測點進行測值，從而提高反演的效率和精度。同時在減振隔振措施中，可針對性地對高敏感度參數進行處理，為控制地面振動提供指導。

圖1 測點布置示意（單位：m）

1 依據現場試驗建立數值分析模型

考察地面振動的現場測試位于寶蘭客運專線榆中站附近，軌道型式為CRTS型雙塊式無砟軌道，線路兩側地勢平坦開闊，便于布設測點（如圖1所示）。地基土層由素填土和粉狀黃土組成。試驗儀器為東方振動研究所INV3062C型24位采集儀，配套加速度傳感器，采樣頻率為1 024 Hz?？紤]到地面振動評價多以垂向響應為指標，在垂向設置加速度傳感器。列車車速利用記錄頭車通過測段的時間求得，現場對測試區段通過的10列動車引起的地面振動進行測試。本文選取速度250 km/h的CRH2G型8節編組動車組測試數據進行分析。

根據現場試驗工況，利用有限元軟件ABAQUS建立了車輛-軌道-路基-地基非線性耦合系統的三維數值模型，車輛為8節編組的整車模型，各部件采用剛體模擬，一系、二系懸掛和車廂連接采用彈簧阻尼系數單元模擬?？紤]車輛-軌道-路基-地基系統的耦合動力作用，在車輛輪對與鋼軌之間施加接觸對模擬輪軌動力相互作用，扣件簡化為彈簧-阻尼單元用以實現鋼軌與軌枕的相互作用，其余各層間假設位移連續，采用綁定約束。為了提高計算效率，采用ABAQUS/Explicit顯式動力學求解器，并通過并行技術加速求解。

1.1 模型參數

車輛部分采用CRH2G型動車組模型，車輛系統編組質量為359.7 t，構架質量為2 760 kg，輪對質量為2 000 kg，車輪半徑為860 mm，固定軸距2 500 mm，編組長度為201.4 m，一系懸掛剛度為1.176×106N/m，一系懸掛阻尼系數1.96×104N·s/m，二系懸掛剛度為0.221×106N/m，二系懸掛阻尼系數4.00×104N·s/m。

軌道參數采用CRTS型雙塊式無砟軌道真實值，見表1。

表1 軌道參數

以中國高速鐵路無砟軌道不平順譜［6］的高低及水平不平順作為激勵輸入軌道模型，土體采用Drucker?Prager彈塑性本構模型。采用黏彈性人工邊界［7］消除有限尺寸區域的邊界效應。地基模型材料參數根據該處地勘資料取值，見表2。地基模型沿線路方向長450 m，深15 m，沿線路中心線兩側均寬35 m，車輛-軌道-路基-地基模型如圖2所示。

表2 地基模型材料參數

1.2 數值結果驗證

將實測參數帶入模型進行計算，得到測點處地面垂向振動的時程響應，并由時程響應通過式（1）得到加速度有效值arms。

圖2 有限元模型及測點布置

式中：n為采樣點的個數；ai為采樣點加速度。

地面各測點垂向振動加速度計算值與實測值對比，見圖3。

圖3 各測點垂向振動加速度計算值與實測值對比

由圖3可知，各測點處地面垂向振動加速度有效值的計算值與實測值均呈衰減趨勢，距離10～20 m時下降得較快，距離20～35 m時下降速度減慢，且兩者吻合得很好，驗證了數值模型的合理性與計算的正確性，建立的車輛-軌道-路基-地基系統數值模型可適應地面振動分析。

2 建立神經網絡

2.1 列車致地面振動參數集及取值范圍

在車輛-軌道-路基-地基耦合系統中，影響地面振動的地基參數主要有密度ρ、彈性模量E、阻尼系數μ、泊松比ξ，內摩擦角c、黏聚力φ。以垂向加速度有效值arms為指標衡量地面振動水平，則系統中其他條件一定時，距線路中心線某處地面振動垂向加速度有效值可表示為與參數集相關的多元函數，即

由測段地勘資料可得式（2）中所需土體參數集的取值范圍，見表3。

表3 土體參數取值范圍

2.2 獲取學習樣本

神經網絡學習樣本是將每個參數在其取值范圍內均作密集離散，利用正交設計原理確定參數集樣本，樣本數量越多，神經網絡逼近非線性問題的能力也越強，但網絡學習收斂速度成倍減慢。數值分析時，獲取樣本的計算代價也隨之增大，對于車輛-軌道-路基-地基這種復雜耦合系統顯然是不經濟的，因此本文采用拉丁超立方抽樣（Latin Hypercube Sampling，LHS）方法對離散后的參數集進行抽樣。LHS方法是一種均勻空間抽樣方法［8］，具有抽樣次數較少、可有效避免重復抽樣等優點，適用于對復雜的多維參數空間抽樣。

LHS方法可按如下步驟實施［9］：

1）定義有p個參數構成的參數集｛x1，x2，…，xp｝，將每個參數的取值范圍等概率地分為q個互不重疊的區間。

2）將每個參數xi在各個區間內隨機抽取1個樣本代表xiq，即每個參數xi含q個樣本值。

3）將每個參數xi的q個樣本值隨機排列，共隨機生成p個排列。

4）每次從p個排列中順序提取各參數的1個樣本值，即形成1個樣本參數集Xq；依次提取q次，即得q個樣本參數集｛x1，x2，…，xq｝。

本文通過LHS方法得到50組參數集，將參數集的各參數代入數值模型進行計算，得到系統的輸出集，從而獲取神經網絡的輸入-輸出學習樣本。

2.3 神經網絡結構的建立

采用MATLAB平臺編制神經網絡程序，將參數集樣本作為輸入向量，各測點振動加速度有效值作為輸出向量，利用模塊化網絡結構，分別建立參數集在6個測點處垂向振動加速度有效值的映射關系。設計6個子神經網絡模塊，如圖4所示。輸入層含6個神經元，對應6個輸入參數，輸出層含1個神經元，對應該測點處加速度有效值。隱含層神經元數量采用公式（3）計算［10］。

式中：t為輸入層神經元數；s為輸出層神經元數；a為1～10的常數；n1為隱含層神經元數。

在本文網絡結構中，t=6，s=1，依據式（3）計算得該網絡模型的隱含層神經元數應在4～13。分別選取不同隱含層神經元數進行網絡訓練，可得n1=10時系統誤差最小，訓練效果最為理想。每個模塊都是6?10?1的BP神經網絡，將學習樣本學習至收斂，得到訓練成熟的神經網絡模型。

圖4 單隱含層的BP神經網絡模型

3 參數敏感度分析

借鑒文獻［11］，建立S型函數的BP網絡，其參數敏感度可定義為網絡輸入對輸出的偏導數，可由式（4）得出。

式中：Ok為輸出層單元k的輸出值；Oi為輸入層單元i的輸入值；Wjnk為中間層單元jn與輸出層單元k的連接權；Wjn-1jn為2個中間層單元jn-1與jn的連接權；Wij1為輸入層單元i與中間層單元j1的連接權；G(ek)為輸出層單元k與中間層單元jn的影響函數；G(ejn)為中間層單元jn與jn-1的影響函數；G(ej1)為中間層單元j1與輸入層單元i的影響函數。

由式（4）可得各測點參數敏感度，見圖5。

圖5 各測點參數敏感度

由圖5可見，各測點處影響車致地面振動的地基土質參數敏感度排序為彈性模量＞阻尼系數＞泊松比＞密度＞內摩擦角＞黏聚力。10～30 m處彈性模量、阻尼系數、泊松比、密度的敏感度顯著大于內摩擦角和黏聚力。這是由于鐵路交通引起的地面振動微小，在線彈性階段［12］黏聚力、內摩擦角屬于材料的塑性參數，故其對地面振動的敏感度低。因此可以認為彈性模量、阻尼系數、泊松比和密度為主要力學參數，而內摩擦角、黏聚力為次要力學參數。隨著距離的增大，彈性模量、泊松比、密度、內摩擦角、黏聚力的敏感度隨距離變化不大，而阻尼系數的敏感度隨距離增大有上升趨勢。這是因為阻尼對地面振動的影響為過程累積性因素，隨振動傳播其對振動的影響越來越凸顯。

4 參數反演

以上確定了影響地面振動水平的主要力學參數和次要力學參數。根據以上結果，對不同的參數分別進行敏感度評估。由于土體的內摩擦角、黏聚力對地面振動影響很小，將各測點處地面振動響應作為輸入向量，土體力學參數（密度、彈性模量、阻尼系數、泊松比）作為輸出向量，建立BP神經網絡模型，經過若干次迭代達到收斂。

將現場實測的各測點處地面振動響應輸入到訓練好的神經網絡中，得到參數反演值及其與實際值的相對誤差，見表4。

表4 參數反演值與實際值相對誤差

由表4可知，與地勘資料中現場實際參數相比，基于神經網絡的參數反演得到參數值相對誤差均較?。?%以下），基本控制在工程可接受的范圍內，證明了神經網絡在通過高速鐵路引起的地面振動響應來反演地基參數的可行性。

5 結論

本文依據現場試驗建立了車輛-軌道-路基-地基非線性耦合系統的三維數值模型，結合現場試驗數據，驗證了數值模型，建立了反映土質參數和測點地面振動加速度有效值之間映射關系的神經網絡，基于敏感度分析理論對地基土質參數進行了敏感度分析與反演，得出以下結論：

1）數值模型采集到地面垂向振動加速度響應與現場試驗真實數值吻合得很好，證明了建模方法及有限元模型的合理性。

2）影響車致地面振動的地基土質參數敏感度的順序依次為彈性模量＞阻尼系數＞泊松比＞密度＞內摩擦角＞黏聚力；其中內摩擦角、黏聚力的敏感度低，因此在解決地面振動問題時僅能用作次要力學參數，其余4項參數可作為主要力學參數；隨著距離的增大，阻尼系數敏感度有上升趨勢，其余5項參數變化不大。

3）由于內摩擦角、黏聚力敏感度低，在反演計算時僅選取彈性模量、阻尼系數、泊松比、密度作為輸出向量，得到反演土質參數與現場實測土質參數相對誤差均小于5%，在工程可接受的范圍內。因此通過神經網絡在高速鐵路引起的地面振動響應問題中來反演地基參數是可行的。