【JY】主成分分析與振型分解

不等待


即關(guān)注

【JY】主成分分析與振型分解的圖1

簡介

   世界的事物原本統(tǒng)一,是因?yàn)閷I(yè)賦予我們看待世界的角度,今天分別從經(jīng)濟(jì)學(xué)方法(主成分分析)和結(jié)構(gòu)分析方法(振型分解法)來討論下這個(gè)東西。

【JY】主成分分析與振型分解的圖2

主成分分析法&振型分解法

    首先小談下主成分分析法(principal components analysis),也稱主分量分析,是利用降維的思想,在損失很少信息的前提下,把多個(gè)指標(biāo)轉(zhuǎn)化為幾個(gè)綜合指標(biāo)的多元統(tǒng)計(jì)方法。通常把轉(zhuǎn)化生成的綜合指標(biāo)稱為主成分,其中每個(gè)主成分都是原始變量的線性組合,且各個(gè)主成分之間互不相關(guān),使得主成分比原始變量具有某些更優(yōu)越的性能。這樣在研究復(fù)雜問題時(shí)就可以只考慮少數(shù)幾個(gè)主成分而不至于損失太多信息,從而更容易抓住主要矛盾,揭示事物內(nèi)部變量之間的規(guī)律性,同時(shí)使問題得到簡化,提高分析效率。

    再談下 振型分解法,在討論多自由度體系的強(qiáng)迫振動(dòng)時(shí),如采用質(zhì)點(diǎn)位移作為坐標(biāo)(稱為幾何坐標(biāo)),則所得到的振動(dòng)方程為耦聯(lián)微分方程,因而必須聯(lián)立求解。對(duì)于無阻尼簡諧強(qiáng)迫振動(dòng),在平穩(wěn)階段,由于各質(zhì)點(diǎn)都作同步振動(dòng),利用這一特性可將微分方程轉(zhuǎn)化為代數(shù)方程,故求解沒有困難。然而,當(dāng)考慮阻尼影響或者在一般動(dòng)力荷載作用下時(shí),求解聯(lián)立的微分方程組就會(huì)比較困難。按振型分解的計(jì)算方法就是針對(duì)這一問題提出來的。振型分解法是 基于坐標(biāo)變換,把原來耦聯(lián)的微分方程組變?yōu)閚個(gè)互相獨(dú)立的微分方程,從而使原來多自由度體系的動(dòng)力計(jì)算變?yōu)橐幌盗袉巫杂啥润w系的問題,當(dāng)然這一方法只限于線性體系的應(yīng)用。下面介紹振型分解法。
    從上可以看出,兩種方法的本質(zhì)都是通過從結(jié)構(gòu)內(nèi)部(事物之間)錯(cuò)綜復(fù)雜的關(guān)系中找出一些振型(主要成分),從而得到結(jié)構(gòu)振動(dòng)的動(dòng)力特性(或從而能有效利用大量統(tǒng)計(jì)數(shù)據(jù)進(jìn)行定量分析),揭示結(jié)構(gòu)在不同地震動(dòng)下各振型的參與成份(變量之間的內(nèi)在關(guān)系),得到對(duì)結(jié)構(gòu)(事物)特征及其發(fā)展規(guī)律的一些深層次的啟發(fā),把研究工作引向深入。
    接下來分別淺析下主成分分析法和振型分解法的原理。
【JY】主成分分析與振型分解的圖3
【JY】主成分分析與振型分解的圖4

【主成分分析法】

    (詳見《多元統(tǒng)計(jì)分析》)設(shè)對(duì)某一事物的研究涉及p個(gè)指標(biāo),分別用X1,X2,…,X。表示,這p個(gè)指標(biāo)構(gòu)成的p維隨機(jī)向量為X=(X1,X2,…,Xp)'。設(shè)隨機(jī)向量X的均值為μ。

    對(duì)X進(jìn)行線性變換,可以形成新的綜合變量,用Y表示,也就是說,新的綜合變量可以由原來的變量線性表示,即滿足下式:

