，

(沈陽(yáng)建筑大學(xué)機(jī)械工程學(xué)院，遼寧 沈陽(yáng) 110168)

0 引言

管道清淤機(jī)器人在極端復(fù)雜的排水管道進(jìn)行清淤作業(yè)，其清淤裝置的剛度和強(qiáng)度對(duì)清淤能力起到至關(guān)重要的作用，為驗(yàn)證該管道機(jī)器人的清淤能力，對(duì)清淤裝置進(jìn)行流固耦合強(qiáng)度分析，并應(yīng)用參數(shù)優(yōu)化的方法對(duì)刮刀、扇形盤(pán)和清淤裝置整體進(jìn)行有限元分析，對(duì)外載荷可行域進(jìn)行優(yōu)化計(jì)算。

機(jī)械結(jié)構(gòu)的有限元強(qiáng)度和剛度分析在機(jī)械設(shè)計(jì)、結(jié)構(gòu)優(yōu)化、研究機(jī)械結(jié)構(gòu)可靠性等方面應(yīng)用廣泛。劉平將有限元分析與經(jīng)典史密斯算法進(jìn)行對(duì)比,驗(yàn)證了有限元分析的可靠性[1]；張庶等建立了裝配狀態(tài)下的有限元模型，得到了裝配體之間剛度與強(qiáng)度的相互影響規(guī)律[2]；丁玉慶等通過(guò)對(duì)轎車(chē)車(chē)身的強(qiáng)度、剛度分析，提出改進(jìn)方案，降低了最大應(yīng)力[3]。

參數(shù)優(yōu)化的方法大多應(yīng)用于機(jī)械結(jié)構(gòu)的設(shè)計(jì)和結(jié)構(gòu)優(yōu)化領(lǐng)域，極少應(yīng)用于剛度和強(qiáng)度的分析中，彭歡[4],胡丹梅[5]等均應(yīng)用參數(shù)優(yōu)化的方法對(duì)相應(yīng)結(jié)構(gòu)進(jìn)行優(yōu)化研究。本文應(yīng)用參數(shù)優(yōu)化的方法，將清淤裝置的外載荷、應(yīng)力和應(yīng)變等參數(shù)化，并以應(yīng)力和應(yīng)變作為約束條件，應(yīng)用ANSYS進(jìn)行了參數(shù)優(yōu)化計(jì)算。

1 管道機(jī)器人及清淤裝置

管道機(jī)器人模型見(jiàn)圖1，由推進(jìn)系統(tǒng)和清淤裝置組成，推進(jìn)系統(tǒng)為整機(jī)提供前進(jìn)動(dòng)力。

圖1 管道機(jī)器人模型

在推進(jìn)系統(tǒng)和主軸旋轉(zhuǎn)時(shí)實(shí)現(xiàn)刮削—攪拌—過(guò)濾—推進(jìn)—自流沖刷的五位一體清淤動(dòng)作。其中，刮刀的刮削，刮刀、扇形盤(pán)和葉片的攪拌對(duì)后續(xù)的清淤效果起到至關(guān)重要的作用，清淤裝置剖視圖見(jiàn)圖2，該清淤裝置具有過(guò)載保護(hù)和越障的能力，通過(guò)刮刀前錐面、扇形盤(pán)(滑塊)、彈簧以及滑道實(shí)現(xiàn)。在高速旋轉(zhuǎn)清淤時(shí)，刮刀、扇形盤(pán)及葉片的強(qiáng)度和剛度是衡量該清淤裝置工作能力的重要參數(shù)，需對(duì)其進(jìn)行研究。

圖2 清淤裝置剖視圖

2 流固耦合強(qiáng)度分析

對(duì)清淤裝置進(jìn)行流固耦合分析，得到流體作用下清淤裝置的應(yīng)力分布狀態(tài)，驗(yàn)證流體對(duì)清淤裝置強(qiáng)度的影響。

通過(guò)對(duì)污泥的流變特性進(jìn)行分析，確定清淤攪拌轉(zhuǎn)速范圍為20～100 r/min時(shí)污泥的黏度趨向于極限黏度[6]，查閱城市排污管道相關(guān)資料[7]，設(shè)定清淤工況，參數(shù)如表1。

表1 清淤裝置的工況參數(shù)

