趙 軍,孟欣佳,張立香
(1.山西機電職業技術學院,山西 長治046011;2.中國農業大學工學院,北京100083)
電動沖擊扳手是用于螺栓、螺母裝卸機械化的高效工具,廣泛應用于機車、建筑、化工、冶金、坦克等行業。芯軸作為電動沖擊扳手的核心部件,芯軸的失效將直接導致扳手無法工作。芯軸的傳統設計是采用材料力學的方法來完成的,而優化和可靠性分析都是靠一些經驗來完成。為此在芯軸的設計過程中不免要造成尺寸的偏大及優質材料的浪費。隨著計算機仿真技術的迅速發展,采用計算機進行合理的可靠性分析便顯得十分方便[1]。
以軌枕螺栓預緊用電動扳手的沖擊結構中的芯軸為研究對象。電動扳手的沖擊結構如圖1所示。

該電動沖擊扳手是通過電動機、減速器帶動沖擊機構運作,從而完成鐵路軌枕螺栓的裝卸工作。它將行星齒輪減速機構作為主傳動機構,可保證沖擊扳手體積小、重量輕,使電動扳手傳遞功率大、扭矩輸出較穩定、沖擊小,并能較好地實現對電動扳手扭矩和轉速的控制。如圖1所示,該沖擊機構主要由芯軸、沖擊頭、滾珠、主壓力彈簧、推力滾珠軸承、沖擊桿等組成。其中芯軸采用V型槽結構,沖擊頭采用人字槽結構并帶有兩個凸爪。沖擊機構的作用主要有兩個,一是將芯軸旋轉的動能轉換成沖擊頭的沖擊力,二是確保芯軸上的負載力矩在規定的范圍內。
軌枕螺栓電動扳手芯軸結構如圖2所示。

本文采用ANSYS有限元分析軟件直接建立芯軸的實體模型,為求解方便,在保證不影響分析結果的前提下,忽略倒角、圓角以及人字槽以簡化模型,提高模型質量,減少計算時間。建模時,芯軸的回轉軸線與ANSYS笛卡爾坐標系的Z軸重合。運用ANSYS參數化建模模塊,將芯軸的各軸段直徑按參數處理,并對其進行初始化,采用的單位制為國際單位制。
根據芯軸的實體模型特性和受載情況,定義材料的彈性模量為210 GPa,泊松比為0.3,由于不考慮芯軸的慣性力,在此不定義材料的密度。單元選用SOLID185單元,該單元是3維8節點固體結構單元,每個節點有3個沿著x、y、z方向平移的自由度。劃分網格時,采用自由網格劃分的方式,單元邊長為0.002,整個模型共劃分159 885個單元,32 218個節點。建立的有限元模型如圖3所示。

芯軸受到的載荷有:驅動力矩T1、阻力矩T2和彈簧力F。加載時將力F轉化為面上的壓力,施加在芯軸直徑為D6至D7的環形端面上。驅動力矩T1的加載是利用“Create Component”命令對輸入軸段D7圓周面上的所有節點創建節點組員,然后在節點組元上施加繞Z軸旋轉的轉矩MZ=T1。同理可在輸出軸段D4上施加阻力矩T2,在此要保證驅動力矩與阻力矩方向相反,即所加阻力矩為MZ=-T2。
本模擬中只允許芯軸繞回轉軸線Z軸旋轉,因此對芯軸兩個端面施加除ROTZ外的所有自由度。
當施加完載荷與約束后,即可對芯軸進行靜力學仿真求解,得到芯軸的Von Mises等效應力(如圖4所示)。由圖4可以看出芯軸的最大應力發生在輸入軸端與彈簧的接觸區域。

將節點的等效應力取絕對值,然后按升序排列,使用*get命令提取應力最大值并賦值給Smax;使用*get命令建立極限狀態方程:G=S-Smax。當G>0時,芯軸處于安全狀態;G≤0時,芯軸失效。
整理命令流,生成概率分析文件(.txt文件),進入ANSYS/PDS分析模塊進行可靠性分析。具體步驟為:(1)指定分析文件;(2)設定輸入變量和輸出變量;(3)選擇分析方法;(4)模擬求解。
本模擬中設定的隨機輸入變量如表1所示,輸出變量為G和Smax。

