基于有限段方法的SMA扭轉(zhuǎn)驅(qū)動器柔性機(jī)構(gòu)動力分析

鄒秀清 董二寶 張世武 許 旻 楊 杰

中國科學(xué)技術(shù)大學(xué)，合肥，230027

0 引言

柔性結(jié)構(gòu)［1］是一種依靠材料的彈性變形來傳輸運(yùn)動或力的裝置，相比傳統(tǒng)的剛性鉸鏈機(jī)構(gòu)，柔性機(jī)構(gòu)在簡化結(jié)構(gòu)、減少磨損和提高精度上都有顯著的優(yōu)勢。柔性機(jī)構(gòu)通常被應(yīng)用在精密工程場合以及仿生機(jī)械及機(jī)器人等領(lǐng)域［1-2］。

在高速、精密以及高性能MEMS應(yīng)用場合，柔性機(jī)構(gòu)的動力學(xué)問題十分突出。因此，柔性機(jī)構(gòu)的動力學(xué)特性分析方法的研究是柔性機(jī)構(gòu)學(xué)研究的一個重要方向。

文獻(xiàn)［3］考慮柔性體變形的動力剛化現(xiàn)象，提出了在有限元方法中引入單元耦合形函數(shù)，再根據(jù)Kane方程得到包含動力剛度項(xiàng)的一致線性化動力學(xué)方程，但是該方法是在結(jié)構(gòu)小變形的前提下推導(dǎo)得出的，對大變形問題并不適用。

文獻(xiàn)［4］針對全柔性機(jī)器人的動力學(xué)分析總結(jié)了幾種方法:一種是偽剛體法，這種方法雖然可以使柔性體的動力學(xué)建模得到很大簡化，但是簡化造成的誤差較大，且復(fù)雜柔性機(jī)構(gòu)的偽剛體模型的建立也非常困難;一種是利用有限元軟件建立柔性系統(tǒng)的有限元模型，這種方法只限于小變形尺度，可以忽略位形變化對機(jī)構(gòu)動態(tài)特性的影響;還有一種就是應(yīng)用傳統(tǒng)的試驗(yàn)方法進(jìn)行分析，這種方法效率低、耗時長、成本高。

柔性體動力學(xué)的有限段方法是基于多剛體動力學(xué)理論發(fā)展起來的柔性體動力學(xué)分析方法，自然計及幾何非線性的影響［5］。對于梁式柔性機(jī)構(gòu)，可以方便地應(yīng)用有限段方法進(jìn)行離散建模，再利用原有的多剛體動力學(xué)理論進(jìn)行計算分析。

本文采用有限段模型［6］對形狀記憶合金(SMA)扭轉(zhuǎn)驅(qū)動器的柔性結(jié)構(gòu)部分進(jìn)行離散建模，結(jié)合ADAMS動力學(xué)分析軟件在系統(tǒng)動力學(xué)分析上的優(yōu)勢，對SMA扭轉(zhuǎn)驅(qū)動器的柔性結(jié)構(gòu)進(jìn)行動力學(xué)分析。

1 有限段原理介紹

有限段方法的基本思想是:將梁式柔性構(gòu)件離散為有限的剛性梁段，段間用3個扭轉(zhuǎn)彈簧、3個線性彈簧及與之相應(yīng)的6個阻尼器相連接，如圖1所示。在小應(yīng)變有限變形的條件下，可用離散的梁段來描述柔性體的慣性特征，用段間的彈簧和阻尼器代表柔性體的彈性和阻尼特性［7］。

圖1 梁式柔性機(jī)構(gòu)的有限段模型

在小應(yīng)變有限變形的情況下，雖然整體的變形比較大，但由于是小應(yīng)變，每一小段的變形仍然符合小變形條件，所以有限段模型在小變形情況和小應(yīng)變大變形情況都是適用的。

為了推導(dǎo)梁段之間彈簧的剛度特性，我們首先考察懸臂梁自由端的位移與作用其上的力(矩)的關(guān)系。

根據(jù)線性梁的理論［8］，在懸臂梁的自由端施加力和力矩，在變形較小的情況下，懸臂梁自由端的位移和轉(zhuǎn)動分量與施加的力和力矩之間的關(guān)系如下。

首先，左端固定情況如圖2所示，懸臂梁自由端上施加的力F和力矩M與位移u和轉(zhuǎn)動分量θ之間的關(guān)系可以表示為

式中，l為梁長度;A為梁截面積;J為梁截面積對X軸的形心二次矩;Iy為梁截面積對Y軸的形心二次矩;Iz為梁截面積對Z軸的形心二次矩;E為材料彈性模量;G為材料剪切模量。

圖2 右端點(diǎn)2作用一般載荷的懸臂梁

同理可以得出右端固定情況下的力和力矩與位移和轉(zhuǎn)動分量之間的關(guān)系，限于篇幅，這里不再一一列出。

