三臂機器人在軌動(dòng)力學(xué)仿真

中國機械工程第27卷第3期2016年2月上半月臂機器人在軌動(dòng)力學(xué)仿真黃晉英1:2杜為民3李輝蔣志宏楊歡21.北京理工大學(xué),北京,10100812.中北大學(xué),太原,0300513.中國民用航空河南空管分局,鄭州,451162摘要:根據航天機器人太空作業(yè)的需要,設計了三臂機器人簡(jiǎn)化模型,應用DH方法建立了坐標變換矩陣,推導了該機器人的運動(dòng)學(xué)方程?；?ADAMS虛擬樣機技術(shù),分別在三種重力場(chǎng)環(huán)境下,對機器人的爬桿過(guò)程進(jìn)行仿真,得到了機械手攀爬過(guò)程受力曲線(xiàn),分析了影響機械手夾持力的主要因素,從而為這類(lèi)機器人的設計提供參考關(guān)鍵詞:機器人;在軌;仿真;夾持力中圖分類(lèi)號:TP241.2DOI:10.3969/1issn.1004-132X.2016.03.008Dynamic Simulation for a Three-arm Robot on OrbitHuang Jinying Du Weimin' Li Hui Jiang Zhihong Yang Huan1. Beijing Institute of Tecnology, Beijing, 1000812. North University of China, Taiyuan, 0300513. Henan Branch of Air Traffic Control of Civil Aviation of China, Zhengzhou, 451162Abstract: A simplified robot model with three arms was designed according to the needs of spaceoperations. The coordinate transformation matrix was established based on the D-H method and thekinematic equations of the robot were derived. Based on ADAMS virtual prototyping technology, theclimbing process of the robot on a pole was simulated for three kinds of gravity, respectively. Theforce curves and effect factors were obtained so as to offer the design references for the robotKey words: robot; on orbit; simulation; clamping force0引言其軌道形成2個(gè)移動(dòng)自由度,3條手臂共計18個(gè)太空環(huán)境極其惡劣,用機器人代替宇航員出自由度。工作時(shí),由3條手臂交替攀爬航天器上艙完成裝配維修和搬運等工作十分必要121?？赵O置的扶手到達指定位置,然后其中一條手臂固間環(huán)境與地面環(huán)境最大的區別是微重力,在微重定在扶手上,其余兩條手臂完成指定工作。這種力環(huán)境下,機器人處于漂浮狀態(tài),因此地面看似簡(jiǎn)三臂機構可使機器人具有多種靈活的運動(dòng)步態(tài)和單的工作將變得難以實(shí)現,其動(dòng)力學(xué)行為與地面操作方式以適應不同的工況,極大地滿(mǎn)足了空間有較大差別,其運動(dòng)學(xué)控制較為困難。多臂機器機器人太空中工作的要求人與雙臂機器人相比,多余的手臂可作為工作中小臂的基座,從而克服了上述難題。三臂機器人可在小臂關(guān)節機身結構上克服腿、臂獨立結構構型機器人本體機構質(zhì)量大、結構復雜等缺點(diǎn),適用于太空作業(yè),正日腰關(guān)節大臂關(guān)節益引起人們的關(guān)注34。然而,這類(lèi)空間機器人的肘關(guān)節大臂研究仍處于起步階段,本研究的目的就是通過(guò)機機械手器人的在軌仿真,分析機器人在軌運動(dòng)中的動(dòng)力機械手軌道學(xué)行為,為其設計和應用提供參考圖1三臂機器人簡(jiǎn)化幾何模型1三臂機器人的系統模型1.1機器人空間坐標系的建立按照D-H方法建立機器人系統坐標系,如圖空間三臂機器人的簡(jiǎn)化模型如圖1所示,其2所示。其中,O1XnY1ZB是航天器的坐標系本體結構為一等邊三角形,擁有3條完全相同的手臂。每條手臂由一組(2只)機械手、小臂和大dYH中國煤化工構件的動(dòng)坐標系,機器CNMHG陣用T(=1,2臂組成,其上的腰關(guān)節、大臂關(guān)節、肘關(guān)節和小臂n)表示,其中,是機器人關(guān)節序號,前置序號關(guān)節形成4個(gè)轉動(dòng)自由度,機械手的2只手掌與1表示目的坐標系的序號,-1=0表示機身坐標收稿日期:2014-08-15修回日期:2015-12-07系。對于本文所設計的機器人,n=18·328三臂機器人在軌動(dòng)力學(xué)仿真—黃晉英杜為民李輝等因此推導出三臂機器人每條臂的坐標變化矩陣為T(mén)=TBT1T2…T則機械臂末端在航天器坐標系下的位姿坐標可表示為(3)式中,P為航天器的位姿坐標。三臂機器人的仿真模型Z2.1機器人的簡(jiǎn)化幾何模型級在ProE軟件中建立三臂機器人的簡(jiǎn)化模型,直接將裝配好的模型導入 ADAMS軟件中進(jìn)行仿圖2三臂機器人空間坐標系真,其簡(jiǎn)化模型如圖1所示。該模型僅包含實(shí)現機1.2機器人運動(dòng)學(xué)方程的建立器人運動(dòng)所必需的機械構件,其余附件及電機等未機器人是由一系列具有空間彎曲軸線(xiàn)的廣義涉及,簡(jiǎn)化模型的主要幾何參數如圖4所示連桿連接在一起構成的。根據三臂機器人的實(shí)際參數,可確定該模型的D-H參數,推導出三臂機器人的運動(dòng)學(xué)方程方程。記機器人的任意兩中100個(gè)關(guān)節序號分別為i和i-1,如圖3所示,固連在該2個(gè)關(guān)節上的動(dòng)坐標系Z軸間的空間夾角(連桿扭角)為a,,空間最小距離(連桿長(cháng)度)為b;,2個(gè)關(guān)節公垂線(xiàn)間的距離為d2,2個(gè)關(guān)節X軸間的空間夾角(連桿夾角)為圖4機器人簡(jiǎn)化模型的主要幾何參數2.2仿真模型建立關(guān)節2.2.1模型質(zhì)量設從ProE中導入的幾何模型是沒(méi)有質(zhì)量的,關(guān)節-1在仿真前必須對各部件進(jìn)行質(zhì)量設置。本研究選用 ADAMS軟件材料庫中的鈦合金,并在重力圖3廣義連桿變換的4個(gè)特征參數D-H參數( Gravity)選項中將重力設置為無(wú)重力(模擬微重上述D-H參數可通過(guò)如下兩次旋轉和兩次力狀態(tài))或g/6(模擬月球環(huán)境重力)1),g為重平移坐標變換從坐標系OXYZ變換到坐標系力加速度O,-1X1Y-1Z:1810:①將坐標系OX,YZ;繞2.2.2運動(dòng)約束設置X1旋轉a:角,使Z:與Z1軸位于同一平面內在 ADAMS/View環(huán)境中,模型構件之間的②將坐標系OXYZ;沿X1軸平移b,將坐標相互運動(dòng)是靠運動(dòng)副來(lái)實(shí)現的,在腰關(guān)節、大臂關(guān)系OXYZ原點(diǎn)移到Z,軸線(xiàn)上;③將坐標系節、腕關(guān)節和小臂關(guān)節之間添加旋轉副,在機械手OXYZ;繞Z軸旋轉θ角,使X軸與X1軸平添加移動(dòng)副,其余構件添加固定副約束行;④將坐標系OXY,Z;沿Z1軸平移d2.2.3接觸定上述坐標變換可用如下坐標變換矩陣當兩個(gè)構件表面間發(fā)生接觸并存在相對運動(dòng)描述:H中國煤化工觸區產(chǎn)生摩擦力,摩擦cose: cosIna按數少CNM口而增大,在 ADAMS中二宵饑奮人攀爬仿真過(guò)程中,當ine sing cos0: sinaCosacosa機械手抓住航天器上的桿狀扶手時(shí),需要定義機械手與扶手之間的接觸。本研究選用 IMPACT中國機械工程第27卷第3期2016年2月上半月型接觸力和庫侖摩擦,其參數如表1所示。其中材料剛度、阻尼和摩擦因數由材料性質(zhì)所決定。表1仿真過(guò)程中接觸參數設置人人參數類(lèi)型參數值材料剛度1.00×10接觸力指數(b)阻尼(N·s/mm)穿透深度(mm)態(tài)摩擦因數動(dòng)態(tài)摩擦因數靜態(tài)滑移速度(mm/s)動(dòng)態(tài)滑移速度(mm/s)2.2.4運動(dòng)軌跡設計與驅動(dòng)添加d)對機器人模型施加驅動(dòng)就是讓其未受約束的圖6機器人仿真步態(tài)示意圖運動(dòng)副按照某種規律變化。根據機器人的幾何結成120°,桿置于1號機械手內,見(jiàn)圖6a。構特點(diǎn),在仿真過(guò)程中只需要添加旋轉驅動(dòng)與滑(2)1號手臂的兩只機械手沿桿軸線(xiàn)方向運移驅動(dòng),將旋轉驅動(dòng)添加到旋轉副和圓柱副,將滑動(dòng)(X向)閉合抓桿,見(jiàn)圖6b,運動(dòng)時(shí)間0.1s移驅動(dòng)添加到滑移副和圓柱副。每個(gè)運動(dòng)副需要(3)除1號手臂外,其余部件繞1號手臂腰關(guān)編寫(xiě)獨立的驅動(dòng)函數,以完成規定的運動(dòng)軌跡。節旋轉90°,機器人在 YBOBZB平面內運動(dòng),運動(dòng)在本研究中,仿真運動(dòng)軌跡規定為三臂輪流抓桿時(shí)間1s,見(jiàn)圖6c。向前移動(dòng),起始動(dòng)作為1號機械手抓桿,2號機械(4)2、3號手臂和機身繞1號手臂大臂關(guān)節在手開(kāi)始運動(dòng)抓桿;當2號機槭手抓穩杄后,1號機YBOεzε平面內旋轉60°,然后,2號手臂繞其大臂械手脫離桿,同時(shí)3號機械手重復前述動(dòng)作準備關(guān)節再旋轉60°,2號手臂到達抓桿前狀態(tài),運動(dòng)抓桿;如此往復,實(shí)現機器人沿攀爬桿向前運動(dòng)。時(shí)間共計4s,見(jiàn)圖6d在整個(gè)攀爬過(guò)程中,機器人始終與 YROBZE平面(5)2號手臂的兩只機械手閉合抓桿,其運動(dòng)平行,機械手的夾持力方向平行于XB軸。時(shí)間為1s;1、2號手臂同時(shí)抓桿時(shí)間持續1s,本文仿真中實(shí)現上述運動(dòng)的驅動(dòng)函數采用(6)1號手臂的兩只機械手向桿的反方向運p函數編寫(xiě)。step函數是 ADAMS函數庫中動(dòng),松開(kāi)爬桿,運動(dòng)時(shí)間1s。常用的函數,是三次多項式逼近階躍函數,由5個(gè)(7)1號手臂繞其大臂關(guān)節旋轉一60°,1、3號參數控制,分別為自變量(x)、自變量的初值(x0)手臂和機身繞2號大臂關(guān)節在 YBOBZE平面內旋和終值(x1),函數值的初值(h)和終值(h1)。在轉-120°,3號手臂繞其大臂旋轉60,3號機械手仿真中,從初值到終值的過(guò)程由軟件按照三次樣到達抓桿前位置,見(jiàn)圖6e和圖6條擬合自動(dòng)計算給出。step函數的具體格式為(8)3號手臂的兩只機槭手閉合抓桿,其運動(dòng)step(x,x0,h。,x1,h1)時(shí)間和驅動(dòng)函數同過(guò)程(2)中設置,兩只手同時(shí)抓圖5為由s數擬合得到的運動(dòng)軌跡曲桿時(shí)間持續1s線(xiàn)13。本研究中采用時(shí)間位移驅動(dòng),即要求各關(guān)(9)重復上述過(guò)程,實(shí)現機器人的整個(gè)攀爬節在規定的時(shí)間內運動(dòng)到規定的位置,自變量為過(guò)程時(shí)間坐標。圖6為爬桿過(guò)程步態(tài)圖不失一般性,本文僅對1號手臂抓兩次桿,23號手臂各抓一次桿進(jìn)行仿真。1號機械手驅動(dòng)函數為tep(time,0,0,0.1,-17.2)+step(time,9,0,10,17.2)+其余驅動(dòng)函數在此不再列出。為減小同時(shí)動(dòng)作帶來(lái)的振動(dòng)影響和前一次動(dòng)作對后一次動(dòng)作的圖5step函數擬合軌跡曲線(xiàn)中國煤化工聽(tīng)有運動(dòng)副的動(dòng)作時(shí)間基于上述驅動(dòng)函數,設計運動(dòng)過(guò)程如下CNMHG(1)初始位置設為1號手臂抓桿前的位置,3條手臂伸直與機身位于同一平面內,即 XBOBY3仿真結果及分析平面,1號手臂平行于YB軸,其余兩條手臂與之本文在無(wú)重力和不同方向g/6重力環(huán)境下三臂機器人在軌動(dòng)力學(xué)仿真—黃晉英杜為民李輝等分別對三臂機器人的爬桿過(guò)程進(jìn)行動(dòng)力學(xué)仿真,觸,在所允許的切入深度下緊握爬桿,但實(shí)際中會(huì )獲得不同環(huán)境下機械手受力情況和大臂關(guān)節運動(dòng)出現圖10b~圖10d所示的情況。圖10b兩側手過(guò)程扭矩曲線(xiàn),為機器人的設計提供參考掌切人深度不一致(左側大于右側),圖10c和圖3.1機械手受力情況分析l0d接觸面變?yōu)椴糠纸佑|,這種接觸方式之間的圖7~圖9分別為三臂機器人在無(wú)重力、g/6變換導致接觸力發(fā)生波動(dòng)重力環(huán)境下3條手臂的機械手受力曲線(xiàn),其中,CONTACT_1.FX為1號手臂的機械手夾持力桿桿(XB向)曲線(xiàn), CONTACT5FX為2號手臂的機械手夾持力(X1向)曲線(xiàn), CONTACT8.FX為3號手臂的機械手夾持力(Xε向)曲線(xiàn)。重心位于機器人幾何中心-CONTACT-1FX桿CONTACT-8FX(c)圖10機械手抓桿接觸狀例如,2~3s時(shí),1號機械手抓緊爬桿,機身時(shí)間t/s繞Y軸旋轉,轉動(dòng)過(guò)程對機械手的受力影響較圖7無(wú)重力環(huán)境下3條手臂機械手受力情況大,造成夾持力發(fā)生波動(dòng),如圖11所示。一側手CONTACT-1FX掌夾持力增大的同時(shí)另一側手掌夾持力減小,為CONTACT-8FX防止機械手脫離爬桿,控制兩側手掌的接觸力在某一數值(如-4000N)上下波動(dòng)。圖11中2條曲線(xiàn)分別為1號機械手在無(wú)重力情況下兩只手掌15的夾持力曲線(xiàn),夾持力峰值最大達到-5700N,波動(dòng)幅度與平均夾持力相比增大30%。11~17時(shí)間t/s時(shí)1號手臂脫離爬桿,3號手臂向抓桿位置擺動(dòng),圖8重力平行于Xs軸環(huán)境下夾持力峰值在12~14s期間波動(dòng)較大,達到3條手臂機械手受力情況5785N,如圖12所示。機體轉動(dòng)慣性會(huì )使機體CONTACT-1FX發(fā)生晃動(dòng),嚴重的晃動(dòng)不僅造成夾持力波動(dòng),若控CONTACT-8FX制不當甚至會(huì )使機械手滑落爬桿,因此,應進(jìn)一步考察機身轉動(dòng)速度對夾持力的影響情況左側手掌--右側手掌時(shí)間t/s圖9重力平行于一Y軸環(huán)境下3條手臂機械手受力情況時(shí)間t/s由上述仿真結果可見(jiàn),雖然1號手臂首次抓圖11無(wú)重力情況下1號機械手夾持力桿持續時(shí)間遠遠短于其他兩條手臂和1號手臂第次抓桿持續時(shí)間,但3條手臂的機械手在抓桿圖13所示為三種重力情況下2號機械手夾過(guò)程中受力基本一致。由此可知,抓桿過(guò)程的速對中:圖中曲線(xiàn)1為無(wú)重力情況下仿真中國煤化工度變化基本不影響夾持力的大小。在整個(gè)爬桿過(guò)重力方向為XB向和CNMH程中,機體運動(dòng)會(huì )引起機械手與桿之間的相對運脂出的是,重力對夾持動(dòng),如圖10所示。圖10a所示為理論抓桿狀態(tài),力有較大影響。當重力方向與夾持力方向一致(即重力方向為Xs向)時(shí),重力在造成夾持力增機械手兩側手掌和軌道與桿的3個(gè)平面完全接大的同時(shí),還會(huì )使得兩只機械手夾持力不相等,在331中國機械工程第27卷第3期2016年2月上半月左側手掌3右側手掌2000(12.4,-2387.4)CONTACT-3.FY- CONTACT-6F600(9.2,-5048.2)(13.1,-5785.1)CONTACT-9.FY68101214161820時(shí)間t/圖12無(wú)重力情況下2號機械手夾持力圖15無(wú)重力環(huán)境下3條手臂機械手軌道受力情況快速抓桿時(shí)出現多次碰撞(圖10),導致抓桿困CONTACT-3.FY難。當重力方向與夾持力方向垂直(即重力方向CONTACT-6F為—YB向)時(shí),重力使得夾持力波動(dòng)幅度和頻率CONTACT-9F減小。以上情況表明,無(wú)重力狀態(tài)的運動(dòng)姿態(tài)較重力場(chǎng)中的運動(dòng)姿態(tài)更難控制無(wú)重力重力方向為XB時(shí)間t/s重力方向為一YB(18.1,20.5)圖16重力平行于XB軸環(huán)境下(12.3,-2377.2)3條手臂機械手軌道受力情況z-2000(10.3,(136-2377.2)M一6000(8.x03.0)·(13.1.-6739.8)4F— CONTACT-3FY8000}(8.0.,-7143.3)(13.6.8576.0)皋-6-c0NmA:6F100006810121416182022時(shí)間t/s時(shí)間t/s圖17重力平行于一YB軸環(huán)境下圖13三種重力情況下2號機械手夾持力3條手臂機械手軌道受力情況3.2機械手軌道受力情況分析由上述曲線(xiàn)看出,在機械手抓緊爬桿后手掌對于上述設計步態(tài),在機械手與爬桿接觸之和桿接觸面之間的相對運動(dòng)比較頻繁,導致軌道前,機被手軌道的一端首先與桿發(fā)生碰撞;此后逐接觸力頻繁波動(dòng)。從碰撞力來(lái)看,重力平行于漸碰撞到達指定部位(圖14中從右至左)。此時(shí)x軸環(huán)境下碰撞力最大。表2為各種環(huán)境下機軌道平面與桿平面接觸或平行,該過(guò)程會(huì )給大臂械手及其軌道最大受力匯總表關(guān)節形成沖擊,造成大臂關(guān)節逆向轉動(dòng)。如果驅表2三種環(huán)境下機械手受力動(dòng)力矩不夠大,則機械手將無(wú)法到達抓桿前指定機械手臂重力平行重力平行無(wú)重力位置。圖15~圖17為三種環(huán)境下軌道的受力情編號于XB軸于Yn軸況,其中 CONTACT6.FY為1號機械手軌道的4006100~7500機械手2號4000受力曲線(xiàn), CONTACT3.FY為2號機械手軌道的受力曲線(xiàn)(-YB向), CONTACT9.FY為3號1990機械手4945機械手軌道的受力曲線(xiàn),2號機械手軌道方向與5000軌道其他二只機械手軌道的坐標方向相反4結論(1)對于本文設計的步態(tài),為了驗證抓桿速度和多關(guān)節同步運動(dòng)對機械手夾持力的影響,對1TV凵中國煤化工別采用不同抓桿速度CNMHG桿速度對夾持力的影不腰部叫轉時(shí)采用多關(guān)節聯(lián)動(dòng)、其他關(guān)節單獨運動(dòng)的方式進(jìn)行仿真分析,結果表明,多圖14機械手軌道與桿碰撞過(guò)程關(guān)節聯(lián)動(dòng)會(huì )造成機器人振動(dòng)加劇,夾持力增大,夾·332三臂機器人在軌動(dòng)力學(xué)仿真—黃晉英杜為民李輝等持力波動(dòng)加大,故在設計機器人控制策略時(shí)應盡ing Its Manipulation[J]. Journal of Beijing Univers量避免多關(guān)節聯(lián)動(dòng)of Posts and Telecommunications. 2002. 25(1)(2)重力對機器人的姿態(tài)影響較大,在設計重34-38,力環(huán)境下機器人攀爬路線(xiàn)時(shí),應盡量讓重力和夾「6]葉平,孫漢旭.具有冗余度的三分支空間機器人逆持力與桿長(cháng)度方向垂直(重力為一YB向),這樣可動(dòng)學(xué)分析冂].機械工程學(xué)報,2005,41(11):58使姿態(tài)容易保持穩定。Ye Ping, Sun Hanxu. Analysis of Inverse Kinema(3)應改進(jìn)設計機械手軌道,盡量減少其邊緣for Three-branch Space Robot with Redundancy對桿的碰撞;同時(shí)改變關(guān)節運動(dòng)順序,將邊緣首先LJ. Chinese Journal of Mechanical Engineering碰撞改為整個(gè)軌道平面與桿全接觸,這樣可減小2005,41(11):58-62碰撞力引起的機器人振動(dòng)。L7] Suga T G, Kumar V. Control of Cooperating Mo(4)在仿真過(guò)程中,仿真參數設置對仿真結果bile Manipulators[J]. IEEE Transaction on Robot影響較大,本文所設置的參數主要以經(jīng)驗值為依and Automation, 2002, 18(1):94-103據,仿真結果可為這類(lèi)空間機器人的設計與控制81方深瑋.基于 ADAMS機器人動(dòng)力學(xué)仿真研究D)提供定性參考北京:北京郵電大學(xué),2009[9 Zhuang Y F, Liu DQ, Wang J G. Dynamic Modeling and Analyzing of a Walking Robot [J].The參考文獻Journal of China Universities of Posts and Telecom[1 Rastegari R, Moosavian S AA Multiple Impedancemunications,2014,21(1):122-128Control of Space Free- flying Robots via VirtualLinkages[J]. Acta astronautica,2010,6:748759[10]甘屹,王均壘,孫福佳.基于給定工作空間的6R2]梁捷,陳力.雙臂空間機器人基于高斯型函數的姿型機器人DH參數優(yōu)化設計[].中國機械工程2014,25(22):3003-3007,3011態(tài)、關(guān)節運動(dòng)模糊自適應補償控制[J.中國機械工Gang Yi. Wang Junlei, Sun Fujia. Optimal Design程,2010,21(3):330-336of dh parameters of a 6r robot for a prescribedLiang Jie, Chen Li. Fuzzy Adaptive CompensatioControl for Space-based Robot System with DualWorkspace [J]. China Mecha2014,25(22):3003-3007,3011arm Based on Gaussian- type Function to Trac[1]姬廣勛.履帶式載人月球車(chē)移動(dòng)系統的設計與運rectory in Joint Space[J. China Mecha動(dòng)性能分析[D].哈爾濱:哈爾濱工業(yè)大學(xué),20gineering,2010,21(3):330-336[3]慧記莊,魏芳勝,高凱,等.基于 ADAMs的冗余驅[12] Sharf1, Zhang YN. A Contact Force Solution forNon-colliding Contact Dynamics Simulation [J].動(dòng)并聯(lián)機器人動(dòng)力學(xué)仿真研究[].工程設計學(xué)報Multibody System Dynamics, 2006, 16(3): 2632012,19(5):362-365Hui Jizhuang, Wei Fangsheng, Gao kai,etal. Study[13]阮鵬,俞志偉,張昊,等基于 ADAMS的仿壁虎n Dynamics Simulation of Redundantly Actuated機器人步態(tài)規劃及仿真[J].機器人,2010,32(4)Parallel Robot Based on ADAMS[J. Chinese Jor499-509nal of Engineering Design, 2012, 19(5): 362-365[4]孫漢旭,賈慶軒,張秋豪,等.基于三分支機器人關(guān)Zhiwei, Zhang Ha節空間軌跡規劃的研究[].北京郵電大學(xué)學(xué)報,Planning and Simulation of Gecko Inspired RobotBased on ADAMS[J. Robot, 2010, 32(4):4992006,29(3):81-85Sun Hanxu, Jia Qingxuan, Zhang Qiuhao, et al.(編輯蘇衛國)Trajectory Planning in Joint Space Based on TriBranch Robot[J]. Journal of Beijing University of作者簡(jiǎn)介:黃晉英,女,1971年生。北京理工大學(xué)宇航學(xué)院博士Posts and Telecommunications, 2006, 29(3):81后研究人員,中北大學(xué)機械與動(dòng)力工程學(xué)院教授。主要研究方向為動(dòng)力學(xué)分析與控制。杜為民,女,1966年生。中國民用航空河[5]孫漢旭,姬羅栓,徐揚生.空間三臂機器人在操作中南空管分局高級工程師。李輝,男,1982年生。北京理工大學(xué)的平衡問(wèn)題[.北京郵電大學(xué)學(xué)報,2002,25(1)中國煤化工生。北京理工大學(xué)機電學(xué)院Sun Hanxu, Ji Luoshuang, Xu Yangsheng Study oCNM云北大學(xué)機械與動(dòng)力工程學(xué)院Balance Problem for a Three-arm Space Robot dur-333