余興勝

(中鐵第四勘察設(shè)計院集團有限公司,武漢 430063)

截面特性計算是結(jié)構(gòu)有限元分析和截面檢算的首要工作。以往矩形、圓形、圓端形等規(guī)則截面多采用各自特定的截面特性計算公式直接計算[1-2]，復(fù)雜異形截面多采用網(wǎng)格分塊模擬的辦法間接計算[3-6]。特定截面特性公式計算效率高，但僅適用于特定形式，分塊模擬算法適用面廣，但精度和效率卻低下，且需要借助一些通用的制圖或分析軟件，如AutoCAD、Ansys、Midas等[7-10]。為找到一種適用于任意截面形式的截面特性計算方法，開發(fā)通用高效的計算軟件，亟需對截面特性計算原理和程序開發(fā)架構(gòu)進行深入研究。

鐵四院在開發(fā)軌道交通橋梁協(xié)同設(shè)計系統(tǒng)時，為了實現(xiàn)通用截面檢算功能，通過對截面特性中面積、靜矩、慣性矩和慣性積的數(shù)學(xué)表達式進行深入研究，找到了一種利用數(shù)值積分直接求解截面特性的通用算法[11-12]。即運用格林公式將截面特性的數(shù)學(xué)表達式中對面積的二重積分轉(zhuǎn)換為對截面邊界曲線的積分，通過代入曲線參數(shù)方程將邊界曲線的二元積分轉(zhuǎn)換為一元積分，進而利用牛頓-萊布尼茲公式精確求解截面特性[13]。該方法直接采用數(shù)值積分算法求解截面特性數(shù)學(xué)表達式，不存在任何假定或模擬[14]，與網(wǎng)格分塊模擬的方法相比，計算效率大大提高，并且計算結(jié)果精確，不受網(wǎng)格分塊程度影響。基于該數(shù)值積分算法，采用面向?qū)ο蟮能浖O(shè)計方法開發(fā)通用截面特性計算軟件，提供靈活方便的數(shù)據(jù)輸入組織方式，適用于實體或空心，整體或分離，單一或組合材料的任意組合截面形式[15]。圖1為通用截面構(gòu)建邏輯圖，其中X、Y代表輪廓頂點的X、Y坐標，F(xiàn)代表圓弧凸度或矢高，E代表橢圓離心率。下面將分別就其中關(guān)鍵的整體技術(shù)路線、數(shù)值積分求解和軟件設(shè)計實現(xiàn)三部分進行介紹。

圖1 通用截面構(gòu)建邏輯圖

1 整體技術(shù)路線

截面特性中面積、靜矩、慣性矩和慣性積的數(shù)學(xué)表達式均為在截面平面閉區(qū)域上對面積的二重積分。如面積是截面內(nèi)各微元面積的積分[1]

(1)

靜矩是截面內(nèi)各微元面積與各微元至指定坐標系X軸或Y軸距離乘積的積分[1]

(2)

(3)

慣性矩是截面內(nèi)各微元面積與各微元至指定坐標系X軸或Y軸距離二次方乘積的積分[1]

(4)

(5)

慣性積是截面內(nèi)各微元面積與各微元至指定坐標系X軸距離和Y軸距離乘積的積分[1]

(6)

由格林公式可知，在平面閉區(qū)域D上的二重積分可以通過沿閉區(qū)域D的邊界曲線L上的曲線積分來表達[14]。設(shè)閉區(qū)域D由分段光滑的曲線L圍成，函數(shù)P(x，y)及Q(x，y)在D上具有一階連續(xù)偏導(dǎo)數(shù)，則有

(7)

式中，L是D的取正向邊界曲線，L的正向如下，當觀察者沿L方向行走時，D內(nèi)在他近處的那一部分總在他的左邊。

利用格林公式可以將上述截面特性計算中對面積的二重積分轉(zhuǎn)換為對輪廓邊界曲線的積分。采用輪廓邊界曲線頂點的凸度坐標或矢高坐標的方式組織數(shù)據(jù)，利用常見的直線段、圓弧段或橢圓弧段組合而成的復(fù)合曲線來分段模擬任意輪廓邊界形狀，可以建立輪廓邊界曲線的參數(shù)化方程。但是格林公式右端對輪廓邊界曲線的積分是含有x，y兩個參數(shù)的二元函數(shù)積分，依然不好求解。可以利用直線段、圓(橢圓)弧段的參數(shù)方程將格林公式右端轉(zhuǎn)化為一元函數(shù)積分。