然后，考察柔性梁的彈簧連接有限段模型。為了建立精確的梁模型，彈簧必須精確表示相鄰段間的彈性。為推導(dǎo)這個表達(dá)式，考察圖3所示的兩個相鄰段Bj和Bk。其中，連接彈簧分別由兩個代表相應(yīng)段彈性的部分組成，同時，連接彈簧代表相應(yīng)段的伸張、彎曲和扭轉(zhuǎn)特性。

圖3 兩個典型相鄰有限段

每段有12個與其相關(guān)的位移和轉(zhuǎn)動分量(與每個端點(diǎn)相關(guān)的有6個)。為了定義這些分量，考察如圖4所示的任一段Bj。令Gj為Bj的質(zhì)心。令Oj和Ok為Bj的彈簧與相鄰段Bi的彈簧的連接點(diǎn)以及Bk的彈簧與相鄰段Bj的彈簧的連接點(diǎn)。令OI、OJ和 OK分別為以 Oj、Gj和 Ok為原點(diǎn)的參考系。令RJ的軸與Bj的主軸連續(xù)共線。RI和RK的軸當(dāng)彈簧未變形時與RJ的軸平行。則OI相對于OJ的位移和RI相對于RJ的轉(zhuǎn)動共有6個自由度。同理，OK相對于OJ的位移和RK相對于RJ的轉(zhuǎn)動共有6個自由度，即與Bj段相關(guān)的位移和轉(zhuǎn)動分量有12個。這些位移和轉(zhuǎn)動變量的定義如表1所示。其中，Rk、Rj分別為以O(shè)K和OI為原點(diǎn)的參考系。

圖4 典型有限段模型

表1 位移和轉(zhuǎn)動變量

Rk相對于Rj的位移和轉(zhuǎn)動變量可表示為

根據(jù)式(1)、式(2)，可將 ukx、uky和 ukz表示為

因?yàn)樵贠k節(jié)點(diǎn)處沒有質(zhì)量，也不受外力，所以有:

所以，ukx、uky和 ukz可以表示為

同樣，相對轉(zhuǎn)動分量可以表示成

由上可解出 Fj2x、Fj2y、Fj2z、Mj2x、Mj2y和 Mj2z:

其中，系數(shù) a、b、c與相應(yīng)段的結(jié)構(gòu)、材料參數(shù)有關(guān)，這里不一一列出，可參考文獻(xiàn)［7］。

到此，我們得出了各個段之間的彈簧剛度特性。

2 基于ADAMS宏命令的有限段模型建立與數(shù)值計算

為了探討有限段模型的準(zhǔn)確性，我們需要建立梁的有限段模型，為了分析梁的變形情況，需要建立有限段模型的多體系統(tǒng)運(yùn)動學(xué)及動力學(xué)方程，如文獻(xiàn)［9-10］所描述的，需要計算出每個有限段的廣義速率，整理出整個系統(tǒng)的動力學(xué)方程組再求解微分方程組而得到系統(tǒng)的特性，模型的建立及求解過程非常復(fù)雜而且對于復(fù)雜模型更是無能為力。

本文應(yīng)用ADAMS軟件建立有限段多剛體模型進(jìn)行分析，把動力學(xué)方程的建立和求解過程交給軟件處理，我們需要做的就是在ADAMS中準(zhǔn)確地建立系統(tǒng)的有限段模型。

要在ADAMS中建立有限段模型，關(guān)鍵是創(chuàng)建各個有限段以及在各有限段之間施加的力和力矩。根據(jù)有限段模型的特點(diǎn)，有限剛體段和力與力矩的建立是一個重復(fù)性的過程，我們利用ADAMS宏命令的循環(huán)語句建立模型。先用循環(huán)命令建立所有的有限段，然后再用一個循環(huán)命令建立各段之間的力和力矩。具體的建立方法如下:如要在第i個有限段part i和第i+1個有限段part i+1之間建立有限段的彈簧力和彈簧力矩，先分別在part i的右端面和part i+1的左端面創(chuàng)建標(biāo)記點(diǎn) marker i和 marker i+1，再在標(biāo)記點(diǎn)marker i和marker i+1之間創(chuàng)建 VFORCE(X、Y、Z三向的直線力)和VTORQUE(X、Y、Z三向的扭轉(zhuǎn)力)。其中，段 part i和段 part i+1之間的力VFORCE的表達(dá)式根據(jù)式(11)、式(13)、式(15)建立，ukx、uky、ukz分別是 marker i+1 相對于marker i的X、Y、Z軸的移動位移。同樣，段part i和段part i+1之間的力VTORQUE的表達(dá)式根據(jù)式(12)、式(14)、式(16)建立，θkx、θky、θkz分別是marker i+1相對于marker i的X、Y、Z軸的轉(zhuǎn)動位移。

