基于絕對(duì)節(jié)點(diǎn)坐標(biāo)法的柔性多體系統(tǒng)靈敏度分析

王鐵成， 陳國(guó)平， 孫東陽

(1.南京航空航天大學(xué) 機(jī)械結(jié)構(gòu)力學(xué)及控制國(guó)家重點(diǎn)實(shí)驗(yàn)室,南京　210016; 2.重慶大學(xué) 航空航天學(xué)院 ，重慶　400044)

?

基于絕對(duì)節(jié)點(diǎn)坐標(biāo)法的柔性多體系統(tǒng)靈敏度分析

王鐵成1， 陳國(guó)平1， 孫東陽2

(1.南京航空航天大學(xué) 機(jī)械結(jié)構(gòu)力學(xué)及控制國(guó)家重點(diǎn)實(shí)驗(yàn)室,南京210016; 2.重慶大學(xué) 航空航天學(xué)院 ，重慶400044)

機(jī)械系統(tǒng)優(yōu)化設(shè)計(jì)可以提高機(jī)械系統(tǒng)的精度和穩(wěn)定性。為了避免優(yōu)化參數(shù)選擇的盲目性，提高設(shè)計(jì)效率，進(jìn)行近似分析和預(yù)測(cè)設(shè)計(jì)趨勢(shì)，首先就需要對(duì)系統(tǒng)參數(shù)進(jìn)行靈敏度分析。因此，靈敏度分析已經(jīng)成為機(jī)械系統(tǒng)動(dòng)力學(xué)分析與優(yōu)化設(shè)計(jì)的橋梁，并成為機(jī)械系統(tǒng)優(yōu)化設(shè)計(jì)的核心問題。

靈敏度計(jì)算方法主要有直接微分法和伴隨變量法。直接微分法是通過計(jì)算系統(tǒng)的運(yùn)動(dòng)學(xué)或動(dòng)力學(xué)方程對(duì)設(shè)計(jì)變量的偏導(dǎo)數(shù)，得到目標(biāo)函數(shù)對(duì)設(shè)計(jì)變量的靈敏度。潘振寬等[1]采用直接微分法對(duì)多剛體系統(tǒng)進(jìn)行了靈敏度分析。針對(duì)遞推形式建立的鏈狀多體系統(tǒng)，Haug等[2]采用直接微分法對(duì)其進(jìn)行了靈敏度分析。直接微分法的優(yōu)點(diǎn)是直觀，易于理解，但隨著系統(tǒng)設(shè)計(jì)變量的增多，求解規(guī)模成倍增長(zhǎng)，導(dǎo)致靈敏度求解的效率降低。伴隨變量法是通過引入伴隨變量，使多設(shè)計(jì)變量系統(tǒng)靈敏度分析計(jì)算效率得到提高，該方法已經(jīng)在多體系統(tǒng)的靈敏度分析中得到了重視和應(yīng)用[3]。Etman等[4]采用伴隨變量法對(duì)多體系統(tǒng)動(dòng)態(tài)特性進(jìn)行了靈敏度分析。Li等[5]和Maly等[6]基于隱式微分代數(shù)方程，推導(dǎo)了多體系統(tǒng)靈敏度分析的伴隨變量法。目前，多體系統(tǒng)靈敏度分析大多是針對(duì)剛體系統(tǒng)，柔性多體系統(tǒng)靈敏度分析的研究還比較少。主要原因是柔性多體系統(tǒng)具有高維數(shù)和強(qiáng)非線性特性，導(dǎo)致分析相對(duì)困難。雖然浮動(dòng)坐標(biāo)法是最常用的一種柔性多體系統(tǒng)建模方法，但是該方法基于小變形假設(shè)，用于分析存在大變形的柔性多體系統(tǒng)會(huì)帶來較大誤差。Shabana等[7]提出了多柔性體絕對(duì)節(jié)點(diǎn)坐標(biāo)建模方法，該方法的理論基礎(chǔ)是有限元和連續(xù)介質(zhì)力學(xué)[8-9]，由該方法推導(dǎo)的多體系統(tǒng)方程的質(zhì)量陣為常數(shù)矩陣，且不存在科氏力和離心力項(xiàng)。為了對(duì)存在大變形的柔性多體系統(tǒng)進(jìn)行靈敏度分析，同時(shí)考慮到多變量情況下伴隨變量法有更高的計(jì)算效率，本文采用伴隨變量法對(duì)基于絕對(duì)節(jié)點(diǎn)坐標(biāo)法建立的柔性多體系統(tǒng)進(jìn)行了靈敏度分析，并將其計(jì)算結(jié)果與直接法的計(jì)算結(jié)果進(jìn)行了比較。

