摘 要：數(shù)控仿真系統(tǒng)運(yùn)行過(guò)程中圖形的真實(shí)感和刷新速率是衡量該系統(tǒng)優(yōu)劣的兩個(gè)重要指標(biāo)。作者從工件模型的幾何表示和仿真過(guò)程中使用的局部繪制技術(shù)兩個(gè)方面來(lái)研究和實(shí)現(xiàn)三軸數(shù)控銑削加工仿真系統(tǒng)，在真實(shí)感圖形的效果和刷新速率方面均取得了理想的效果。 關(guān)鍵詞：數(shù)控仿真，真實(shí)感圖形，局部繪制，CAD/CAM 1 引言 　　隨著機(jī)械零件制造復(fù)雜程度的提高，對(duì)數(shù)控代碼正確性驗(yàn)證的需求越來(lái)越迫切，數(shù)控仿真系統(tǒng)利用計(jì)算機(jī)圖形顯示技術(shù)模擬實(shí)際加工過(guò)程，是驗(yàn)證數(shù)控加工程序正確性的有力工具之一[1]。國(guó)外已有比較成熟的數(shù)控加工仿真系統(tǒng)，如美國(guó)的CGTech公司的Vericut系統(tǒng)等[2];國(guó)內(nèi)眾多研究機(jī)構(gòu)也在該方面做了大量的研究工作，推出了比較成熟的數(shù)控仿真系統(tǒng) ，但相比國(guó)外同類產(chǎn)品而言，在仿真過(guò)程動(dòng)畫的效果方面仍然存在一定的差距。因此作者在研究過(guò)程中，借鑒部分已有成果 [3][4]，對(duì)傳統(tǒng)的Z-MAP方法進(jìn)行改進(jìn)，設(shè)計(jì)并實(shí)現(xiàn)了一種基于OpenGL的模板緩存機(jī)制、繪制范圍二次確定的局部繪制算法的三軸數(shù)控銑削加工仿真系統(tǒng)。 2 改進(jìn)的Z-MAP方法 　　假設(shè)刀具平行于z軸，三軸數(shù)控銑削加工過(guò)程存在兩個(gè)明顯的特點(diǎn)[5]：1）只有工件的上表面才是加工表面;2）以工件底平面上任意一點(diǎn)為起點(diǎn)且平行于z軸的一條射線與工件的上表面有且僅有一個(gè)交點(diǎn)。 　　正是在這樣的條件下，將Z-MAP方法應(yīng)用于三軸數(shù)控銑削仿真系統(tǒng)的工件模型的幾何表示才成為可能。 　　2.1 傳統(tǒng)的Z-MAP方法 　　Hsu和Yang[6]提出的Z-MAP方法是一種特殊的基于離散模型的表示方法。該方法的基本思想是將原始的工件模型假定為長(zhǎng)方體，將該長(zhǎng)方體在XOY平面的投影得到的長(zhǎng)方形按照一定精度離散為m×n個(gè)均勻分布的正方形網(wǎng)格，從而可以將整個(gè)工件模型離散為m×n個(gè)以各個(gè)正方形為底面的小長(zhǎng)方體集合。 　　在引入Z-MAP方法建立工件的幾何表示后，刀具對(duì)工件實(shí)際切削過(guò)程的幾何表示就被近似表示為離散小長(zhǎng)方體高度連續(xù)更新（降低）的過(guò)程，因此該方法被形象地稱為“割草法”。 　　2.2 改進(jìn)后的Z-MAP方法 　　作者在研究過(guò)程中采用了Z-MAP方法建立工件的幾何表示，發(fā)現(xiàn)該方法仍然存在著不足?；赯-MAP方法建立的工件模型僅僅是一系列不同高度的小長(zhǎng)方體單元的集合，這樣處理的后果是嚴(yán)重破壞了工件模型上表面的幾何連續(xù)性，特別是當(dāng)離散精度較低時(shí)，難以表達(dá)切削過(guò)程中工件上表面的局部曲面特征。因此作者在傳統(tǒng)的Z-MAP方法的基礎(chǔ)上進(jìn)行了改進(jìn)，設(shè)計(jì)了改進(jìn)的Z-MAP方法，算法思想如下： 　　1）在按照一定的離散精度將工件模型的底平面離散為一系列的正方形網(wǎng)格后，以各網(wǎng)格點(diǎn)為起點(diǎn)、工件模型頂面為終點(diǎn)可以得到一系列的z向線段; 　　2）將所有的z向線段與工件模型上表面的交點(diǎn)按照一定的規(guī)則以三角形的形式連接構(gòu)成的三角形網(wǎng)格面即可被視為工件模型的近似表示。 3 局部繪制技術(shù) 　　數(shù)控仿真過(guò)程中工件模型形狀特征的改變僅在當(dāng)前切削段的局部區(qū)域中，其他場(chǎng)景保持不變，如果僅對(duì)當(dāng)前切削段局部區(qū)域進(jìn)行局部繪制將會(huì)極大地減少參與消隱、光照和圖形生成的數(shù)據(jù)量，并最終提高圖形的刷新速率，獲得每秒至少24幀的動(dòng)畫顯示效果。 　　作者在數(shù)控仿真的研究過(guò)程中，設(shè)計(jì)并實(shí)現(xiàn)了復(fù)雜場(chǎng)景下基于Z-MAP方法的工件模型的局部繪制方法，其基本流程如下： 　　1）在設(shè)備坐標(biāo)系下確定切削段刀具運(yùn)動(dòng)包絡(luò)體在屏幕上的顯示區(qū)域 Q; 　　2）世界坐標(biāo)系下在切削最低平面上初選工件模型的矩形區(qū)域F，F(xiàn)是圖形顯示結(jié)果和Q發(fā)生重疊關(guān)系的最小區(qū)域; 　　3）在切削最低平面上精確排除上述矩形區(qū)域F中不與Q發(fā)生重疊關(guān)系的部分區(qū)域，并最終精確確定局部區(qū)域R; 　　4）使用模板緩存技術(shù)在屏幕顯示區(qū)域Q中顯示由精選得到的局部繪制區(qū)域R。 　　3.1 切削段刀具運(yùn)動(dòng)包絡(luò)體的確定 　　世界坐標(biāo)系中工件模型三角網(wǎng)格面局部繪制范圍的確定過(guò)程，從本質(zhì)上講，是在工件模型空間中尋求可能會(huì)被某切削段的刀具運(yùn)動(dòng)包絡(luò)體所遮擋的最大三角網(wǎng)格范圍;而設(shè)備坐標(biāo)系中圖形局部繪制區(qū)域的確定也是尋求某切削段的刀具運(yùn)動(dòng)包絡(luò)體的最大圖形顯示范圍。精確意義上的刀具運(yùn)動(dòng)包絡(luò)體的確定需要考慮刀具類型的影響，特別是球頭刀包絡(luò)體的組成比較復(fù)雜。但在實(shí)際確定該模型時(shí)，基于以下兩個(gè)原因，可對(duì)模型的建立進(jìn)行了適當(dāng)?shù)暮?jiǎn)化處理：1）由于世界坐標(biāo)系和設(shè)備坐標(biāo)系中所確定的都是最大局部區(qū)域或范圍，因此將任何類型的刀具都簡(jiǎn)化為長(zhǎng)方體模型;2）在實(shí)際切削過(guò)程中，工件外形的變化僅局限在某一切削段局部范圍內(nèi)，而且工件上表面是工件外形在高度方向上變化的最大范圍。 　　根據(jù)上述分析，對(duì)切削段刀具運(yùn)動(dòng)包絡(luò)體的確定本質(zhì)上就是對(duì)刀具運(yùn)動(dòng)包絡(luò)體頂點(diǎn)位置的確定。如圖3-1所示，根據(jù)這八個(gè)頂點(diǎn)的坐標(biāo)依次確定刀具運(yùn)動(dòng)包絡(luò)體的六個(gè)表面就確定了該刀具運(yùn)動(dòng)包絡(luò)體，而刀具運(yùn)動(dòng)包絡(luò)體的確定是局部繪制算法的基礎(chǔ)。 [align=center] 圖3-1 切削段刀具運(yùn)動(dòng)包絡(luò)體[/align] 　　3.2 世界坐標(biāo)系中局部繪制區(qū)域的確定 　　世界坐標(biāo)系中工件模型上表面三角網(wǎng)格面局部繪制區(qū)域確定的目的是查找出所有可能被刀具運(yùn)動(dòng)包絡(luò)體所遮擋的最大范圍。為了體現(xiàn)“最大范圍”的原則，作者在算法設(shè)計(jì)中引入了“切削最低平面”的概念， “切削最低平面”的定義如下： 　　 [align=center] 圖3-2 局部繪制區(qū)域的確定示意圖[/align] 　　其中Z[sub]MIN[/sub]是指所有切削過(guò)程全部結(jié)束后工件模型上表面三角網(wǎng)格面頂點(diǎn)中z坐標(biāo)的最小值，而且工件模型空間中存在這樣一個(gè)明確的規(guī)律：對(duì)于在“切削最低平面”上沒(méi)有被刀具運(yùn)動(dòng)包絡(luò)體所遮擋的三角網(wǎng)格面部分而言，雖然實(shí)際情況下該區(qū)域的三角網(wǎng)格面往往高于“切削最低平面”，但在這種情況下這些區(qū)域也絕對(duì)不可能被刀具運(yùn)動(dòng)包絡(luò)體所遮擋，因此這部分區(qū)域在局部繪制時(shí)可以不予計(jì)算和顯示。如圖3-2中的由四條白色直線所確定的平面為“切削最低平面”，黑色長(zhǎng)方體為刀具運(yùn)動(dòng)包絡(luò)體。 　　為了確定局部繪制區(qū)域，最簡(jiǎn)單的方法是依次判斷平行于X（或Y）方向的各條網(wǎng)格線的假想顯示圖形是否和刀具運(yùn)動(dòng)包絡(luò)體的假想顯示圖形相交，如果相交則精確確定每條網(wǎng)格線上和刀具運(yùn)動(dòng)包絡(luò)體顯示圖形的相貫區(qū)間，所有相貫區(qū)間的集合即為局部繪制區(qū)域。但是這種方法存在的問(wèn)題是計(jì)算量大、計(jì)算效率低。為此作者采用了初選和精選兩個(gè)階段來(lái)提高計(jì)算效率。初選的目的是確定工件上表面和刀具運(yùn)動(dòng)包絡(luò)體發(fā)生遮擋關(guān)系的最小矩形區(qū)域，圖3-2中所示的四條白色粗線所包圍的矩形區(qū)域即為初選結(jié)果;精選的目的是在初選所得矩形區(qū)域內(nèi)排除絕對(duì)不發(fā)生遮擋關(guān)系的范圍，而其余部分即為由于可能發(fā)生遮擋關(guān)系而需要重新繪制的區(qū)域，圖3-2中所示的在刀具運(yùn)動(dòng)包絡(luò)體中的白色區(qū)域部分即為精選結(jié)果。 　　3.3 屏幕圖形局部繪制區(qū)域的確定 　　根據(jù)得到的刀具運(yùn)動(dòng)包絡(luò)體，可以確定在屏幕中所改變的區(qū)域僅在該刀具運(yùn)動(dòng)包絡(luò)體范圍內(nèi)，而在刀具運(yùn)動(dòng)包絡(luò)體以外的區(qū)域?qū)⒉话l(fā)生改變，因此作者在研究和實(shí)現(xiàn)的過(guò)程中使用了OpenGL所提供的模板緩存機(jī)制來(lái)限制屏幕的繪制范圍。 　　OpenGL所提供的模板緩存將繪制限制在屏幕的某些部分進(jìn)行。模板緩存可以保持屏幕上某些部分的圖形不變化，而其他部分仍然可以進(jìn)行正常的圖形顯示[7]。在三軸數(shù)控銑削加工仿真中，局部繪制的區(qū)域只可能在當(dāng)前切削段刀具運(yùn)動(dòng)包絡(luò)體中，因此將由當(dāng)前切削段刀具運(yùn)動(dòng)包絡(luò)體在屏幕上的顯示區(qū)域設(shè)置為模板。 　　3.4 局部繪制的完整流程 　　步驟1：清空模板緩存; 　　步驟2：?jiǎn)⒂蒙疃染彺婧湍０寰彺娈嫷毒哌\(yùn)動(dòng)包絡(luò)體，并且設(shè)置顏色緩存中的信息為不能修改，不顯示所畫的刀具運(yùn)動(dòng)包絡(luò)體，但相應(yīng)的模板緩存中的信息設(shè)置為1; 　　步驟3：開(kāi)啟深度緩存，繪制一個(gè)遠(yuǎn)平面; 　　步驟4：設(shè)置顏色緩存為可以寫入，繪制工件模型改變部分。 　　在數(shù)控加工圖形仿真過(guò)程中，要使局部重繪區(qū)域的圖形有正確的消隱效果，就要將工件模型在屏幕所確定區(qū)域的深度信息設(shè)置為最大值。OpenGL沒(méi)有提供修改局部區(qū)域深度信息的功能，因此作者使用繪制遠(yuǎn)平面的方法修改由工件模型所確定的包絡(luò)體的深度信息。由于該遠(yuǎn)平面只是用于修改深度信息的一個(gè)媒介，因此不需要顯示在屏幕上。 4 結(jié)論 　　作者根據(jù)工件模型的幾何表示和局部繪制技術(shù)實(shí)現(xiàn)了三軸數(shù)控銑削仿真系統(tǒng)，現(xiàn)以半徑為4的球頭刀銑削一個(gè)半圓柱形的零件為例在三軸數(shù)控銑削加工仿真系統(tǒng)進(jìn)行了實(shí)驗(yàn)。當(dāng)工件模型的離散精度為0.8時(shí)，工件模型上表面三角網(wǎng)格頂點(diǎn)個(gè)數(shù)為187×300個(gè)，使用全局繪制方法刷新速率為每秒（28.723±1）幀，而使用局部繪制方法刷新速率為每秒（67.703±1）幀，效果圖如圖4-1所示。當(dāng)工件模型的離散精度為0.5時(shí)，工件模型上表面三角網(wǎng)格頂點(diǎn)個(gè)數(shù)為300×480個(gè)，使用局部繪制方法刷新速率為每秒（10.264±1）幀，使用局部繪制方法刷新速率為每秒（64.170±1）幀，效果圖如圖4-2所示。 [align=center] 圖4-1 離散精度為0.8時(shí)仿真效果圖 圖4-2 離散精度為0.5時(shí)仿真效果圖[/align] 　　本文作者創(chuàng)新點(diǎn)：對(duì)傳統(tǒng)的Z-MAP方法改進(jìn)的基礎(chǔ)上，設(shè)計(jì)并實(shí)現(xiàn)了一種基于OpenGL的模板緩存機(jī)制、繪制范圍二次確定的局部繪制算法的三軸數(shù)控銑削加工仿真系統(tǒng)。本系統(tǒng)使用的數(shù)控代碼來(lái)源于鄭州輕工業(yè)學(xué)院機(jī)電工程學(xué)院數(shù)控加工中心FANUC數(shù)控機(jī)床，Pentium® 4 CPU 1.49GHZ,512MB內(nèi)存的測(cè)試環(huán)境，直接讀取NC文件，得到以上測(cè)試數(shù)據(jù)，該系統(tǒng)的應(yīng)用每年將節(jié)約開(kāi)發(fā)經(jīng)費(fèi)約50萬(wàn)元。 參考文獻(xiàn) 　　[1]. 應(yīng)珂.基于OpenGL的數(shù)控仿真三維圖形的研究[J].微計(jì)算機(jī)信息，2006，22（9-1）:226-228. 　　[2]. 崔洪斌，方憶湘，張嘉鈺等.計(jì)算機(jī)輔助設(shè)計(jì)基礎(chǔ)及應(yīng)用[M]. 北京：清華大學(xué)出版社，2002. 　　[3]. 姜曉峰. 數(shù)控加工仿真關(guān)鍵技術(shù)研究[D]. 南京：南京航空航天大學(xué)，1999. 　　[4]. 伍鐵軍， 周來(lái)水， 周儒榮 . 數(shù)控仿真的實(shí)時(shí)真實(shí)感圖形顯示[J]. 計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)報(bào) ，2000，12（4）： 291—293. 　　[5]. 羅堃. 三角片離散法實(shí)現(xiàn)數(shù)控銑床加工仿真[J]. 計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào),2001，13（11）：1024—1028. 　　[6]. Sang-Kyu Lee, Sung-Lim Ko .Development of simulation system for machining process using enhance Z map model [J]. Journal of Materials Processing Technology，2002，（ 130-131）：608-617. 　　[7]. Mason Woo， Jackie Neider， Tom Davis， Dave Shreiner 著 吳斌 ，段海波，薛鳳武 譯. OpenGL 編程權(quán)威指南（第三版）[M].北京：中國(guó)電力出版社，2003.