基于ANSYS-PDS的數(shù)控導(dǎo)軌磨床主軸可靠性分析

胡軍

?

基于ANSYS-PDS的數(shù)控導(dǎo)軌磨床主軸可靠性分析

胡軍

（泰州市產(chǎn)品質(zhì)量監(jiān)督檢驗(yàn)院，江蘇 泰州 225300）

根據(jù)可靠性分析理論，運(yùn)用有限元分析軟件的概率設(shè)計(jì)功能（模塊），對數(shù)控導(dǎo)軌磨床的主軸進(jìn)行可靠性分析。分析時(shí)按照工程實(shí)際，將主軸的幾何尺寸、載荷、強(qiáng)度極限等參數(shù)作為隨機(jī)輸入變量，采用蒙特卡羅法進(jìn)行抽樣計(jì)算，可以得到比傳統(tǒng)設(shè)計(jì)方法精度更高的主軸失效概率。以MK5220數(shù)控導(dǎo)軌磨床主軸為例進(jìn)行計(jì)算，結(jié)果表明利用ANSYS/PDS模塊對磨床主軸進(jìn)行可靠性分析是可行的。

數(shù)控導(dǎo)軌磨床；主軸；蒙特卡羅法；可靠性

數(shù)控導(dǎo)軌磨床的關(guān)鍵結(jié)構(gòu)是傳動(dòng)結(jié)構(gòu)，主軸是傳動(dòng)結(jié)構(gòu)的關(guān)鍵零部件，它的工作狀況直接影響數(shù)控導(dǎo)軌磨床的加工質(zhì)量，制約著磨床的發(fā)展。主軸直接承受來自電機(jī)的扭矩，而且在實(shí)際工況下，磨床的加工環(huán)境是惡劣的，主軸很容易產(chǎn)生各種故障，甚至損壞。因此，為提高磨床主軸的性能，可以對其進(jìn)行可靠性分析?？煽啃苑治龇椒梢詫⒂绊懥慵ぷ餍阅艿淖兞浚ㄈ绯叽?、載荷）視為不確定量，從而分析這些變量對零件可靠性的影響，進(jìn)而提高零部件設(shè)計(jì)的準(zhǔn)確性和可靠性[1-2]。本文就是根據(jù)文獻(xiàn)[3]運(yùn)用ANSYS-PDS模塊對實(shí)際工況下的磨床主軸進(jìn)行可靠性分析。

1 磨床主軸可靠性分析過程及理論

1.1 基本原理

可靠性工程是綜合型的系統(tǒng)工程，產(chǎn)品的可靠性與其設(shè)計(jì)、制造、使用等環(huán)節(jié)密切相關(guān)[4]。其可靠性設(shè)計(jì)是以應(yīng)力—強(qiáng)度干涉模型為基礎(chǔ)。一般情況下，用表示施加在產(chǎn)品上的物理量，如載荷，稱為產(chǎn)品所受的應(yīng)力；用表示產(chǎn)品能夠承受這種應(yīng)力的能力，稱為產(chǎn)品的強(qiáng)度。要使產(chǎn)品在規(guī)定的時(shí)間內(nèi)可靠的工作，就必須滿足＝－≥0。其中，應(yīng)力、強(qiáng)度都是某些變量的函數(shù)。產(chǎn)品的可靠度就是值大于等于零時(shí)的概率。在磨床主軸可靠性分析中，將實(shí)際工況下的許用應(yīng)力值與最大應(yīng)力值的差值定義為，那么，主軸的可靠度為＝（≥0)[5]。

1.2 主軸可靠性分析基本過程

本文運(yùn)用ANSYS-PDS模塊對磨床主軸進(jìn)行可靠性分析，根據(jù)GB/T 9061－2006中的相關(guān)要求，對磨床主軸在最大扭矩工作狀況下作可靠性分析。具體分析流程如下[6]：

（1）創(chuàng)建主軸的可靠性分析文件，包括建立主軸模型，施加載荷、約束及求解等過程；

（2）對主軸進(jìn)行可靠性分析，其中包括定義隨機(jī)變量，選擇仿真循環(huán)次數(shù)等過程；

（3）可靠度的求取。

1.3 主軸可靠性分析方法