直線段的參數(shù)方程

(8)

(9)

式中，a、c分別為直線段方位角的余弦值、正弦值；b、d分別為起點的X、Y坐標；t[0，len]為直線段上某點距離起點的長度變量；len為直線段長度。

圓弧段的參數(shù)方程

(10)

(11)

式中，r為圓弧半徑；a、b分別為圓心的X、Y坐標；t[sa，ea]為直線段上某點距離起點的弧長變量；sa、ea分別為圓弧的起點、終點方位角。

橢圓弧段的參數(shù)，只需乘以橢圓弧Y軸長與X軸長的比值k即可，圓弧段k=1，是橢圓弧段的特殊情況，為減少變量，本文接下來均以圓弧段代表，公式中將不再出現(xiàn)參數(shù)k，橢圓弧的公式讀者可自行推導(dǎo)。

在一元函數(shù)積分學(xué)中，牛頓-萊布尼茲公式

(12)

表示：F′(x)在區(qū)間[a，b]上的積分可以通過它的原函數(shù)F(x)在這個區(qū)間端點上的值來表達[14]，故格林公式右端代入曲線的參數(shù)方程后的一元積分可以通過尋找其原函數(shù)來求解。在接下來的數(shù)值積分求解部分將分別推導(dǎo)面積、靜矩、慣性矩和慣性積的曲線積分函數(shù)的原函數(shù)。

2 數(shù)值積分求解

2.1 面積

格林公式中函數(shù)P(x，y)和Q(x，y)可以有很多種形式，對應(yīng)的格林公式右端對邊界曲線的積分公式就會不同，雖然單個曲線段的積分結(jié)果不一樣，但是無論函數(shù)P(x，y)和Q(x，y)采用哪一種形式，整個閉合邊界曲線的各曲線段積分的代數(shù)和一定是一樣的。如面積積分時函數(shù)P(x，y)和Q(x，y)分別取(Q=2x，P=y)和(Q=0.5x，P=-0.5y)，格林公式左端偏微分計算的結(jié)果是一樣的，雖然對應(yīng)的單個曲線段的積分結(jié)果不一樣，但整個閉合邊界曲線各曲線段積分的代數(shù)和是一樣的。為降低后續(xù)數(shù)值積分求解的難度，提高計算效率，要力求使右端曲線積分形式最簡。如直接假定函數(shù)P(x，y)和Q(x，y)中的一個為零可以大大降低求解的難度[12]。

令Q=x、P=0，則

因此可將面積的二重積分公式用格林公式轉(zhuǎn)換成對邊界曲線積分公式

(13)

將直線段的參數(shù)方程式(8)、式(9)代入式(13)，可得面積的直線段定積分公式

利用牛頓-萊布尼茲公式求解上式時，推導(dǎo)面積A的直線段定積分函數(shù)的原函數(shù)

FAL=0.5act2+bct

(14)

同理，將圓弧段的參數(shù)方程式(10)、式(11)代入式(13)，可得面積的圓弧段定積分公式

利用牛頓-萊布尼茲公式求解上式時，推導(dǎo)面積A的圓弧段定積分函數(shù)的原函數(shù)

FAA=0.5r2t+0.25r2sin2t+arsint

(15)

2.2 靜矩

與面積的二重積分公式一樣，靜矩的二重積分用格林公式轉(zhuǎn)換成邊界曲線積分時P(x，y)和Q(x，y)也可以有很多種形式，可假定其中一個為零來降低求解的難度。

2.2.1靜矩Sy

令Q=0.5x2、P=0，則

(16)

將直線段的參數(shù)方程式(8)、式(9)代入公式(16)，可得靜矩Sy的直線段定積分公式

利用牛頓-萊布尼茲公式求解上式時，推導(dǎo)靜矩Sy的直線段定積分函數(shù)的原函數(shù)

FSyL=a2ct3/6+abct2/2+0.5b2ct

(17)

同理，將圓弧段的參數(shù)方程式(10)、式(11)代入公式(16)，可得靜矩Sy的圓弧段定積分公式

