孫　奇，王　丙，李聚軒，王　洋，吉磊杰

(1.軍事交通學(xué)院 國防交通系，天津300161； 2.75660部隊，廣西 桂林 541002；3．62101部隊，武漢 430010)

?

● 基礎(chǔ)科學(xué)與技術(shù)Basic Science & Technology

單個移動荷載作用下浮橋的動態(tài)響應(yīng)

孫奇1，王丙1，李聚軒1，王洋2，吉磊杰3

(1.軍事交通學(xué)院 國防交通系，天津300161； 2.75660部隊，廣西 桂林 541002；3．62101部隊，武漢 430010)

為研究浮橋在單個移動荷載作用下的動態(tài)響應(yīng)，將浮橋的連接接頭進行線性處理并簡化為彈性地基梁，將平均剛度法求得的平均剛度作為浮橋的整體剛度，并考慮荷載偏心距離的影響，應(yīng)用有限元方法對單個移動荷載作用下的動態(tài)響應(yīng)進行了計算分析。結(jié)果表明：荷載的移動速度和偏心距離對浮橋的運動和受力均有不同程度的影響，當(dāng)車輛通過浮橋時，應(yīng)控制車輛的移動速度和偏心距離，減小浮橋的變形和受力，以保證車輛能夠安全通過。

浮橋；動態(tài)響應(yīng)；單個移動荷載；移動速度；偏心距離

當(dāng)浮橋的連接接頭采用線性處理時，即可應(yīng)用梁模型來對浮橋進行簡化，浮橋下沉產(chǎn)生的浮力可用彈簧模擬，于是浮橋的力學(xué)模型便可采用彈性地基梁，在土木工程和水利工程中彈性地基梁有著相當(dāng)廣泛的應(yīng)用。

到目前為止，已經(jīng)有很多學(xué)者提出了不同的解析方法和數(shù)值解法來對彈性地基梁進行計算。Lee等[1]對彈性地基梁的靜力、動力及穩(wěn)定問題進行了分析。Gutierrez等[2]和Laura等[3]用簡單多項式的組合來表示梁的模態(tài)，對非均勻單跨梁在各種邊界條件下的動態(tài)響應(yīng)進行了研究。Li等[4]采用雙層彈簧阻尼單元來模擬彈性地基，對彈性地基上Timoshenko梁在移動荷載作用下的動力響應(yīng)進行了分析。Khavassefat等[5]構(gòu)建一種準靜力計算準則來計算黏彈性基礎(chǔ)上柔性結(jié)構(gòu)在移動荷載作用下的應(yīng)力和應(yīng)變。文獻[6—8]對梁在移動荷載作用下的動力響應(yīng)進行了較深入的研究。

與普通梁相比，彈性地基梁是無窮多次超靜定結(jié)構(gòu)，計算彈性地基梁時必須同時考慮地基的變形。本文在對移動荷載作用下浮橋的動態(tài)響應(yīng)進行研究時采用Winkler地基模型，Winkler地基模型假設(shè)作用在地基表面任一點單位面積上的壓力與該點的沉陷量成正比。

1　控制方程

根據(jù)有限元法的標準過程，即可得到浮橋系統(tǒng)的動態(tài)響應(yīng)控制方程：

MD"(t)+CD'(t)+KD(t)=Fo(t)

(1)

式中：M為系統(tǒng)質(zhì)量矩陣；C為阻尼矩陣；K為剛度矩陣；Fo(t)為外荷載向量；D(t)為系統(tǒng)的節(jié)點自由度向量。

式(1)中的Fo(t)包含了表面力、體積力和集中力，均是時間的函數(shù)，通常體積力與浮體的靜態(tài)浮力相平衡，因此只需考慮式(1)中與結(jié)構(gòu)體積力不平衡的荷載部分。本文應(yīng)用有限元軟件進行計算分析。

2　計算模型

浮橋整體由15只全形舟組成，舟兩端連接處的標號為1—16，每個全形舟的長為L=6．7 m，不考慮岸邊舟和跳板對浮橋變形的影響。浮橋全長Ltotal=100．5 m，寬度B=8．082 m，型深D=1．07 m，水密度ρ=1 000 kg/m3，重力加速度g=9．81 m/s2，彈性地基的地基系數(shù)K=ρ×g×B，彈性模量E=2．1 GPa，浮橋的截面為規(guī)則的矩形，兩端簡支。浮橋的平面圖和正面圖如圖1、2所示，本文將浮橋簡化為連接接頭線性處理的彈性地基梁進行計算，兩端簡支。

圖1　浮橋平面

圖2　浮橋正面

由于舟與舟之間的接頭間隙會大大降低浮橋的整體抗彎剛度，因此將單個舟的抗彎剛度代入彈性地基梁運動方程進行計算與實際不符，可用平均剛度法進行求解。當(dāng)一常量力P0作用在浮橋中點時，可以求得浮橋在靜載作用下的最大位移ymax，從而反算得到浮橋整體的抗彎剛度[9]:

(2)

圖3　靜載作用下連接接頭線性處理浮橋的位移比較

3　數(shù)值計算結(jié)果分析

本文對單個移動荷載作用下連接接頭線性處理的浮橋的動態(tài)響應(yīng)進行計算，并分析浮橋的受力情況，考慮單個移動荷載作用時偏心距離對浮橋動態(tài)響應(yīng)的影響，在所有工況中均假定荷載為勻速直線運動且始終與浮橋表面相接觸。在應(yīng)用軟件進行計算過程中，將動荷載簡化為靜荷載并將浮橋所受荷載轉(zhuǎn)化為節(jié)點荷載，當(dāng)節(jié)點荷載以不同速度在浮橋上移動時，通過控制節(jié)點荷載作用的時間間隔，完成不同速度的荷載通過浮橋時的動態(tài)響應(yīng)計算。

圖4(a)—(d)分別為無偏心的單個移動荷載以不同速度移動至浮橋跨中時浮橋的垂向變形圖、縱向轉(zhuǎn)角變形圖、應(yīng)力最大值沿橋長方向的分布圖(指各截面中的最大應(yīng)力，其所在點的位置不盡相同，故未在圖中標出，下同)和跨中節(jié)點的垂向位移時程圖。由圖4可知，荷載作用點附近浮橋的垂向位移、縱向轉(zhuǎn)角和所受應(yīng)力最大，浮橋的垂向變形、縱向轉(zhuǎn)角和所受應(yīng)力隨著荷載移動速度的增加而增加；當(dāng)荷載在浮橋上移動的過程中，浮橋的位移變形波向兩側(cè)傳播，荷載移動速度越高，波動越劇烈，同時由于波動使得浮橋的動力響應(yīng)關(guān)于荷載不對稱；當(dāng)荷載移動速度較高時，浮橋的縱向轉(zhuǎn)角變形比較大，荷載在浮橋上移動的過程中會出現(xiàn)較為明顯的“爬坡效應(yīng)”[13]，因此需要控制荷載在浮橋上的移動速度以減小浮橋的變形和受力。

(a)垂向變形

(b)縱向轉(zhuǎn)角變形

(c)最大應(yīng)力分布

(d)跨中節(jié)點垂向位移圖4　單個移動荷載無偏心時浮橋的動力響應(yīng)

圖5(a)—(d)分別為偏心0.5 m的單個移動荷載以不同速度移動至浮橋跨中時浮橋的垂向變形圖、縱向轉(zhuǎn)角變形圖、應(yīng)力最大值沿橋長方向的分布圖和橫向轉(zhuǎn)角變形圖。由圖5可知，荷載作用點附近浮橋的垂向位移、縱向轉(zhuǎn)角和橫向轉(zhuǎn)角最大，浮橋的垂向變形、縱向轉(zhuǎn)角和橫向轉(zhuǎn)角總體上隨著荷載移動速度的增加而增加，荷載移動速度越快，浮橋波動越劇烈；浮橋各截面的最大應(yīng)力隨著荷載移動速度的增加而減小，這可能與浮橋自身的振動特性有關(guān)，由于浮橋兩端簡支，荷載的偏心影響使得浮橋兩端約束處承受較大的應(yīng)力，而荷載移動速度v=60 km/h時的橫向轉(zhuǎn)角小于v=50 km/h時的橫向轉(zhuǎn)角，這可能是由于浮橋橫向轉(zhuǎn)動的相位變化所導(dǎo)致。

(a)垂向變形

(b)縱向轉(zhuǎn)角變形

(c)最大應(yīng)力分布

(d)橫向轉(zhuǎn)角變形圖5　單個移動荷載偏心0.5 m時浮橋的動力響應(yīng)

圖6(a)—(d)分別為偏心1.0 m的單個移動荷載以不同速度移動至浮橋跨中時浮橋的垂向變形圖、縱向轉(zhuǎn)角變形圖、應(yīng)力最大值沿橋長方向的分布圖和橫向轉(zhuǎn)角變形圖。由圖6可知，偏心1.0 m的單個移動荷載作用下浮橋的縱向轉(zhuǎn)角、各截面的最大應(yīng)力和橫向轉(zhuǎn)角的變化規(guī)律與偏心0.5 m的單個移動荷載作用時的情況相同，但浮橋的垂向變形不隨著荷載移動速度的變化而單調(diào)變化，這可能與浮橋自身的振動特性有關(guān)。

(a)垂向變形

(b)縱向轉(zhuǎn)角變形

(c)最大應(yīng)力分布

(d)橫向轉(zhuǎn)角變形圖6　單個移動荷載偏心1.0 m時浮橋的動力響應(yīng)