在ANSYS-PDS分析過程中，常用的概率設(shè)計(jì)方法有蒙特卡羅法和響應(yīng)面法。蒙特卡羅法對有限元有非常好的適應(yīng)性，能模擬實(shí)際問題的真實(shí)特征[7]。本文采用蒙特卡羅法進(jìn)行可靠性分析。

2 磨床主軸的可靠性分析

2.1 主軸受力分析

主傳動(dòng)系統(tǒng)工作情況如下：電機(jī)（額定功率為30 kW，額定轉(zhuǎn)速為1470 r/min，經(jīng)傳動(dòng)比為150:290的帶輪帶動(dòng)主軸運(yùn)動(dòng)，主軸另一端連接砂輪，砂輪尺寸為600×150×305 mm。數(shù)控導(dǎo)軌磨床在運(yùn)行過程中，其主軸轉(zhuǎn)速范圍為600～1100 r/min，磨削深度為0.02～0.1 mm。

主軸模型簡化圖如圖1所示。主軸材料為45#鋼，許用扭轉(zhuǎn)切應(yīng)力為30 MPa，安全系數(shù)?。?。

圖1 主軸模型簡化圖

根據(jù)文獻(xiàn)[3]中的有關(guān)規(guī)定，機(jī)床在檢驗(yàn)時(shí)需做負(fù)荷試驗(yàn)，包括機(jī)床主傳動(dòng)系統(tǒng)扭矩試驗(yàn)，主傳動(dòng)系統(tǒng)最大功率試驗(yàn)等。由于磨床主軸工作時(shí)主要受扭矩作用，且主軸在最大功率下運(yùn)行時(shí)所受扭矩小于最大扭矩試驗(yàn)中主軸承受的扭矩，因此本文主要研究最大扭矩條件下主軸的可靠性。

按照該標(biāo)準(zhǔn)中最大扭矩的近似計(jì)算法，最大扭矩的公式為：

式中：為扭矩，N·m；為切削時(shí)電動(dòng)機(jī)的輸入功率，kW；0為機(jī)床裝有工件時(shí)的空運(yùn)轉(zhuǎn)功率，kW；為機(jī)床主軸轉(zhuǎn)速，r/min。

計(jì)算得max＝540 N·m。按125%計(jì)算得出進(jìn)行可靠性分析時(shí)＝680 N·m。

磨床在一般情況下磨削時(shí)，磨削深度a＝0.05 mm，主軸轉(zhuǎn)速為750 r/min，工件速度為20 m/s。由于磨削力受多方面因素影響，沒有準(zhǔn)確的公式計(jì)算，在本文中運(yùn)用平面磨削力的經(jīng)驗(yàn)公式計(jì)算得出切向磨削力[8]：

按照F/F＝3.8求出磨削時(shí)所受的徑向磨削力。經(jīng)計(jì)算得：

F＝282 N

F＝1071.6 N

2.2 主軸建模、網(wǎng)格劃分及加載求解

采用ANSYS參數(shù)化設(shè)計(jì)語言APDL對主軸進(jìn)行建模，建模過程中在不影響主軸受力的條件下對其進(jìn)行了一定的簡化，忽略了主軸的倒角等[9]。選用4節(jié)點(diǎn)的平面單元模型Quad 4 node 182先對主軸截面進(jìn)行網(wǎng)格劃分，再旋轉(zhuǎn)得到主軸的網(wǎng)格劃分圖，如圖2所示。然后對主軸分別施加邊界條件和載荷，再進(jìn)行求解。其中扭矩及徑向力均施加在對應(yīng)節(jié)點(diǎn)上。最大扭矩下，施加約束及載荷后求得的主軸應(yīng)力云圖如圖3、圖4所示。

圖2 主軸有限元模型

圖3 主軸最大扭矩試驗(yàn)應(yīng)力圖

圖4 主軸正常切削應(yīng)力圖

2.3 可靠性計(jì)算與分析

分別整理保存分析求解過程的命令流。寫入可靠性分析文件。將等效應(yīng)力按降序排列并提取最大等效應(yīng)力值“maxstress”，然后定義極限狀態(tài)方程，即求出強(qiáng)度與最大應(yīng)力值的差值，將之與之前ANSYS應(yīng)力分析命令一起構(gòu)成概率分析文件。進(jìn)入PDS模塊，讀取保存的可靠性分析宏文件，選擇、定義隨機(jī)輸入變量（見表1）和輸出變量服從的分布類型及其參數(shù)[10]。

