高 鵬，黃珊珊，歐陽峰，陳 達

（河海大學港口海岸與近海工程學院，南京 210098）

寬支座連續(xù)梁彎矩簡化計算方法研究

高 鵬，黃珊珊，歐陽峰，陳 達

（河海大學港口海岸與近海工程學院，南京 210098）

高樁碼頭縱橫梁正負彎矩的最大值是工程設計中的控制因素，文章在分布彈性支承連續(xù)梁模型的基礎上，提出寬支座連續(xù)梁跨中和支座的正負彎矩最大值的簡化計算公式。對于不同懸跨比的連續(xù)梁，通過ANSYS對不同的支座相對寬度和相對剛性系數(shù)進行有限元模擬，根據(jù)數(shù)值分析結(jié)果確定公式中的系數(shù)。對某一高樁碼頭工程縱梁，將簡化公式和有限元的計算結(jié)果進行對比，兩者誤差小于5%，說明簡化公式具有一定參考意義。

寬支座；縱橫梁；彎矩最大值；簡化計算；ANSYS

高樁碼頭是普遍應用的碼頭結(jié)構(gòu)型式之一。隨著裝卸設備的大型化以及碼頭向外海的發(fā)展，縱橫梁節(jié)點處基樁的數(shù)量和樁徑越來越大，樁帽也相應的越來越大，有些碼頭的縱橫梁支座寬度甚至達到梁跨度的0.6倍左右［1］，于是就產(chǎn)生了支座寬度過大情況下梁的內(nèi)力應該如何計算的新問題。

高樁碼頭一般采用平面計算方法，對于支座寬度不大的梁，以往直接將支座視作點支承進行內(nèi)力計算［2］。對于支座寬度較大的梁，2010版的《高樁碼頭設計與施工規(guī)范》［3］提出當支座相對寬度和相對剛性系數(shù)滿足表1條件時，將連續(xù)梁結(jié)構(gòu)簡化為彈性點支承連續(xù)梁計算，對支座處的負彎矩和剪力進行削峰處理，但未考慮寬支座對跨中正彎矩的影響；當不滿足表1條件時，寬支座對連續(xù)梁的支承作用可以用支座寬度范圍內(nèi)近似均勻分布的彈簧來替代，類似于彈性地基，支座反力按文克爾地基模型考慮，這種假定接近于真實情況，能得到具有一定精度的計算結(jié)果，但計算過程較復雜。

另外，根據(jù)梁與支座連接處剛度的大小，有學者提出了節(jié)點域剛性連續(xù)梁法［4］和節(jié)點域彈性連續(xù)梁法［5］，以兩端帶節(jié)點域桿梁為計算單元；還有學者提出一種“寬支承方式”［6］，以支座范圍內(nèi)的梁段和跨中梁段為計算單元，這三種方法均是通過推導單元剛度矩陣，并借助計算機計算梁的內(nèi)力，計算過程繁瑣且合理性需做進一步探討。

連續(xù)梁彎矩的最大值是結(jié)構(gòu)設計的主要控制條件，為了尋求更加簡便和實用的彎矩最大值的計算方法，對于寬支座結(jié)構(gòu)，本文基于分布彈性支承連續(xù)梁模型，提出均布荷載作用下跨中正彎矩和支座負彎矩最大值的簡化計算公式，并基于數(shù)學模型分析確定公式參數(shù)，經(jīng)與實際工程案例有限元計算結(jié)果對比，驗證簡化計算公式的可靠性。

表1 按點支承的連續(xù)梁計算的條件Tab.1 Conditions of calculation based on point supporting continuous beam

1 計算模型

根據(jù)2010版的《高樁碼頭設計與施工規(guī)范》，對于寬支座高樁碼頭縱橫梁，將其簡化為帶有懸臂的分布彈性支承連續(xù)梁，如圖1所示，l為跨長；m為懸跨比，即懸臂長與跨長之比；β為支座寬度與跨長之比；K為支座的剛度系數(shù)；E和I分別為連續(xù)梁的彈性模量和截面慣性距。若在支座范圍以有限個等間距布置的彈簧代替，則單個彈簧的剛度系數(shù)K0=K/n，n為支座范圍內(nèi)的彈簧個數(shù)。