1絕對(duì)節(jié)點(diǎn)坐標(biāo)建模的動(dòng)力學(xué)方程

基于絕對(duì)節(jié)點(diǎn)坐標(biāo)法的一維兩節(jié)點(diǎn)梁?jiǎn)卧鐖D1所示，單元上任意點(diǎn)的位置矢量可表示為：

(1)

式中：e為單元節(jié)點(diǎn)坐標(biāo)，S是定義在總體坐標(biāo)系上的形函數(shù)，分別表示如下[7]

(2)

圖1　平面梁?jiǎn)卧狥ig.1Plane beam element

(3)

式中:s1=1-3ξ2+2ξ3，s2=ξ-3ξ2+2ξ3，s3=3ξ2-2ξ3，s4=l(ξ3-ξ2)，這里ξ=x/l。

根據(jù)式(1)，單元的動(dòng)能可表示為

(4)

式中:Me=∫ρSTSdV，為單元的常數(shù)質(zhì)量陣，ρ和V分別是材料的密度和單元的體積。

根據(jù)連續(xù)介質(zhì)力學(xué)理論，單元的總應(yīng)變能包含彎曲應(yīng)變能Uel和軸向拉伸應(yīng)變能Uet，表示為

(5)

式中：εi和κ分別為單元應(yīng)變和曲率。

基于虛功原理建立單元的動(dòng)力學(xué)方程表示為

(6)

其中:Qe為單元受到的廣義外力，Qk為單元廣義彈性力。

單元廣義彈性力由單元應(yīng)變能對(duì)單元坐標(biāo)求偏導(dǎo)獲得，即

(7)

式中：Kl和Kt分別表示單元的軸向拉伸非線性剛度陣和彎曲線性剛度陣。

則含約束的柔性體k的動(dòng)力學(xué)方程為

(8)

其中:Be為布爾矩陣。

2靈敏度分析

伴隨變量法計(jì)算靈敏度，需要先求解動(dòng)力學(xué)方程，然后根據(jù)動(dòng)力學(xué)計(jì)算結(jié)果，逆向求解伴隨變量為未知量的微分方程組，其中，動(dòng)力學(xué)計(jì)算結(jié)果需要擬合函數(shù)，可使在伴隨變量微分方程的數(shù)值計(jì)算時(shí)，能取得在所計(jì)算時(shí)間內(nèi)的任意時(shí)刻的函數(shù)值。

本文選用積分型目標(biāo)函數(shù)對(duì)絕對(duì)節(jié)點(diǎn)坐標(biāo)法建模的柔性多體系統(tǒng)進(jìn)行靈敏度分析。目標(biāo)函數(shù)設(shè)為

(9)

狀態(tài)變量和設(shè)計(jì)變量分別為

q=[q1q2…qn]T

(10)

b=[b1b2…bl]T

(11)

多體系統(tǒng)動(dòng)力學(xué)方程和約束方程分別為

(12)

Φ(q,b,t)=0

(13)

為了計(jì)算靈敏度，先用萊布尼茲法對(duì)目標(biāo)函數(shù)求導(dǎo)，得到

(14)

(15)

(16)

引入伴隨變量μ和ν，分別乘以式(12)和式(13)并積分得

(17)

(18)

式(17)和式(18)對(duì)設(shè)計(jì)變量求導(dǎo)得

ΦTqbλ+ΦTqqλqb-Qb-Qqqb)dt=0

(19)

(20)

(21)

其中

將式(15)，式(20)和式(21)相加得

(22)

選取伴隨變量使得積分式中qb和λb前的系數(shù)分別等于零，得如下關(guān)于伴隨變量的微分代數(shù)方程

(23)

Φqμ-Fλ=0

(24)

(25)

再由式(25)推出下面兩式

(26)

(27)

采用絕對(duì)節(jié)點(diǎn)坐標(biāo)法建立的柔性多體系統(tǒng)動(dòng)力學(xué)方程，質(zhì)量矩陣M為滿秩常數(shù)陣，進(jìn)一步可得到

(28)

(29)

根據(jù)初始條件式(28)、(29)，逆向求解關(guān)于伴隨變量的微分代數(shù)方程組(23)、(24)，將計(jì)算結(jié)果代入式(16)，得目標(biāo)函數(shù)靈敏度

