• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    四元數(shù)與歐拉角剛體動力學(xué)數(shù)值積分算法及其比較

    2014-03-13 08:43徐小明鐘萬勰
    計算機(jī)輔助工程 2014年1期

    徐小明 鐘萬勰

    摘要:為推廣四元數(shù)保辛積分在工程中的應(yīng)用,對歐拉角表示的狀態(tài)方程數(shù)值積分與四元數(shù)的保辛積分進(jìn)行比較.重陀螺的數(shù)值仿真結(jié)果表明四元數(shù)保辛積分的數(shù)值結(jié)果明顯優(yōu)于歐拉角狀態(tài)方程積分.與歐拉角狀態(tài)方程積分相比,四元數(shù)保辛積分在剛體動力學(xué)的數(shù)值仿真中更具優(yōu)勢.

    關(guān)鍵詞:四元數(shù); 歐拉角; 剛體動力學(xué); 保辛; 重陀螺

    中圖分類號: TP391.9; O313.3

    文獻(xiàn)標(biāo)志碼: A

    0引言

    四元數(shù)、歐拉角和方向余弦[1]是描述剛體旋轉(zhuǎn)最主要的3種坐標(biāo)形式.方向余弦法需要9個參量,應(yīng)用較少;而另外2種則應(yīng)用廣泛,如航天器姿態(tài)控制和捷聯(lián)式慣性導(dǎo)航系統(tǒng)[1]等,對于兩者的研究也卷帙浩繁,但對兩者優(yōu)劣的評價卻褒貶不一.

    描述剛體在三維空間中的運動姿態(tài)可采用2類12種歐拉角系統(tǒng),分別對應(yīng)于不同的旋轉(zhuǎn)軸先后次序.目前公認(rèn)的用歐拉角描述旋轉(zhuǎn)的固有缺陷是奇異性問題[2],即:無論采用哪種歐拉角系統(tǒng),都不可避免地會含有奇異點,使得在該點附近區(qū)域進(jìn)行的數(shù)值積分精度不高.對于小角度旋轉(zhuǎn),只要采用適當(dāng)?shù)臍W拉角系統(tǒng)便可避開奇異點;而在大角度旋轉(zhuǎn)時,若想避開奇異點,必須提供2套歐拉角系統(tǒng)交替進(jìn)行計算.據(jù)此,黃雪樵[3]提出一種“雙歐法”的計算方法;近幾年仍有學(xué)者[4]在繼續(xù)研究該方法.雙歐法雖然解決奇異性問題,但是計算過程較復(fù)雜.

    四元數(shù)用于描述剛體旋轉(zhuǎn),沒有奇異性問題,可很好地描述剛體的全角度旋轉(zhuǎn).然而,四元數(shù)需要滿足長度等于1的單位約束,這成為制約其應(yīng)用的限制.在實際應(yīng)用中,經(jīng)常采用的正交化修正等方法只能緩解長度的偏移,無法從根本上解決問題;黃雪樵[3]在其雙歐法中也有所討論.目前,對于單位約束最佳的解決方案是將四元數(shù)表示的剛體動力學(xué)方程導(dǎo)入微分代數(shù)方程范疇,近年來逐漸有學(xué)者[57]展開相關(guān)問題的研究.徐小明等[8]提出一種基于四元數(shù)理論描述剛體旋轉(zhuǎn)的保辛數(shù)值積分方法.該方法先將問題導(dǎo)入微分代數(shù)方程系統(tǒng),然后利用分析結(jié)構(gòu)力學(xué)理論[9]進(jìn)行逐步積分,該積分保辛.該方法在積分點上嚴(yán)格滿足四元數(shù)長度等于1的約束條件,而在區(qū)段內(nèi)部則由插值近似,屬于祖沖之類方法[10].數(shù)值算例表明仿真效果優(yōu)異.

    本文簡要介紹四元數(shù)和歐拉角的基本理論,以重陀螺為例對2種表示形式的數(shù)值積分進(jìn)行比較.對于歐拉角表述,應(yīng)用比較普遍的狀態(tài)方程表述.研究表明,以歐拉角和角速度為狀態(tài)變量形成的1階微分方程在使用差分近似積分時,精度損失很大,能量不守恒;該現(xiàn)象被周江華等[11]稱為“睡陀螺”.這是一個值得注意的問題,卻未得到廣泛關(guān)注;而采用文獻(xiàn)[8]給出的保辛格式,四元數(shù)單位長度約束條件得到滿足,仿真結(jié)果優(yōu)異,能量也達(dá)到守恒.

    1剛體旋轉(zhuǎn)及其運動學(xué)表示

    剛體運動由質(zhì)心平動和繞質(zhì)心轉(zhuǎn)動組成.如果剛體上有一固定點,則稱為剛體定點轉(zhuǎn)動問題.假設(shè)Oxyz為系統(tǒng)的慣性坐標(biāo)系,O為原點,亦為固定點.Ox′y′z′為隨體坐標(biāo)系,固定于剛體上.若將剛體的隨體坐標(biāo)軸取為與固定點有關(guān)的主軸,則剛體定點轉(zhuǎn)動可由式(1)描述.

    2種數(shù)值積分的系統(tǒng)能量隨時間變化情況見圖4.在歐拉角表示中,雖然采用辛歐拉格式進(jìn)行數(shù)值積分,但是能量卻保持得不好,這驗證對狀態(tài)方程應(yīng)用辛歐拉格式并不能保辛.與之相反,四元數(shù)的數(shù)值積分保辛,其能量保持得很好,這也是保辛積分的優(yōu)勢所在.四元數(shù)保辛積分的約束誤差見圖5,表明在時間積分過程中單位長度約束條件滿足得很好.

    4結(jié)束語

    介紹2種剛體旋轉(zhuǎn)的數(shù)值積分,一種基于歐拉角表示,另一種基于四元數(shù)表示.以重陀螺的高速旋轉(zhuǎn)為例,對2種數(shù)值積分進(jìn)行比較發(fā)現(xiàn):以歐拉角、角速度組成狀態(tài)變量,然后直接使用辛歐拉格式不

    能保辛,能量衰減很快,數(shù)值積分存在缺陷;與之相

    反,采用四元數(shù)表示,根據(jù)分析結(jié)構(gòu)力學(xué)的保辛積分方法,并結(jié)合祖沖之類方法的思想,可以很好地避免約束違約,仿真效果令人滿意,可作為陀螺等仿真分析的有力工具.

    本文僅對以歐拉角、角速度組成狀態(tài)變量的數(shù)值積分進(jìn)行研究,對其他形式并未涉及,對其積分效果不佳的成因亦未研究,還有很多工作有待展開.

    參考文獻(xiàn):

    [1]張樹俠, 孫靜. 捷聯(lián)式慣性導(dǎo)航系統(tǒng)[M]. 北京: 國防工業(yè)出版社, 1992: 4880.

    [2]趙曉穎, 溫立書, 么彩蓮. 歐拉角參數(shù)表示下姿態(tài)的2階運動奇異性[J]. 科學(xué)技術(shù)與工程, 2012, 12(3): 634637.

    [3]黃雪樵. 克服歐拉方程奇異性的雙歐法[J]. 飛行力學(xué), 1994, 12(4): 2837.

    [4]李躍軍, 閻超. 飛行器姿態(tài)角解算的全角度雙歐法[J]. 北京航空航天大學(xué)學(xué)報, 2007, 33(5): 505508.

    [5]NIKRAVESH P E, WEHAGE R A, KWON K. Euler parameters in computational kinematics and dynamics: Part 1[J]. J Mechanisms, Transmissions & Automation Des, 1985, 107(3): 358365.

    [6]BETSCH P, SIEBERT R. Rigid body dynamics in terms of quaternions: Hamiltonian formulation and conserving numerical integration[J]. Int J Numer Methods Eng, 2009, 79(4): 444473.

    [7]UDWADIA F E, SCHUTTE A D. An alternative derivation of the quaternion equations of motion for rigidbody rotational dynamics[J]. J Appl Mech, 2010, 77(4): 44505.

    [8]徐小明, 鐘萬勰. 剛體動力學(xué)的四元數(shù)表示及保辛積分[J]. 應(yīng)用數(shù)學(xué)和力學(xué), 2014, 35(1): 111.

    [9]鐘萬勰, 高強(qiáng). 約束動力系統(tǒng)的分析結(jié)構(gòu)力學(xué)積分[J]. 動力學(xué)與控制學(xué)報, 2006, 4(3): 193200.

    [10]鐘萬勰, 高強(qiáng), 彭海軍. 經(jīng)典力學(xué)辛講[M]. 大連: 大連理工大學(xué)出版社, 2013: 202241.

    [11]周江華, 苗育紅, 李宏, 等. 四元數(shù)在剛體姿態(tài)仿真中的應(yīng)用研究[J].飛行力學(xué), 2000, 18(4): 2832.

    [12]HAIRER E, LUBICH C, WANNER G. Geometric numerical integration: structurepreserving algorithms for ordinary differential equations[M]. Berlin: Springer, 2006: 189.

    (編輯于杰)

    摘要:為推廣四元數(shù)保辛積分在工程中的應(yīng)用,對歐拉角表示的狀態(tài)方程數(shù)值積分與四元數(shù)的保辛積分進(jìn)行比較.重陀螺的數(shù)值仿真結(jié)果表明四元數(shù)保辛積分的數(shù)值結(jié)果明顯優(yōu)于歐拉角狀態(tài)方程積分.與歐拉角狀態(tài)方程積分相比,四元數(shù)保辛積分在剛體動力學(xué)的數(shù)值仿真中更具優(yōu)勢.

    關(guān)鍵詞:四元數(shù); 歐拉角; 剛體動力學(xué); 保辛; 重陀螺

    中圖分類號: TP391.9; O313.3

    文獻(xiàn)標(biāo)志碼: A

    0引言

    四元數(shù)、歐拉角和方向余弦[1]是描述剛體旋轉(zhuǎn)最主要的3種坐標(biāo)形式.方向余弦法需要9個參量,應(yīng)用較少;而另外2種則應(yīng)用廣泛,如航天器姿態(tài)控制和捷聯(lián)式慣性導(dǎo)航系統(tǒng)[1]等,對于兩者的研究也卷帙浩繁,但對兩者優(yōu)劣的評價卻褒貶不一.

    描述剛體在三維空間中的運動姿態(tài)可采用2類12種歐拉角系統(tǒng),分別對應(yīng)于不同的旋轉(zhuǎn)軸先后次序.目前公認(rèn)的用歐拉角描述旋轉(zhuǎn)的固有缺陷是奇異性問題[2],即:無論采用哪種歐拉角系統(tǒng),都不可避免地會含有奇異點,使得在該點附近區(qū)域進(jìn)行的數(shù)值積分精度不高.對于小角度旋轉(zhuǎn),只要采用適當(dāng)?shù)臍W拉角系統(tǒng)便可避開奇異點;而在大角度旋轉(zhuǎn)時,若想避開奇異點,必須提供2套歐拉角系統(tǒng)交替進(jìn)行計算.據(jù)此,黃雪樵[3]提出一種“雙歐法”的計算方法;近幾年仍有學(xué)者[4]在繼續(xù)研究該方法.雙歐法雖然解決奇異性問題,但是計算過程較復(fù)雜.

    四元數(shù)用于描述剛體旋轉(zhuǎn),沒有奇異性問題,可很好地描述剛體的全角度旋轉(zhuǎn).然而,四元數(shù)需要滿足長度等于1的單位約束,這成為制約其應(yīng)用的限制.在實際應(yīng)用中,經(jīng)常采用的正交化修正等方法只能緩解長度的偏移,無法從根本上解決問題;黃雪樵[3]在其雙歐法中也有所討論.目前,對于單位約束最佳的解決方案是將四元數(shù)表示的剛體動力學(xué)方程導(dǎo)入微分代數(shù)方程范疇,近年來逐漸有學(xué)者[57]展開相關(guān)問題的研究.徐小明等[8]提出一種基于四元數(shù)理論描述剛體旋轉(zhuǎn)的保辛數(shù)值積分方法.該方法先將問題導(dǎo)入微分代數(shù)方程系統(tǒng),然后利用分析結(jié)構(gòu)力學(xué)理論[9]進(jìn)行逐步積分,該積分保辛.該方法在積分點上嚴(yán)格滿足四元數(shù)長度等于1的約束條件,而在區(qū)段內(nèi)部則由插值近似,屬于祖沖之類方法[10].數(shù)值算例表明仿真效果優(yōu)異.

    本文簡要介紹四元數(shù)和歐拉角的基本理論,以重陀螺為例對2種表示形式的數(shù)值積分進(jìn)行比較.對于歐拉角表述,應(yīng)用比較普遍的狀態(tài)方程表述.研究表明,以歐拉角和角速度為狀態(tài)變量形成的1階微分方程在使用差分近似積分時,精度損失很大,能量不守恒;該現(xiàn)象被周江華等[11]稱為“睡陀螺”.這是一個值得注意的問題,卻未得到廣泛關(guān)注;而采用文獻(xiàn)[8]給出的保辛格式,四元數(shù)單位長度約束條件得到滿足,仿真結(jié)果優(yōu)異,能量也達(dá)到守恒.

    1剛體旋轉(zhuǎn)及其運動學(xué)表示

    剛體運動由質(zhì)心平動和繞質(zhì)心轉(zhuǎn)動組成.如果剛體上有一固定點,則稱為剛體定點轉(zhuǎn)動問題.假設(shè)Oxyz為系統(tǒng)的慣性坐標(biāo)系,O為原點,亦為固定點.Ox′y′z′為隨體坐標(biāo)系,固定于剛體上.若將剛體的隨體坐標(biāo)軸取為與固定點有關(guān)的主軸,則剛體定點轉(zhuǎn)動可由式(1)描述.

    2種數(shù)值積分的系統(tǒng)能量隨時間變化情況見圖4.在歐拉角表示中,雖然采用辛歐拉格式進(jìn)行數(shù)值積分,但是能量卻保持得不好,這驗證對狀態(tài)方程應(yīng)用辛歐拉格式并不能保辛.與之相反,四元數(shù)的數(shù)值積分保辛,其能量保持得很好,這也是保辛積分的優(yōu)勢所在.四元數(shù)保辛積分的約束誤差見圖5,表明在時間積分過程中單位長度約束條件滿足得很好.

    4結(jié)束語

    介紹2種剛體旋轉(zhuǎn)的數(shù)值積分,一種基于歐拉角表示,另一種基于四元數(shù)表示.以重陀螺的高速旋轉(zhuǎn)為例,對2種數(shù)值積分進(jìn)行比較發(fā)現(xiàn):以歐拉角、角速度組成狀態(tài)變量,然后直接使用辛歐拉格式不

    能保辛,能量衰減很快,數(shù)值積分存在缺陷;與之相

    反,采用四元數(shù)表示,根據(jù)分析結(jié)構(gòu)力學(xué)的保辛積分方法,并結(jié)合祖沖之類方法的思想,可以很好地避免約束違約,仿真效果令人滿意,可作為陀螺等仿真分析的有力工具.

    本文僅對以歐拉角、角速度組成狀態(tài)變量的數(shù)值積分進(jìn)行研究,對其他形式并未涉及,對其積分效果不佳的成因亦未研究,還有很多工作有待展開.

    參考文獻(xiàn):

    [1]張樹俠, 孫靜. 捷聯(lián)式慣性導(dǎo)航系統(tǒng)[M]. 北京: 國防工業(yè)出版社, 1992: 4880.

    [2]趙曉穎, 溫立書, 么彩蓮. 歐拉角參數(shù)表示下姿態(tài)的2階運動奇異性[J]. 科學(xué)技術(shù)與工程, 2012, 12(3): 634637.

    [3]黃雪樵. 克服歐拉方程奇異性的雙歐法[J]. 飛行力學(xué), 1994, 12(4): 2837.

    [4]李躍軍, 閻超. 飛行器姿態(tài)角解算的全角度雙歐法[J]. 北京航空航天大學(xué)學(xué)報, 2007, 33(5): 505508.

    [5]NIKRAVESH P E, WEHAGE R A, KWON K. Euler parameters in computational kinematics and dynamics: Part 1[J]. J Mechanisms, Transmissions & Automation Des, 1985, 107(3): 358365.

    [6]BETSCH P, SIEBERT R. Rigid body dynamics in terms of quaternions: Hamiltonian formulation and conserving numerical integration[J]. Int J Numer Methods Eng, 2009, 79(4): 444473.

    [7]UDWADIA F E, SCHUTTE A D. An alternative derivation of the quaternion equations of motion for rigidbody rotational dynamics[J]. J Appl Mech, 2010, 77(4): 44505.

    [8]徐小明, 鐘萬勰. 剛體動力學(xué)的四元數(shù)表示及保辛積分[J]. 應(yīng)用數(shù)學(xué)和力學(xué), 2014, 35(1): 111.

    [9]鐘萬勰, 高強(qiáng). 約束動力系統(tǒng)的分析結(jié)構(gòu)力學(xué)積分[J]. 動力學(xué)與控制學(xué)報, 2006, 4(3): 193200.

    [10]鐘萬勰, 高強(qiáng), 彭海軍. 經(jīng)典力學(xué)辛講[M]. 大連: 大連理工大學(xué)出版社, 2013: 202241.

    [11]周江華, 苗育紅, 李宏, 等. 四元數(shù)在剛體姿態(tài)仿真中的應(yīng)用研究[J].飛行力學(xué), 2000, 18(4): 2832.

    [12]HAIRER E, LUBICH C, WANNER G. Geometric numerical integration: structurepreserving algorithms for ordinary differential equations[M]. Berlin: Springer, 2006: 189.

    (編輯于杰)

    摘要:為推廣四元數(shù)保辛積分在工程中的應(yīng)用,對歐拉角表示的狀態(tài)方程數(shù)值積分與四元數(shù)的保辛積分進(jìn)行比較.重陀螺的數(shù)值仿真結(jié)果表明四元數(shù)保辛積分的數(shù)值結(jié)果明顯優(yōu)于歐拉角狀態(tài)方程積分.與歐拉角狀態(tài)方程積分相比,四元數(shù)保辛積分在剛體動力學(xué)的數(shù)值仿真中更具優(yōu)勢.

    關(guān)鍵詞:四元數(shù); 歐拉角; 剛體動力學(xué); 保辛; 重陀螺

    中圖分類號: TP391.9; O313.3

    文獻(xiàn)標(biāo)志碼: A

    0引言

    四元數(shù)、歐拉角和方向余弦[1]是描述剛體旋轉(zhuǎn)最主要的3種坐標(biāo)形式.方向余弦法需要9個參量,應(yīng)用較少;而另外2種則應(yīng)用廣泛,如航天器姿態(tài)控制和捷聯(lián)式慣性導(dǎo)航系統(tǒng)[1]等,對于兩者的研究也卷帙浩繁,但對兩者優(yōu)劣的評價卻褒貶不一.

    描述剛體在三維空間中的運動姿態(tài)可采用2類12種歐拉角系統(tǒng),分別對應(yīng)于不同的旋轉(zhuǎn)軸先后次序.目前公認(rèn)的用歐拉角描述旋轉(zhuǎn)的固有缺陷是奇異性問題[2],即:無論采用哪種歐拉角系統(tǒng),都不可避免地會含有奇異點,使得在該點附近區(qū)域進(jìn)行的數(shù)值積分精度不高.對于小角度旋轉(zhuǎn),只要采用適當(dāng)?shù)臍W拉角系統(tǒng)便可避開奇異點;而在大角度旋轉(zhuǎn)時,若想避開奇異點,必須提供2套歐拉角系統(tǒng)交替進(jìn)行計算.據(jù)此,黃雪樵[3]提出一種“雙歐法”的計算方法;近幾年仍有學(xué)者[4]在繼續(xù)研究該方法.雙歐法雖然解決奇異性問題,但是計算過程較復(fù)雜.

    四元數(shù)用于描述剛體旋轉(zhuǎn),沒有奇異性問題,可很好地描述剛體的全角度旋轉(zhuǎn).然而,四元數(shù)需要滿足長度等于1的單位約束,這成為制約其應(yīng)用的限制.在實際應(yīng)用中,經(jīng)常采用的正交化修正等方法只能緩解長度的偏移,無法從根本上解決問題;黃雪樵[3]在其雙歐法中也有所討論.目前,對于單位約束最佳的解決方案是將四元數(shù)表示的剛體動力學(xué)方程導(dǎo)入微分代數(shù)方程范疇,近年來逐漸有學(xué)者[57]展開相關(guān)問題的研究.徐小明等[8]提出一種基于四元數(shù)理論描述剛體旋轉(zhuǎn)的保辛數(shù)值積分方法.該方法先將問題導(dǎo)入微分代數(shù)方程系統(tǒng),然后利用分析結(jié)構(gòu)力學(xué)理論[9]進(jìn)行逐步積分,該積分保辛.該方法在積分點上嚴(yán)格滿足四元數(shù)長度等于1的約束條件,而在區(qū)段內(nèi)部則由插值近似,屬于祖沖之類方法[10].數(shù)值算例表明仿真效果優(yōu)異.

    本文簡要介紹四元數(shù)和歐拉角的基本理論,以重陀螺為例對2種表示形式的數(shù)值積分進(jìn)行比較.對于歐拉角表述,應(yīng)用比較普遍的狀態(tài)方程表述.研究表明,以歐拉角和角速度為狀態(tài)變量形成的1階微分方程在使用差分近似積分時,精度損失很大,能量不守恒;該現(xiàn)象被周江華等[11]稱為“睡陀螺”.這是一個值得注意的問題,卻未得到廣泛關(guān)注;而采用文獻(xiàn)[8]給出的保辛格式,四元數(shù)單位長度約束條件得到滿足,仿真結(jié)果優(yōu)異,能量也達(dá)到守恒.

    1剛體旋轉(zhuǎn)及其運動學(xué)表示

    剛體運動由質(zhì)心平動和繞質(zhì)心轉(zhuǎn)動組成.如果剛體上有一固定點,則稱為剛體定點轉(zhuǎn)動問題.假設(shè)Oxyz為系統(tǒng)的慣性坐標(biāo)系,O為原點,亦為固定點.Ox′y′z′為隨體坐標(biāo)系,固定于剛體上.若將剛體的隨體坐標(biāo)軸取為與固定點有關(guān)的主軸,則剛體定點轉(zhuǎn)動可由式(1)描述.

    2種數(shù)值積分的系統(tǒng)能量隨時間變化情況見圖4.在歐拉角表示中,雖然采用辛歐拉格式進(jìn)行數(shù)值積分,但是能量卻保持得不好,這驗證對狀態(tài)方程應(yīng)用辛歐拉格式并不能保辛.與之相反,四元數(shù)的數(shù)值積分保辛,其能量保持得很好,這也是保辛積分的優(yōu)勢所在.四元數(shù)保辛積分的約束誤差見圖5,表明在時間積分過程中單位長度約束條件滿足得很好.

    4結(jié)束語

    介紹2種剛體旋轉(zhuǎn)的數(shù)值積分,一種基于歐拉角表示,另一種基于四元數(shù)表示.以重陀螺的高速旋轉(zhuǎn)為例,對2種數(shù)值積分進(jìn)行比較發(fā)現(xiàn):以歐拉角、角速度組成狀態(tài)變量,然后直接使用辛歐拉格式不

    能保辛,能量衰減很快,數(shù)值積分存在缺陷;與之相

    反,采用四元數(shù)表示,根據(jù)分析結(jié)構(gòu)力學(xué)的保辛積分方法,并結(jié)合祖沖之類方法的思想,可以很好地避免約束違約,仿真效果令人滿意,可作為陀螺等仿真分析的有力工具.

    本文僅對以歐拉角、角速度組成狀態(tài)變量的數(shù)值積分進(jìn)行研究,對其他形式并未涉及,對其積分效果不佳的成因亦未研究,還有很多工作有待展開.

    參考文獻(xiàn):

    [1]張樹俠, 孫靜. 捷聯(lián)式慣性導(dǎo)航系統(tǒng)[M]. 北京: 國防工業(yè)出版社, 1992: 4880.

    [2]趙曉穎, 溫立書, 么彩蓮. 歐拉角參數(shù)表示下姿態(tài)的2階運動奇異性[J]. 科學(xué)技術(shù)與工程, 2012, 12(3): 634637.

    [3]黃雪樵. 克服歐拉方程奇異性的雙歐法[J]. 飛行力學(xué), 1994, 12(4): 2837.

    [4]李躍軍, 閻超. 飛行器姿態(tài)角解算的全角度雙歐法[J]. 北京航空航天大學(xué)學(xué)報, 2007, 33(5): 505508.

    [5]NIKRAVESH P E, WEHAGE R A, KWON K. Euler parameters in computational kinematics and dynamics: Part 1[J]. J Mechanisms, Transmissions & Automation Des, 1985, 107(3): 358365.

    [6]BETSCH P, SIEBERT R. Rigid body dynamics in terms of quaternions: Hamiltonian formulation and conserving numerical integration[J]. Int J Numer Methods Eng, 2009, 79(4): 444473.

    [7]UDWADIA F E, SCHUTTE A D. An alternative derivation of the quaternion equations of motion for rigidbody rotational dynamics[J]. J Appl Mech, 2010, 77(4): 44505.

    [8]徐小明, 鐘萬勰. 剛體動力學(xué)的四元數(shù)表示及保辛積分[J]. 應(yīng)用數(shù)學(xué)和力學(xué), 2014, 35(1): 111.

    [9]鐘萬勰, 高強(qiáng). 約束動力系統(tǒng)的分析結(jié)構(gòu)力學(xué)積分[J]. 動力學(xué)與控制學(xué)報, 2006, 4(3): 193200.

    [10]鐘萬勰, 高強(qiáng), 彭海軍. 經(jīng)典力學(xué)辛講[M]. 大連: 大連理工大學(xué)出版社, 2013: 202241.

    [11]周江華, 苗育紅, 李宏, 等. 四元數(shù)在剛體姿態(tài)仿真中的應(yīng)用研究[J].飛行力學(xué), 2000, 18(4): 2832.

    [12]HAIRER E, LUBICH C, WANNER G. Geometric numerical integration: structurepreserving algorithms for ordinary differential equations[M]. Berlin: Springer, 2006: 189.

    (編輯于杰)

    国产精品亚洲av一区麻豆| av欧美777| 嫩草影视91久久| 免费在线观看影片大全网站| 亚洲全国av大片| 成年av动漫网址| 天天影视国产精品| 欧美黄色片欧美黄色片| 99精国产麻豆久久婷婷| 香蕉国产在线看| 伦理电影免费视频| 狂野欧美激情性bbbbbb| 1024香蕉在线观看| 汤姆久久久久久久影院中文字幕| 人人澡人人妻人| 欧美成人午夜精品| 国产日韩欧美在线精品| 亚洲av欧美aⅴ国产| 各种免费的搞黄视频| 久久精品国产综合久久久| 蜜桃国产av成人99| 欧美黑人欧美精品刺激| 国产在线视频一区二区| 老司机亚洲免费影院| 色婷婷久久久亚洲欧美| 美国免费a级毛片| 下体分泌物呈黄色| 亚洲国产精品成人久久小说| 狂野欧美激情性bbbbbb| 婷婷成人精品国产| tube8黄色片| kizo精华| 欧美日韩中文字幕国产精品一区二区三区 | 色播在线永久视频| 伊人久久大香线蕉亚洲五| 亚洲精品一二三| 精品国产一区二区三区久久久樱花| 日韩欧美免费精品| 日韩人妻精品一区2区三区| 成人手机av| 精品一区二区三区四区五区乱码| 国产男人的电影天堂91| 国产精品 欧美亚洲| 亚洲国产精品一区三区| 午夜福利视频精品| 五月天丁香电影| 午夜精品国产一区二区电影| 国产亚洲av片在线观看秒播厂| 一区福利在线观看| 最近最新中文字幕大全免费视频| 国产精品亚洲av一区麻豆| 久久青草综合色| 欧美性长视频在线观看| 国产成人系列免费观看| 中国美女看黄片| 熟女少妇亚洲综合色aaa.| 中文字幕精品免费在线观看视频| 99九九在线精品视频| 国产精品一区二区精品视频观看| 婷婷成人精品国产| 亚洲 国产 在线| 亚洲成人免费av在线播放| 男女无遮挡免费网站观看| 欧美在线黄色| 国产免费av片在线观看野外av| 亚洲视频免费观看视频| 日韩,欧美,国产一区二区三区| 久久人妻熟女aⅴ| 久久性视频一级片| 桃花免费在线播放| 亚洲欧美日韩高清在线视频 | 亚洲精品在线美女| av天堂在线播放| 日韩欧美一区视频在线观看| 黄色视频不卡| 成人国产av品久久久| 国产精品免费视频内射| 亚洲成人国产一区在线观看| 国产极品粉嫩免费观看在线| 亚洲成人手机| 日本91视频免费播放| 国产亚洲欧美在线一区二区| 美女大奶头黄色视频| 老司机福利观看| 国产一区二区 视频在线| 亚洲欧美精品综合一区二区三区| 亚洲精品一二三| 欧美精品一区二区免费开放| 国产片内射在线| 欧美乱码精品一区二区三区| 男女无遮挡免费网站观看| 亚洲精华国产精华精| 国产精品一区二区在线不卡| 每晚都被弄得嗷嗷叫到高潮| 久久精品人人爽人人爽视色| 久久久国产一区二区| 一边摸一边抽搐一进一出视频| 午夜福利乱码中文字幕| 欧美亚洲 丝袜 人妻 在线| 黑人欧美特级aaaaaa片| 精品亚洲乱码少妇综合久久| 母亲3免费完整高清在线观看| 后天国语完整版免费观看| 少妇的丰满在线观看| 多毛熟女@视频| 一级毛片精品| 一二三四在线观看免费中文在| 十八禁高潮呻吟视频| 午夜福利,免费看| 久久香蕉激情| 欧美精品啪啪一区二区三区 | av片东京热男人的天堂| 黄色视频在线播放观看不卡| 一本—道久久a久久精品蜜桃钙片| 精品久久久精品久久久| 韩国精品一区二区三区| 久久久水蜜桃国产精品网| 交换朋友夫妻互换小说| 99精品久久久久人妻精品| 国产精品亚洲av一区麻豆| 国产一区二区 视频在线| 欧美性长视频在线观看| 成人国产一区最新在线观看| 他把我摸到了高潮在线观看 | 老司机午夜十八禁免费视频| 久久久久网色| 精品乱码久久久久久99久播| 99热网站在线观看| 亚洲成人免费电影在线观看| 俄罗斯特黄特色一大片| 欧美大码av| 国产精品一二三区在线看| 欧美激情极品国产一区二区三区| 人妻久久中文字幕网| 黑人巨大精品欧美一区二区mp4| 捣出白浆h1v1| 国产成人一区二区三区免费视频网站| 两个人看的免费小视频| 精品国产国语对白av| 91精品伊人久久大香线蕉| 欧美日韩av久久| 亚洲精品一区蜜桃| 亚洲精品美女久久av网站| 亚洲av国产av综合av卡| 大片电影免费在线观看免费| 午夜福利视频精品| 飞空精品影院首页| 欧美97在线视频| 亚洲人成77777在线视频| 香蕉国产在线看| www.熟女人妻精品国产| 国产无遮挡羞羞视频在线观看| 国产欧美亚洲国产| 国产亚洲午夜精品一区二区久久| 男人舔女人的私密视频| 男女边摸边吃奶| xxxhd国产人妻xxx| 极品少妇高潮喷水抽搐| 大香蕉久久成人网| 免费在线观看视频国产中文字幕亚洲 | 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩精品网址| 90打野战视频偷拍视频| 一本—道久久a久久精品蜜桃钙片| 搡老熟女国产l中国老女人| 性色av一级| 在线观看免费视频网站a站| 亚洲欧美精品自产自拍| 一区二区三区四区激情视频| 久热这里只有精品99| 国产欧美日韩一区二区精品| 女性被躁到高潮视频| a级片在线免费高清观看视频| 美女主播在线视频| 黄色视频在线播放观看不卡| 啦啦啦 在线观看视频| 狂野欧美激情性xxxx| 99热全是精品| 狠狠狠狠99中文字幕| 黑人欧美特级aaaaaa片| 国产亚洲av片在线观看秒播厂| 中国美女看黄片| 国产免费视频播放在线视频| kizo精华| av网站在线播放免费| 蜜桃国产av成人99| 日韩人妻精品一区2区三区| 亚洲精品av麻豆狂野| 久久国产精品男人的天堂亚洲| 在线观看免费日韩欧美大片| av不卡在线播放| 精品乱码久久久久久99久播| 久久久国产一区二区| 国产成人欧美| 人成视频在线观看免费观看| 国产免费一区二区三区四区乱码| 中文字幕人妻丝袜一区二区| 一级片'在线观看视频| 999精品在线视频| 成年美女黄网站色视频大全免费| 亚洲色图综合在线观看| 国产精品久久久人人做人人爽| 18禁黄网站禁片午夜丰满| 亚洲午夜精品一区,二区,三区| 日本黄色日本黄色录像| 免费久久久久久久精品成人欧美视频| 极品少妇高潮喷水抽搐| 五月开心婷婷网| 麻豆国产av国片精品| 青春草亚洲视频在线观看| 亚洲欧美清纯卡通| 欧美精品高潮呻吟av久久| 国产男女超爽视频在线观看| 99久久国产精品久久久| 王馨瑶露胸无遮挡在线观看| 久久精品aⅴ一区二区三区四区| 国产熟女午夜一区二区三区| 久久国产精品影院| 久久精品国产亚洲av高清一级| 欧美一级毛片孕妇| 国产免费一区二区三区四区乱码| 女人被躁到高潮嗷嗷叫费观| 午夜福利,免费看| 汤姆久久久久久久影院中文字幕| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品国产av成人精品| av天堂在线播放| 亚洲精品粉嫩美女一区| 欧美日韩福利视频一区二区| 丰满少妇做爰视频| 少妇精品久久久久久久| 交换朋友夫妻互换小说| 久久影院123| 国产精品久久久久久精品古装| 国产欧美日韩一区二区三 | 欧美日韩视频精品一区| 美女福利国产在线| 国产亚洲av片在线观看秒播厂| 精品福利永久在线观看| 午夜激情久久久久久久| 99久久综合免费| 欧美亚洲日本最大视频资源| 80岁老熟妇乱子伦牲交| 久久人妻福利社区极品人妻图片| 老司机午夜十八禁免费视频| 精品一区二区三区av网在线观看 | 成人国语在线视频| 亚洲欧美日韩高清在线视频 | 99香蕉大伊视频| 一区二区三区四区激情视频| 欧美成人午夜精品| 国产免费av片在线观看野外av| 女性生殖器流出的白浆| 欧美在线一区亚洲| 中文字幕色久视频| 啦啦啦在线免费观看视频4| 欧美亚洲日本最大视频资源| 99国产精品免费福利视频| 后天国语完整版免费观看| 久久ye,这里只有精品| 久久国产精品影院| 美女高潮喷水抽搐中文字幕| 亚洲激情五月婷婷啪啪| 亚洲欧洲精品一区二区精品久久久| 久久久久国产精品人妻一区二区| 少妇被粗大的猛进出69影院| 青草久久国产| 丁香六月欧美| 极品少妇高潮喷水抽搐| 国产成人精品无人区| 亚洲欧洲精品一区二区精品久久久| 另类精品久久| www.av在线官网国产| 欧美日韩福利视频一区二区| 国产成人av激情在线播放| 狂野欧美激情性bbbbbb| 亚洲男人天堂网一区| 中文字幕人妻熟女乱码| 在线观看舔阴道视频| 日韩视频在线欧美| 男人操女人黄网站| 丝袜脚勾引网站| 日韩欧美一区二区三区在线观看 | 国产精品久久久久久精品古装| av超薄肉色丝袜交足视频| 夜夜骑夜夜射夜夜干| 成人免费观看视频高清| 看免费av毛片| 久久天堂一区二区三区四区| 国产成人啪精品午夜网站| 涩涩av久久男人的天堂| 菩萨蛮人人尽说江南好唐韦庄| 国产精品国产三级国产专区5o| 99国产精品一区二区蜜桃av | 一级毛片女人18水好多| 黄片小视频在线播放| 久久久国产一区二区| av在线app专区| 欧美日韩一级在线毛片| 建设人人有责人人尽责人人享有的| 欧美精品人与动牲交sv欧美| 下体分泌物呈黄色| 国产福利在线免费观看视频| 亚洲专区字幕在线| 人人妻人人澡人人爽人人夜夜| 国产日韩一区二区三区精品不卡| 欧美黑人精品巨大| 国产精品亚洲av一区麻豆| www.自偷自拍.com| 18禁黄网站禁片午夜丰满| 婷婷色av中文字幕| 香蕉国产在线看| 亚洲精品成人av观看孕妇| 亚洲精品中文字幕一二三四区 | 一个人免费看片子| av片东京热男人的天堂| 涩涩av久久男人的天堂| 十分钟在线观看高清视频www| 亚洲全国av大片| 免费观看a级毛片全部| 天堂俺去俺来也www色官网| 老鸭窝网址在线观看| 欧美性长视频在线观看| 在线看a的网站| 又黄又粗又硬又大视频| 国产精品av久久久久免费| 美女福利国产在线| 国产一区二区三区在线臀色熟女 | 亚洲七黄色美女视频| 一区二区三区乱码不卡18| 欧美日韩中文字幕国产精品一区二区三区 | 大片电影免费在线观看免费| 国产在线观看jvid| 婷婷丁香在线五月| av欧美777| 老司机福利观看| 在线观看一区二区三区激情| 老司机福利观看| 国产成人系列免费观看| 成人三级做爰电影| 国产精品香港三级国产av潘金莲| 国产成人啪精品午夜网站| 精品国产一区二区三区四区第35| 久久久国产一区二区| 国产日韩欧美在线精品| 高清av免费在线| 亚洲第一av免费看| 午夜激情久久久久久久| 人妻 亚洲 视频| 国产精品免费大片| 他把我摸到了高潮在线观看 | 麻豆乱淫一区二区| 国产精品成人在线| 999久久久精品免费观看国产| 大香蕉久久成人网| 免费在线观看影片大全网站| 国产极品粉嫩免费观看在线| 老熟女久久久| 少妇人妻久久综合中文| 午夜福利视频在线观看免费| 中文精品一卡2卡3卡4更新| 在线观看一区二区三区激情| 欧美另类亚洲清纯唯美| 考比视频在线观看| 两个人看的免费小视频| 欧美日本中文国产一区发布| 窝窝影院91人妻| 岛国毛片在线播放| 99九九在线精品视频| 亚洲性夜色夜夜综合| 欧美日韩国产mv在线观看视频| 黑丝袜美女国产一区| 王馨瑶露胸无遮挡在线观看| 免费观看人在逋| 免费在线观看视频国产中文字幕亚洲 | 久久99热这里只频精品6学生| 一二三四社区在线视频社区8| 国产成人av激情在线播放| www.999成人在线观看| 色视频在线一区二区三区| 美女高潮喷水抽搐中文字幕| 午夜免费鲁丝| 国产精品九九99| 精品高清国产在线一区| 91av网站免费观看| 9热在线视频观看99| 亚洲精华国产精华精| 少妇 在线观看| a在线观看视频网站| 建设人人有责人人尽责人人享有的| 大香蕉久久成人网| 午夜精品国产一区二区电影| 日本撒尿小便嘘嘘汇集6| 国产日韩欧美亚洲二区| 中文精品一卡2卡3卡4更新| 久久精品熟女亚洲av麻豆精品| 国产男女内射视频| 亚洲,欧美精品.| 色视频在线一区二区三区| 欧美日本中文国产一区发布| 亚洲av成人一区二区三| 久久精品国产亚洲av高清一级| 成年人午夜在线观看视频| 亚洲国产欧美一区二区综合| 动漫黄色视频在线观看| 别揉我奶头~嗯~啊~动态视频 | av网站在线播放免费| h视频一区二区三区| 一级毛片电影观看| 亚洲精品美女久久久久99蜜臀| 精品国内亚洲2022精品成人 | 一区二区三区四区激情视频| 女性生殖器流出的白浆| 亚洲专区字幕在线| av超薄肉色丝袜交足视频| 秋霞在线观看毛片| 亚洲天堂av无毛| 天天躁夜夜躁狠狠躁躁| 亚洲精品乱久久久久久| av视频免费观看在线观看| 9色porny在线观看| 国产免费现黄频在线看| 超碰成人久久| 丝袜人妻中文字幕| 亚洲精品国产av蜜桃| 性少妇av在线| 国产淫语在线视频| 久久影院123| 成年人免费黄色播放视频| 大片免费播放器 马上看| 久久中文字幕一级| 97人妻天天添夜夜摸| 欧美精品一区二区大全| 精品国产超薄肉色丝袜足j| 日韩有码中文字幕| 国产成人系列免费观看| 欧美另类亚洲清纯唯美| 丝袜喷水一区| av福利片在线| 少妇粗大呻吟视频| 操美女的视频在线观看| 日韩精品免费视频一区二区三区| 9色porny在线观看| 国产一区二区三区在线臀色熟女 | 老鸭窝网址在线观看| 久久久久久亚洲精品国产蜜桃av| 中文字幕制服av| 久久精品亚洲av国产电影网| av免费在线观看网站| 精品免费久久久久久久清纯 | 99国产精品99久久久久| 免费观看av网站的网址| 不卡av一区二区三区| 1024视频免费在线观看| 国产高清视频在线播放一区 | 国产不卡av网站在线观看| 久久青草综合色| 亚洲性夜色夜夜综合| 美女国产高潮福利片在线看| 日本av免费视频播放| 男人舔女人的私密视频| 老司机福利观看| 久久亚洲国产成人精品v| 久久99热这里只频精品6学生| 久久国产精品大桥未久av| 色老头精品视频在线观看| 国产成人影院久久av| 天堂俺去俺来也www色官网| 超色免费av| 侵犯人妻中文字幕一二三四区| 免费女性裸体啪啪无遮挡网站| 国产亚洲精品一区二区www | 久久久久国产精品人妻一区二区| 亚洲av欧美aⅴ国产| 欧美性长视频在线观看| 亚洲性夜色夜夜综合| netflix在线观看网站| 成人亚洲精品一区在线观看| 久久国产亚洲av麻豆专区| 女人久久www免费人成看片| 国产1区2区3区精品| av有码第一页| 91精品伊人久久大香线蕉| 岛国毛片在线播放| 免费一级毛片在线播放高清视频 | 日本五十路高清| av福利片在线| 亚洲熟女毛片儿| 男女免费视频国产| 亚洲欧美一区二区三区久久| 黄色视频不卡| 国产一区二区 视频在线| 在线永久观看黄色视频| 欧美国产精品一级二级三级| 一区二区日韩欧美中文字幕| 国产成人免费观看mmmm| 久久人人97超碰香蕉20202| 中文欧美无线码| 国产视频一区二区在线看| 国产成人精品无人区| 欧美 亚洲 国产 日韩一| 99re6热这里在线精品视频| 国产精品一区二区在线不卡| 啦啦啦啦在线视频资源| 亚洲av国产av综合av卡| 老司机影院成人| 免费观看人在逋| 搡老熟女国产l中国老女人| 国产三级黄色录像| av国产精品久久久久影院| 亚洲精品国产精品久久久不卡| 亚洲中文日韩欧美视频| 亚洲精品粉嫩美女一区| 两个人看的免费小视频| 国产伦理片在线播放av一区| www.精华液| 午夜福利一区二区在线看| 免费久久久久久久精品成人欧美视频| 无遮挡黄片免费观看| 成人国语在线视频| 午夜福利视频在线观看免费| 欧美黄色淫秽网站| 欧美精品高潮呻吟av久久| 日韩中文字幕视频在线看片| 亚洲国产毛片av蜜桃av| 亚洲第一欧美日韩一区二区三区 | 午夜免费成人在线视频| av天堂久久9| 国产三级黄色录像| 九色亚洲精品在线播放| 午夜老司机福利片| 国产又色又爽无遮挡免| 久久精品国产综合久久久| 国产黄色免费在线视频| 精品亚洲成a人片在线观看| 久久影院123| 国产有黄有色有爽视频| 久久久久久免费高清国产稀缺| 亚洲精品乱久久久久久| 老汉色av国产亚洲站长工具| 色老头精品视频在线观看| 99香蕉大伊视频| 在线观看免费视频网站a站| 真人做人爱边吃奶动态| 日韩人妻精品一区2区三区| 国产深夜福利视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产成人精品在线电影| 国产精品国产三级国产专区5o| 久久久久久亚洲精品国产蜜桃av| 亚洲男人天堂网一区| 99国产精品一区二区三区| 久9热在线精品视频| av欧美777| 久久99热这里只频精品6学生| 人人妻人人澡人人爽人人夜夜| 国产99久久九九免费精品| www日本在线高清视频| 久久久欧美国产精品| 成人三级做爰电影| 在线观看免费视频网站a站| 黄色怎么调成土黄色| 国产av精品麻豆| 精品卡一卡二卡四卡免费| 国产精品九九99| 一本综合久久免费| 永久免费av网站大全| 久久精品亚洲熟妇少妇任你| 欧美在线一区亚洲| 精品福利观看| 国产成人影院久久av| 国产一级毛片在线| 欧美日韩一级在线毛片| 一区福利在线观看| 精品卡一卡二卡四卡免费| 免费人妻精品一区二区三区视频| 搡老岳熟女国产| 日日摸夜夜添夜夜添小说| 久久天躁狠狠躁夜夜2o2o| 精品亚洲成a人片在线观看| 精品一区二区三区av网在线观看 | 午夜福利在线免费观看网站| 国产在线免费精品| 日韩 亚洲 欧美在线| 一区二区av电影网| 亚洲第一青青草原| 岛国在线观看网站| 99久久国产精品久久久| 亚洲,欧美精品.| 女性被躁到高潮视频| 国产成人一区二区三区免费视频网站| 久久中文看片网| 欧美激情久久久久久爽电影 | 黄色毛片三级朝国网站| 国产黄色免费在线视频| 国产一区有黄有色的免费视频| 啦啦啦视频在线资源免费观看| 久久久久国内视频| 天天躁狠狠躁夜夜躁狠狠躁| 啦啦啦在线免费观看视频4| 国产黄色免费在线视频| 免费少妇av软件| 99国产精品99久久久久| 国产精品偷伦视频观看了| 精品视频人人做人人爽| 咕卡用的链子| 99香蕉大伊视频| 国产91精品成人一区二区三区 | 亚洲精品久久成人aⅴ小说| netflix在线观看网站| 成人黄色视频免费在线看| 国产av一区二区精品久久| 国产精品久久久久成人av| 人人妻人人澡人人看| 美女高潮喷水抽搐中文字幕|