基于VB語言的鋼管及鑄鐵管壁厚的程序計算

李 煒

(太原供水設計研究院有限公司，山西 太原 030009)

?

·計算機技術及應用·

基于VB語言的鋼管及鑄鐵管壁厚的程序計算

李 煒

(太原供水設計研究院有限公司，山西 太原 030009)

對埋地給排水管道進行了受力與變形分析，結合計算機編程語言VB，得到了關于計算管道壁厚的應用程序，并將該程序的計算結果與手冊計算結果作了對比，驗證了計算程序的適用性。

管道，VB程序語言，管壁厚度，彈性模量

0 引言

城鎮基礎設施中的給水、排水及中水等管道工程是城市建設必不可少的基礎性工程。管道在實際應用中，往往處于復雜的受力狀態。有內壓作用的埋地管道的管材應為管壁材質一致且厚度相等的內外壁均勻平滑的平壁管管材。

埋地給水排水管道結構上作用，有管自重、管內水重、管道外壁上的豎向和側向土壓力和水壓力、地面車輛荷載、堆積荷載和人群荷載等。此外給水壓力管道上尚有作用在管道內壁的運行內水壓力，作用在管道系統縱向由環境溫度變化產生的溫度作用，管道縱向出現不均勻沉降時產生的縱向力和管道在轉彎、三通、堵頭等部位由內水產生的軸向推力等。

對埋地剛性管道，管頂豎向土壓力標準值(每延米)，Fsv,k=CcγsHsBc，其中，Cc為填埋式土壓力系數，在我國，對開槽施工管道，開槽的槽寬都比較大，都屬于填埋式埋管情況。其受力形式如圖1,圖2所示。

1 基于力學理論的管道受力分析

1.1 強度分析

1)荷載組合作用下管壁截面上最大彎矩計算：

式中：φ——彎矩折減系數，可取0.7～1.0；Ep——管材彈性模量；kgm,kvm,kwm——管道自重、豎向土壓力、管內水重作用下管壁截面的最大彎矩系數。

2)設計內水壓力作用下，管壁截面上的拉力設計值計算：

N=ψcγQFwd,kr0b0。

式中：ψc——可變荷載組合系數，可取0.90；γQ——設計內水壓力及地面車輛荷載或堆積荷載的分項系數，可取1.4；

r0,b0——計算半徑、計算寬度。

3)計算鋼管管壁截面的最大環向應力：

其中，t0為計算壁厚，當管材為鋼管時，可取t0=t-2；當管材為球墨鑄鐵管時，可取t0=0.975t-1.5。

1.2 穩定驗算

鋼管管壁截面臨界壓力計算：

式中：vs——回填土的泊松比;n——鋼管管壁失穩時的折縐波數，其取值應使Fcr,k為最小值并為不小于2的正整數。

1.3 變形驗算

在準永久組合作用下管道的豎向最大變形為：

2 回填土的綜合彈性模量

在強度計算中，管側回填土的綜合彈性模量的計算，則是通過計算x=Br/D1和y=Ee/En，由x和y的值共同查表得到z=ξ。但在程序中要實現這個二元函數z=f(x,y)，必須將《給水排水工程結構設計手冊》中表8.4.3-3的數據進行曲線擬合，得到函數z=f(x,y)的表達式，然后便可在程序中進行代碼的編寫。

利用Excel對數據處理的功能，將表8.4.3-3的數據輸入到Excel中，形成光滑的曲線，并選擇添加趨勢線的功能，得到擬合曲線的圖形及函數表達式，如圖3,圖4所示。

從圖3，圖4中可發現，有下列關系表達式存在：

z=A(x)lny+B(x)。

其中，A(x)和B(x)是關于x的函數，關系如圖5所示。

根據上述分析，可以得到z=f(x,y)的函數表達式為：

z=f(x,y)=[0.545 5ln(x)-0.892 5]ln(y)+0.049 9x2-0.369 8x+1.587 2。

3 基于VB的編程實現

通過設計程序的界面，優化各功能按鍵的布置及邏輯順序，采用VB程序語言編寫程序代碼。

應用程序的Visual Basic代碼：

Private Sub Command1_Click()

If Option1.Value = True Then

t0 = Val(Text2.Text) - 2

