2130龍門加工中心有限元模型的建立
2130龍門加工中心有限元模型的建立2.2.1有限元法簡介有限元法是求解負責(zé)微分方程近似解的一種常用的且非常有效的工具，是現(xiàn)代 數(shù)字化科技的一種重要的基礎(chǔ)性原理，也是現(xiàn)代工程數(shù)值分析中應(yīng)用廣泛的一種方 法。有限元分析主要包含三個方面：（1)有限元分析的基本數(shù)學(xué)力學(xué)原理；（2) 基于原理所形成的實用軟件；（3)使用時的計算機硬件。有限元法常應(yīng)用于電磁力學(xué)計算、結(jié)構(gòu)力學(xué)計算、流體力學(xué)計算，使用有限元 軟件ANSYS、ANSYS Workbench等進行有限元模擬，能求解由桿、梁、板殼、塊 體等各類單元構(gòu)成的彈性、彈塑性或塑性問題[25]。本文主要采用ANSYS Workbench 分析軟件對橫梁系統(tǒng)進行靜動熱態(tài)特性分析。2_2_2有限元軟件ANSYS Workbench簡介自 ANSYS7.0 開始，ANSYS 公司推出了 ANSYS 經(jīng)典版（Mechanical APDL) 和ANSYS Workbench版兩個版本，并且目前均己開發(fā)至15.0版本。ANSYS Workbench是求解實際問題的新一代產(chǎn)品，把ANSYS系列產(chǎn)品融合在仿真平臺，使數(shù)據(jù)無縫實現(xiàn)傳遞以及共享，提高了仿真的效率，并且保證了仿真模 擬的通用性和精確性[26]。2_2_4有限元模型的簡化和建立XH2130龍門加工中心結(jié)構(gòu)十分復(fù)雜，若完全根據(jù)實際情況建立模型，則網(wǎng)格 劃分的難度會增強許多，在有限元軟件中產(chǎn)生過多有限單元，從而導(dǎo)致計算量過于 龐大，影響計算效率，甚至?xí)斐汕蠼庥嬎憬Y(jié)果無法收斂的情況。有許多零件對于 有限元計算和分析是沒有作用的，為減少有限元網(wǎng)格劃分時間，提高網(wǎng)格質(zhì)量，因 此，在利用Pro/E建立機床三維模型時，必須對整機模型進行適當(dāng)?shù)暮喕ＤＰ秃喕恼_性與有限元計算結(jié)果的準確性有著直接緊密的聯(lián)系，當(dāng)然，網(wǎng)格劃分的精細程度也會影響模型計算結(jié)果的精度。在建立機床的整機模型的過程中， 我們需要在保證分析的準確度和精確度的前提下，對結(jié)構(gòu)進行合適、有效的簡化以 獲取可供分析的結(jié)構(gòu)和可應(yīng)用的計算結(jié)果。XH2130龍門結(jié)構(gòu)主要由立柱、動橫梁、 滑座、滑枕四大部分以及連接各個結(jié)構(gòu)件導(dǎo)軌等組成。同時，各個主要結(jié)構(gòu)件上均 有很多微小的結(jié)構(gòu)，這些結(jié)構(gòu)對于計算結(jié)果的影響很小，但會大大降低計算效率， 因此需要將這些小的結(jié)構(gòu)進行簡化，對模型的簡化按以下原則進行：（1)去除結(jié)構(gòu)中 諸如倒角、圓角、凸緣、溝槽等工藝特征；（2)忽略結(jié)構(gòu)中尺寸較小的孔；（3)除去 定位裝配結(jié)構(gòu)，如不直接受力部位的空洞等；（4)平面化和直接化處理模型中的小 錐度、小曲率曲面；（5)若零部件結(jié)構(gòu)對整機的動態(tài)、靜態(tài)特性影響很小，則忽略 不計[29]。簡化后的模型如圖2-3所示。圖2-3整機簡化模型2.2.5接觸關(guān)系有限元軟件中可以設(shè)置5種接觸關(guān)系[3°]:(1)Bonded (綁定）：為接觸的默認設(shè)置。接觸面或邊之間無滑移、無分離， 忽略了其中的間隙和穿透，屬于線性求解。(2)No Separation (不分離接觸）：類似于綁定接觸。它只適用于面。不允許 接觸區(qū)域的面分離，但是沿著接觸面可以有小的無摩擦滑動。綁定的和不分離的接 觸式最基礎(chǔ)的線性行為，僅僅需要一次迭代。(3)Frictionless (無摩擦接觸）：只適用于面接觸。為標(biāo)準的單邊接觸，也就 是如果出現(xiàn)分離則法向壓力為零，摩擦系數(shù)為零，允許自由滑動，裝配體中會施加 弱彈簧以幫助穩(wěn)定求解，因此，根據(jù)不同的載荷，模型間可以出現(xiàn)間隙。它是非線性求解，因為在載荷施加過程中接觸面積可能會發(fā)生改變。(4)Rough (粗糙接觸）：只適用于面接觸。這種接觸方式與無摩擦接觸類型 相似，但其表現(xiàn)為完全的摩擦接觸即沒有相對滑動。在預(yù)設(shè)情況下，不會自動消除 間隙。(5)Frictional (有摩擦接觸）：只適用于面接觸。摩擦系數(shù)可以是任意非負值。 在發(fā)生相對滑動前，兩接觸面之間可以通過接觸區(qū)域傳遞一定數(shù)量的剪應(yīng)力。這種 情況有點像膠水。模型在滑動發(fā)生前定義一個等效的剪應(yīng)力，作為接觸壓力的一部 分。一旦剪應(yīng)力超過此值，兩面將發(fā)生相對滑動。對機床的各個結(jié)構(gòu)件簡化之后，通過有限元分析，可以建立各個結(jié)構(gòu)件之間的 接觸關(guān)系，為了簡化計算，將動橫梁與立柱之間的導(dǎo)軌接觸以及滑座與動橫梁之間 的導(dǎo)軌接觸關(guān)系設(shè)置成bonded關(guān)系，而為了保證滑枕的計算精度，滑枕與滑座之間 采用的是貼塑導(dǎo)軌，且滑枕處的變形量也較大，將滑枕與滑座之間的貼塑導(dǎo)軌設(shè)置 為 Frictionless。2_2_6定義材料屬性和網(wǎng)格劃分為減少有限元網(wǎng)格劃分時間，提高網(wǎng)格質(zhì)量，在建模時對整機模型進行適當(dāng)?shù)?簡化，然后導(dǎo)入ANSYS有限元析軟件中，材料屬性如表2-2所示。表2-2主要材料參數(shù) 材料 彈性模量（GPa) 密度（g/cm3) 泊松比 HT300(橫梁、立柱） 130 7.35 0.27 QT600-3 (滑枕） 169 7.12 0.28 Bearing (導(dǎo)軌、滑塊） 280 7.82 0.30 運用自由網(wǎng)格劃分類型對整機進行劃分網(wǎng)格，設(shè)置相對大小為Coarse，劃分網(wǎng) 格結(jié)果如圖2-4所示，共有212219個節(jié)點，110945個網(wǎng)格單元[31]。2.2.7約束及載入條件(1)約束條件①機床立柱固定在地面上，故將兩個立柱的地面設(shè)置成固定約束；②動橫梁通過立柱上的滾珠絲杠實現(xiàn)上下移動，設(shè)置絲母座沿導(dǎo)軌Y方向位移 約束；滑座絲母座沿X方向位移約束；滑枕絲母座沿Y方向位移約束。(2)受力情況①自身重力：重量加速度設(shè)置為Y方向，重力加速度大小為9806.6mm/s2。②切削力：設(shè)計切削******扭矩為1200Nm，切削刀盤的直徑為300mm。因此，切削力為：1200/0.15=8000N。切削扭矩作用在滑枕頂部電主軸與滑枕相連接的表面 上，切削載荷作用在滑枕內(nèi)部用于固定銑頭的圓形固定座上，采用軸承載荷的方式 進行加載 Bearing Load=8000N。③在橫梁和滑枕上增加平衡油缸時受力情況由于重型機床的結(jié)構(gòu)件的自身重量的影響，為了降低驅(qū)動電機的功率，減小絲 杠的負擔(dān)，在垂直移動部件上采用了平衡油缸來平衡結(jié)構(gòu)的自身重力。在XH2130 型龍門加工中心中，上下運動（沿Y軸方向）的活動部件包括：動橫梁和滑枕，因 此，在動橫梁和滑枕上各設(shè)計了兩個平衡油缸，用于平衡結(jié)構(gòu)的自身重力。動橫梁平衡油缸固定在兩個立柱上，用拉桿拉住動橫梁的兩端的吊塊，因此， 作用在立柱上的反作用力是未簡化情況下，動橫梁上所有零件重量的總和：即動橫 梁質(zhì)量（9300kg)、滑座質(zhì)量（2220kg)、滑枕質(zhì)量（2108kg)的總和為13628kg， 平均每一個立柱上的反作用力為66822N。而實際載入在動橫梁吊塊上的平衡拉力為 簡化后的模型質(zhì)量總和113570N，平均在每個吊塊上的拉力即為56785N。滑枕平衡油缸是為克服滑枕上所有零件的自身重力的，兩個滑枕平衡油缸分別 位于滑枕的兩側(cè)，并固定在滑座上，拉桿末端固定在滑枕的頂端。作用在滑枕上的 拉力為簡化后的模型重量：94849.8N，而作用在滑塊上的反作用力為滑枕結(jié)構(gòu)的實 際重量：20672N[32]。具體受力和約束情況如圖2-5所示。