大型龍門銑床主軸滑枕有限元分析說明書,大型,龍門,銑床,主軸,有限元分析,說明書,仿單 大型龍門銑床主軸滑枕有限元分析 摘 要 隨著計算機的普及和計算速度的提高，有限元分析已越來越廣泛地應用到工程分析計算問題，其分析手段已成為現代機械設備CAE分析的主流方法，也是解決我國機械產品設計分析水平落后的關鍵手段之一。主軸滑枕是機床中一種必不可少的通用零部件，亦是影響機床工作性能的薄弱環(huán)節(jié)。 本論文的主要目的就是針對機床主軸關鍵部件——滑枕進行相關分析。采用Ansys命令流方式建立高速數控機床主軸滑枕實體分析模型（其中加強筋厚度給出了三種，分別為15mm，20mm，25mm），并對他們在滑動軸承支撐、高速旋轉條件下動力學行為進行有限元分析，求解出系統(tǒng)模態(tài)頻率、模態(tài)振型、瞬態(tài)動力響應和周期激勵下諧響應譜，在分析的基礎上進行比較，確定系統(tǒng)的薄弱環(huán)節(jié)，對大型龍門銑床主軸滑枕的結構變形規(guī)律，給出了定性分析結論，選擇出了三種加強筋厚度中最優(yōu)的一種，同時，在分析的基礎上提出了一些改進措施，為工作中進行機械系統(tǒng)結構CAE分析打下良好的基礎，為滑枕結構優(yōu)化設計提供理論依據和指導。 關鍵詞：CAE；主軸滑枕；有限元分析；變形分析；改進措施 Abstract With the development of computers and the speed of calculation, finite element analysis has become increasingly widely in the calculation of engineering analysis, of which analysis methods has become the mainstream of CAE analysis methods, and the analytical tools has become one of the key method to solve the problem of the level of product design behind about our country mechanical analysis. Spindle ram is the necessary and general components of machine tools, which is also the weak link with frequent failures. The main purpose of this paper is to analysis the spindle - ram which is the key components for machine tool. Ansys command is used to create high-speed flow axis CNC machine tool ram model entities, and do the finite element analysis under the conditions of the sliding bearing support and the high-speed rotary dynamic behavior.Their kinetic analysis, a system for solving modal frequencies, modal vibration mode, the transient dynamic response and the cycle of harmonic excitation response spectrum. In based of the analysis to make sure the weak links in the system,and find out the large gantry milling machine spindle ram deformation of the structure for giving the conclusions of the qualitative analysis, which carry out the work of CAE analysis of the structure of mechanical systems and lay a good foundation for the structural optimization ram designed to provide the theoretical basis and guidance. Keywords: CAE；Spindle ram；Finite element analysis；Deformation analysis；Improvement measures 目 錄 1 緒論................ 1 1.1 有限元分析的發(fā)展現狀 1 1.2 有限元分析的發(fā)展趨勢 2 2 有限元分析理論及ANSYS軟件的介紹………………………………………....5 2.1 有限元分析理論 5 2.2 ANSYS軟件的介紹 7 3 模型及系統(tǒng)分析 10 3.1 模型描述 11 3.2 主軸滑枕靜力分析 12 3.2.1加強筋厚度為15mm時主軸滑枕靜力分析 13 3.2.2加強筋厚度為20mm時主軸滑枕靜力分析 14 3.2.3加強筋厚度為25mm時主軸滑枕靜力分析 15 3.2.4結論 ……………………….. 15 3.3 主軸滑枕動力學分析 16 3.3.1主軸滑枕模態(tài)頻率 16 3.3.2主軸滑枕模態(tài)振型 18 3.3.3主軸滑枕瞬態(tài)動力學分析 21 3.3.4主軸滑枕周期激勵下諧響應 24 4 總結 29 謝 辭……………………………………………………………………………….30 參考文獻…………………………………………………………………………...31 1 緒論 有限元方法的計算思想早在 20 世紀 40 年代就提出了。1943 年 Courant 在一篇 論文中就提出用一組三角形單元和最小位能原理研究了 St.Venant 的扭轉問題。但真 正用以解決工程中的數值計算問題是在上世紀 50 年代電子計算機出現以后，1956 年 Truner、 Clough、 Martin、 和 Topp 發(fā)表了在結構力學采用有限元方法的第一篇論文，但當時還沒有“有限元方法”這個名稱。“有限元方法” （Finite Element Method）這個名稱，是 1960 年 Clough在一篇結構分析的計算論文中首次明確提出的。在我國，上世紀 60 年代初期馮康等人獨立地建立了有限元方法的數學理論并應用于工程實際問題。自從那時候起，不論國際還是國內，隨著高速電子電子計算機的日益廣泛的應用，有限元方法得到廣泛的發(fā)展。近年來隨著計算機技術的普及和計算速度的不斷提高，有限元分析在工程設計和分析中得到了越來越廣泛的重視， 已經成為解決復雜的工程分析計算問題的有效途徑。從“有限元”這個名詞第一次出現，到今天有限元在工程上得到廣泛應用，經歷了三十多年的發(fā)展歷史，理論和算法都已經日趨完善。有限元的核心思想是結構的離散化，就是將實際結構假想地離散為有限數目的規(guī)則或不規(guī)則的單元組合體，實際結構的物理性能可以通過對離散體進行分析，得出滿足工程精度的近似結果來替代對實際結構的分析，這樣可以解決很多實際工程需要解決而理論分析又無法解決的復雜問題。 1.1 有限元分析的發(fā)展現狀 隨著有限元方法研究的深入，過去不能解決或能解決但求解精度不高的問題，都 得到了新的解決方案。 傳統(tǒng)的有限元方法假設：分析域是無限的；材料是同質的，甚至在大部分的分析中認為材料是各向同性的；對邊界條件簡化處理。但實際問題往往是分析域有限、材料各向異性或邊界條件難以確定等。 為解決這類問題， 美國的 Heofanis St Rouboulis、Lin Zhang 等人提出用 GFEM(Generalized Finite Element Method) 解決分析域內含有大量孔洞特征的問題；比利時的 Nguyen Dang Hung 和越南的 Tran Thanh Ngoc提出用 HSM(the Hybrid metis Singular element of Membrane plate) 解決實際開裂問題(結構尺寸有限，形狀任意，邊界條件復雜，材料特性任意)。 傳統(tǒng)的有限元斷裂力學技術(the ?nite element fracture mechanics techniques)在解決零件中出現裂縫這類問題時，需要在曲線型裂紋前緣附近的區(qū)域細分網格。這樣無論是從網格生成的角度看還是從求解的角度看，都需要花費大量的時間。而且在循環(huán)加載的情況下產生的次裂紋將會使分析變得更加復雜。為此，美國的 Daniel SPipkinsay、Satya N Atlurib 提出了 FEAM(Finite Element Alternating Method) 。該方法在求解應力集中因子時，可在不犧牲精度的情況下節(jié)省時間，用它分析具有橢圓裂紋或部分橢圓裂紋的結構是很有用的。此外，西班牙的 Onate E 和波蘭的 Rojek J 將 DEM(Discrete Element Method) 和有限元方法結合解決地質力學中的動態(tài)分析問題；瑞典的 Birgersson F 和英國的 Finnveden S 針對有限元方法在頻域中的應用提出了 SFEM(Spectral Finite Element Method)。 在有限元方法應用領域不斷擴展、求解精度不斷提高的同時，有限元方法也從分析比較向優(yōu)化設計方向發(fā)展。印度 Mahanty 博士用 ANSYS 對拖拉機前橋進行優(yōu)化設計，結果不但降低了約 40% 的前橋自重，還避免了在制造過程中的大量焊接工藝，降低了生產成本。 有限元方法在國內的應用也十分廣泛。20 世紀 80 年代我國大連理工大學工程力學研究所開發(fā)成功了國內第一個通用有限元程序系統(tǒng) JIGFEX，并在 1983 年開發(fā)出了它的微機版 JIG2FEX-W。目前有限元方法已滲透到工程分析的各個領域，從大型的三峽工程到微米級器件都采用有限元進行分析。它在我國經濟發(fā)展中擁有廣闊的發(fā)展前景。 1.2 有限元分析的發(fā)展趨勢 近年來隨著計算機技術的普及和計算速度的不斷提高，有限元分析在工程設計和分析中得到了越來越廣泛的重視， 已經成為解決復雜的工程分析計算問題的有效途徑，現在從汽車到航天飛機幾乎所有的設計制造都已離不開有限元分析計算，其在機械制 造、材料加工、航空航天、汽車、土木建筑、電子電器、國防軍工、船舶、鐵道、石化、 能源、科學研究等各個領域的廣泛使用已使設計水平發(fā)生了質的飛躍，主要表現在以下幾個方面： ? 增加產品和工程的可靠性 ? 在產品的設計階段發(fā)現潛在的問題 ? 經過分析計算，采用優(yōu)化設計方案，降低原材料成本 ? 縮短產品投向市場的時間 ? 模擬試驗方案，減少試驗次數，從而減少試驗經費 國際上早在上世紀 60 年代初就開始投入大量的人力和物力開發(fā)有限元分析程序，但真正的 CAE 軟件是誕生于上世紀 70 年代初期，而近 15 年則是 CAE 軟件商品化的發(fā)展階段，CAE 開發(fā)商為滿足市場需求和適應計算機硬、軟件技術的迅速發(fā)展，在大力推銷其軟件產品的同時，對軟件的功能、性能，用戶界面和前、后處理能力，都進行了大幅度的改進與擴充。這就使得目前市場上知名的 CAE 軟件，在功能、性能、易用性、可靠性以及對運行環(huán)境的適應性方面，基本上滿足了用戶的當前需求，從而幫助用戶解決了成千上萬個工程實際問題，同時也為科學技術的發(fā)展和工程應用做出了不可磨滅的貢獻。目前流行的 CAE 分析軟件主要有 NASTRAN、 ADINA 、 ANSYS、ABAQUS、MARC、MAGSOFT、COSMOS 等。MSC-NASTRAN 軟件因為和 NASA的特殊關系，在航空航天領域有著很高的地位，它以最早期的主要用于航空航天方面的線性有限元分析系統(tǒng)為基礎，兼并了 PDA 公司的 PATRAN，又在以沖擊、接觸為特長的 DYNA3D 的基礎上組織開發(fā)了 DYTRAN。近來又兼并了非線性分析軟件MARC， 成為目前世界上規(guī)模最大的有限元分析系統(tǒng)。 ANSYS 軟件致力于耦合場的分析計算， 能夠進行結構、 流體、 熱、 電磁四種場的計算， 已博得了世界上數千家用戶的鐘愛。ADINA 非線性有限元分析軟件由著名的有限元專家、麻省理工學院的 K.J.Bathe教授領導開發(fā)，其單一系統(tǒng)即可進行結構、流體、熱的耦合計算。并同時具有隱式和顯式兩種時間積分算法。由于其在非線性求解、流固耦合分析等方面的強大功能，迅速成為有限元分析軟件的后起之秀，現已成為非線性分析計算的首選軟件。 縱觀當今國際上 CAE 軟件的發(fā)展情況，可以看出有限元分析方法的一些發(fā)展趨 勢： 1、與 CAD 軟件的無縫集成 當今有限元分析軟件的一個發(fā)展趨勢是與通用 CAD 軟件的集成使用，即在用CAD 軟件完成部件和零件的造型設計后，能直接將模型傳送到 CAE 軟件中進行有限元網格劃分并進行分析計算，如果分析的結果不滿足設計要求則重新進行設計和分析，直到滿意為止，從而極大地提高了設計水平和效率。為了滿足工程師快捷地解決復雜工程問題的要求，許多商業(yè)化有限元分析軟件都開發(fā)了和著名的 CAD 軟件（例如 Pro/ENGINEER、Unigraphics、SolidEdge、SolidWorks、IDEAS、Bentley 和AutoCAD 等）的接口。有些 CAE 軟件為了實現和 CAD 軟件的無縫集成而采用了CAD 的建模技術，如 ADINA 軟件由于采用了基于 Parasolid 內核的實體建模技術，能和以 Parasolid 為核心的 CAD 軟件（如 Unigraphics、 SolidEdge、 SolidWorks）實現真正無縫的雙向數據交換。 2、更為強大的網格處理能力 有限元法求解問題的基本過程主要包括：分析對象的離散化、有限元求解、計算結果的后處理三部分。由于結構離散后的網格質量直接影響到求解時間及求解結果的正確性與否，近年來各軟件開發(fā)商都加大了其在網格處理方面的投入，使網格生成的質量和效率都有了很大的提高，但在有些方面卻一直沒有得到改進，如對三維實體模型進行自動六面體網格劃分和根據求解結果對模型進行自適應網格劃分，除了個別商業(yè)軟件做得較好外，大多數分析軟件仍然沒有此功能。自動六面體網格劃分是指對三維實體模型程序能自動的劃分出六面體網格單元，現在大多數軟件都能采用映射、拖拉、掃略等功能生成六面體單元，但這些功能都只能對簡單規(guī)則模型適用，對于復雜的三維模型則只能采用自動四面體網格劃分技術生成四面體單元。對于四面體單元，如果不使用中間節(jié)點，在很多問題中將會產生不正確的結果，如果使用中間節(jié)點將會引起求解時間、收斂速度等方面的一系列問題，因此人們迫切的希望自動六面體網格功能的出現。自適應性網格劃分是指在現有網格基礎上，根據有限元計算結果估計計算誤差、重新劃分網格和再計算的一個循環(huán)過程。對于許多工程實際問題，在整個求解過程中，模型的某些區(qū)域將會產生很大的應變，引起單元畸變，從而導致求解不能進行下去或求解結果不正確，因此必須進行網格自動重劃分。自適應網格往往是許多工程問題如裂紋擴展、薄板成形等大應變分析的必要條件。 3、由求解線性問題發(fā)展到求解非線性問題 隨著科學技術的發(fā)展，線性理論已經遠遠不能滿足設計的要求，許多工程問題如材料的破壞與失效、裂紋擴展等僅靠線性理論根本不能解決，必須進行非線性分析求解，例如薄板成形就要求同時考慮結構的大位移、大應變（幾何非線性）和塑性（材料非線性）； 而對塑料、 橡膠、 陶瓷、 混凝土及巖土等材料進行分析或需考慮材料的塑性、蠕變效應時則必須考慮材料非線性。眾所周知，非線性問題的求解是很復雜的，它不僅涉及到很多專門的數學問題，還必須掌握一定的理論知識和求解技巧，學習起來也較為困難。為此國外一些公司花費了大量的人力和物力開發(fā)非線性求解分析軟件，如ADINA、ABAQUS 等。它們的共同特點是具有高效的非線性求解器、豐富而實用的非線性材料庫，ADINA 還同時具有隱式和顯式兩種時間積分方法。 4、由單一結構場求解發(fā)展到耦合場問題的求解 有限元分析方法最早應用于航空航天領域，主要用來求解線性結構問題，實踐證明這是一種非常有效的數值分析方法。而且從理論上也已經證明，只要用于離散求解對象的單元足夠小，所得的解就可足夠逼近于精確值?，F在用于求解結構線性問題的有限元方法和軟件已經比較成熟，發(fā)展方向是結構非線性、流體動力學和耦合場問題的求解。例如由于摩擦接觸而產生的熱問題，金屬成形時由于塑性功而產生的熱問題,需要結構場和溫度場的有限元分析結果交叉迭代求解，即“熱力耦合”的問題。當流體在彎管中流動時，流體壓力會使彎管產生變形，而管的變形又反過來影響到流體的流動…… 這就需要對結構場和流場的有限元分析結果交叉迭代求解，即所謂“流固耦合”的問題。由于有限元的應用越來越深入，人們關注的問題越來越復雜，耦合場的求解必定成為 CAE 軟件的發(fā)展方向。 5、程序面向用戶的開放性 隨著商業(yè)化的提高，各軟件開發(fā)商為了擴大自己的市場份額， 滿足用戶的需求， 在軟件的功能、易用性等方面花費了大量的投資，但由于用戶的要求千差萬別，不管他們怎樣努力也不可能滿足所有用戶的要求，因此必須給用戶一個開放的環(huán)境，允許用戶根據自己的實際情況對軟件進行擴充，包括用戶自定義單元特性、用戶自定義材料本構（結構本構、熱本構、流體本構） 、用戶自定義流場邊界條件、用戶自定義結構斷裂判據和裂紋擴展規(guī)律等等。 2 有限元分析理論及ANSYS軟件的介紹 有限元分析，即使用有限元方法來分析靜態(tài)或動態(tài)的物體或系統(tǒng)。在這種方法中一個物體或系統(tǒng)被分解為由多個相互聯結的、簡單、獨立的點組成的幾何模型。在這種方法中這些獨立的點的數量是有限的，因此被稱為有限元。ANSYS軟件是融結構、流體、電場、磁場、聲場分析于一體的大型通用有限元分析軟件。它能與多數CAD軟件接口，實現數據的共享和交換，如Pro/Engineer, NASTRAN, Alogor, I－DEAS, AutoCAD等， 是現代產品設計中的高級CAD工具之一。 2.1 有限元分析理論 有限元分析(FEA，Finite Element Analysis)的基本概念是用較簡單的問題代替復雜問題后再求解，用有限元方法來分析靜態(tài)或動態(tài)的物體或系統(tǒng)。在這種方法中一個物體或系統(tǒng)被分解為由多個相互聯結的、簡單、獨立的點組成的幾何模型。在這種方法中這些獨立的點的數量是有限的，因此被稱為有限元。它將求解域看成是由許多稱為有限元的小的互連子域組成，對每一單元假定一個合適的(較簡單的)近似解，然后推導求解這個域總的滿足條件(如結構的平衡條件)，從而得到問題的解。這個解不是準確解，而是近似解，因為實際問題被較簡單的問題所代替。由于大多數實際問題難以得到準確解，而有限元不僅計算精度高，而且能適應各種復雜形狀，因而成為行之有效的工程分析手段。 由實際的物理模型中推導出來的平衡方程式被使用到每個點上，由此產生了一個方程組。這個方程組可以用線性代數的方法來求解。有限元分析的精確度無法無限提高。元的數目到達一定高度后解的精確度不再提高，只有計算時間不斷提高。 有限元分析可被用來分析比較復雜的、用一般地代數方法無法足夠精確地分析的系統(tǒng)，它可以提供使用其它方法無法提供的結果。在實踐中一般使用電腦來求解在分析時出現的巨量的方程組。 有限元的概念早在幾個世紀前就已產生并得到了應用，例如用多邊形(有限個直線單元)逼近圓來求得圓的周長，但作為一種方法而被提出，則是最近的事。有限元法最初被稱為矩陣近似方法，應用于航空器的結構強度計算，并由于其方便性、實用性和有效性而引起從事力學研究的科學家的濃厚興趣。經過短短數十年的努力，隨著計算機技術的快速發(fā)展和普及，有限元方法迅速從結構工程強度分析計算擴展到幾乎所有的科學技術領域，成為一種豐富多彩、應用廣泛并且實用高效的數值分析方法。 有限元方法與其他求解值問題近似方法的根本區(qū)別在于它的近似性僅限于相對小的子域中。20世紀60年代初首次提出結構力學計算有限元概念的克拉夫(Clough)教授形象地將其描繪為：“有限元法=Rayleigh Ritz法+分片函數”，即有限元法是Rayleigh Ritz法的一種局部化情況。不同于求解(往往是困難的)滿足整個定義域邊界條件的允許函數的Rayleigh Ritz法，有限元法將函數定義在簡單幾何形狀(如二維問題中的三角形或任意四邊形)的單元域上(分片函數)，且不考慮整個定義域的復雜邊界條件，這是有限元法優(yōu)于其他近似方法的原因之一。 　　對于不同物理性質和數學模型的問題，有限元求解法的基本步驟是相同的，只是具體公式推導和運算求解不同。有限元求解問題的基本步驟通常為： 　　第一步：問題及求解域定義：根據實際問題近似確定求解域的物理性質和幾何區(qū)域。 　　第二步：求解域離散化：將求解域近似為具有不同有限大小和形狀且彼此相連的有限個單元組成的離散域，習慣上稱為有限元網絡劃分。顯然單元越小(網絡越細)則離散域的近似程度越好，計算結果也越精確，但計算量及誤差都將增大，因此求解域的離散化是有限元法的核心技術之一。 　　第三步：確定狀態(tài)變量及控制方法：一個具體的物理問題通?？梢杂靡唤M包含問題狀態(tài)變量邊界條件的微分方程式表示，為適合有限元求解，通常將微分方程化為等價的泛函形式。 　　第四步：單元推導：對單元構造一個適合的近似解，即推導有限單元的列式，其中包括選擇合理的單元坐標系，建立單元試函數，以某種方法給出單元各狀態(tài)變量的離散關系，從而形成單元矩陣(結構力學中稱剛度陣或柔度陣)。 　　為保證問題求解的收斂性，單元推導有許多原則要遵循。對工程應用而言，重要的是應注意每一種單元的解題性能與約束。例如，單元形狀應以規(guī)則為好，畸形時不僅精度低，而且有缺秩的危險，將導致無法求解。 　　第五步：總裝求解：將單元總裝形成離散域的總矩陣方程(聯合方程組)，反映對近似求解域的離散域的要求，即單元函數的連續(xù)性要滿足一定的連續(xù)條件?？傃b是在相鄰單元結點進行，狀態(tài)變量及其導數(可能的話)連續(xù)性建立在結點處。 　　第六步：聯立方程組求解和結果解釋：有限元法最終導致聯立方程組。聯立方程組的求解可用直接法、迭代法和隨機法。求解結果是單元結點處狀態(tài)變量的近似值。對于計算結果的質量，將通過與設計準則提供的允許值比較來評價并確定是否需要重復計算。 簡言之，有限元分析可分成三個階段，前處理、處理和后處理。前處理是建立有限元模型，完成單元網格劃分;后處理則是采集處理分析結果，使用戶能簡便提取信息，了解計算結果。 在分析一個物體或系統(tǒng)中的壓力和變形時有限元分析是一種常用的手段，此外它還被用來分析許多其它問題如熱傳導、流體力學和電力學。 有限元分析通常借助計算機軟件完成，著名工程軟件有 MSC NASTRAN，ADINA，LS-DYNA，ANSYS，ABAQUS，2D-sigma等。 2.2 ANSYS軟件的介紹 ANSYS軟件是融結構、流體、電場、磁場、聲場分析于一體的大型通用有限元分析軟件。由世界上最大的有限元分析軟件公司之一的美國ANSYS開發(fā)，它能與多數CAD軟件接口，實現數據的共享和交換，如Pro/Engineer, NASTRAN, Alogor, I－DEAS, AutoCAD等， 是現代產品設計中的高級CAD工具之一。軟件主要包括三個部分：前處理模塊，分析計算模塊和后處理模塊。前處理模塊提供了一個強大的實體建模及網格劃分工具，用戶可以方便地構造有限元模型；分析計算模塊包括結構分析（可進行線性分析、非線性分析和高度非線性分析）、流體動力學分析、電磁場分析、聲場分析、壓電分析以及多物理場的耦合分析，可模擬多種物理介質的相互作用，具有靈敏度分析及優(yōu)化分析能力；后處理模塊可將計算結果以彩色等值線顯示、梯度顯示、矢量顯示、粒子流跡顯示、立體切片顯示、透明及半透明顯示（可看到結構內部）等圖形方式顯示出來，也可將計算結果以圖表、曲線形式顯示或輸出。軟件提供了100種以上的單元類型，用來模擬工程中的各種結構和材料。該軟件有多種不同版本，可以運行在從個人機到大型機的多種計算機設備上，如PC，SGI，HP，SUN，DEC，IBM，CRAY等。 啟動ANSYS，進入歡迎畫面以后，程序停留在開始平臺。從開始平臺（主菜單）可以進入各處理模塊：PREP7（通用前處理模塊），SOLUTION（求解模塊），POST1（通用后處理模塊），POST26（時間歷程后處理模塊）。ANSYS用戶手冊的全部內容都可以聯機查閱。 　用戶的指令可以通過鼠標點擊菜單項選取和執(zhí)行，也可以在命令輸入窗口通過鍵盤輸入。命令一經執(zhí)行，該命令就會在.LOG文件中列出，打開輸出窗口可以看到.LOG文件的內容。如果軟件運行過程中出現問題，查看.LOG文件中的命令流及其錯誤提示，將有助于快速發(fā)現問題的根源。.LOG 文件的內容可以略作修改存到一個批處理文件中，在以后進行同樣工作時，由ANSYS自動讀入并執(zhí)行，這是ANSYS軟件的第三種命令輸入方式。這種命令方式在進行某些重復性較高的工作時，能有效地提高工作速度。 ANSYS的各處理模塊具體內容有： ①前處理模塊PREP7： ●實體建模 ANSYS程序提供了兩種實體建模方法：自頂向下與自底向上。 自頂向下進行實體建模時，用戶定義一個模型的最高級圖元，如球、棱柱，稱為基元，程序則自動定義相關的面、線及關鍵點。用戶利用這些高級圖元直接構造幾何模型，如二維的圓和矩形以及三維的塊、球、錐和柱。無論使用自頂向下還是自底向上方法建模，用戶均能使用布爾運算來組合數據集，從而“雕塑出”一個實體模型。ANSYS程序提供了完整的布爾運算，諸如相加、相減、相交、分割、粘結和重疊。在創(chuàng)建復雜實體模型時，對線、面、體、基元的布爾操作能減少相當可觀的建模工作量。ANSYS程序還提供了拖拉、延伸、旋轉、移動、延伸和拷貝實體模型圖元的功能。附加的功能還包括圓弧構造、切線構造、通過拖拉與旋轉生成面和體、線與面的自動相交運算、自動倒角生成、用于網格劃分的硬點的建立、移動、拷貝和刪除。 自底向上進行實體建模時，用戶從最低級的圖元向上構造模型，即：用戶首先定義關鍵點，然后依次是相關的線、面、體。 ●網格劃分 ANSYS程序提供了使用便捷、高質量的對CAD模型進行網格劃分的功能。包括四種網格劃分方法：延伸劃分、映像劃分、自由劃分和自適應劃分。延伸網格劃分可將一個二維網格延伸成一個三維網格。映像網格劃分允許用戶將幾何模型分解成簡單的幾部分，然后選擇合適的單元屬性和網格控制，生成映像網格。ANSYS程序的自由網格劃分器功能是十分強大的，可對復雜模型直接劃分，避免了用戶對各個部分分別劃分然后進行組裝時各部分網格不匹配帶來的麻煩。自適應網格劃分是在生成了具有邊界條件的實體模型以后，用戶指示程序自動地生成有限元網格，分析、估計網格的離散誤差，然后重新定義網格大小，再次分析計算、估計網格的離散誤差，直至誤差低于用戶定義的值或達到用戶定義的求解次數。 ②求解模塊SOLUTION 前處理階段完成建模以后，用戶可以在求解階段獲得分析結果。 點擊快捷工具區(qū)的SAVE_DB將前處理模塊生成的模型存盤，退出 Preprocessor，點擊實用菜單項中的Solution，進入分析求解模塊。在該階段，用戶可以定義分析類型、分析選項、載荷數據和載荷步選項，然后開始有限元求解。 ANSYS軟件提供的分析類型如下： 1.結構靜力分析 用來求解外載荷引起的位移、應力和力。靜力分析很適合求解慣性和阻尼對結構的影響并不顯著的問題。ANSYS程序中的靜力分析不僅可以進行線性分析，而且也可以進行非線性分析，如塑性、蠕變、膨脹、大變形、大應變及接觸分析。 2.結構動力學分析 結構動力學分析用來求解隨時間變化的載荷對結構或部件的影響。與靜力分析不同，動力分析要考慮隨時間變化的力載荷以及它對阻尼和慣性的影響。ANSYS可進行的結構動力學分析類型包括：瞬態(tài)動力學分析、模態(tài)分析、諧波響應分析及隨機振動響應分析。 3.結構非線性分析 結構非線性導致結構或部件的響應隨外載荷不成比例變化。ANSYS程序可求解靜態(tài)和瞬態(tài)非線性問題，包括材料非線性、幾何非線性和單元非線性三種。 4.動力學分析 ANSYS程序可以分析大型三維柔體運動。當運動的積累影響起主要作用時，可使用這些功能分析復雜結構在空間中的運動特性，并確定結構中由此產生的應力、應變和變形。 5.熱分析 程序可處理熱傳遞的三種基本類型：傳導、對流和輻射。熱傳遞的三種類型均可進行穩(wěn)態(tài)和瞬態(tài)、線性和非線性分析。熱分析還具有可以模擬材料固化和熔解過程的相變分析能力以及模擬熱與結構應力之間的熱－結構耦合分析能力。 6.電磁場分析 主要用于電磁場問題的分析，如電感、電容、磁通量密度、渦流、電場分布、磁力線分布、力、運動效應、電路和能量損失等。還可用于螺線管、調節(jié)器、發(fā)電機、變換器、磁體、加速器、電解槽及無損檢測裝置等的設計和分析領域。 7.流體動力學分析 ANSYS流體單元能進行流體動力學分析，分析類型可以為瞬態(tài)或穩(wěn)態(tài)。分析結果可以是每個節(jié)點的壓力和通過每個單元的流率。并且可以利用后處理功能產生壓力、流率和溫度分布的圖形顯示。另外，還可以使用三維表面效應單元和熱－流管單元模擬結構的流體繞流并包括對流換熱效應。 8.聲場分析 程序的聲學功能用來研究在含有流體的介質中聲波的傳播，或分析浸在流體中的固體結構的動態(tài)特性。這些功能可用來確定音響話筒的頻率響應，研究音樂大廳的聲場強度分布，或預測水對振動船體的阻尼效應。 9.壓電分析 用于分析二維或三維結構對AC（交流）、DC（直流）或任意隨時間變化的電流或機械載荷的響應。這種分析類型可用于換熱器、振蕩器、諧振器、麥克風等部件及其它電子設備的結構動態(tài)性能分析?？蛇M行四種類型的分析：靜態(tài)分析、模態(tài)分析、諧波響應分析、瞬態(tài)響應分析 ③后處理模塊POST1和POST26 ANSYS軟件的后處理過程包括兩個部分：通用后處理模塊POST1和時間歷程后處理模塊POST26。通過友好的用戶界面，可以很容易獲得求解過程的計算結果并對其進行顯示。這些結果可能包括位移、溫度、應力、應變、速度及熱流等，輸出形式可以有圖形顯示和數據列表兩種。 ●通用后處理模塊POST1 點擊實用菜單項中的“General Postproc”選項即可進入通用后處理模塊。這個模塊對前面的分析結果能以圖形形式顯示和輸出。例如，計算結果（如應力）在模型上的變化情況可用等值線圖表示，不同的等值線顏色，代表了不同的值（如應力值）。濃淡圖則用不同的顏色代表不同的數值區(qū)（如應力范圍），清晰地反映了計算結果的區(qū)域分布情況。 ●時間歷程響應后處理模塊POST26 點擊實用菜單項中的TimeHist Postpro選項即可進入時間歷程響應后處理模塊。這個模塊用于檢查在一個時間段或子步歷程中的結果，如節(jié)點位移、應力或支反力。這些結果能通過繪制曲線或列表查看。繪制一個或多個變量隨頻率或其它量變化的曲線，有助于形象化地表示分析結果。另外，POST26還可以進行曲線的代數運算。 ④ANSYS軟件的質量認證 ANSYS軟件是第一個通過ISO9001質量認證的大型分析設計類軟件，是美國機械工程師協會（ASME）、美國核安全局（NQA）及近二十種專業(yè)技術協會認證的標準分析軟件。在國內第一個通過了中國壓力容器標準化技術委員會認證并在國務院十七個部委推廣使用。 3 模型及系統(tǒng)分析 本論文所要研究的是大型龍門銑床主軸滑枕有限元分析。首先是介紹ANSYS實體建模的基本方法。有限元分析是針對特定的模型而進行的，因此，必須建立一個有物理原型的準確的數學模型。一般情況下，有限元模型包含所有的節(jié)點、單元、材料屬性、實常數、邊界條件、荷載條件以及用這些條件來定義的物理場特征。建模與網格劃分是有限元分析的基礎。而后是對該模型進行動力學分析。結構動力學分析是用來求解隨時間變化的荷載對結構或部件的影響。與靜力學分析不同，動力學分析要考慮隨時間變化的力荷載以及它對阻尼和慣性的影響。ANSYS可進行的結構動力學分析類型包括：模態(tài)分析、瞬態(tài)動力學分析、諧波響應分析及隨機振動響應分析。在實際工程中，結構除了受永久性荷載作用外，還會受到或大或小的動荷載作用。通常，當荷載變化緩慢，變化周期遠大于結構的自振周期時，其動力響應是很小的，可將它作為靜力荷載處理。反之，對于那些變化激烈，動力作用顯著的荷載必須考慮結構的動力響應。常見的動力荷載有簡諧荷載、沖擊荷載、突加荷載、快速移動荷載和隨機荷載等。在本文中分別對實體模型的模態(tài)頻率、模態(tài)振型、瞬態(tài)動力分析及簡諧激勵下的諧響應進行分析，在每一項中，又分別對不同參數的模型（即模型加強筋厚度分別為15mm、20mm、25mm）進行分析和比較，得出結論。 3.1 模型描述 首先介紹ANSYS建模的基本方法。由節(jié)點和單元構成的有限元模型與結構系統(tǒng)的幾何外形是基本一致的。有限元模型的建立可分為直接法和間接法。直接法是根據結構的幾何外形建立節(jié)點和單元，因此它只適用于簡單的結構系統(tǒng)。間接法適用于節(jié)點及單元數目較多的復雜幾何外形的結構系統(tǒng)，該方法通過點、線、面和體，先建立有限元模型，再進行實體網格劃分，以完成有限元模型的建立。 實體模型的建立有下列3種方法： ●自底向上法（Bottom-up Method） 由建立最低單元的點到最高單元的體，即先建立點，再由點連成線，然后由線組合成面，最后由面組合成體。 ●自頂向下（Top-down Method） 此方法直接建立較高單元對象，同時將一起產生所對應的較低單元對象，對象單元商低順序依次為體、面、線及點。所謂布爾運算就是對對象進行相互加、減或者組合等操作。 ●混合使用前兩種方法 通常，可結合前兩種方法進行綜合運用，但應考慮到要獲得什么樣的有限元模型，即在進行網格劃分時，是要產生自由網格劃分還是對應網格劃分。自由網格劃分時，實體模型的建立比較簡單，只要所有的面或體能接合成一個體就可以；對應網格劃分時，平面結構一定要由四邊形或三邊形的面相接而成，立體結構一定要由六面體相接而成。 然后說明本論文中龍門銑床主軸滑枕模型的設計參數及建立方法。 龍門銑床主軸滑枕工作狀態(tài)為： 1、 輸入電機（交流伺服）：功率37kW，額定轉速1500rpm，額定扭矩235.5Nm； 2、 傳動箱傳動比：高擋1.4，低擋6.776； 3、 主軸轉速220rpm，最大輸出扭矩1595Nm； 4、 使用刀具直徑200mm，刀盤齒數10，刀盤厚度70mm； 5、 主軸直徑128.57mm； 6、 滑枕與滑座之間靠六面導軌、壓板約束，長度為1500mm； 7、 滑枕材料為球墨鑄鐵QT600-3A; 8、 滑枕全長2500mm，滑枕伸出滑座1000mm，使用直徑200mm的刀盤，承受切削扭矩1595Nm，環(huán)境溫度變化15℃，電機與傳動箱發(fā)熱的溫升為30℃。 有限元模型建立的好壞關系到分析計算的準確性和計算成本。一般而言，可采用三維CAD軟件建立待分析結構的三維實體模型，然后通過初始圖形交換規(guī)范 (Initial Graphics Exchange Specification，IGES) 導入有限元分析軟件中，然而，對復雜模型而言，該方法基本上不能成功。因此，使用較多的利用有限元分析軟件本身的建模功能，忽略細節(jié)信息，獲得具有幾何近似性的求解模型。本論文采用ANSYS軟件，用自底向上法（Bottom-up Method），即利用點、線、面、體的模型生成方法，編寫求解命令流文件，生成滑枕模型，并求出滑枕的變形量及其分布情況,對滑枕的結構進行動力學分析。 滑枕主軸組合模型：有限元分析實體模型如圖1(a)所示，網格剖分如圖1(b)所示。用SOLID45單元求解，劃分單元數目68939個。 （a）實體模型 （b）網格剖分 圖1 滑枕主軸組合模型 3.2 主軸滑枕靜力分析 靜力分析計算在固定不變的荷載作用下結構的效應，它不考慮慣性和阻尼的影響，例如結構受隨時間變化荷載的情況?？墒?，靜力分析可以計算那些固定不變的慣性荷載對結構影響（例如重力和離心力），以及那些可以近似為等價靜力作用的隨時間變化荷載（例如通常在許多建筑規(guī)范中所定義的等價靜力風載和地震荷載）。 靜力分析用于計算由那些不包括慣性和阻尼效應的荷載作用于結構或部件上所引起的位移、應力、應變和力。固定不變的荷載和響應是一種假定，即假定荷載和結構的響應隨時間的變化非常緩慢。靜力分析所施加的荷載包括：外部施加的作用力和壓力、穩(wěn)態(tài)的慣性力（如重力和離心力）、位移荷載及溫度荷載。靜力分析既可以是線性的也可以是非線性的。非線性靜力分析包括所有的非線性類型：大變形、塑性、蠕變、應力剛化、接觸（間隙）單元以及超彈性單元等。 3.2.1 加強筋厚度為15mm時主軸滑枕靜力分析 工況：滑枕伸出滑座1000mm，約束住滑槽3個面兩個方向移動，同時約束住滑槽沿Z方向的移動，約束加在側面上一個小區(qū)域。 采用SOLID45單元求解，劃分單元數目68939個。將銑削力作為動態(tài)力，乘以動載荷系數10，則實加力為：10×F1＝7975N×10＝79750N。 圖2(a)、2(c)所示為變形圖，最大變形為：0.14122mm。圖2(b)、2(d)所示為應力云圖；最大應力為集中力作用點，而滑塊上應力不超過1.402MPa，該結構強度足夠。 (b) X方向應力云圖 (a) X方向變形圖 (d) 總應力云圖 (c) 總變形圖 圖2 變形和應力云圖 結論：工況求解獲得最大變形為0.14122mm，應力求解結果遠遠小于鋼材屈服極限306Mpa，結果強度足夠。 3.2.2 加強筋厚度為20mm時主軸滑枕靜力分析 工況：滑枕伸出滑座1000mm，約束住滑槽3個面兩個方向移動，同時約束住滑槽沿Z方向的移動，約束加在側面上一個小區(qū)域。 采用SOLID45單元求解，劃分單元數目68939個。將銑削力作為動態(tài)力，乘以動載荷系數10，則實加力為：10×F1＝7975N×10＝79750N。 圖3(a)、3(c)所示為變形圖，最大變形為：0.14027mm。圖3(b)、3(d)所示為應力云圖；最大應力為集中力作用點，而滑塊上應力不超過1.412MPa，該結構強度足夠。 (b) X方向應力云圖 (a) X方向變形圖 (d) 總應力云圖 (c) 總變形圖 圖3 變形和應力云圖 結論：工況求解獲得最大變形為0.14027mm,應力求解結果遠遠小于鋼材屈服極限306Mpa，結果強度足夠。 3.2.3 加強筋厚度為25mm時主軸滑枕靜力分析 工況：滑枕伸出滑座1000mm，約束住滑槽3個面兩個方向移動，同時約束住滑槽沿Z方向的移動，約束加在側面上一個小區(qū)域。 采用SOLID45單元求解，劃分單元數目68939個。將銑削力作為動態(tài)力，乘以動載荷系數10，則實加力為：10×F1＝7975N×10＝79750N。 圖4(a)、4(c)所示為變形圖，最大變形為：1.735mm。圖4(b)、4(d)所示為應力云圖；最大應力為集中力作用點，而滑塊上應力達到576.261MPa，該結構強度不夠。 (b) X方向應力云圖 (a) X方向變形圖 (d) 總應力云圖 (c) 總變形圖 圖4 變形和應力云圖 結論：工況求解獲得最大變形為1.735mm，應力求解結果大于鋼材屈服極限306Mpa，結果強度不夠。 3.2.4 結論 由以上3點可以看出，在相同工況下，加強筋厚度為15mm時，求解獲得最大變形為0.14122mm,應力求解結果1.402MPa遠遠小于鋼材屈服極限306Mpa，強度足夠；加強筋厚度為20mm時，求解獲得最大變形為0.14027mm,應力求解結果1.412MPa遠遠小于鋼材屈服極限306Mpa，強度足夠；加強筋厚度為25mm時，求解獲得最大變形為1.735mm,應力求解結果576.261MPa大于鋼材屈服極限306Mpa，強度不夠。 相比較可以看出，加強筋厚度為15mm、20mm時，滑枕強度都足夠。而且可以發(fā)現，一方面，變形最大的位置在滑塊上下底面靠近刀盤的部位，應該采取結構措施，減少變形量；另一方面，當單側力加載時，變形量普遍偏大，其原因在于滑塊為彎扭組合變形。同時，假如結構載荷對稱性好，對減少滑塊變形有利。 3.3 主軸滑枕動力學分析 針對該結構具有穩(wěn)定的激勵源（電機），對結構進行模態(tài)分析、瞬態(tài)動力學分析和諧響應分析，求出系統(tǒng)的模態(tài)頻率、模態(tài)振型、瞬態(tài)動力響應及諧響應，考察結構抗振動性。針對該結構主軸剛性足夠的特點，并結合排查問題的需要，在對結構進行合理簡化的基礎上，分析主軸滑枕在力矩作用下的變形和應力發(fā)布狀況。 3.3.1 主軸滑枕模態(tài)頻率 模態(tài)分析主要用于確定設計結構或機械部件的振動特性，即結構的固有頻率和各階振型，它們是承受動態(tài)荷載結構設計中的重要參數，也是其他各類型動力學分析（如諧響應分析、瞬態(tài)動力分析等）的基礎。 ANSYS的模態(tài)分析是一個線性分析，因此在模態(tài)分析中，任何非線性特性，例如塑性和接觸單元，即使定義了也會被忽略。ANSYS求解模態(tài)分析的方法有子空間法、分塊Lanczos法、縮減法、非對稱法、阻尼法、QR阻尼法以及PowerDynamics法等。在很多場合，模態(tài)分析都起到舉足輕重的作用，例如結構都必須避免共振。進行模態(tài)分析之后，可以了解結構的固有頻率和振動型式，并對此采取必要的措施，從而避免在使用中由于共振的因素造成不必要的損失。 工況：滑枕伸出滑座1000mm，約束住滑槽3個面兩個方向移動，同時約束住滑槽沿Z方向的移動，約束加在側面上一個小區(qū)域。 (1)加強筋厚度為15mm時主軸滑枕模態(tài)頻率 前10階模態(tài)頻率求解結果如表1所示： 表1 前10階頻率(Hz) 階數 頻率（HZ） 階數 頻率（HZ） 1 163.65 6 807.20 2 168.27 7 841.89 3 392.98 8 1173.1 4 418.31 9 1194.9 5 766.59 10 1256.9 結論：基頻(第1階固有頻率)為163.65Hz。遠高于電機工頻激勵（1500rpm/60m=25Hz）。因此，滑枕本身不可能由于傳遞系統(tǒng)激勵作用而產生過大的振動。 (2)加強筋厚度為20mm時主軸滑枕模態(tài)頻率 前10階模態(tài)頻率求解結果如表2所示： 表2 前10階頻率(Hz) 階數 頻率（HZ） 階數 頻率（HZ） 1 163.90 6 808.94 2 168.69 7 842.57 3 392.83 8 1174.3 4 421.14 9 1201.2 5 768.55 10 1266.8 結論：基頻(第1階固有頻率)為163.90Hz。遠高于電機工頻激勵（1500rpm/60m=25Hz）。因此，滑枕本身不可能由于傳遞系統(tǒng)激勵作用而產生過大的振動。 (3)加強筋厚度為25mm時主軸滑枕模態(tài)頻率 前10階模態(tài)頻率求解結果如表3所示： 階數 頻率（HZ） 階數 頻率（HZ） 1 163.98 6 810.27 2 168.70 7 843.29 3 392.83 8 1171.9 4 423.89 9 1205.6 5 768.43 10 1275.8 表3 前10階頻率(Hz) 結論：基頻(第1階固有頻率)為163.98Hz。遠高于電機工頻激勵（1500rpm/60m=25Hz）。因此，滑枕本身不可能由于傳遞系統(tǒng)激勵作用而產生過大的振動。 (4)結論 比較表1、表2、表3中所得模態(tài)頻率，可以看出，工況基頻分別為163.65Hz、163.90 Hz、163.98 Hz，并且3個基頻相差不大，都遠高于電機工頻激勵頻率（1500rpm/60m=25Hz）。在相同工況下，隨著加強筋厚度的增大，相同階數下主軸滑枕的固有頻率也在相應地增大。通過模態(tài)分析可知：該結構不會因為電機工頻率激勵而產生共振。 3.3.2 主軸滑枕模態(tài)振型 模態(tài)分析是用來確定結構的振動特性的一種技術，通過它可以確定自然頻率、振型和振型參與系數（即在特定方向上某個振型在多大程度上參與了振動）。 進行模態(tài)分析可以使結構設計避免共振或以特定頻率進行振動（例如揚聲器），可以使工程師認識到結構對于不同類型的動力載荷是如何響應的，也有助于在其他動力分析中估算求解控制參數（如時間步長）。由于結構的振動特性決定結構對于各種動力載荷的響應情況，所以在準備進行其他動力分析之前首先要進行模態(tài)分析。 工況：滑枕伸出滑座1000mm，約束住滑槽3個面兩個方向移動，同時約束住滑槽沿Z方向的移動，約束加在側面上一個小區(qū)域。 (1)加強筋厚度為15mm時主軸滑枕模態(tài)振型 前6階振型如圖5所示，以彎扭組合變形為主，表明：一旦發(fā)生結構上共振或自激振動，將產生過大的彎扭組合變形。 (b) 第2階振型 (a) 第1階振型 (d) 第4階振型 (c) 第3階振型 (e) 第5階振型 (f) 第6階振型 圖5 懸伸1000mm時滑枕主軸組合模型前6階振型 (2)加強筋厚度為20mm時主軸滑枕模態(tài)振型 前6階振型如圖6所示，以彎扭組合變形為主，表明：一旦發(fā)生結構上共振或自激振動，將產生過大的彎扭組合變形。 (b) 第2階振型 (a) 第1階振型 (d) 第4階振型 (c) 第3階振型 (e) 第5階振型 (f) 第6階振型 圖6 懸伸1000mm時滑枕主軸組合模型前6階振型 (3)加強筋厚度為25mm時主軸滑枕模態(tài)振型 前6階振型如圖7所示，以彎扭組合變形為主，表明：一旦發(fā)生結構上共振或自激振動，將產生過大的彎扭組合變形。 (b) 第2階振型 (a) 第1階振型 (d) 第4階振型 (c) 第3階振型 (e) 第5階振型 (f) 第6階振型 圖7 懸伸1000mm時滑枕主軸組合模型前6階振型 (4)結論 由前面模態(tài)頻率分析可知，不同加強筋厚度下，工況基頻分別為163.65Hz、163.90 Hz、163.98Hz，遠高于電機工頻激勵頻率（1500rpm/60m=25Hz）。通過模態(tài)分析可知：該結構不會因為電機工頻率激勵而產生共振，而且隨著加強筋厚度的增大，主軸滑枕的基頻也相應增大，而相應的振型變形卻相應減小。 3.3.3 主軸滑枕瞬態(tài)動力學分析 瞬態(tài)動力學分析完全不同于靜力學分析，它分析結構的運動特性，包括各種運動或爆破問題的分析。動力學分析是結構分析的一大分支，有其自身的一套理論體系和求解方法，ANSYS的瞬態(tài)動力學分析中除了其自身的分析模塊外，還有LS-DYNA模塊。 瞬態(tài)動力分析又稱為時間歷程分析，用來確定結構在隨時間變化的荷載作用下的動力響