下面，我們探討有限段模型的準(zhǔn)確性。首先，比較直梁的有限段模型和有限元模型在恒力作用下的變形。對厚度h=4mm、寬度b=2mm、長度L=80mm的懸臂直梁，在ADAMS中建立直梁的有限段離散模型，如圖5所示，各剛體段間施加的力和力矩按以上有限段模型計算的剛度乘以變形時段間的相對位置變化得到。分析結(jié)果如圖6所示，從圖中可以看出，隨著有限段模型不斷細(xì)化，有限段模型的分析結(jié)果不斷靠近有限元模型的分析結(jié)果。

圖5 直梁的有限段模型

圖6 直梁靜力分析比較

然后，比較曲梁的有限段模型和有限元模型在恒力作用下的變形。在曲梁的有限段模型中，相鄰兩有限段之間的中性軸并不是平行的，而是有一定曲率夾角存在的，在推導(dǎo)各段之間的彈簧剛性時存在更復(fù)雜的耦合情況。但是，隨著有限段的不斷細(xì)化，相鄰兩段之間的夾角不斷減小，對彈簧剛度的影響也不斷減小甚至可以忽略。因此，在曲梁應(yīng)用上，為了簡化有限段模型的段間彈簧剛度的推導(dǎo)，在計算段間剛度時，忽略相鄰兩段之間曲率的影響，將相鄰兩段近似為直梁段，段間剛度仍然按照直梁的方法推導(dǎo)。任意兩小段之間近視為直線段，軸向方向近視為兩段相鄰點(diǎn)處的切線方向。隨著有限段離散的細(xì)化，兩小段之間的傾角不斷減小，這種由近似造成的誤差也將不斷減小。對厚度 h=4mm、寬度 b=2mm、半徑R=50mm的1/4圓弧，在ADAMS中的模型如圖7所示，建模方法與直梁一致。分析結(jié)果如圖8所示，從圖中可以看出，隨著有限段模型不斷細(xì)化，有限段模型的分析結(jié)果也不斷靠近有限元模型的分析結(jié)果，且誤差與直梁的有限段模型相差不大。

基于以上比較結(jié)果可知，盡管有限段模型在曲梁分析上的應(yīng)用存在忽略梁段之間曲率的簡化，但是對有限段離散足夠細(xì)的模型，其分析結(jié)果仍是可靠的。所以，我們認(rèn)為有限段模型在一定范圍內(nèi)的分析應(yīng)用是可靠的。

圖7 曲梁的有限段模型

圖8 曲梁靜力分析比較

3 SMA扭轉(zhuǎn)驅(qū)動器的動態(tài)特性分析

SMA扭轉(zhuǎn)驅(qū)動器的工作原理如圖9所示，SMA合金絲按一定規(guī)則纏繞在柔性結(jié)構(gòu)及外圓環(huán)上，加熱SMA合金絲使其收縮帶動柔性結(jié)構(gòu)發(fā)生變形，從而使中間傳動軸孔轉(zhuǎn)過一定的角度。由于本文的主要工作集中在有限段模型建立和SMA扭轉(zhuǎn)驅(qū)動器動態(tài)特性分析上，SMA扭轉(zhuǎn)驅(qū)動器中柔性結(jié)構(gòu)的優(yōu)化設(shè)計以及SMA絲的布置問題在其他文章中有詳細(xì)介紹，這里不再贅述。

圖9 SMA扭轉(zhuǎn)驅(qū)動器結(jié)構(gòu)圖與實(shí)物

本文結(jié)合有限段模型和ADAMS動力學(xué)分析軟件來分析SMA扭轉(zhuǎn)驅(qū)動器的動態(tài)特性。SMA扭轉(zhuǎn)驅(qū)動器的柔性結(jié)構(gòu)在優(yōu)化設(shè)計的過程中是采用B樣條曲線來描述曲線形狀的，在MATLAB中離散成若干個坐標(biāo)點(diǎn)。在建立有限段模型時，為了方便，直接以這些離散坐標(biāo)點(diǎn)為離散位置進(jìn)行劃分。由式(11)～式(16)計算各段之間的彈性特性，在ADAMS中通過應(yīng)用宏命令建立有限段多剛體模型的方法建立SMA扭轉(zhuǎn)驅(qū)動器的離散有限段模型。

首先分析靜力情況。在ADAMS中建立SMA扭轉(zhuǎn)驅(qū)動器的有限段模型，如圖10所示。在緩慢施加的對稱恒力10N的載荷作用下，柔性梁厚度為0.4mm的扭轉(zhuǎn)驅(qū)動器的輸出轉(zhuǎn)角如圖11所示。圖12顯示不同厚度柔性梁的驅(qū)動器輸出轉(zhuǎn)角。

圖10 SMA扭轉(zhuǎn)驅(qū)動器的有限段模型

