国产成人嫩模一区二区,成人黄色小视频在线观看,国产成人夜夜嗨精品视频,成人无码A级毛片免费

 帶有非線性聯(lián)軸器軸系穩(wěn)態(tài)響應(yīng)計(jì)算方法的研究

5-1帶有非線性聯(lián)軸器軸系力學(xué)和數(shù)學(xué)模型的建立
在第三章扯饶，根據(jù)試驗(yàn)的基礎(chǔ)建立了具有非線性遲滯特性聯(lián)袖器的恢復(fù)力模型荡将。這一章，將研究一個(gè)帶有這種聯(lián)軸器的軸系如圖5-1所示岳悟，該軸系有n個(gè)圓盤，由鋼絲繩聯(lián)軸器與主機(jī)相連接，按以下原則建立力學(xué)模型：
1.每個(gè)圓盤均視為剛性勻質(zhì)，所有圓盤的質(zhì)量mi都有不同程度的偏心距ei颠放。轉(zhuǎn)軸軸線垂直通過各圓盤的幾何中心；
2.設(shè)軸承吭敢、軸承座以及聯(lián)軸器各向同性碰凶，靜坐標(biāo)系如圖5-1所示，支座處理成簡支鹿驼；
3.聯(lián)軸器從動端處理成一集中質(zhì)量m_b欲低，由于聯(lián)軸器主動端與主機(jī)軸相連接，主機(jī)袖相對軸系軸來說畜晰，較短較粗砾莱，剛性較大，變形較小舷蟀，故將聯(lián)軸器主動端恤磷、主機(jī)軸視為一體，位移為零戈勾，聯(lián)軸器從動端與主動端之間由非線性彈簧及非線性阻尼器聯(lián)結(jié)绷煎；
4.軸系為小位移振動，且忽略回轉(zhuǎn)效應(yīng)户痒。

當(dāng)主機(jī)以角速度ω轉(zhuǎn)動時(shí)火毕，軸系在各偏心力的作用下產(chǎn)生穩(wěn)態(tài)振動，其運(yùn)動微分方程可表示為：

式中 分別為軸系的質(zhì)量和左陣赠槽、阻尼矩陣和剛度矩陣剃炬。其中阻尼矩陣，a幢戳，b為比例常數(shù)莫诲。

其中：y_i,z_i分別為軸系中第i個(gè)圓盤在y,z方向的位移分量泰涡； =dy_i/dt， =d²y_i/dt²撇涡， =dz_i/dt像兆， =d²z_i/dt²分別為對應(yīng)于位移分量y_i,z_i的速度和加速度，y_b,z_b為聯(lián)軸器從動端集中質(zhì)量塊m_b在y,z方向的位移糜曲； =dy_b/dt河哑， =d²y_b/dt²， =dz_i/dt龟虎， =d²z_d/dt²分別為對應(yīng)于位移分量yb,zb的速度和加速度璃谨；,為y,z方向的激勵(lì)力向量，其中m_ie_iω²cos（ωt+ ）鲤妥，m_ie_iω²sin（ωt+ ）分別為y,z方向的偏心力分量佳吞； （ ，y_b棉安， 容达，ωt), （ ，z_b垂券， ，ωt）分別為非線性彈性聯(lián)軸器恢復(fù)力 在y,z方向的分量羡滑。 為第i個(gè)圓盤質(zhì)量偏心的初相位角菇爪，即軸系靜止時(shí)，第i個(gè)圓盤質(zhì)心和幾何中心連線與水平軸的夾角柒昏。
由于假設(shè)軸承凳宙、軸承座以及聯(lián)軸器各向同性，式（5-1〕轻调、式（5-2）解法相同尽诀，只需要討論二者之一即可。