利用牛頓-萊布尼茲公式求解上式時，推導(dǎo)靜矩Sy的圓弧段定積分函數(shù)的原函數(shù)

(18)

2.2.2靜矩Sx

令Q=0、P=-0.5y2，則

(19)

將直線段的參數(shù)方程式(8)、式(9)代入式(19)，可得靜矩Sx的直線段定積分公式

利用牛頓-萊布尼茲公式求解上式時，推導(dǎo)靜矩Sx的直線段定積分函數(shù)的原函數(shù)

FSxL=-ac2t3/6-0.5acdt2-0.5ad2t

(20)

同理，將圓弧段的參數(shù)方程式(10)、式(11)代入式(19)，可得靜矩Sx的圓弧段定積分公式

利用牛頓-萊布尼茲公式求解上式時，推導(dǎo)靜矩Sx的圓弧段定積分函數(shù)的原函數(shù)

(21)

2.3 慣性矩

慣性矩的二重積分用格林公式轉(zhuǎn)換成邊界曲線積分時P(x，y)和Q(x，y)也可以有很多種形式，可假定其中一個為零來降低求解的難度。

2.3.1慣性矩Iy

令Q=x3/3，P=0，則

(22)

將直線段的參數(shù)方程式(8)、式(9)代入式(22)，可得慣性矩Iy的直線段定積分公式

利用牛頓-萊布尼茲公式求解上式時，推導(dǎo)慣性矩Iy的直線段定積分函數(shù)的原函數(shù)

(23)

同理，將圓弧段的參數(shù)方程式(10)、式(11)代入式(22)，可得慣性矩Iy的圓弧段定積分公式

利用牛頓-萊布尼茲公式求解上式時，推導(dǎo)慣性矩Iy的圓弧段定積分函數(shù)的原函數(shù)

(24)

2.3.2慣性矩Ix

令Q=0、P=-y3/3，則

(25)

將直線段的參數(shù)方程式(8)、式(9)代入式(25)，可得慣性矩Ix的直線段定積分公式

利用牛頓-萊布尼茲公式求解上式時，推導(dǎo)慣性矩Ix的直線段定積分函數(shù)的原函數(shù)

(26)

同理，將圓弧段的參數(shù)方程式(10)、式(11)代入式(25)，可得慣性矩Ix的圓弧段定積分公式

利用牛頓-萊布尼茲公式求解上式時，推導(dǎo)慣性矩Ix的圓弧段定積分函數(shù)的原函數(shù)

(27)

2.4 慣性積

慣性積的二重積分用格林公式轉(zhuǎn)換成邊界曲線積分時P(x，y)和Q(x，y)也可以有很多種形式，可假定其中一個為零來降低求解的難度。

令Q=0.5x2y、P=0，則

(28)

將直線段的參數(shù)方程式(8)、式(9)代入式(28)，可得慣性積Ixy的直線段定積分公式

(0.5b2c2+abcd)t+0.5b2cd)dt

利用牛頓-萊布尼茲公式求解上式時，推導(dǎo)慣性積Ixy的直線段定積分函數(shù)的原函數(shù)

(29)

同理圓弧段的參數(shù)方程式(10)、式(11)代入式(28)，可得慣性積Ixy的圓弧段定積分公式

0.5a2r2sintcost+0.5br3cos3t+

abr2cos2t+0.5a2brcost)dt

利用牛頓-萊布尼茲公式求解上式時，推導(dǎo)慣性積Ixy的圓弧段定積分函數(shù)的原函數(shù)

(30)

3 軟件設(shè)計實現(xiàn)

3.1 程序結(jié)構(gòu)

通用截面特性計算軟件根據(jù)圖1通用截面構(gòu)建邏輯圖設(shè)計，采用面向?qū)ο蟮拈_發(fā)方法實現(xiàn)。程序結(jié)構(gòu)包括任意截面、輪廓形狀(含構(gòu)造器)、邊(直線邊和圓(橢)弧邊)和頂點，其相互之間的包含、繼承關(guān)系如圖2所示。

圖2 程序類結(jié)構(gòu)