得到流體作用于清淤裝置的最大應(yīng)力分布，需在清淤盤(pán)的極限轉(zhuǎn)速為100 r/min，軸向前進(jìn)速度為2.5 m/s和流體流速為2.5 m/s時(shí)，進(jìn)行流固耦合仿真分析。得到流體作用下清淤裝置的表面壓力分布圖見(jiàn)圖3，清淤裝置基于壓力的流體運(yùn)動(dòng)曲線見(jiàn)圖4。

圖3 清淤裝置的表面壓力分布

圖4 基于壓力的流體運(yùn)動(dòng)曲線

仿真結(jié)果分析：

清淤裝置最大表面壓力和流體最大壓力均出現(xiàn)在清淤裝置貼近管道下方處，有利于刮刀將管道底部的淤泥攪拌至推移狀態(tài)。

清淤裝置最大壓力為68.5 Pa，流體最大壓力為0.1 MPa，清淤裝置的屈服應(yīng)力為280 MPa，故流體對(duì)結(jié)構(gòu)作用的壓力遠(yuǎn)遠(yuǎn)小于屈服極限，故流體對(duì)結(jié)構(gòu)的強(qiáng)度影響很小。

(三)目前我國(guó)經(jīng)濟(jì)發(fā)展速率已經(jīng)出現(xiàn)了一定的下滑，且隨著人民幣貶值、國(guó)內(nèi)經(jīng)濟(jì)轉(zhuǎn)型等客觀壓力。經(jīng)濟(jì)發(fā)展中的不穩(wěn)定因素得到進(jìn)一步凸顯。而形成有效的國(guó)際財(cái)經(jīng)合作體系則是共同抵抗風(fēng)險(xiǎn)的有效路徑。通過(guò)建立穩(wěn)步的國(guó)際財(cái)經(jīng)合作體系擴(kuò)大有效的基本盤(pán)，綜合提升各個(gè)地區(qū)的經(jīng)濟(jì)穩(wěn)定程度是保障我國(guó)經(jīng)濟(jì)持續(xù)發(fā)展的根本。

如圖4，流體在進(jìn)入清淤裝置時(shí)，隨著清淤裝置進(jìn)行高速旋轉(zhuǎn)；通過(guò)清淤裝置后，流體處于高速無(wú)規(guī)則運(yùn)動(dòng)，清淤效果顯著。

3 外載荷參數(shù)優(yōu)化

清淤裝置在清淤過(guò)程中難免會(huì)遇到較大阻力的工況，其彈簧系統(tǒng)具有越障和過(guò)載保護(hù)功能，該功能主要依靠刮刀的軸向和切向力過(guò)大時(shí)實(shí)現(xiàn)， 經(jīng)過(guò)第2節(jié)的流固耦合分析，可忽略流體對(duì)清淤裝置的強(qiáng)度分析。對(duì)清淤裝置的受力分析見(jiàn)圖5，其中，F(xiàn)τ為切向刮削力，F(xiàn)n為軸向淤積物沖擊力，分析可知刮刀、扇形盤(pán)是發(fā)生最大彎矩的零件。

圖5 受力簡(jiǎn)圖

3.1 建立優(yōu)化模型

參數(shù)優(yōu)化的目的是，在清淤裝置滿足預(yù)定的剛度和強(qiáng)度的基礎(chǔ)上，尋找刮刀錐面處可承受的最大外載荷。參數(shù)化優(yōu)化數(shù)學(xué)模型為：

(1)

xi為不同大小外載荷；Fi,Fj分別為切向刮削力和軸向沖擊力(根據(jù)有限元靜力學(xué)分析試探性設(shè)置外載荷參數(shù)變化范圍：100 N≤Fτ≤5 000 N，100 N≤Fn≤5 000 N)；σmax為清淤裝置所受最大應(yīng)力；[σ]基盤(pán)許用應(yīng)力(查閱資料參照水輪機(jī)葉片的安全系數(shù)最大為2.5，結(jié)構(gòu)的屈服應(yīng)力為280 MPa，故設(shè)定許用應(yīng)力為112 MPa);δmax為基盤(pán)最大變形；[δ]為基盤(pán)許用形變(3 mm，按照管徑直徑1%要求)。

3.2 參數(shù)優(yōu)化計(jì)算

將簡(jiǎn)化的清淤裝置三維模型導(dǎo)入到ANSYS中，因六面體網(wǎng)格相比四面體網(wǎng)格單元數(shù)量少，能降低分析計(jì)算誤差，適用于網(wǎng)格質(zhì)量高的結(jié)構(gòu)，所以采用六面體網(wǎng)格劃分方法進(jìn)行網(wǎng)格劃分，得到有限元模型如圖6，劃分網(wǎng)格單元27 325個(gè)，節(jié)點(diǎn)數(shù)17 566個(gè)，網(wǎng)格排列緊湊規(guī)整，質(zhì)量高。