式（5-l）表示一個(gè)局部非線性彈性和阻尼元件的振動系統(tǒng)皮匪，對這種系統(tǒng)陆宝，按現(xiàn)有常規(guī)方法來求解是非常困難的，為此涮观，本文以GILM為基礎(chǔ)發(fā)展了一種稱為SSGILM(Separate System-Gal-erkin and Improved Levenbery-Marquar-dt）的方法來求解此類微分方程組坑状。

5-2  SSGILM法

一.振動微分方程組的改寫和解耦
首先把有局部非線性系統(tǒng)振動微分方程組（5-1）式改寫成只有線性常系數(shù)的微分方程組和一個(gè)具有非線性變參數(shù)的微分方程兩部分：

式中：

式（5-4）的P_y中含有y_b和 而（5-5）的F_j中含有y_j，（j=1梳让，2却坦，… n）。這樣喻名，軸系分成運(yùn)動微分方程耦聯(lián)的兩個(gè)子系統(tǒng)一線性軸系子系統(tǒng)和非線性聯(lián)軸器子系統(tǒng)殃练。

式（5-4）為線性方程薇雳，按常規(guī)方法可求出無阻尼的各階固有頻率p_i以及對應(yīng)的主振型Y_i向量（i=1，2扫皱，…足绅，n），Y_i分別除以相應(yīng)廣義質(zhì)量的平方根 得到正則振型Y_Ni向量啸罢。引入正則振型坐標(biāo)W_Ni编检，對（5-4）式進(jìn)行坐標(biāo)變換；

式中扰才，W_Ni為正則振型向量W_N中的第i個(gè)分量允懂， 和 ，分別為它對時(shí)間t的一次導(dǎo)數(shù)和二次導(dǎo)數(shù)衩匣；I為單位矩陣蕾总； s_i=（a+ ）/2P_i=c_ii/2p_im_i為振型比例阻尼比；P_Ni為P_N= P_y激勵(lì)力向量中的第i個(gè)分量琅捏。
經(jīng)正則坐標(biāo)變換生百，雖然得到互不耦合的線性微分方程組（5-9)，但是（5-9）式仍然不能像單自由度振動系統(tǒng)那樣求解柄延，因?yàn)椋?-9）式中的激勵(lì)力包含未知的振動位移y_b和振動速度 蚀浆。如果y_b已知，就可以從（5-9）式中解得W_Ni搜吧，進(jìn)而可以由（5-7）式得到各y_i箍伏。因此需要求得y_b。
圖5-1所示軸系在主機(jī)帶動下轉(zhuǎn)動時(shí)乒踢，各圓盤質(zhì)量偏心將產(chǎn)生周期偏心激勵(lì)力揪馅，根據(jù)第二章的試驗(yàn)結(jié)果，聯(lián)軸器非線性恢復(fù)力Q是時(shí)間的周期函數(shù)捡霹，因此當(dāng)軸系中有聯(lián)軸器這種局部非線性元件時(shí)瞳州，可以設(shè)它的位移響應(yīng)是周期性的，即假設(shè)y_b有下面的形式：

由式（5-6）卧袄、（5-8）和（5-9）可知寒护，激勵(lì)力由二類力構(gòu)成，一類是質(zhì)量偏心力m_ie_iω²cos（ωt+ 岖疲，另一類是恢復(fù)力k_i（n+1）y_b-c_i（n+1） 咏摔，為求解（5-8）式方便，將（5-10）代入（5-8）式极阴，并將正則激勵(lì)力分解成：

P_N= P_y=P_N1+P_N2                                        （5-11）

式中：

由（5-11）和（5-12）式可知昙百，若{a}已知，就可以由（5-12）式求出各W_Ni碟狞，代W_Ni入（5-7）式啄枕，可以得到各y_i婚陪。由此可以把求y_b的問題轉(zhuǎn)化為求{a}=[a₀  a₁  的問題。

二频祝、Galerkin法的應(yīng)用

為了求出{a}泌参，引入無量綱時(shí)間τ= ，則（5-10）式可以寫成：

式中[q₀  q₁   … q_N  ]^T的各分量為聯(lián)軸器恢復(fù)力 的傅立葉級數(shù)各簡諧分量的系數(shù)分量常空。[f₁₀  f₁₁    …  f_1N    ]^T沽一，…[f_jo  f_j1   … f_jN  ]^T，…[f_n0  f_n1   … f_nN  ]^T分別為彈性恢復(fù)力和阻尼恢復(fù)力的F₁=k_（n+1） +c_（n+1） …漓糙，F(xiàn)_j=k_（n+1）+ c_（n+1） …铣缠，F(xiàn)_n=k_（n+1）+_（n+1） 的傅立葉級數(shù)各簡諧分量的系數(shù)向量。

由以上推導(dǎo)可知昆禽，求解帶有非線性聯(lián)軸器軸系的穩(wěn)態(tài)響應(yīng)蝗蛙，可以歸結(jié)為聯(lián)合求解2N+1個(gè)非線性代數(shù)方程組（5-15）和求解n個(gè)微分方程組（5-4)而重點(diǎn)又在于由（5-15）式求得2N

