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

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

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

本文編號(hào)：3877730