【正文】 式（284）的條件都成立。 ()則存在滿足下式（285）的拉格朗日乘子矢量。 ()其中為任意的。在式()和()中，，，運用拉格朗日乘子定理于式()和()，則存在拉格朗日乘子矢量，對于任意的應(yīng)滿足： ()由此得到運動方程的拉格朗日乘子形式： ()式()還必須滿足式()、()和()表示的位置約束方程、速度約束方程及加速度約束方程，如下： ()， ()， ()以上三式其維數(shù)同式()。式()、()、()和()組成約束機械系統(tǒng)的完整的運動方程。將式()與()聯(lián)立表示為矩陣形式： ()式()即為多體系統(tǒng)動力學(xué)中最重要的動力學(xué)運動方程，被稱為歐拉拉格朗日方程（EulerLagrange Equation），式()還必須滿足式()和()。它是一個微分代數(shù)方程組（Differential Algebraic Equations DAEs），不同于單純的常微分方程組（Ordinary Differential Equations ODEs）問題，其求解關(guān)鍵在于避免積分過程中的違約現(xiàn)象，此外，還要注意DAE問題的剛性問題。顯然，式()有且僅有唯一解的充要條件是其系數(shù)矩陣非奇異，但這一條件不利于實際中的判斷，可以給出更為實用的判斷。如果式()滿足如下條件：a. ，；b. 對任意且。則式()中的系數(shù)矩陣是非奇異的，且和是唯一確定的。這是多體系統(tǒng)運動方程解的存在定理。可以據(jù)此判斷，如果系統(tǒng)質(zhì)量矩陣是正定的，并且約束獨立，那么運動方程就有唯一解。實際中的系統(tǒng)質(zhì)量矩陣通常是正定的，只要保證約束是獨立的，運動方程就會有解。在實際數(shù)值迭代求解過程中，需要給定初始條件，包括位置初始條件和速度初始條件。此時，如果要使運動方程有解，還需要滿足初值相容條件，也就是要使位置初始條件滿足位置約束方程，速度初始條件滿足速度約束方程。對于由式()及()、()確定的系統(tǒng)動力學(xué)方程，初值相容條件為： () ()（3）正向動力學(xué)分析、逆向動力學(xué)分析與靜平衡分析對于一個確定的約束多體系統(tǒng)，其動力學(xué)分析不同于運動學(xué)分析，并不需要系統(tǒng)約束方程的維數(shù)等于系統(tǒng)廣義坐標(biāo)的維數(shù)。在給定外力的作用下，從初始的位置和速度，求解滿足位置約束式()及速度約束式()的運動方程式()，就可得到系統(tǒng)的加速度和相應(yīng)的速度、位置響應(yīng)，以及代表約束反力的拉格朗日乘子，這種已知外力求運動及約束反力的動力學(xué)分析，稱為正向動力學(xué)分析。如果約束多體系統(tǒng)約束方程的維數(shù)與系統(tǒng)廣義坐標(biāo)的維數(shù)相等，也就是對系統(tǒng)施加與系統(tǒng)自由度相等的驅(qū)動約束，那么該系統(tǒng)在運動學(xué)上就被完全確定，、速度方程和加速度方程可求解系統(tǒng)運動。在此情況下，式()的雅可比矩陣是非奇異方陣，即： ()展開式()的運動方程，為： () ()由式()可解得，再由式()可求得，拉格朗日乘子就唯一地確定了作用在系統(tǒng)上的約束力和力矩（主要存在于運動副中）。這種由確定的運動求系統(tǒng)約束反力的動力學(xué)分析就是逆向動力學(xué)分析。如果一個系統(tǒng)在外力作用下保持靜止?fàn)顟B(tài)，也就是說，如果： ()那么，就說該系統(tǒng)處于平衡狀態(tài)。將式()代入運動方程式()，得到平衡方程： ()由平衡方程式()及約束方程式()可求出狀態(tài)和拉格朗日乘子。這種求系統(tǒng)的平衡狀態(tài)及在平衡狀態(tài)下的約束反力的動力學(xué)分析稱為（靜）平衡分析。（4）約束反力對于約束機械系統(tǒng)中的構(gòu)件，設(shè)其與系統(tǒng)中某構(gòu)件存在運動學(xué)約束或驅(qū)動約束，約束編號為。除連體坐標(biāo)系外，再在構(gòu)件上以某點為原點建立一個新的固定于構(gòu)件上的坐標(biāo)系，稱為運動副坐標(biāo)系，設(shè)從坐標(biāo)系到坐標(biāo)系的變換矩陣為，從坐標(biāo)系到坐標(biāo)系的變換矩陣為，則可導(dǎo)出由約束產(chǎn)生的反作用力和力矩分別為： () ()以上兩式中，為約束對應(yīng)的拉格朗日乘子，反作用力和力矩均為運動副坐標(biāo)系中的量。2．三維系統(tǒng)的廣義坐標(biāo)比二維系統(tǒng)復(fù)雜得多，使得問題規(guī)模更大。這里討論的是與二維多體系統(tǒng)動力學(xué)分析相應(yīng)的內(nèi)容，包括微分代數(shù)混合方程組的建立，三種類型的動力學(xué)分析等。只是根據(jù)三維系統(tǒng)情況的不同，給出了角速度表示和歐拉參數(shù)表示的兩種不同形式的運動方程。（1）空間自由剛體的變分運動方程對于空間任意剛體構(gòu)件，令其連體坐標(biāo)系原點固定于剛體質(zhì)心，此時連體坐標(biāo)系也稱為質(zhì)心坐標(biāo)系，設(shè)剛體質(zhì)量為，其相對于質(zhì)心坐標(biāo)系的慣性張量為，再設(shè)作用在剛體上的總外力，外力相對于質(zhì)心坐標(biāo)系原點的力矩為，則相對于剛體質(zhì)心坐標(biāo)系的剛體牛頓歐拉變分運動方程為： ()其中，為剛體質(zhì)心的虛位移，為剛體的虛轉(zhuǎn)動，為剛體質(zhì)心位移，為剛體的在坐標(biāo)系中表示的角速度。（2）空間約束機械系統(tǒng)的運動方程—角加速度形式考慮由個剛體組成的空間約束機械系統(tǒng)，系統(tǒng)的廣義坐標(biāo)選為： () ()定義： () () () () () () ()110） ()運用上述符號將系統(tǒng)中每個構(gòu)件的牛頓歐拉方程式()總和為： ()同樣地，在理想約束情況下，對于一個系統(tǒng)只要考慮外力或外力矩，如此，由式()可得到約束機械系統(tǒng)的變分運動方程： ()此式必須適用于所有運動學(xué)上所容許的虛位移和虛轉(zhuǎn)動。系統(tǒng)中的運動學(xué)約束和驅(qū)動約束如式()及()，聯(lián)立表示為： ()與式()和式()不同的是，動力學(xué)約束中允許總體上大于零的自由度出現(xiàn)，設(shè)其維數(shù)為。此外，采用歐拉參數(shù)還須滿足歐拉參數(shù)歸一化約束（其維數(shù)為nb）： ()對式()微分，得到用虛位移和虛轉(zhuǎn)動表示的變分形式： ()式()即為系統(tǒng)虛位移和虛轉(zhuǎn)動在運動學(xué)上所容許的條件。式中不包括式()所表示的歐拉參數(shù)歸一化約束，因為如式()所表明的，當(dāng)以為變量時，與歐拉參數(shù)歸一化約束相關(guān)的速度方程是自動得到滿足的，同樣，以為變量時，與歐拉參數(shù)歸一化約束相關(guān)的變分方程也是自動得滿足的。對式()和()應(yīng)用拉格朗日乘子定理，則存在一個拉格朗日乘子矢量，滿足： ()式中和是任意的，由此導(dǎo)出空間約束機械系統(tǒng)的牛頓歐拉運動方程為： () ()對式()分別求一次及兩次導(dǎo)數(shù)，： () .()式中速度右項和加速度右項與式()和()中定義相同。式()和()維數(shù)與式()相同，為。由式()、()及()可以得到方陣形式的系統(tǒng)運動方程為： ()式()與約束方程式()及速度方程式()一起組成描述nb個剛體系統(tǒng)運動的微分代數(shù)方程組。式()的維數(shù)關(guān)系滿足：。由于角速度是不可積的，因此應(yīng)把式()看作是速度變量及代數(shù)變量的一階微分代數(shù)混合方程。與二維約束機械系統(tǒng)類似，如果滿足：a. 矩陣行滿秩；b. 對任意且，存在。則式()中系數(shù)矩陣非奇異，系統(tǒng)運動方程有唯一解。在實際迭代計算中，需要給定初始位置條件和，以及初始速度條件和，初始條件須分別滿足位置約束方程和速度約束方程： () ()（3）空間約束機械系統(tǒng)的運動方程—歐拉參數(shù)形式，可以得出用歐拉參數(shù)導(dǎo)數(shù)表示的空間機械系統(tǒng)運動方程。，可進一步導(dǎo)出歐拉參數(shù)二階導(dǎo)數(shù)與角加速度關(guān)系，對于單個剛體的這些運動量間關(guān)系，在這里一并給出： () () () () ()則用歐拉參數(shù)表示式()所示的系統(tǒng)速度方程為： ()130）用歐拉參數(shù)表示式()所示的系統(tǒng)加速度方程為： () 式（）和（）維數(shù)同樣為。對歐拉參數(shù)式()微分，得到歐拉參數(shù)速度方程： ()其維數(shù)關(guān)系為，且其中： ()對式()微分，得到歐拉參數(shù)加速度方程： ()其維數(shù)關(guān)系同式()。約束方程變分式()用歐拉參數(shù)表示為： ()歐拉參數(shù)約束方程()的變分式為： ()對于用角速度表示的約束機械系統(tǒng)的變分運動方程式()，根據(jù)式() 替換，根據(jù)式()替換，根據(jù)式()替換，可得到用歐拉參數(shù)表示的系統(tǒng)變分運動方程： ()其中，： ()式()必須對滿足式()和()的所有和成立。應(yīng)用拉格朗日乘子定理于式()及()和()，則有拉格朗日乘子和使下式對任意和成立。 ()從而得到歐拉參數(shù)形式的約束機械系統(tǒng)變分運動方程： () ()由式()和()，再加上運動學(xué)約束和驅(qū)動約束的加速度方程式()及歐拉參數(shù)約束的加速度方程式()，可得到方陣形式的空間約束機械系統(tǒng)運動方程為： ()式()與式()和式()的運動學(xué)約束和歐拉參數(shù)歸一化約束以及式()和式()的相應(yīng)速度方程一起，構(gòu)成歐拉參數(shù)形式表達的系統(tǒng)微分代數(shù)運動方程的混合方程組。其維數(shù)關(guān)系為：式()的系數(shù)矩陣非奇異的條件同本節(jié)二維多體系統(tǒng)的運動方程解存在條件。實際求解中需要給出初始位置和及初始速度和，且初始位置條件必須滿足位置約束方程，初始速度條件必須滿足速度約束方程，如下：