文獻［7］的研究表明彈性支承連續(xù)梁的跨數(shù)在5跨及以上時，中間幾跨（跨3、4…）之間、中間支座（支座3、4…）之間的彎矩相差很小。因此，為簡化研究，本文基于5跨連續(xù)梁模型（圖2）進行簡化計算公式推導，所得結(jié)論同樣適用于5跨以上連續(xù)梁。圖2所示的超靜定結(jié)構(gòu)，連續(xù)梁的內(nèi)力與其本身的絕對剛度無關，而與梁體和樁的相對剛度有關，可用無量綱的相對剛性系數(shù)α=6EI/Kl3［8］表示。因此，該連續(xù)梁模型的內(nèi)力大小與懸跨比m、支座相對寬度β、相對剛性系數(shù)α以及荷載q有關。

圖1 分布彈性支承連續(xù)梁模型Fig.1 Continuous beam with uniformly distribution elastic supporting

圖2 5跨分布彈性支承連續(xù)梁模型Fig.2 5-span continuous beam with uniformly distribution elastic supporting

2 彎矩最大值的簡化公式

在進行連續(xù)梁的斷面和配筋設計時，正負彎矩最大值是主要控制條件（下拉為正，上拉為負），而均布荷載是控制縱橫梁內(nèi)力的主要因素，因此本文主要針對均布荷載作用下連續(xù)梁的正負彎矩最大值開展研究。懸臂長度對均布荷載下連續(xù)梁彎矩最大值出現(xiàn)的位置有重要的影響，在實際高樁碼頭縱橫梁結(jié)構(gòu)中，懸臂長度通常都較?。?］，懸臂部分荷載產(chǎn)生的負彎矩及其對其他跨的彎矩值影響均較小，負彎矩的最大值出現(xiàn)在中間的某支座，正彎矩的最大值出現(xiàn)在邊跨跨中（圖3）。

圖3 彎矩最大值出現(xiàn)的位置Fig.3 Position of the maximum bending moment

支座寬度會對支座負彎矩產(chǎn)生明顯的削峰效果，跨中正彎矩也會因?qū)挾鹊脑黾佣鴾p小［8］，故對寬支座連續(xù)梁的簡化計算采用折減的方法，即先忽略支座寬度的影響，按照彈性點支承進行計算，再對計算結(jié)果進行相應的折減，即

式中：Mmax為寬支座彎矩最大值；Mi為假定彈性點支承情況對應支座或跨的彎矩極大值；M'為折減彎矩。

2.1 負彎矩最大值

圖4為5跨連續(xù)梁（左半部分）在均布荷載q作用下的計算圖示。其中，圖4-a為彈性點支承的情況，圖中Ni為點支承支座反力；圖4-b為考慮支座寬度的情況，由分布彈性支承模型，假設每個支座下的反力由2n+1根彈簧代替，每根彈簧間距d相等，每根彈簧反力 fij相等，fij為第i個支座下的第j根彈簧，最靠近懸臂端的支座i＝1，其他依序編號。

圖4 5跨連續(xù)梁簡化計算圖Fig.4 Simplified calculation diagram for 5-span continuous beam

假定負彎矩最大值出現(xiàn)在支座3的中心，利用疊加法可分別求出其在彈性點支承和寬支座支承下的負彎矩和（順時針為正）

式中：Mq為均布荷載對支座3中心的力矩大?。籰ij為彈簧反力對支座3中心的力臂。若假設相同荷載下，圖4-a和圖4-b的各支座總反力對應相等，即

由式（7）可知，支座寬度增加，與點支承相比，分布彈性支承對支座3中心產(chǎn)生了額外的正彎矩N3βl/8（即折減彎矩M'），且這一正彎矩是由支座3本身貢獻的，與其余支座無關。

若實際負彎矩最大值出現(xiàn)在支座i處，同樣能得到上述結(jié)論。并且考慮到實際上支座范圍內(nèi)的反力并非均勻分布以及其他因素的影響，引入系數(shù)k1對求得的折減彎矩進行調(diào)整，合并常系數(shù)后得到M'=k1Niβl，從而得到支座負彎矩最大值的計算公式，關鍵在于確定系數(shù)k1。