【JY】主成分分析與振型分解的圖5

    而它的幾何意義在于,在處理涉及多個(gè)指標(biāo)問題的時(shí)候,為了提高分析的效率,可以不直接對(duì)p個(gè)指標(biāo)構(gòu)成的p維隨機(jī)向量X=(X1,X2,…,Xp)'進(jìn)行分析,而是先對(duì)向量X進(jìn)行線性變換,形成少數(shù)幾個(gè)新的綜合變量Y1,Y2,…,Yp,使得各綜合變量之間相互獨(dú)立且能解釋原始變量盡可能多的信息。這樣,在以損失很少部分信息為代價(jià)的前提下,達(dá)到簡化數(shù)據(jù)結(jié)構(gòu),提高分析效率的目的。

【JY】主成分分析與振型分解的圖6

主成分分析方法的幾何意義

    進(jìn)行主成分分析的目的之一是減少變量的個(gè)數(shù),所以一般不會(huì)取p個(gè)主成分,而是取m<p個(gè)主成分。m取多少比較合適,是一個(gè)很實(shí)際的問題,通常以所取m使得累積貢獻(xiàn)率達(dá)到85%以上為宜( 這就等同振型分解法中的質(zhì)量參與系數(shù),其中特征值λ等同振型分解法的圓頻率w),即:
【JY】主成分分析與振型分解的圖7
【JY】主成分分析與振型分解的圖8

【振型分解法】

   首先第一步建立動(dòng)力方程,此處不做太多解釋,可看《往期  結(jié)構(gòu)動(dòng)力學(xué)初步——單質(zhì)點(diǎn)結(jié)構(gòu)的瞬態(tài)動(dòng)力學(xué)分析》,也可通過基本的力學(xué)原理進(jìn)行建立動(dòng)力方程。

    如:通過微分方法——非變分的方法(如牛頓第二定律、達(dá)朗貝爾原理、拉格朗日方程);變分的方法(如虛功原理)。

    或通過積分的方法——非變分的方法(如能量守恒原理、動(dòng)量變化定理);變分的方法(如哈密頓原理)。

    將一個(gè)n維自由度的結(jié)構(gòu),動(dòng)力方程可寫成:
【JY】主成分分析與振型分解的圖9
矩陣形式即: 【JY】主成分分析與振型分解的圖10
給一個(gè)不改變結(jié)構(gòu)特性的振動(dòng),則可以計(jì)算得到:
【JY】主成分分析與振型分解的圖11

矩陣形式即:

【JY】主成分分析與振型分解的圖12

其中
【JY】主成分分析與振型分解的圖13
將上述求解方程帶入運(yùn)動(dòng)方程中可得到一個(gè)齊次方程組:
【JY】主成分分析與振型分解的圖14
矩陣形式為: 【JY】主成分分析與振型分解的圖15
即可得到: 【JY】主成分分析與振型分解的圖16
若要使振型 【JY】主成分分析與振型分解的圖17恒不為零,則要滿足:   【JY】主成分分析與振型分解的圖18
    通過上式即可求得結(jié)構(gòu)的 圓頻率和振型。對(duì)于特別復(fù)雜或則超多自由度的結(jié)構(gòu),通常會(huì)取質(zhì)量參與系數(shù)>90%或95%。
【JY】主成分分析與振型分解的圖19
【JY】主成分分析與振型分解的圖20

【算例】

計(jì)算串模型8層樓,每層質(zhì)量m=10t,每層剛度 10^8 N/m ,計(jì)算前六階振型和周期。

根據(jù)結(jié)構(gòu)振型分解法計(jì)算分析可得:
Sap2000結(jié)果呈現(xiàn):
【JY】主成分分析與振型分解的圖21

Etbas結(jié)果呈現(xiàn):

【JY】主成分分析與振型分解的圖22

