隨著國民經(jīng)濟和國防建設(shè)的飛速發(fā)展,對工程應(yīng)用中機械系統(tǒng)產(chǎn)品的動態(tài)性能要求在不斷提高,因此需要對一些較為復(fù)雜的機械系統(tǒng)動力學(xué)特性進(jìn)行準(zhǔn)確且高效地分析和預(yù)測,設(shè)計更加高效穩(wěn)定的數(shù)值算法來滿足系統(tǒng)動力學(xué)的數(shù)值仿真需求,針對多體系統(tǒng)動力學(xué)微分-代數(shù)方程形式,在時間區(qū)間上構(gòu)造A-、L-穩(wěn)定方法,分別基于等距節(jié)點和切比雪夫節(jié)點、勒讓德節(jié)點等非等距節(jié)點建立求解格式,依據(jù)Ehle定理及猜想,與帕德逼近式進(jìn)行對比得到待定系數(shù)矩陣和向量,從而獲得A-、L-穩(wěn)定求解公式,進(jìn)一步將該格式推廣得到任意節(jié)點的待定矩陣和向量,循環(huán)求解過程采用牛頓迭代法計算。以平面雙連桿機械臂系統(tǒng)為例,得到各指標(biāo)的微分-代數(shù)方程組,采用基于均勻節(jié)點的L-穩(wěn)定方法進(jìn)行數(shù)值仿真,對不同時間區(qū)間節(jié)點數(shù)和步長以及各個指標(biāo)的結(jié)果進(jìn)行比較,并與經(jīng)典Runge-Kutta法以及基于等距節(jié)點的A-穩(wěn)定方法、基于非等距節(jié)點的L-穩(wěn)定方法、B-穩(wěn)定的Gauss方法、L-穩(wěn)定的Radau IA、Radau IIA、Lobatto IIIC方法進(jìn)行對比,同時將構(gòu)造的基于等距節(jié)點的A-穩(wěn)定方法與A-穩(wěn)定的Lobatto IIIA、Lobatto IIIB方法...

【文章頁數(shù)】：54 頁

【學(xué)位級別】：碩士

本文編號：3877730