為了得到普遍規(guī)律，式（8）兩邊同時除以ql2，得到準數(shù)

2.2 正彎矩最大值

對于正彎矩最大值，假定其出現(xiàn)在跨i的跨中，文獻［10］是采用寬支座與點支承兩種情況下的彎矩比值來折減的，即，給出了折減系數(shù)K'的計算表。K'綜合考慮了支座相對寬度β和相對剛性系數(shù)α的影響。本文采用類似的方法，令k2=K'-1，則，關鍵在于確定系數(shù)k2。寫成準數(shù)形式

2.3 k1和k2的確定

對于圖2中的5跨分布彈性支承連續(xù)梁，本文利用ANSYS建立其有限元模型（圖5）。梁體采用beam4單元（跨長l取1），分布彈簧采用Combin14彈簧單元，懸臂長ml=m，支座寬度βl=β，彈簧間距取1/120l=0.008 33。梁體的抗彎剛度取EI，相對剛性系數(shù)取α，則支座的豎向剛度系數(shù)K=6EI/αl3，單根彈簧的剛度系數(shù)。豎向荷載q取1。

考慮到豎向荷載下梁不會產(chǎn)生水平位移，故為梁體施加了水平約束；彈簧單元與梁單元連接處采用節(jié)點耦合；彈簧單元下端采用全約束。

考慮到實際工程中縱橫梁的懸臂長度都較小，本文模擬了懸跨比m=0.17、0.25、0.33三種情況。結(jié)合現(xiàn)有碼頭工程實際情況，考慮支座相對寬度β為0、0.17、0.33、0.50、0.60和相對剛性系數(shù)α為0.08、0.10、0.20、0.40、0.60，得到不同組合條件下連續(xù)梁的彎矩圖。

數(shù)值分析發(fā)現(xiàn)，連續(xù)梁跨中正彎矩最大值始終出現(xiàn)在邊跨1的跨中，而支座負彎矩最大值出現(xiàn)的位置與支座相對寬度β和相對剛性系數(shù)α有關，可以分為兩種情況：當α=0.08及α=0.1、β≤0.17，出現(xiàn)在支座2處（圖6-a）；其余情況，則出現(xiàn)在支座3處（圖6-b）。

圖5 5跨分布彈性支承連續(xù)梁有限元模型Fig.5 Finite element model of 5-span continuous beam with uniformly distribution elastic bearings

圖6 負彎矩最大值出現(xiàn)位置的兩種情況Fig.6 Two different positions of the maximum negative moment

由圖6-a還可發(fā)現(xiàn)，當負彎矩最大值出現(xiàn)在支座2處時，中間幾個支座的負彎矩極大值相差不大，故簡單起見，進行寬支座連續(xù)梁負彎矩最大值的計算時，可以假定其始終出現(xiàn)在支座3處。

k1和k2的值可根據(jù)數(shù)值分析的結(jié)果確定，令。對于式（9），；對于式（10），。

（1）負彎矩最大值。

表2和圖7-a～7-c為不同懸跨比m下負彎矩最大值的有限元計算結(jié)果。

表2 負彎矩最大值的計算結(jié)果（m=0.33）Tab.2 Calculation results of the maximum negative moment

由圖7-a，k1與β近似成一次函數(shù)關系，且不同相對剛性系數(shù)下（β，k1）幾乎落在同一直線上，用最小二乘法擬合后得到k1=-0.05β+0.131。

由圖7-b和7-c可知，m不同的取值對k1的影響不大，（β，k1）也近似落在k1=-0.05β+0.131上。從而，負彎矩最大值簡化計算公式為

（2）正彎矩最大值。

表3和圖8-a～8-c為不同m取值正彎矩最大值的有限元計算結(jié)果：