自編Matlab結(jié)果:

【JY】主成分分析與振型分解的圖23

根據(jù)主成分分析方法計(jì)算分析可得:

我們將需要處理的“大數(shù)據(jù)”定為[M]×[K],則有下列數(shù)據(jù):
【JY】主成分分析與振型分解的圖24

利用Matlab做主成分分析可得:

(各特征值和累計(jì)貢獻(xiàn)率)

【JY】主成分分析與振型分解的圖25

將特征值繪制成圖可得到:

利用統(tǒng)計(jì)學(xué)軟件SAS做主成分分析可得:
【JY】主成分分析與振型分解的圖26
【JY】主成分分析與振型分解的圖27

(注:求特征值方法有很多種,Matlab的princomp命令和SAS內(nèi)置的計(jì)算方法選擇有點(diǎn)差別,但結(jié)果基本一致)

【JY】主成分分析與振型分解的圖28

 (注:藍(lán)色為Matlab的princomp命令,紅色為SAS內(nèi)置計(jì)算)

它們分別計(jì)算得到的周期為:
【JY】主成分分析與振型分解的圖29
【JY】主成分分析與振型分解的圖30
【JY】主成分分析與振型分解的圖31

總結(jié)

【JY】主成分分析與振型分解的圖32

    同一種事物在不同身份的人看來是不一樣的,數(shù)學(xué)的人稱它為特征向量,土木的人稱它為振型,經(jīng)濟(jì)學(xué)的人稱它為主成分,當(dāng)然還有累計(jì)貢獻(xiàn)率-振型參與系數(shù)…… 

【JY】主成分分析與振型分解的圖33
【JY】主成分分析與振型分解的圖34
【JY】主成分分析與振型分解的圖35
【JY】主成分分析與振型分解的圖36

    所以說世界消解統(tǒng)一的宇宙秩序,解構(gòu)了由這種宇宙秩序所設(shè)定的統(tǒng)一的價(jià)值世界,價(jià)值世界出現(xiàn)了多元分化。世界原本是沒有意義的,人們?yōu)榱藬[脫和超越當(dāng)下世界的短暫性和有限性,通過預(yù)設(shè)與想象,以各種身份(工程師、經(jīng)濟(jì)學(xué)家、醫(yī)生)的角度賦予世界以某種意義, 也就是說世界原本是統(tǒng)一的,但同時(shí)也是多元的,對(duì)于同一個(gè)事物,從不同的角度理解會(huì)得到不一樣的結(jié)論,世界不是非黑即白、非對(duì)即錯(cuò)。

概念為先,機(jī)理為本,下期再見!

建源之光——工程俠

 往期推薦 ·


往期精彩

#性能分析

【JY】基于性能的抗震設(shè)計(jì)淺析(一)

【JY】基于性能的抗震設(shè)計(jì)淺析(二)

【JY】淺析消能附加阻尼比

【JY】近斷層結(jié)構(gòu)設(shè)計(jì)策略分析與討論

#概念機(jī)理

【JY】基于Ramberg-Osgood本構(gòu)模型的雙線性計(jì)算分析

【JY】結(jié)構(gòu)瑞利阻尼與經(jīng)濟(jì)訂貨模型

【JY】結(jié)構(gòu)動(dòng)力學(xué)初步-單質(zhì)點(diǎn)結(jié)構(gòu)的瞬態(tài)動(dòng)力學(xué)分析

【JY】從一根懸臂梁說起

【JY】反應(yīng)譜的詳解與介紹

#其他

【JY】位移角還是有害位移角?

【JY】如何利用python來編寫GUI?


【JY】主成分分析與振型分解的圖37

~路過別擦肩~


登錄后免費(fèi)查看全文
立即登錄
App下載
技術(shù)鄰APP
工程師必備
  • 項(xiàng)目客服
  • 培訓(xùn)客服
  • 平臺(tái)客服

TOP

4
1
1