其中構(gòu)造器ShapeMaker類的主要功能是基于矩形切割快速自動化生成常規(guī)輪廓形狀頂點集，在數(shù)據(jù)組織部分將對其原理作詳細介紹。頂點Vertex為結(jié)構(gòu)體，有X、Y、F三個字段，X、Y為頂點的x、y坐標，F(xiàn)為該頂點與下一頂點圓弧段的矢高或凸度，若為直線段則F=0。邊Edge類為抽象基類，包含起點、終點兩個頂點Vertex對象，以及Area、Sx、Sy、Ix、Iy、Ixy六個抽象方法。直線邊LineEdge類繼承自邊Edge類，根據(jù)起點、終點坐標計算出直線段長度len，直線段方位角的余弦值a和正弦值c，建立直線段的參數(shù)方程。將參數(shù)t的上下限[0，len]分別代入面積、靜矩、慣性矩、慣性積的直線段定積分函數(shù)的原函數(shù)求解，實現(xiàn)直線邊的Area、Sx、Sy、Ix、Iy、Ixy方法。圓弧邊ArcEdge類繼承自邊Edge類，根據(jù)起點、終點坐標以及矢高或凸度計算出圓弧半徑r、圓心X坐標a、圓心Y坐標b、圓弧段起點方位家sa和終點方位角ea，建立圓弧段的參數(shù)方程。將參數(shù)t的上下限[sa，ea]分別代入面積、靜矩、慣性矩、慣性積的圓弧段定積分函數(shù)的原函數(shù)原求解，實現(xiàn)圓弧邊的Area、Sx、Sy、Ix、Iy、Ixy方法。

平面單連通區(qū)域就是不含有“洞”(包括點“洞”)的區(qū)域，復(fù)連通區(qū)域是含有“洞”(包括點“洞”)的區(qū)域[14]。輪廓形狀Shape類用于描述一個平面單連通區(qū)域或者復(fù)連通區(qū)域的一個“洞”，基本屬性包括彈模系數(shù)，是否外輪廓(復(fù)連通區(qū)域中的“洞”以內(nèi)輪廓表示)，以及封閉邊界的頂點Vertex集。通過封閉邊界的頂點Vertex集創(chuàng)建邊Edge集，當前頂點F為0則創(chuàng)建直線段，若F不為0則創(chuàng)建圓弧段。匯總各個邊Edge的Area、Sx、Sy、Ix、Iy、Ixy，并考慮彈模系數(shù)計算出整個輪廓形狀的換算面積、靜矩、慣性矩和慣性積。任意截面Section類用于描述單個或組合的復(fù)連通區(qū)域，包括各個不同材料的內(nèi)外輪廓形狀Shape集，外輪廓加，內(nèi)輪廓減，逐個匯總各個輪廓形狀Shape的Area、Sx、Sy、Ix、Iy、Ixy，計算出整個截面的換算面積、靜矩、慣性矩和慣性積[15]。

3.2 數(shù)據(jù)組織

任意截面的內(nèi)部數(shù)據(jù)組織采用內(nèi)外輪廓形狀集的形式。不同材質(zhì)的內(nèi)外輪廓形狀均采用統(tǒng)一的封閉邊界坐標點集的方式定義。通過設(shè)定各個輪廓形狀的彈模系數(shù)及是否外輪廓，實現(xiàn)整體或分離、空心或?qū)嶓w的多種材料任意組合截面。

直接輸入輪廓形狀的邊界坐標點集的方式較為繁瑣且容易出錯，基于對矩形輪廓的任意切割，可以實現(xiàn)一種簡潔的參數(shù)化輪廓構(gòu)造器，通過對矩形四角分別設(shè)置倒角(斜角、圓角或內(nèi)直角)及四邊分別設(shè)置開槽可實現(xiàn)絕大多數(shù)常見輪廓形式。調(diào)整圓倒角的位置和半徑實現(xiàn)圓形、圓端形(圓倒角等于順橋向長一半)、弧端形(圓倒角大于順橋向長一半)、圓角矩形等。調(diào)整斜倒角的長度、寬度和位置實現(xiàn)八角形、菱形、梯形、平行四邊形等。調(diào)整內(nèi)直倒角的長度、寬度和位置實現(xiàn)T形、L形、啞鈴形等。如圖3所示，取消規(guī)則截面復(fù)選框，即可自動生成輪廓邊界坐標點集[15]。