(30)

3算例分析

本文以受重力作用的矩形截面柔性單擺系統(tǒng)為研究對(duì)象，如圖2所示。

圖2　柔性單擺Fig.2 Flexible simple pendulum

表1　兩種方法的靈敏度計(jì)算結(jié)果

表1給出了分別采用伴隨變量法和直接微分法得到的目標(biāo)函數(shù)對(duì)單擺的密度、彈性模量、長(zhǎng)度、截面長(zhǎng)度和寬度的靈敏度。同時(shí)對(duì)比了兩種方法計(jì)算結(jié)果的相對(duì)誤差，可以發(fā)現(xiàn)，其相對(duì)誤差都在4%以內(nèi)。而且可以發(fā)現(xiàn)，目標(biāo)函數(shù)對(duì)單擺長(zhǎng)度、截面長(zhǎng)度和寬度的變化比較敏感，而密度和彈性模量對(duì)目標(biāo)函數(shù)的影響很小。采用直接微分法計(jì)算的單擺末端橫向位移對(duì)五個(gè)設(shè)計(jì)變量的靈敏度如圖3~圖7所示，由圖可知五個(gè)設(shè)計(jì)變量中長(zhǎng)度的變化對(duì)末端橫向位移的影響最明顯，這與基于目標(biāo)函數(shù)的靈敏度分析結(jié)論一致。

圖3 x對(duì)ρ的靈敏度Fig.3Thesensitivityofxwithrespecttoρ圖4 x對(duì)E的靈敏度Fig.4ThesensitivityofxwithrespecttoE圖5 x對(duì)l的靈敏度Fig.5Thesensitivityofxwithrespecttol

圖6　x對(duì)a的靈敏度Fig.6 The sensitivity of x with respect to a

圖7　x對(duì)b的靈敏度Fig.7 The sensitivity of x with respect to b

為了分析伴隨變量法和直接微分法的計(jì)算效率，對(duì)比了五個(gè)設(shè)計(jì)變量情況下兩種方法的計(jì)算時(shí)間，如表2所示。 可以看出，考慮一個(gè)設(shè)計(jì)變量時(shí)，伴隨變量法所用計(jì)算時(shí)間多于直接微分法，約為直接微分法計(jì)算時(shí)間的2倍；考慮兩個(gè)設(shè)計(jì)變量時(shí)，兩種方法的計(jì)算時(shí)間基本相當(dāng)；三個(gè)設(shè)計(jì)變量時(shí)，伴隨變量法所用時(shí)間約為直接微分法的75%，當(dāng)五個(gè)設(shè)計(jì)變量的情況時(shí)，伴隨變量法的計(jì)算時(shí)間約為直接法計(jì)算時(shí)間一半。由此可見，隨著設(shè)計(jì)變量的增加，伴隨法的計(jì)算效率將高于直接微分法。這主要是因?yàn)殡S著設(shè)計(jì)變量的增多，伴隨變量法隨之增加的計(jì)算時(shí)間主要就在靈敏度積分函數(shù)上，在此算例中大概耗時(shí)25 s左右，而直接微分法隨著設(shè)計(jì)變量增加，需要再計(jì)算一組微分方程，勢(shì)必增加求解時(shí)間。

表2　兩種方法的計(jì)算時(shí)間

4結(jié)論

本文對(duì)基于絕對(duì)節(jié)點(diǎn)坐標(biāo)法建立的柔性多體系統(tǒng)進(jìn)行了靈敏度分析。以矩形截面的柔性單擺系統(tǒng)為例，以工程上常用的積分型函數(shù)為目標(biāo)函數(shù)，使用伴隨變量法和直接微分法分別計(jì)算了目標(biāo)函數(shù)對(duì)柔性體長(zhǎng)度、彈性模量、密度、截面的長(zhǎng)度和寬度的靈敏度。通過分析發(fā)現(xiàn)，兩種方法計(jì)算結(jié)果的相對(duì)誤差均在4%以內(nèi)，因此，計(jì)算結(jié)果是可靠的。同時(shí)可以發(fā)現(xiàn)隨著設(shè)計(jì)變量的增加，伴隨變量法比直接法有更高的計(jì)算效率。為工程上，采用伴隨變量法對(duì)存在大變形的柔性多體系統(tǒng)的進(jìn)行靈敏度分析提供了理論依據(jù)。

