連載-“老朱陪您學(xué)Civil”之梁橋抗震4

2016-10-29 by:CAE仿真在線來源:互聯(lián)網(wǎng)

連載-“老朱陪您學(xué)Civil”之梁橋抗震,你們的最愛-朱老師,小邁致力于與您分享一手技術(shù)資料,用技術(shù)創(chuàng)造幸福,你我一起!

連載-“老朱陪您學(xué)Civil”之梁橋抗震4Midas Civil分析案例圖片1

3.Civil程序計(jì)算振型的三種方法

在一般的有限元分析中,由于系統(tǒng)的自由度很多,同時(shí)在研究系統(tǒng)的響應(yīng)時(shí),往往只需要了解少數(shù)較低的特征值及相應(yīng)的特征向量,因此在有限元分析中,發(fā)展了一些適應(yīng)上述特點(diǎn)而效率較高的解法(子空間迭代法、lanczos)。

midas Civil中除了提供精確的特征向量法分析外,還提供了與荷載相關(guān)的Ritz向量分析法。多重Ritz向量能用于線性和非線性結(jié)構(gòu)的動(dòng)力分析。與精確特征向量法相比,多重Ritz向量法用更少的時(shí)間可產(chǎn)生更精確的結(jié)果。

3.1子空間迭代法(wilson著作《結(jié)構(gòu)靜力與動(dòng)力分析》)

連載-“老朱陪您學(xué)Civil”之梁橋抗震4Midas Civil分析案例圖片2

子空間迭代法是求解大型矩陣特征值問題的最常用最有效的方法之一,它適合于求解部分特征值解,被廣泛應(yīng)用于結(jié)構(gòu)動(dòng)力學(xué)的有限元分析中。

子空間迭代法是假設(shè)r個(gè)起始向量(采用移頻法,通過特征值的移動(dòng)和已收斂的特征向量的移出,使r保持在較小的數(shù)值,從而顯著提高計(jì)算效率和改進(jìn)收斂速度)同時(shí)進(jìn)行迭代(通過求解減縮廣義特征值問題)以求得矩陣的前p(<r)個(gè)特征值和特征向量。(如果r不是足夠大,一方面可能漏掉可能激起的振型;另一方面又可能引入不可能激起的振型)。

3.2 Lanczos方法

連載-“老朱陪您學(xué)Civil”之梁橋抗震4Midas Civil應(yīng)用技術(shù)圖片3

Lanczos方法和Ritz向量法的共同特點(diǎn)是直接生成一組Lanczos向量或者Ritz向量,對(duì)運(yùn)動(dòng)方程進(jìn)行縮減,然后求解縮減了的運(yùn)動(dòng)方程的特征值問題,避免了迭代步驟(采用直接疊加法),從而具有更高的計(jì)算效率。

Lanczos法和Ritz向量法本質(zhì)上一致,但是在實(shí)際計(jì)算中,由于計(jì)算機(jī)的截?cái)嗾`差和舍入誤差,導(dǎo)致數(shù)值上的不穩(wěn)定性(例如虛假的多重特征值現(xiàn)象),因此妨礙了Lanczos方法的實(shí)際應(yīng)用。20世紀(jì)70年代以后,很多研究工作者提出了不少Lanczos向量的重正交技術(shù)以調(diào)高其算法的穩(wěn)定性,Ritz向量法從這個(gè)意義上說可以是這種,但由于他改變了成Lanczos向量的算法公式,導(dǎo)致以后求解的不是對(duì)角矩陣的特征值問題,而是一般矩陣特征值問題。

3.3多重Ritz向量法(求解的是一般矩陣特征值問題)

連載-“老朱陪您學(xué)Civil”之梁橋抗震4Midas Civil培訓(xùn)教程圖片4

多重Ritz向量法認(rèn)可結(jié)構(gòu)動(dòng)態(tài)響應(yīng)是空間荷載分布的函數(shù),考慮動(dòng)力荷載的空間分布(當(dāng)定義了初始向量后,第一個(gè)向量塊的靜態(tài)響應(yīng)就來源于該初始荷載向量),可以避免漏掉可能激起的振型和引入不可能激起的振型,能夠顯著提高計(jì)算效率。

3.4 子空間迭代法、Lanczos法和多重Ritz向量法算例比較(荷載作用是對(duì)稱)

連載-“老朱陪您學(xué)Civil”之梁橋抗震4Midas Civil培訓(xùn)教程圖片5

由圖可以清楚的發(fā)現(xiàn):多重Ritz向量法和另外2種分析方法得到的振型圖是大不相同的。通過模態(tài)查看,可以發(fā)現(xiàn)多重Ritz向量法所有振型都是對(duì)稱的(荷載作用是對(duì)稱的),因?yàn)樗紤]了空間荷載分布狀態(tài)及動(dòng)力貢獻(xiàn),所以他忽略了所有反對(duì)稱振型。對(duì)于反對(duì)稱振型,并不是由荷載激發(fā)的,荷載在這些振型的動(dòng)力貢獻(xiàn)為零。故:反應(yīng)譜分析時(shí)較好的方法是采用多重Ritz向量法。

4.已知振型(模態(tài))情況下如何計(jì)算反應(yīng)譜下的結(jié)構(gòu)受力情況

4.1單振型施加慣性力的方法