+1個(gè)未知量［a₀  a₁  ］^T。由(5-15）式可知醉鳖，它不同于一般的非線性方程組捡硅，在一般的非線性方程組中，未知量是顯含的剪炮，而（5-15）式的未知量{a}是隱含的讶瘩，而且{q｝是{a｝的非線件函數(shù)。不能通過將y_b的傅立葉展開式代人恢復(fù)力O的表達(dá)式中來直接求得藤门，這是因?yàn)镼表達(dá)式中的振幅 不能表示成{a}的顯函數(shù)的緣故俭整。因此要想由（5-1）式解得｛a}，還需要解決兩個(gè)問題：一是如何求 ,二是用什么來解隱含{a}的方程組（5-15)仙衩。下面就決這些問題。

三熔厌、ILM法

將非線性方程組簡寫成：

由于向量函數(shù){r（{a}）}中非線性項(xiàng){q（{a}）}的復(fù)雜性篇挡，可行的求解方法是采用數(shù)值法。

對于式（5-17）的數(shù)值法求解問題誉煎，我們將其轉(zhuǎn)化為與之等效的目標(biāo)函數(shù)：

的極小值問題报葬。即非線笥方程組（5-17）的解{a}可以通過求此極小值問題的最小二乘解得到。在此采用一種改進(jìn)的LM算法（ILM）來求解這一問題痛贬。ILM法的基本思想簡述如下：

當(dāng)LM法用于（5-18）問題時(shí)挺久，迭代公式為：

式中μ_k是LM參數(shù)，I為單位矩陣驻债，{p^（k）（μ_k）}為搜索逼近（5-17）式解的校正向量乳规，它依賴于參數(shù)μ_k值，控制搜索方向和長度合呐。

LM算法的優(yōu)點(diǎn)是暮的，在計(jì)算（5-20）式時(shí)笙以，可以不考慮[J（{a^（k）}）]是否滿秩，μ_k可以起到 （{a}）下降參數(shù)的作用冻辩，{p^（k）（μ_k）}是 （{a}）在{a^（k）}處的下降方向猖腕，而且總存在這樣的μ≥0，使 （{a^（k）}+ p^（k））< {a^（k）}）（對所有k）恨闪，隨著μ值的增大倘感，{p^（k）（μ_k）}的方向越來越接近 （{a^（k）}的負(fù)梯度方向，此時(shí)咙咽，收斂速度變慢老玛，但放寬了對初始近似{a^（0）}的要求。缺點(diǎn)是選擇μ滿足 （{a^（k+1）}）< （{a^（k）}）時(shí)犁珠，在一個(gè)迭代步驟中需要多次求解方程組（5-20）逻炊，大大增加了計(jì)算量，因此有必要對LM法進(jìn)行改進(jìn)雇牍，以使選擇μ_k時(shí)即不增加太多的計(jì)算量嗅呻，又能保證滿足下降性質(zhì)，即有一定的收斂速度乌心，又能改善（5-20）方程組的病態(tài)性导劝。改進(jìn)后的迭代公式為：