參 考 文 獻(xiàn)

[1] 潘振寬，丁潔玉，王鈺. 基于隱式微分/代數(shù)方程的多體系統(tǒng)動(dòng)力學(xué)設(shè)計(jì)靈敏度分析方法[J]. 動(dòng)力學(xué)與控制,2004, 2(2): 66-69.

PAN Zhen-kuan, DING Jie-yu, WANG Yu. Design sensitivity analysis of multibody system dynamics descrided by implicit differential/algebraic[J]. Journal of Dynamics and Control,2004, 2(2): 66-69.

[2] Haug E J,Neel K M, Krishnasawami P. Design sensitivity analysis and optimization of dynamically driven systems[M]. Computer Aided Analysis and Optim-ization of Mechanical System Dynamics.Springer Berlin Heidelberg,1983:555-636.

[3] 康新忠, 王寶元. 機(jī)械系統(tǒng)動(dòng)態(tài)優(yōu)化設(shè)計(jì)的靈敏度方法[J]. 機(jī)械工程學(xué)報(bào), 1990, 13(1): 18-23.

KANG Xin-zhong, WANG Bao-yuan. Sensitivity analysis method of the dynamic optimal design for mechanical systems[J]. Journal of Vibration Engineering,1990, 13(1): 18-23.

[4] Etman L F P, Van Campen D H, Schoofs A J G. Optimization of multibody systems using approximation concepts[C]//IUTAM.Symposium on Optimization of Mechanical Systems. Netherlands：Springer，1996：81-88.

[5] Li S, Petzold L. Software and algorithms for sensitivity analysis of large-scale differential-algebraic systems[J]. Journal of Computational and Applied Mathematics, 2000, 125:131-145.

[6] Maly T, Pctzold L R. Numerical methods and software for sensitivity analysis of differential-algebraic systems[J]. Applied Numerical Mathematics, 1996, 20: 57-59.

[7] Shabana A. An absolute nodal coordinates formulation for the large rotation and deformation analysis of flexible bodies[R]. University of Illionis at Chicago, 1996.

[8] Shabana A. Computational continuum mechanics[M]. New York: Cambridge University Press, 2011.

[9] Bonet J, Wood R. Nonlinear continuum mechanics for finite element analysis[M]. Cambridge: Cambridge University Press, 1997.

第一作者 王鐵成 男，博士生，1979年生

摘要：針對(duì)多個(gè)設(shè)計(jì)變量情況下的柔性多體系統(tǒng)靈敏度分析，采用伴隨變量法對(duì)基于絕對(duì)節(jié)點(diǎn)坐標(biāo)法建立的柔性多體系統(tǒng)進(jìn)行了研究。為了驗(yàn)證的該方法的計(jì)算效率，分別采用直接微分法和伴隨變量法對(duì)受重力作用的柔性單擺進(jìn)行了研究, 結(jié)果表明：這兩種方法計(jì)算結(jié)果的誤差很小，隨著設(shè)計(jì)變量數(shù)量的增加，伴隨變量法有更高的計(jì)算效率。

關(guān)鍵詞：多體系統(tǒng);絕對(duì)節(jié)點(diǎn)坐標(biāo)法;靈敏度;伴隨變量法

Sensitivity analysis of flexible multibody systems based on absolute nodal coordinate formulation

WANGTie-cheng1,CHENGuo-ping1,SUNDong-yang2(1. The State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China;2. College of Aerospace Engineering, Chongqing University, Chongqing 400044, China)

Abstract:For multiple design variables, an adjoint variable method was applied in sensitivity analysis of flexible multibody systems based on the absolute node coordinate formulation. In order to verify the computational efficiency of the method, a flexible pendulum under gravity was analysed by using the direct differentiation method and the adjoint variable method respectively. The results show that the errors of the two methods are both small, and the adjoint variable method has higher computational efficiency with the increase of design variables.

Key words:multibody system; sensitivity; absolute node coordinates formulation; adjoint variable method

中圖分類號(hào):O313.7

文獻(xiàn)標(biāo)志碼：A DOI:10.13465/j.cnki.jvs.2015.24.015

通信作者陳國(guó)平 男，教授，博士生導(dǎo)師，1956年生

收稿日期：2014-09-04修改稿收到日期：2014-11-26

基金項(xiàng)目：江蘇高校優(yōu)勢(shì)學(xué)科建設(shè)工程資助項(xiàng)目