喜歡這套資料就充值下載吧。。。資源目錄里展示的都可在線預(yù)覽哦。。。下載后都有，，請(qǐng)放心下載，，文件全都包含在內(nèi)，，【有疑問咨詢QQ:1304139763或 1064457796】
============================================喜歡這套資料就充值下載吧。。。資源目錄里展示的都可在線預(yù)覽哦。。。下載后都有，，請(qǐng)放心下載，，文件全都包含在內(nèi)，，【有疑問咨詢QQ:1304139763或 1064457796】
============================================
外文翻譯資料
機(jī)電一體化技術(shù)及其應(yīng)用研究
1 機(jī)電一體化技術(shù)發(fā)展
機(jī)電一體化是機(jī)械、微、控制、機(jī)、信息處理等多學(xué)科的交叉融合，其發(fā)展和進(jìn)步有賴于相關(guān)技術(shù)的進(jìn)步與發(fā)展，其主要發(fā)展方向有數(shù)字化、智能化、模塊化、化、人性化、微型化、集成化、帶源化和綠色化。
1.1 數(shù)字化
微控制器及其發(fā)展奠定了機(jī)電產(chǎn)品數(shù)字化的基礎(chǔ)，如不斷發(fā)展的數(shù)控機(jī)床和機(jī)器人；而計(jì)算機(jī)網(wǎng)絡(luò)的迅速崛起，為數(shù)字化設(shè)計(jì)與制造鋪平了道路，如虛擬設(shè)計(jì)、計(jì)算機(jī)集成制造等。數(shù)字化要求機(jī)電一體化產(chǎn)品的軟件具有高可靠性、易操作性、可維護(hù)性、自診斷能力以及友好人機(jī)界面。數(shù)字化的實(shí)現(xiàn)將便于遠(yuǎn)程操作、診斷和修復(fù)。
1.2 智能化
即要求機(jī)電產(chǎn)品有一定的智能,使它具有類似人的邏輯思考、判斷推理、自主決策等能力。例如在CNC數(shù)控機(jī)床上增加人機(jī)對(duì)話功能，設(shè)置智能I/O接口和智能工藝數(shù)據(jù)庫(kù)，會(huì)給使用、操作和維護(hù)帶來極大的方便。隨著模糊控制、神經(jīng)網(wǎng)絡(luò)、灰色、小波理論、混沌與分岔等人工智能技術(shù)的進(jìn)步與發(fā)展，為機(jī)電一體化技術(shù)發(fā)展開辟了廣闊天地。
1.3 模塊化
由于機(jī)電一體化產(chǎn)品種類和生產(chǎn)廠家繁多，研制和開發(fā)具有標(biāo)準(zhǔn)機(jī)械接口、動(dòng)力接口、環(huán)境接口的機(jī)電一體化產(chǎn)品單元模塊是一項(xiàng)復(fù)雜而有前途的工作。如研制具有集減速、變頻調(diào)速電機(jī)一體的動(dòng)力驅(qū)動(dòng)單元；具有視覺、圖像處理、識(shí)別和測(cè)距等功能的電機(jī)一體控制單元等。這樣，在產(chǎn)品開發(fā)設(shè)計(jì)時(shí)，可以利用這些標(biāo)準(zhǔn)模塊化單元迅速開發(fā)出新的產(chǎn)品。
1.4 網(wǎng)絡(luò)化
由于網(wǎng)絡(luò)的普及，基于網(wǎng)絡(luò)的各種遠(yuǎn)程控制和監(jiān)視技術(shù)方興未艾。而遠(yuǎn)程控制的終端設(shè)備本身就是機(jī)電一體化產(chǎn)品，現(xiàn)場(chǎng)總線和局域網(wǎng)技術(shù)使家用電器網(wǎng)絡(luò)化成為可能，利用家庭網(wǎng)絡(luò)把各種家用電器連接成以計(jì)算機(jī)為中心的計(jì)算機(jī)集成家用電器系統(tǒng)，使人們?cè)诩依锟沙浞窒硎芨鞣N高技術(shù)帶來的好處，因此，機(jī)電一體化產(chǎn)品無疑應(yīng)朝網(wǎng)絡(luò)化方向發(fā)展。
1.5 人性化
機(jī)電一體化產(chǎn)品的最終使用對(duì)象是人，如何給機(jī)電一體化產(chǎn)品賦予人的智能、情感和人性顯得愈來愈重要，機(jī)電一體化產(chǎn)品除了完善的性能外，還要求在色彩、造型等方面與環(huán)境相協(xié)調(diào)，使用這些產(chǎn)品，對(duì)人來說還是一種享受，如家用機(jī)器人的最高境界就是人機(jī)一體化。
1.6 微型化
微型化是精細(xì)加工技術(shù)發(fā)展的必然，也是提高效率的需要。微機(jī)電系統(tǒng)(Micro Electronic Mechanical Systems，簡(jiǎn)稱MEMS)是指可批量制作的，集微型機(jī)構(gòu)、微型傳感器、微型執(zhí)行器以及信號(hào)處理和控制電路，直至接口、通信和電源等于一體的微型器件或系統(tǒng)。自1986年美國(guó)斯坦福大學(xué)研制出第一個(gè)醫(yī)用微探針，1988年美國(guó)加州大學(xué)Berkeley分校研制出第一個(gè)微電機(jī)以來，國(guó)內(nèi)外在MEMS工藝、材料以及微觀機(jī)理方面取得了很大進(jìn)展，開發(fā)出各種MEMS器件和系統(tǒng)，如各種微型傳感器（壓力傳感器、微加速度計(jì)、微觸覺傳感器），各種微構(gòu)件（微膜、微粱、微探針、微連桿、微齒輪、微軸承、微泵、微彈簧以及微機(jī)器人等）。
1.7 集成化
集成化既包含各種技術(shù)的相互滲透、相互融合和各種產(chǎn)品不同結(jié)構(gòu)的優(yōu)化與復(fù)合，又包含在生產(chǎn)過程中同時(shí)處理加工、裝配、檢測(cè)、管理等多種工序。為了實(shí)現(xiàn)多品種、小批量生產(chǎn)的自動(dòng)化與高效率，應(yīng)使系統(tǒng)具有更廣泛的柔性。首先可將系統(tǒng)分解為若干層次，使系統(tǒng)功能分散，并使各部分協(xié)調(diào)而又安全地運(yùn)轉(zhuǎn)，然后再通過軟、硬件將各個(gè)層次有機(jī)地聯(lián)系起來，使其性能最優(yōu)、功能最強(qiáng)。
1.8 帶源化
是指機(jī)電一體化產(chǎn)品自身帶有能源，如太陽能電池、燃料電池和大容量電池。由于在許多場(chǎng)合無法使用電能，因而對(duì)于運(yùn)動(dòng)的機(jī)電一體化產(chǎn)品，自帶動(dòng)力源具有獨(dú)特的好處。帶源化是機(jī)電一體化產(chǎn)品的發(fā)展方向之一。
1.9 綠色化
技術(shù)的發(fā)展給人們的生活帶來巨大變化，在物質(zhì)豐富的同時(shí)也帶來資源減少、生態(tài)環(huán)境惡化的后果。所以，人們呼喚保護(hù)環(huán)境，回歸，實(shí)現(xiàn)可持續(xù)發(fā)展，綠色產(chǎn)品概念在這種呼聲中應(yīng)運(yùn)而生。綠色產(chǎn)品是指低能耗、低材耗、低污染、舒適、協(xié)調(diào)而可再生利用的產(chǎn)品。在其設(shè)計(jì)、制造、使用和銷毀時(shí)應(yīng)符合環(huán)保和人類健康的要求，機(jī)電一體化產(chǎn)品的綠色化主要是指在其使用時(shí)不污染生態(tài)環(huán)境，產(chǎn)品壽命結(jié)束時(shí)，產(chǎn)品可分解和再生利用。
2 機(jī)電一體化技術(shù)在鋼鐵中應(yīng)用
在鋼鐵企業(yè)中，機(jī)電一體化系統(tǒng)是以微處理機(jī)為核心，把微機(jī)、工控機(jī)、數(shù)據(jù)通訊、顯示裝置、儀表等技術(shù)有機(jī)的結(jié)合起來，采用組裝合并方式，為實(shí)現(xiàn)工程大系統(tǒng)的綜合一體化創(chuàng)造有力條件，增強(qiáng)系統(tǒng)控制精度、質(zhì)量和可靠性。機(jī)電一體化技術(shù)在鋼鐵企業(yè)中主要應(yīng)用于以下幾個(gè)方面：
2.1 智能化控制技術(shù)(IC)
由于鋼鐵具有大型化、高速化和連續(xù)化的特點(diǎn)，傳統(tǒng)的控制技術(shù)遇到了難以克服的困難，因此非常有必要采用智能控制技術(shù)。智能控制技術(shù)主要包括專家系統(tǒng)、模糊控制和神經(jīng)等，智能控制技術(shù)廣泛于鋼鐵的產(chǎn)品設(shè)計(jì)、生產(chǎn)、控制、設(shè)備與產(chǎn)品質(zhì)量診斷等各個(gè)方面，如高爐控制系統(tǒng)、電爐和連鑄車間、軋鋼系統(tǒng)、煉鋼———連鑄———軋鋼綜合調(diào)度系統(tǒng)、冷連軋等。
2.2 分布式控制系統(tǒng)(ＤCS)
分布式控制系統(tǒng)采用一臺(tái)中央機(jī)指揮若干臺(tái)面向控制的現(xiàn)場(chǎng)測(cè)控計(jì)算機(jī)和智能控制單元。分布式控制系統(tǒng)可以是兩級(jí)的、三級(jí)的或更多級(jí)的。利用計(jì)算機(jī)對(duì)生產(chǎn)過程進(jìn)行集中監(jiān)視、操作、管理和分散控制。隨著測(cè)控技術(shù)的，分布式控制系統(tǒng)的功能越來越多。不僅可以實(shí)現(xiàn)生產(chǎn)過程控制，而且還可以實(shí)現(xiàn)在線最優(yōu)化、生產(chǎn)過程實(shí)時(shí)調(diào)度、生產(chǎn)計(jì)劃統(tǒng)計(jì)管理功能，成為一種測(cè)、控、管一體化的綜合系統(tǒng)。ＤCS具有特點(diǎn)控制功能多樣化、操作簡(jiǎn)便、系統(tǒng)可以擴(kuò)展、維護(hù)方便、可靠性高等特點(diǎn)。ＤCS是監(jiān)視集中控制分散，故障面小，而且系統(tǒng)具有連鎖保護(hù)功能，采用了系統(tǒng)故障人工手動(dòng)控制操作措施，使系統(tǒng)可靠性高。分布式控制系統(tǒng)與集中型控制系統(tǒng)相比，其功能更強(qiáng)，具有更高的安全性。是當(dāng)前大型機(jī)電一體化系統(tǒng)的主要潮流。
2.3 開放式控制系統(tǒng)(OCS)
開放控制系統(tǒng)(Open Control System)是計(jì)算機(jī)技術(shù)發(fā)展所引出的新的結(jié)構(gòu)體系概念?！伴_放”意味著對(duì)一種標(biāo)準(zhǔn)的信息交換規(guī)程的共識(shí)和支持，按此標(biāo)準(zhǔn)設(shè)計(jì)的系統(tǒng)，可以實(shí)現(xiàn)不同廠家產(chǎn)品的兼容和互換，且資源共享。開放控制系統(tǒng)通過工業(yè)通信網(wǎng)絡(luò)使各種控制設(shè)備、管理計(jì)算機(jī)互聯(lián)，實(shí)現(xiàn)控制與經(jīng)營(yíng)、管理、決策的集成，通過現(xiàn)場(chǎng)總線使現(xiàn)場(chǎng)儀表與控制室的控制設(shè)備互聯(lián)，實(shí)現(xiàn)測(cè)量與控制一體化。
2.4 計(jì)算機(jī)集成制造系統(tǒng)(CIMS)
鋼鐵企業(yè)的CIMS是將人與生產(chǎn)經(jīng)營(yíng)、生產(chǎn)管理以及過程控制連成一體，用以實(shí)現(xiàn)從原料進(jìn)廠，生產(chǎn)加工到產(chǎn)品發(fā)貨的整個(gè)生產(chǎn)過程全局和過程一體化控制。目前鋼鐵企業(yè)已基本實(shí)現(xiàn)了過程自動(dòng)化，但這種“自動(dòng)化孤島”式的單機(jī)自動(dòng)化缺乏信息資源的共享和生產(chǎn)過程的統(tǒng)一管理，難以適應(yīng)鋼鐵生產(chǎn)的要求。未來鋼鐵企業(yè)競(jìng)爭(zhēng)的焦點(diǎn)是多品種、小批量生產(chǎn)，質(zhì)優(yōu)價(jià)廉，及時(shí)交貨。為了提高生產(chǎn)率、節(jié)能降耗、減少人員及現(xiàn)有庫(kù)存，加速資金周轉(zhuǎn)，實(shí)現(xiàn)生產(chǎn)、經(jīng)營(yíng)、管理整體優(yōu)化，關(guān)鍵就是加強(qiáng)管理，獲取必須的效益，提高了企業(yè)的競(jìng)爭(zhēng)力。美國(guó)、日本等一些大型鋼鐵企業(yè)在20世紀(jì)80年代已廣泛實(shí)現(xiàn)CIMS化。
2.5 現(xiàn)場(chǎng)總線技術(shù)(ＦBT)
現(xiàn)場(chǎng)總線技術(shù)(Ｆieｄ Bus Technology)是連接設(shè)置在現(xiàn)場(chǎng)的儀表與設(shè)置在控制室內(nèi)的控制設(shè)備之間的數(shù)字式、雙向、多站通信鏈路。采用現(xiàn)場(chǎng)總線技術(shù)取代現(xiàn)行的信號(hào)傳輸技術(shù)(如4～20mA，ＤC直流傳輸)就能使更多的信息在智能化現(xiàn)場(chǎng)儀表裝置與更高一級(jí)的控制系統(tǒng)之間在共同的通信媒體上進(jìn)行雙向傳送。通過現(xiàn)場(chǎng)總線連接可省去66%或更多的現(xiàn)場(chǎng)信號(hào)連接導(dǎo)線?，F(xiàn)場(chǎng)總線的引入導(dǎo)致ＤCS的變革和新一代圍繞開放自動(dòng)化系統(tǒng)的現(xiàn)場(chǎng)總線化儀表，如智能變送器、智能執(zhí)行器、現(xiàn)場(chǎng)總線化檢測(cè)儀表、現(xiàn)場(chǎng)總線化PLC(Programmable Logic Controller)和現(xiàn)場(chǎng)就地控制站等的發(fā)展。
2.6 交流傳動(dòng)技術(shù)
傳動(dòng)技術(shù)在鋼鐵工業(yè)中起作至關(guān)重要的作用。隨著電力技術(shù)和微電子技術(shù)的發(fā)展，交流調(diào)速技術(shù)的發(fā)展非常迅速。由于交流傳動(dòng)的優(yōu)越性，電氣傳動(dòng)技術(shù)在不久的將來由交流傳動(dòng)全面取代直流傳動(dòng)，數(shù)字技術(shù)的發(fā)展，使復(fù)雜的矢量控制技術(shù)實(shí)用化得以實(shí)現(xiàn)，交流調(diào)速系統(tǒng)的調(diào)速性能已達(dá)到和超過直流調(diào)速水平。現(xiàn)在無論大容量電機(jī)或中小容量電機(jī)都可以使用同步電機(jī)或異步電機(jī)實(shí)現(xiàn)可逆平滑調(diào)速。交流傳動(dòng)系統(tǒng)在軋鋼生產(chǎn)中一出現(xiàn)就受到用戶的歡迎，應(yīng)用不斷擴(kuò)大。
4
夾具夾緊力的優(yōu)化及對(duì)工件定位精度的影響
B.Li 和 S.N.Mellkote
布什伍德拉夫機(jī)械工程學(xué)院，佐治亞理工學(xué)院，格魯吉亞，美國(guó)研究所
由于夾緊和加工，在工件和夾具的接觸部位會(huì)產(chǎn)生局部彈性變形，使工件尺寸發(fā)生變化，進(jìn)而影響工件的最終加工質(zhì)量。這種效應(yīng)可通過最小化夾具設(shè)計(jì)優(yōu)化，夾緊力是一個(gè)重要的設(shè)計(jì)變量，可以得到優(yōu)化，以減少工件的位移。本文提出了一種確定多夾緊夾具受到準(zhǔn)靜態(tài)加工部位的最佳夾緊力的新方法。該方法采用彈性接觸力學(xué)模型代表夾具與工件接觸，并涉及制定和解決方案的多目標(biāo)優(yōu)化模型的約束。夾緊力的最優(yōu)化對(duì)工件定位精度的影響通過3-2-1式銑夾具的例子進(jìn)行了分析。
關(guān)鍵詞：彈性 接觸 模型 夾具 夾緊力 優(yōu)化
前言
定位和夾緊的工件加工中的兩個(gè)關(guān)鍵因素。要實(shí)現(xiàn)夾具的這些功能，需將工件定位到一個(gè)合適的基準(zhǔn)上并夾緊，采用的夾緊力必須足夠大，以抑制工件在加工過程中產(chǎn)生的移動(dòng)。然而，過度的夾緊力可誘導(dǎo)工件產(chǎn)生更大的彈性變形 ，這會(huì)影響它的位置精度，并反過來影響零件質(zhì)量。所以有必要確定最佳夾緊力，來減小由于彈性變形對(duì)工件的定位誤差，同時(shí)滿足加工的要求。在夾具分析和綜合領(lǐng)域上的研究人員使用了有限元模型的方法或剛體模型的方法。大量的工作都以有限元方法為基礎(chǔ)被報(bào)道[參考文獻(xiàn)1-8]。隨著得墨忒耳[8]，這種方法的限制是需要較大的模型和計(jì)算成本。同時(shí)，多數(shù)的有限元基礎(chǔ)研究人員一直重點(diǎn)關(guān)注的夾具布局優(yōu)化和夾緊力的優(yōu)化還沒有得到充分討論，也有少數(shù)的研究人員通過對(duì)剛性模型[9-11]對(duì)夾緊力進(jìn)行了優(yōu)化，剛型模型幾乎被近似為一個(gè)規(guī)則完整的形狀。得墨忒耳[12，13]用螺釘理論解決的最低夾緊力，總的問題是制定一個(gè)線性規(guī)劃，其目的是盡量減少在每個(gè)定位點(diǎn)調(diào)整夾緊力強(qiáng)度的法線接觸力。接觸摩擦力的影響被忽視，因?yàn)樗^法線接觸力相對(duì)較小，由于這種方法是基于剛體假設(shè)，獨(dú)特的三維夾具可以處理超過6個(gè)自由度的裝夾，復(fù)和倪[14]也提出迭代搜索方法，通過假設(shè)已知摩擦力的方向來推導(dǎo)計(jì)算最小夾緊力，該剛體分析的主要限制因素是當(dāng)出現(xiàn)六個(gè)以上的接觸力是使其靜力不確定，因此，這種方法無法確定工件移位的唯一性。
這種限制可以通過計(jì)算夾具——工件系統(tǒng)[15]的彈性來克服，對(duì)于一個(gè)相對(duì)嚴(yán)格的工件，該夾具在機(jī)械加工工件的位置會(huì)受夾具點(diǎn)的局部彈性變形的強(qiáng)烈影響。Hockenberger和得墨忒耳[16]使用經(jīng)驗(yàn)的接觸力變形的關(guān)系（稱為元功能），解決由于夾緊和準(zhǔn)靜態(tài)加工力工件剛體位移。同一作者還考察了加工工件夾具位移對(duì)設(shè)計(jì)參數(shù)的影響[17]。桂 [18] 等 通過工件的夾緊力的優(yōu)化定位精度彈性接觸模型對(duì)報(bào)告做了改善，然而，他們沒有處理計(jì)算夾具與工件的接觸剛度的方法，此外，其算法的應(yīng)用沒有討論機(jī)械加工刀具路徑負(fù)載有限序列。李和Melkote [19]和烏爾塔多和Melkote [20]用接觸力學(xué)解決由于在加載夾具夾緊點(diǎn)彈性變形產(chǎn)生的接觸力和工件的位移，他們還使用此方法制定了優(yōu)化方法夾具布局[21]和夾緊力[22]。但是，關(guān)于multiclamp系統(tǒng)及其對(duì)工件精度影響的夾緊力的優(yōu)化并沒有在這些文件中提到 。
本文提出了一種新的算法，確定了multiclamp夾具工件系統(tǒng)受到準(zhǔn)靜態(tài)加載的最佳夾緊力為基礎(chǔ)的彈性方法。該法旨在盡量減少影響由于工件夾緊位移和加工荷載通過系統(tǒng)優(yōu)化夾緊力的一部分定位精度。接觸力學(xué)模型，用于確定接觸力和位移，然后再用做夾緊力優(yōu)化，這個(gè)問題被作為多目標(biāo)約束優(yōu)化問題提出和解決。通過兩個(gè)例子分析工件夾緊力的優(yōu)化對(duì)定位精度的影響，例子涉及的銑削夾具3-2-1布局。
1． 夾具——工件聯(lián)系模型
1．1 模型假設(shè)
該加工夾具由L定位器和帶有球形端的c形夾組成。工件和夾具接觸的地方是線性的彈性接觸，其他地方完全剛性。工件——夾具系統(tǒng)由于夾緊和加工受到準(zhǔn)靜態(tài)負(fù)載。夾緊力可假定為在加工過程中保持不變，這個(gè)假設(shè)是有效的，在對(duì)液壓或氣動(dòng)夾具使用。在實(shí)際中，夾具工件接觸區(qū)域是彈性分布，然而，這種模式的發(fā)展，假設(shè)總觸剛度（見圖1）第i夾具接觸力局部變形如下：
(1) 其中（j=x，y，z）表示，在當(dāng)?shù)刈幼鴺?biāo)系切線和法線方向的接觸剛度
第 19 頁 共 15 頁
圖1 彈簧夾具——
工件接觸模型。
表示在第i個(gè)
接觸處的坐標(biāo)系
（j=x，y，z）是對(duì)應(yīng)沿著xyz方向的彈性變形，分別 （j= x，y，z）的代表和切向力接觸 ，法線力接觸。
1．2 工件——夾具的接觸剛度模型
集中遵守一個(gè)球形尖端定位，夾具和工件的接觸并不是線性的，因?yàn)榻佑|半徑與隨法線力呈非線性變化 [23]。由于法線力接觸變形作用于半徑和平面工件表面之間，這可從封閉赫茲的辦法解決縮進(jìn)一個(gè)球體彈性半空間的問題。對(duì)于這個(gè)問題， 是法線的變形，在[文獻(xiàn)23 第93頁]中給出如下：
（2）
其中式中 和是工件和夾具的彈性模量，、分別是工件和材料的泊松比。
切向變形沿著和切線方向）硅業(yè)切力距有以下形式[文獻(xiàn)23第217頁]
（3）
其中、 分別是工件和夾具剪切模量
一個(gè)合理的接觸剛度的線性可以近似從最小二乘獲得適合式 （2），這就產(chǎn)生了以下線性化接觸剛度值：在計(jì)算上述的線性近似，
（4）
(5)
正常的力被假定為從0到1000N，且最小二乘擬合相應(yīng)的R2值認(rèn)定是0.94。
2．夾緊力優(yōu)化
我們的目標(biāo)是確定最優(yōu)夾緊力，將盡量減少由于工件剛體運(yùn)動(dòng)過程中，局部的夾緊和加工負(fù)荷引起的彈性變形，同時(shí)保持在準(zhǔn)靜態(tài)加工過程中夾具——工件系統(tǒng)平衡，工件的位移減少，從而減少定位誤差。實(shí)現(xiàn)這個(gè)目標(biāo)是通過制定一個(gè)多目標(biāo)約束優(yōu)化問題的問題，如下描述。
2.1 目標(biāo)函數(shù)配方
工件旋轉(zhuǎn)，由于部隊(duì)輪換往往是相當(dāng)小[17]的工件定位誤差假設(shè)為確定其剛體翻譯基本上，其中 、、和 是 沿，和三個(gè)正交組件（見圖2）。
圖2 工件剛體平移和旋轉(zhuǎn)
工件的定位誤差歸于裝夾力，然后可以在該剛體位移的范數(shù)計(jì)算如下：
（6）
其中表示一個(gè)向量二級(jí)標(biāo)準(zhǔn)。
但是作用在工件的夾緊力會(huì)影響定位誤差。當(dāng)多個(gè)夾緊力作用于工件，由此產(chǎn)生的夾緊力為,有如下形式：
（7）
其中夾緊力是矢量，夾緊力的方向矩陣，是夾緊力是矢量的方向余弦，、和 是第i個(gè)夾緊點(diǎn)夾緊力在、和方向上的向量角度（i=1、2、3...,C）。
在這個(gè)文件中，由于接觸區(qū)變形造成的工件的定位誤差，被假定為受的作用力是法線的，接觸的摩擦力相對(duì)較小，并在進(jìn)行分析時(shí)忽略了加緊力對(duì)工件的定位誤差的影響。意指正常接觸剛度比，是通過（i=1，2…L）和最小的所有定位器正常剛度相乘，并假設(shè)工件、、取決于、、的方向，各自的等效接觸剛度可有下式計(jì)算得出（見圖3），工件剛體運(yùn)動(dòng)，歸于夾緊行動(dòng)現(xiàn)在可以寫成：
(8)
工件有位移，因此，定位誤差的減小可以通過盡量減少產(chǎn)生的夾緊力向量 范數(shù)。因此，第一個(gè)目標(biāo)函數(shù)可以寫為：
最小化 （9）
要注意，加權(quán)因素是與等效接觸剛度成正比的在、和 方向上。通過使用最低總能量互補(bǔ)參考文獻(xiàn)[15，23]的原則求解彈性力學(xué)接觸問題得出A的組成部分是唯一確定的，這保證了夾緊力和相應(yīng)的定位反應(yīng)是“真正的”解決方案，對(duì)接觸問題和產(chǎn)生的“真正”剛體位移，而且工件保持在靜態(tài)平衡，通過夾緊力的隨時(shí)調(diào)整。因此，總能量最小化的形式為補(bǔ)充的夾緊力優(yōu)化的第二個(gè)目標(biāo)函數(shù)，并給出：
最小化 （10）
其中代表機(jī)構(gòu)的彈性變形應(yīng)變能互補(bǔ)，代表由外部力量和力矩配合完成，是遵守對(duì)角矩陣的， 和是所有接觸力的載體。
如圖3 加權(quán)系數(shù)計(jì)算確定的基礎(chǔ)
內(nèi)蒙古科技大學(xué)本科生畢業(yè)設(shè)計(jì)（外文翻譯）
2.2 摩擦和靜態(tài)平衡約束
在（10）式優(yōu)化的目標(biāo)受到一定的限制和約束，他們中最重要的是在每個(gè)接觸處的靜摩擦力約束。庫(kù)侖摩擦力的法律規(guī)定（是靜態(tài)摩擦系數(shù)）,這方面的一個(gè)非線性約束和線性化版本可以使用，并且[19]有：
（11）
假設(shè)準(zhǔn)靜態(tài)載荷，工件的靜力平衡由下列力和力矩平衡方程確保（向量形式）：
(12)
其中包括在法線和切線方向的力和力矩的機(jī)械加工力和工件重量。
2.3界接觸力
由于夾具——工件接觸是單側(cè)面的，法線的接觸力只能被壓縮。這通過以下的的約束表（i=1，2…,L+C） （13）
它假設(shè)在工件上的法線力是確定的，此外，在一個(gè)法線的接觸壓力不能超過壓工件材料的屈服強(qiáng)度（）。這個(gè)約束可寫為：
（i=1，2，…,L+C） (14)
如果是在第i個(gè)工件——夾具的接觸處的接觸面積，完整的夾緊力優(yōu)化模型，可以寫成：最小化 (15)
3．模型算法求解
式（15）多目標(biāo)優(yōu)化問題可以通過求解約束[24]。這種方法將確定的目標(biāo)作為首要職能之一，并將其轉(zhuǎn)換成一個(gè)約束對(duì)。該補(bǔ)充（）的主要目的是處理功能，并由此得到夾緊力（）作為約束的加權(quán)范數(shù)最小化。對(duì)為主要目標(biāo)的選擇，確保選中一套獨(dú)特可行的夾緊力，因此，工件——夾具系統(tǒng)驅(qū)動(dòng)到一個(gè)穩(wěn)定的狀態(tài)（即最低能量狀態(tài)），此狀態(tài)也表示有最小的夾緊力下的加權(quán)范數(shù)。 的約束轉(zhuǎn)換涉及到一個(gè)指定的加權(quán)范數(shù)小于或等于，其中是 的約束，假設(shè)最初所有夾緊力不明確，要確定一個(gè)合適的。在定位和夾緊點(diǎn)的接觸力的計(jì)算只考慮第一個(gè)目標(biāo)函數(shù)（即）。雖然有這樣的接觸力，并不一定產(chǎn)生最低的夾緊力，這是一個(gè)“真正的”可行的解決彈性力學(xué)問題辦法，可完全抑制工件在夾具中的位置。這些夾緊力的加權(quán)系數(shù)，通過計(jì)算并作為初始值與比較，因此，夾緊力式（15）的優(yōu)化問題可改寫為:
最小化 （16）
由： (11)–(14) 得。
類似的算法尋找一個(gè)方程根的二分法來確定最低的上的約束, 通過盡可能降低上限，由此產(chǎn)生的最小夾緊力的加權(quán)范數(shù)。 迭代次數(shù)K，終止搜索取決于所需的預(yù)測(cè)精度和，有參考文獻(xiàn)[15]：
（17）
其中表示上限的功能，完整的算法在如圖4中給出。
圖4 夾緊力的優(yōu)化算法（在示例1中使用）。 圖5 該算法在示例2使用
4． 加工過程中的夾緊力的優(yōu)化及測(cè)定
上一節(jié)介紹的算法可用于確定單負(fù)載作用于工件的載體的最佳夾緊力，然而，刀具路徑隨磨削量和切割點(diǎn)的不斷變化而變化。因此，相應(yīng)的夾緊力和最佳的加工負(fù)荷獲得將由圖4算法獲得，這大大增加了計(jì)算負(fù)擔(dān)，并要求為選擇的夾緊力提供標(biāo)準(zhǔn)， 將獲得滿意和適宜的整個(gè)刀具軌跡 ，用保守的辦法來解決下面將被討論的問題，考慮一個(gè)有限的數(shù)目（例如m）沿相應(yīng)的刀具路徑設(shè)置的產(chǎn)生m個(gè)最佳夾緊力，選擇記為， ， …,在每個(gè)采樣點(diǎn)，考慮以下四個(gè)最壞加工負(fù)荷向量：
(18)、和表示在、和方向上的最大值，、和上的數(shù)字1，2，3分別代替對(duì)應(yīng)的和另外兩個(gè)正交切削分力，而且有：
雖然4個(gè)最壞情況加工負(fù)荷向量不會(huì)在工件加工的同一時(shí)刻出現(xiàn)，但在每次常規(guī)的進(jìn)給速度中，刀具旋轉(zhuǎn)一次出現(xiàn)一次，負(fù)載向量引入的誤差可忽略。因此，在這項(xiàng)工作中，四個(gè)載體負(fù)載適用于同一位置，（但不是同時(shí)）對(duì)工件進(jìn)行的采樣 ，夾緊力的優(yōu)化算法圖4，對(duì)應(yīng)于每個(gè)采樣點(diǎn)計(jì)算最佳的夾緊力。夾緊力的最佳形式有：
(i=1,2,…,m) (j=x,y z,r) (19)
其中是最佳夾緊力的四個(gè)情況下的加工負(fù)荷載體，(C=1，2，…C)是每個(gè)相應(yīng)的夾具在第i個(gè)樣本點(diǎn)和第j負(fù)荷情況下力的大小。是計(jì)算每個(gè)負(fù)載點(diǎn)之后的結(jié)果，一套簡(jiǎn)單的“最佳”夾緊力必須從所有的樣本點(diǎn)和裝載條件里發(fā)現(xiàn)，并在所有的最佳夾緊力中選擇。這是通過在所有負(fù)載情況和采樣點(diǎn)排序，并選擇夾緊點(diǎn)的最高值的最佳的夾緊力，見于式 （20）：
（k=1，2，…，C） （20）
只要這些具備，就得到一套優(yōu)化的夾緊力，驗(yàn)證這些力，以確保工件夾具系統(tǒng)的靜態(tài)平衡。否則，會(huì)出現(xiàn)更多采樣點(diǎn)和重復(fù)上述程序。在這種方式中，可為整個(gè)刀具路徑確定“最佳”夾緊力 ，圖5總結(jié)了剛才所描述的算法。請(qǐng)注意，雖然這種方法是保守的，它提供了一個(gè)確定的夾緊力，最大限度地減少工件的定位誤差的一套系統(tǒng)方法。
5．影響工件的定位精度
它的興趣在于最早提出了評(píng)價(jià)夾緊力的算法對(duì)工件的定位精度的影響。工件首先放在與夾具接觸的基板上，然后夾緊力使工件接觸到夾具，因此，局部變形發(fā)生在每個(gè)工件夾具接觸處，使工件在夾具上移位和旋轉(zhuǎn)。隨后，準(zhǔn)靜態(tài)加工負(fù)荷應(yīng)用造成工件在夾具的移位。工件剛體運(yùn)動(dòng)的定義是由它在、和方向上的移位和自轉(zhuǎn)（見圖2），
如前所述，工件剛體位移產(chǎn)生于在每個(gè)夾緊處的局部變形，假設(shè)為相對(duì)于工件的質(zhì)量中心的第i個(gè)位置矢量定位點(diǎn)，坐標(biāo)變換定理可以用來表達(dá)在工件的位移，以及工件自轉(zhuǎn)如下: (21)
其中表示旋轉(zhuǎn)矩陣,描述當(dāng)?shù)卦诘趇幀相聯(lián)系的全球坐標(biāo)系和是一個(gè)旋轉(zhuǎn)矩陣確定工件相對(duì)于全球的坐標(biāo)系的定位坐標(biāo)系。假設(shè)夾具夾緊工件旋轉(zhuǎn)，由于旋轉(zhuǎn)很小，故也可近似為：
（22）
方程（21）現(xiàn)在可以改寫為： （23）
其中是經(jīng)方程（21）重新編排后變換得到的矩陣式，是夾緊和加工導(dǎo)致的工件剛體運(yùn)動(dòng)矢量。工件與夾具單方面接觸性質(zhì)意味著工件與夾具接觸處沒有拉力的可能。因此，在第i裝夾點(diǎn)接觸力可能與的關(guān)系如下：
（24）
其中是在第i個(gè)接觸點(diǎn)由于夾緊和加工負(fù)荷造成的變形，意味著凈壓縮變形，而負(fù)數(shù)則代表拉伸變形; 是表示在本地坐標(biāo)系第i個(gè)接觸剛度矩陣，是單位向量. 在這項(xiàng)研究中假定液壓/氣動(dòng)夾具，根據(jù)對(duì)外加工負(fù)荷，故在法線方向的夾緊力的強(qiáng)度保持不變，因此，必須對(duì)方程（24）的夾緊點(diǎn)進(jìn)行修改為：
（25）
其中是在第i個(gè)夾緊點(diǎn)的夾緊力，讓表示一個(gè)對(duì)外加工力量和載體的6×1矢量。并結(jié)合方程（23）—（25）與靜態(tài)平衡方程，得到下面的方程組：
（26）
其中，其中表示相乘。由于夾緊和加工工件剛體移動(dòng)，q可通過求解式（26）得到。工件的定位誤差向量， （見圖6），
現(xiàn)在可以計(jì)算如下： （27）
其中是考慮工件中心加工點(diǎn)的位置向量，且
6．模擬工作
較早前提出的算法是用來確定最佳夾緊力及其對(duì)兩例工件精度的影響例如：
1．適用于工件單點(diǎn)力。
2．應(yīng)用于工件負(fù)載準(zhǔn)靜態(tài)銑削序列
如左圖7 工件夾具配置中使用的模擬研究 工件夾具定位聯(lián)系; 、和全球坐標(biāo)系。
3-2-1夾具圖7所示，是用來定位并控制7075 - T6鋁合金（127毫米×127毫米×38.1毫米）的柱狀塊。假定為球形布局傾斜硬鋼定位器/夾具在表1中給出。工件——夾具材料的摩擦靜電對(duì)系數(shù)為0.25。使用伊利諾伊大學(xué)開發(fā)EMSIM程序[參考文獻(xiàn)26] 對(duì)加工瞬時(shí)銑削力條件進(jìn)行了計(jì)算，如表2給出例（1），應(yīng)用工件在點(diǎn)（109.2毫米，25.4毫米，34.3毫米）瞬時(shí)加工力，圖4中表3和表4列出了初級(jí)夾緊力和最佳夾緊力的算法 。該算法如圖5所示 ，一個(gè)25.4毫米銑槽使用EMSIM進(jìn)行了數(shù)值模擬，以減少起步（0.0毫米，25.4毫米，34.3毫米）和結(jié)束時(shí)（127.0毫米，25.4毫米，34.3毫米）四種情況下加工負(fù)荷載體，
（見圖8）。模擬計(jì)算銑削力數(shù)據(jù)在表5中給出。
圖8最終銑削過程模擬例如2。
表6中5個(gè)坐標(biāo)列出了為模擬抽樣調(diào)查點(diǎn)。最佳夾緊力是用前面討論過的排序算法計(jì)算每個(gè)采樣點(diǎn)和負(fù)載載體最后的夾緊力和負(fù)載。
7．結(jié)果與討論
例如算法1的繪制最佳夾緊力收斂圖9，
圖9
對(duì)于固定夾緊裝置在圖示例假設(shè)（見圖7），由此得到的夾緊力加權(quán)范數(shù)有如下形式：.結(jié)果表明，最佳夾緊力所述加工條件下有比初步夾緊力強(qiáng)度低得多的加權(quán)范數(shù)，最初的夾緊力是通過減少工件的夾具系統(tǒng)補(bǔ)充能量算法獲得。由于夾緊力和負(fù)載造成的工件的定位誤差，如表7。結(jié)果表明工件旋轉(zhuǎn)小，加工點(diǎn)減少錯(cuò)誤從13.1％到14.6％不等。在這種情況下，所有加工條件改善不是很大，因?yàn)閺淖畛跬ㄟ^互補(bǔ)勢(shì)能確定的最小化的夾緊力值已接近最佳夾緊力。圖5算法是用第二例在一個(gè)序列應(yīng)用于銑削負(fù)載到工件，他應(yīng)用于工件銑削負(fù)載一個(gè)序列。最佳的夾緊力，，對(duì)應(yīng)列表6每個(gè)樣本點(diǎn)，隨著最后的最佳夾緊力，在每個(gè)采樣點(diǎn)的加權(quán)范數(shù)和最優(yōu)的初始夾緊力繪圖10，在每個(gè)采樣點(diǎn)的加權(quán)范數(shù)的，，和繪制。
結(jié)果表明，由于每個(gè)組成部分是各相應(yīng)的最大夾緊力，它具有最高的加權(quán)范數(shù)。如圖10所示，如果在每個(gè)夾緊點(diǎn)最大組成部分是用于確定初步夾緊力，則夾緊力需相應(yīng)設(shè)置，有比相當(dāng)大的加權(quán)范數(shù)。故是一個(gè)完整的刀具路徑改進(jìn)方案。上述模擬結(jié)果表明，該方法可用于優(yōu)化夾緊力相對(duì)于初始夾緊力的強(qiáng)度，這種做法將減少所造成的夾緊力的加權(quán)范數(shù)，因此將提高工件的定位精度。
圖10
8．結(jié)論
該文件提出了關(guān)于確定多鉗夾具，工件受準(zhǔn)靜態(tài)加載系統(tǒng)的優(yōu)化加工夾緊力的新方法。夾緊力的優(yōu)化算法是基于接觸力學(xué)的夾具與工件系統(tǒng)模型，并尋求盡量減少應(yīng)用到所造成的工件夾緊力的加權(quán)范數(shù)，得出工件的定位誤差。該整體模型，制定一個(gè)雙目標(biāo)約束優(yōu)化問題，使用-約束的方法解決。該算法通過兩個(gè)模擬表明，涉及3-2-1型，二夾銑夾具的例子。今后的工作將解決在動(dòng)態(tài)負(fù)載存在夾具與工件在系統(tǒng)的優(yōu)化，其中慣性，剛度和阻尼效應(yīng)在確定工件夾具系統(tǒng)的響應(yīng)特性具有重要作用。
9．參考資料：
1、J. D. Lee 和L. S. Haynes .《柔性夾具系統(tǒng)的有限元分析》交易美國(guó)ASME，工程雜志工業(yè) ：134-139頁。
2、W. Cai, S. J. Hu 和J. X. Yuan .“柔性鈑金夾具：原理，算法和模擬”，交易美國(guó)ASME，制造科學(xué)與工程雜志 ：1996 318-324頁。
3、P. Chandra, S. M. Athavale, R. E. DeVor 和S. G. Kapoor.“負(fù)載對(duì)表面平整度的影響”工件夾具制造科學(xué)研討會(huì)論文集1996，第一卷：146-152頁。
4、R. J. Menassa 和V. R. DeVries.“適用于選拔夾具設(shè)計(jì)與優(yōu)化方法，美國(guó)ASME工業(yè)工程雜志：113 、 412-414，1991。
5、A. J. C. Trappey, C. Su 和J. Hou.《計(jì)算機(jī)輔助夾具分析中的應(yīng)用有限元分析和數(shù)學(xué)優(yōu)化模型》， 1995 ASME程序，MED： 777-787頁。
6、 S. N. Melkote, S. M. Athavale, R. E. DeVor, S. G. Kapoor 和J. Burkey .“基于加工過程仿真的加工裝置作用力系統(tǒng)研究”， NAMRI/SME：207–214頁, 1995
7、“考慮工件夾具，夾具接觸相互作用布局優(yōu)化模擬的結(jié)果” 341-346，1998。
8、E. C. DeMeter. 《快速支持布局優(yōu)化》，國(guó)際機(jī)床制造， 碩士論文 1998。
9、Y.-C. Chou, V. Chandru, M. M. Barash .《加工夾具機(jī)械構(gòu)造的數(shù)學(xué)算法：分析和合成》，美國(guó)ASME，工程學(xué)報(bào)工業(yè)“：1989 299-306頁。
10、S. H. Lee 和 M. R. Cutkosky. 《具有摩擦性的夾具規(guī)劃》 美國(guó)ASME，工業(yè)工程學(xué)報(bào)：1991，320–327頁。
11、S. Jeng, L. Chen 和W. Chieng.“最小夾緊力分析”，國(guó)際機(jī)床制造，碩士論文 1995年。
12、E. C. DeMeter.《加工夾具的性能的最小——最大負(fù)荷標(biāo)準(zhǔn)》 美國(guó)ASME，工業(yè)工程雜志 ：1994
13、E. C. DeMeter .《加工夾具最大負(fù)荷的性能優(yōu)化模型》 美國(guó)ASME，工業(yè)工程雜志 1995。
14、JH復(fù)和AYC倪.“核查和工件夾持的夾具設(shè)計(jì)”方案優(yōu)化，設(shè)計(jì)和制造，4，碩士論文： 307-318，1994。
15、T. H. Richards、埃利斯 霍伍德.1977,《應(yīng)力能量方法分析》，1977。
16、M. J. Hockenberger and E. C. DeMeter. 對(duì)工件準(zhǔn)靜態(tài)分析功能位移在加工夾具的應(yīng)用程序，制造科學(xué)雜志與工程： 325–331頁, 1996。
Int J Adv Manuf Technol (2001) 17:104113 2001 Springer-Verlag London Limited Fixture Clamping Force Optimisation and its Impact on Workpiece Location Accuracy B. Li and S. N. Melkote George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Georgia, USA Workpiece motion arising from localised elastic deformation at fixtureworkpiece contacts owing to clamping and machining forces is known to affect significantly the workpiece location accuracy and, hence, the final part quality. This effect can be minimised through fixture design optimisation. The clamping force is a critical design variable that can be optimised to reduce the workpiece motion. This paper presents a new method for determining the optimum clamping forces for a multiple clamp fixture subjected to quasi-static machining forces. The method uses elastic contact mechanics models to represent the fixtureworkpiece contact and involves the formulation and solution of a multi-objective constrained optimisation model. The impact of clamping force optimisation on workpiece location accuracy is analysed through examples involving a 32-1 type milling fixture. Keywords: Elastic contact modelling; Fixture clamping force; Optimisation 1. Introduction The location and immobilisation of the workpiece are two critical factors in machining. A machining fixture achieves these functions by locating the workpiece with respect to a suitable datum, and clamping the workpiece against it. The clamping force applied must be large enough to restrain the workpiece motion completely during machining. However, excessive clamping force can induce unacceptable level of workpiece elastic distortion, which will adversely affect its location and, in turn, the part quality. Hence, it is necessary to determine the optimum clamping forces that minimise the workpiece location error due to elastic deformation while satisfying the total restraint requirement. Previous researchers in the fixture analysis and synthesis area have used the finite-element (FE) modelling approach or Correspondence and offprint requests to: Dr S. N. Melkote, George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0405, USA. E-mail: shreyes.melkoteme.gatech.edu the rigid-body modelling approach. Extensive work based on the FE approach has been reported 18. With the exception of DeMeter 8, a common limitation of this approach is the large model size and computation cost. Also, most of the FE- based research has focused on fixture layout optimisation, and clamping force optimisation has not been addressed adequately. Several researchers have addressed fixture clamping force optimisation based on the rigid-body model 911. The rigid body modelling approach treats the fixture-element and work- piece as perfectly rigid solids. DeMeter 12, 13 used screw theory to solve for the minimum clamping force. The overall problem was formulated as a linear program whose objective was to minimise the normal contact force at each locating point by adjusting the clamping force intensity. The effect of the contact friction force was neglected because of its relatively small magnitude compared with the normal contact force. Since this approach is based on the rigid body assumption, it can uniquely only handle 3D fixturing schemes that involve no more than 6 unknowns. Fuh and Nee 14 also presented an iterative search-based method that computes the minimum clamping force by assuming that the friction force directions are known a priori. The primary limitation of the rigid-body analysis is that it is statically indeterminate when more than six contact forces are unknown. As a result, workpiece displace- ments cannot be determined uniquely by this method. This limitation may be overcome by accounting for the elasticity of the fixtureworkpiece system 15. For a relatively rigid workpiece, the location of the workpiece in the machining fixture is strongly influenced by the localised elastic defor- mation at the fixturing points. Hockenberger and DeMeter 16 used empirical contact force-deformation relations (called meta- functions) to solve for the workpiece rigid-body displacements due to clamping and quasi-static machining forces. The same authors also investigated the effect of machining fixture design parameters on workpiece displacement 17. Gui et al 18 reported an elastic contact model for improving workpiece location accuracy through optimisation of the clamping force. However, they did not address methods for calculating the fixtureworkpiece contact stiffness. In addition, the application of their algorithm for a sequence of machining loads rep- resenting a finite tool path was not discussed. Li and Melkote 19 and Hurtado and Melkote 20 used contact mechanics to Fixture Clamping Force Optimisation 105 solve for the contact forces and workpiece displacement pro- duced by the elastic deformation at the fixturing points owing to clamping loads. They also developed methods for optimising the fixture layout 21 and clamping force using this method 22. However, clamping force optimisation for a multiclamp system and its impact on workpiece accuracy were not covered in these papers. This paper presents a new algorithm based on the contact elasticity method for determining the optimum clamping forces for a multiclamp fixtureworkpiece system subjected to quasi- static loads. The method seeks to minimise the impact of workpiece motion due to clamping and machining loads on the part location accuracy by systematically optimising the clamping forces. A contact mechanics model is used to deter- mine a set of contact forces and displacements, which are then used for the clamping force optimisation. The complete prob- lem is formulated and solved as a multi-objective constrained optimisation problem. The impact of clamping force optimis- ation on workpiece location accuracy is analysed via two examples involving a 32-1 fixture layout for a milling oper- ation. 2. FixtureWorkpiece Contact Modelling 2.1 Modelling Assumptions The machining fixture consists of L locators and C clamps with spherical tips. The workpiece and fixture materials are linearly elastic in the contact region, and perfectly rigid else- where. The workpiecefixture system is subjected to quasi- static loads due to clamping and machining. The clamping force is assumed to be constant during machining. This assumption is valid when hydraulic or pneumatic clamps are used. In reality, the elasticity of the fixtureworkpiece contact region is distributed. However, in this model development, lumped contact stiffness is assumed (see Fig. 1). Therefore, the contact force and localised deformation at the ith fixturing point can be related as follows: F i j = k i j d i j (1) where k i j (j = x,y,z) denotes the contact stiffness in the tangential and normal directions of the local x i ,y i ,z i coordinate frame, d i j Fig. 1. A lumped-spring fixtureworkpiece contact model. x i , y i , z i , denote the local coordinate frame at the ith contact. (j = x,y,z) are the corresponding localised elastic deformations along the x i ,y i , and z i axes, respectively, F i j (j = x,j,z) represents the local contact force components with F i x and F i y being the local x i and y i components of the tangential force, and F i z the normal force. 2.2 WorkpieceFixture Contact Stiffness Model The lumped compliance at a spherical tip locator/clamp and workpiece contact is not linear because the contact radius varies nonlinearly with the normal force 23. The contact deformation due to the normal force P i acting between a spherical tipped fixture element of radius R i and a planar workpiece surface can be obtained from the closed-form Hertz- ian solution to the problem of a sphere indenting an elastic half-space. For this problem, the normal deformation D i n is given as 23, p. 93: D i n = S 9(P i ) 2 16R i (E*) 2 D 1/3 (2) where 1 E* = 1 - n 2 w E w + 1 - n 2 f E f E w and E f are Youngs moduli for the workpiece and fixture materials, respectively, and n w and n f are Poisson ratios for the workpiece and fixture materials, respectively. The tangential deformation D i t (= D i tx or D i ty in the local x i and y i tangential directions, respectively) due to a tangential force Q i (= Q i x or Q i y ) has the following form 23, p. 217: D ti t = Q i 8a i S 2 - n f G f + 2 - n w G w D (3) where a i = S 3P i R i 4 S 1 - n f E f + 1 - n w E w DD 1/3 and G w and G f are shear moduli for the workpiece and fixture materials, respectively. A reasonable linear approximation of the contact stiffness can be obtained from a least-squares fit to Eq. (2). This yields the following linearised contact stiffness values: k i z = 8.82 S 16R i (E*) 2 9 D 1/3 (4) k i x = k i y = 4 E* S 2 - n j G f + 2 - n w G w D - 1 k i z (5) In deriving the above linear approximation, the normal force P i was assumed to vary from 0 to 1000 N, and the correspond- ing R 2 value of the least-squares fit was found to be 0.94. 3. Clamping Force Optimisation The goal is to determine the set of optimal clamping forces that will minimise the workpiece rigid-body motion due to 106 B. Li and S. N. Melkote localised elastic deformation induced by the clamping and machining loads, while maintaining the fixtureworkpiece sys- tem in quasi-static equilibrium during machining. Minimisation of the workpiece motion will, in turn, reduce the location error. This goal is achieved by formulating the problem as a multi- objective constrained optimisation problem, as described next. 3.1 Objective Function Formulation Since the workpiece rotation due to fixturing forces is often quite small 17 the workpiece location error is assumed to be determined largely by its rigid-body translation Dd w = DX w DY w DZ w T , where DX w , DY w , and DZ w are the three orthogonal components of Dd w along the X g , Y g , and Z g axes (see Fig. 2). The workpiece location error due to the fixturing forces can then be calculated in terms of the L 2 norm of the rigid-body displacement as follows: iDd w i = (DX w ) 2 + (DY w ) 2 + (DZ w ) 2 ) (6) where ii denotes the L 2 norm of a vector. In particular, the resultant clamping force acting on the workpiece will adversely affect the location error. When mul- tiple clamping forces are applied to the workpiece, the resultant clamping force, P R C = P R X P R y P R Z T , has the form: P R C = R C P C (7) where P C = P L+1 .P L+C T is the clamping force vector, R C = n L+1 .n L+C T is the clamping force direction matrix, n L+i = cosa L+i cosb L+i cosg L+i T is the clamping force direction cosine vector, and a L+i , b L+i , and g L+i are angles made by the clamping force vector at the ith clamping point with respect to the X g , Y g , Z g coordinate axes (i = 1,2,. . .,C). In this paper, the workpiece location error due to contact region deformation is assumed to be influenced only by the normal force acting at the locatorworkpiece contacts. The frictional force at the contacts is relatively small and is neg- lected when analysing the impact of the clamping force on the workpiece location error. Denoting the ratio of the normal contact stiffness, k i z , to the smallest normal stiffness among all locators, k s z ,byj i (i = 1,. . .,L), and assuming that the workpiece rests on N X , N Y , and N Z number of locators oriented in the X g , Fig. 2. Workpiece rigid body translation and rotation. Y g , and Z g directions, the equivalent contact stiffness in the X g , Y g , and Z g directions can be calculated as k s zSO N X i=1 j iD , k s zSO N Y i=1 j iD , and k s zSO N Z i=1 j iD respectively (see Fig. 3). The workpiece rigid-body motion, Dd w , due to clamping action can now be written as: Dd w = 3 P R X k s zSO N X i=1 j iD P R Y k s zSO N Y i=1 j iD P R Z k s z SO N Z i=1 j iD 4 T (8) The workpiece motion, and hence the location error can be reduced by minimising the weighted L 2 norm of the resultant clamping force vector. Therefore, the first objective function can be written as: Minimize iP R C i w = ! 11 P R X O N X i=1 j i 2 2 + 1 P R Y O N Y i=1 j i 2 2 + 1 P R Z O N Z i=1 j i 2 2 2 (9) Note that the weighting factors are proportional to the equival- ent contact stiffnesses in the X g , Y g , and Z g directions. The components of P R C are uniquely determined by solving the contact elasticity problem using the principle of minimum total complementary energy 15, 23. This ensures that the clamping forces and the corresponding locator reactions are “true” solutions to the contact problem and yield “true” rigid- body displacements, and that the workpiece is kept in static equilibrium by the clamping forces at all times. Therefore, the minimisation of the total complementary energy forms the second objective function for the clamping force optimisation and is given by: Minimise (U* - W*) = 1 2 FO L+C i=1 (F i x ) 2 k i x + O L+C i=1 (F i y ) 2 k i y + O L+C i=1 (F i z ) 2 k i z G (10) = .l T Ql Fig. 3. The basis for the determination of the weighting factor for the L 2 norm calculation. Fixture Clamping Force Optimisation 107 where U* represents the complementary strain energy of the elastically deformed bodies, W* represents the complementary work done by the external force and moments, Q = diag c 1 x c 1 y c 1 z .c L+C x c L+C y c L+C z is the diagonal contact compliance matrix, c i j = (k i j ) - 1 , and l = F 1 x F 1 y F 1 z .F L+C x F L+C y F L+C z T is the vector of all contact forces. 3.2 Friction and Static Equilibrium Constraints The optimisation objective in Eq. (10) is subject to certain constraints and bounds. Foremost among them is the static friction constraint at each contact. Coulombs friction law states that (F i x ) 2 +(F i y ) 2 ) #m i s F i z (m i s is the static friction coefficient). A conservative and linearised version of this nonlinear con- straint can be used and is given by 19: uF i x u + uF i y u #m i s F i z (11) Since quasi-static loads are assumed, the static equilibrium of the workpiece is ensured by including the following force and moment equilibrium equations (in vector form): O F = 0 (12) O M = 0 where the forces and moments consist of the machining forces, workpiece weight and the contact forces in the normal and tangential directions. 3.3 Bounds Since the fixtureworkpiece contact is strictly unilateral, the normal contact force, P i , can only be compressive. This is expressed by the following bound on P i : P i $ 0(i = 1, . . ., L + C) (13) where it is assumed that normal forces directed into the workpiece are positive. In addition, the normal compressive stress at a contact cannot exceed the compressive yield strength (S y ) of the workpiece material. This upper bound is written as: P i # S y A i (i = 1, . . .,L+C) (14) where A i is the contact area at the ith workpiecefixture con- tact. The complete clamping force optimisation model can now be written as: Minimize f = H f 1 f 2 J = H .l T Ql iP R C i w J (15) subject to: (11)(14). 4. Algorithm for Model Solution The multi-objective optimisation problem in Eq. (15) can be solved by the e-constraint method 24. This method identifies one of the objective functions as primary, and converts the other into a constraint. In this work, the minimisation of the complementary energy (f 1 ) is treated as the primary objective function, and the weighted L 2 norm of the resultant clamping force (f 2 ) is treated as a constraint. The choice of f 1 as the primary objective ensures that a unique set of feasible clamping forces is selected. As a result, the workpiecefixture system is driven to a stable state (i.e. the minimum energy state) that also has the smallest weighted L 2 norm for the resultant clamping force. The conversion of f 2 into a constraint involves specifying the weighted L 2 norm to be less than or equal to e, where e is an upper bound on f 2 . To determine a suitable e,itis initially assumed that all clamping forces are unknown. The contact forces at the locating and clamping points are computed by considering only the first objective function (i.e. f 1 ). While this set of contact forces does not necessarily yield the lowest clamping forces, it is a “true” feasible solution to the contact elasticity problem that can completely restrain the workpiece in the fixture. The weighted L 2 norm of these clamping forces is computed and taken as the initial value of e. Therefore, the clamping force optimisation problem in Eq. (15) can be rewritten as: Minimize f 1 = .l T Ql (16) subject to: iP R C i w $e, (11)(14). An algorithm similar to the bisection method for finding roots of an equation is used to determine the lowest upper bound for iP R C i w . By decreasing the upper bound e as much as possible, the minimum weighted L 2 norm of the resultant clamping force is obtained. The number of iterations, K, needed to terminate the search depends on the required prediction accuracy d and ueu, and is given by 25: K = F log 2 S ueu d DG (17) where I denotes the ceiling function. The complete algorithm is given in Fig. 4. 5. Determination of Optimum Clamping Forces During Machining The algorithm presented in the previous section can be used to determine the optimum clamping force for a single load vector applied to the workpiece. However, during milling the magnitude and point of cutting force application changes continuously along the tool path. Therefore, an infinite set of optimum clamping forces corresponding to the infinite set of machining loads will be obtained with the algorithm of Fig. 4. This substantially increases the computational burden and calls for a criterion/procedure for selecting a single set of clamping forces that will be satisfactory and optimum for the entire tool path. A conservative approach to addressing these issues is discussed next. Consider a finite number (say m) of sample points along the tool path yielding m corresponding sets of optimum clamp- ing forces denoted as P 1 opt , P 2 opt ,.,P m opt . At each sampling 108 B. Li and S. N. Melkote Fig. 4. Clamping force optimisation algorithm (used in example 1). point, the following four worst-case machining load vectors are considered: F X max = F max X F 1 Y F 1 Z T F Y max = F 2 X F max Y F 2 Z T F Z max = F 3 X F 3 Y F max Z T (18) F r max = F 4 X F 4 Y F 4 Z T where F max X , F max Y , and F max Z are the maximum X g , Y g , and Z g components of the machining force, the superscripts 1, 2, 3 of F X , F Y , and F Z stand for the other two orthogonal machining force components corresponding to F max X , F max Y , and F max Z , respectively, and iF r max i = max(F X ) 2 +(F Y ) 2 +(F Z ) 2 ). Although the four worst-case machining load vectors will not act on the workpiece at the same instant, they will occur once per cutter revolution. At conventional feedrates, the error introduced by applying the load vectors at the same point would be negligible. Therefore, in this work, the four load vectors are applied at the same location (but not simultaneously) on the workpiece corresponding to the sam- pling instant. The clamping force optimisation algorithm of Fig. 4 is then used to calculate the optimum clamping forces corresponding to each sampling point. The optimum clamping forces have the form: P i jmax = C i 1j C i 2j .C i Cj T (i = 1, . . .,m)(j = x,y,z,r) (19) where P i jmax is the vector of optimum clamping forces for the four worst-case machining load vectors, and C i kj (k = 1,. . .,C) is the force magnitude at each clamp corresponding to the ith sample point and the jth load scenario. After P i jmax is computed for each load application point, a single set of “optimum” clamping forces must be selected from all of the optimum clamping forces found for each clamp from all the sample points and loading conditions. This is done by sorting the optimum clamping force magnitudes at a clamping point for all load scenarios and sample points and selecting the maximum value, C max k , as given in Eq. (20): C max k # C i kj (k =