式中才 為矩陣 ^（k）的元素。 匹羹， 分別為矩陣L^（k）齐皂， ^（k）的元素。公式（5-23）和（5-27）和（5-28）中不含參數(shù)聲μ_k糕但，因此改變μ_k時(shí)仅汰，只需重新計(jì)算(5-24),而不需要重新進(jìn)行三角分解。顯然唠魏，改進(jìn)后的算法減少了由于選擇合適的μ_k而帶來的計(jì)算量蝉站。這種算法與LM法的不同之處在于用正定矩陣L_k 代替了單位矩陣I，它不僅調(diào)整了矩陣 ^（k）的主對角元素例吹，而且對整個(gè)矩陣進(jìn)行了調(diào)整捧颅，它比LM法有更好的收斂速度。
為了得到合適的μ_k值较雕，使得搜索向量{p^（k）（μ_k）}更快更好地滿足 （{a^（k）}+ p^（k））< {a^（k）}）碉哑，加速收斂，改善LM法對初始值{a^（0）}要求過高的缺點(diǎn)亮蒋，本文采用一種選擇合適μ_k值的新算法扣典，對LM法進(jìn)行了改進(jìn)，選擇計(jì)算步驟如下：

①.取定初始近似{a^（0）}和初始值 ，為了調(diào)整μ_k值激捏，引人調(diào)整系數(shù)v, 可以取得小一些设塔，取 =10^-2，v=5远舅。

