黃君寶,鄧成發(fā),2
(1. 浙江省水利河口研究院,杭州 310020;2. 浙江省水利防災(zāi)減災(zāi)重點實驗室,杭州 310018)
渡槽是一種輸送渠道水流而跨越河渠、溪谷、洼地、道路等的架空輸水建筑物。世界上最早的渡槽誕生于中東和西亞地區(qū)。中國修建渡槽也有悠久的歷史,據(jù)記載,西漢時修渠所建的“飛渡”即為渡槽,距今約2000年[1]。
我國20世紀(jì)50-80年代修建了大量渡槽,目前仍有相當(dāng)一部分渡槽仍在承擔(dān)輸水任務(wù)。由于建設(shè)年代久,加之當(dāng)時設(shè)計標(biāo)準(zhǔn)低,渡槽運行中出現(xiàn)了不同程度的開裂、滲漏及變形等現(xiàn)象[2,3],嚴(yán)重影響渡槽的運行安全,因此對這些老舊渡槽進(jìn)行安全評價是目前急需開展的工作。由于目前尚缺乏渡槽安全評價的相關(guān)標(biāo)準(zhǔn)或規(guī)范,同時考慮到渡槽本身的結(jié)構(gòu)特性,較多學(xué)者采用有限元法對渡槽進(jìn)行結(jié)構(gòu)安全分析[4-6]。汪興萌等[7]對運行28年的U型渡槽進(jìn)行了外觀質(zhì)量、混凝土強度及混凝土碳化深度等檢測,并建立了ANSYS有限元計算模型,分析了渡槽各控制截面內(nèi)力值。宋志鵬[8]等以建成57年的拱形渡槽為研究對象,采用現(xiàn)場質(zhì)量檢測、現(xiàn)場載荷試驗和三維有限元方法,對渡槽進(jìn)行質(zhì)量檢測和整體結(jié)構(gòu)的靜力、動力安全性評價。鄭重陽[9]等采用三維有限元軟件對矩形渡槽的受力性能進(jìn)行了模擬,研究了渡槽槽身各控制斷面的應(yīng)力及變形分布規(guī)律。
現(xiàn)階段,渡槽有限元計算時基本采用原設(shè)計資料進(jìn)行建模,由于渡槽在運行多年后槽身及支承結(jié)構(gòu)可能已發(fā)生變形,因此采用原設(shè)計資料建模顯然不能反映渡槽實際現(xiàn)狀,同時還有部分渡槽設(shè)計資料存在遺失的現(xiàn)象。本文以運行40年的雙曲拱渡槽為例,采用三維激光掃描技術(shù)對渡槽進(jìn)行實景復(fù)制,建立渡槽三維點云模型,經(jīng)數(shù)據(jù)處理后形成實體模型再進(jìn)行三維數(shù)值計算,實現(xiàn)了渡槽結(jié)構(gòu)的安全性評估,為類似工程提供參考。
大塘口渡槽是一座修建于1978年的空腹式雙曲拱渡槽。槽身為有拉桿加肋矩形槽,凈寬3.82 m,凈深2.5 m,拉桿間距1.95 m。渡槽拱架采用三肋二波,主拱圈采用等截面懸鏈線,凈跨22.0 m,共3跨;腹拱采用連續(xù)微彎拱,凈跨2.0 m。渡槽全長108.9 m,設(shè)計流量16.5 m3/s,加大流量為18.0 m3/s。渡槽現(xiàn)狀照片見圖1所示。
圖1 渡槽現(xiàn)狀Fig.1 Present situation of aqueduct
通過渡槽的外觀質(zhì)量檢測,目前渡槽存在槽身接縫處滲漏較嚴(yán)重,槽板表面裂縫以及橫隔板與邊肋連接處開裂破損等現(xiàn)象。
根據(jù)文獻(xiàn)[10]對渡槽進(jìn)行混凝土抗壓強度(回彈法)檢測,結(jié)果見表1所示。由于渡槽運行年代較久,回彈值齡期修正系數(shù)取0.9[11],由表可見,各構(gòu)件現(xiàn)齡期混凝土抗壓強度推定值均符合設(shè)計要求。
表1 混凝土抗壓強度推定值檢測結(jié)果Tab.1 Test results of concrete strength
碳化深度檢測表明,各構(gòu)件均存在一定程度碳化現(xiàn)象,其中槽板碳化深度相對較大(36.6~46.1mm),接近鋼筋保護(hù)層厚度(52~68 mm),其他構(gòu)件炭化深度在3.1~9.7 mm,碳化深度均明顯小于鋼筋保護(hù)層厚度。
由于渡槽原設(shè)計資料缺失,無法進(jìn)行結(jié)構(gòu)建模,故采用德國Z+F 5010C激光掃描儀進(jìn)行三維掃描,該儀器特性指標(biāo)主要為:最大測程187.3 m,最小測程0.3 m,分辨率0.1 mm,垂直視野范圍320°,水平視野范圍360°。三維激光掃描儀見圖2所示,第1跨三維掃描結(jié)果見圖3所示。
圖2 三維激光掃描儀Fig.2 Three-dimensional laser scanner
圖3 第1跨三維點云圖Fig.3 3d point Cloud of first span aqueduct
三維激光掃描的數(shù)據(jù)處理流程見圖4所示。首先對不同站點的原始點云數(shù)據(jù)進(jìn)行配準(zhǔn),經(jīng)精簡、去噪、濾波等數(shù)據(jù)處理后,再利用逆向工程軟件,如Geomagic對點云進(jìn)行特征分析和數(shù)據(jù)處理,實現(xiàn)離散點云數(shù)據(jù)的NURBS曲面重構(gòu)和三維重建。
圖4 三維激光掃描數(shù)據(jù)處理流程Fig.4 3-D laser scanning data processing flow
根據(jù)已知的設(shè)計資料以及現(xiàn)場測量結(jié)果,模型精度檢測成果見表2所示。由表可以看出,模型尺寸與設(shè)計或?qū)崪y值基本相當(dāng),模型精度總體較高,滿足計算要求。
4.1.1 有限元模型
將重建后三維模型導(dǎo)入有限元軟件后進(jìn)行網(wǎng)格劃分。由于渡槽在結(jié)構(gòu)上具有對稱性,為縮短計算時間,同時基于偏安全角度考慮,本次僅對第1跨進(jìn)行建模分析,共劃分4節(jié)點四面體單元143 582個,見圖5。各構(gòu)件現(xiàn)齡期混凝土抗壓強度推定值均滿足設(shè)計要求,因此計算模型中彈性模量均采用設(shè)計彈性模量。
表2 模型精度檢測成果Tab.2 Achievements of model accuracy detection
圖5 渡槽有限元分析模型Fig.5 Finite element analysis model of aqueduct
4.1.2 計算結(jié)果分析
為評價渡槽的最大受力性能,選擇滿槽工況進(jìn)行計算分析。各控制截面應(yīng)力計算結(jié)果見表3所示,最大主應(yīng)力、最小主應(yīng)力云圖見圖6、圖7所示。從圖和表中看出:滿槽工況下主拱圈最大主壓應(yīng)力為5.185 MPa,出現(xiàn)在中肋上游拱腳,小于混凝土軸心抗壓強度標(biāo)準(zhǔn)(18.7 MPa);主拱圈最小主拉應(yīng)力為0.895 MPa,出現(xiàn)在中肋拱頂,小于混凝土的軸心抗拉強度標(biāo)準(zhǔn)值(1.9 MPa)。
表3 各控制截面應(yīng)力計算值Tab.3 Stress calculated values of control sections
圖6 最大主應(yīng)力云圖Fig.6 Maximum principal stress nephogram
圖7 最小主應(yīng)力云圖Fig.7 Minimum principal stress nephogram
圖8為滿槽工況下渡槽撓度云圖。從圖8中可以看出,滿槽工況下主拱圈撓度自拱腳至跨中逐漸增大,最大撓度出現(xiàn)在拱頂中肋,為5.71 mm,遠(yuǎn)小于規(guī)范允許最大撓度27.5 mm(L/800)。
圖8 渡槽撓度云圖Fig.8 Deflection nephogram of aqueduct
該渡槽為雙曲拱渡槽,結(jié)構(gòu)的整體穩(wěn)定性主要取決于主拱圈,因此有必要對主拱圈進(jìn)行穩(wěn)定性復(fù)核。主拱圈軸線長度為s,換算為直桿時,長度l0=0.36s,直桿兩端承受的平均軸向壓力Nm計算如下[12]:
(1)
(2)
式中:H為拱腳水平推力,kN;φm為半跨拱的弦與水平線的夾角,°;f為拱圈的計算矢高,m;l為拱圈的計算跨度,m。
計算出的Nm即為換算直桿的軸向力設(shè)計值,根據(jù)不同計算工況,再按軸心受壓構(gòu)件進(jìn)行正截面受壓承載力計算,結(jié)果見表4所示。由表4可見,主拱圈縱向穩(wěn)定性滿足規(guī)范要求。
表4 主拱圈縱向穩(wěn)定性安全系數(shù)Tab.4 Safety factor for longitudinal stability of main arch ring
(1)通過三維激光掃描技術(shù)建立渡槽三維點云模型,經(jīng)數(shù)據(jù)處理后形成實體模型再進(jìn)行三維數(shù)值計算是一種切實可行的方法。
(2)通過三維有限元計算及結(jié)構(gòu)穩(wěn)定性復(fù)核,渡槽主拱圈最大主應(yīng)力、最小主應(yīng)力及最大撓度均小于設(shè)計標(biāo)準(zhǔn),且主拱圈縱向穩(wěn)定性滿足規(guī)范要求,因此認(rèn)為渡槽支承結(jié)構(gòu)整體仍處于安全狀態(tài)。
□