Text8 = 0.018 * Val(Text3.Text) * Val(Text1.Text)

Text9 = 0.0002466 * D0 * Val(Text2.Text)

Text10 = 7.85 * Val(Text1.Text) * D / 1000000

x = Val(Text4.Text) * 1000 / Val(Text1.Text)

y = 0.8

z = (2.892 * Log(x) - 2.055) * Log(y) + 0.0499 * x ^ 2 - 0.3698 * x + 1.5872

Ed = 4 * z

m1 = 0.2424 * Val(Text9.Text) + 0.19939 * Val(Text8.Text) + 0.1224 * Val(Text10.Text) + 1.9782 * Val(Text1.Text) / 1000

m2 = 1 + 0.732 * Ed * (r0 / t0) ^ 3 / 206000

M = 0.75 * m1 * r0 / m2

N = 1.26 * Val(Text6.Text) * r0

Text11 = N / t0 + 6 * M / t0 ^ 2

……

借助VB的強大的程序開發功能，結合上述關于計算管道管壁厚度的力學理論，基于Visual Basic 6.0編碼平臺，開發出了適用于計算管道管壁的.exe文件。

4 計算結果的比較

結合《給水排水工程結構設計手冊》中給出的計算實例，將手冊中相關參數輸入到應用程序中，得到軟件的計算結果。

4.1 軟件計算結果

通過應用程序的分析計算，得計算結果如下所示：

管壁截面最大環向應力：σθ=231.327 MPa。

降溫時截面最大組合應力為182.531 MPa。

升溫時截面最大組合應力為212.109 MPa。

4.2 手冊計算結果

算例中，計算得到的結果如下：

管壁截面最大環向應力：σθ=233.23 MPa。

降溫時截面最大組合應力為183.96 MPa。

升溫時截面最大組合應力為213.57 MPa。

4.3 計算結果對比

截面的最大環向應力：

降溫時截面最大組合應力：

升溫時截面最大組合應力：

結合上述分析得出，兩組計算結果非常接近，最大誤差僅為0.8%。同時，應用程序在編寫設計時，各參數的選取均以考慮最不利情況為前提，因而其計算結果是偏于安全的。

5 結語

綜合上述各節分析，對埋地管道進行了受力與變形的分析研究，結合計算機編程語言VB，得到的關于計算管道壁厚的應用程序，能夠適應實際工程在不同條件下的管道壁厚的計算，考慮的影響因素比較全面，具有實際的應用意義。

[1] 給水排水工程結構設計手冊[M].第2版.北京：中國建筑工業出版社，2007.

[2] GB 50015—2003,建筑給水排水設計規范[S].

[3] 丁皓江,何福保,謝貽權,等.彈性力學中的有限單元法[M].第2版.北京：機械工業出版社，1989.

[4] 徐芝綸.彈性力學[M].第3版.北京：高等教育出版社，1990.

[5] 閆 安.巧用EXCEL處理統計數據表[J].中國統計，2007(9)：99-100.

[6] 唐世雄．EXCEL統計功能的應用研究[J].成都信息工程學院學報，2002，17(4)：160-161.

[7] 顧建軍.Visual Basic程序設計[M].南寧：廣西科學技術出版社，2009.

[8] 段玉琴.VB程序設計課程教學實踐的幾點體會[J].才智，2009(3)：120-121.

[9] 湯蘭芳.《VB程序設計》學習興趣培養探索[J].電腦知識與技術，2009(19)：69-70.

[10] 薛曉萍.VB程序設計教學方法幾點體會[J].電腦知識與技術(學術交流)，2007(12)：83-84.

On steel pipes based on VB language and program calculation of cast iron pipe’s wall thickness

Li Wei

(TaiyuanWater-SupplyDesignResearchInstitute,Taiyuan030009,China)

The paper undertakes the stress and deformation analysis of the buried water-supply and drainage pipelines, achieves the application program in calculating the pipeline’s wall thickness by combining with the computer program language VB, has the comparison between the calculation result and booklet calculation result of the program, and proves the adaptability of the calculation program.

pipeline, VB program language, pipe’s wall thickness, elastic modulus

1009-6825(2017)09-0255-03

2017-01-16

李 煒(1985- )，男，助理工程師

TP319

A