②.按（5-21) 式首先算得［J（{a^（0）}）]闰蛔，然后按（5-26）式算得0⁽⁰⁾ ，再根據(jù)（5-27）图柏、（5-28）式算得矩陣L^(0）和 ^（0）中的各元素序六，按(5-25）、(5-23）和（5-24）算得{P^（0）（ ）}蚤吹；按（5-17）算得{r（{a^（0）}）}拒吧；按（5-18）算 （{a^（0）}）。并計(jì)算{a^（1）}={a^（0）}+{P^（0）（ ）}和 （{a^（1）}）
③.(l）如果少 （{a^（1）}）≤ （{a^（0）}）泡募，則以 /v代替 重新計(jì)算{P^（0）（ /v）}和新的{ ^（1）}描蹦， （{ ^（1）}）。（a）如果有 （{ ^（1）}）≤ （{a^（0）}）票虎，則取定μ₀= /v仆阶；（b）如果 （{a^（1）}）> （{a^（0）}），則取定μ₀= 摔色，（2） （{a^（1）}）> （{a^（0）}）驯祖，則需增大μ值，取 =v^m 数辱，以 代替 彪御，逐次對m=1，2初藐，…時(shí)的 代入（5-24）進(jìn)行計(jì)算筋顽，最后取定使不等式 （{a^（0）}）+{ P^（0）（ ）}）< （{a^（0）}）成立的最小正整數(shù) 對應(yīng)的 作為μ₀以及{a^（1）}+{P（ ）}。

④.假定已求得第k次近似{ a^（k）}和相應(yīng)的 =μ_k-1乍恐，算出[J（{ a^（k）}）]哪自， ^（k），L^（k）禁熏， ^（k），r（{ a^（k）}）}邑彪， （{ a^（k）}）瞧毙，然后算得{P（ ）}，并計(jì)算{ a^（k+1）}={ a^（k）}+{P（ ）}和 （{ a^（k+1）}）寄症。

⑤.（1）如果 （{ a^（k+1）}）≤ （{ a^（k）}）宙彪，則減少 以/V代替 重新計(jì)算{P（/V ）}和對應(yīng)的{ ^（k+1）}及 （{ ^（k+1）}）（a）如果 （{ ^（k+1）}）≤ （{ a^（k）}），則取定μ_k=/V ；（b）如果 （{ ^（k+1）}）> （{ a^（k）}）释漆，則取定μ_k= （2）如果 （{ ^（k+1）}）> （{ a^（k）}）悲没，則增大 值，可取 =v^m ,以 代替 ,逐次時(shí)m=1男图，2示姿，…時(shí)的 代入（5-24）式進(jìn)行計(jì)算，最后取定使不定式 ({ a^（k）}+{P（ ）}< （{ a^（k）}）揭凭。

⑥.目標(biāo)函數(shù) （{a}）的k次迭代和k+1次迭代如果滿足絕對誤差和相對誤差精度源快，即可將{ a^（k+1）}作為近似解。如果不滿足俯重，則以{ a^（k+1）}代替{ a^（k）}炫瘤，μ_k代替μ_k-1去執(zhí)行步驟⑤，直到滿足要求為止夕锹。

至此枷配，（5-17）式的求解方法已經(jīng)建立起來，日（5-17）式求解{a}還需要求出 檬桅，{q}和（f_j）（j=1物少，2，…绿捶，n）棚疏。另外，在用（5-21）式求[J（{a}）]時(shí)遇到r_i（{a}）對隱含向量{a}求偏導(dǎo)數(shù)的麻煩婆精，因此在計(jì)算[J（{a}）]時(shí)继溯，需要做些處理。由（5-15）式可知沈条，它是頻域內(nèi)的2N+1個(gè)非線性代數(shù)方程組需忿，而非線性遲滯恢復(fù)力 是時(shí)域函數(shù)，因此在求解（5-15）時(shí)需要進(jìn)行頻域與時(shí)域相互變換蜡歹，并做一些數(shù)學(xué)上的處理屋厘，下面逐一敘述。

5-3一步用SSGILM法計(jì)算帶有非線笥聯(lián)軸器軸系穩(wěn)態(tài)響應(yīng)

為了求得未知向量{a}月而，還需做以下工作：

一汗洒、向量函數(shù){r（{ a^（k）}）}的計(jì)算

要計(jì)算{r（{ a^（k）}）}為聯(lián)軸器恢復(fù)力 的傅立葉級數(shù)展開式各簡諧分量的系數(shù)向量。對于給定的{ a^（k）}父款，不能由時(shí)域表達(dá)式求得頻域內(nèi)的{q（{ a^（k）}）}溢谤，而必須按以下步驟計(jì)算。

首先憨攒，由給定的頻域向量{q（{ a^（k）}）}世杀，根據(jù) 的N個(gè)時(shí)域離散值，記為Y（1）～Y（N）和DY（1）～DY（N），其次瞻坝，由下式求得 幅值

然后蛛壳，將 以及 時(shí)域離散值代入第三章中得到的恢復(fù)力 的表達(dá)式（3-8）、（3-12）牢簸、（3-27）和（3-41）式骨矗，算得對應(yīng)于{ a^（k）}的恢復(fù)力Q在一個(gè)周期內(nèi)N個(gè)Q的時(shí)域離散值進(jìn)行FFT變換，即可得到非線性恢復(fù)力Q在頻域內(nèi)的各諧波分量的系數(shù)向量{q（{ a^（k）}）}

2. {f_j（{ a^（k）}）}的計(jì)算

{f_j（{ a^（k）}）}是式（5-6）中恢復(fù)力F_j的傅立葉級數(shù)展開式各諧波分量的系數(shù)向量碘淘。當(dāng)向量{ a^（k）}給定時(shí)刮血，首先由5-2節(jié)中所述方法算得各個(gè) 和 （j=1，2惰采，…颂鞭，n），其次由（5-6）式中恢復(fù)力F_j的各諧波分量的系數(shù)向量{f_j（{ a^（k）}）}栋锣。

算得{q（{ a^（k）}）}和{f_j（{ a^（k）}）}后刑袒，代入（5-15）式即可得到{r（{ a^（k）}）}。

二谍竿、雅可比矩陣[J（{ a^（k）}）]的計(jì)算

為計(jì)算方便起見程悼，將雅可比矩陣式（5-21）分解成三個(gè)矩陣：

[J（{ a^（k）}）]=[J₁]+[J₂]+[J₃]                  （5-30）

1.[J₁]的計(jì)算

由（5-15）和（5-21）式可得：

[J₁]=

2. [J₂]的計(jì)算

[J₂]的計(jì)算實(shí)際上就是對 {q（{ a^（k）}）}/ { a^（k）}的計(jì)算。由于恢復(fù)力Q的函數(shù)表達(dá)式是{ a^（k）}的隱函數(shù)路影，因此不可能由Q的表達(dá)式直接進(jìn)行FFT變換來求得{q（{ a^（k）}）}诬像，也可能由{q（{ a^（k）}）}直接對{ a^（k）}；求偏導(dǎo)數(shù)來得到[J₂]：

由第三章中的式（3-8）闸婴、（3-12）坏挠、（3-41）和本章中的式（5-13）可得：

式中， 邪乍，降狠， ，庇楞， 可以直接由Q的表達(dá)式和（5-13）進(jìn)行解析計(jì)算得到榜配， / { a^（k）}， / { a^（k）}也可以由（5-13）式直接進(jìn)行解{ a^（k）}吕晌，選取增量△{ a^（k）}蛋褥，然后計(jì)算一個(gè)周期內(nèi)對應(yīng)的 （{ a^（k）}+△{ a^（k）}）和 （{ a^（k）}）值，得：

由（5-32）式算得 / { a^（k）}的時(shí)域離散值后睛驳，進(jìn)行FFT變換即可得到頻域中的 {q（{ a^（k）}）}/ { a^（k）}烙心。

3. [J₃]的計(jì)算

[J₃]的計(jì)算實(shí)際上就是偏導(dǎo)數(shù) {f_j（{ a^（k）}）}/ { a^（k）}的計(jì)算。計(jì)算方法與計(jì)算 {q（{ a^（k）}）}的方法類似柏靶，故不再贅述。對于[J₂]，[J₃]铛田、分別有式（5-34）挂蹦、（5-35）。

在計(jì)算出{r（{ a^（k）}）}和{J（{ a^（k）}）}之后嚣历，就可以按5-2節(jié)中發(fā)展的ILM法進(jìn)行搜索迭代計(jì)算叼叙，由初始向量{a^（0）}不斷搜索逼近滿足精度要求的{ a^（k+1）}，然后代入（5-10）盟挤、（5-11）么酪、（5-12）和（5-7）即可求得聯(lián)軸器從動動端質(zhì)理塊m_b的穩(wěn)態(tài)響應(yīng)y_b以及軸系中各圓盤轉(zhuǎn)子穩(wěn)態(tài)應(yīng)Y。同時(shí)還可以得到對應(yīng)的幅頻圖腻俏。

5-4小結(jié)
本章對帶有非線性聯(lián)抽器的n+1個(gè)自由度軸系在四盤質(zhì)量偏心產(chǎn)生的激勵(lì)力作用下穩(wěn)態(tài)響應(yīng)計(jì)算方法進(jìn)行了研究长尼，以GLM法為基礎(chǔ)，提出了一種稱為SSG

ILM (Separate System-Ga1erkin and Imp-r0ved Leenberg-Marquardt的方法目尘，用于計(jì)算局部具有非線性特性軸系的穩(wěn)態(tài)響應(yīng)瞻替，
主要工作如下：
1.建立帶有非線性聯(lián)鈾器n+1個(gè)自由度軸系的止）學(xué)模型和數(shù)學(xué)模型，為了求解這種局部含有非線性元件的軸系的穩(wěn)態(tài)響應(yīng)诺秒，將軸系分成兩個(gè)子系統(tǒng)遥倦，即含有n個(gè)自由度的線性軸系和含有非線性變參數(shù)的聯(lián)軸器子系統(tǒng)。
2.對線性軸系子系統(tǒng)運(yùn)動微分方程組用正則坐標(biāo)變換進(jìn)行解藕處理占锯；對非線性變參數(shù)的聯(lián)軸器子系統(tǒng)運(yùn)動微分方程用Galerkin法袒哥，變成一組隱含系數(shù)向量{a}的非線性代數(shù)方程組。
3.針對LM法的缺點(diǎn)消略，采用已有的選擇合適μ值的方法和加快收斂速度的算法對LM法進(jìn)行了改進(jìn)堡称，用改進(jìn)算法對上面得到的非線生方程組中的未知向量{a}進(jìn)行搜索計(jì)算，最后解出局都具有非線性特性軸系的穩(wěn)態(tài)響應(yīng)疑俭。