• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于數(shù)值解的非線性振動系統(tǒng)定性分析

    2022-10-27 09:11:32張紅兵李萬祥
    振動與沖擊 2022年20期
    關鍵詞:龐加萊不動點平衡點

    陳 寧, 張紅兵, 李萬祥

    (蘭州交通大學 機電工程學院,蘭州 730070)

    在非線性振動系統(tǒng)的定性研究中,龐加萊映射是極其重要的思想工具,可以將連續(xù)微分動力學系統(tǒng)轉(zhuǎn)化離散差分動力學系統(tǒng),使之深化對分叉行為和過程的研究。在非線性動力學研究中,定量分析是定性分析的基礎,定量分析的結(jié)果是定性分析中龐加萊映射得以構(gòu)建的基礎。但是在當前的研究中,對復雜的、難以解析或不可解析的振動系統(tǒng),定性研究的深度遠不及對其的定量分析。這是由于對這類系統(tǒng)無法得到精確的解析解,而在現(xiàn)行的定性分析方法中,一個可以被分析的龐加萊映射必須依賴精確的解析解,而非數(shù)值解。

    非線性振動系統(tǒng)的定量分析方法可劃分為解析法和數(shù)值法兩個范疇,與之相對應的在定性分析過程中,理論上也具有依靠解析解和數(shù)值解兩種構(gòu)建龐加萊映射的基本方式。為便于下文闡述,將映射的可視化結(jié)果簡稱映像,將映射的迭代不動點簡稱為平衡點,分別將通過精確解析解、近似解析解和數(shù)值解構(gòu)架起來的龐加萊映射簡稱為解析映射、近似解析映射和數(shù)值映射。

    龐加萊映射只是一個特殊的運算,從本質(zhì)上來說依靠精確的解都可以構(gòu)建起系統(tǒng)的龐加萊映射。Luo等[1-3]考慮了分段線性系統(tǒng)可以分段解析的特點,利用解析解以及周期運動的假設構(gòu)建系統(tǒng)的周期映射,深刻地揭示了含間隙振動系統(tǒng)的分叉過程和機理。其中得出擦邊可能是導致分叉奇異的發(fā)現(xiàn)和結(jié)論,這既是非線性動力學領域持續(xù)關注的熱點之一[4-6],也是機械振動系統(tǒng)中普遍需要解決的問題之一。但是大多數(shù)的系統(tǒng)都不能利用解析法進行定量分析,如文獻[7-9]中針對含庫倫摩擦振動系統(tǒng)做出的研究,雖然可以分段解析,但是系統(tǒng)中由于間隙和摩擦力的兩個相互不干涉的因素存在,無法完全以解析的方法構(gòu)建系統(tǒng)的龐加萊映射,必須采用分支程序以輔助完成對系統(tǒng)分叉的分析。類似這種結(jié)合解析解和數(shù)值方法的定量分析方法可以統(tǒng)稱為半解析法,是針對復雜分段線性系統(tǒng)切實有效的方法,如李萬祥等[10]就單自由度非光滑系統(tǒng)采用半解析法進行研究,得出單自由度碰撞振動系統(tǒng)也可以存在Hopf分叉的發(fā)現(xiàn)。

    面對更加復雜的系統(tǒng),如分段光滑系統(tǒng)中,是難以分段精確解析的,陳洪月等[11]在對采煤機這類復雜設備的動力學響應研究中,就只能以數(shù)值模擬的方法求得系統(tǒng)的龐加萊截面圖??梢哉f,基于精確的解析解、半解析解和數(shù)值解都是可以構(gòu)建起系統(tǒng)的龐加萊映射。

    對于解析映射可以方便地利用不動點處雅可比矩陣的特征值信息,對系統(tǒng)中發(fā)生的分叉等做出精確的判別,如文獻[12]中對內(nèi)伊馬克沙克-音叉分岔點處的兩參數(shù)開折問題的研究。但是對于數(shù)值映射,由于缺少后續(xù)對應的分析方法,對其研究只能止步于觀察數(shù)值映像的分叉圖和結(jié)合經(jīng)驗做出分析,然后利用經(jīng)驗對系統(tǒng)中的分叉行為等做出解釋,如文獻[13]中關于采煤機和復雜齒輪系統(tǒng)的動力學特性研究,就以數(shù)值模擬的方法結(jié)合龐加萊截面法求得系統(tǒng)的數(shù)值映像及分叉圖,但是并未做到文獻[14]中對分叉分析的深刻成度,僅依靠觀察分岔圖從而得出部分規(guī)律。上述現(xiàn)象反映出數(shù)值映射的可分析性不足,不像解析映射那樣可以對龐加萊映射平衡點做出分析,這也是數(shù)值映射不被廣大學者所公認的原因。

    數(shù)值求解適合于所有振動系統(tǒng)的求解,這是解析分析所不具備的優(yōu)勢。雖然數(shù)值映射的概念尚未被認可,但是已有大部分學者在變相地使用,如文獻[15-19]中針對含摩擦系統(tǒng)、分段光滑系統(tǒng)、含間隙系統(tǒng)等復雜多體動力學系統(tǒng)的研究,通過大量龐加萊截面法來分析問題,其過程的本質(zhì)是龐加萊映射的數(shù)值化實現(xiàn)方法,其結(jié)果和解析映射的結(jié)果一樣為龐加萊映像圖。

    數(shù)值計算的適用性更強,應用面更廣,但是數(shù)值映射在定性研究中發(fā)揮的作用卻相對更小,這其中也包含學者對數(shù)值解精確性有所懷疑的因素。大量的理論分析和試驗證明,只要嚴格控制算法穩(wěn)定性,必定可以保證系統(tǒng)的積分精度。鑒于其求解精度,數(shù)值解是目前大多數(shù)近似解析方法精度驗證的標準。如劉汝逾等[20-22]在解決非線性振動問題時,就以數(shù)值解為標準對各自所提出的近似解析法的求解精度進行驗證。大量拓展近似解析法的論文,都以數(shù)值解驗證近似解析解精度,這說明數(shù)值解具有比近似解析解更高的可靠性。

    針對數(shù)值映射尚不能被分析研究的問題,部分學者嘗試采用近似解析解替代精確解析解建立系統(tǒng)的近似解析映射,但是目前還沒有實質(zhì)性的突破以便在實際中廣泛應用。雖然一個成熟的近似解析方法足以對系統(tǒng)做出有效的定量研究,但是函數(shù)映射的特點和近似解析解的本質(zhì)特征也決定利用近似解析映射可能導致定性分析結(jié)果失真的必然性。

    因此,考慮數(shù)值映射的可靠性和廣泛的適用性,如果能夠?qū)?shù)值映射找到與之匹配的分析過程和方法,必定可以深化對不可解析甚至難以解析振動系統(tǒng)的定性分析。

    分叉點處龐加萊映射迭代不動點的雅可比矩陣特征值是研究高余維分叉的關鍵,目前成熟的方法中缺乏對數(shù)值映射迭代不動點極其雅可比矩陣特征值的求解方法,這是導致對復雜系統(tǒng)定性分析深入度不足的主要原因,也將是本課題研究的主要工作。

    1 近似迭代不動點的數(shù)值計算方法

    系統(tǒng)龐加萊映射的迭代不動點是分叉行為分析的關鍵,因此利用數(shù)值映射的本質(zhì)特點尋找分叉點附近系統(tǒng)的迭代不動點是首要工作。

    1.1 關于分叉點處龐加萊映射的分析

    在系統(tǒng)分叉點的某個去心鄰域內(nèi)無論迭代過程是漸近收斂還是發(fā)散狀態(tài),其在時間進程中向穩(wěn)定狀態(tài)過渡的過程都十分緩慢。由于這種過程的緩慢性和穩(wěn)定性,通過自然迭代的方法求解不動點是不現(xiàn)實的,但是也為和其他方法的結(jié)合提供可能性和前提條件。

    在此定義某種基于數(shù)值解建立的離散動力系統(tǒng),由映射Py描述,如式(1)所示。

    ynyn+1=Py(μ,yn)yn

    (1)

    式中:yn+1為系統(tǒng)參數(shù)yn在映射Py作用下的象; 映射Py由原像yn和參數(shù)μ共同決定。

    設映射Py的迭代不動點為Y,根據(jù)不動點的映射條件應該滿足如式(2)所示的代數(shù)方程,則Y同樣的和參數(shù)μ相關。

    Y=Py(μ,Y)Y

    (2)

    參數(shù)yn可以表示為關于不動點Y和不動點的誤差εn之間的關系,如式(3)所示。將式(3)代入式(1),并將Py(μ,yn)按Taylor級數(shù)在不動點Y處展開,代入式(2)可以重構(gòu)與映射Py等價的映射Pε如式(4)所示, 映射Pε可以描述誤差參數(shù)εn和εn+1之間的映射關系,如式(5)所示。

    yn=Y+εn

    (3)

    (4)

    εnεn+1=Pε(μ,εn)εn

    (5)

    式中:εn為參數(shù)yn相對于不動點Y的誤差; 原點O即映射Pε的不動點,滿足式(6)所示的平衡關系。

    O=Pε(μ,O)O

    (6)

    當系統(tǒng)的不動點穩(wěn)定時,映射Pε在原點O處的譜半徑ρ[Pε(μ,O)]<1。假設系統(tǒng)在參數(shù)μ=μ0處發(fā)生某種基于不動點的經(jīng)典分叉,此時必然有ρ[Pε(μ0,O)]=1,即存在Pε的雅可比矩陣的特征值λ滿足如式(7)所示的形式,和分叉參數(shù)μ0有關。

    λj(μ0)=cosθ±sinθi∈diag[Pε(μ0,O)]

    (7)

    利用映射關系式(5)分析擾動εn的傳遞情況

    (8)

    假設存在轉(zhuǎn)化矩陣Qε(μ0,O)使Pε(μ0,O)可以對角化為Λε(μ0,O),滿足式(9)的形式。

    Qε(μ0,O)εn+1=Λε(μ0,O)Qε(μ0,O)εn

    (9)

    令Qε(μ0,O)εn=Zn,則映射Λε可以表示誤差εn正則化后的誤差參數(shù)Zn和Zn+1之間的映射關系,如式(10)所示。

    Zn+1=Λε(μ0,O)Zn

    (10)

    當條件式(7)成立時,初始擾動εn在某個空間方向或者空間曲面內(nèi)臨界穩(wěn)定,既不會收斂,也不會發(fā)散,保持相對初值的穩(wěn)定性,則式(10)在迭代過程中必然滿足式(11)的關系。

    ‖Zn+1‖=‖Zn‖

    (11)

    但是在一般情況下,由于計算的精度以及其他原因,只能夠確定分叉參數(shù)所在的某個鄰域,即僅僅可以找到某個μ∈U(μ0),可以表示為

    μ=μ0+μΔ

    (12)

    則特征值式(7)可表示為含有誤差的形式

    λj(μ)=[cos(θ+θΔ)±sin(θ+θΔ)i][1+Δλ]∈
    diag[Pε(μ,εn)]

    (13)

    式中: Δλ為由于參數(shù)μ的誤差引起的譜半徑的變化,Δλ=ρ[Pε(μ,εn)]-1;θΔ為特征值輻角的誤差; 此時若有μΔ→0, 則Δλ→0,(1+Δλ)→1,且使條件式(14)成立。

    ‖Zn+m-Zn‖=
    ‖Zn‖[1+Δλ]m-‖Zn‖?‖Zn‖

    (14)

    說明在有限次數(shù)的迭代之后,由于系統(tǒng)的發(fā)散和收斂造成的攝動的變化量要遠遠小于初始攝動量。在系統(tǒng)的某個維度面上,如系統(tǒng)的正則范式所對應的基平面,狀態(tài)參數(shù)的某些分量繞特定的點發(fā)生旋轉(zhuǎn),此時其平衡點近似在旋轉(zhuǎn)中心位置所在,但是由于其衰減指數(shù)十分接近1,故而其旋轉(zhuǎn)半徑的攝動量十分小。利用映射關系式(5)可知即便是收斂情況,通過持續(xù)迭代求取迭代不動點也是極其困難的,當?shù)l(fā)散時,必然是一種不可取的方案。

    1.2 算法原理和過程

    在自然迭代過程中,如果初始狀態(tài)的攝動量‖Zn‖變化速度較為緩慢,則在有限的時間內(nèi),‖Zn‖的值會保持相對的穩(wěn)定,可以借助這種相對穩(wěn)定特征來加速迭代收斂的過程。

    在此,可以借助分叉點處映射的特征進行迭代不動點的求解,如圖1所示。圖1中:O為其理論平衡點在系統(tǒng)特定低維度子空間中的投影;Zn,Zn+1,Zn+2為連續(xù)相鄰三次狀態(tài)構(gòu)成的龐加萊映像,近似滿足式(15)、式(16)的關系;Oμ為其外接圓心。

    (15)

    ∠ZnOZn+1=∠Zn+1OZn+2=θ+θΔ

    (16)

    可設龐加萊映像中狀態(tài)參數(shù)的坐標如式(17)所示。

    Zn[‖εn‖cos 0,‖εn‖sin 0],
    Zn+1[‖εn+1‖cos(θ+θΔ),‖εn+1‖sin(θ+θΔ)],
    Zn+2[‖εn+2‖cos(2θ+2θΔ),
    ‖εn+2‖sin(2θ+2θΔ)]

    (17)

    利用幾何關系可求得外接圓圓心Oμ坐標

    (18)

    計算的外接圓心的理論攝動半徑和初始狀態(tài)攝動半徑之比如式(19)所示。

    (19)

    當sinθ→0時,代入式(19)可得

    (20)

    由式(20)可以看出,當|sinθ|>|Δλ|時,經(jīng)過該方法計算后的攝動量是絕對小于自然迭代過程的,對平衡點失穩(wěn)的狀態(tài)同樣適用。初始狀態(tài)Zn經(jīng)過兩次自然迭代后的攝動半徑為‖Zn+2‖=‖Zn‖(1+Δλ)2,和自然迭代法相比,其收斂速度之比如式(21)所示。

    (21)

    當sinθ=0時,

    (22)

    同樣滿足當|sinθ|>|Δλ|時,該迭代計算過程是絕對收斂的,且效率絕對優(yōu)越于自然迭代。

    1.3 算法穩(wěn)定分析和修正

    式(22)也反映出對于sinθ=0的情況該算法效果不太理想,極有可能導致算法失穩(wěn)失效。對于這種情況,可以給出算法失效的判決方式和替代計算方法。

    在此以運算前的初始狀態(tài)和運算后的初始狀態(tài),分別經(jīng)過m次的自然迭代,則可以輕易求得以原初始狀態(tài)Zn為起點的迭代序列有限集合s1,以及以運算后的結(jié)果Zoμ為初始條件的迭代序列s2。

    s1: {Zn,Zn+1,Zn+2,Zn+3,…,Zn+m};

    s2: {Zoμ,Zoμ+1,Zoμ+2,Zoμ+3,…,Zoμ+m},Zoμ=Oμ。

    其集合S1和S2的構(gòu)建以外接矩形為例,如圖2所示。圖2和圖3中:黑點表示有限的龐加萊映像;矩形

    表示按照特定方向最小外接矩形構(gòu)造的凸集;圓點表示以原初始條件映射序列的有限元素集合s1;三角形表示運算后的值為初始狀態(tài)的映射序列有限元素集合s2;菱形表示最小凸集的幾何中心。

    以迭代計算前后的結(jié)果為起始點所構(gòu)造的兩個凸集的集合關系,如圖3所示。此時:若有S2∈S1,則算法必然穩(wěn)定;若有S1?S2,則可以判定為算法失穩(wěn)。當結(jié)果如圖3(a)所示時,運算可以連續(xù)壓縮平衡點所在區(qū)間,多次迭代后可得到有效的近似平衡點;當集合包含關系如圖3(b)所示時,算法處于失穩(wěn)狀態(tài);圖3(c)和圖3(d)所示情況說明,已經(jīng)達到計算機的或者程序求解的精度極限,或者原初始狀態(tài)已經(jīng)超出原平衡點的吸引域,并呈現(xiàn)指數(shù)發(fā)散狀態(tài)。對高維度的系統(tǒng),可以利用迭代映像在不同相平面上的投影來構(gòu)建所需要的兩個凸集,并結(jié)合圖3所對應的集合關系判斷算法的有效性。

    基于上述的特點,無論龐加萊映射發(fā)散與否只要轉(zhuǎn)換矩陣的特征值不存在大于1的正實數(shù),可以在判斷原算法失效的條件下,改用多次迭代的凸集幾何中心或者若干次龐加萊映射的重心替代原算法中的圓心Oμ,作為迭代運算的結(jié)果,如式(23)或式(24)所示。

    (23)

    (24)

    當離散樣本數(shù)量m較大時,理論上對各種基于周期運動發(fā)生的分叉都具有較好的收斂性和較為普遍的適應性。

    2 近似特征值的數(shù)值計算

    對平衡點的分析主要分析其雅可比矩陣的特征值,確定近似平衡點特征值也是十分重要的工作。可以采用第1章提供的迭代方法求得在系統(tǒng)Λε中的近似平衡點Zn。

    此時,可以保證‖Zn‖足夠小,假設存在若干初始狀態(tài)Zkn屬于列向量組Zsn,Zsn及其Zkn滿足式(25)所示的條件。

    (25)

    式中,E0=e0I,e0為對每個初始狀態(tài)Zn預設的初始誤差量,I為單位矩陣,且有rank(I)=rank(Zsn)。

    對初始狀態(tài)組Zsn和Zn取相同次數(shù)的龐加萊映射,可以得到關于近似不動點和初始狀態(tài)組的映射結(jié)果,Zsn+m及Zn+m,誤差關系可表示為

    (26)

    式中,Em為m次復合映射后的誤差矩陣,根據(jù)離散系統(tǒng)的穩(wěn)定性理論,可以利用式(27)確定雅克比矩陣的所有特征值。

    (27)

    方程γm=1在復數(shù)域內(nèi)具有m個解,其解集γ的表達式為

    (28)

    假設近似平衡點Zn處的特征值λ的真解如式(29)所示。

    λ=ρcosΦ±jρsinΦ

    (29)

    式中:ρ為特征值的模長,ρ=‖λ‖;Φ為特征值在復平面內(nèi)的角度。根據(jù)復數(shù)的乘積運算,特征值的m次冪可以表示為

    λm=ρm{cos[mod(mΦ,2π)]+
    jsin[mod(mΦ,2π)]}

    (30)

    通過對復數(shù)取根式可以得出的模接近真值,但是其角度信息并不一定正確。其中,當m=1時可以近似確定輻角位置,但是當m=1時的精確性較低,可以利用復數(shù)的指數(shù)運算規(guī)律,通過迭代法,分別運算m=[1,2,3,…]時的情況,利用漸近性的原理逐步逼近系統(tǒng)雅可比矩陣特征值的輻角真解,根據(jù)式(31)的迭代過程算法可以求得特征值輻角逼近解。

    (31)

    式中:Φm為m次誤差迭代對應的特征值輻角;Φ(m)為逐步逼近的不動點雅可比舉證的特征值輻角。

    3 算例分析

    3.1 算例模型

    以下以一種含間隙兩自由度彈簧非線性柔性碰撞振動系統(tǒng)在含三倍諧波的情況下的動力學研究分析為例,對其中某個類似Flip-NS型余維二的分叉過程進行研究。

    選取算例的力學模型如圖4所示。圖4中:M1和M2為振子的質(zhì)量;C1,C2和C為阻尼器的阻尼系數(shù);K1,K2和K為彈簧的剛度系數(shù);D為靜止狀態(tài)時的振子間的間隙;g為振子之間的實時動態(tài)間隙,負值表示壓縮量;F為外激勵;FD和-FD為由于彈簧的含間隙接觸緣故造成的作用與反作用力。由此可以建立系統(tǒng)的微分方程式(32)。

    (32)

    式中:彈簧剛度系數(shù)Kn具有非線性表達式(33);外激勵F含有高次諧波成分如式(34)所示;FD表達式如式(35)所示。

    (33)

    式中,χ為非線性剛度比例系數(shù)。

    (34)

    式中:P為無量綱的比例系數(shù);Pn為諧波強度;φn為n次諧波的相位。

    (35)

    由于間隙造成系統(tǒng)中含有一個非光滑切變界面ΣC,如式(36)所示,在其兩側(cè)系統(tǒng)的微分方程不連續(xù)。為保證系統(tǒng)狀態(tài)穿越界面ΣC的必然性,參考線性系統(tǒng)對間隙D做出約束,在此將彈簧力學性能曲線進行近似線性替代,并假設外激勵只存在主諧波,即理想化假設K1,2,3(Δ)=kn,F(xiàn)=P1cos(ω1t+φ1)。

    ΣC:X1-X2=D

    (36)

    假設系統(tǒng)位移的穩(wěn)定待定解如式(37)所示。

    (37)

    式中,A1,A2,B1和B1為解的待定系數(shù),將其代入方程并分離三角函數(shù)的相關項可得系數(shù)的平衡表達式(38)。

    (38)

    其中,

    在線性近似條件下系統(tǒng)振動時可穿越切變截面的充要條件如式(39)所示。

    (A1-A2)2+(B1-B2)2>D2

    (39)

    根據(jù)表達式(38)并定義間隙系數(shù)δ表征靜態(tài)間隙量D,滿足如式(40)所示的關系。

    (40)

    取模型參數(shù):M1=50;M2=4.5;k1=500;k2=1 244;k=2 244;χn=500;C1=0.3;C2=20.03;C=5.4;δ=0.007;ω=15.23;P1=0.128 0;P2=0;P3=0.154 4;φ1=0.24;φ2=0;φ3=0.63。

    以RK-4數(shù)值積分法進行仿真,采用穩(wěn)定的二分法進行求解穿過非光滑轉(zhuǎn)折界面的時間間斷點,求解精度取1×10-13。以系統(tǒng)主諧波外激勵周期定義龐加萊截面Σ如式(41)所示,在數(shù)值積分的過程中,利用龐加萊截面法構(gòu)建離散映射。

    Σ:mod(ωt,2π)=φ1

    (41)

    3.2 分叉分析

    通過數(shù)值仿真可以得到如圖5所示的系統(tǒng)分叉圖,具有雙穩(wěn)態(tài)的特征。圖5中:灰色分支的分叉過程比較明確,是單周期運動失穩(wěn)以激變的方式演化為擬周期運動,在擬周期區(qū)間發(fā)生的鎖相等行為也極其明確;黑色分支在P=12.5附近出現(xiàn)從擬周期演化為二倍周期運動的分叉過程,從分叉過程推斷,該分叉具有Flip-NS分叉的趨勢,但是由于本系統(tǒng)得不到系統(tǒng)的解析映射,故而無法利用解析映射的特征值來定義該分叉過程的類型。

    利用第1章提供的方法,可求得圖5中所對應兩支分叉過程的平衡點變化規(guī)律,如圖6所示。通過對比圖5和圖6可以看出,第1章所提供的方法對基于單周期運動的分叉具有良好的適應性,平衡點的激變過程和圖5中的激變過程相一致。圖中求解結(jié)果顯示,對基于多周期運動的分叉過程適應性不足,這是由于方法機制所導致,對于多周期運動,可以通過構(gòu)建復合映射的方式將其轉(zhuǎn)化為單周期運動,以便于研究。但是對于混沌運動,由于混沌運動的瀝遍性,該方法已經(jīng)失效。

    同時可以看出,對于圖5中存在疑點的分叉區(qū)間,圖6中對平衡點的求解是收斂且穩(wěn)定的,可以繼續(xù)利用第2章提出的迭代方法求得系統(tǒng)雅可比矩陣特征值的近似值。對平衡點的求解精度取1×10-13,然后對該近似平衡點的特征值進行37次迭代運算,求得系統(tǒng)的雅可比矩陣特征值變化如圖7所示。

    圖7中,當P=12.5時,不動點的雅可比矩陣其中一對共軛特征值處于單位圓外部,具有明顯擬周期運動的特點,和理論分析的結(jié)果相一致。隨著參數(shù)的遞增,其雅可比矩陣特征值的虛部分量逐漸減小為0,直接導致系統(tǒng)的龐加萊映像由擬周期狀態(tài)遷變?yōu)槎吨芷跔顟B(tài)??梢钥闯?,始終僅僅有一對特征值在單位圓附近變化,并不符合發(fā)生Flip-NS組合分叉時兩對特征值同時穿越單位圓的特點。

    上述分叉過程從分叉圖的變化規(guī)律分析,基本符合Flip-NS分叉的特征,但是從圖7內(nèi)特征值的變化規(guī)律中可以看出,其中一對復共軛特征值在單位圓外側(cè)退化為重實特征值,然后其中一個特征值從-1處穿越單位圓,該變化規(guī)律符合Hopf分叉周期二鎖相狀態(tài)的特征,其演化過程如圖8所示。

    4 結(jié) 論

    本文總結(jié)和分析了振動系統(tǒng)中不能利用數(shù)值解進行定性分析的現(xiàn)象及其原因。本質(zhì)原因是現(xiàn)行的定性分析方法對解析解的依賴;關鍵原因是無法利用數(shù)值解求得龐加萊映射的平衡點。文中對上述問題做出了針對性研究,總結(jié)了基于數(shù)值解求解平衡點及其龐加萊映射特征值的方法,并針對一個存在疑點的分叉行為成功地進行了分析。

    文中提到的方法在算例模型上的實踐被證明是可行有效的,非線性系統(tǒng)的數(shù)值解和非線性定性理論之間可以實現(xiàn)對接,利用差分和迭代原理基于數(shù)值解可以完成非線性系統(tǒng)的定性分析,擺脫在定性分析時對解析解的依賴。

    猜你喜歡
    龐加萊不動點平衡點
    龐加萊偵察術
    中外文摘(2022年7期)2022-05-17 09:36:42
    龐加萊偵察術
    青年文摘(2021年20期)2021-12-11 18:45:12
    龐加萊偵查術
    一類抽象二元非線性算子的不動點的存在性與唯一性
    活用“不動點”解決幾類數(shù)學問題
    探尋中國蘋果產(chǎn)業(yè)的產(chǎn)銷平衡點
    煙臺果樹(2019年1期)2019-01-28 09:34:58
    電視庭審報道,如何找到媒體監(jiān)督與司法公正的平衡點
    傳媒評論(2018年7期)2018-09-18 03:45:52
    在給專車服務正名之前最好找到Uber和出租車的平衡點
    IT時代周刊(2015年7期)2015-11-11 05:49:56
    不動點集HP1(2m)∪HP2(2m)∪HP(2n+1) 的對合
    推挽式直流變換器的龐加萊映射圖分析
    亚洲精品粉嫩美女一区| 欧美3d第一页| 欧洲精品卡2卡3卡4卡5卡区| 欧美潮喷喷水| 搞女人的毛片| 不卡一级毛片| 欧美日韩在线观看h| 久久精品国产亚洲av香蕉五月| 久久精品国产清高在天天线| 国产精品爽爽va在线观看网站| 日日摸夜夜添夜夜爱| 欧美日韩乱码在线| 国产精品日韩av在线免费观看| 欧美性猛交╳xxx乱大交人| 日韩成人av中文字幕在线观看 | 淫秽高清视频在线观看| 久久久国产成人免费| 免费看日本二区| 国产乱人偷精品视频| 看片在线看免费视频| 两性午夜刺激爽爽歪歪视频在线观看| 99久久九九国产精品国产免费| 婷婷精品国产亚洲av| 12—13女人毛片做爰片一| 色在线成人网| 床上黄色一级片| 观看美女的网站| 日韩高清综合在线| 成人av在线播放网站| 国产精品久久久久久久久免| 亚洲精品456在线播放app| 最近最新中文字幕大全电影3| 18+在线观看网站| 一个人看视频在线观看www免费| 九九在线视频观看精品| 国产午夜精品久久久久久一区二区三区 | 18+在线观看网站| 成年版毛片免费区| 无遮挡黄片免费观看| 搡女人真爽免费视频火全软件 | 欧美日韩乱码在线| 一区二区三区四区激情视频 | 一本久久中文字幕| 亚洲av电影不卡..在线观看| 国产精品野战在线观看| 伊人久久精品亚洲午夜| 久久久久久九九精品二区国产| 神马国产精品三级电影在线观看| 国产成年人精品一区二区| 乱人视频在线观看| 高清毛片免费观看视频网站| 国产精品伦人一区二区| 午夜免费男女啪啪视频观看 | a级毛片a级免费在线| 亚洲婷婷狠狠爱综合网| 久久午夜福利片| 国产美女午夜福利| 人妻丰满熟妇av一区二区三区| 亚洲欧美日韩东京热| av.在线天堂| 一a级毛片在线观看| 91在线精品国自产拍蜜月| 国产69精品久久久久777片| 尾随美女入室| 村上凉子中文字幕在线| 一区二区三区四区激情视频 | 丝袜喷水一区| 人妻制服诱惑在线中文字幕| 国产av在哪里看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲内射少妇av| 成人午夜高清在线视频| 国产亚洲欧美98| 淫妇啪啪啪对白视频| 国产精品国产高清国产av| 此物有八面人人有两片| 久久精品影院6| 久久韩国三级中文字幕| 99久国产av精品| 最近最新中文字幕大全电影3| 伦精品一区二区三区| 久久精品91蜜桃| 亚洲国产精品sss在线观看| 91av网一区二区| 精品午夜福利在线看| 亚洲国产精品成人综合色| 亚洲国产色片| 国产成人a区在线观看| 欧美不卡视频在线免费观看| av国产免费在线观看| 色哟哟·www| 性插视频无遮挡在线免费观看| 中文字幕精品亚洲无线码一区| 午夜久久久久精精品| 看免费成人av毛片| 久久精品夜夜夜夜夜久久蜜豆| 国产一区二区亚洲精品在线观看| 亚洲av免费高清在线观看| 亚洲久久久久久中文字幕| 国产午夜精品论理片| 两个人的视频大全免费| 长腿黑丝高跟| 美女cb高潮喷水在线观看| 晚上一个人看的免费电影| 看片在线看免费视频| 你懂的网址亚洲精品在线观看 | 午夜老司机福利剧场| 99热6这里只有精品| 啦啦啦观看免费观看视频高清| 一个人看的www免费观看视频| 午夜福利18| 欧美日本亚洲视频在线播放| 久久精品综合一区二区三区| 国产成人影院久久av| 一级av片app| 小蜜桃在线观看免费完整版高清| 成年女人看的毛片在线观看| 亚洲一级一片aⅴ在线观看| 成人午夜高清在线视频| 成人亚洲精品av一区二区| 亚洲无线观看免费| 不卡视频在线观看欧美| 搡老岳熟女国产| 有码 亚洲区| 最近手机中文字幕大全| 蜜桃久久精品国产亚洲av| 欧美色视频一区免费| 欧美国产日韩亚洲一区| 搞女人的毛片| 亚洲自拍偷在线| 精品久久久久久久久av| 午夜老司机福利剧场| 国产精华一区二区三区| 淫妇啪啪啪对白视频| 精品不卡国产一区二区三区| 亚洲美女视频黄频| 久久韩国三级中文字幕| 国内久久婷婷六月综合欲色啪| 亚洲欧美成人精品一区二区| 久久热精品热| 国国产精品蜜臀av免费| 人妻少妇偷人精品九色| 精品久久久久久久久久免费视频| 日本-黄色视频高清免费观看| 日韩制服骚丝袜av| 春色校园在线视频观看| 1024手机看黄色片| 国产探花在线观看一区二区| 久久久久免费精品人妻一区二区| 久久精品久久久久久噜噜老黄 | 噜噜噜噜噜久久久久久91| 久久草成人影院| 欧美区成人在线视频| 久久精品影院6| a级毛色黄片| 99久久九九国产精品国产免费| 丰满人妻一区二区三区视频av| 成人av在线播放网站| 亚洲性久久影院| 亚洲精品日韩av片在线观看| 亚洲aⅴ乱码一区二区在线播放| 蜜桃久久精品国产亚洲av| 女同久久另类99精品国产91| 国产一区二区在线av高清观看| 最近2019中文字幕mv第一页| 色5月婷婷丁香| 3wmmmm亚洲av在线观看| 亚洲四区av| 日韩一区二区视频免费看| 男人舔奶头视频| 激情 狠狠 欧美| 韩国av在线不卡| 伦精品一区二区三区| 少妇被粗大猛烈的视频| 国产成人影院久久av| 国产色爽女视频免费观看| 亚洲一区二区三区色噜噜| 日本精品一区二区三区蜜桃| 国产精品美女特级片免费视频播放器| 搞女人的毛片| 亚洲精品色激情综合| 久99久视频精品免费| 午夜免费男女啪啪视频观看 | 美女高潮的动态| 天天一区二区日本电影三级| 国产高清视频在线观看网站| 国产亚洲欧美98| 少妇高潮的动态图| 波多野结衣高清作品| 国产淫片久久久久久久久| 99热精品在线国产| 久久精品国产自在天天线| 一级毛片我不卡| 日韩一区二区视频免费看| 18禁黄网站禁片免费观看直播| 久久九九热精品免费| 亚洲,欧美,日韩| 高清日韩中文字幕在线| 干丝袜人妻中文字幕| 免费av不卡在线播放| 男女边吃奶边做爰视频| 精品人妻一区二区三区麻豆 | 99久国产av精品国产电影| 尾随美女入室| 精品久久久久久久久av| 一本精品99久久精品77| 3wmmmm亚洲av在线观看| 亚洲在线自拍视频| 日韩一本色道免费dvd| 免费人成视频x8x8入口观看| 午夜日韩欧美国产| 国产一区亚洲一区在线观看| 内地一区二区视频在线| 日韩欧美 国产精品| 国产黄色小视频在线观看| 国产色爽女视频免费观看| 亚洲乱码一区二区免费版| 中国美女看黄片| 国产伦精品一区二区三区四那| 美女大奶头视频| 亚洲av免费高清在线观看| 国产探花极品一区二区| 99热只有精品国产| 午夜激情欧美在线| 变态另类成人亚洲欧美熟女| av在线播放精品| 青春草视频在线免费观看| 亚洲欧美日韩无卡精品| 久久久色成人| av福利片在线观看| 成人午夜高清在线视频| 看片在线看免费视频| 久久婷婷人人爽人人干人人爱| 精品无人区乱码1区二区| 国产精品亚洲一级av第二区| 欧美日韩乱码在线| 久久久欧美国产精品| 久久久成人免费电影| 午夜激情欧美在线| 悠悠久久av| 国产成人a∨麻豆精品| 看黄色毛片网站| 精品一区二区免费观看| 午夜福利成人在线免费观看| 最好的美女福利视频网| 乱系列少妇在线播放| 免费无遮挡裸体视频| 国产精品三级大全| 国产在线精品亚洲第一网站| 成人国产麻豆网| 久久精品夜色国产| 麻豆久久精品国产亚洲av| 一个人观看的视频www高清免费观看| 精品久久久久久久人妻蜜臀av| 日韩制服骚丝袜av| 午夜精品国产一区二区电影 | 精品国产三级普通话版| 又爽又黄无遮挡网站| 国产高潮美女av| 卡戴珊不雅视频在线播放| 97超碰精品成人国产| 美女黄网站色视频| 99久久成人亚洲精品观看| 国产精品无大码| 国产精品一区www在线观看| 亚洲av二区三区四区| 乱人视频在线观看| 欧美潮喷喷水| 亚洲熟妇中文字幕五十中出| 亚洲不卡免费看| av在线观看视频网站免费| 国产欧美日韩一区二区精品| 亚洲第一区二区三区不卡| 国产精品国产三级国产av玫瑰| 啦啦啦啦在线视频资源| 人人妻,人人澡人人爽秒播| 亚洲精品国产av成人精品 | 日本欧美国产在线视频| 高清午夜精品一区二区三区 | 中文字幕免费在线视频6| 日日干狠狠操夜夜爽| 在线播放国产精品三级| 欧美bdsm另类| 国产单亲对白刺激| 亚洲激情五月婷婷啪啪| 日韩 亚洲 欧美在线| 欧美日韩一区二区视频在线观看视频在线 | 国产私拍福利视频在线观看| 午夜精品国产一区二区电影 | 深夜a级毛片| 国产片特级美女逼逼视频| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩乱码在线| 久久久欧美国产精品| 国产精品一及| 蜜桃久久精品国产亚洲av| 老熟妇仑乱视频hdxx| 春色校园在线视频观看| 99久久九九国产精品国产免费| 99九九线精品视频在线观看视频| 亚洲精华国产精华液的使用体验 | 麻豆国产av国片精品| 日本成人三级电影网站| а√天堂www在线а√下载| 91在线观看av| 噜噜噜噜噜久久久久久91| 99精品在免费线老司机午夜| 欧美精品国产亚洲| 久久久精品大字幕| 十八禁网站免费在线| 午夜久久久久精精品| 亚洲美女视频黄频| 色尼玛亚洲综合影院| 少妇被粗大猛烈的视频| 成人亚洲精品av一区二区| 久久人人精品亚洲av| 国产成人aa在线观看| 免费在线观看影片大全网站| 国产欧美日韩一区二区精品| 我的老师免费观看完整版| 日韩高清综合在线| 亚洲中文字幕日韩| 亚洲精品国产av成人精品 | 国产成人aa在线观看| 不卡一级毛片| 麻豆成人午夜福利视频| 成人毛片a级毛片在线播放| 18禁裸乳无遮挡免费网站照片| 一个人看的www免费观看视频| 日日干狠狠操夜夜爽| 国产精品亚洲美女久久久| 久久久精品94久久精品| 欧美高清成人免费视频www| 国产黄a三级三级三级人| 国产在线精品亚洲第一网站| 日本精品一区二区三区蜜桃| 黄色欧美视频在线观看| 国产亚洲91精品色在线| 国产精品电影一区二区三区| 久久久久久国产a免费观看| 精品国内亚洲2022精品成人| 亚洲熟妇中文字幕五十中出| 黄色欧美视频在线观看| 亚洲精品粉嫩美女一区| 波多野结衣巨乳人妻| 国产一区二区激情短视频| 久久久精品大字幕| 亚洲精品粉嫩美女一区| 亚洲人成网站在线观看播放| 成人特级av手机在线观看| 99热精品在线国产| 国产精品日韩av在线免费观看| 99精品在免费线老司机午夜| 国国产精品蜜臀av免费| 国产亚洲精品综合一区在线观看| 国产精品女同一区二区软件| 国产成人一区二区在线| 最近中文字幕高清免费大全6| 欧美人与善性xxx| 一级毛片久久久久久久久女| 久久鲁丝午夜福利片| 免费观看人在逋| 久久久久久久久久成人| 女人十人毛片免费观看3o分钟| 成人永久免费在线观看视频| 69人妻影院| 亚洲欧美日韩东京热| 精品免费久久久久久久清纯| 国内精品一区二区在线观看| 六月丁香七月| 国产亚洲精品久久久久久毛片| 精品久久久久久久末码| 国产熟女欧美一区二区| 欧美色视频一区免费| 我的老师免费观看完整版| 最近在线观看免费完整版| 国产极品精品免费视频能看的| 最近手机中文字幕大全| 国产午夜福利久久久久久| 美女cb高潮喷水在线观看| 欧美三级亚洲精品| 1024手机看黄色片| 日韩精品青青久久久久久| 狂野欧美白嫩少妇大欣赏| 男人舔奶头视频| 麻豆乱淫一区二区| 亚洲乱码一区二区免费版| 欧美中文日本在线观看视频| 91av网一区二区| 精品午夜福利视频在线观看一区| 九九久久精品国产亚洲av麻豆| 亚洲国产精品久久男人天堂| 在线免费观看的www视频| 久久人妻av系列| 亚州av有码| 亚洲在线自拍视频| 色综合站精品国产| 中国美白少妇内射xxxbb| 我要搜黄色片| 网址你懂的国产日韩在线| 欧美成人一区二区免费高清观看| 熟女电影av网| 亚洲国产色片| 国产男靠女视频免费网站| 性欧美人与动物交配| 国产一区二区三区av在线 | 婷婷亚洲欧美| 一本久久中文字幕| ponron亚洲| 国产中年淑女户外野战色| 亚洲真实伦在线观看| 夜夜夜夜夜久久久久| 一级a爱片免费观看的视频| 亚洲av电影不卡..在线观看| 国产精品一及| 成人欧美大片| 超碰av人人做人人爽久久| 国内精品一区二区在线观看| 97热精品久久久久久| 全区人妻精品视频| 国产精品一区www在线观看| 国产精品久久久久久亚洲av鲁大| 男人狂女人下面高潮的视频| 日本 av在线| 久久精品国产鲁丝片午夜精品| 丰满人妻一区二区三区视频av| 国产色婷婷99| 成人特级av手机在线观看| 亚洲性久久影院| 女人十人毛片免费观看3o分钟| av福利片在线观看| 精品一区二区三区视频在线| 自拍偷自拍亚洲精品老妇| 国产v大片淫在线免费观看| 日韩精品有码人妻一区| 国产aⅴ精品一区二区三区波| 亚州av有码| 免费一级毛片在线播放高清视频| 成人漫画全彩无遮挡| 看免费成人av毛片| 插阴视频在线观看视频| 大型黄色视频在线免费观看| 成人二区视频| 久久久久久久久大av| 亚洲三级黄色毛片| 黄色视频,在线免费观看| 久久久久国产网址| 插阴视频在线观看视频| 中文字幕av成人在线电影| 久久99热6这里只有精品| 小说图片视频综合网站| 日韩高清综合在线| 日韩一本色道免费dvd| 干丝袜人妻中文字幕| 国内揄拍国产精品人妻在线| 又黄又爽又刺激的免费视频.| 亚洲熟妇熟女久久| 亚洲成av人片在线播放无| 国产高清有码在线观看视频| 国产精品久久久久久久电影| 狂野欧美白嫩少妇大欣赏| 亚洲国产精品sss在线观看| 日本在线视频免费播放| 国产欧美日韩精品亚洲av| 欧美日本视频| 少妇猛男粗大的猛烈进出视频 | 久久精品国产清高在天天线| 看片在线看免费视频| 欧美日韩乱码在线| 九九在线视频观看精品| 国产中年淑女户外野战色| 精品久久久久久久久亚洲| 国产综合懂色| 狂野欧美激情性xxxx在线观看| 免费看光身美女| 高清毛片免费观看视频网站| 亚洲激情五月婷婷啪啪| 精品99又大又爽又粗少妇毛片| 免费不卡的大黄色大毛片视频在线观看 | 亚洲av.av天堂| 国内精品久久久久精免费| 乱人视频在线观看| 在线观看66精品国产| 精品一区二区三区视频在线| 国产 一区精品| 乱码一卡2卡4卡精品| 九色成人免费人妻av| 熟妇人妻久久中文字幕3abv| 久久久久国内视频| 亚洲av中文av极速乱| 亚洲欧美精品综合久久99| 97超视频在线观看视频| 亚洲成av人片在线播放无| 夜夜爽天天搞| 国产高潮美女av| 男人狂女人下面高潮的视频| 免费不卡的大黄色大毛片视频在线观看 | 99久久中文字幕三级久久日本| a级一级毛片免费在线观看| 综合色丁香网| 插逼视频在线观看| 三级毛片av免费| 午夜a级毛片| 免费人成在线观看视频色| 日韩欧美免费精品| av.在线天堂| 久久99热6这里只有精品| 国产蜜桃级精品一区二区三区| 成人一区二区视频在线观看| 日日干狠狠操夜夜爽| 亚洲中文字幕一区二区三区有码在线看| 亚洲成a人片在线一区二区| 国产高清激情床上av| 联通29元200g的流量卡| 两个人视频免费观看高清| 日本黄大片高清| 一级av片app| 国产成年人精品一区二区| 成人永久免费在线观看视频| 亚洲精品国产av成人精品 | 精品人妻一区二区三区麻豆 | 亚洲三级黄色毛片| 狠狠狠狠99中文字幕| 国产极品精品免费视频能看的| 麻豆乱淫一区二区| 我的女老师完整版在线观看| 在线播放无遮挡| 日韩欧美精品免费久久| 亚洲真实伦在线观看| 一本精品99久久精品77| 嫩草影院精品99| 亚洲专区国产一区二区| 麻豆乱淫一区二区| 国产精品,欧美在线| 尾随美女入室| 久久国内精品自在自线图片| 久久人人精品亚洲av| 观看免费一级毛片| 免费观看人在逋| 午夜福利在线观看免费完整高清在 | 国产久久久一区二区三区| 久久久久久久午夜电影| 亚洲专区国产一区二区| 欧美色欧美亚洲另类二区| 一进一出好大好爽视频| 欧美日韩国产亚洲二区| 一本久久中文字幕| 欧美日韩在线观看h| 亚洲在线观看片| 嫩草影院新地址| 国内精品美女久久久久久| 欧美成人精品欧美一级黄| 久久人人爽人人爽人人片va| 婷婷亚洲欧美| 少妇人妻精品综合一区二区 | 久久久精品94久久精品| 在线a可以看的网站| 国产亚洲欧美98| 亚洲18禁久久av| 国产精品乱码一区二三区的特点| 好男人在线观看高清免费视频| 国产一区二区在线观看日韩| 啦啦啦观看免费观看视频高清| 女人十人毛片免费观看3o分钟| 午夜福利高清视频| 高清毛片免费看| 久久久精品欧美日韩精品| 在线天堂最新版资源| 久久精品国产99精品国产亚洲性色| 国产精品人妻久久久影院| 国内精品一区二区在线观看| 亚洲人与动物交配视频| 高清日韩中文字幕在线| 久久精品国产亚洲av涩爱 | 成人精品一区二区免费| 最后的刺客免费高清国语| 国产真实伦视频高清在线观看| 亚洲av成人精品一区久久| 亚洲欧美日韩卡通动漫| 97碰自拍视频| av在线老鸭窝| 日日干狠狠操夜夜爽| 久久精品人妻少妇| 最近最新中文字幕大全电影3| 噜噜噜噜噜久久久久久91| 全区人妻精品视频| 精品日产1卡2卡| 淫秽高清视频在线观看| 精品熟女少妇av免费看| 如何舔出高潮| 成年女人永久免费观看视频| 亚洲成人av在线免费| 欧美中文日本在线观看视频| 亚洲精品粉嫩美女一区| 国产老妇女一区| 中出人妻视频一区二区| 久久久精品94久久精品| av.在线天堂| 国产91av在线免费观看| 综合色av麻豆| 国产成人精品久久久久久| 欧美高清性xxxxhd video| 国产av不卡久久| av天堂中文字幕网| 欧美丝袜亚洲另类| 国产成年人精品一区二区| 老熟妇乱子伦视频在线观看| 亚洲成a人片在线一区二区| 国产国拍精品亚洲av在线观看| 变态另类成人亚洲欧美熟女| 男人狂女人下面高潮的视频| 久久精品综合一区二区三区| 久久精品91蜜桃| 欧美成人精品欧美一级黄| 麻豆国产97在线/欧美|