螺旋刀具切削土壤的三維有限元分析

螺旋刀具切削土壤的三維有限元分析

螺旋式車床在切割和研磨過程中完成了切割、改進和傳輸，具有結構簡單、拆卸緊湊、能耗低等特點。這是開溝機技術的一個重要開發(fā)方向。螺旋刀具是開溝機的關鍵部件,研究螺旋刀具對土壤的破壞過程,對優(yōu)化螺旋式開溝機的運動參數(shù)和結構參數(shù),降低開溝能耗,具有重要的作用。傳統(tǒng)的耕作機具與土壤相互作用的試驗方法,不僅耗費大量人力、財力和時間,且易受客觀條件的限制,使得試驗結果的精度受到影響。通過計算機模擬,可以快捷地虛擬研究螺旋刀具與土壤的切削過程,更好地揭示螺旋刀具-土壤的工作機理和指導螺旋式開溝機的設計。因此,運用ANSYS/LS_DYNA(顯式非線性動力分析有限元軟件)對螺旋刀具切削土壤進行數(shù)值模擬具有重要的現(xiàn)實意義。目前針對螺旋刀具切削仿真的研究文獻很少,主要困難在于螺旋結構為空間曲面三維結構,難以在有限元軟件中建立模型,但根據(jù)切削機理,與之相仿的研究有刀具切削金屬、耕耘機刀具切削土壤的模擬仿真等。國內的研究,有利用數(shù)學模擬方法探討土壤的切削過程,包括傳統(tǒng)分析方法、離散單元法(DEM)和有限單元法(FEM)。在離散單元法方面,吉林大學的張銳及中國農業(yè)大學的徐泳等分別介紹了目前地面力學研究領域中離散單元法在土壤機械特性動態(tài)仿真中的應用。在有限單元法方面,從2001年開始,國內學者使用有限元法研究土壤切削問題迅速增加,主要包括研究進展的綜合論述以及方法的探討等。從1977年開始,國外學者分別采用二維有限元和三維有限元法對寬、窄齒耕作部件進行土壤切削模擬,并采用有限元法研究深耕鏟的土壤切削性能,有的還研究了動力影響時土壤切削的有限元分析方法。近年來用有限元方法研究土壤切削特征取得了較大的進展。有限元法已成為研究耕作部件對土壤高速切削問題的有效工具,但運用有限元法在螺旋刀具切屑土壤的研究還鮮有報道。筆者運用ANSYS/LS_DYNA進行螺旋刀具切削土壤的顯示動力學分析,探討了其切削過程,分析了切削功耗的大小及螺旋刀具力學特性,旨在揭示螺旋刀具-土壤的工作機理,為螺旋式開溝機的動力選型和結構參數(shù)的優(yōu)化設計提供科學依據(jù)。1有限預算建模1.1維圖的繪制螺旋刀具是螺旋式開溝機的關鍵部件,根據(jù)實際的工況要求,設計螺旋刀具由雙頭變螺距螺旋葉片和空心圓柱軸組成,其中變螺距螺旋葉片的表面為變螺距圓柱正螺旋面,變螺距圓柱正螺旋面的形成是以變螺距螺旋線T為導線,使母線與螺旋線軸線成90°沿導線運動形成。螺旋線的參數(shù)方程為?????x=rcos2πky=rsin2πkz=e(2πr)2k2+(2πrtgα0)k{x=rcos2πky=rsin2πkz=e(2πr)2k2+(2πrtgα0)k式中e=h1?2πrk0tgα0(2πrk0)2e=h1-2πrk0tgα0(2πrk0)2tgαk=4eπrk+tgα0其中k為螺旋線上升圈數(shù),α0為初始螺旋角,αk為螺旋線上升圈時的螺旋角,h1為螺旋上升k0圈時沿螺旋軸線上升的高度,r為圓柱半徑。設螺旋葉片的厚度為n,空心圓柱軸內半徑為r0。參考國內學者的研究,設計r=0.029m,r0=0.021m,α0=10°,k0=2圈,h1=0.3m,n=0.003m。根據(jù)實際的開溝要求確定螺旋葉片的外徑為0.16m,螺旋葉片的高度為0.37m,空心圓柱軸的高度為0.45m。利用Pro/E繪制三維圖(圖1)。工作時,動力傳輸給空心圓柱軸,使其帶動螺旋葉片做直線和旋轉運動,切削土壤。1.2彈性材料的選擇由于切削過程的復雜性,伴隨著應力、應變的急劇變化,包括彈、塑性和斷裂的變形,因此,為準確反映切削的過程,假設如下:Ⅰ.螺旋葉片和空心圓柱軸的材料為各向同性的線彈性材料;Ⅱ.如果讓螺旋葉片同時做水平和旋轉運動,則螺旋葉片做復雜的三維運動,在有限元中難以實現(xiàn)在彈性體上同時加載兩種運動,因此將模型簡化為螺旋刀具的旋轉運動和土壤的水平運動;Ⅲ.根據(jù)土壤材料模型發(fā)展的現(xiàn)狀和計算機運算的水平,在數(shù)值模擬中,只考慮對土壤切削的過程,不考慮開溝過程中土壤的升運過程。1.3維模型的建立螺旋刀具由螺旋葉片和空心圓柱軸組成,而螺旋葉片為空間復雜的三維實體,無法在ANSYS中直接建模,因此采用在Pro/E中由方程建立空間三維變螺距圓柱螺旋線,再利用混合掃描的方法建立變螺距螺旋葉片的三維模型,通過Pro/E與ANSYS的無縫連接,將三維實體模型導入ANSYS軟件。將螺旋刀具的單元定義為LS-DYNAExplicit單元SHELL163,為減小沙漏能,采用全積分的殼單元算法(fullyintegratedshellelement)。考慮到應力、應變數(shù)據(jù)的采集,材料采用線性各向同性(linearisotropic),材料的彈性模量為0.21×1012Pa,密度為7.8×103kg/m3,泊松比為0.288。利用掃掠方法來劃分有限元網(wǎng)格,劃分單元數(shù)為3892個。1.4mohpa模型對土壤進行切削模擬,土壤本構關系的選擇對模擬的準確性影響很大。本研究采用的材料為MAT147(MAT_FHWA_SOIL),材料采用修正的Mohr_Cloulomb屈服準則。此土壤材料對屈服面進行了修正,當剪切力很小時,屈服面是一個光滑的平面,并且垂直于壓力軸,修正后的Mohr_Cloulomb屈服表面表達式為F=?psinφ+J2K(θ)2+Ahypsin2φ?????????????????√?ccosφ=0F=-psinφ+J2Κ(θ)2+Ahypsin2φ-ccosφ=0式中p為壓力;φ為內摩擦角;J2為應力偏張量的第2不變量;K(θ)為張量平面角的函數(shù);c為粘聚力;Ahyp為確定修正后的Mohr_Cloulomb屈服面和標準的Mohr_Cloulomb屈服面相似程度。當Ahyp=0時,上式表示的是標準Mohr_Cloulomb屈服面。該材料同時考慮了塑性硬化、塑性軟化、應變速率效應和孔隙水壓力效應,相對LS_DYNA提供的其他土壤模型,其仿真的結果更加準確。結合武漢地區(qū)的土壤特性(黃棕壤),在仿真中采用的土壤模型的參數(shù):土壤中水分密度為1.0×103kg/m3,土壤含水率為3.4%,土壤密度為2.082×103kg/m3,土壤體積模量為0.35×108Pa,土壤剪切模量為0.2×108Pa,土壤內聚力為0.22×105Pa,內摩擦角為0.436rad,土粒相對密度為2.68,Ahyp為0.997×103,粘塑性參數(shù)為1.1,單元破壞的水平為0.99,第3不變量效應偏心參數(shù)Eccen為0.7。考慮到螺旋刀具對土壤的切削方式,建立土壤的模型為空心的半圓弧狀,中心內半圓半徑為0.161m,外半圓半徑為0.6m,厚度為0.4m。土壤單元定義為LS-DYNAExplicit單元SOLID164,為減少沙漏能,采用全積分算法?？紤]仿真的精度以及計算機的運算能力,對參與切削的土壤進行網(wǎng)格細化,其余土壤網(wǎng)格劃分相對較粗。采用掃掠方法劃分有限元網(wǎng)格,劃分單元數(shù)為14282個。螺旋刀具切削土壤的有限元模型如圖2所示。2結果與分析2.1添加邊界條件利用LS_DYNA對模型進行仿真分析。在仿真分析中,螺旋刀具做旋轉運動,土壤向螺旋葉片做直線運動,兩者接觸后,螺旋葉片不斷的銑削土壤。筆者主要研究在一定的工況下螺旋刀體對土壤的切削過程及完成一個切削量的作業(yè)時,整個切削模型最大的功耗以及螺旋葉片的應力等的變化,為螺旋刀具的運動參數(shù)及結構參數(shù)的優(yōu)化打下基礎。在切削過程中,添加如下邊界條件:Ⅰ.螺旋刀具添加繞z軸旋轉的動力載荷,旋轉速度為480r/min;Ⅱ.土壤添加0.0278m/s的直線位移初速度;Ⅲ.空心圓柱軸的上端面和螺旋葉片的內端面添加z向約束;Ⅳ.土壤外圓弧端面添加了除y向位移外的其它約束,考慮切削過程中,距離切削面較遠的土壤不受到任何擾動;Ⅴ.添加螺旋葉片與土壤之間的接觸為面面侵蝕接觸選項(*ERODING_SURFACE_TO_SURFACE)。設置求解時間為0.25s,在ANSYS前處理器中生成K文件,經修改材料、載荷等參數(shù)后,提交DYNA求解器進行求解。2.2模擬分析的結果1間線間的運動圖3為土壤因切削力的作用而破壞的過程。因螺旋葉片半徑為0.160m,土壤內圓弧初始半徑為0.161m,故兩者之間存在0.001m的間距。經過約0.03s的運動,螺旋葉片開始和土壤接觸,土壤受力被破壞。之后隨著螺旋葉片與土壤接觸范圍的增加,土壤被破壞的面積逐漸增加。土壤因受剪切和擠壓而被破壞,被破壞土壤的邊界線形狀和螺旋葉片外端形狀相似,因土壤的直線位移和螺旋葉片的旋轉運動,使銑削不斷向葉片之間尚未被銑削的土壤推進,實現(xiàn)了螺旋葉片對土壤的連續(xù)切削。2不同單元的土壤能量在切削過程中,總的能耗包括維持螺旋葉片的旋轉和土壤模型的運動所需的動能和兩者相互作用的內能,動能基本保持不變,總的能耗隨著內能的不斷增加而增大。由圖4可知,剛開始切削時,總能耗的變化最強烈,隨著切削的進行,能量的變化逐漸趨于平緩,這主要是因為土壤顆粒在初始被破壞后相互之間的結合力減小,而螺旋葉片還未接觸到新的土壤層,總能量的變化放緩。將總的能耗對時間求導,得到切削功率如圖5所示。當土壤單元失效后,螺旋刀具未能接觸到新的土壤單元,因此會出現(xiàn)功率的波動。由圖5可得最大的功率約為3.7kW,這與根據(jù)理論公式計算所得的4kW基本相符。3螺旋葉片與土壤接觸單元s066-ss65由螺旋刀具在切削土壤過程中某一時刻等效應力的分布,可知Mises等效應力在螺旋葉片的內外端面較大,這主要是因為螺旋葉片的外端面與土壤之間的切削作用以及螺旋葉片內端面的曲率較大,成為應力集中區(qū)域。較大的Mises等效應力值集中于約1.5×108Pa,接近材料的屈服應力,因此在實際切削時可以考慮以螺旋葉片為基體,添加物理特性好的刀片在螺旋葉片端部,進行切削。圖6為螺旋葉片端面兩單元S966與S361的等效應力隨時間變化的曲線,其中S966為螺旋與土壤接觸的單元(A線),S361為螺旋葉片內部單元(B線),因S966為螺旋端面與土壤接觸單元,故在同一時刻其等效應力較S361單元大。由曲線A的變化趨勢和切削過程進行對比,可以得出該單元或其附近單元接觸新的土壤單元時刻及應力的大小,曲線尖峰時刻即是該單元或其附近單元與新的土壤單元接觸時刻。分析螺旋葉片單元的受力,為螺旋葉片參數(shù)的優(yōu)化、刀片參數(shù)的選擇打下基礎。3基于ansys/lsdynas的螺旋刀具切削土壤仿真本試驗結果表明:用所構建的螺旋刀具切削土壤的有限元模型應用于螺旋開溝機螺旋刀具與土壤之間相互作用研究是可行的;采用MAT147(MAT_FHWA_SOIL)的土壤模型,該模型采用修正的Mohr_Cloulomb屈服準則,可同時考慮了塑性硬化、塑性軟化、應變速率效應和孔隙水壓力效應,相對LS_DYNA提供的其他土壤模型,其仿真的結果更加準確,通過數(shù)值模擬標明其適合于土體切削的數(shù)