彭 旋,唐華平
(中南大學(xué)機電工程學(xué)院,湖南長沙 410083)
后裝壓縮式垃圾車是城市中最常見的環(huán)衛(wèi)車型,在垃圾車裝卸載垃圾的過程中,垃圾及填料器產(chǎn)生的載荷會加載到車廂上,從而影響垃圾車車廂的使用壽命,但由于垃圾成分復(fù)雜使得垃圾車在進行加卸載時的受力也很復(fù)雜導(dǎo)致對垃圾車進行力學(xué)分析比較困難,筆者以建立壓縮垃圾本構(gòu)模型為基礎(chǔ),使用ANSYS軟件對后裝壓縮式垃圾車滿載擠入工況進行有限元仿真分析,得到相應(yīng)工況下垃圾車車廂的應(yīng)力分布情況并與現(xiàn)場試驗結(jié)果相對比,驗證了有限元模型的正確性,為以后進行壓縮式垃圾車車廂的應(yīng)力分析提出了一種新方法。
后裝壓縮式式圾車主要運載對象為城市生活垃圾。城市生活壓縮垃圾成分復(fù)雜力學(xué)特性難以確定,經(jīng)研究表明[1]壓縮垃圾的力學(xué)特性與土壤的力學(xué)特性非常的相似,垃圾在被壓縮的過程中,其應(yīng)力應(yīng)變曲線表現(xiàn)為雙曲線,且具有加工硬化特性,符合一般土壤的應(yīng)力應(yīng)變曲線,故可以使用土體力學(xué)中的鄧肯-張本構(gòu)模型模擬垃圾的應(yīng)力應(yīng)變關(guān)系。鄧肯-張模型是國內(nèi)外廣泛采用的巖土模型,其應(yīng)力應(yīng)變曲線為雙曲線,主要優(yōu)點是可以利用常規(guī)三軸固結(jié)排水剪試驗確定模型參數(shù),參數(shù)測試簡單,應(yīng)用方便。
Duncan和Chang根據(jù)常規(guī)三軸試驗的結(jié)果提出了在某一圍壓σ3下,在變形過程中t時刻土體體積模量和彈性模量可以表示為:
由摩爾庫倫定律有:
式中:Rf為材料破壞比;Ei為初始彈性模量;Pa為大氣壓;k為初始切線模量數(shù);Kb為體積模量系數(shù);S為當前材料的應(yīng)力水平[2]。
對于三維情況[3],在應(yīng)力水平 S的計算中引入中間主應(yīng)力,用廣義剪應(yīng)力q代替(σ1-σ3),用平均主應(yīng)力 p代替 σ3,用三維 Mohr-Coulomb準則 qf代替(σ1-σ3)f,則有:
將式( 5) ~ ( 9) 代入原鄧肯- 張模型中,則土體的切線模量為:
根據(jù)實驗結(jié)果,垃圾車車廂裝滿垃圾時,車廂前后的上部區(qū)域都未填充垃圾,如圖1所示,但由于后上部未充滿垃圾區(qū)域很小可視為垃圾已充滿該區(qū)域,ANSYS中建立垃圾有限元模型如圖2所示。
圖1 垃圾未充滿區(qū)域
圖2 壓縮垃圾模型
將式(5)~(9)代入原鄧肯-張模型中,則土體的切線模量為:
垃圾在被壓縮過程中,隨著壓縮次數(shù)的增加垃圾的密度也會隨之增加當壓縮次數(shù)達到一定數(shù)值時垃圾密度的變化很?。?],如圖3所示。
圖3 垃圾壓實密度與裝載次數(shù)關(guān)系曲線
垃圾車在每次裝載過程中推鏟會往后移動一段固定距離,假設(shè)每一次裝載壓縮后的垃圾密度視為均勻分布,進行完第三次壓縮后垃圾密度隨裝載次數(shù)變化很小可視為密度不發(fā)生變化,根據(jù)圖3曲線,設(shè)斜線斜率為k。垃圾的總體積為11.5 m3,實驗中垃圾總質(zhì)量為9 t,第一次壓縮的壓縮垃圾位于車廂最前端,該區(qū)域上部未填滿垃圾,測得體積v1=1.68 m3,第2~5次壓縮的垃圾體積分別相等,可以測得體積分別為v2=1.74 m3,第六次壓縮垃圾的體積v3=2.86 m3,根據(jù)體積與密度的計算有:
式中:ρ1、ρ2、ρ3、ρ4、ρ5、ρ6分別為第 1 次到第 6 次加載的垃圾密度,ρ0為垃圾的初始密度,v1為第一次加載的垃圾體積,v2為第2次到第5次分別加載的垃圾體積,v3為第六次加載的垃圾體積,k為圖3曲線中的斜率,m為壓縮垃圾的總質(zhì)量。
計算得到垃圾根據(jù)壓縮次數(shù)得到的垃圾密度分布圖如圖4所示。鄧肯-張E-B本構(gòu)模型在進行計算時需要用到八個參數(shù),所有參數(shù)都可以通過對垃圾進行三軸試驗獲得。取垃圾初始孔隙比為2.0,壓縮垃圾的鄧肯-張E-B參數(shù)如表1[5]。
圖4 垃圾密度分布圖
表1 壓縮垃圾鄧肯-張E-B參數(shù)孔隙比e=2.0
建立模型使用的單元及材料如表2,3所列。
表2 模型中使用的單元
表3 模型中定義的材料屬性
垃圾車車廂廂壁以及加強筋為薄壁殼體類,加載后受平面法向壓力,故使用殼體(Shell)單元,箱體的鋼板和型鋼使用相同的鋼鐵材料。垃圾與廂體的相互作用通過在垃圾與廂體之間施加接觸單元來實現(xiàn)[6]。
對于車廂模型根據(jù)模型的尺寸使用面單元自由網(wǎng)格劃分,對于垃圾模型考慮到材料的特性以及接觸的收斂,采用六面體網(wǎng)格進行SWEEP網(wǎng)格劃分。
垃圾車車廂通過底架與主車體相連接,在車廂與汽車大梁接觸處設(shè)置全約束,推鏟與支撐架接觸處設(shè)置全約束,約束如圖5所示。
在滿載擠入工況下垃圾車車廂受到垃圾的擠入力F1、壓縮力F2、重力G2以及填裝器重力G1的作用力[7],如圖 6 所示。
圖5 有限元分析模型及邊界約束
圖6 車廂滿載擠入工況受力情況
垃圾車在壓縮的過程中通過刮板對垃圾產(chǎn)生擠入力,擠入力可分解為水平分量 Fb1和豎直分量,其中豎直分量Fa1克服區(qū)域C垃圾自重Gc后最終傳遞到頂板,如圖7所示。假定擠入力均勻地施加到頂板區(qū)域A,則頂板垃圾擠入載荷計算如圖7。
圖7 頂板擠入載荷受力區(qū)域
式中:F1為垃圾擠入載荷,Gc為C區(qū)域垃圾重量,Sa為A區(qū)域頂板面積。
鄧肯-張模型是基于軸向加載的三軸實驗提出來的[9],需要對垃圾進行軸向加載,軸向載荷Fb1大小為擠入力F1的水平分量。同時填料器對車廂的載荷通過鉸點、密封條和鎖緊機構(gòu)傳遞給車廂,將填料器對密封條的壓力力視為均勻分布其大小為重力G1在垂直于密封條上的分量,滿載擠入工況受力情況如表4所列。
表4 滿載擠入工況車廂受力情況
根據(jù)表4中的工況對車廂加載計算,得出應(yīng)力與位移云圖,如圖8所示,從計算結(jié)果可以看出來車廂所受最大應(yīng)力在車廂尾部頂板處,車廂尾部受力比較大,前部受力較小,這與車廂最后一次加載尾部受力較大的實際情況是相符的。
圖8 車廂應(yīng)力云圖
為了驗證計算結(jié)果以及有限元模型的準確性,需要對垃圾車在滿載擠入工況下的車廂進行應(yīng)力測試,并將測試結(jié)果與計算結(jié)果進行對比驗證有限元模型以及壓縮垃圾本構(gòu)模型的準確性。
試驗測試時,在垃圾站垃圾車連續(xù)完成裝載工況刮板裝載垃圾至極限位置,在車廂頂板、側(cè)板、一共選取15個主要點進行應(yīng)力測試,測試三次取平均值。測試點位置如圖9,各測點均布粘貼45°應(yīng)變花。
圖9 車廂側(cè)板與車廂頂板應(yīng)力片位置
在滿載擠入工況下車廂各點所測得的工作應(yīng)力如表5、6所列,測點1和測點7處于垃圾未充滿區(qū)域,未測到應(yīng)力。
從表5與表6中可知,仿真與表對比可以看出誤差基本控制在10%左右,證明有限元模型以及壓縮垃圾本構(gòu)模型是有效的。
表5 車廂頂板實驗與仿真應(yīng)力比 /MPa
表6 車廂側(cè)板試驗與仿真應(yīng)力比較 /MPa
通過實驗測試與計算仿真相互驗證的結(jié)果表明:用基于鄧肯-張模型建立的壓縮垃圾本構(gòu)模型對垃圾車車廂和壓縮垃圾進行有限元分析是可行的。通過實驗驗證計算精度滿足工程要求,可代替繁瑣低效的試驗工作。通過數(shù)值分析,可以對后裝壓縮式垃圾車車廂結(jié)構(gòu)進行優(yōu)化設(shè)計。
[1] 張振營,楊云芳,陳云敏.城市生活垃圾的應(yīng)力壓縮曲線及壓縮參數(shù)[J].浙江理工大學(xué)學(xué)報,2008(1):119 -122.
[2] 劉 招,苗隆德.基于APDL的混凝土面板堆石壩三維非線形有限元分析[J].西北水力發(fā)電,2004:17 -20.
[3] 盧廷浩,邵松桂.天生橋一級水電站面板堆石壩三維非線性有限元分析[J].紅水河,1996:20-23.
[4] 黃仁華,趙由才.城市生活垃圾壓實的工程試驗研究[J].同濟大學(xué)學(xué)報(自然科學(xué)版),2000(3):376-378.
[5] 張振營,陳云敏.固體廢棄物壓縮特性及填埋場沉降研究[J].巖土工程學(xué)報,2005(1):116-120.
[6] 徐延海,賈麗萍.基于接觸面幾何的接觸問題解法[J].機械工程學(xué)報,2002,38(10):131 -134.
[7] 丁繼斌.后裝壓縮式垃圾車實驗研究[J].機械工程師,2005(10):65-67.
[8] 陳樹勛,王素暖,白 斌,等.壓縮垃圾車結(jié)構(gòu)的載荷描述與優(yōu)化設(shè)計[J].機械工程學(xué)報,2008(3):213-219.
[9] 孔德志,朱俊高.鄧肯-張模型幾種改進方法的比較[J].巖土力學(xué),2004(6):971-974.