【正文】 再次感謝所有對(duì)我支持和幫助的人。 在近 1年的論文設(shè)計(jì)和學(xué)習(xí)中，感謝同我一項(xiàng)目組中共同學(xué)習(xí)的陳井春和徐旭陽 的幫助，也由衷感謝和我一起朝夕相處的 125宿舍的舍友以及其他的師兄師姐們，他們?cè)o我的學(xué)習(xí)和生活提供了許多幫助，在我一度頹廢不振的時(shí)候給予我精神上的支持和學(xué)術(shù)上的一并探討，使我順利走出困境，柳暗花明又一村。從方案的提出以及確定，再到方案相關(guān)設(shè)計(jì)理論 以及后來理論方面的闡述，得到了王老師的悉心指導(dǎo)，沒有她的幫助，這篇論文是沒有辦法如此順利的完成的。在有限元分析中，實(shí)現(xiàn)對(duì)機(jī)床主軸系統(tǒng)的全參數(shù)化建模，將可更加全面的對(duì)其進(jìn)行優(yōu)化設(shè)計(jì)。如設(shè)計(jì)冷卻系統(tǒng)控制軸承的溫升，采用熱誤差補(bǔ)償降低主軸頭的跳動(dòng)量等。 1. 在對(duì)機(jī)床主軸系統(tǒng)進(jìn)行溫度場(chǎng)分析時(shí)，關(guān)鍵問題是部件間的接觸傳導(dǎo)率、零部件表面與流體間的對(duì)流換熱系數(shù)等理論計(jì)算還不具有普遍性 ，對(duì)這些理論的完善，使得機(jī)床主軸系統(tǒng)熱邊界條件的精確計(jì)算成為可能，在以后的設(shè)計(jì)中，將能夠?qū)崿F(xiàn)對(duì)機(jī)床主軸系統(tǒng)更為合理的設(shè)計(jì)。 本論文最大特點(diǎn)是 理論和實(shí)際相結(jié)合，對(duì) CK6140數(shù)控車床主軸的溫度分布有了客觀的描述，并為以后提高車床的精度，效率奠定了堅(jiān)實(shí)的基礎(chǔ)。 2. 基于有限元法建立了溫度場(chǎng)數(shù)學(xué)模型， 對(duì) CK6140數(shù)控車床主軸在高速運(yùn)動(dòng)過程中的溫度有了直觀的了解 。 32 結(jié)論與展望 本 論文 針對(duì)主軸在高速運(yùn)動(dòng)的過程中產(chǎn)生大量的熱量使主軸溫度升高，導(dǎo)致主軸變形，從而影響加工精度這一不利因素，因此有必要對(duì)主軸熱變形進(jìn)行大量研究，提出補(bǔ)償措施，提高數(shù)控機(jī)床的加工精度。 對(duì)主軸軸線呈對(duì)稱分布的溫度場(chǎng)使主軸產(chǎn)生偏斜的可能性較小，反之，即使溫升不是很高，也有可能使主軸產(chǎn)生較大的偏斜。 對(duì)主軸加載過程中要注意系統(tǒng)單位的統(tǒng)一， ANSYS 系統(tǒng)中用的單位為國際單位，要注意將機(jī)械設(shè)計(jì)中通用的毫米（ mm）單位轉(zhuǎn)化為米（ m） 。 如圖 46所示。選擇主菜單 → General Postproc→ Plot Results→ Contour Plot→ Nodal Solu，出現(xiàn) Contour Nodal Solution Data 對(duì)話框，如圖45所示，選擇 Nodal Temperature. 30 圖 45 查看結(jié)果 6) 最終計(jì)算得到主軸溫度場(chǎng)分布云圖。選擇主菜單選擇主菜單 → Solution→ Solve→ Current 此時(shí)通過計(jì)算及分析，對(duì)車床主軸進(jìn)行了溫度場(chǎng)的分析。 輸入相關(guān)參數(shù)。再 重復(fù)以上步驟，依次選擇 B和 C表面的所有節(jié)點(diǎn) 施加同樣的載荷 。再選擇 Utility Menu→ Select→ Entities 命令，出現(xiàn) Select Entities 對(duì)話框，在第一個(gè)列表框中選擇 Nodes 選項(xiàng) ,在第二個(gè)下拉列表中選擇 Attached to 選項(xiàng)，在選擇 Areas All 單選按鈕，單擊 OK。如圖 42 所示。軸承 A 發(fā)熱量 ? ， B 發(fā)熱量 ? ， C 發(fā)熱量 ? 。( m m ) 9 0 . 78 7 0/)1 0 3161 9 5251 2 0121 1 2101 1 5272 0 0 5511 0 51 1 2901 1 5861 1 3801 1 021 7 0753074447035(???????????????????????????????sdsr ad / 21200 ??? ?? 查閱得：空氣普朗特?cái)?shù)： Pr= 空氣導(dǎo)熱系數(shù)： ?? 空氣運(yùn)動(dòng)粘度： smf / 26???? 得： 雷諾數(shù)： 46232 )( ??? ????? ??fsD v d? 則，努謝爾特?cái)?shù)： 0 )107(1 3 3 3/13/243/13/2 ?????? DDNu 所以有： 3???? h 計(jì)算得，對(duì)流換熱系數(shù) h= W/(m2 ?℃ ) 主軸系統(tǒng)溫度場(chǎng)分布 以環(huán)境溫度 40℃ ， CA6140 數(shù)控車床主軸導(dǎo)熱系數(shù) k 為 (m?℃ )。在強(qiáng)迫對(duì)流條件下，當(dāng)主軸以一定的轉(zhuǎn)速旋轉(zhuǎn)時(shí)，與空氣間的對(duì)流換熱系數(shù)可按下式計(jì)算： 3/13/2 DDNu ? （ , 5 ????D ） ? sD dhNu ?? 式中， DNu —— 努謝爾特?cái)?shù)； DRe —— 雷諾數(shù)，fsD v d 2Re ?? ? ； Pr —— 普朗特?cái)?shù)； h —— 對(duì)流換熱系數(shù)； ? —— 空氣導(dǎo)熱系數(shù)； sd —— 主軸當(dāng)量直徑， l ldldldd nns ???? ?2211， nllll ?21 ?? ； ? —— 主軸角速度； fv —— 空氣運(yùn)動(dòng)粘度。 1．導(dǎo)熱系數(shù)：根據(jù)傳熱學(xué)理論， CA6140 主軸材料為 45 鋼，導(dǎo)熱系數(shù)為(m?℃ )。 需要確定的有限元模型的邊界條件和初始條件有：主軸前后體上的熱流率、主軸系統(tǒng)與空氣的對(duì)流換熱、溫度分布初始條件以及機(jī)床主軸部件的熱傳導(dǎo)。 將其作為熱載荷施加于主軸系統(tǒng)有限元模型的軸承體上，就可以進(jìn)一步計(jì)算主軸箱系統(tǒng)的溫度場(chǎng) 。 主軸扭矩計(jì)算公式如下 [12]： nPT /9550? 其中 T 為主軸傳動(dòng)的扭矩，單位 mN? ， P 表示軸的計(jì)算功率，單位 KW， n表示主軸的計(jì)算轉(zhuǎn)速， CA6140 數(shù)控車床以 1200r/min 的正常轉(zhuǎn)速計(jì)算。 其中 Fr， Ft分別表示主軸上大齒輪所受軸向力和切向力。 傳遞全功率的最低轉(zhuǎn)速處于低速傳動(dòng)組內(nèi)，所有動(dòng)力通過小齒輪以 26： 58的比例傳給主軸。 mdPfM 111 ? （ N?mm） 式中， 1f —— 與軸承類型和所受負(fù)荷有關(guān)的經(jīng)驗(yàn)常數(shù)； 1P —— 確定軸承摩擦力矩的計(jì)算負(fù)荷（ N）； md —— 軸承中徑（ mm）。 21 MMM ?? 負(fù)荷項(xiàng) 1M 與上述第 (1)、（ 2）、（ 3）項(xiàng)有關(guān)，它決定了起動(dòng)摩擦力矩和低速運(yùn)轉(zhuǎn)時(shí)的摩擦力矩的大小，反映了彈性滯后和局部差動(dòng)滑動(dòng)的摩擦損耗。這樣，就可以把熱變形的計(jì)算問題，轉(zhuǎn)化為靜剛度的計(jì)算問題。因此，在單元與單元聯(lián)接的節(jié)點(diǎn)上引起了熱應(yīng)力。由于溫度的變化，每個(gè)單元都要產(chǎn)生熱變形。這樣就可以用差分的原理，把描述熱傳導(dǎo)的復(fù)雜的偏微分方程簡(jiǎn)化為若干個(gè)線性代數(shù)方程，并且利用電子計(jì)算機(jī)，解出各個(gè)單元的溫度分布。 熱變形有限元法 熱變形計(jì)算的有限元法的基本思路是，把部件分割為有限（一定）數(shù)量的單元。例如傳動(dòng)軸的兩端如果裝有圓錐滾子軸承，冷態(tài)時(shí)無 間隙 ，工作一段時(shí)間后，由于箱體的散熱條件比較好，軸的溫度將高于 箱體 的溫度； 熱膨脹使軸承內(nèi)產(chǎn)生軸向附加載荷，這個(gè)載荷又將使軸承進(jìn)一步發(fā)熱，嚴(yán)重時(shí)會(huì)損壞軸承。如果前后軸承的溫升相同，則上升后的主軸軸線仍平行于溫升前的軸線；如果前軸承的溫度高于后軸承，則主軸軸線將向上傾斜。不均勻的熱膨脹對(duì)精度的影響比均勻熱膨脹更大。 熱變形和熱應(yīng)力 零件受熱膨脹有兩種可能：一種是均勻的熱膨脹，一種是不均勻的熱膨脹。 機(jī)床的熱量主要是從某幾個(gè)熱源發(fā)出的，所以熱源處溫度最高，離熱源處越遠(yuǎn)則溫度越低，這就形成了溫度場(chǎng)。機(jī)床主軸滑動(dòng)軸承的間隙很小，容許的間隙變化量更小，故滑動(dòng)軸承允許的溫升比滾動(dòng)軸承要低一些。這是因?yàn)闇厣龑⒏淖兓瑒?dòng)軸承的間隙，溫升過高將導(dǎo)致間隙消失而“抱軸”。對(duì)于一般機(jī)床的主軸箱來說就是主軸箱軸承處的溫度和溫升。機(jī)床上的熱源主要是電動(dòng)機(jī)、主軸軸承、高速運(yùn)動(dòng)的導(dǎo)軌、液壓油池、液壓缸等。一般機(jī)床約需運(yùn)轉(zhuǎn)二至三小時(shí)，才能達(dá)到熱平衡。為方便起見，如果每小時(shí)的溫升不超過 5℃ ，就認(rèn)為達(dá)到了熱平衡。 如果熱平衡溫度為 0t ℃，在升溫過程中，經(jīng)過時(shí)間 ? （小時(shí)）與達(dá)到的溫度t℃之間的關(guān)系為 [10]： ???????? ?? ? tkett ?10（℃） kAcmkt ? （ s） 式中 tk —— 時(shí)間常數(shù)（ s）； k —— 向周圍環(huán)境的散熱系數(shù)（ W/ ?2m ℃）； A —— 散熱面積（ 2m ）； c —— 比熱容（ J/kg.℃ ）； m —— 發(fā)熱物質(zhì)的質(zhì)量（ kg）。隨著機(jī)床溫度的升高，溫差加大，散熱也增加。完成網(wǎng)格的劃分如圖 315所示： 圖 315 網(wǎng)格劃分結(jié)果 21 第四章 主軸的溫度場(chǎng)分析 機(jī)床熱變形理論 熱平衡和溫度場(chǎng) 機(jī)床一方面在工作時(shí)產(chǎn)生熱量，另一方面又向周圍發(fā)散熱量。 圖 314 CK6140主軸 實(shí)體模型 劃分網(wǎng)格 由于主軸的各節(jié)之間的尺寸過渡問題，選用自由網(wǎng)格方式對(duì)軸體劃分網(wǎng)格。 選擇 Main Menu: Preprocessor→ Modeling→ Operate→ Booleans→ Subtract→ Volumes，選擇圖 314 建立的幾何體，使其成為一個(gè)整體。 圖 313 回到初始界面圖 7）重復(fù)步驟 2)、 3)、 4)，根據(jù)數(shù)據(jù)完成 主軸內(nèi)孔的繪制 。 圖 312 主軸的外輪廓 19 6）點(diǎn)擊應(yīng)用菜單 WorkPlane→ offset WP to→ Global Origin。然后根據(jù)圖形要求重復(fù)步驟 4。完成圖形如圖 311 右圖所示。在WP X 與 WP Y 分別輸入 0、 0，在 Radius 中輸入 35，在 Depth 中輸入 28。 4） 從主菜單中選擇 Main Menu: Preprocessor→ Modeling→ Create→Volumes→ Cylinder→ Solid cylinder。 圖 310 軸端面視圖 3） 點(diǎn)擊應(yīng)用菜單 WorkPlane→ Offset WP by increments，彈出 Offset WP對(duì)話框，在 Snaps X,Y,Z offsets 輸入 0， 0， 3 點(diǎn)擊“ OK”按鈕。打開界面如圖 310所示。 建立主軸三維實(shí)體模型 1） 點(diǎn)擊右側(cè)工具欄 圖標(biāo)，打開等軸測(cè)視界面。如圖 39所示。 設(shè)置步驟： 1）從主菜單中選擇 Main Menu :Preprocessor → Material Props→ Material Model 命令，將打開 Define Material Model Behavior(定義材料16 模型屬性 )窗口，如圖 38所示。在對(duì)機(jī)床主軸系統(tǒng)進(jìn)行熱特性分析中，進(jìn)行溫度場(chǎng)分析時(shí)必須確定下列參數(shù)導(dǎo)熱系數(shù)、對(duì)流系數(shù)、密度、比熱；熱變形分析時(shí)則必須確定彈性