圖3 規(guī)則輪廓構(gòu)造器及切割形狀示例

鋼筋混凝土截面中的鋼筋可看作是復(fù)連通區(qū)域中的點“洞”，采用鋼筋集的形式定義其鋼筋信息。鋼筋集定義截面每根縱向受力鋼筋的坐標、直徑、類型、成束根數(shù)等信息，手工輸入鋼筋集繁瑣且容易出錯。通過數(shù)學(xué)算法將內(nèi)外輪廓邊界向截面內(nèi)偏移保護層厚度生成布筋中線輪廓，再根據(jù)設(shè)定的直徑、成束根數(shù)、束數(shù)或間距等參數(shù)自動生成鋼筋集，可快速解決鋼筋集的輸入問題。

圖4 通用截面檢算軟件界面

4 應(yīng)用情況

基于該通用截面特性計算方法編寫的通用截面檢算軟件具備靈活方便的數(shù)據(jù)輸入組織方式，適用于實體或空心，整體或分離，單一或組合材料的任意組合截面形式，計算結(jié)果精確，計算效率高。廣泛應(yīng)用在各種橋梁墩臺、梁部、橋塔等的設(shè)計檢算工作中。特別是在其被整合進軌道交通橋梁協(xié)同設(shè)計系統(tǒng)以后，已先后在武漢、昆明、無錫、尼日利亞拉各斯等國內(nèi)外城市軌道交通，穂莞深、廣肇等城際軌道交通，溫州、南京等市域鐵路，蘇州有軌道電車，長沙中低速磁浮，柳州跨座式單軌等項目的初步設(shè)計、施工圖設(shè)計階段成功應(yīng)用[15]。系統(tǒng)成果顯著提高橋梁設(shè)計效率，目前已被鐵四院所有正在開展的橋梁設(shè)計項目所采用。

[1] 孫訓(xùn)方，方孝淑，關(guān)來泰.材料力學(xué)[M].3版.北京：高等教育出版社，1994.

[2] 童森林.橋梁設(shè)計算法新解[M].北京：中國鐵道出版社，1998.

[3] 余倫明，徐鋒.橋梁結(jié)構(gòu)截面特性計算的梯形分塊法[J].科技信息,2008(2)：154，138.

[4] 項貽強，余建華.橋梁結(jié)構(gòu)分析的數(shù)值方法及其程序——在正交橋、斜彎橋設(shè)計中的應(yīng)用[M].北京：人民交通出版社，1993.

[5] 陳常松，陽斌，顏東煌.任意梁截面特性值的有限元計算方法[J].長沙理工大學(xué)學(xué)報，2013(3)：41-46.

[6] 周凌遠，李喬，張士中.采用有限單元法計算梁任意形狀截面特性[J].計算力學(xué)學(xué)報,2008(5)：634-639.

[7] 胡曉倫，王慧萍，陳立山.兩類橋梁構(gòu)件截面的幾何特性計算[J].交通與計算機，2005(3)：86-89.

[8] 智耕.截面特性參數(shù)的計算方法探討[J].價值工程，2012(3)：35-35.

[9] 劉志紅.平面特殊圖形幾何性質(zhì)計算方法研究[J].安陽工學(xué)院學(xué)報，2011(4)：57-59：

[10] 曲春升.計算軟件求解截面特性的實用方法介紹[J].四川建筑，2007(8)：56-58

[11] 李娜，周小勇.任意截面幾何特性研究[J]．山西建筑，2007(7)：58-60

[12] 羅賢海，湯國興，溫革.二維梁截面特性計算原理及軟件[J]．機械設(shè)計，2003(5)：57-59.

[13] Li Bin， Yu Xingsheng， et al. Stress & Strain Analysis of Specially Shaped RC Column in Elastic-Plastic Stage[C]∥Proceedings of The Twelfth International Symposium on Structural Engineering. Beijing: Science Press，2013.

[14] 同濟大學(xué)數(shù)學(xué)教研室.高等數(shù)學(xué)[M].3版.北京：高等教育出版社，1988.

[15] 余興勝.軌道交通橋梁協(xié)同設(shè)計系統(tǒng)成果報告[D].武漢：中鐵第四勘察設(shè)計院集團有限公司，2015.