表1 隨機變量統計表Tab.1 The statistical table of random variables
ANSYS/PDS概率設計模塊提供了兩種可靠性分析方法:蒙特卡羅(Monte-carlo)模擬法和響應面法[2]。蒙特卡羅模擬技術是概率分析中最常用的方法,它清楚地模擬實際問題的真實行為特征。一個仿真循環就代表一個加工制造的零件,該零件承受一個特定系列的載荷和邊界條件的作用[3-4]。模擬次數越多,說明抽樣數目越大,精度就越高[5]。在ANSYS中,蒙特卡羅模擬技術可以選擇直接抽樣或拉丁超立方抽樣(簡稱LHS抽樣法)進行分析。LHS抽樣法具有抽樣“記憶”功能,可以避免直接抽樣法數據點集中而導致的仿真循環重復問題[6]。因此,本文選擇 Monte-carlo中的LHS抽樣法進行仿真分析。
為獲得較高的模擬精度,取模擬樣本次數為1 000次,得最大等效應力的抽樣過程如圖5所示。圖5中,最大等效應力隨著抽樣次數的增加逐漸趨于平穩,說明采用的循環次數足夠。

可靠度是產品在規定的時間內和規定的條件下,完成預定功能的概率,是衡量可靠性水平高低的重要指標之一[7]。由可靠性理論可知,求一個結構的可靠度就是求極限狀態方程大于零的概率。因此,極限狀態方程G>0時的概率,即為芯軸的可靠度R,即R=P(G>0)。
經過計算分析后,可以查看置信度為0.95時,極限狀態方程G>0的概率,計算結果如下:

結果表明芯軸在模擬樣本為1 000,置信度為0.95 時的可靠度為0.999 996,接近于1,說明該設計具有很高的可靠性。
概率靈敏度是一個非常重要的參數,通過該參數可以查看隨機變量對可靠度的影響程度。在設計和制造過程中,應嚴格控制這些參數,得到更可靠、質量更好的產品。
圖6是由“Spearman Rank”表示的芯軸可靠度對隨機輸入變量的靈敏度圖。

條狀圖按靈敏度的大小,從左到右依次排列,0刻度以上的表示可靠度與該參數成正比,0刻度以下的表示可靠度與該參數成反比。餅狀圖從90°軸線開始,按靈敏度的大小,順時針旋轉排列,餅狀圖只能反映隨機變量對可靠度的影響程度,而不能反映可靠度的變化規律。通過ANSYS/PDS輸出文件可查看各隨機輸入參數的靈敏度數值(如圖7所示)。

圖7中,正值表示芯軸可靠度隨該參數的增大而增大,負值表示芯軸可靠度隨該參數的減小而減小,數值大小表示該參數對芯軸可靠度的影響程度。參數S的改變對芯軸可靠度的影響比例為0.993,參數F的改變對芯軸可靠度的影響比例為-0.157,說明隨機變量S對芯軸的可靠度影響程度最大,且隨參數S的增大而增大。
當設計無法滿足可靠度要求時,可根據各參數對芯軸可靠度的影響規律對各參數進行修正,使芯軸的可靠度達到設計的要求。本研究為改進軌枕螺栓電動扳手芯軸的設計,提高芯軸的可靠性提供了參考,具有一定的現實意義。
[1] 高 暉,王延遐,董 敏,等.基于ANSYS活塞頭部的可靠性分析[J] .農業裝備與車輛工程,2011(3):15-17.
[2] 宮恩祥,黃銘科,葉 娟,等.基于ANSYS的液力透平軸可靠性分析[J] .石油機械,2011,39(1):39- 42.
[3] 施衛東,顏品蘭,蔣小平,等.基于ANSYS/PDS的泵軸可靠性分析[J] .排灌機械,2008,26(5):1-4.
[4] 張愛華,任工昌.基于ANSYS的概率設計的高速電主軸抗共振的可靠性分析[J] .機械設計與制造,2010(7):112-114.
[5] 葉勇,郝艷華,張昌漢.基于ANSYS的機構可靠性分析[J] .機械工程與自動化,2004(6):63-65.
[6] 張洪才,何 波.有限元分析 -ANSYS13.0從入門到實戰構[M] .北京:機械工業出版社,2011.
[7] 孫志禮,陳良玉.實用機械可靠性設計理論與方法[M] .北京:科學出版社,2003.