劉勇 高軍

(安徽省交通規(guī)劃設計研究院，安徽合肥 230088)

鋼管混凝土拱橋是一種新型橋型，它充分發(fā)揮鋼與混凝土兩種材料的受力特點，且又給人們視覺上美的感受，使得它在橋梁工程中迅速發(fā)展。拱肋結構的形式多種多樣，有啞鈴式、四肢式、多肋式。拱肋結構是主要承重結構，不同的矢跨比，拱軸系數，結構受力是各不相同的。作為工程人員，為了對比、分析其結構受力特性，需要利用有限元工具對其分析。結構有限元分析，首先要明確結構的節(jié)點坐標，因此，對拱肋的節(jié)點準確確定非常重要。

1 問題的提出及分析

本文以湖南省南縣茅草街大橋為工程背景，對主拱肋的坐標確定進行分析。茅草街大橋為(80+368+80)m的中承式系桿拱橋，拱肋斷面形式為四肢式，如圖1所示。顯然，拱肋結構在空間上是三維坐標體系，在計算跨度恒定的情況下，若變化拱軸系數、矢高等參數，則三維坐標在變動，如何準確的確定節(jié)點坐標。下面，就以茅草街大橋主拱橋拱肋坐標確定進行理論上探討。

茅草街大橋拱軸線方程為:

大橋主拱肋的各參數為:f=71.2m，k=0.999 93，m=1.543，l=368m。

如圖2所示，對于B點，將橫坐標x之值代入公式(1)，能求得縱坐標y之值。從圖1可知，B1，B2兩點關于B點對稱(自身坐標系下)，同時根據李特(Ritter)公式:

其中，n為拱肋變厚系數，本橋取值為0.248 4;β為拱軸水平傾角;Ⅰd為拱頂處慣性矩。

對拱軸線方程一階求導，得:

以上是求解B1B2斷面上B1，B2坐標的過程，但一般工程中，我們關心的是要知道如圖2所示的A，C兩點坐標。很顯然，完全用上述方法求解是無法得出我們想要的結果。但是我們可以從如圖1所示上述方法中得到啟示，如圖3所示。

假設拱軸上一點P，AP垂直于P1P，根據幾何關系，我們可利用上述方法得出A點坐標，這樣問題歸結為確定P點坐標。因為我們已知的只有B點坐標，從而后面的工作是利用B點坐標求出P 點坐標。假設點 B(x，f(x))，點 P(x+Δx，f(x+Δx))，若使得 P點滿足AP垂直于PP1，根據幾何關系及拱軸線性質，必有Δx=P2P=APsinβp。這樣問題就是尋找出一個Δx值，使得此處的P點滿足Δx=P2P=APsinβp。顯然，我們需要經過反復試算才能找出這個Δx，實踐經驗告訴我們，這個工作比較繁瑣，我們可借助編制程序幫我們完成這項工作。

2 程序實現

根據本文第二節(jié)的分析，采用試算法求出Δx之值。問題歸結為選擇一個合適的步長，求出合適的Δx。

我們可采用對分法搜索，求出P點。初試Δx0，我們可利用已知點B求得，設B點處橫坐標為x0。如圖2所示，由幾何關系及拱軸性質可知Δx0≤B1B，因此，不妨可以(x0，x0+B1B)為區(qū)間，采用對分法搜索:1)取值x1=1/2(x0+x0+B1B)，令此處為點P1。2)代入到式(4)，計算出拱軸水平傾角β1，再由式(2)，即李特公式算出Ⅰ1，從而算出AP1長度，判別AP1×sinβ1是否大于B1B值。若小于B1B的話，則以(x1，x1+B1B)為區(qū)間重新搜索。若大于B1B的話，則以(x0，x1)為區(qū)間重新搜索。3)對于第二個步驟，作出判別后，重新選擇區(qū)間，重復第1，2步驟，直至算出適合的Δx，這種解法是數值分析法，無法得到解析解，一般 只要滿足一定精度條件，可認為該點是所求的點。可設定的精度條件取為，對于以上設置的精度條件完全滿足工程精度要求。

其中，Δx(n+1)為第n+1次迭代后的數值;Δx(n)為第n次迭代后的數值。以上的分析步驟能求解出上弦形心處A的坐標，采用同樣的分析方法，能求解出下弦形心處C點的坐標。

表1 湖南省交通勘察設計院與本文程序計算結果對比 m

3 程序實踐

根據以上分析，采用Fortran語言編制計算程序，以茅草街大橋為分析對象，用本文的程序確定拱肋節(jié)點坐標。

取幾個斷面比較如下，A，B，C的點位如圖1，圖2所示。

湖南省交通勘察設計院與本文程序計算結果對比見表1。

4 結語

本文程序的結果與湖南省交通勘察設計院提供的結果相差較小，說明本文的分析方法、程序是合理可行的。對工程技術人員來說，如果需要改變拱軸系數、拱肋變厚系數、矢跨等參數作對比分析的話，可以利用本文的方法和程序較容易確定節(jié)點坐標劃分單元進行有限元分析。

[1]范立礎.橋梁工程(下冊)[M].北京:人民交通出版社，1996.

[2]陳寶春.鋼管混凝土拱橋設計與施工[M].北京:人民交通出版社，2000.

[3]顏東煌.橋梁結構電算程序設計[M].長沙:湖南大學出版社，1999.

[4]譚浩強，田淑清.FORTRAN語言程序設計[M].北京:高等教育出版社，1986.