圖11 SMA扭轉(zhuǎn)驅(qū)動器輸出轉(zhuǎn)角曲線

圖12 不同厚度下的輸出轉(zhuǎn)角

然后，在動載荷情況下觀察扭轉(zhuǎn)驅(qū)動器的扭轉(zhuǎn)情況。分別施加脈沖載荷和階躍載荷，觀察驅(qū)動器的響應(yīng)情況。

圖13顯示柔性梁厚度為0.5mm時，在脈沖載荷作用下的時域和頻域響應(yīng)結(jié)果。圖14顯示不同厚度的柔性梁在脈沖載荷下的響應(yīng)頻率，從圖14中可以看出厚度越大，響應(yīng)頻率越高。

在階躍載荷作用下，柔性梁厚度為0.5mm時的扭轉(zhuǎn)響應(yīng)如圖15所示。

4 結(jié)論

(1)有限段方法是對梁式柔性機(jī)構(gòu)進(jìn)行動力學(xué)分析的有效方法。

(2)利用ADAMS宏命令在ADAMS中創(chuàng)建有限段模型系統(tǒng)，可以大大簡化多剛體系統(tǒng)動力學(xué)方程的建立過程。

圖13 脈沖載荷下扭轉(zhuǎn)響應(yīng)

圖14 不同厚度下的響應(yīng)頻率

圖15 階躍載荷下扭轉(zhuǎn)響應(yīng)

(3)將有限段方法應(yīng)用于SMA柔性扭轉(zhuǎn)驅(qū)動器的動力學(xué)分析，獲得了扭轉(zhuǎn)驅(qū)動器的相關(guān)動態(tài)特性。

［1］Her I，Midha A.A Compliance Number Concept for Compliant Mechanisms and Type Synthesis［J］.ASME Journal of Mechanisms，Transmissions，and Automation in Design，1987，109:348-355.

［2］王振海，張衛(wèi)紅.柔性結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計發(fā)展概況［J］.機(jī)械設(shè)計，2004，21(3):1-4.Wang Zhenhai，Zhang Weihong.Topological Optimization Design for Flexible Structure［J］.Journal of Machine Design，2004，21(3):1-4.

［3］劉又午，王樹新，王建明，等.復(fù)雜結(jié)構(gòu)航天器的柔性多體動力學(xué)分析［J］.航天器工程，1998，7(3):45-52.Liu Youwu，Wang Shuxin，Wang Jianming，et al.Flexible MultibodyDynamicAnalysisofComplex Structural Spacecraft［J］. Spacecraft Engineering，1998，7(3):45-52.

［4］于靖軍，畢樹生，宗光華，等.全柔性機(jī)器人機(jī)構(gòu)結(jié)構(gòu)動力學(xué)分析方法研究［J］.機(jī)械工程學(xué)報，2004，40(8):54-57.Yu Jingjun，Bi Shusheng，Zong Guanghua，et al.Research on Structural Dynamics of a General Fully Compliant Robotic Mechanism［J］.Journal of Mechanical Engineering，2004，40(8):54-57.

［5］閆紹澤，劉冰清，劉又午，等.柔性多體系統(tǒng)動力學(xué)——有限段方法［J］.河北工業(yè)大學(xué)學(xué)報.1997，26(3):1-8.Yan Shaoze，Liu Bingqing，Liu Youwu，et al.Dynamics of Flexible Multibody Systems－The Finite Segment Method［J］.Journal of Hebei University of Technology，1997，26(3):1-8.

［6］Huston R L.Multibody Dynamics Including the Effects of Flexibility and Compliance［J］.Computer and Science，1981，14(5/6):443-451.

［7］李東旭.撓性航天器結(jié)構(gòu)動力學(xué)［M］.北京:科學(xué)出版社，2010.

［8］劉鴻文.材料力學(xué)［M］.北京:高等教育出版社，2004.

［9］金國光，劉又午，張明成，等.含柔性梁機(jī)械系統(tǒng)動力學(xué)分析的有限段方法［J］.鄭州工業(yè)大學(xué)學(xué)報，2000，21(1):23-25.Jin Guoguang，Liu Youwu，Zhang Mingcheng，et al.Finite Segment Method of Dynamic Analysis of Mechanical System with Beam－like Flexible Bodies［J］.Journal of Zhengzhou University of Technology，2000，21(1):23-25.

［10］董龍雷，閆桂榮，余建軍，等.帶中心剛體的旋轉(zhuǎn)柔性件有限段建模［J］.西安交通大學(xué)學(xué)報，2001，35(9):977-981.Dong Longlei，Yan Guirong，Yu Jianjun，et al.Finite Segment Modeling of a Rotating Flexible Body with Central Rigid Body［J］.Journal of Xi’an Jiaotong University，2001，35(9):977-981.