圖6 有限元模型

對(duì)有限元設(shè)定約束如圖7，在清淤裝置旋轉(zhuǎn)前進(jìn)過(guò)程中，刮刀的錐面可實(shí)現(xiàn)越障和過(guò)載保護(hù)功能,是受外載荷的部位。

圖7 有限元約束

進(jìn)行有限元靜力學(xué)計(jì)算，通過(guò)Parameter Set參數(shù)連接Response Surface Optimization模塊建立參數(shù)優(yōu)化模塊，并計(jì)算。

3.3 參數(shù)優(yōu)化后處理

圖8 外載荷與刮刀關(guān)系

圖9 外載荷與扇形盤(pán)關(guān)系

圖10 外載荷與清淤裝置關(guān)系

結(jié)果分析見(jiàn)表2，軸向力參數(shù)區(qū)間100～5 000 N，切向力參數(shù)區(qū)間1 000～5 000 N。

仿真結(jié)果表明：

在軸向力參數(shù)區(qū)間100～5 000 N內(nèi)，刮刀、扇形盤(pán)和清淤裝置整體的應(yīng)力均滿足設(shè)定安全系數(shù)最大為2.5時(shí)，材料的許用應(yīng)力為112 MPa的要求，其中清淤裝置最大應(yīng)力達(dá)到110 MPa接近許用應(yīng)力值；變形均滿足許用形變3 mm的要求。

表2 仿真結(jié)果分析

在切向力參數(shù)區(qū)間100～5 000 N內(nèi)，刮刀、扇形盤(pán)和清淤裝置整體的變形均在許用形變內(nèi)；刮刀和清淤裝置的最大應(yīng)力大于許用應(yīng)力，扇形盤(pán)的最大應(yīng)力小于許用應(yīng)力，所以切向刮削力的選取需要進(jìn)一步優(yōu)化。

3.4 響應(yīng)面優(yōu)化分析

進(jìn)入Response Surface Optimization項(xiàng)目的 Optimization模塊，建立優(yōu)化參數(shù)的優(yōu)化目標(biāo)，即設(shè)定應(yīng)力上限為112 MPa，進(jìn)行優(yōu)化計(jì)算，得到在應(yīng)力約束的條件下軸向沖擊力與切向刮削力的許用值見(jiàn)圖11；變形上限為3 mm，進(jìn)行優(yōu)化計(jì)算，在變形約束條件下的軸向沖擊力與切向刮削力的許用值見(jiàn)圖12。

圖11 應(yīng)力約束下軸、切向力許用關(guān)系

由圖12分析，在變形約束下，軸向沖擊力和切向刮削力在其變化范圍內(nèi)均可行。

方形為可行點(diǎn)，圓形為不可行點(diǎn)，由圖11分析：在應(yīng)力約束下，切向力最大允許3 500 N，軸向力在100～5 000 N內(nèi)均可行。

圖12 變形約束下軸、切向力許用關(guān)系

4 結(jié)束語(yǔ)

本文應(yīng)用有限元方法，從流體沖擊和固體沖擊兩方面對(duì)清淤裝置的外載荷進(jìn)行了計(jì)算，得到結(jié)論如下：

①應(yīng)用流固耦合的方法對(duì)清淤裝置整體進(jìn)行有限元計(jì)算，驗(yàn)證在極限流體沖擊下，流體對(duì)結(jié)構(gòu)的強(qiáng)度影響不大，可忽略不計(jì)。

②應(yīng)用參數(shù)優(yōu)化的方法，以強(qiáng)度和剛度為約束條件，以外載荷為參數(shù)變量，對(duì)清淤裝置的刮刀端部外載荷進(jìn)行有限元優(yōu)化分析，通過(guò)軸向沖擊力和切向刮削力的可行域，確定清淤裝置的作業(yè)能力，為清淤裝置的過(guò)載保護(hù)和越障系統(tǒng)的設(shè)計(jì)計(jì)算提供理論參考。

③切向力對(duì)結(jié)構(gòu)的應(yīng)力影響顯著，對(duì)變形影響呈線性增加；軸向力對(duì)變形影響顯著，對(duì)應(yīng)力呈線性增加。