表1 隨機(jī)輸入變量統(tǒng)計(jì)表

2.4 可靠性分析結(jié)果

選擇蒙特卡羅法進(jìn)行概率設(shè)計(jì)，選擇拉丁抽樣法，模擬樣本次數(shù)為60次，通過計(jì)算可得到主軸在許用扭轉(zhuǎn)切應(yīng)力為20 MPa，置信度為95%時(shí)，主軸在最大扭矩試驗(yàn)時(shí)可靠度為85.97%。

3 結(jié)論

本文采用可靠性分析理論及ANSYS/PDS（上接第26頁）模塊，根據(jù)工程實(shí)際及國家標(biāo)準(zhǔn)的相關(guān)規(guī)定，采用蒙特卡羅法，以MK5220數(shù)控導(dǎo)軌磨床主軸為例，進(jìn)行可靠性分析，分析計(jì)算得到了其在最大扭矩下的可靠度。通過ANSYS/PDS定量的計(jì)算出磨床主軸在惡劣工況下的可靠性。同時(shí)這種分析方法具有普遍適用性，可以對類似機(jī)械零件進(jìn)行可靠性分析，在工程應(yīng)用中，具有一定的實(shí)用價(jià)值。

[1]余偉煒，高炳軍. ANSYS在機(jī)械與化工裝備中的應(yīng)用[M]. 北京：中國水利水電出版社，2006.

[2]Reh S，Lethbridge P，Ostergaard D. Quality-based design with probabilistic methods[J]. ANSYS Solutions，2002，2（2）：20-22.

[3]GB/T 9061－2006，金屬切削機(jī)床通用技術(shù)條件[S].

[4]劉惟信. 機(jī)械可靠性設(shè)計(jì)[M]. 北京：清華大學(xué)出版社，1996.

[5]施衛(wèi)東，顏品蘭，蔣小平，等. 基于ANSYS/PDS的泵軸可靠性分析[J]. 排灌機(jī)械，2008，26（5）：1-4.

[6]博弈創(chuàng)作室. ANSYS 9.0經(jīng)典產(chǎn)品高級(jí)分析技術(shù)與實(shí)例詳解 [M]. 北京：中國水利水電出版社，2006.

[7]小颯工作室. 最新經(jīng)典ANSYS及Workbench教程[M]. 北京：電子工業(yè)出版社，2004.

[8]陳章燕. 磨削力的試驗(yàn)研究[J]. 北京工業(yè)大學(xué)學(xué)報(bào)，1990，16（4）：98-106.

[9]劉永峰，孫宇，張新洲. 基于ANSYS-PDS的高速?zèng)_床曲軸可靠性分析[J]. 鍛壓技術(shù)，2012，37（5）：155-157，184.

[10]何雪浤，蔣彥收，孫志禮. 基于ANSYS/PDS模塊的高速電主軸剛度可靠性分析[J]. 機(jī)電工程，2011，28（6）：645-647，69.

CNC Guideway Grinding Spindle Reliability Analysis Based on ANSYS/PDS Module

HU Jun

( Product Quality Supervising and Inspecting Institute of Taizhou City, Taizhou 225300, China)

According to the theory of reliability analysis, using probability design function of ANSYS (ANSYS/PDS module) to analyze reliability of CNC guideway grinding spindle. According to the engineering practice, geometry size, load, ultimate strength of the spindle as the random input variables, the spindle failure probability can be calculated by the Monte Carlo method. This paper take MK5220 CNC guideway grinding spindle as example, the result indicates that it is feasible that grinding spindle reliability analysis by ANSYS/PDS.

CNC guideway grinding；spindle；monte carlo method；ANSYS-PDS；reliability

TP202+.1

A

1006－0316 (2018) 04－0024－04

10.3969/j.issn.1006-0316.2018.04.006

2017－07－31

胡軍（1989－），男，江蘇靖江人，碩士研究生，助理工程師，主要從事機(jī)床等機(jī)械產(chǎn)品的檢驗(yàn)工作。