由圖8-a～8-c，對于m不同的取值，（β，k2）均大致落在經(jīng)過原點的拋物線附近，且不同相對剛性系數(shù)下拋物線的“陡峭”程度不同，α越大，越平坦。根據(jù)此特征，假設k2與β的平方成正比，用比例系數(shù)A反映不同相對剛性系數(shù)的影響，即k2=Aβ2。懸跨比m不同，A的取值不同。表4為不同懸跨比下由最小二乘法擬合得到的A的取值，對于其他的懸跨比和相對剛性系數(shù)，可由表中數(shù)據(jù)線性插值求得。從而，正彎矩最大值的簡化公式為

圖7 調(diào)整系數(shù)k1與β關系曲線Fig.7 Curve k1-β

圖8 折減系數(shù)k2與β關系曲線Fig.8 Curve k2-β

表3 正彎矩最大值的計算結(jié)果（m=0.33）Tab.3 Calculation results of the maximum positive moment （m=0.33）

為了進一步研究連續(xù)梁兩端懸臂對正負彎矩最大值的影響，本文補充模擬了懸跨比m=0.42、0.50和0.67三種情況，研究發(fā)現(xiàn)，隨著懸臂的增長，懸臂端支座1的負彎矩迅速增大，靠近懸臂的邊跨1的正彎矩則逐漸減小，負彎矩最大值出現(xiàn)在支座1處，正彎矩最大值出現(xiàn)在邊跨1的跨中。因此，本文的簡化公式僅適用于懸臂較短的情況，即懸跨比m≤0.33，高樁碼頭連續(xù)縱橫梁大多滿足這個條件。

表4 不同懸跨比、相對剛性系數(shù)A的取值Tab.4 Value of A for different ratios of cantilever to span and relative rigidity coefficients

3 算例

計算結(jié)構(gòu)選取某典型寬樁帽高樁碼頭縱梁部分（圖9）?？v梁總長為60 m，排架間距l(xiāng)為7 m，共有8跨，兩端懸臂長為2 m，支座寬度為1.5 m?？v梁為鋼筋混凝土結(jié)構(gòu)EI=3.0×106kN·m2，樁的剛性系數(shù)為473 849 kN/m，相對寬度β=0.21，相對剛性系數(shù)α=6EI/Kl3=0.11，懸跨比m=0.29。考慮豎向均布荷載q=150 kN/m。

圖9 某典型寬支座高樁碼頭縱梁示意圖Fig.9 Longitudinal beam of a certain high-piled wharf engineering

建立8跨連續(xù)梁有限元模型，不考慮支座寬度時，計算結(jié)果如表5；考慮支座寬度時，計算結(jié)果如表6（自左往右分別為支座1、跨1、支座2、跨2…，對稱的支座和跨從略）。

由表5和表6可知，8跨連續(xù)梁支座3、支座4、支座5處的負彎矩極值相差甚小，跨3和跨4的正彎矩極值之差也是如此，說明本文基于5跨連續(xù)梁進行分析是合理的。另外，負彎矩最大值出現(xiàn)在支座3處，正彎矩最大值出現(xiàn)在邊跨1的跨中，也與前面的分析一致。

（1）負彎矩最大值。

與表6中負彎矩最大值-430.18 kN·m相比，兩者誤差0.67%，具有較高精度，滿足工程要求。

（2）正彎矩最大值。

α=0.11，m=0.29，由表4用內(nèi)插法求得A=-1.009。

表5 8跨連續(xù)梁有限元計算結(jié)果（點支承）Tab.5 Calculation results of 8-span continuous beam by FEM（point supporting）

表6 8跨連續(xù)梁有限元計算結(jié)果（寬支承）Tab.6 Calculation results of 8-span continuous beam by FEM（wide supporting）

與表6中正彎矩最大值442.92 kN·m相比，兩者誤差0.56%，具有較高精度，滿足工程要求。

4 結(jié)論

