【正文】 e21TVK??（）式中 ——單元剛度矩陣，且2eK ????0000[ ]LLLLTTTTzypBDACENdxIdxINdxGINdx???????（） 單元動力學方程將計入幾何非線性因素的單元動能()、變形能()帶入 Lagrange 方程 ?dTVFt???????????（）得到單元動力學方程為 ?*eeeMKPQ（） ?*112TeeK???（） ??18Tdiag????（） 121178eeeK??????????????（）將未計入幾何非線性因素的單元動能()、變形能()帶入 Lagrange 方程得到單元動力學方程為山東科技大學工程碩士學位論文 高速空間并聯(lián)坐標測量機的彈性動力學建模20 ?2eeeMKFPQ????（） 式中——單元外加載荷的廣義力列陣；eF——是其它單元給予所分析單元的作用力列陣；P，是系統(tǒng)單元的剛體慣性力列陣，對系統(tǒng)作剛體運動分析后即可得出。e=rQM?? 下面求解系統(tǒng)坐標系下的單元彈性動力學方程。首先將單元局部坐標系 中的oxyz?廣義坐標列陣 轉化為系統(tǒng)坐標系下（這里指定系 中） 。引入定系 下的單元廣義δ??A??A坐標列陣 ij? （,Tijijij???????）引入以下坐標旋轉變換矩陣 0000ijijijij ijijijRR?????????（） 其中 12122122iiiii i iij iiiii i icscscsR?????? ?? ?? ???????由此得到廣義坐標列陣 和 之間的關系表達式δij ijR??（） 將式()對時間 t 求一、二階導數(shù)，得到 ijij?????（）山東科技大學工程碩士學位論文 高速空間并聯(lián)坐標測量機的彈性動力學建模21 2ijijijR??????（）將（） （）兩式代入（）得到在定坐標系下考慮幾何非線性因素時的梁單元的彈性運動微分方程： ij ij ijijij ij ijMCKP??????（）式中 TijijeijMR? 2ijijijC? *1TTijijeijijijeKK? ()ijijPRFQ?定坐標系下未計入非線性因素時梁單元彈性動力學方程和上面的相同只是 不同，ijK且 2TTijijeijijijeRMR???（）由以上得到了考慮幾何非線性因素和未考慮幾何非線性因素影響的定系下的單元動力學方程。在下文的算例中將給出計入幾何非線性因素和不計入幾何非線性因素之間的比對圖。 支鏈彈性動力學方程的建立4UPS/UPU 高速空間并聯(lián)坐標測量機是由定平臺、動平臺以及連接它們的分支等組成的。其中，定平臺通過一個驅動分支 UPU（虎克鉸 —移動副—虎克鉸）和四個結構完全相同的驅動分支 UPS（虎克鉸—移動副—球鉸） ，與動平臺相連接；要實現(xiàn)并聯(lián)坐標測量機動平臺的位姿變化，可通過改變其驅動桿的桿長值來實現(xiàn)。假定擺動桿 為剛性桿，伸縮桿 為柔性桿，動平臺 為iPUiSP)5,21(??54321S剛性架，關節(jié)柔性忽略不計（由于并聯(lián)機構的支鏈一般較短，而關節(jié)變形引起的系統(tǒng)累積誤差也比較小，且并聯(lián)機構多為閉環(huán)系統(tǒng)，對系統(tǒng)的關節(jié)變形具有一定的約束作用，故關節(jié)變形可以忽略） 。這里將伸縮桿 等分為 個梁單元，并將這些單元依次編號iSm山東科技大學工程碩士學位論文 高速空間并聯(lián)坐標測量機的彈性動力學建模22為 ，其中 表示擺動桿和伸縮桿在第 個單元接觸。1,2m? nn 由以上假定建立支鏈的彈性動力學方程，只需對伸縮桿進行分析即可。 當 j= 時，由于這些單元被剛性體包圍其廣義坐標都為零，即,1?? 0ij?? 當 j= 時，由于第 個單元的左端被剛體包圍，此單元的前九個廣義坐n?n標為零。即 。此單元的右端始終重合于下一個單元的左端，從此1290iini????單元以后的單元也是如此。則有： （）10(1)1(1)212(1)313(1)414(1)515(1)616(1)717(1)818(1)9ij ijij ijij ijij ijij ijij ijij ijij ijij ij??????????????????????????對于 單元考慮到動平臺鉸鏈對單元 右端點的約束，不同支鏈 端點處的曲率不同。mmi當 i=1 時 （）1780iim??當 i=2,3,4,5 時 （）6178iii??也就是說當 i=1 時， 是獨立廣義坐標，但在 i=1,2,3,4,5 時為 0。i由以上分析可以從各個單元中得到支鏈的廣義坐標。當 時，支鏈廣義坐標個數(shù)為 個。=1i ??9+7n具體表示為： （）(1)()2(1)9(2)1916,,Tiniiniimiiq?????? ?????? ?時，支鏈廣義坐標個數(shù)為有 個。=345??+n具體表示為：山東科技大學工程碩士學位論文 高速空間并聯(lián)坐標測量機的彈性動力學建模23 （）(1)()2(1)9(2)1915,,Tiniiniimiiq?????? ?????? ?為得出支鏈的彈性動力學方程，需要先建立 與 的表達式間的關系，iij?根據(jù)測量機支鏈的約束方程組，將列陣 與 關系表達式表示為iij （）ijijAq??不同的支鏈轉換關系矩陣也不同。當 1i? （）??8(97)1(92)89()18921()60(1,3,))0(1,)0 )mnijjnmjmjnBnAEjmG??????????????????? ???? ??? ? ?當 時=2,345i （）??18(97)(92)189()18925(1,3,)0)(1,))mnijjnmjmjnBnAEjmG?????????? ????????????????? ? ?式中矩陣 ，將 時的矩陣 去掉最后一列，即為 時的矩陣 ，i?186? 2,345i 185G?矩陣 、 在各個支鏈上都相同。 表示 1818 階的單位矩陣， 、 矩陣表示如EBEB186?下：山東科技大學工程碩士學位論文 高速空間并聯(lián)坐標測量機的彈性動力學建模24 000000000101001000B? ?? ?? ???????? ?189??????（）山東科技大學工程碩士學位論文 高速空間并聯(lián)坐標測量機的彈性動力學建模25 （）10000010000010000G?1860100000?? ?? ?? ?? ?? ?? ?? ?? ?將 代入方程 ()兩邊并左乘 得：ijijAq??ijA A???ijijijijMq+CKq=P??（）式中 ?TijijMA?（） ?TijijC（） ?TijijKA?（） TijijP（）上式是對一個單元列出的方程，機構支鏈 中有多少個單元，就可以列出多少個形如i（）的方程，將這些方程疊加，得到山東科技大學工程碩士學位論文 高速空間并聯(lián)坐標測量機的彈性動力學建模26 （） ?MqCKqPiiii????式中 A1mijijM??（ ） （）A1mijijC?? （）1ijijK? （）A1miijjP?由此得到該測量機支鏈 的彈性動力學方程。 系統(tǒng)彈性動力學方程的建立 系統(tǒng)運動學約束動平臺與支鏈連接時要滿足以下條件：①各個支鏈與動平臺連接點的位移必須與動平臺上與各個支鏈連接點的位移相一致；②各個支鏈對動平臺的作用力之和應與作用于動平臺的外力和慣性力相平衡。下面根據(jù)條件① ，進行系統(tǒng)運動學約束的推導。由圖 可知，在動系 下，各個鉸鏈與動平臺的連接點 的坐標是??B )5,4321(?iS定值。這里動平臺的六個自由度定義為動系 沿著系統(tǒng)定系 的BBZYXO?AAZYXO?變化。顯然，鉸鏈點 處的坐標值分別是上述六個獨立變量的函數(shù)。因此，)5,4321(?iS動平臺的位移量可以表示為 。??Tqq6543210,由第二章可知，動系 相對于定系 的變換矩陣為 ，當給出動平臺的??BARAB型歐拉角 和動系 原點 在定系 下的坐標XYZ?),(???BO??時，由式()，以歐拉角的形式表示動系 相對于定系 的變??TBBOAZ?P ??A換矩陣： ???????????? 100),( BZYXAB Zcscs YsXc?????????T山東科技大學工程碩士學位論文 高速空間并聯(lián)坐標測量機的彈性動力學建模27（） 假設動平臺的名義運動位姿在點 處，由于機構支鏈的彈性變形使動平臺的實際BO運動位姿發(fā)生了微小的變動（即 ） ，最終移動到了點 處。動平BZYX?????, BO?臺與支鏈的約束關系如圖 所示。則由坐標系 到坐標系 的??? ZYX?變換矩陣為 T? （）?????????? 10)cos()sin()sin(i ico00)si(co)sin( )cos(n)si(i iicoc BZYX????????????? ??? 動平臺與支鏈的約束關系Fig. Constraint between the moving platform and branch因為六個變換參數(shù) 都是微小的變量，由泰勒展開式以及麥克Bzyx?????,勞林公式，可以得到 的近似表達式T?山東科技大學工程碩士學位論文