由圖4—6可知，當(dāng)荷載的移動速度v=30～60 km/h，荷載移動至浮橋跨中時，浮橋的垂向變形、各截面的最大應(yīng)力和橫向轉(zhuǎn)角隨著荷載偏心距離的增加而增加，荷載的偏心距離對浮橋縱向轉(zhuǎn)角的影響不大，浮橋在無偏心的單個移動荷載作用下各截面最大應(yīng)力隨著荷載移動速度的增加而增加，與有偏心的單個移動荷載作用時的變化規(guī)律相反，且由于荷載偏心的影響，浮橋兩端的受力明顯大于無偏心荷載作用時的情況。因此，當(dāng)車輛、坦克等荷載在浮橋上移動時應(yīng)盡量減小荷載的偏心距離以使浮橋承受較小的變形和受力。

4　結(jié)　語

本文將連接接頭線性處理的浮橋簡化為彈性地基梁，應(yīng)用有限元方法對單個移動荷載作用下浮橋的動態(tài)響應(yīng)進行了計算分析，得到了浮橋運動和受力的變化規(guī)律與荷載移動速度和偏心距離之間的關(guān)系，為浮橋以后的研究設(shè)計提供了一定的技術(shù)支撐，同時也為下一步研究單列或雙列移動荷載對浮橋動態(tài)響應(yīng)的影響打下了基礎(chǔ)。

[1]LEE S Y, KOU Y H, LIN F Y. Stability of a Timoshenko beam resting on a Winkler foundation[J]. Journal of Sound and Vibration, 1992, 153(2): 193-202.

[2]GUTIERREZ R H, LAURA P A A. Transverse vibrations of beams traversed by point masses: a general approximate solution[J]. Journal of Sound and Vibration, 1996, 195(2):353-358.

[3]GUTIERREZ R H, LAURA P A A. Vibrations of a beam of non-uniform cross-section traversed by a time varying concentrated force[J]. Journal of Sound and Vibration, 1997, 207(3):419-425.

[4]LI K F, LIU W N, MARKINE V, et al. Analytical study on the dynamic displacement response of a curved track subjected to moving loads[J]. Applied Physics and Engineering, 2013, 14(12):867-879.

[5]KHAVASSEFAT P, JELAGIN D, BIRGISSON B. A computational framework for viscoelastic analysis of flexible pavements under moving loads[J]. Materials and Structures, 2012,45(11):1655-1671.

[6]周孫基,張?zhí)焐?魏建友.長期超限移動荷載作用下T梁承載能力試驗研究[J].施工技術(shù),2013,42(23):30-33.

[7]孫明杰,許燕,張沫.移動荷載對橋梁結(jié)構(gòu)位移的影響分析[J]. 建筑技術(shù)開發(fā),2013,40(3):17-19.

[8]ZhAO Y H, ZHANG J Y, LI T, et al. Aerodynamic performances and vehicle dynamic response of high-speed trains passing each other[J]. Journal of Modern Transportation, 2012,20(1):36-43.

[9]RANJIT A T. Analytical expressions for the transfer matrix to standardize scalp potentials to infinity reference[J]. Journal of Computational Neuroscience, 2011, 31(3): 609-613.

[10]江召兵.快速重載作用下長浮橋多體系統(tǒng)波動堆積效應(yīng)研究[D].南京：解放軍理工大學(xué),2008.

[11]付世曉.錨泊彈性浮體及通載浮橋的非線性水彈性響應(yīng)研究[D].上海:上海交通大學(xué),2005.

[12]陳徐均,林鑄明,吳廣懷,等.通載浮橋動態(tài)響應(yīng)位移的測試方法與數(shù)值分析[J].振動測試與診斷,2006,26(2):97-101.

[13]TSAI C Y, HUANG S H. Transfer matrix for rotor coupler with parallel misalignment[J]. Journal of Mechanical Science and Technology, 2009,23(5):1383-1395.

(編輯：張峰)

Dynamic Response of Floating Bridge To Single Moving Load

SUN Qi1, WANG Bing1, LI Juxuan1, WANG Yang2， JI Leijie3

(1.National Defense Traffic Department, Military Transportation University, Tianjin 300161, China；2.Unit 75660, Guilin 541002, China； 3.Unit 62101, Wuhan 430010, China)

To study the dynamic response of a floating bridge under the action of a single moving load, the joints of the floating bridge are processed linearly and floating bridge is simplified into elastic foundation beam, the average rigidity method is used to obtain the average rigidity instead of the total rigidity of floating bridge, and the dynamic response of floating bridge to a single moving load is calculated and analyzed by finite element method according to the eccentric distance of the load. The results show that the speed and the eccentric distance of the load have major impact on the motion and stress of floating bridge. Therefore, the speed and the eccentric distance of vehicle should be controlled to decrease the deformation and the stress of floating bridge so as to ensure the safe pass of vehicles.

floating bridge; dynamic response; single moving load; moving speed; eccentric distance

2015-10-21；

2015-12-30．

天津市科技計劃項目(14ZCZDSF00024).

孫奇(1984—)，男，碩士，助教．

10．16807/j.cnki.12-1372/e.2016.05.020

U448．19

A

1674-2192(2016)05- 0080- 05