對于高樁碼頭寬支座縱橫梁，一般其懸臂長度較短，本文基于5跨分布彈性支承連續(xù)梁模型，在一定的假設條件下，導出了均布荷載下寬支座連續(xù)梁彎矩最大值的簡化計算公式。經(jīng)過有限元模擬發(fā)現(xiàn)，當支座相對寬度β在0～0.6，相對剛性系數(shù)α在0.08～0.6變化時，負彎矩的最大值總出現(xiàn)中間支座附近，正彎矩的最大值相應地出現(xiàn)在邊跨的跨中。根據(jù)數(shù)值分析結(jié)果，得到適用于懸臂長度較短連續(xù)梁（懸跨比小于0.33）的正負彎矩最大值的無量綱計算公式：。針對某典型寬支座高樁碼頭縱梁，分別用簡化公式和有限元對彎矩最大值進行計算，結(jié)果表明兩者誤差小于5%，計算精度滿足工程要求，具有一定的參考價值。

本文只考慮了豎向均布荷載，對于其他荷載情況有待于進一步討論。

［1］周偉才.大樁帽碼頭排架計算方法及模型試驗研究［D］.南京：河海大學，2005.

［2］任忠.樁帽對橫梁內(nèi)力計算的影響［J］.水運工程，1989（7）：52-53.

［3］JTS 167-1-2010，高樁碼頭設計與施工規(guī)范［S］.

［4］連競，宋向群.結(jié)點剛性域?qū)Ω邩洞a頭排架內(nèi)力的影響［J］.水運工程，1990（9）：12-18.

［5］周偉才，陶桂蘭.考慮節(jié)點域和剪切變形影響的高樁碼頭橫向排架的內(nèi)力計算［J］.水利水運工程學報，2015（1）：56-59. ZHOU W C,TAO G L.Internal force calculation of bents of standing-pile wharfs considering influences of joint-panels and shear deformation［J］.Hydro-Science and Engineering,2005（1）：56-59.

［6］陸東漢.用寬支承方式計算高樁梁板式碼頭上的軌道梁［J］.水運工程，2004（12）：36-39. LU D H.Calculating the Rail Bearer on High-piled Slab&Girder Wharf in the Wide Bearing Form［J］.Port and Waterway Engineering,2004（12）：36-39.

［7］汪基偉.五跨彈性支承連續(xù)深受彎梁的內(nèi)力計算［J］.河海大學學報（自然科學版），1998，26（3）：109-111. WANG J W.Calculation of Internal Forces of Elasticly Supported Continuous Deep Beams［J］.Journal of Hohai University:Natural Sciences,1998,26（3）：109-111.

［8］張華平，羅文華，張金華.支座寬度及相對剛性系數(shù)對連續(xù)梁彎矩的影響［J］.重慶交通大學學報（自然科學版），2010，29 （1）：8-10. ZHANG H P,LUO W H,ZHANG J H.Influence of Abutment Width and Relative Rigidity Coefficients on Moment of Continuous Beam［J］.Journal of Chongqing Jiaotong University:Natural Sciences,2010,29（1）：8-10.

［9］王元戰(zhàn).港口與海岸水工建筑物［M］.北京：人民交通出版社，2013.

［10］吉明.寬支座連續(xù)梁模型試驗及簡化計算方法研究［D］.南京：河海大學，2006.

Research on simplified calculation method of bending moment of continuous beam with wide bearing

GAO Peng，HUANG Shan?shan，OU Yang?feng，CHEN Da
（College of Harbor，Coastal and Offshore Engineering，Hohai University，Nanjing 210098，China）

In the design of high-piled wharfs,the maximum value of bending moment of longitudinal and transverse beams is the leading factor.Based on the model of continuous beam with uniformly distribution elastic bearing,a simplified formula of the maximum value of positive and negative bending moment was put forward in this paper.For continuous beam with different ratios of cantilever to span,simulations of bearings with different relative width and relative rigidity coefficient were carried out through ANSYS finite element software.The coefficient in the simplified formula was determined by the numerical analysis results.Taking the longitudinal beam of a certain highpiled wharf engineering as an example,calculation results of the maximum value of bending moments were compared between the simplified formula and ANSYS.The error between them is less than 5%.

wide bearing;longitudinal and transverse beam;maximum value of bending moment;simplified calculation;ANSYS

TV 698

A

1005-8443（2016）01-0089-07

2015-05-20；

2015-09-15

高鵬(1992-)，男，江蘇南通人，碩士研究生，主要從事港口與近海工程結(jié)構(gòu)計算研究。

Biography：GAO Peng（1992-），male，master student.