• 一種基于宏纖維複合材料的拉索附加彎矩減振系統的制作方法

    文檔序號:11229615
    一種基于宏纖維複合材料的拉索附加彎矩減振系統的制造方法與工藝

    本發明涉及結構工程、自動控制技術領域,尤其涉及一種基于宏纖維複合材料的拉索附加彎矩減振系統。



    背景技術:

    隨著交通的日益發達,人們對于橋梁的跨度與舒適度也越來越高,因此斜拉橋與懸索橋對拉索的減振要求也越來越高。作爲斜拉橋與懸索橋的主要承重構件之一,拉索起著傳遞橋面荷載的重要作用,並且拉索的特點是阻尼小,柔度大和固有頻率頻段較寬。這些特點使得拉索極易在外荷載或橋面橋塔的作用下産生多種振動形式。由于拉索長期處于不同振動形式作用的組合或單獨作用下,拉索極可能發生疲勞破壞或者發生拉索錨固端破壞。拉索的疲勞破壞對于斜拉橋與懸索橋是致命的,這將造成嚴重的人員傷亡和經濟損失。

    關于拉索的減振策略,目前主要分爲三種,空氣動力學法、改變拉索自身特性與安裝機械阻尼器的方法。其中,通過改變拉索自身特性來減振的方法是一種易于實現的方法,一般爲輔助索法。雖然輔助索法易于布置,但輔助索減振效果受環境影響較大且不美觀,因此減小外界環境對拉索的影響與改變拉索自身特性相結合的方法是拉索振動控制的重要問題之一。

    宏纖維複合材料(MFC)是新型智能壓電纖維複合材料,主要由三層組成,包括上下兩層的交叉指型電極層與中間的壓電陶瓷纖維和聚合物基質。MFC具有質量輕、韌性大、出力大、耐腐蝕性好、易于粘貼與布置等特點,具有很好的發展前景。因此結合宏纖維複合材料的韌性大、出力大、耐腐蝕性好,拉索-MFC系統將會對拉索的振動控制起到極大的作用。



    技術實現要素:

    基于背景技術存在的技術問題,本發明目的之一在于提供一種方法:通過對拉索軸向同位置處對稱粘貼的MFC施加電壓,産生與外界荷載作用引起的彎矩方向相反的附加彎矩,減小拉索承受的彎矩,從而減小拉索的振幅。

    一種宏纖維複合材料的拉索,所述拉索的軸向設有若幹MFC組和位移傳感器,所述位移傳感器連接dSPACE實時仿真系統的輸入端,所述dSPACE實時仿真系統的輸出端連接高壓放大器的輸入端,所述高壓放大器的輸出端連接MFC組。

    優選的,所述MFC組是將若幹大小相同的MFC對稱地粘貼在拉索鋼絞線的周圍。

    優選的,所述dSPACE實時仿真系統與計算機連接。

    優選的,采用對同位置處對稱粘貼的MFC施加相反方向電壓的方法,從而使各MFC組産生相反方向的應變,並通過控制算法使MFC産生的附加彎矩與外界荷載産生的彎矩方向相反,構成拉索-MFC減振系統用于實時減小拉索的振幅

    一種基于宏纖維複合材料的拉索附加彎矩減振系統,方法步驟如下:

    S1:首先通過模態截斷法大致地截取所要控制的低階模態,以此作爲確定粘貼MFC與傳感器位置的預先估計;

    S2:建立拉索動力學方程與狀態空間方程組;

    S3:由于拉索是一種頻率較爲密集的柔性結構,然後選用內平衡降階方法,構造對角塊占優的系統可控、可觀Gram矩陣,求解系統的奇異值,對比奇異值大小,若沒有出現明顯大于其它奇異值的情況,就增加MFC數量,直到出現遠大于其它奇異值的奇異值時停止,保留奇異值較大的內平衡模態作爲系統保留模態。再通過步驟二的狀態空間方程與保留模態修正爲基于模態坐標的狀態空間方程;

    S4:編寫模糊PID控制算法,提取所截取模態的模態函數矩陣在傳感器位置處的模態函數建立新的模態函數矩陣,並將傳感器得到的物理空間信號與模態坐標信號和新的模態函數矩陣聯立得到方程組,這樣可以實時地求出模態坐標信號。由此實現了物理空間與模態空間的轉換。通過dSPACE實時仿真系統輸出控制信號,通過高壓放大器放大電壓作用于MFC,使MFC出力抵抗外界荷載作用,從而起到減振的作用。

    進一步的,所述S2中拉索動力學方程與狀態空間方程組的建立方法,方法步驟入下:

    設E爲拉索的彈性模量,I爲慣性矩,S0爲拉索初始索力,S爲振動引起的索力增量,η爲y方向的撓度,w爲振動引起的y方向撓度,qx爲忽略垂度影響的拉索沿水平跨長分布重量,t爲時間,T爲系統動能,V爲系統勢能,W爲外力做功,M(x)爲MFC産生的附加彎矩。

    因爲MFC對拉索的附加彎矩在拉索與MFC的橫截面邊界處有類似于突變的增長,且粘貼MFC的面積遠小于未粘貼位置,所以MFC在除粘貼位置外的部位産生的彎矩較小,對未粘貼部位振幅的不利影響可忽略不計。

    1)首先建立恒載作用下的靜力平衡方程

    EIη(4)-S0η″-qx=0

    2)推導拉索在附加彎矩時的動力方程

    根據哈密頓原理有

    所以

    上式的後兩項與邊界條件有關,其中邊界的δy與δy′均爲0。

    所以

    所以拉索的動力學方程爲

    3)將拉索的動力學方程寫爲狀態空間方程的形式。

    一種基于宏纖維複合材料的拉索附加彎矩減振系統統用于控制拉索的振幅。與現有技術相比,本發明具有的有益效果在于:

    本發明的系統提出的基于MFC的拉索附加彎矩減振方法,將多組MFC在拉索軸向同位置處對稱粘貼,通過對同位置處對稱粘貼的MFC施加方向相反的電壓,從而使各組MFC産生相反方向的應變,並通過控制算法使MFC産生的附加彎矩與外界荷載産生的彎矩方向相反,構成拉索-MFC減振系統以此實時地減小拉索的振幅。

    本發明提出的一種基于宏纖維複合材料的拉索附加彎矩減振系統,方便實施、便于觀測,提高了工作效率;該系統測試靈敏度高,能夠精准調節拉索的頻率。

    附圖說明

    附圖用來提供對本發明的進一步理解,並且構成說明書的一部分,與本發明的實施例一起用于解釋本發明,並不構成對本發明的限制。在附圖中:

    圖1爲本發明的基于MFC的拉索附加彎矩減振系統的布置示意圖,以斜拉索爲例;

    圖2爲本發明的MFC的粘貼位置橫截面示意圖;

    圖3爲本發明的基于MFC的拉索附加彎矩減振系統的流程示意圖,以斜拉索爲例;

    圖4爲基于模糊PID的拉索振動控制系統;

    圖5爲模糊PID與PID的控制效果對比圖;

    圖6爲一階頻率下拉索的控制效果圖;

    圖7爲二階頻率下拉索的控制效果圖。

    圖中:1-MFC組、2-拉索、3-拉索固定端、4-橋塔、5-橋面、6-位移傳感器、7-高壓放大器、8-計算機、9-dSPACE實時仿真系統、10-拉索鋼絞線、11-套管、12-填料、13-MFC。

    具體實施方式

    下面結合具體實施例對本發明作進一步解說。

    實施例1

    圖1表示的是基于MFC的拉索附加彎矩減振系統,在橋塔4和橋面5之間通過拉索固定端3連接拉索2,在拉索2上設有MFC組1和位移傳感器6,位移傳感器6連接dSPACE實時仿真系統9的輸入端,dSPACE實時仿真系統9的輸出端連接高壓放大器7的輸入端,高壓放大器7的輸出端連接MFC組1,dSPACE實時仿真系統9與計算機連接8。

    實施例2

    圖2表示的是MFC的橫向粘貼位置。在拉索鋼絞線10的周圍兩兩對稱粘貼四個MFC 13,拉索鋼絞線10的外部設有套管11,在拉索鋼絞線10與套管11之間還設有填料12。

    實施例3

    圖3爲基于MFC的拉索附加彎矩減振系統的流程示意圖,(1)首先要確定MFC與傳感器的粘貼位置,先通過模態截斷法大致地截取所要控制的低階模態,以此作爲確定粘貼MFC與傳感器位置的預先估計;(2)建立基于MFC的拉索附加彎矩的動力方程和狀態空間方程;(3)由于拉索是一種頻率較爲密集的柔性結構,然後選用內平衡降階方法,構造對角塊占優的系統可控、可觀Gram矩陣,求解系統的奇異值,對比奇異值大小,若沒有出現明顯大于其它奇異值的情況,就增加MFC與傳感器數量或改變MFC與傳感器的粘貼位置,直到出現遠大于其它奇異值的奇異值時停止,保留奇異值較大的內平衡模態作爲系統保留模態,修正由初始估計得到的狀態方程組,得到最終使用的狀態方程組,此方法解決了頻率密集的問題並確定了MFC與傳感器的最終粘貼位置;(4)最後編寫模糊PID控制算法,當有外荷載或橋塔橋面作用時,將傳感器得到的物理空間信號轉化爲模態坐標信號,再通過dSPACE實時仿真系統輸出控制信號,通過高壓放大器放大電壓作用于MFC,使MFC出力抵抗外界荷載作用,從而起到減振的作用。

    實施例4

    1)指定相關參數

    以如圖1所示的斜拉索爲模型,E爲拉索的彈性模量,I爲慣性矩,S0爲拉索初始索力,S爲振動引起的索力增量,η爲y方向的撓度,w爲振動引起的y方向撓度,qx爲忽略垂度影響的拉索沿水平跨長分布重量,t爲時間,M(x)爲MFC産生的附加彎矩,p(x,t)爲幹擾力,P0爲激勵力的幅值,θ爲激勵力的圓頻率,q(t)爲模態位移,V(x)爲振型曲線。且粘貼MFC的面積遠小于未粘貼面積,所以MFC在除粘貼位置外的部位産生的彎矩較小,對未粘貼部位振幅的不利影響可忽略不計。選用的MFC爲M8514-P1型MFC,可知其長度爲85mm,寬度爲14mm,厚度爲0.3mm,彈性模量爲30.336GPa,壓電常數d33爲400×10-12C/N,指叉電極間距爲0.35mm。本實例爲粘貼5對MFC,拉索的參數爲:長20m,直徑25mm,初始索力爲2000N,單位質量3.17kg/m,彈性模量爲130GPa。

    2)建立拉索動力學方程

    此處使用激振器進行單點激勵,假設P0爲100N,激振位置爲距原點位置5m處,所以取

    p(x,t)=P0δ(x-0.25L)sinθt

    3)列出初始條件與邊界條件並求解固有頻率

    假設初始條件與邊界條件爲

    v(x,0)=v0(x)

    v(0,t)=0

    v(L,t)=0

    所以可得

    w(x,t)=V(x)q(t)

    分離變量得

    S0V″-EIV(4)2mV=0

    求解可得

    V=A1sin(αx)+A2cos(αx)+A3sh(βx)+A4ch(βx)

    其中

    代入初始條件與邊界條件可得下式,求解下式得到各階頻率,由于下式爲超越方程,所以只能通過數值求解解出。

    2αβ(1-cos(αL)ch(βL))+(β22)sin(αL)sh(βL)=0

    計算得出前十階頻率

    f1=0.7176Hz、f2=1.4972Hz、f3=2.3900Hz、f4=3.4323Hz、f5=4.6484Hz

    f6=6.0537Hz、f7=7.6580Hz、f8=9.4679Hz、f9=11.4874Hz、f10=13.7195Hz

    4)求解各階振型

    由此可以求解拉索的前十階振型,可以看出拉索具有密頻的特性,即一階與二階頻率大小相差不大,因此選用內平衡降階技術,由此可得前十階振型。

    V1(x)=X1(sin(0.176x)-0.193sh(0.913x))-0.193X1(cos(0.176x)-ch(0.913x))

    V2(x)=X2(sin(0.349x)-0.363sh(0.961x))-0.363X2(cos(0.349x)-ch(0.961x))

    V3(x)=X3(sin(0.518x)-0.500sh(1.035x))-0.500X3(cos(0.518x)-ch(1.035x))

    V4(x)=X4(sin(0.683x)-0.606sh(1.126x))-0.606X4(cos(0.683x)-ch(1.126x))

    V5(x)=X5(sin(0.846x)-0.686sh(1.232x))-0.686X5(cos(0.846x)-ch(1.232x))

    V6(x)=X6(sin(1.007x)-0.747sh(1.348x))-0.747X6(cos(1.007x)-ch(1.348x))

    V7(x)=X7(sin(1.167x)-0.793sh(1.471x))-0.793X7(cos(1.167x)-ch(1.471x))

    V8(x)=X8(sin(1.326x)-0.829sh(1.600x))-0.829X8(cos(1.326x)-ch(1.600x))

    V9(x)=X9(sin(1.485x)-0.856sh(1.734x))-0.856X9(cos(1.485x)-ch(1.734x))

    V10(x)=X10(sin(1.643x)-0.878sh(1.871x))-0.878X10(cos(1.643x)-ch(1.871x))

    其中

    X1=0.1840、X2=0.1740、X3=0.1637、X4=0.1551、X5=0.1486

    X6=0.1438、X7=0.1402、X8=0.1376、X9=0.1354、X10=0.1337

    5)計算各階模態阻尼比

    由此可得到模態方程組,引入模態阻尼比

    寫成狀態空間方程

    其中

    B=diag(B1 B2 B3 … B10)

    模態阻尼依靠Rayleigh阻尼考慮

    假設前兩階模態阻尼比爲0.22%,則可得出a0與a1的值,進而求出剩余各階阻尼比,步驟如下

    由此得

    可得

    ξn=0.22%,0.22%,0.28%,0.37%0.48%,0.62%,0.77%,0.95%,1.15%1.37%

    6)拉索的模型降階

    通過內平衡降階技術降階,首先求解兩個李雅普諾夫方程中的Wc與W0,方程如下兩式

    AWc+WcAT+BBT=0

    ATW0+W0A+CTC=0

    然後在求解WcW0的特征值即爲奇異值的平方,奇異值表征了系統模態的可控性與可觀性,以此可以根據奇異值的大小來截取模態,截取奇異值較大的模態。根據奇異值的大小,這裏截取前兩階模態。奇異值如下:

    δ=5.602、5.577、1.287、1.281、0.394、0.392、0.145、0.144…

    7)計算模態幹擾力與模態控制力的大小

    a)模態幹擾力

    0.7176Hz的模態幹擾力爲

    F1(t)=35.868sin(62.832t)F2(t)=54.849sin(62.832t)

    1.4972Hz的模態幹擾力爲

    F1(t)=35.868sin(125.66t)F2(t)=54.849sin(125.66t)

    2Hz的模態幹擾力爲

    F1(t)=35.868sin(188.50t)F2(t)=54.849sin(188.50t)

    b)模態控制力

    宏纖維複合材料的力電公式爲

    由于宏纖維複合材料邊界處橫截面的正應力一定爲0,所以此處可將正應力的分布形式近似爲抛物線的形式,如下所示

    其中a與b分別爲MFC的坐標值,d爲拉索直徑。

    MFC附加彎矩的等效分布荷載公式如下:

    由此可知彎矩二階導數與電壓之間關系爲M″=3.967U,並可以計算模態控制力

    其中ai1與ai2爲第i片MFC的端部與尾部坐標值。

    本實例爲粘貼5對MFC爲了有效的控制這兩階模態,將其分別粘貼在距坐標原點5m、9m、10m、11m和15m處,其中5m、9m、10m、11m處MFC的出力與15m處MFC相反,這樣可以使兩個模態控制都不會過小。經計算如下二式,以此作爲對模態控制力的仿真限值。下面兩式爲

    FC1=0.19122U FC2=0.24611U

    傳感值的求取是通過將兩階模態位移值分別乘以各自所屬的振型函數得到拉索位移響應公式,再求傳感器坐標下的函數值。在模糊PID控制中,通過對傳感值與傳感值導數的模糊化,可以調節PID參數。拉索的位移響應公式i如下:

    w=V1(x)×q1(t)+V2(x)×q2(t)

    8)模糊PID控制計算

    本實例使用如下所示的控制系統求出拉索模態坐標下的響應,再將其轉化爲物理空間下的響應,本實例給出了在拉索坐標原點15m處,兩個模態頻率下的控制效果圖與2Hz作用下的模糊PID算法與PID算法的對比圖,如圖4所示。

    從圖4可以看出模糊PID控制優于PID控制,模糊PID控制控制效果更好且平穩,說明選擇模糊PID算法是合適的,從圖5、6可以看出MFC附加彎矩的方法有著良好的控制效果。

    以上所述,僅爲本發明較佳的具體實施方式,但本發明的保護範圍並不局限于此,任何熟悉本技術領域的技術人員在本發明揭露的技術範圍內,根據本發明的技術方案及其發明構思加以等同替換或改變,都應涵蓋在本發明的保護範圍之內。

    再多了解一些
    當前第1頁1 2 3 
    網友詢問留言 已有0條留言
    • 還沒有人留言評論。精彩留言會獲得點贊!
    1
    X-POWER-BY FNC Vbeta FROM ZZtest