喜歡這套資料就充值下載吧。。。資源目錄里展示的都可在線預(yù)覽哦。。。下載后都有,,請放心下載,,文件全都包含在內(nèi),,【有疑問咨詢QQ:414951605 或 1304139763】
========================================喜歡這套資料就充值下載吧。。。資源目錄里展示的都可在線預(yù)覽哦。。。下載后都有,,請放心下載,,文件全都包含在內(nèi),,【有疑問咨詢QQ:414951605 或 1304139763】
========================================
黃河科技學(xué)院本科畢業(yè)設(shè)計(論文)任務(wù)書
工 學(xué)院 機 械 系 機械設(shè)計制造及其自動化 專業(yè) 08 級 1 班
學(xué)號 080105007 學(xué)生 劉 文 超 指導(dǎo)教師 王 良 文
畢業(yè)設(shè)計(論文)題目 GQ50型鋼筋切斷機的結(jié)構(gòu)設(shè)計與運動仿真
畢業(yè)設(shè)計(論文)工作內(nèi)容與基本要求(目標(biāo)、任務(wù)、途徑、方法,應(yīng)掌握的原始資料(數(shù)據(jù))、參考資料(文獻(xiàn))以及設(shè)計技術(shù)要求、注意事項等)
一、設(shè)計技術(shù)要求、原始資料(數(shù)據(jù))、參考資料(文獻(xiàn))
1、設(shè)計要求:
1)建立運動分析模型 ,分析曲柄滑塊機構(gòu)的位移、速度、加速度。
2)完成50型鋼筋切斷機的結(jié)構(gòu)設(shè)計。
3)完成傳動系統(tǒng)的計算機的仿真分析(采用三維軟件建模并仿真)。
2、原始資料(數(shù)據(jù)):
切斷機的電機,4kw;切斷次數(shù),28(次/分);切斷行程:40mm。
3、參考資料
1)筑路機械工程著作;
2)建筑機械雜志。
二、設(shè)計目標(biāo)與任務(wù)
設(shè)計出滿足要求的鋼筋切斷機,并完成其裝配圖與部分零件圖,查閱文獻(xiàn)資料不少于12篇,其中外文資料不少于2篇。
1、文獻(xiàn)綜述一篇,不少于3000字,與專業(yè)相關(guān)的英文翻譯一篇,不少于3000漢字。
2、畢業(yè)設(shè)計說明書一份,內(nèi)容與字?jǐn)?shù)都不少于規(guī)定的任務(wù)量。
3、圖紙若干(折合后不少于A1圖紙3張,可以用計算機繪圖)。
4、包含本次設(shè)計的所有內(nèi)容的光盤一張。
畢業(yè)設(shè)計(論文)撰寫規(guī)范及有關(guān)要求,請查閱《黃河科技學(xué)院本科畢業(yè)設(shè)計(論文)指導(dǎo)手冊》。
三、時間安排
1-4 周 完成開題報告、文獻(xiàn)翻譯、文獻(xiàn)綜述及總體方案設(shè)計。
5-11 周 完成總體設(shè)計、完成部分機構(gòu)的裝配圖及部分零件圖并撰寫說明書。
12-13 周 修改論文、資格審查等。
14 周 畢業(yè)答辯。
畢業(yè)設(shè)計(論文)時間: 年 月 日至 年 月 日
計 劃 答 辯 時 間: 年 月 日
專業(yè)(教研室)審批意見:
審批人簽名:
課題名稱
GQ50型鋼筋切斷機的結(jié)構(gòu)設(shè)計與運動仿真
課題來源
教師擬訂
課題類型
AX
指導(dǎo)教師
王良文
學(xué)生姓名
劉文超
專 業(yè)
機械設(shè)計制造及其自動化
學(xué) 號
080105007
一、調(diào)研資料的準(zhǔn)備
1. 參觀黃河金剛石有限公司的GQ40型和GQ50型鋼筋切斷機,了解鋼筋切斷機的工作原理。
2. 利用圖書館及網(wǎng)絡(luò)資源,查閱相關(guān)資料。
二、設(shè)計的目的與要求
1. 建立運動分析模型。
2. 建立機械的動力學(xué)模型。
3. 設(shè)計整體結(jié)構(gòu)。
三、思路與預(yù)期成果
1、設(shè)計思路:
(1)首先查詢大量的相關(guān)資料,并熟悉AUTO CAD及Inventor2011繪圖軟件,達(dá)到能夠較為熟練的使用,以提高設(shè)計及繪圖效率。
(2)根據(jù)畢業(yè)時間的安排,逐步推進(jìn)設(shè)計進(jìn)程,有不懂之處就盡快與老師進(jìn)行溝通,得到老師指導(dǎo)并解決問題。
2、預(yù)期的成果
(1)完成文獻(xiàn)綜述一篇,不少與3000字,與專業(yè)相關(guān)的英文翻譯一篇,不少于3000字。
(2)完成內(nèi)容與字?jǐn)?shù)都不少于規(guī)定量的畢業(yè)設(shè)計說明書一份,繪制出各主要零件的工程圖。
(3)建立GQ50型鋼筋切斷機的三維模型,并進(jìn)行運動仿真。
四、任務(wù)完成的階段內(nèi)容及時間安排
1-4 周 完成開題報告、文獻(xiàn)翻譯及文獻(xiàn)綜述。
5-11 周 完成總體設(shè)計,基本完成機構(gòu)的裝配圖及零件圖,并撰寫說明書。
12-13 周 修正并完善論文終稿,進(jìn)行資格審查。
14 周 畢業(yè)答辯。
五、完成設(shè)計(論文)所具備的條件因素
修完機械設(shè)計,工程材料等課程,借助圖書館的相關(guān)文獻(xiàn)資料,以及網(wǎng)絡(luò)等資源。
指導(dǎo)教師簽名: 日期:
黃河科技學(xué)院畢業(yè)設(shè)計(論文)開題報告表
課題來源:(1)教師擬訂;(2)學(xué)生建議;(3)企業(yè)和社會征集;(4)科研單位提供
課題類型:(1)A—工程設(shè)計(藝術(shù)設(shè)計);B—技術(shù)開發(fā);C—軟件工程;D—理論研究;E—調(diào)查報告
(2)X—真實課題;Y—模擬課題;Z—虛擬課題
要求(1)、(2)均要填,如AY、BX等
目 錄
1.任務(wù)書………………………………………………………1
2.開題報告……………………………………………………2
3.指導(dǎo)教師評閱表……………………………………………3
4.主審教師評審表……………………………………………4
5.畢業(yè)設(shè)計答辯評審與總成績評定表………………………5
6.畢業(yè)設(shè)計說明書……………………………………………6
7.文獻(xiàn)綜述……………………………………………………58
8.文獻(xiàn)翻譯……………………………………………………64
9.光盤
10.設(shè)計圖紙
單位代碼 02
學(xué) 號 080105007
分 類 號 TH
密 級
畢業(yè)設(shè)計
文獻(xiàn)綜述
院(系)名稱
工學(xué)院機械系
專業(yè)名稱
機械設(shè)計制造及其自動化
學(xué)生姓名
劉文超
指導(dǎo)教師
王良文
2012年 03 月 20
黃河科技學(xué)院畢業(yè)設(shè)計(文獻(xiàn)綜述) 第 8 頁
鋼筋切斷機運動仿真綜述
摘要
本文調(diào)研了最近幾年國內(nèi)外有關(guān)鋼筋切斷機的專利、論文,以及現(xiàn)在的研發(fā)現(xiàn)狀。分析研究鋼筋切斷機的動力學(xué)性能是進(jìn)行鋼筋切斷機設(shè)計、改良的理論基礎(chǔ)。為了推動鋼筋切斷機設(shè)計創(chuàng)新,縮短鋼筋切斷機的設(shè)計周期,在產(chǎn)品試制及完成測試之前明確了解鋼筋切斷機的動力學(xué)參數(shù),及時發(fā)現(xiàn)并消除設(shè)計過程中存在的缺陷,最大程度地節(jié)省資金和時間。采用計算機仿真技術(shù)結(jié)合動力學(xué)理論,對鋼筋切斷機進(jìn)行零件三維設(shè)計、裝配,建立鋼筋切斷機三維仿真分析模型,模擬鋼筋切斷機運行狀況,對鋼筋切斷機進(jìn)行仿真分析研究, 以加快產(chǎn)品技術(shù)更新。
關(guān)鍵詞: 鋼筋切斷機,動力學(xué)模型,結(jié)構(gòu)設(shè)計,運動仿真
一.引言
鋼筋切斷機是一種各類工程建設(shè)領(lǐng)域廣泛使用的設(shè)備,在近年來的生產(chǎn)、使用中呈現(xiàn)快速增長的趨勢。國內(nèi)各相關(guān)廠家先后開發(fā)了GQ32型、GQ40型、GQ50型、GQ60型、GQ75型等多個品種的切斷機。隨著建設(shè)工程規(guī)模的擴大,建筑質(zhì)量的要求提高,鋼筋切斷機的規(guī)格有逐漸增大的趨勢。各生產(chǎn)廠家為提高自己的生產(chǎn)產(chǎn)品的技術(shù)性能、增加產(chǎn)品的競爭優(yōu)勢,不斷優(yōu)化結(jié)構(gòu),努力降低生產(chǎn)成本,并不斷開發(fā)新的品種。但由于未見國內(nèi)外有關(guān)鋼筋切斷機動力學(xué)研究的文獻(xiàn),因而在改良產(chǎn)品設(shè)計及開發(fā)新的產(chǎn)品過程中,在產(chǎn)品試制及完成測試之前尚不能對結(jié)構(gòu)的動力學(xué)參數(shù)及驅(qū)動電機在工作過程中的真實運動情況有比較明確的了解,而電機在工作過程中的速度降是衡量鋼筋切斷機設(shè)計質(zhì)量的重要指標(biāo)之一。為此,我們首先建立了一種鋼筋切斷機的通用動力學(xué)計算模型,使我們在產(chǎn)品的改良設(shè)計及研發(fā)新產(chǎn)品的過程中通過AutoCAD完成產(chǎn)品結(jié)構(gòu)圖的設(shè)計后,通過聯(lián)機計算的方法就可以清楚的求得所開發(fā)產(chǎn)品的動力學(xué)參數(shù)及驅(qū)動電機的真實運動情況,可以極大的加快產(chǎn)品的開發(fā)速度及方便的修改產(chǎn)品的結(jié)構(gòu)參數(shù)。相關(guān)工作對于提高國內(nèi)產(chǎn)品技術(shù)開發(fā)工作的進(jìn)度,提高產(chǎn)品的競爭力有一定的積極作用。[6][7][8]
二.鋼筋切斷機簡介
鋼筋切斷機是一種在工程及建筑領(lǐng)域中廣泛使用的設(shè)備,由于近年來房地產(chǎn)的迅速發(fā)展,鋼筋切斷機的生產(chǎn)、使用呈現(xiàn)快速增長的趨勢,隨著建筑工程規(guī)模的擴大和對建筑質(zhì)量要求的提高,對鋼筋切斷機的性能也提出了更高的要求。[9]
目前,國產(chǎn)鋼筋切斷機的典型結(jié)構(gòu)下圖所示,電機通過一級帶傳動、三級齒輪傳動來驅(qū)動曲柄滑塊機構(gòu)帶動活動刀做往復(fù)運動,固定在活動刀座上的活動刀片和固定在機體上的固定刀片一起作用完成對鋼筋的切斷。[10]
三.相關(guān)專利產(chǎn)品介紹(專利[1][2][3])
專利[1]的產(chǎn)品是一種以柴油機為動力的鋼筋切斷機,其特征在于將鋼筋切斷機與柴油機相組合、中間用皮帶傳動相聯(lián)接構(gòu)成一柴油鋼筋切斷機。其優(yōu)點和積極效果是:不需要電力,使沒有電力或者電力不足的地方,也夠使用鋼筋切斷機完成鋼筋的切斷。
右圖所示的是專利[1]產(chǎn)品的結(jié)構(gòu)示意圖:1為鋼筋切斷機機構(gòu),2為底架、3為皮帶傳動、4為柴油機、5 和6分別是柴油機4和切斷機1與底架2相連接的連接螺栓。柴油機4和鋼筋切斷機構(gòu)1都被連接螺栓5和6固定在底架2上,發(fā)動柴油機4通過皮帶傳動3帶動鋼筋切斷機構(gòu)1的刀具來完成切斷鋼筋的工作要求。
專利[2]的產(chǎn)品是一種改進(jìn)的鋼筋切斷機。它的箱體由側(cè)面箱板和箱底板組裝成封閉型式:其傳動系統(tǒng)采用三軸(六齒輪)三級齒輪傳動結(jié)構(gòu);而其曲柄桿與滑塊采用鉤型連接結(jié)構(gòu)。這種改進(jìn),使箱體制造工藝簡化。由于是三級齒輪傳動節(jié)省了一根中間傳動軸,故而使整機體積縮小了許多,且降低了成本。而曲柄桿與滑塊的鉤型連接方式,不僅結(jié)構(gòu)簡單,而且提高了機械強度,不易損壞,使整機延長了壽命。
專利[3]的產(chǎn)品是一種改進(jìn)的鋼筋切斷機。包括電機、飛輪、第一級至第四級齒輪傳動裝置、離合器、定動切刀和哈夫式齒輪箱體,第一級齒輪傳動裝置的傳動軸兩端分別套裝有飛輪;哈夫式齒輪箱體以曲軸中心水平線為依據(jù),由鑄成上、下各一半的兩塊箱體組合而成。本實用新型由于在現(xiàn)有四級齒輪傳動鋼筋切斷機的基礎(chǔ)上,采用哈夫式箱體并在第一級齒輪傳動裝置的傳動軸兩端分別套裝有飛輪,故其故障率低,使用壽命長,阻力大;且安裝、維修方便該專利產(chǎn)品是對國內(nèi)外產(chǎn)品的一種改進(jìn),是一種實用新型的鋼筋切斷機,該產(chǎn)品主要包括電動機,飛輪,第一級至第四級齒輪傳動裝置、離合器、定動刀片和哈夫式齒輪箱體,起主要特別之處在于:第一級齒輪傳動裝置的傳動軸兩端分別裝有飛輪。具體結(jié)構(gòu)如上圖。
四. 鋼筋切斷機幾個主要部件
(1) 刀口裝置(專利[4])
專利[4]指推動刀口軸桿的偏心軸外圍裝設(shè)有方形滑塊,于偏心軸轉(zhuǎn)動時能使方形滑塊上下移動,進(jìn)而使刀片軸桿以平面方式向外推動或拉回,使鋼筋的切口平直,并于切斷鋼筋后,偏心軸拉回軸桿撞擊偏心軸的力,形成面的使受力點平均,降低偏心軸的耗損,且該偏心軸與推動齒輪間設(shè)有一固定隔板,以避免撞擊的同時造成齒輪的偏移而損壞,該刀口軸桿的另一側(cè)設(shè)有一活動隔板,能增進(jìn)維修時的便利性。
(2) 離合器組件(專利[5])
專利[5]鋼筋切斷機的曲軸離合器。該器包括曲軸、內(nèi)齒輪、大齒輪和安裝在兩輪內(nèi)接之間的活動銷齒機構(gòu)。該機構(gòu)又由轉(zhuǎn)銷軸、轉(zhuǎn)銷齒、撥動擋板、拉簧和伸縮擋銷組成。通過操縱轉(zhuǎn)銷齒接觸于大齒輪的轉(zhuǎn)銷齒推槽內(nèi)而使大齒輪帶動內(nèi)齒輪和曲軸轉(zhuǎn)動。或操縱轉(zhuǎn)銷齒接觸于內(nèi)齒輪的轉(zhuǎn)銷軸槽內(nèi)而使內(nèi)齒輪脫離大齒輪,使大齒輪空轉(zhuǎn)。該離合器可適用于切斷機、沖床、剪切機床上,其具有離合平穩(wěn)、體積小、結(jié)構(gòu)緊湊的優(yōu)點。
五. 鋼筋切斷機存在的缺點和問題
目前,國內(nèi)外現(xiàn)有的鋼筋切斷機,就其傳動方式分為兩種類型:一種是開式傳動,一種是閉式傳動。[11]開式傳動的形式,其軸、齒輪等主要傳動件都要暴露在外,各傳動件均要靠人工加油,極為不方便。閉式傳動的形式,主要傳動件均安裝在機器箱體內(nèi),形成整體密封,其機體為整體構(gòu)件,加工工藝復(fù)雜。國內(nèi)外的鋼筋切斷機主要存在以下幾個問題:
(1) 用電動機作原動機的鋼筋切斷機和用柴油機作原動機的鋼筋切斷機,所用的原動機不能代換,這給使用帶來了一定的不方便。
(2) 國內(nèi)外使用的全封閉結(jié)構(gòu)的鋼筋切斷機的活動刀處采用上開蓋的機體,大都是采用的是減速器結(jié)構(gòu),用螺栓連接,而這種結(jié)構(gòu)的剛性較差,而且在結(jié)合面處加工困難,并且常常漏油。而機體采用整體式結(jié)構(gòu)的鋼筋切斷機,活動刀處采用的是側(cè)開蓋結(jié)構(gòu),由于切斷機工作時所受的側(cè)向力很大,常常損壞。
(3) 國內(nèi)外所使用的鋼筋切斷機在空套齒輪與軸結(jié)合處,常常容易損壞。
(4) 國內(nèi)外所使用的鋼筋切斷機,其離合器操縱機構(gòu)只有一種方式,即腳踏式,這不適應(yīng)于不同習(xí)慣的操作者和不利于緊急情況的操作。
(5) 國內(nèi)外所使用的鋼筋切斷機,大都采用的是小功率的電動機,工作主軸的偏心距也很小,這樣在剪切公稱直徑范圍內(nèi)的不同直徑的鋼筋時,就要更換刀片或者在刀片的后面加減刀墊,才能完成剪切不同尺寸鋼筋的切斷要求。
(6) 國內(nèi)外所使用的鋼筋切斷機,其固定刀片側(cè)隙的調(diào)整是靠人工選配尺寸固定的刀墊來實現(xiàn)的。
六.鋼筋切斷機的改良方向
真對鋼筋切斷機存在的這些問題,提出以下改良方向:
(1)結(jié)構(gòu)與造型的改良。
目前的鋼筋切斷機有閉式與開式兩大類。開式由于體積大,搬移不便、潤滑差等,已經(jīng)很少生產(chǎn)使用。而閉式結(jié)構(gòu)的體積雖較為緊湊,但存在的主要問題是在傳動系統(tǒng)出現(xiàn)故障時,維修不太方便。目前在市場上出現(xiàn)了一種在整體閉式結(jié)構(gòu)上的改良設(shè)計。[12][13]即將機體一側(cè)的1/2~1/3設(shè)計為可拆卸結(jié)構(gòu),這種結(jié)構(gòu)雖然在加工上增加了工序,但鑄造環(huán)節(jié)更簡單,安裝比較方便,尤其是便于售后服務(wù)。
(2)輕量化設(shè)計
輕量化設(shè)計可節(jié)省資源、減少生產(chǎn)、使用、回收等環(huán)節(jié)中了浪費,而且就機器本身而言,也可以通過輕量化設(shè)計,在最低成本下達(dá)到機器的使用性能。就目前來說,可通過建立虛擬樣機、并對其進(jìn)行有限元分析實現(xiàn)輕量化設(shè)計及制造。[14][15]
七. 其他
機械動力學(xué)的相關(guān)內(nèi)容:
機械動力學(xué)主要研究的是機械的振動和平衡的問題。目前,我們正在高速發(fā)展的階段,各個部門迫切需要大量的新的高效率,高速度,高精度和高自動化的機械設(shè)備。隨著機械速度的提高,機械平衡和震動的問題成為機械設(shè)計中的一個問題,在設(shè)計高精度的機器時,就必然涉及各種動力學(xué)的因素,需要精確計算各部件的真實運動情況以及考慮部件的彈性、運動副中的磨察等因素對構(gòu)件運動的影響,才能使各部件的動作協(xié)調(diào),機械運轉(zhuǎn)正常。[17][18][19][20]
八. 結(jié)束語
產(chǎn)品的設(shè)計開發(fā)過程是一個不斷改良完善、技術(shù)優(yōu)化發(fā)展的過程。唯有通過長期的市場考驗,且有優(yōu)良的性價比,才能真正占領(lǐng)市場。
鋼筋切斷機的運動仿真的設(shè)計必將更加有利于鋼筋切斷機產(chǎn)品的開發(fā),可以加快產(chǎn)品的開發(fā)速度以及開發(fā)出更加高效的產(chǎn)品。
參考文獻(xiàn)
專利[1] 孟憲民,劉鴻鷹,章友文,郭俊平,付竹祥 柴油鋼筋切斷機
專利號:91209535.0
專利[2] 陶昆 , 鋼筋切斷機 專利號:94217280.9
專利[3] 蔣家善 改進(jìn)的鋼筋切斷機 專利號:00240526.1
專利[4] 邱文聰 鋼筋切斷機刀口裝置 專利號:96205392.9
專利[5] 全承浩 鋼筋切斷機的曲軸離合器 專利號:95214902.8
[6] 胡永華;劉中 《鋼筋切斷機參數(shù)化設(shè)計系統(tǒng)開發(fā)》 2011年4月
[7] 陶浩;段紅杰 《鋼筋切斷機的動力學(xué)性能》 2008年 5月
[8] 段紅杰;陶浩 《鋼筋切斷機動力學(xué)仿真設(shè)計》 2010年4月
[9] 王文良,李長詩 《鋼筋切斷機通用動力學(xué)》 2008年8月
[10] 張金芳 袁俊華 陳 偉 《人機界面型全自動鋼筋切斷機》 2010年8月
[11] 劉瑞東 《液壓隨動切斷、液壓馬達(dá)送料鋼筋定長切斷機》 2010年10月
[12] 王良文 《國產(chǎn)鋼筋切斷機的生產(chǎn)現(xiàn)狀與改良方向》 2009年1月
[13] 樊智慧 《鋼筋切斷機曲軸連桿機構(gòu)的改進(jìn)》 1995年3月
[14] 孟慶華 《鋼筋切斷機液壓系統(tǒng)的設(shè)計研究》 2010年8月
[15] 賀生華 《普通鋼筋切斷機技術(shù)改造》 2008年一月
[16] 車仁煒 陸念力 王樹春 《一種新型鋼筋切斷機的設(shè)計研究》 2004年
[17] 唐錫寬《機械動力學(xué)》 高等教育出版社 1984.5
[18] 張策 《機械動力學(xué)》 北京:高等教育出版社 2001年4月
[19] Bevefidge R,Bolme D,Teixeira M,et a1.The CSU Face Iden—tification Evaluation System[R].Users’Guide Version 5.0,technical report,Computer Science Dept.,Colorado State Univ.,May 2003.
[20] Huang L—L,Shimizu A,and Kobatake H.Classification—based face detection using gabor filter features[A].The 6th IEEE Internati0nal Conference on Automatic Face and Gesture Recognition(FG2004)[C].2OO4:397-402.
單位代碼 02
學(xué) 號 080105007
分 類 號 TH
密 級
畢業(yè)設(shè)計
文獻(xiàn)翻譯
院(系)名稱
工學(xué)院機械系
專業(yè)名稱
機械設(shè)計制造及其自動化
學(xué)生姓名
劉文超
指導(dǎo)教師
王良文
2012 年 3月 20日
黃河科技學(xué)院畢業(yè)設(shè)計(文獻(xiàn)翻譯) 第 11 頁
利用離線仿真結(jié)果和3D模型變形以實現(xiàn)實時遙控結(jié)構(gòu)變形的方法
摘要
DTP2,一個為了演示和改進(jìn)遠(yuǎn)程操作設(shè)備ITER全面的物理測試設(shè)備,已經(jīng)在芬蘭建立起來。首個裝備有SCEE的RH設(shè)備原型CMM已經(jīng)于2008年10月移交給DTP2。其目的是為了驗證CMM/SCEE原型可以被成功的應(yīng)用于第二個暗盒的RH的運作。在 F4E 授與 " DTP2 測試設(shè)備運行和升級準(zhǔn)備 " 結(jié)束的時候,第二個暗盒的 RH 的運行成功地為F4E代表做了證明。
得益于CMM/SCEE機器人的設(shè)計,所以當(dāng)它在3.6米長的控制桿上運載9噸重的第二暗盒時,具有相當(dāng)大的機械彈性。這也就導(dǎo)致數(shù)據(jù)不精確,并且用于控制系統(tǒng)的3D模型也不能準(zhǔn)確的反映CMM/SCEE機器人的變化狀態(tài)。 為了提高其精確度,已經(jīng)發(fā)展出了一種在虛擬環(huán)境中控制其彈性的方法。加載在CMM/SCEE上載荷的作用大小被測量并且最小化到由控制系統(tǒng)軟件執(zhí)行的載荷補償模型上。這種優(yōu)化的方法利用有限元分析,通過3D模型的變形解釋了控制系統(tǒng)機器的結(jié)果變形。 這將促使CMM/SCEE的絕對精度和3D模型的適合性有一個相當(dāng)大的改進(jìn),這對RH應(yīng)用程序是至關(guān)重要的,因為控制裝置的視覺信息是受周圍環(huán)境的限制的。
關(guān)鍵詞:國際熱核實驗反應(yīng)堆 遠(yuǎn)程控制 偏濾器測試平臺2 虛擬工程 彈性變形 虛擬現(xiàn)實
1.引言
這篇論文展示了系統(tǒng)控制軟件執(zhí)行的負(fù)載補償功能是怎樣改進(jìn)DTP2控制系統(tǒng)的絕對精度和可視化精度的。同時也找到了一種通過3D模型變形來解釋DTP2的結(jié)構(gòu)變形的新方法,利用有限元分析來導(dǎo)出變化范圍。除此之外,真實的組件變形,2D結(jié)構(gòu)變形會被用于顯示每單位結(jié)構(gòu)負(fù)荷的運行評估。
在第二暗盒安裝程序的時候,CMM在電機傳動裝置的協(xié)作下,沿著射線方向行進(jìn)到維持隧道的頂端。(圖1)。在垂直面上的上升和傾斜運動可以用來憑借向上維持隧道來控制暗盒的方位。當(dāng)熱運動到達(dá)第二暗盒時,由CRO和HRO回轉(zhuǎn)連接的可以用來改變暗盒的方位。
Fig. 1. CMM and SCEE structural representation.
2. DTP2的偏差研究目錄
2.1 CMM/SCEE檢驗
在CMM/SCEE傳遞到DTP2后,系統(tǒng)綜合階段開始啟動,以為實際測試做系統(tǒng)準(zhǔn)備。這個測試在工廠啟動實施,在RH控制室中結(jié)束測試運行。[1]
最初用于暗盒運行的的程序是被用來教學(xué)的,這與暗盒有持續(xù)性的視覺聯(lián)系。CMM/SCEE良好的重復(fù)精度(3mm)保證了運行程序成功重復(fù)。然而,由于CMM/SCEE的完全精度不夠準(zhǔn)確,導(dǎo)致靜態(tài)的三維模型不能很好的支持運行。三維模型有時候會出現(xiàn)暗盒與DRM沖突的情況,但實際情況是一切運行良好。很明顯,在遠(yuǎn)程操作之前,系統(tǒng)的絕對精度需要改進(jìn)。
2.2 載荷補償
在運行程序的時候,載荷對CMM-SCEE運動鏈的影響會被測量。并且知道位于暗盒尖端的定位誤差最大接近80mm。這些測量數(shù)據(jù)被用來生成載荷補償以改進(jìn)絕對精度。這種解決方法對于RH維持通道的運行是十分普遍的,但是對于負(fù)載補償模型,確實一個基于CRO蓮價值的平臺。由于CMM在將來普遍支持其他終端執(zhí)行器,所以這種方法簡單,易用。具體的查表只應(yīng)用于在特殊的環(huán)形SCEE軌道中的操作。補償功能還可以明顯改善設(shè)備性能。因此,暗盒尖端的最大誤差由80mm下降到了5mm。[1]
載荷補償?shù)膶嵤﹥r值可以參考圖2中的笛卡爾坐標(biāo)系。根據(jù)暗盒是否加載到HRO上,解決的方法也分為兩個階段?!袄硐朐O(shè)備的笛卡爾參考系”(圖2)表達(dá)了與軸相連的HRO鏈的位置坐標(biāo)。因此,HRO鏈僅僅被用在改變縱軸周圍暗盒的方向,并且,之后CRO鏈可以應(yīng)用于y軸參考數(shù)據(jù)的評估。因此,在熱運動時,負(fù)載補償?shù)墓δ芤蕾囉贑RO鏈。
如果在HRO鏈上沒有負(fù)載,一個逆運動學(xué)解可以直接用于解決聯(lián)合相應(yīng)的數(shù)值參考。解決方法是使用包括基于簽名修正的CMM / SCEE校準(zhǔn)的Denavit-Hartenberg參數(shù)計算。
Fig. 2. Left: load compensation in cartesian space. Right: implementation of load effect to the joint data of the real device.
當(dāng)載荷是連接到HRO關(guān)節(jié),在這種情況下,由于笛卡爾參考也會受CMM/SCEE的撓度影響,所以機器人的逆運動學(xué)解并不可直接使用。當(dāng)產(chǎn)生的作用是已知的,正確的評估CRO鏈的價值可以迭代利用負(fù)載補償或者定義了并列價值與CRO鏈價值之間的適應(yīng)性。迭代解和7th多項式都能很好的應(yīng)用于實踐中。CRO價值鏈被定義后,在x,y,z方向的位置補償和圓周與徑向的定向補償可以做成笛卡爾參考系。由于CMM/SCEE缺乏在yaw方向上移動的能力,所以無法做運動補償。
Fig. 3. Ansys FEM result (DCM lifted from RH interface).
Fig. 4. CATIA FEM result (DCM lifted from RH interface).
2.3 改進(jìn)遙控裝置的可視化精度
當(dāng)增加一個鏈接到CMM/SCEE的三維模型上時,暗盒在yaw方向的傾斜是可視化的。這個鏈接已經(jīng)被放置在勾板和暗盒之間。因此,操作者可以看見暗盒傾斜的作用 ,它在垂直方向上最大有效運動距離是10mm。
為了增加可視化精度,當(dāng)暗盒連接DRM通道內(nèi)部與外部的時候,壓力差超過上升油缸提供的載荷,暗盒的重量也轉(zhuǎn)化都勾板或者DRM通道或者其他別的地方。(圖2)
2.4 偏濾暗盒模型的偏差計算結(jié)果
在真實的運行環(huán)境中,暗盒的三維模型是不能完全反映其模型形狀的。當(dāng)DCM處理終端感應(yīng)器并停留在環(huán)形通道上時,它會傾斜。(圖3-5)
DCM的形變分別用分析軟件和CATIA有限元建模工具來計算。這兩種結(jié)果會被比較。如果限定條件比較正確、全面,那么兩種工具的分析結(jié)果是相似的。
在下一階段,有限元分析結(jié)果會被分解。然后勾板的水平和垂直偏轉(zhuǎn)會與DTP2實驗室中真實的DCM測試結(jié)果比較(平臺1)。這種測試裝置是Sokkia NET05高精度三維調(diào)試系統(tǒng)。
Fig. 5. Vertical and horizontal deflections in respect to cassette structures.
在有限元分析結(jié)果和Sokkia NET05測試結(jié)果比較后,得出分析結(jié)論。
2.5 偏濾暗盒的偏差的可視化
根據(jù)機器的適用性和標(biāo)準(zhǔn)性原則設(shè)計了DCM。結(jié)果,其壓力總是在建筑材料的比例極限之下,并且具有一定的線彈性。初次測是在胡可定律的線彈性假設(shè)下進(jìn)行的。
因此,有限元分析結(jié)果是可以應(yīng)用于DCM形變的可視化的。加載裝置的形狀必然會反映到遠(yuǎn)程觀測系統(tǒng)的數(shù)據(jù)中。計算變形的遠(yuǎn)程可視化可以在兩種不同的方法下進(jìn)行。傳統(tǒng)的方法是把一個整體分成碎片,并在這些碎片間建立聯(lián)系[4]。這種方法需要大量的分析工作。鏈接的位置和最大鏈接就是這種分析的結(jié)果。基于這種分析結(jié)果,它被DCM分為三段鏈,并用兩個回轉(zhuǎn)節(jié)鏈接起來。圖6
本文提到的方法是運用3D變形——3D模型逐漸改變的過程——基于有限元分析結(jié)果去描述形變。形變,或者是3D變形,是物體從一種形狀變?yōu)榱硪环N形狀的過程[2]。這種技術(shù)可以直接使用有限元分析結(jié)果而不用麻煩的求的近似值。另外,這種方法能夠運用有限元分析出的每單位范圍內(nèi)的作用結(jié)果(圖7)。這就提供了一個更高層次的應(yīng)用能力,以適用于那些接受多個外力影響的復(fù)雜系統(tǒng)。
Fig. 6. The body of the cassette is divided into three rigid links connected with two rotational joints to approximate mechanical flexibilities.
Fig. 7. Simplified example of 3 links deformed by 9 individual morph targets (forces).
為了改變模型,我們使用直線切削沒變形的三維模型和有限元殘缺模型的高點。如果直線切削的精度不足的話,一個更先進(jìn)的變形算法是應(yīng)變場插入法。
運用達(dá)索系統(tǒng)可視化工具5.0來觀察變形(圖8)。這種虛擬環(huán)境是由ITER CATIA與有限元模型連接起來,用于創(chuàng)建變形范圍。
這種推薦方法的好處有以下幾點:
運用未加載荷狀態(tài)和變形狀態(tài)間的完全彈性變形,來直接使用有限元結(jié)果在每單位范圍內(nèi)的最大壓力。
更容易利用真實系統(tǒng)中離線和在線的變形結(jié)果。
更加精確的展示復(fù)雜系統(tǒng)的各個環(huán)節(jié),并且能完全控制連續(xù)變形的點,而不粗略的接近。
使從復(fù)雜系統(tǒng)中分離出單個外力因素引起的變形成為可能。
2.6 控制三維模型的變形
控制三維模型的變形意味著虛擬環(huán)境中的變形必須遵循現(xiàn)實環(huán)境中的變形。變形信息可以依據(jù)提前測量的運行狀態(tài)或者運用液壓系統(tǒng)的壓力來估測外力,因此,運用了現(xiàn)有的傳感器信息。
考慮到機器人的操作,一個更精確的方法可以通過采用應(yīng)變規(guī)來測量機器人鏈接的實際變形來達(dá)到。在實驗室的實驗中,應(yīng)變規(guī)將被安裝到DCM上。
應(yīng)變規(guī)的優(yōu)勢:
變形范圍方法的互補性,每一個應(yīng)力都可以通過專用的應(yīng)變規(guī)來單獨測量其范圍。
具有即時測量精確應(yīng)變的能力,不用依賴于提前測量的靜態(tài)變形或者不能對每一個應(yīng)力不能直接使用的液壓壓力
3. 未來工作
DTP2的三維虛擬樣機組件描述如下。最初,DCM彈性研究集中在解決子系統(tǒng)的彈性問題。然后,研究會包括整個CMM機器結(jié)構(gòu),并且鏈接結(jié)構(gòu)接近變形過程。
更多復(fù)雜的變形途徑,如三維領(lǐng)域里的線性插值法將會被使用。這將會聯(lián)合由真實DCM變形控制的三維樣機的彈性變形。
我們將進(jìn)一步分析柔性機器人關(guān)節(jié)的鏈接和分離。這將有助于創(chuàng)建更精確的有限元模型。最后,彈性鏈接的作用和影響會反映到整個系統(tǒng)中。
免責(zé)聲明
本文觀點不代表歐洲委員會的觀點。
致謝
這項工作是在EURATOM和TEKES的聯(lián)合契約下,由歐洲委員會支持,在EFDA工作框架下完成的。
參考文獻(xiàn)
[1] Internal Reports of Grant “DTP2 test facility operation and upgrade preparation”, 2010.
[2] J. Gomes, B. Costa, L. Darsa, L. Velho, Graphical objects, Visual Computer 12 (1996) 269–282.
[3] Han-Bing Yan, Shi-Min Hu, Ralph RMartin, 3D morphing using strain field interpolation, Journal of Computer Science and Technology archive 22 (1) (2007) 147–155.
[4] A.D. Luca, W. Book, Robots with flexible elements, in: Siciliano, Khatib (Eds.), Springer Handbook of Robotics, Springer, 2008, pp. 287–319.
Fusion Engineering and Design 86 (2011) 19581962 Contents lists available at ScienceDirect Fusion Engineering and Design journal homepage: A method for enabling real-time structural control results and 3D model morphing Sauli Krassi VTT article Article Available Keywords: ITER Remote Divertor Virtual Flexibility Deformation Morphing Virtual reality facility, DTP2 (Divertor Test Platform 2) has been established in Finland for the Remote Handling (RH) equipment designs for ITER. The first prototype the Cassette Multifunctional Mover (CMM) equipped with Second Cassette End to DTP2 in October 2008. The purpose is to prove that CMM/SCEE prototype can 2nd cassette RH operations. At the end of F4E grant “DTP2 test facility oper- the RH operations of the 2nd cassette were successfully demonstrated Fusion For Energy (F4E). CMM/SCEE robot has relatively large mechanical flexibilities when the robot 2nd Cassette on the 3.6-m long lever. This leads into a poor absolute where the 3D model, which is used in the control system, does not reflect the actual deformed state of the CMM/SCEE robot. To improve the accuracy, the new method has been developed in order to handle the flexibilities within the control systems virtual environment. The effect of the load on the CMM/SCEE has been measured and minimized in the load compensation model, which is implemented in the control system software. The proposed method accounts for the structural deformations of the robot in the control system through the 3D model morphing by utilizing the finite 1. been the system. deformations ing, morph 2D load radial nel and position of 0920-3796/$ doi: element method (FEM) analysis for morph targets. This resulted in a considerable improvement of the CMM/SCEE absolute accuracy and the adequacy of the 3D model, which is crucially important in the RH applications, where the visual information of the controlled device in the surrounding environment is limited. 2010 Elsevier B.V. All rights reserved. Introduction The paper presents how the load compensation functions have implemented in the control system software to improve absolute accuracy and visualization accuracy of DTP2 control It also proposes a new method for accounting structural in DTP2 control system through 3D model morph- utilizing the finite element method (FEM) analysis results for targets. In addition to the actual component morphing, the texture morphing will be utilized for representing the structural per each component for better operator evaluation. During the 2nd cassette installation process, CMM travels into direction, towards the reactor, on top the maintenance tun- rails with the aid of an electric motor drive, Fig. 1. The lifting tilting motions in the vertical plane are used for controlling the and orientation of the cassette according to uphill profile the maintenance tunnel. The SCEE, which consists of the can- Corresponding author. Tel.: +358 504118792. E-mail address: sauli.kivirantavtt.fi (S. Kiviranta). Fig. 1. CMM and SCEE structural representation. tilever (CRO) and the hook-plate (HRO) rotational joints, is devoted to change the position and orientation of the cassette during the toroidal motion towards the place of the 2nd cassette. see front matter 2010 Elsevier B.V. All rights reserved. 10.1016/j.fusengdes.2010.11.015 system by utilizing offline simulation Kiviranta , Hannu Saarinen, Harri Mkinen, Boris Technical Research Centre of Finland, P.O. Box 1300, FI-33101 Tampere, Finland info history: online 30 December 2010 handling Test Platform 2 engineering abstract A full scale physical test demonstrating and refining RH equipment at DTP2 is Effector (SCEE) delivered be used successfully for the ation and upgrade preparation”, to the representatives of Due to its design, the carries the nine-ton-weighting accuracy and into the situation deformation in remote handling S. Kiviranta et al. / Fusion Engineering and Design 86 (2011) 19581962 1959 2. Content of DTP2 deflection studies 2.1. CMM/SCEE trials After delivery of CMM/SCEE to DTP2, the system integration phase was started in order to prepare the system for the actual test- ing. This testing was done in two phases starting from the factory floor and ending to the RH operated tests from the control room 1. The initial motion programs for the cassette operations were done by teaching, while having a continuous visual contact with the cassette. Good repetition accuracy (3 mm) of the CMM/SCEE guaranteed successful repetition of the motion programs. However, static 3D model could not support operations properly, because of poor absolute accuracy of CMM/SCEE. On the grounds of 3D model, it seemed that the cassette was colliding with Divertor Region Mock-up (DRM) structure although in practice everything was fine. It was very clear, that absolute accuracy of the sys- tem should be improved before the remote operations could be started. 2.2. Load compensation The effect of load to the CMM-SCEE kinematic chain (body, wheels, links and joints) was measured during the motion pro- grams. And it was realized that the positioning error at the tip of the cassette was in the worst case close to 80 mm. The measure- ment data was utilized for creating load compensation functions to improve the absolute accuracy. The solution is general for the RH maintenance tunnel operations, but for the toroidal operations the load compensation model is a look-up table based on the values of the CRO joint. The compensation approach is simple and well reasoned because of generality for CMM to support future CMM operations with other end effectors. The specific look-up tables are used only when operating with SCEE for performing a spe- cific toroidal trajectory in and out. The compensation functions helped to improve the performance of the equipment considerably. Thus, the positioning error at the furthest point of the cassette was reduced from almost 80 mm to about 5 mm 1. Theimplementationofloadcompensationintothecartesianref- erence values can be seen in Fig. 2. The solution is divided into two phases depending on whether the cassette is loaded to the HRO or not. Cartesian reference for ideal equipment (Fig. 2) is expressing the location of a coordinate system (Fig. 1) which coincides with the axis of the HRO joint. Thereby, HRO joint can be used only for changing the orientation of the cassette around the vertical axis and only the CRO joint can be used to reach a y-coordinate value of the reference data. For this reason, the load compensation functions during the toroidal motions depend on the CRO joint. If there is no load in HRO joint, an inverse kinematics solution can be used directly to solve corresponding values for the joint ref- erences. The solution is calculated using the DenavitHartenberg parameters, which include corrections based on the signature cal- ibration of the CMM/SCEE 1. When the load is attached to the HRO joint, the inverse kine- matics solution cannot be used directly because in this case the cartesian reference includes also the components that rep- resent the deflections of CMM/SCEE. When the effect of load is known, the correct value for CRO joint can be found either iteratively utilizing the inverse of the load compensation or by defining the least-square polynomial fit between the measured y- coordinate values and the corresponding CRO joint values. Both iterative solution and 7th order polynomial fit are working well in practice. After the CRO joint value is defined, the position compensation in the x-, y- and z-directions and the orientation compensation in the Roll- (R) and Pitch- (P) directions can be made with respect to the cartesian reference. The compensation movement in the Yaw (W) direction cannot be done because of lacking the ability to move in the yaw-direction with the CMM/SCEE. Fig. 2. Left: load compensation in cartesian space. Right: implementation of load effect to the joint data of the real device. 1960 S. Kiviranta et al. / Fusion Engineering and Design 86 (2011) 19581962 2.3. an This Because which direction. tacting over weight to 2.4. Mock-up Cassette CAD end-effectors ware were results Then and that The dinate Fig. 3. Ansys FEM result (DCM lifted from RH interface). Fig. 4. CATIA FEM result (DCM lifted from RH interface). Improving teleoperator visualization accuracy Tilting the cassette in the yaw-direction can be visualized when additional joint is added to the 3D model of CMM/SCEE, (Fig. 2). joint has been placed between the hook plate and cassette. of that, the operator can see the effect of the cassette tilting, is 10 mm at the end of toroidal movement in the vertical To increase the visualization accuracy, when the cassette is con- the inner and outer rails of DRM, the pressure difference the Lift cylinder provides estimation about the load, as the of the cassette is gradually transferred from the hook plate the DRM rails or the other way around, (Fig. 2). Calculation of the deflections of the Divertor Cassette In the real operation environment, the shape of the Divertor Mock-up (DCM) is never equally represented by the 3D- model. The DCM deflects, when it is handled with the CMM and when it rests on the toroidal rails (Figs. 35). The deformations of the DCM were calculated using Ansys soft- and CATIA FEM-tools. The results of these two calculations compared. In conclusion, both FEM tools provide similar if the restraints are specified correctly. In the next phase, the FEM results were divided to components. horizontal and vertical deflections of the hook plate handling resting on the toroidal rails were compared to measurements had been done for real DCM in the DTP2 laboratory (Table 1). measurement device is Sokkia NET05 high precision 3D coor- measuring system (theodolite). Fig. 5. Vertical and horizontal deflections in respect to cassette structures. Comparison between the FEM results and the Sokkia measure- ments showed that the real DCM behaves as it was analyzed. 2.5. Visualization of the deflections of the Divertor Cassette Mock-up Design of the DCM has been made according to applicable design rules and standards of the machine design. As a result, the stresses are always below the proportional limit of the construction mate- rial and the behavior of material is linearly elastic. The initial tests in this study were conducted under the Hookes law assumption for linear deformations. Hence the results of the FEM analysis can be utilized for the visualization of the DCM deformations. Problem of loaded device shape not being reflected to the teleoperator view makes accurate controlofthesystemnearimpossible.Teleoperatorvisualizationby accounting deformations can be carried out in two different ways. The traditional method is to divide a body into pieces and to create the link mechanisms between the pieces 4. This approach requires a lot of analysis work. The position of the joints and the maximum joint values are the result of these analyses. Based on the analyses, it was recommended to divide DCM into three links, which were connected with two rotational joints, Fig. 6. Method proposed by this paper is to use the 3D morphing the process of gradual transformation between 3D bodies to describe the deformations of the body based on the FEM analysis results. The metamorphosis or the (3D) morphing of the 3D graphical objects, also known as shape interpolation, is the process of transform- ing one shape into another 2. This technique allows utilization of the FEM analysis results directly without laborious link-joint approximations. In addition, this method enables the use of sep- Table 1 FEM results compared to Sokkia measurements. DCM deflections Criteria Horizontal Vertical Unit Total deformation based on FEM between hook plate handling and resting on toroidal rails 7.2 9.2 mm Sokkia measurements between hook plate handling and resting on the toroidal rails 6.6 7.3 mm S. Kiviranta et al. / Fusion Engineering and Design 86 (2011) 19581962 1961 Fig. rotational Fig. (forces). arate target ( accounting multiple tion deformed strain tion Systems built the 6. The body of the cassette is divided into three rigid links connected with two joints to approximate mechanical flexibilities. 7. Simplified example of 3 links deformed by 9 individual morph targets deformation results by utilizing FEM results for each morph per force applicable to the component for a given scenario Fig. 7). This provides a high degree of adaptation capabilities for a large variety of flexibilities in complex systems where sources of forces can affect each part of the system. For morphing the model, we have used the linear interpola- between the vertices in the non-deformed 3D model and the FEM model. A more advanced morphing algorithm is the field interpolation 3 if the accuracy of the linear interpola- is not sufficient for the given application. For visualizing the deformations to the teleoperator, Dassault Virtools 5.0 was used (Fig. 8). The virtual environment is by directly utilizing ITER CATIA models in conjunction with FEM models that are utilized for creating morph targets. The benefits of the proposed method are the following: Applicationofthefullyflexiblemeshmorphingbetweenthecom- ponents unloaded neutral states and the deformed states for a given maximum force per morph target thus directly utilizing the FEM results, Fig. 4. Easier reuse of the existing deformation data obtained by the offline and online analysis of real systems. More accurate representation of each section of the complex sys- tem components and full control over the continuum possible deformation points, instead of rough estimations gained trough joint-link approximation. Possibility to combine multiple deformations separated by indi- vidual forces in complex system. Fig. 8. Example of DCM morph targets within virtual environment. 2.6. Controlling of the 3D model deformations Controlling the 3D model deformations means that the defor- mations in the virtual environment have to follow the actual deformations. The deformation information can be determined based on the previously measured deformations for a given operation state or by using the hydraulic system pressures to estimate the force, hence utilizing existing sensor informa- tion. In the case of robot operations, a more accurate solution can be achieved by employing strain gauges to measure the actual defor- mations of the robot links. In the laboratory tests, the strain gauges will be installed to the DCM. The advantage of the strain gauges includes: Complementarity to the morph target method, where for each force one can have an individual morph target controlled by a dedicated strain gauge. Ability to measure the exact deformations immediately, not relying on the previously measured static deformation data or hydraulic pressures that may not be available for all force direc- tions. 3. Future work The continuation regarding the fully flexible 3D virtual proto- typing of the DTP2 robot components will be as follows. Initially, the DCM flexibility studies will focus on solving the flexibility prob- lems of subsystems. Later, the studies will cover the whole CMM robot structures and will result in appropriate joint structures of the rigid objects to reflect the deformed states of the preceding objects in the link chain. More sophisticated morphing methods such as the 3D strain fields in replacement of the linear interpolation will be applied. This will be combined with the flexible 3D virtual prototype to be controlled by the strain gauges measurements of the real DCM deformations. We will further analyze the combined and separate effect of the flexibility of the robot links and the joints. This will help to create a more precise FEM models. Finally, the influence of the link flexibilities on the dynamic response of the system will be explored. 1962 S. Kiviranta et al. / Fusion Engineering and Design 86 (2011) 19581962 Disclaimer The views and opinions expressed herein do not necessarily reflect those of the European Commission. Acknowledgements This work is supported by the European Commission under the contract of association between EURATOM/TEKES, was carried out within the frame-work of the European Fusion Development Agreement. References 1 Internal Reports of Grant “DTP2 test facility operation and upgrade preparation”, 2010. 2 J.Gomes,B.Costa,L.Darsa,L.Velho,Graphicalobjects,VisualComputer12(1996) 269282. 3 Han-Bing Yan, Shi-Min Hu, Ralph R Martin, 3D morphing using strain field inter- polation, Journal of Computer Science and Technology archive 22 (1) (2007) 147155. 4 A.D. Luca, W. Book, Robots with flexible elements, in: Siciliano, Khatib (Eds.), Springer Handbook of Robotics, Springer, 2008, pp. 287319.