喜歡就充值下載吧，，資源目錄下展示的全都有，，下載后全都有，所見即所得，歡迎充值購買哦======================喜歡就充值下載吧，，資源目錄下展示的全都有，，下載后全都有，所見即所得，歡迎充值購買哦======================喜歡就充值下載吧，，資源目錄下展示的全都有，，下載后全都有，所見即所得，歡迎充值購買哦====================== 夾具夾緊力的優(yōu)化及對工件定位精度的影響 B.Li 和 S.N.Mellkote 布什伍德拉夫機械工程學院，佐治亞理工學院，格魯吉亞，美國研究所 由于夾緊和加工，在工件和夾具的接觸部位會產(chǎn)生局部彈性變形，使工件尺寸發(fā)生變化，進而影響工件的最終加工質(zhì)量。這種效應可通過最小化夾具設計優(yōu)化，夾緊力是一個重要的設計變量，可以得到優(yōu)化，以減少工件的位移。本文提出了一種確定多夾緊夾具受到準靜態(tài)加工部位的最佳夾緊力的新方法。該方法采用彈性接觸力學模型代表夾具與工件接觸，并涉及制定和解決方案的多目標優(yōu)化模型的約束。夾緊力的最優(yōu)化對工件定位精度的影響通過3-2-1式銑夾具的例子進行了分析。 關鍵詞：彈性 接觸 模型 夾具 夾緊力 優(yōu)化 前言 定位和夾緊的工件加工中的兩個關鍵因素。要實現(xiàn)夾具的這些功能，需將工件定位到一個合適的基準上并夾緊，采用的夾緊力必須足夠大，以抑制工件在加工過程中產(chǎn)生的移動。然而，過度的夾緊力可誘導工件產(chǎn)生更大的彈性變形 ，這會影響它的位置精度，并反過來影響零件質(zhì)量。所以有必要確定最佳夾緊力，來減小由于彈性變形對工件的定位誤差，同時滿足加工的要求。在夾具分析和綜合領域上的研究人員使用了有限元模型的方法或剛體模型的方法。大量的工作都以有限元方法為基礎被報道[參考文獻1-8]。隨著得墨忒耳[8]，這種方法的限制是需要較大的模型和計算成本。同時，多數(shù)的有限元基礎研究人員一直重點關注的夾具布局優(yōu)化和夾緊力的優(yōu)化還沒有得到充分討論，也有少數(shù)的研究人員通過對剛性模型[9-11]對夾緊力進行了優(yōu)化，剛型模型幾乎被近似為一個規(guī)則完整的形狀。得墨忒耳[12，13]用螺釘理論解決的最低夾緊力，總的問題是制定一個線性規(guī)劃，其目的是盡量減少在每個定位點調(diào)整夾緊力強度的法線接觸力。接觸摩擦力的影響被忽視，因為它較法線接觸力相對較小，由于這種方法是基于剛體假設，獨特的三維夾具可以處理超過6個自由度的裝夾，復和倪[14]也提出迭代搜索方法，通過假設已知摩擦力的方向來推導計算最小夾緊力，該剛體分析的主要限制因素是當出現(xiàn)六個以上的接觸力是使其靜力不確定，因此，這種方法無法確定工件移位的唯一性。 這種限制可以通過計算夾具——工件系統(tǒng)[15]的彈性來克服，對于一個相對嚴格的工件，該夾具在機械加工工件的位置會受夾具點的局部彈性變形的強烈影響。Hockenberger和得墨忒耳[16]使用經(jīng)驗的接觸力變形的關系（稱為元功能），解決由于夾緊和準靜態(tài)加工力工件剛體位移。同一作者還考察了加工工件夾具位移對設計參數(shù)的影響[17]。桂 [18] 等 通過工件的夾緊力的優(yōu)化定位精度彈性接觸模型對報告做了改善，然而，他們沒有處理計算夾具與工件的接觸剛度的方法，此外，其算法的應用沒有討論機械加工刀具路徑負載有限序列。李和Melkote [19]和烏爾塔多和Melkote [20]用接觸力學解決由于在加載夾具夾緊點彈性變形產(chǎn)生的接觸力和工件的位移，他們還使用此方法制定了優(yōu)化方法夾具布局[21]和夾緊力[22]。但是，關于multiclamp系統(tǒng)及其對工件精度影響的夾緊力的優(yōu)化并沒有在這些文件中提到 。 本文提出了一種新的算法，確定了multiclamp夾具工件系統(tǒng)受到準靜態(tài)加載的最佳夾緊力為基礎的彈性方法。該法旨在盡量減少影響由于工件夾緊位移和加工荷載通過系統(tǒng)優(yōu)化夾緊力的一部分定位精度。接觸力學模型，用于確定接觸力和位移，然后再用做夾緊力優(yōu)化，這個問題被作為多目標約束優(yōu)化問題提出和解決。通過兩個例子分析工件夾緊力的優(yōu)化對定位精度的影響，例子涉及的銑削夾具3-2-1布局。 1． 夾具——工件聯(lián)系模型 1．1 模型假設 該加工夾具由L定位器和帶有球形端的c形夾組成。工件和夾具接觸的地方是線性的彈性接觸，其他地方完全剛性。工件——夾具系統(tǒng)由于夾緊和加工受到準靜態(tài)負載。夾緊力可假定為在加工過程中保持不變，這個假設是有效的，在對液壓或氣動夾具使用。在實際中，夾具工件接觸區(qū)域是彈性分布，然而，這種模式的發(fā)展，假設總觸剛度（見圖1）第i夾具接觸力局部變形如下： (1) 其中（j=x，y，z）表示，在當?shù)刈幼鴺讼登芯€和法線方向的接觸剛度 第 19 頁 共 15 頁 圖1 彈簧夾具—— 工件接觸模型。 表示在第i個 接觸處的坐標系 （j=x，y，z）是對應沿著xyz方向的彈性變形，分別 （j= x，y，z）的代表和切向力接觸 ，法線力接觸。 1．2 工件——夾具的接觸剛度模型 集中遵守一個球形尖端定位，夾具和工件的接觸并不是線性的，因為接觸半徑與隨法線力呈非線性變化 [23]。由于法線力接觸變形作用于半徑和平面工件表面之間，這可從封閉赫茲的辦法解決縮進一個球體彈性半空間的問題。對于這個問題， 是法線的變形，在[文獻23 第93頁]中給出如下： （2） 其中式中 和是工件和夾具的彈性模量，、分別是工件和材料的泊松比。 切向變形沿著和切線方向）硅業(yè)切力距有以下形式[文獻23第217頁] （3） 其中、 分別是工件和夾具剪切模量 一個合理的接觸剛度的線性可以近似從最小二乘獲得適合式 （2），這就產(chǎn)生了以下線性化接觸剛度值：在計算上述的線性近似， （4） (5) 正常的力被假定為從0到1000N，且最小二乘擬合相應的R2值認定是0.94。 2．夾緊力優(yōu)化 我們的目標是確定最優(yōu)夾緊力，將盡量減少由于工件剛體運動過程中，局部的夾緊和加工負荷引起的彈性變形，同時保持在準靜態(tài)加工過程中夾具——工件系統(tǒng)平衡，工件的位移減少，從而減少定位誤差。實現(xiàn)這個目標是通過制定一個多目標約束優(yōu)化問題的問題，如下描述。 2.1 目標函數(shù)配方 工件旋轉(zhuǎn)，由于部隊輪換往往是相當小[17]的工件定位誤差假設為確定其剛體翻譯基本上，其中 、、和 是 沿，和三個正交組件（見圖2）。 圖2 工件剛體平移和旋轉(zhuǎn) 工件的定位誤差歸于裝夾力，然后可以在該剛體位移的范數(shù)計算如下： （6） 其中表示一個向量二級標準。 但是作用在工件的夾緊力會影響定位誤差。當多個夾緊力作用于工件，由此產(chǎn)生的夾緊力為,有如下形式： （7） 其中夾緊力是矢量，夾緊力的方向矩陣，是夾緊力是矢量的方向余弦，、和 是第i個夾緊點夾緊力在、和方向上的向量角度（i=1、2、3...,C）。 在這個文件中，由于接觸區(qū)變形造成的工件的定位誤差，被假定為受的作用力是法線的，接觸的摩擦力相對較小，并在進行分析時忽略了加緊力對工件的定位誤差的影響。意指正常接觸剛度比，是通過（i=1，2…L）和最小的所有定位器正常剛度相乘，并假設工件、、取決于、、的方向，各自的等效接觸剛度可有下式計算得出（見圖3），工件剛體運動，歸于夾緊行動現(xiàn)在可以寫成： (8) 工件有位移，因此，定位誤差的減小可以通過盡量減少產(chǎn)生的夾緊力向量 范數(shù)。因此，第一個目標函數(shù)可以寫為： 最小化 （9） 要注意，加權(quán)因素是與等效接觸剛度成正比的在、和 方向上。通過使用最低總能量互補參考文獻[15，23]的原則求解彈性力學接觸問題得出A的組成部分是唯一確定的，這保證了夾緊力和相應的定位反應是“真正的”解決方案，對接觸問題和產(chǎn)生的“真正”剛體位移，而且工件保持在靜態(tài)平衡，通過夾緊力的隨時調(diào)整。因此，總能量最小化的形式為補充的夾緊力優(yōu)化的第二個目標函數(shù)，并給出： 最小化 （10） 其中代表機構(gòu)的彈性變形應變能互補，代表由外部力量和力矩配合完成，是遵守對角矩陣的， 和是所有接觸力的載體。 如圖3 加權(quán)系數(shù)計算確定的基礎 內(nèi)蒙古科技大學本科生畢業(yè)設計（外文翻譯） 2.2 摩擦和靜態(tài)平衡約束 在（10）式優(yōu)化的目標受到一定的限制和約束，他們中最重要的是在每個接觸處的靜摩擦力約束。庫侖摩擦力的法律規(guī)定（是靜態(tài)摩擦系數(shù)）,這方面的一個非線性約束和線性化版本可以使用，并且[19]有： （11） 假設準靜態(tài)載荷，工件的靜力平衡由下列力和力矩平衡方程確保（向量形式）： (12) 其中包括在法線和切線方向的力和力矩的機械加工力和工件重量。 2.3界接觸力 由于夾具——工件接觸是單側(cè)面的，法線的接觸力只能被壓縮。這通過以下的的約束表（i=1，2…,L+C） （13） 它假設在工件上的法線力是確定的，此外，在一個法線的接觸壓力不能超過壓工件材料的屈服強度（）。這個約束可寫為： （i=1，2，…,L+C） (14) 如果是在第i個工件——夾具的接觸處的接觸面積，完整的夾緊力優(yōu)化模型，可以寫成：最小化 (15) 3．模型算法求解 式（15）多目標優(yōu)化問題可以通過求解約束[24]。這種方法將確定的目標作為首要職能之一，并將其轉(zhuǎn)換成一個約束對。該補充（）的主要目的是處理功能，并由此得到夾緊力（）作為約束的加權(quán)范數(shù)最小化。對為主要目標的選擇，確保選中一套獨特可行的夾緊力，因此，工件——夾具系統(tǒng)驅(qū)動到一個穩(wěn)定的狀態(tài)（即最低能量狀態(tài)），此狀態(tài)也表示有最小的夾緊力下的加權(quán)范數(shù)。 的約束轉(zhuǎn)換涉及到一個指定的加權(quán)范數(shù)小于或等于，其中是 的約束，假設最初所有夾緊力不明確，要確定一個合適的。在定位和夾緊點的接觸力的計算只考慮第一個目標函數(shù)（即）。雖然有這樣的接觸力，并不一定產(chǎn)生最低的夾緊力，這是一個“真正的”可行的解決彈性力學問題辦法，可完全抑制工件在夾具中的位置。這些夾緊力的加權(quán)系數(shù)，通過計算并作為初始值與比較，因此，夾緊力式（15）的優(yōu)化問題可改寫為: 最小化 （16） 由： (11)–(14) 得。 類似的算法尋找一個方程根的二分法來確定最低的上的約束, 通過盡可能降低上限，由此產(chǎn)生的最小夾緊力的加權(quán)范數(shù)。 迭代次數(shù)K，終止搜索取決于所需的預測精度和，有參考文獻[15]： （17） 其中表示上限的功能，完整的算法在如圖4中給出。 圖4 夾緊力的優(yōu)化算法（在示例1中使用）。 圖5 該算法在示例2使用 4． 加工過程中的夾緊力的優(yōu)化及測定 上一節(jié)介紹的算法可用于確定單負載作用于工件的載體的最佳夾緊力，然而，刀具路徑隨磨削量和切割點的不斷變化而變化。因此，相應的夾緊力和最佳的加工負荷獲得將由圖4算法獲得，這大大增加了計算負擔，并要求為選擇的夾緊力提供標準， 將獲得滿意和適宜的整個刀具軌跡 ，用保守的辦法來解決下面將被討論的問題，考慮一個有限的數(shù)目（例如m）沿相應的刀具路徑設置的產(chǎn)生m個最佳夾緊力，選擇記為， ， …,在每個采樣點，考慮以下四個最壞加工負荷向量： (18)、和表示在、和方向上的最大值，、和上的數(shù)字1，2，3分別代替對應的和另外兩個正交切削分力，而且有： 雖然4個最壞情況加工負荷向量不會在工件加工的同一時刻出現(xiàn)，但在每次常規(guī)的進給速度中，刀具旋轉(zhuǎn)一次出現(xiàn)一次，負載向量引入的誤差可忽略。因此，在這項工作中，四個載體負載適用于同一位置，（但不是同時）對工件進行的采樣 ，夾緊力的優(yōu)化算法圖4，對應于每個采樣點計算最佳的夾緊力。夾緊力的最佳形式有： (i=1,2,…,m) (j=x,y z,r) (19) 其中是最佳夾緊力的四個情況下的加工負荷載體，(C=1，2，…C)是每個相應的夾具在第i個樣本點和第j負荷情況下力的大小。是計算每個負載點之后的結(jié)果，一套簡單的“最佳”夾緊力必須從所有的樣本點和裝載條件里發(fā)現(xiàn)，并在所有的最佳夾緊力中選擇。這是通過在所有負載情況和采樣點排序，并選擇夾緊點的最高值的最佳的夾緊力，見于式 （20）： （k=1，2，…，C） （20） 只要這些具備，就得到一套優(yōu)化的夾緊力，驗證這些力，以確保工件夾具系統(tǒng)的靜態(tài)平衡。否則，會出現(xiàn)更多采樣點和重復上述程序。在這種方式中，可為整個刀具路徑確定“最佳”夾緊力 ，圖5總結(jié)了剛才所描述的算法。請注意，雖然這種方法是保守的，它提供了一個確定的夾緊力，最大限度地減少工件的定位誤差的一套系統(tǒng)方法。 5．影響工件的定位精度 它的興趣在于最早提出了評價夾緊力的算法對工件的定位精度的影響。工件首先放在與夾具接觸的基板上，然后夾緊力使工件接觸到夾具，因此，局部變形發(fā)生在每個工件夾具接觸處，使工件在夾具上移位和旋轉(zhuǎn)。隨后，準靜態(tài)加工負荷應用造成工件在夾具的移位。工件剛體運動的定義是由它在、和方向上的移位和自轉(zhuǎn)（見圖2）， 如前所述，工件剛體位移產(chǎn)生于在每個夾緊處的局部變形，假設為相對于工件的質(zhì)量中心的第i個位置矢量定位點，坐標變換定理可以用來表達在工件的位移，以及工件自轉(zhuǎn)如下: (21) 其中表示旋轉(zhuǎn)矩陣,描述當?shù)卦诘趇幀相聯(lián)系的全球坐標系和是一個旋轉(zhuǎn)矩陣確定工件相對于全球的坐標系的定位坐標系。假設夾具夾緊工件旋轉(zhuǎn)，由于旋轉(zhuǎn)很小，故也可近似為： （22） 方程（21）現(xiàn)在可以改寫為： （23） 其中是經(jīng)方程（21）重新編排后變換得到的矩陣式，是夾緊和加工導致的工件剛體運動矢量。工件與夾具單方面接觸性質(zhì)意味著工件與夾具接觸處沒有拉力的可能。因此，在第i裝夾點接觸力可能與的關系如下： （24） 其中是在第i個接觸點由于夾緊和加工負荷造成的變形，意味著凈壓縮變形，而負數(shù)則代表拉伸變形; 是表示在本地坐標系第i個接觸剛度矩陣，是單位向量. 在這項研究中假定液壓/氣動夾具，根據(jù)對外加工負荷，故在法線方向的夾緊力的強度保持不變，因此，必須對方程（24）的夾緊點進行修改為： （25） 其中是在第i個夾緊點的夾緊力，讓表示一個對外加工力量和載體的6×1矢量。并結(jié)合方程（23）—（25）與靜態(tài)平衡方程，得到下面的方程組： （26） 其中，其中表示相乘。由于夾緊和加工工件剛體移動，q可通過求解式（26）得到。工件的定位誤差向量， （見圖6）， 現(xiàn)在可以計算如下： （27） 其中是考慮工件中心加工點的位置向量，且 6．模擬工作 較早前提出的算法是用來確定最佳夾緊力及其對兩例工件精度的影響例如： 1．適用于工件單點力。 2．應用于工件負載準靜態(tài)銑削序列 如左圖7 工件夾具配置中使用的模擬研究 工件夾具定位聯(lián)系; 、和全球坐標系。 3-2-1夾具圖7所示，是用來定位并控制7075 - T6鋁合金（127毫米×127毫米×38.1毫米）的柱狀塊。假定為球形布局傾斜硬鋼定位器/夾具在表1中給出。工件——夾具材料的摩擦靜電對系數(shù)為0.25。使用伊利諾伊大學開發(fā)EMSIM程序[參考文獻26] 對加工瞬時銑削力條件進行了計算，如表2給出例（1），應用工件在點（109.2毫米，25.4毫米，34.3毫米）瞬時加工力，圖4中表3和表4列出了初級夾緊力和最佳夾緊力的算法 。該算法如圖5所示 ，一個25.4毫米銑槽使用EMSIM進行了數(shù)值模擬，以減少起步（0.0毫米，25.4毫米，34.3毫米）和結(jié)束時（127.0毫米，25.4毫米，34.3毫米）四種情況下加工負荷載體， （見圖8）。模擬計算銑削力數(shù)據(jù)在表5中給出。 圖8最終銑削過程模擬例如2。 表6中5個坐標列出了為模擬抽樣調(diào)查點。最佳夾緊力是用前面討論過的排序算法計算每個采樣點和負載載體最后的夾緊力和負載。 7．結(jié)果與討論 例如算法1的繪制最佳夾緊力收斂圖9， 圖9 對于固定夾緊裝置在圖示例假設（見圖7），由此得到的夾緊力加權(quán)范數(shù)有如下形式：.結(jié)果表明，最佳夾緊力所述加工條件下有比初步夾緊力強度低得多的加權(quán)范數(shù)，最初的夾緊力是通過減少工件的夾具系統(tǒng)補充能量算法獲得。由于夾緊力和負載造成的工件的定位誤差，如表7。結(jié)果表明工件旋轉(zhuǎn)小，加工點減少錯誤從13.1％到14.6％不等。在這種情況下，所有加工條件改善不是很大，因為從最初通過互補勢能確定的最小化的夾緊力值已接近最佳夾緊力。圖5算法是用第二例在一個序列應用于銑削負載到工件，他應用于工件銑削負載一個序列。最佳的夾緊力，，對應列表6每個樣本點，隨著最后的最佳夾緊力，在每個采樣點的加權(quán)范數(shù)和最優(yōu)的初始夾緊力繪圖10，在每個采樣點的加權(quán)范數(shù)的，，和繪制。 結(jié)果表明，由于每個組成部分是各相應的最大夾緊力，它具有最高的加權(quán)范數(shù)。如圖10所示，如果在每個夾緊點最大組成部分是用于確定初步夾緊力，則夾緊力需相應設置，有比相當大的加權(quán)范數(shù)。故是一個完整的刀具路徑改進方案。上述模擬結(jié)果表明，該方法可用于優(yōu)化夾緊力相對于初始夾緊力的強度，這種做法將減少所造成的夾緊力的加權(quán)范數(shù)，因此將提高工件的定位精度。 圖10 8．結(jié)論 該文件提出了關于確定多鉗夾具，工件受準靜態(tài)加載系統(tǒng)的優(yōu)化加工夾緊力的新方法。夾緊力的優(yōu)化算法是基于接觸力學的夾具與工件系統(tǒng)模型，并尋求盡量減少應用到所造成的工件夾緊力的加權(quán)范數(shù)，得出工件的定位誤差。該整體模型，制定一個雙目標約束優(yōu)化問題，使用-約束的方法解決。該算法通過兩個模擬表明，涉及3-2-1型，二夾銑夾具的例子。今后的工作將解決在動態(tài)負載存在夾具與工件在系統(tǒng)的優(yōu)化，其中慣性，剛度和阻尼效應在確定工件夾具系統(tǒng)的響應特性具有重要作用。 9．參考資料： 1、J. D. Lee 和L. S. Haynes .《柔性夾具系統(tǒng)的有限元分析》交易美國ASME，工程雜志工業(yè) ：134-139頁。 2、W. Cai, S. J. Hu 和J. X. Yuan .“柔性鈑金夾具：原理，算法和模擬”，交易美國ASME，制造科學與工程雜志 ：1996 318-324頁。 3、P. Chandra, S. M. Athavale, R. E. DeVor 和S. G. Kapoor.“負載對表面平整度的影響”工件夾具制造科學研討會論文集1996，第一卷：146-152頁。 4、R. J. Menassa 和V. R. DeVries.“適用于選拔夾具設計與優(yōu)化方法，美國ASME工業(yè)工程雜志：113 、 412-414，1991。 5、A. J. C. Trappey, C. Su 和J. Hou.《計算機輔助夾具分析中的應用有限元分析和數(shù)學優(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)化》，國際機床制造， 碩士論文 1998。 9、Y.-C. Chou, V. Chandru, M. M. Barash .《加工夾具機械構(gòu)造的數(shù)學算法：分析和合成》，美國ASME，工程學報工業(yè)“：1989 299-306頁。 10、S. H. Lee 和 M. R. Cutkosky. 《具有摩擦性的夾具規(guī)劃》 美國ASME，工業(yè)工程學報：1991，320–327頁。 11、S. Jeng, L. Chen 和W. Chieng.“最小夾緊力分析”，國際機床制造，碩士論文 1995年。 12、E. C. DeMeter.《加工夾具的性能的最小——最大負荷標準》 美國ASME，工業(yè)工程雜志 ：1994 13、E. C. DeMeter .《加工夾具最大負荷的性能優(yōu)化模型》 美國ASME，工業(yè)工程雜志 1995。 14、JH復和AYC倪.“核查和工件夾持的夾具設計”方案優(yōu)化，設計和制造，4，碩士論文： 307-318，1994。 15、T. H. Richards、埃利斯 霍伍德.1977,《應力能量方法分析》，1977。 16、M. J. Hockenberger and E. C. DeMeter. 對工件準靜態(tài)分析功能位移在加工夾具的應用程序，制造科學雜志與工程： 325–331頁, 1996。 Int J Adv Manuf Technol (2001) 17:104113 2001 Springer-Verlag London LimitedFixture Clamping Force Optimisation and its Impact onWorkpiece Location AccuracyB. Li and S. N. MelkoteGeorge W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Georgia, USAWorkpiece motion arising from localised elastic deformationat fixtureworkpiece contacts owing to clamping and machiningforces is known to affect significantly the workpiece locationaccuracy and, hence, the final part quality. This effect can beminimised through fixture design optimisation. The clampingforce is a critical design variable that can be optimised toreduce the workpiece motion. This paper presents a newmethod for determining the optimum clamping forces for amultiple clamp fixture subjected to quasi-static machiningforces. The method uses elastic contact mechanics modelsto represent the fixtureworkpiece contact and involves theformulation and solution of a multi-objective constrainedoptimisation model. The impact of clamping force optimisationon workpiece location accuracy is analysed through examplesinvolving a 32-1 type milling fixture.Keywords: Elasticcontactmodelling;Fixtureclampingforce; Optimisation1.IntroductionThe location and immobilisation of the workpiece are twocritical factors in machining. A machining fixture achievesthese functions by locating the workpiece with respect to asuitable datum, and clamping the workpiece against it. Theclamping force applied must be large enough to restrain theworkpiece motion completely during machining. However,excessive clamping force can induce unacceptable level ofworkpiece elastic distortion, which will adversely affect itslocation and, in turn, the part quality. Hence, it is necessaryto determine the optimum clamping forces that minimise theworkpiece location error due to elastic deformation whilesatisfying the total restraint requirement.Previous researchers in the fixture analysis and synthesisarea have used the finite-element (FE) modelling approach orCorrespondenceandoffprintrequeststo:DrS.N.Melkote,George W. Woodruff School of Mechanical Engineering, GeorgiaInstitute of Technology, Atlanta, Georgia 30332-0405, USA. E-mail:shreyes.melkoteme.gatech.eduthe rigid-body modelling approach. Extensive work based onthe FE approach has been reported 18. With the exceptionof DeMeter 8, a common limitation of this approach is thelarge model size and computation cost. Also, most of the FE-based research has focused on fixture layout optimisation, andclamping force optimisation has not been addressed adequately.Several researchers have addressed fixture clamping forceoptimisation based on the rigid-body model 911. The rigidbody modelling approach treats the fixture-element and work-piece as perfectly rigid solids. DeMeter 12, 13 used screwtheory to solve for the minimum clamping force. The overallproblem was formulated as a linear program whose objectivewas to minimise the normal contact force at each locatingpoint by adjusting the clamping force intensity. The effect ofthe contact friction force was neglected because of its relativelysmall magnitude compared with the normal contact force. Sincethis approach is based on the rigid body assumption, it canuniquely only handle 3D fixturing schemes that involve nomore than 6 unknowns. Fuh and Nee 14 also presentedan iterative search-based method that computes the minimumclamping force by assuming that the friction force directionsare known a priori. The primary limitation of the rigid-bodyanalysis is that it is statically indeterminate when more thansix 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 theelasticity of the fixtureworkpiece system 15. For a relativelyrigid workpiece, the location of the workpiece in the machiningfixture is strongly influenced by the localised elastic defor-mation at the fixturing points. Hockenberger and DeMeter 16used empirical contact force-deformation relations (called meta-functions) to solve for the workpiece rigid-body displacementsdue to clamping and quasi-static machining forces. The sameauthors also investigated the effect of machining fixture designparameters on workpiece displacement 17. Gui et al 18reported an elastic contact model for improving workpiecelocation accuracy through optimisation of the clamping force.However, they did not address methods for calculating thefixtureworkpiece contact stiffness. In addition, the applicationof their algorithm for a sequence of machining loads rep-resenting a finite tool path was not discussed. Li and Melkote19 and Hurtado and Melkote 20 used contact mechanics toFixture Clamping Force Optimisation105solve for the contact forces and workpiece displacement pro-duced by the elastic deformation at the fixturing points owingto clamping loads. They also developed methods for optimisingthe fixture layout 21 and clamping force using this method22. However, clamping force optimisation for a multiclampsystem and its impact on workpiece accuracy were not coveredin these papers.This paper presents a new algorithm based on the contactelasticity method for determining the optimum clamping forcesfor a multiclamp fixtureworkpiece system subjected to quasi-static loads. The method seeks to minimise the impact ofworkpiece motion due to clamping and machining loads onthe part location accuracy by systematically optimising theclamping forces. A contact mechanics model is used to deter-mine a set of contact forces and displacements, which are thenused for the clamping force optimisation. The complete prob-lem is formulated and solved as a multi-objective constrainedoptimisation problem. The impact of clamping force optimis-ation on workpiece location accuracy is analysed via twoexamples involving a 32-1 fixture layout for a milling oper-ation.2.FixtureWorkpiece Contact Modelling2.1Modelling AssumptionsThe machining fixture consists of L locators and C clampswith spherical tips. The workpiece and fixture materials arelinearly 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 forceis assumed to be constant during machining. This assumption isvalid when hydraulic or pneumatic clamps are used.In reality, the elasticity of the fixtureworkpiece contactregion is distributed. However, in this model development,lumped contact stiffness is assumed (see Fig. 1). Therefore, thecontact force and localised deformation at the ith fixturingpoint can be related as follows:Fij= kijdij(1)where kij(j = x,y,z) denotes the contact stiffness in the tangentialand normal directions of the local xi,yi,zicoordinate frame, dijFig. 1. A lumped-spring fixtureworkpiece contact model. xi, yi, zi,denote the local coordinate frame at the ith contact.(j = x,y,z) are the corresponding localised elastic deformationsalong the xi,yi, and ziaxes, respectively, Fij(j = x,j,z) representsthe local contact force components with Fixand Fiybeing thelocal xiand yicomponents of the tangential force, and Fizthenormal force.2.2WorkpieceFixture Contact Stiffness ModelThe lumped compliance at a spherical tip locator/clamp andworkpiece contact is not linear because the contact radiusvaries nonlinearly with the normal force 23. The contactdeformation due to the normal force Piacting between aspherical tipped fixture element of radius Riand a planarworkpiece surface can be obtained from the closed-form Hertz-ian solution to the problem of a sphere indenting an elastichalf-space. For this problem, the normal deformation Dinisgiven as 23, p. 93:Din=S9(Pi)216Ri(E*)2D1/3(2)where1E*=1 n2wEw+1 n2fEfEwand Efare Youngs moduli for the workpiece and fixturematerials, respectively, and nwand nfare Poisson ratios forthe workpiece and fixture materials, respectively.The tangential deformation Dit(= Ditxor Dityin the local xiand yitangential directions, respectively) due to a tangentialforce Qi(= Qixor Qiy) has the following form 23, p. 217:Dtit=Qi8aiS2 nfGf+2 nwGwD(3)whereai=S3PiRi4S1 nfEf+1 nwEwDD1/3and Gwand Gfare shear moduli for the workpiece and fixturematerials, respectively.A reasonable linear approximation of the contact stiffnesscan be obtained from a least-squares fit to Eq. (2). This yieldsthe following linearised contact stiffness values:kiz= 8.82S16Ri(E*)29D1/3(4)kix= kiy=4E*S2 njGf+2 nwGwD1kiz(5)In deriving the above linear approximation, the normal forcePiwas assumed to vary from 0 to 1000 N, and the correspond-ing R2value of the least-squares fit was found to be 0.94.3.Clamping Force OptimisationThe goal is to determine the set of optimal clamping forcesthat will minimise the workpiece rigid-body motion due to106B. Li and S. N. Melkotelocalised elastic deformation induced by the clamping andmachining loads, while maintaining the fixtureworkpiece sys-tem in quasi-static equilibrium during machining. Minimisationof 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.1Objective Function FormulationSince the workpiece rotation due to fixturing forces is oftenquite small 17 the workpiece location error is assumed to bedetermined largely by its rigid-body translation Ddw= DXwDYwDZwT, where DXw, DYw, and DZware the three orthogonalcomponents of Ddwalong the Xg, Yg, and Zgaxes (see Fig. 2).The workpiece location error due to the fixturing forces canthen be calculated in terms of the L2norm of the rigid-bodydisplacement as follows:iDdwi =(DXw)2+ (DYw)2+ (DZw)2)(6)where i i denotes the L2norm of a vector.In particular, the resultant clamping force acting on theworkpiece will adversely affect the location error. When mul-tiple clamping forces are applied to the workpiece, the resultantclamping force, PRC= PRXPRyPRZT, has the form:PRC= RCPC(7)wherePC= PL+1. . .PL+CTistheclampingforcevector,RC= nL+1. . .nL+CTis the clamping force direction matrix,nL+i= cosaL+icosbL+icosgL+iTis the clamping force directioncosine vector, and aL+i, bL+i, and gL+iare angles made by theclamping force vector at the ith clamping point with respectto the Xg, Yg, Zgcoordinate axes (i = 1,2,. . .,C).In this paper, the workpiece location error due to contactregion deformation is assumed to be influenced only by thenormal force acting at the locatorworkpiece contacts. Thefrictional force at the contacts is relatively small and is neg-lected when analysing the impact of the clamping force on theworkpiece location error. Denoting the ratio of the normalcontact stiffness, kiz, to the smallest normal stiffness among alllocators, ksz, by ji(i = 1,. . .,L), and assuming that the workpiecerests on NX, NY, and NZnumber of locators oriented in the Xg,Fig. 2. Workpiece rigid body translation and rotation.Yg, and Zgdirections, the equivalent contact stiffness in theXg, Yg, and Zgdirections can be calculated askszSONXi=1jiD, kszSONYi=1jiD, and kszSONZi=1jiDrespectively (see Fig. 3). The workpiece rigid-body motion,Ddw, due to clamping action can now be written as:Ddw=3PRXkszSONXi=1jiDPRYkszSONYi=1jiDPRZkszSONZi=1jiD4T(8)The workpiece motion, and hence the location error can bereduced by minimising the weighted L2norm of the resultantclamping force vector. Therefore, the first objective functioncan be written as:Minimize iPRCiw=!11PRXONXi=1ji22+1PRYONYi=1ji22+1PRZONZi=1ji222(9)Note that the weighting factors are proportional to the equival-ent contact stiffnesses in the Xg, Yg, and Zgdirections.The components of PRCare uniquely determined by solvingthe contact elasticity problem using the principle of minimumtotal complementary energy 15, 23. This ensures that theclamping 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 staticequilibrium by the clamping forces at all times. Therefore, theminimisation of the total complementary energy forms thesecond objective function for the clamping force optimisationand is given by:Minimise (U* W*) =12FOL+Ci=1(Fix)2kix+OL+Ci=1(Fiy)2kiy+OL+Ci=1(Fiz)2kizG(10)= .lTQlFig. 3. The basis for the determination of the weighting factor for theL2norm calculation.Fixture Clamping Force Optimisation107where U* represents the complementary strain energy of theelastically deformed bodies, W* represents the complementarywork done by the external force and moments, Q = diagc1xc1yc1z. . . cL+CxcL+CycL+Cz is the diagonal contact compliancematrix, cij= (kij)1, and l = F1xF1yF1z. . . FL+CxFL+CyFL+CzTis thevector of all contact forces.3.2Friction and Static Equilibrium ConstraintsThe optimisation objective in Eq. (10) is subject to certainconstraints and bounds. Foremost among them is the staticfriction constraint at each contact. Coulombs friction law statesthat(Fix)2+ (Fiy)2) # misFiz(misis the static friction coefficient).A conservative and linearised version of this nonlinear con-straint can be used and is given by 19:uFixu + uFiyu # misFiz(11)Since quasi-static loads are assumed, the static equilibriumof the workpiece is ensured by including the following forceand moment equilibrium equations (in vector form):OF = 0(12)OM = 0where the forces and moments consist of the machining forces,workpiece weight and the contact forces in the normal andtangential directions.3.3BoundsSince the fixtureworkpiece contact is strictly unilateral, thenormal contact force, Pi, can only be compressive. This isexpressed by the following bound on Pi:Pi$ 0(i = 1, . . ., L + C)(13)where it is assumed that normal forces directed into theworkpiece are positive.In addition, the normal compressive stress at a contact cannotexceed the compressive yield strength (Sy) of the workpiecematerial. This upper bound is written as:Pi# SyAi(i = 1, . . .,L+C)(14)where Aiis the contact area at the ith workpiecefixture con-tact.The complete clamping force optimisation model can nowbe written as:Minimize f =Hf1f2J=H.lTQliPRCiwJ(15)subject to: (11)(14).4.Algorithm for Model SolutionThe multi-objective optimisation problem in Eq. (15) can besolved by the e-constraint method 24. This method identifiesone of the objective functions as primary, and converts theother into a constraint. In this work, the minimisation of thecomplementary energy (f1) is treated as the primary objectivefunction, and the weighted L2norm of the resultant clampingforce (f2) is treated as a constraint. The choice of f1as theprimary objective ensures that a unique set of feasible clampingforces is selected. As a result, the workpiecefixture system isdriven to a stable state (i.e. the minimum energy state) thatalso has the smallest weighted L2norm for the resultantclamping force.The conversion of f2into a constraint involves specifyingthe weighted L2norm to be less than or equal to e, where eis an upper bound on f2. To determine a suitable e, it isinitially assumed that all clamping forces are unknown. Thecontact forces at the locating and clamping points are computedby considering only the first objective function (i.e. f1). Whilethis set of contact forces does not necessarily yield the lowestclamping forces, it is a “true” feasible solution to the contactelasticity problem that can completely restrain the workpiecein the fixture. The weighted L2norm of these clamping forcesis computed and taken as the initial value of e. Therefore,the clamping force optimisation problem in Eq. (15) can berewritten as:Minimize f1= .lTQl(16)subject to: iPRCiw$ e, (11)(14).An algorithm similar to the bisection method for findingroots of an equation is used to determine the lowest upperbound for iPRCiw. By decreasing the upper bound e as muchas possible, the minimum weighted L2norm of the resultantclamping force is obtained. The number of iterations, K, neededto terminate the search depends on the required predictionaccuracy d and ueu, and is given by 25:K =Flog2SueudDG(17)where I denotes the ceiling function. The complete algorithmis given in Fig. 4.5.Determination of Optimum ClampingForces During MachiningThe algorithm presented in the previous section can be usedto determine the optimum clamping force for a single loadvector applied to the workpiece. However, during millingthe magnitude and point of cutting force application changescontinuously along the tool path. Therefore, an infinite set ofoptimum clamping forces corresponding to the infinite set ofmachining loads will be obtained with the algorithm of Fig. 4.This substantially increases the computational burden and callsfor a criterion/procedure for selecting a single set of clampingforces that will be satisfactory and optimum for the entire toolpath. A conservative approach to addressing these issues isdiscussed next.Consider a finite number (say m) of sample points alongthe tool path yielding m corresponding sets of optimum clamp-ing forces denoted as P1opt, P2opt, . . ., Pmopt. At each sampling108B. Li and S. N. MelkoteFig. 4. Clamping force optimisation algorithm (used in example 1).point, the following four worst-case machining load vectorsare considered:FXmax= FmaxXF1YF1ZTFYmax= F2XFmaxYF2ZTFZmax= F3XF3YFmaxZT(18)Frmax= F4XF4YF4ZTwhere FmaxX, FmaxY, and FmaxZare the maximum Xg, Yg, and Zgcomponents of the machining force, the superscripts 1, 2, 3 ofFX, FY, and FZstand for the other two orthogonal machiningforcecomponentscorrespondingtoFmaxX, FmaxY, and FmaxZ,respectively, and iFrmaxi = max(FX)2+ (FY)2+ (FZ)2).Although the four worst-case machining load vectors willnot act on the workpiece at the same instant, they will occuronce per cutter revolution. At conventional feedrates, the errorintroduced by applying the load vectors at the same pointwould be negligible. Therefore, in this work, the four loadvectorsareappliedatthesamelocation(butnotsimultaneously) on the workpiece corresponding to the sam-pling instant.The clamping force optimisation algorithm of Fig. 4 is thenused to calculate the optimum clamping forces correspondingto each sampling point. The optimum clamping forces havethe form:Pijmax= Ci1jCi2j. . . CiCjT(i = 1, . . .,m)(j = x,y,z,r)(19)where Pijmaxis the vector of optimum clamping forces for thefour worst-case machining load vectors, and Cikj(k = 1,. . .,C)is the force magnitude at each clamp corresponding to the ithsample point and the jth load scenario.After Pijmaxis computed for each load application point, asingle set of “optimum” clamping forces must be selected fromall of the optimum clamping forces found for each clamp fromall the sample points and loading conditions. This is done bysorting the optimum clamping force magnitudes at a clampingpoint for all load scenarios and sample points and selectingthe maximum value, Cmaxk, as given in Eq. (20):Cmaxk# Cikj(k = 1,. . .,C)(20)Once this is complete, a set of optimised clamping forcesPopt= Cmax1Cmax2. . . CmaxCTis obtained. These forces must beverified for their ability to ensure static equilibrium of theworkpiecefixture system. Otherwise, more sampling points areselected and the aforementioned procedure repeated. In thisfashion, the “optimum” clamping force, Popt, can be determinedfor the entire tool path. Figure 5 summarises the algorithm justdescribed. Note that although this approach is conservative, itprovides a systematic way of determining a set of clampingforces that minimise the workpiece location error.6.Impact on Workpiece LocationAccuracyIt is of interest to evaluate the impact of the clamping forcealgorithm presented earlier on the workpiece location accuracy.The workpiece is first placed on the fixture baseplate in contactwith the locators. Clamping forces are then applied to pushthe workpiece against the locators. Consequently, localiseddeformations occur at each workpiecefixture contact, causingthe workpiece to translate and rotate in the fixture. Sub-sequently, the quasi-static machining load is applied causingadditional motion of the wor