混合監測受損索集中載荷廣義位移識別方法

            文檔序號:6187423閱讀:120來源:國知局
            混合監測受損索集中載荷廣義位移識別方法
            【專利摘要】混合監測受損索集中載荷廣義位移識別方法基于混合監測,通過監測索結構溫度和環境溫度來決定是否需要更新索結構的力學計算基準模型,得到計入索結構溫度和環境溫度的索結構的力學計算基準模型,在此模型的基礎上計算獲得單位損傷被監測量數值變化矩陣。依據被監測量當前數值向量同被監測量當前初始數值向量、單位損傷被監測量數值變化矩陣和待求的被評估對象當前名義損傷向量間存在的近似線性關系算出被評估對象當前名義損傷向量的非劣解,據此可以在有溫度變化時,識別出支座廣義位移、集中載荷變化量和受損索。
            【專利說明】混合監測受損索集中載荷廣義位移識別方法
            【技術領域】
            [0001]斜拉橋、懸索橋、桁架結構等結構有一個共同點,就是它們有許多承受拉伸載荷的部件,如斜拉索、主纜、吊索、拉桿等等,該類結構的共同點是以索、纜或僅承受拉伸載荷的桿件為支承部件,為方便起見,本方法將該類結構表述為“索結構”,并將索結構的所有承載索、承載纜,及所有僅承受軸向拉伸或軸向壓縮載荷的桿件(又稱為二力桿件),為方便起見統一稱為“索系統”,本方法中用“支承索”這一名詞指稱承載索、承載纜及僅承受軸向拉伸或軸向壓縮載荷的桿件,有時簡稱為“索”,所以在后面使用“索”這個字的時候,對桁架結構實際就是指二力桿件。在結構服役過程中,對支承索或索系統的健康狀態的正確識別關系到整個索結構的安全。在環境溫度發生變化時,索結構的溫度一般也會隨著發生變化,在索結構溫度發生變化時,在索結構服役過程中,索結構支座可能發生廣義位移,索結構承受的集中載荷也可能發生變化,同時索結構的健康狀態也可能在發生變化,在這種復雜條件下,本方法基于混合監測(本方法通過對本節前述不同類型的索結構的可測量參數的變化的混合監測來判斷索結構的健康狀態,本方法將所有被監測的索結構特征參量統稱為“被監測量”,由于此時被監測量是由索結構的不同類型的可測量參數混合組成,本方法稱此為混合監測)來識別支座廣義位移、受損索和索結構承受的集中載荷的變化量,屬工程結構健康監測領域。
            【背景技術】
            [0002]剔除載荷變化、索結構支座廣義位移和結構溫度變化對索結構健康狀態識別結果的影響,從而準確地識別結構的健康狀態的變化,是目前迫切需要解決的問題;剔除載荷變化、索結構健康狀態變化和結構溫度變化對索結構支座廣義位移識別結果的影響,從而準確地識別索結構支座廣義位移,也是目前迫切需要解決的問題;同樣的,剔除結構溫度變化、索結構支座廣義位移和結構健康狀態變化對結構承受的集中載荷的變化量的識別結果的影響,對結構安全同樣具有重要意義。基于結構健康監測技術,本方法公開了一種解決這三個問題的有效方法。
            [0003]當索結構承受的集中載荷出現變化時、或索結構發生支座廣義位移、或索結構的溫度發生變化時、或索系統的健康狀態發生變化(例如發生損傷)時、或者四種情況同時發生時,會引起索結構的可測量參數的變化,例如會引起索力的變化,會影響索結構的變形或應變,會影響索結構的形狀或空間坐標,會引起過索結構的每一點的任意假想直線的角度坐標的變化(例如結構表面任意一點的切平面中的任意一根過該點的直線的角度坐標的變化,或者結構表面任意一點的法線的角度坐標的變化),所有的這些變化都包含了索系統的健康狀態信息,也包含了支座廣義位移和集中載荷的變化量信息,也就是說可以利用索結構的可測量參數來識別支座廣義位移、受損索和集中載荷的變化量。
            [0004]在支座有廣義位移時,目前已公開的技術、方法中,有些僅僅能夠在其它所有條件不變時(僅僅只有結構承受的載荷發生變化)識別結構承受載荷的變化,有些僅僅能夠在其它所有條件不變時(僅僅只有結構健康狀態發生變化)識別結構健康狀態的變化,有些僅僅能夠在其它所有條件不變時(僅僅只有結構溫度和結構健康狀態發生變化)識別結構(環境)溫度和結構健康狀態的變化,目前還沒有一種公開的、有效的方法能夠同時識別結構承受載荷、結構(環境)溫度和結構健康狀態的變化,或者說在結構所承受的載荷和結構(環境)溫度同時變化時,還沒有有效的方法能夠識別結構健康狀態的變化,而結構承受的載荷和結構(環境)溫度是常常變化的,所以如何在結構承受的載荷和結構(環境)溫度變化時,剔除載荷變化和結構溫度變化對索結構健康狀態識別結果的影響,從而準確地識別結構的健康狀態的變化,是目前迫切需要解決的問題,本方法公開了一種方法,在支座有廣義位移時,可以在索結構承受的集中載荷和結構(環境)溫度發生變化時,剔除支座廣義位移、載荷變化和結構溫度變化對索結構健康狀態識別結果的影響,基于被監測量監測來識別受損索,對索結構的安全具有重要的價值。
            [0005]同樣的,在目前公開的方法中,還沒有出現能夠剔除支座廣義位移、結構溫度變化和支承索健康狀態影響的、從而實現集中載荷變化程度的正確識別的方法,而對結構來說,載荷變化的識別也是非常重要的;本方法在識別出受損索的同時,還能同時識別出集中載荷的變化,即本方法能夠剔除支座廣義位移、結構溫度變化和支承索健康狀態變化的影響,實現集中載荷變化程度的正確識別;本方法還能夠剔除載荷變化、索結構健康狀態變化和結構溫度變化對索結構支座廣義位移識別結果的影響,從而準確地識別索結構支座廣義位移。
            [0006]也就是說,本方法實現了已有方法不可能具備的三種功能。

            【發明內容】

            [0007]技術問題:本方法公開了一種方法,實現了已有方法不可能具備的三種功能,分別是,一、剔除支座廣義位移、集中載荷變化和結構溫度變化對索結構健康狀態識別結果的影響,從而準確地識別出支承索的健康狀態;二、本方法在識別出受損索的同時,還能同時識別出集中載荷的變化,即本方法能夠剔除支座廣義位移、結構溫度變化和支承索健康狀態變化的影響,實現集中載荷變化程度的正確識別;三、本方法還能夠剔除載荷變化、索結構健康狀態變化和結構溫度變化對索結構支座廣義位移識別結果的影響,從而準確地識別索結構支座廣義位移。
            [0008]技術方案:本方法由三部分組成。分別是建立結構健康監測系統所需的知識庫和參量的方法、基于知識庫(含參量)和實測被監測量的結構健康狀態評估方法、健康監測系統的軟件和硬件部分。
            [0009]在本方法中,用“支座空間坐標”指稱支座關于笛卡爾直角坐標系的X、Y、Z軸的坐標,也可以說成是支座關于X、Y、Z軸的空間坐標,支座關于某一個軸的空間坐標的具體數值稱為支座關于該軸的空間坐標分量,本方法中也用支座的一個空間坐標分量表達支座關于某一個軸的空間坐標的具體數值;用“支座角坐標”指稱支座關于X、Y、Z軸的角坐標,支座關于某一個軸的角坐標的具體數值稱為支座關于該軸的角坐標分量,本方法中也用支座的一個角坐標分量表達支座關于某一個軸的角坐標的具體數值;用“支座廣義坐標”指稱支座角坐標和支座空間坐標全體,本方法中也用支座的一個廣義坐標分量表達支座關于一個軸的空間坐標或角坐標的具體數值;支座關于X、Y、Z軸的坐標的改變稱為支座線位移,也可以說支座空間坐標的改變稱為支座線位移,本方法中也用支座的一個線位移分量表達支座關于某一個軸的線位移的具體數值;支座關于X、Y、Z軸的角坐標的改變稱為支座角位移,本方法中也用支座的一個角位移分量表達支座關于某一個軸的角位移的具體數值;支座廣義位移指稱支座線位移和支座角位移全體,本方法中也用支座的一個廣義位移分量表達支座關于某一個軸的線位移或角位移的具體數值;支座線位移也可稱為平移位移,支座沉降是支座線位移或平移位移在重力方向的分量。
            [0010]首先確認索結構承受的可能發生變化的集中載荷的數量。根據索結構所承受的集中載荷的特點,確認其中“所有可能發生變化的集中載荷”,或者將所有的集中載荷視為“所有可能發生變化的集中載荷”,設共有JZW個可能發生變化的集中載荷。[0011]集中載荷分為集中力和集中力偶兩種,在坐標系中,例如在笛卡爾直角坐標系中,一個集中力可以分解成三個分量,同樣的,一個集中力偶也可以分解成三個分量,在本方法中將一個集中力分量或一個集中力偶分量稱為一個集中載荷。
            [0012]設索結構的支承索的數量、索結構的支座廣義位移分量的數量和JZW個“所有可能發生變化的集中載荷”的數量之和為N。為敘述方便起見,本方法統一稱被評估的支承索、支座廣義位移和“所有可能發生變化的集中載荷”為“被評估對象”,共有N個被評估對象。給被評估對象連續編號,該編號在后續步驟中將用于生成向量和矩陣。
            [0013]被監測的多類參量可以包括:索力、應變、角度和空間坐標,分別敘述如下:
            [0014]設索系統中共有Q根支承索,索結構的被監測的索力數據由索結構上M1個指定索的M1個索力數據來描述,索結構索力的變化就是所有指定索的索力的變化。每次共有M1個索力測量值或計算值來表征索結構的索力信息。M1是一個不小于0的整數。
            [0015]索結構的被監測的應變數據可由索結構上K2個指定點的、及每個指定點的L2個指定方向的應變來描述,索結構應變數據的變化就是K2個指定點的所有被測應變的變化。每次共有M2 (M2=K2XL2)個應變測量值或計算值來表征索結構應變。M2是一個不小于0的整數。
            [0016]索結構的被監測的角度數據由索結構上K3個指定點的、過每個指定點的L3個指定直線的、每個指定直線的H3個角度坐標分量來描述,索結構角度的變化就是所有指定點的、所有指定直線的、所有指定的角度坐標分量的變化。每次共有M3 (M3=K3XL3XH3)個角度坐標分量測量值或計算值來表征索結構的角度信息。M3是一個不小于0的整數。
            [0017]索結構的被監測的形狀數據由索結構上K4個指定點的、及每個指定點的L4個指定方向的空間坐標來描述,索結構形狀數據的變化就是K4個指定點的所有坐標分量的變化。每次共有M4 (M4=K4XL4)個坐標測量值或計算值來表征索結構形狀。M4是一個不小于0的整數。
            [0018]綜合上述被監測量,整個索結構共有M (M=Mi+M2+M3+M4)個被監測量,定義參量K(K=M1+K2+K3+K4), K 和 M 不得小于 No
            [0019]為方便起見,在本方法中將“索結構的被監測的所有參量”簡稱為“被監測量”。給M個被監測量連續編號,該編號在后續步驟中將用于生成向量和矩陣。本方法用用變量j表示這一編號,j=l, 2, 3,…,M。
            [0020]本方法的第一部分:建立結構健康監測系統所需的知識庫和參量的方法。具體如下:
            [0021]1.首先確定“本方法的索結構的溫度測量計算方法”。由于索結構的溫度可能是變化的,例如索結構的不同部位的溫度是隨著日照強度的變化而變化、隨著環境溫度的變化而變化的,索結構的表面與內部的溫度有時可能是隨時間變化的,索結構的表面與內部的溫度可能是不同的,索結構的表面與內部的溫度差是隨時間變化的,這就使得考慮溫度條件時的索結構的力學計算和監測相當復雜,為簡化問題、減少計算量和降低測量成本,更是為了提高計算精度,本方法提出“本方法的索結構的溫度測量計算方法”,具體如下:
            [0022]第一步,查詢或實測得到索結構組成材料及索結構所處環境的隨溫度變化的傳熱學參數,利用索結構的設計圖、竣工圖和索結構的幾何實測數據,利用這些數據和參數建立索結構的傳熱學計算模型。查詢索結構所在地不少于2年的近年來的氣象資料,統計得到這段時間內的陰天數量記為T個陰天,統計得到T個陰天中每一個陰天的0時至次日日出時刻后30分鐘之間的最高氣溫與最低氣溫,日出時刻是指根據地球自轉和公轉規律確定的氣象學上的日出時刻,可以查詢資料或通過常規氣象學計算得到所需的每一日的日出時亥IJ,每一個陰天的0時至次日日出時刻后30分鐘之間的最高氣溫減去最低氣溫稱為該陰天的日氣溫的最大溫差,有T個陰天,就有T個陰天的日氣溫的最大溫差,取T個陰天的日氣溫的最大溫差中的最大值為參考日溫差,參考日溫差記為Al;。查詢索結構所在地和所在海拔區間不少于2年的近年來的氣象資料或實測得到索結構所處環境的溫度隨時間和海拔高度的變化數據和變化規律,計算得到索結構所在地和所在海拔區間不少于2年的近年來的索結構所處環境的溫度關于海拔高度的最大變化率△ Th,為方便敘述取ATh的單位為。C /m。在索結構的表面上取“R個索結構表面點”,后面將通過實測得到這R個索結構表面點的溫度,稱實測得到的溫度數據為“R個索結構表面溫度實測數據”,如果是利用索結構的傳熱學計算模型,通過傳熱計算得到這R個索結構表面點的溫度,就稱計算得到的溫度數據為“R個索結構表面溫度計算數據”。在索結構的表面上取“R個索結構表面點”時,“R個索結構表面點”的數量與分布必須滿足的條件在后面敘述。從索結構所處的最低海拔到最聞海拔之間,在索結構上均布選取不少于二個不同的海拔聞度,在每一個選取的海拔聞度處、在水平面與索結構表面的交線處至少選取兩個點,從選取點處引索結構表面的外法線,所有選取的外法線方向稱為“測量索結構沿壁厚的溫度分布的方向”,測量索結構沿壁厚的溫度分布的方向與“水平面與索結構表面的交線”相交,在選取的測量索結構沿壁厚的溫度分布的方向中必須包括`索結構的向陽面外法線方向和索結構的背陰面外法線方向,沿每一個測量索結構沿壁厚的溫度分布的方向在索結構中均布選取不少于三個點,測量所有被選取點的溫度,測得的溫度稱為“索結構沿厚度的溫度分布數據”,其中沿與同一“水平面與索結構表面的交線”相交的、“測量索結構沿壁厚的溫度分布的方向”測量獲得的“索結構沿厚度的溫度分布數據”,在本方法中稱為“相同海拔高度索結構沿厚度的溫度分布數據”,設選取了 H個不同的海拔高度,在每一個海拔高度處,選取了 B個測量索結構沿壁厚的溫度分布的方向,沿每個測量索結構沿壁厚的溫度分布的方向在索結構中選取了 E個點,其中H和E都不小于3,B不小于2,設HBE為H與B和E的乘積,對應的共有HBE個“測量索結構沿厚度的溫度分布數據的點”,后面將通過實測得到這HBE個“測量索結構沿厚度的溫度分布數據的點”的溫度,稱實測得到的溫度數據為“HBE個索結構沿厚度溫度實測數據”,如果是利用索結構的傳熱學計算模型,通過傳熱計算得到這HBE個測量索結構沿厚度的溫度分布數據的點的溫度,就稱計算得到的溫度數據為“HBE個索結構沿厚度溫度計算數據”;設BE為B和E的乘積,本方法中在每一個選取的海拔高度處共有BE個“相同海拔高度索結構沿厚度的溫度分布數據”。在索結構所在地按照氣象學測量氣溫要求選取一個位置,將在此位置實測得到符合氣象學測量氣溫要求的索結構所在環境的氣溫;在索結構所在地的空曠無遮擋處選取一個位置,該位置應當在全年的每一日都能得到該地所能得到的該日的最充分的日照,在該位置安放一塊碳鋼材質的平板,稱為參考平板,該參考平板的一面向陽,稱為向陽面,參考平板的向陽面是粗糙的和深色的,參考平板的向陽面應當在全年的每一日都能得到一塊平板在該地所能得到的該日的最充分的日照,參考平板的非向陽面覆有保溫材料,將實時監測得到參考平板的向陽面的溫度。本方法中對同一個量實時監測的任何兩次測量之間的時間間隔不得大于30分鐘,測量記錄數據的時刻稱為實際記錄數據時刻。
            [0023]第二步,實時監測得到上述R個索結構表面點的R個索結構表面溫度實測數據,同時實時監測得到前面定義的索結構沿厚度的溫度分布數據,同時實時監測得到符合氣象學測量氣溫要求的索結構所在環境的氣溫數據;通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的索結構所在環境的氣溫實測數據序列,索結構所在環境的氣溫實測數據序列由當日日出時刻到次日日出時刻后30分鐘之間的索結構所在環境的氣溫實測數據按照時間先后順序排列,找到索結構所在環境的氣溫實測數據序列中的最高溫度和最低溫度,用索結構所在環境的氣溫實測數據序列中的最高溫度減去最低溫度得到索結構所在環境的當日日出時刻到次日日出時刻后30分鐘之間的最大溫差,記為AT6max ;由索結構所在環境的氣溫實測數據序列通過常規數學計算得到索結構所在環境的氣溫關于時間的變化率,該變化率也隨著時間變化;通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的參考平板的向陽面的溫度的實測數據序列,參考平板的向陽面的溫度的實測數據序列由當日日出時刻到次日日出時刻后30分鐘之間的參考平板的向陽面的溫度的實測數據按照時間先后順序排列,找到參考平板的向陽面的溫度的實測數據序列中的最高溫度和最低溫度,用參考平板的向陽面的溫度的實測數據序列中的最高溫度減去最低溫度得到參考平板的向陽面的溫度的當日日出時刻到次日日出時刻后30分鐘之間的最大溫差,記為A Tpfflax ;通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的所有R個索結構表面點的索結構表面溫度實測數據序列,有R個索結構表面點就有R個索結構表面溫度實測數據序列,每一個索結構表面溫度實測數據序列由一個索結構表面點的當日日出時刻到次日日出時刻后30分鐘之間的索結構表面溫度實測數據按照時間先后順序排列,找到每一個索結構表面溫度實測數據序列中的最高溫度和最低溫度,用每一個索結構表面溫度實測數據序列中的最高溫度減去最低溫度得到每一個索結構表面點的溫度的當日日出時刻到次日日出時刻后30分鐘之間的最大溫差,有R個索結構表面點就有R個當日日出時刻到次日日出時刻后30分鐘之間的最大溫差數值,其中的最大值記為ATsmax ;由每一索結構表面溫度實測數據序列通過常規數學計算得到每一個索結構表面點的溫度關于時間的變化率,每一個索結構表面點的溫度關于時間的變化率也隨著時間變化。通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的、在同一時刻、HBE個“索結構沿厚度的溫度分布數據”后,計算在每一個選取的海拔高度處共計BE個“相同海拔高度索結構沿厚度的溫度分布數據”中的最高溫度與最低溫度的差值,這個差值的絕對值稱為“相同海拔高度處索結構厚度方向最大溫差”,選取了 H個不同的海拔高度就有H個“相同海拔高度處索結構厚度方向最大溫差”,稱這H個“相同海拔高度處索結構厚度方向最大溫差”中的最大值為“索結構厚度方向最大溫差”,記為ATtmax。[0024]第三步,測量計算獲得索結構穩態溫度數據;首先,確定獲得索結構穩態溫度數據的時刻,與決定獲得索結構穩態溫度數據的時刻相關的條件有六項,第一項條件是獲得索結構穩態溫度數據的時刻介于當日日落時刻到次日日出時刻后30分鐘之間,日落時刻是指根據地球自轉和公轉規律確定的氣象學上的日落時刻,可以查詢資料或通過常規氣象學計算得到所需的每一日的日落時刻;第二項條件的a條件是在當日日出時刻到次日日出時刻后30分鐘之間的這段時間內,參考平板最大溫差△ Tpmax和索結構表面最大溫差ATsmax都不大于5攝氏度;第二項條件的b條件是在當日日出時刻到次日日出時刻后30分鐘之間的這段時間內,在前面測量計算得到的環境最大誤差AT6max不大于參考日溫差Al;,且參考平板最大溫差A Tpmax減去2攝氏度后不大于A Temax,且索結構表面最大溫差A Tsmax不大于ATpmax ;只需滿足第二項的a條件和b條件中的一項就稱為滿足第二項條件;第三項條件是在獲得索結構穩態溫度數據的時刻,索結構所在環境的氣溫關于時間的變化率的絕對值不大于每小時0.1攝氏度;第四項條件是在獲得索結構穩態溫度數據的時刻,R個索結構表面點中的每一個索結構表面點的溫度關于時間的變化率的絕對值不大于每小時0.1攝氏度;第五項條件是在獲得索結構穩態溫度數據的時刻,R個索結構表面點中的每一個索結構表面點的索結構表面溫度實測數據為當日日出時刻到次日日出時刻后30分鐘之間的極小值;第六項條件是在獲得索結構穩態溫度數據的時刻,“索結構厚度方向最大溫差” ATtmax不大于I攝氏度;本方法利用上述六項條件,將下列三種時刻中的任意一種稱為“獲得索結構穩態溫度數據的數學時刻”,第一種時刻是滿足上述“與決定獲得索結構穩態溫度數據的時刻相關的條件”中的第一項至第五項條件的時刻,第二種時刻是僅僅滿足上述“與決定獲得索結構穩態溫度數據的時刻相關的條件”中的第六項條件的時刻,第三種時刻是同時滿足上述“與決定獲得索結構穩態溫度數據的時刻相關的條件”中的第一項至第六項條件的時刻;當獲得索結構穩態溫度數據的數學時刻就是本方法中實際記錄數據時刻中的一個時,獲得索結構穩態溫度數據的時刻就是獲得索結構穩態溫度數據的數學時刻;如果獲得索結構穩態溫度數據的數學時刻不是本方法中實際記錄數據時刻中的任一個時刻,則取本方法最接近于獲得索結構穩態溫度數據的數學時刻的那個實際記錄數據的時刻為獲得索結構穩態溫度數據的時刻;本方法將使用在獲得索結構穩態溫度數據的時刻測量記錄的量進行索結構相關健康監測分析;本方法近似認為獲得索結構穩態溫度數據的時刻的索結構溫度場處于穩態,即此時刻的索結構溫度不隨時間變化,此時刻就是本方法的“獲得索結構穩態溫度數據的時刻”;然后`,根據索結構傳熱特性,利用獲得索結構穩態溫度數據的時刻的“R個索結構表面溫度實測數據”和“HBE個索結構沿厚度溫度實測數據”,利用索結構的傳熱學計算模型,通過常規傳熱計算得到在獲得索結構穩態溫度數據的時刻的索結構的溫度分布,此時索結構的溫度場按穩態進行計算,計算得到的在獲得索結構穩態溫度數據的時刻的索結構的溫度分布數據包括索結構上R個索結構表面點的計算溫度,R個索結構表面點的計算溫度稱為R個索結構穩態表面溫度計算數據,還包括索結構在前面選定的HBE個“測量索結構沿厚度的溫度分布數據的點”的計算溫度,HBE個“測量索結構沿厚度的溫度分布數據的點”的計算溫度稱為“HBE個索結構沿厚度溫度計算數據”,當R個索結構表面溫度實測數據與R個索結構穩態表面溫度計算數據對應相等時,且“HBE個索結構沿厚度溫度實測數據”與“HBE個索結構沿厚度溫度計算數據”對應相等時,計算得到的在獲得索結構穩態溫度數據的時刻的索結構的溫度分布數據在本方法中稱為“索結構穩態溫度數據”,此時的“R個索結構表面溫度實測數據”稱為“R個索結構穩態表面溫度實測數據”,“HBE個索結構沿厚度溫度實測數據”稱為“HBE個索結構沿厚度穩態溫度實測數據”;在索結構的表面上取“R個索結構表面點”時,“R個索結構表面點”的數量與分布必須滿足三個條件,第一個條件是當索結構溫度場處于穩態時,當索結構表面上任意一點的溫度是通過“R個索結構表面點”中與索結構表面上該任意點相鄰的點的實測溫度線性插值得到時,線性插值得到的索結構表面上該任意點的溫度與索結構表面上該任意點的實際溫度的誤差不大于5% ;索結構表面包括支承索表面;第二個條件是“R個索結構表面點”中在同一海拔高度的點的數量不小于4,且“R個索結構表面點”中在同一海拔高度的點沿著索結構表面均布;“R個索結構表面點”沿海拔高度的所有兩兩相鄰索結構表面點的海拔高度之差的絕對值中的最大值Ah不大于0.2°C除以A Th得到的數值,為方便敘述取A Th的單位為。C/m,為方便敘述取Ah的單位為m ;“R個索結構表面點”沿海拔高度的兩兩相鄰索結構表面點的定義是指只考慮海拔高度時,在“R個索結構表面點”中不存在一個索結構表面點,該索結構表面點的海拔高度數值介于兩兩相鄰索結構表面點的海拔高度數值之間;第三個條件是查詢或按氣象學常規計算得到索結構所在地和所在海拔區間的日照規律,再根據索結構的幾何特征及方位數據,在索結構上找到全年受日照時間最充分的那些表面點的位置,“R個索結構表面點”中至少有一個索結構表面點是索結構上全年受日照時間最充分的那些表面點中的一個點。
            [0025]2.建立索結構的初始力學計算基準模型A。(例如有限元基準模型)和當前初始力學計算基準模型(例如有限元基準模型)的方法,建立與A。對應的被監測量初始數值向量C。的方法,建立與Attj對應的被監測量當前初始數值向量Cttj的方法。在本方法中K、c0、Attj和Cttj是不斷更新的。建立和更新A^C。、#。和Cttj的方法如下。被監測量初始數值向量C0的編號規則與M個被監測量的編號規則相同。
            [0026]建立初始力學計算基準模型A。時,在索結構竣工之時,或者在建立結構健康監測系統前,按照“本方法的索結構的溫度測量計算方法”測量計算得到“索結構穩態溫度數據”(可以用常規溫度測量方法測量,例如使用熱電阻測量),此時的“索結構穩態溫度數據”用向量T。表示,稱為初始索結構穩態溫度數據向量T。。在實測得到T。的同時,也就是在獲得索結構穩態溫度數據的時刻的同一時刻,使用常規方法直接測量計算得到索結構的所有被監測量的初始數。使用常規方法(查資料或實測)得到索結構所使用的各種材料的隨溫度變化的物理參數(例如熱膨脹系數)和力學性能參數(例如彈性模量、泊松比);在實測計算得到初始索結構穩態溫度數據向量T。的同時,也就是在獲得索結構穩態溫度數據的時刻的同一時刻,使用常規方法實測計算得到索結構的實測計算數據。索結構的實測計算數據首先是包括支承索的無損檢測數據在內的能夠表達索的健康狀態的數據,索結構的實測計算數據還是包括索結構支座初始廣義位移測量數據、索結構初始幾何數據、索力數據、拉桿拉力數據、初始索結構支座廣義坐標數據、索結構模態數據、結構應變數據、結構角度測量數據、結構空間坐標測量數據和載荷數據在內的實測數據。索結構的初始幾何數據可以是所有索的端點的空間坐標數據加上結構上一系列的點的空間坐標數據,目的在于根據這些坐標數據確定索結構的幾何特征。索結構支座初始廣義位移測量數據指在建立初始力學計算基準模型A。時,索結構支座相對于索結構設計狀態下的支座所發生的廣義位移。對斜拉橋而言,初始幾何數據可以是所有索的端點的空間坐標數據加上橋梁兩端上若干點的空間坐標數據,這就是所謂的橋型數據。“所有可能發生變化的集中載荷”的變化量在建立初始力學計算基準模型A。時全部為O,也就是說后面識別出的“所有可能發生變化的集中載荷”的變化量是相對于建立初始力學計算基準模型A。時結構所承受的對應集中載荷的變化量。利用支承索的無損檢測數據等能夠表達支承索的健康狀態的數據、索結構支座初始廣義位移測量數據以及“所有可能發生變化的集中載荷”的變化量數據建立被評估對象初始損傷向量d。(如式(I)所示),用d。表示索結構(用初始力學計算基準模型A。表示)的被評估對象的初始健康狀態。如果沒有支承索的無損檢測數據及其他能夠表達支承索的健康狀態的數據時,或者可以認為結構初始狀態為無損傷無松弛狀態時,向量d。的中與支承索相關的各元素數值取O。向量d。中與集中載荷的變化量相關的各元素數值取O。如果沒有索結構支座初始廣義位移測量數據或者可以認為索結構支座初始廣義位移為O時,向量d。的中與索結構支座廣義位移相關的各元素數值取O。利用索結構的設計圖、竣工圖和初始索結構的實測數據、支承索的無損檢測數據、索結構支座初始廣義位移測量數據和索結構所使用的各種材料的隨溫度變化的物理和力學性能參數和初始索結構穩態溫度數據向量T。,利用力學方法(例如有限元法)計入“索結構穩態溫度數據”建立初始力學計算基準模型A。。
            [0027]
            【權利要求】
            1.混合監測受損索集中載荷廣義位移識別方法,其特征在于所述方法包括: a.為敘述方便起見,本方法統一稱被評估的支座廣義位移分量、支承索和集中載荷為被評估對象,設被評估的支座廣義位移分量的數量、支承索的數量和集中載荷的數量之和為N,即被評估對象的數量為N ;確定被評估對象的編號規則,按此規則將索結構中所有的被評估對象編號,該編號在后續步驟中將用于生成向量和矩陣;本方法用變量k表示這一編號,k=l,2,3,…,N ;確定混合監測時指定的將被監測索力的支承索,設索系統中共有Q根支承索,索結構的被監測的索力數據由索結構上M1個指定支承索的M1個索力數據來描述,索結構索力的變化就是所有指定支承索的索力的變化;每次共有M1個索力測量值或計算值來表征索結構的索力信息^是一個不小于O不大于Q的整數;確定混合監測時指定的將被監測應變的被測量點,索結構的被監測的應變數據可由索結構上K2個指定點的、及每個指定點的L2個指定方向的應變來描述,索結構應變數據的變化就是K2個指定點的所有被測應變的變化;每次共有M2個應變測量值或計算值來表征索結構應變,M2為K2和L2之積;M2是不小于O的整數;確定混合監測時指定的將被監測角度的被測量點,索結構的被監測的角度數據由索結構上K3個指定點的、過每個指定點的L3個指定直線的、每個指定直線的H3個角度坐標分量來描述,索結構角度的變化就是所有指定點的、所有指定直線的、所有指定的角度坐標分量的變化;每次共有M3個角度坐標分量測量值或計算值來表征索結構的角度信息,M3為K3、L3和H3之積;M3是一個不小于O的整數;確定混合監測時指定的將被監測的形狀數據,索結構的被監測的形狀數據由索結構上K4個指定點的、及每個指定點的L4個指定方向的空間坐標來描述,索結構形狀數據的變化就是K4個指定點的所有坐標分量的變化;每次共有M4個坐標測量值或計算值來表征索結構形狀,M4為K4和L4之積;M4是一個不小于O的整數;綜合上述混合監測的被監測量,整個索結構共有M個被監測量,M為%、M2'M3和M4之和,定義參量K,K為M1、K2、K3和K4之和,K和M不得小于被評估對象的數量N ;為方便起見,在本方法中將本步所列出的M個被監測量簡稱為“被監測量”;本方法中對同一個量實時監測的任何兩次測量之間的時間間隔不得大于30分鐘,測量記錄數據的時刻稱為實際記錄數據時刻; b.本方法定義“本方法的索結構的溫度測量計算方法”按步驟bl至b3進行; bl:查詢或實測得到索結構組成材料及索結構所處環境的隨溫度變化的傳熱學參數,利用索結構的設計圖、竣工圖和索結構的幾何實測數據,利用這些數據和參數建立索結構的傳熱學計算模型;查詢索結構所在地不少于2年的近年來的氣象資料,統計得到這段時間內的陰天數量記為T個陰天,在本方法中將白天不能見到太陽的一整日稱為陰天,統計得到T個陰天中每一個陰天的0時至次日日出時刻后30分鐘之間的最高氣溫與最低氣溫,日出時刻是指根據地球自轉和公轉規律確定的氣象學上的日出時刻,不表示當天一定可以看見太陽,可以查詢資料或通過常規氣象學計算得到所需的每一日的日出時刻,每一個陰天的0時至次日日出時刻后30分鐘之間的最高氣溫減去最低氣溫稱為該陰天的日氣溫的最大溫差,有T個陰天,就有T個陰天的日氣溫的最大溫差,取T個陰天的日氣溫的最大溫差中的最大值為參考日溫差,參考日溫差記為△!;;查詢索結構所在地和所在海拔區間不少于2年的近年來的氣象資料或實測得到索結構所處環境的溫度隨時間和海拔高度的變化數據和變化規律,計算得到索結構所在地和所在海拔區間不少于2年的近年來的索結構所處環境的溫度關于海拔高度的最大變化率△ Th,為方便敘述取ATh的單位為。C/m;在索結構的表面上取“R個索結構表面點”,取“R個索結構表面點”的具體原則在步驟b3中敘述,后面將通過實測得到這R個索結構表面點的溫度,稱實測得到的溫度數據為“R個索結構表面溫度實測數據”,如果是利用索結構的傳熱學計算模型,通過傳熱計算得到這R個索結構表面點的溫度,就稱計算得到的溫度數據為“R個索結構表面溫度計算數據”;從索結構所處的最低海拔到最高海拔之間,在索結構上均布選取不少于三個不同的海拔高度,在每一個選取的海拔高度處、在水平面與索結構表面的交線處至少選取兩個點,從選取點處引索結構表面的外法線,所有選取的外法線方向稱為“測量索結構沿壁厚的溫度分布的方向”,測量索結構沿壁厚的溫度分布的方向與“水平面與索結構表面的交線”相交,在選取的測量索結構沿壁厚的溫度分布的方向中必須包括索結構的向陽面外法線方向和索結構的背陰面外法線方向,沿每一個測量索結構沿壁厚的溫度分布的方向在索結構中均布選取不少于三個點,測量所有被選取點的溫度,測得的溫度稱為“索結構沿厚度的溫度分布數據”,其中沿與同一“水平面與索結構表面的交線”相交的、“測量索結構沿壁厚的溫度分布的方向”測量獲得的“索結構沿厚度的溫度分布數據”,在本方法中稱為“相同海拔高度索結構沿厚度的溫度分布數據”,設選取了 H個不同的海拔高度,在每一個海拔高度處,選取了 B個測量索結構沿壁厚的溫度分布的方向,沿每個測量索結構沿壁厚的溫度分布的方向在索結構中選取了 E個點,其中H和E都不小于3,B不小于2,設HBE為H與B和E的乘積,對應的共有HBE個“測量索結構沿厚度的溫度分布數據的點”,后面將通過實測得到這HBE個“測量索結構沿厚度的溫度分布數據的點”的溫度,稱實測得到的溫度數據為“HBE個索結構沿厚度溫度實測數據”,如果是利用索結構的傳熱學計算模型,通過傳熱計算得到這HBE個測量索結構沿厚度的溫度分布數據的點的溫度,就稱計算得到的溫度數據為“HBE個索結構沿厚度溫度計算數據”;設BE為B和E的乘積,本方法中在每一個選取的海拔高度處共有BE個“相同海拔高度索結構沿厚度的溫度分布數據”;在索結構所在地按照氣象學測量氣溫要求選取一個位置,將在此位置實測得到符合氣象學測量氣溫要求的索結構所在環境的氣溫;在索結構所在地的空曠無遮擋處選取一個位置,該位置應當在全年的每一日都能得到該地所能得到的該日的最充分的日照,在該位置安放一塊碳鋼材質的平板,稱為參考平板,參考平板與地面不可接觸,參考平板離地面距離不小于1.5米,該參考平板的一面向陽,稱為向陽面,參考平板的向陽面是粗糙的和深色的,參考平板的向陽面應當在全年的每一日都能得到一塊平板在該地所能得到的該日的最充分的日照,參考平板的非向陽面覆有保溫材料,將實時監測得到參考平板的向陽面的溫度; b2:實時監測得到上述R個索結構表面點的R個索結構表面溫度實測數據,同時實時監測得到前面定義的索結構沿厚度的溫度分布數據,同時實時監測得到符合氣象學測量氣溫要求的索結構所在環境的氣溫數據;通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的索結構所在環境的氣溫實測數據序列,索結構所在環境的氣溫實測數據序列由當日日出時刻到次日日出時刻后30分鐘之間的索結構所在環境的氣溫實測數據按照時間先后順序排列,找到索結構所在環境的氣溫實測數據序列中的最高溫度和最低溫度,用索結構所在環境的氣溫實測數據序列中的最高溫度減去最低溫度得到索結構所在環境的當日日出時刻到次日日出時刻后30分鐘之間的最大溫差,稱為環境最大溫差,記為AT6max;由索結構所在環境的氣溫實測數據序列通過常規數學計算得到索結構所在環境的氣溫關于時間的變化率,該變化率也隨著時間變化;通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的參考平板的向陽面的溫度的實測數據序列,參考平板的向陽面的溫度的實測數據序列由當日日出時刻到次日日出時刻后30分鐘之間的參考平板的向陽面的溫度的實測數據按照時間先后順序排列,找到參考平板的向陽面的溫度的實測數據序列中的最高溫度和最低溫度,用參考平板的向陽面的溫度的實測數據序列中的最高溫度減去最低溫度得到參考平板的向陽面的溫度的當日日出時刻到次日日出時刻后30分鐘之間的最大溫差,稱為參考平板最大溫差,記為ATpmax ;通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的所有R個索結構表面點的索結構表面溫度實測數據序列,有R個索結構表面點就有R個索結構表面溫度實測數據序列,每一個索結構表面溫度實測數據序列由一個索結構表面點的當日日出時刻到次日日出時刻后30分鐘之間的索結構表面溫度實測數據按照時間先后順序排列,找到每一個索結構表面溫度實測數據序列中的最高溫度和最低溫度,用每一個索結構表面溫度實測數據序列中的最高溫度減去最低溫度得到每一個索結構表面點的溫度的當日日出時刻到次日日出時刻后30分鐘之間的最大溫差,有R個索結構表面點就有R個當日日出時刻到次日日出時刻后30分鐘之間的最大溫差數值,其中的最大值稱為索結構表面最大溫差,記為ATsmax ;由每一索結構表面溫度實測數據序列通過常規數學計算得到每一個索結構表面點的溫度關于時間的變化率,每一個索結構表面點的溫度關于時間的變化率也隨著時間變化;通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的、在同一時刻、HBE個“索結構沿厚度的溫度分布數據”后,計算在每一個選取的海拔高度處共計BE個“相同海拔高度索結構沿厚度的溫度分布數據”中的最高溫度與最低溫度的差值,這個差值的絕對值稱為“相同海拔高度處索結構厚度方向最大溫差”,選取了 H個不同的海拔高度就有H個“相同海拔高度處索結構厚度方向最大溫差”,稱這H個“相同海拔高度處索結構厚度方向最大溫差”中的最大值為“索結構厚度方向最大溫差”,記為A Ttmax ; b3:測量計算獲得索結構穩態溫度數據;首先,確定獲得索結構穩態溫度數據的時刻,與決定獲得索結構穩態溫度數據的時刻相關的條件有六項,第一項條件是獲得索結構穩態溫度數據的時刻介于當日日落時刻到次日日出時刻后30分鐘之間,日落時刻是指根據地球自轉和公轉規律確定的氣 象學上的日落時刻,可以查詢資料或通過常規氣象學計算得到所需的每一日的日落時刻;第二項條件的a條件是在當日日出時刻到次日日出時刻后30分鐘之間的這段時間內,參考平板最大溫差ATpmaj^P索結構表面最大溫差ATsmax都不大于5攝氏度;第二項條件的b條件是在當日日出時刻到次日日出時刻后30分鐘之間的這段時間內,在前面測量計算得到的環境最大誤差ATraiax不大于參考日溫差A I;,且參考平板最大溫差ATpmax減去2攝氏度后不大于A T6max,且索結構表面最大溫差ATsmax不大于ATpmax;只需滿足第二項的a條件和b條件中的一項就稱為滿足第二項條件;第三項條件是在獲得索結構穩態溫度數據的時刻,索結構所在環境的氣溫關于時間的變化率的絕對值不大于每小時0.1攝氏度;第四項條件是在獲得索結構穩態溫度數據的時刻,R個索結構表面點中的每一個索結構表面點的溫度關于時間的變化率的絕對值不大于每小時0.1攝氏度;第五項條件是在獲得索結構穩態溫度數據的時刻,R個索結構表面點中的每一個索結構表面點的索結構表面溫度實測數據為當日日出時刻到次日日出時刻后30分鐘之間的極小值;第六項條件是在獲得索結構穩態溫度數據的時刻,“索結構厚度方向最大溫差” A Ttmax不大于I攝氏度;本方法利用上述六項條件,將下列三種時刻中的任意一種稱為“獲得索結構穩態溫度數據的數學時刻”,第一種時刻是滿足上述“與決定獲得索結構穩態溫度數據的時刻相關的條件”中的第一項至第五項條件的時刻,第二種時刻是僅僅滿足上述“與決定獲得索結構穩態溫度數據的時刻相關的條件”中的第六項條件的時刻,第三種時刻是同時滿足上述“與決定獲得索結構穩態溫度數據的時刻相關的條件”中的第一項至第六項條件的時刻;當獲得索結構穩態溫度數據的數學時刻就是本方法中實際記錄數據時刻中的一個時,獲得索結構穩態溫度數據的時刻就是獲得索結構穩態溫度數據的數學時刻;如果獲得索結構穩態溫度數據的數學時刻不是本方法中實際記錄數據時刻中的任一個時刻,則取本方法最接近于獲得索結構穩態溫度數據的數學時刻的那個實際記錄數據的時刻為獲得索結構穩態溫度數據的時刻;本方法將使用在獲得索結構穩態溫度數據的時刻測量記錄的量進行索結構相關健康監測分析;本方法近似認為獲得索結構穩態溫度數據的時刻的索結構溫度場處于穩態,即此時刻的索結構溫度不隨時間變化,此時刻就是本方法的“獲得索結構穩態溫度數據的時刻”;然后,根據索結構傳熱特性,利用獲得索結構穩態溫度數據的時刻的“R個索結構表面溫度實測數據”和“HBE個索結構沿厚度溫度實測數據”,利用索結構的傳熱學計算模型,通過常規傳熱計算得到在獲得索結構穩態溫度數據的時刻的索結構的溫度分布,此時索結構的溫度場按穩態進行計算,計算得到的在獲得索結構穩態溫度數據的時刻的索結構的溫度分布數據包括索結構上R個索結構表面點的計算溫度,R個索結構表面點的計算溫度稱為R個索結構穩態表面溫度計算數據,還包括索結構在前面選定的HBE個“測量索結構沿厚度的溫度分布數據的點”的計算溫度,HBE個“測量索結構沿厚度的溫度分布數據的點”的計算溫度稱為“HBE個索結構沿厚度溫度計算數據”,當R個索結構表面溫度實測數據與R個索結構穩態表面溫度計算數據對應相等時,且“HBE個索結構沿厚度溫度實測數據”與“HBE個索結構沿厚度溫度計算數據”對應相等時,計算得到的在獲得索結構穩態溫度數據的時刻的索結構的溫度分布數據在本方法中稱為“索結構穩態溫度數據”,此時的“R個索結構表面溫度實測數據”稱為“R個索結構穩態表面溫度實測數據”,“HBE個索結構沿厚度溫度實測數據”稱為“HBE個索結構沿厚度穩態溫度實測數據”;在索結構的表面上取“R個索結構表面點”時,“R個索結構表面點”的數量與分布必須滿足三個條件,第一個條件是當索結構溫度場處于穩態時,當索結構表面上任意一點的溫度是通過“R個索結構表面點”中與索結構表面上該任意點相鄰的點的實測溫度線性插值得到時,線性插值得到的索結構表面上該任意點的溫度與索結構表面上該任意點的實際溫度的誤差不大于5% ;索結構表面包括支承索表面;第二個條件是“R個索結構表面點”中在同一海拔高度的點的數量不小于4,且“R個索結構表面點”中在同一海拔高度的點沿著索結構表面均布;“R個索結構表面點”沿海拔高度的所有兩兩相鄰索結構表面點的海拔高度之差的絕對值中的最大值Δ h不大于0.2°C除以Δ Th得到的數值,為方便敘述取Δ Th的單位為。C/m,為方便敘述取Δh的單位為m ;“R個索結構表面點”沿海拔高度的兩兩相鄰索結構表面點的定義是指只考慮海拔高度時,在“R個索結構表面點”中不存在一個索結構表面點,該索結構表面點的海拔高度數值介于兩兩相鄰索結構表面點的海拔高度數值之間;第三個條件是查詢或按氣象學常規計算得到索結構所在地和所在海拔區間的日照規律,再根據索結構的幾何特征及方位數據,在索結構上找到全年受日照時間最充分的那些表面點的位置,“R個索結構表面點”中至少有一個索結構表面點是索結構上全年受日照時間最充分的那些表面點中的一個點;c.按照“本方法的索結構的溫度測量計算方法”直接測量計算得到初始狀態下的索結構穩態溫度數據,初始狀態下的索結構穩態溫度數據稱為初始索結構穩態溫度數據,記為“初始索結構穩態溫度數據向量T。” ;實測或查資料得到索結構所使用的各種材料的隨溫度變化的物理和力學性能參數;在實測得到T。的同時,也就是在獲得初始索結構穩態溫度數據向量T。的時刻的同一時刻,直接測量計算得到初始索結構的實測數據,初始索結構的實測數據是包括索結構集中載荷測量數據、索結構分布載荷測量數據、索結構體積載荷測量數據、索結構支座初始廣義位移測量數據、所有被監測量的初始數值、所有支承索的初始索力數據、 初始索結構模態數據、初始索結構應變數據、初始索結構幾何數據、初始索結構支座廣義坐標數據、初始索結構角度數據、初始索結構空間坐標數據在內的實測數據,在得到初始索結構的實測數據的同時,測量計算得到包括支承索的無損檢測數據在內的能夠表達支承索的健康狀態的數據,此時的能夠表達支承索的健康狀態的數據稱為支承索初始健康狀態數據;所有被監測量的初始數值組成被監測量初始數值向量C。,被監測量初始數值向量C。的編號規則與M個被監測量的編號規則相同;利用支承索初始健康狀態數據、索結構支座初始廣義位移測量數據和索結構集中載荷測量數據建立被評估對象初始損傷向量d。,向量d。表示用初始力學計算基準模型A。表示的索結構的被評估對象的初始健康狀態;被評估對象初始損傷向量d。的元素個數等于N,d。的元素與被評估對象是一一對應關系,向量d。的元素的編號規則與被評估對象的編號規則相同;如果d。的某一個元素對應的被評估對象是索系統中的一根支承索,那么d。的該元素的數值代表對應支承索的初始損傷程度,若該元素的數值為O,表示該元素所對應的支承索是完好的,沒有損傷的,若其數值為100%,則表示該元素所對應的支承索已經完全喪失承載能力,若其數值介于0和100%之間,則表示該支承索喪失了相應比例的承載能力;如果d。的某一個元素對應的被評估對象是某一個支座的某一個廣義位移分量,那么d。的該元素的數值代表這個支座的該廣義位移分量的初始數值;如果d。的某一個元素對應的被評估對象是某一個集中載荷,本方法中取d。的該元素數值為0,代表這個集中載荷的變化的初始數值為0 ;如果沒有索結構支座初始廣義位移測量數據或者可以認為索結構支座初始廣義位移為0時,向量d。中與索結構支座廣義位移相關的各元素數值取0 ;如果沒有支承索的無損檢測數據及其他能夠表達支承索的健康狀態的數據時,或者可以認為結構初始狀態為無損傷無松弛狀態時,向量d。中與支承索相關的各元素數值取0 ;初始索結構支座廣義坐標數據指索結構設計狀態下的支座廣義坐標數據,索結構支座初始廣義位移測量數據指在建立初始力學計算基準模型A。時,索結構支座相對于索結構設計狀態下的支座所發生的廣義位移; d.根據索結構的設計圖、竣工圖和初始索結構的實測數據、支承索初始健康狀態數據、索結構支座初始廣義位移測量數據、索結構集中載荷測量數據、索結構分布載荷測量數據、索結構體積載荷測量數據、索結構所使用的各種材料的隨溫度變化的物理和力學性能參數、初始索結構穩態溫度數據向量T。和前面步驟得到的所有的索結構數據,建立計入“索結構穩態溫度數據”的索結構的初始力學計算基準模型A。,基于A。計算得到的索結構計算數據必須非常接近其實測數據,其間的差異不得大于5% ;對應于A。的“索結構穩態溫度數據”就是“初始索結構穩態溫度數據向量T。”;對應于A。的被評估對象健康狀態用被評估對象初始損傷向量d。表示;對應于A。的所有被監測量的初始數值用被監測量初始數值向量C。表示;第一次建立計入“索結構穩態溫度數據”的索結構的當前初始力學計算基準模型#。、被監測量當前初始數值向量Cttj和“當前初始索結構穩態溫度數據向量TtJ ;第一次建立索結構的當前初始力學計算基準模型Attj和被監測量當前初始數值向量Cttj時,索結構的當前初始力學計算基準模型Attj就等于索結構的初始力學計算基準模型A。,被監測量當前初始數值向量Cttj就等于被監測量初始數值向量C。'At0對應的“索結構穩態溫度數據”稱為“當前初始索結構穩態溫度數據”,記為“當前初始索結構穩態溫度數據向量Tt/,第一次建立索結構的當前初始力學計算基準模型Attj時,Tt0就等于T。^t0的被評估對象的初始健康狀態與A0的被評估對象的健康狀態相同,也用被評估對象初始損傷向量d。表示,在后面的循環過程中Attj的被評估對象的初始健康狀態始終用被評估對象初始損傷向量d。表示;T。和d。是A0的參數,由A。的力學計算結果得到的所有被監測量的初始數值與C。表示的所有被監測量的初始數值相同,因此也可以說C。由A。的力學計算結果組成Jt0和d。是Attj的參數,Cttj由At0的力學計算結果組成; e.從這里進入由第e步到第m步的循環;在結構服役過程中,不斷按照“本方法的索結構的溫度測量計算方法”不斷實測計算獲得“索結構穩態溫度數據”的當前數據,“索結構穩態溫度數據”的當前數據稱為“當前索結構穩態溫度數據”,記為“當前索結構穩態溫度數據向量Tt'向量Tt的定義方式與向量T。的定義方式相同; f.根據當前索結構穩態溫度數據向量Tt,按照步驟fl至f3更新當前初始力學計算基準模型A1。、被監測量當前初始數值向量Cttj和當前初始索結構穩態溫度數據向量Tttj ; fl.比較Tt與Tt0,如果Tt等于浐。,則A1。、Ct0和Tttj保持不變;否則需要按下列步驟對AVUt0和Tt0進行更新; f2.計算Tt與T。的差,Tt與T。的差就是當前索結構穩態溫度數據關于初始索結構穩態溫度數據的變化,Tt與T。的差用穩態溫度變化向量S表示,S等于Tt減去T。,S表示索結構穩態溫度數據的變化; f3.對A。中的索結構施加溫度變化,施加的溫度變化的數值就取自穩態溫度變化向量S,對A。中的索結構施加的溫度變化后得到更新的當前初始力學計算基準模型#。,更新Y。的同時,I"。所有元素數值也用Tt的所有元素數值對應代替,即更新了 I"。,這樣就得到了正確地對應于Attj的Tttj ;更新Cttj的方法是:當更新Attj后,通過力學計算得到Attj中所有被監測量的、當前的具體數值,這些具體數值組成Ct0 -X0的支承索的初始健康狀態始終用被評估對象初始損傷向量d。表示; g.在當前初始力學計算基準模型Attj的基礎上按照步驟gl至g4進行若干次力學計算,通過計算獲得索結構單位損傷被監測量數值變化矩陣A C和被評估對象單位變化向量Du ; gl.索結構單位損傷被監測量數值變化矩陣AC是不斷更新的,即在更新當前初始力學計算基準模型#。、被監測量當前初始數值向量Cttj和當前初始索結構穩態溫度數據向量Tt0之后,必須接著更新索結構單位損傷被監測量數值變化矩陣△ C和被評估對象單位變化向量Du ; g2.在索結構的當前初始力學計算基準模型Attj的基礎上進行若干次力學計算,計算次數數值上等于所有被評估對象的數量N,有N個評估對象就有N次計算;依據被評估對象的編號規則,依次進行計算;每一次計算假設只有一個被評估對象在原有損傷或廣義位移或集中載荷的基礎上再增加單位損傷或單位廣義位移或集中載荷單位變化,具體的,如果該被評估對象是索系統中的一根支承索,那么就假設該支承索在向量d。表示的該支承索已有損傷的基礎上再增加單位損傷,如果該被評估對象是一個支座的一個方向的廣義位移分量,就假設該支座在該位移方向再增加單位廣義位移,如果該被評估對象是一個集中載荷,就假設該集中載荷在向量d。表示的該集中載荷已有變化量的基礎上再增加集中載荷單位變化,用Duk記錄這一增加的單位損傷或單位廣義位移或集中載荷單位變化,其中k表示增加單位損傷或單位廣義位移或集中載荷單位變化的被評估對象的編號,Duk是被評估對象單位變化向量Du的一個元素,被評估對象單位變化向量Du的元素的編號規則與向量d。的元素的編號規則相同;每一次計算中增加單位損傷或單位廣義位移或集中載荷單位變化的被評估對象不同于其它次計算中增加單位損傷或單位廣義位移或集中載荷單位變化的被評估對象,每一次計算都利用力學方法計算索結構的所有被監測量的當前計算值,每一次計算得到的所有被監測量的當前計算值組成一個被監測量計算當前向量,被監測量計算當前向量的元素編號規則與被監測量初始數值向量C。的元素編號規則相同; g3.每一次計算得到的被監測量計算當前向量減去被監測量當前初始數值向量Cttj得到一個向量,再將該向量的每一個元素都除以該次計算所假設的單位損傷或單位廣義位移或集中載荷單位變化數值,得到一個被監測量單位變化向量,有N個被評估對象就有N個被監測量單位變化向量; g4.由這N個被監測量單位變化向量按照N個被評估對象的編號規則,依次組成有N列的索結構單位損傷被監測量數值變化矩陣AC ;索結構單位損傷被監測量數值變化矩陣AC的每一列對應于一個被監測量單位變化向量;索結構單位損傷被監測量數值變化矩陣AC的每一行對應于同一個被監測量在不同被評估對象增加單位損傷或單位廣義位移或集中載荷單位變化時的不同的單位變化幅度;索結構單位損傷被監測量數值變化矩陣AC的列的編號規則與向量d。的元素的編號規則相同,索結構單位損傷被監測量數值變化矩陣AC的行的編號規則與M個被監測量的編號規則相同; h.在實測得到當前索結構穩態溫度數據向量Tt的同時,實測得到在獲得當前索結構穩態溫度數據向量Tt的時刻的 同一時刻的索結構的所有被監測量的當前實測數值,組成被監測量當前數值向量C ;被監測量當前數值向量C和被監測量當前初始數值向量Cttj與被監測量初始數值向量C。的定義方式相同,三個向量的相同編號的元素表示同一被監測量在不同時刻的具體數值;i.定義被評估對象當前名義損傷向量d,被評估對象當前名義損傷向量d的元素個數等于被評估對象的數量,被評估對象當前名義損傷向量d的元素和被評估對象之間是一一對應關系,被評估對象當前名義損傷向量d的元素數值代表對應被評估對象的名義損傷程度或名義廣義位移或名義集中載荷變化量;向量d的元素的編號規則與向量d。的元素的編號規則相同; j.依據被監測量當前數值向量C同被監測量當前初始數值向量C1。、索結構單位損傷被監測量數值變化矩陣AC和待求的被評估對象當前名義損傷向量d間存在的近似線性關系,該近似線性關系可表達為式1,式I中除d外的其它量均為已知,求解式I就可以算出被評估對象當前名義損傷向量d ; C=C.+ AC ?d 式 I k.定義被評估對象當前實際損傷向量da,被評估對象當前實際損傷向量da的元素個數等于被評估對象的數量,被評估對象當前實際損傷向量da的元素和被評估對象之間是一一對應關系,被評估對象當前實際損傷向量da的元素數值代表對應被評估對象的實際損傷程度或實際廣義位移或實際集中載荷變化量;向量da的元素的編號規則與向量d。的元素的編號規則相同;l.利用式2表達的被評估對象當前實際損傷向量da的第k個元素d\同被評估對象初始損傷向量d。的第k個元素(I。,和被評估對象當前名義損傷向量d的第k個元素dk間的關系,計算得到被評估對象當前實際損傷向量da的所有元素;

            【文檔編號】G01K13/00GK103616219SQ201310662061
            【公開日】2014年3月5日 申請日期:2013年12月9日 優先權日:2013年12月9日
            【發明者】韓玉林, 韓佳邑 申請人:東南大學
            網友詢問留言 已有0條留言
            • 還沒有人留言評論。精彩留言會獲得點贊!
            1
            婷婷六月激情在线综合激情,亚洲国产大片,久久中文字幕综合婷婷,精品久久久久久中文字幕,亚洲一区二区三区高清不卡,99国产精品热久久久久久夜夜嗨 ,欧美日韩亚洲综合在线一区二区,99国产精品电影,伊人精品线视天天综合,精品伊人久久久大香线蕉欧美
            亚洲精品1区 国产成人一级 91精品国产欧美一区二区 亚洲精品乱码久久久久久下载 国产精品久久久久久久伊一 九色国产 国产精品九九视频 伊人久久成人爱综合网 欧美日韩亚洲区久久综合 欧美日本一道免费一区三区 夜夜爽一区二区三区精品 欧美日韩高清一区二区三区 国产成人av在线 国产精品对白交换绿帽视频 国产视频亚洲 国产在线欧美精品 国产精品综合网 国产日韩精品欧美一区色 国产日韩精品欧美一区喷 欧美日韩在线观看区一二 国产区精品 欧美视频日韩视频 中文字幕天天躁日日躁狠狠躁97 视频一二三区 欧美高清在线精品一区二区不卡 国产精品揄拍一区二区久久 99久久综合狠狠综合久久aⅴ 亚洲乱码视频在线观看 日韩在线第二页 亚洲精品无码专区在线播放 成人亚洲网站www在线观看 欧美三级一区二区 99久久精品免费看国产高清 91麻豆国产在线观看 最新日韩欧美不卡一二三区 成人在线观看不卡 日韩国产在线 在线亚洲精品 亚洲午夜久久久久中文字幕 国产精品成人久久久久久久 精品国产一区二区在线观看 欧美精品国产一区二区三区 中文在线播放 亚洲第一页在线视频 国产午夜精品福利久久 九色国产 精品国产九九 国产永久视频 久久精品人人做人人综合试看 国产一区二区三区免费观看 亚洲精品国产电影 9999热视频 国产精品资源在线 麻豆久久婷婷国产综合五月 国产精品免费一级在线观看 亚洲国产一区二区三区青草影视 中文在线播放 国产成人综合在线 国产在线观看色 国产亚洲三级 国产片一区二区三区 久久99精品久久久久久牛牛影视 亚洲欧美日韩国产 四虎永久免费网站 国产一毛片 国产精品视频在 九九热在线精品 99精品福利视频 色婷婷色99国产综合精品 97成人精品视频在线播放 精品久久久久久中文字幕 亚洲欧美一区二区三区孕妇 亚洲欧美成人网 日韩高清在线二区 国产尤物在线观看 在线不卡一区二区 91网站在线看 韩国精品福利一区二区 欧美日韩国产成人精品 99热精品久久 国产精品免费视频一区 高清视频一区 精品九九久久 欧美日韩在线观看免费 91欧美激情一区二区三区成人 99福利视频 亚洲国产精品91 久热国产在线 精品久久久久久中文字幕女 国产精品久久久久久久久99热 成人自拍视频网 国产精品视频久久久久久 久久影院国产 国产玖玖在线观看 99精品在线免费 亚洲欧美一区二区三区导航 久久久久久久综合 国产欧美日韩精品高清二区综合区 国产精品视频自拍 亚洲一级片免费 久久久久久九九 国产欧美自拍视频 视频一区二区在线观看 欧美日韩一区二区三区久久 中文在线亚洲 伊人热人久久中文字幕 日韩欧美亚洲国产一区二区三区 欧美亚洲国产成人高清在线 欧美日韩国产码高清综合人成 国产性大片免费播放网站 亚洲午夜综合网 91精品久久一区二区三区 国产无套在线播放 国产精品视频网站 国产成人亚洲精品老王 91在线网站 国产视频97 欧美黑人欧美精品刺激 国产一区二区三区免费在线视频 久久久国产精品免费看 99re6久精品国产首页 久久精品91 国产成人一级 国产成人精品曰本亚洲 日本福利在线观看 伊人成综合网 久久综合一本 国产综合久久久久久 久久精品成人免费看 久久福利 91精品国产91久久久久久麻豆 亚洲精品成人在线 亚洲伊人久久精品 欧美日本二区 国产永久视频 国产一区二 一区二区福利 国产一毛片 亚洲精品1区 毛片一区二区三区 伊人久久大香线蕉综合影 国产欧美在线观看一区 亚洲国产欧洲综合997久久 国产一区二区免费视频 国产91精品对白露脸全集观看 久久亚洲国产伦理 欧美成人伊人久久综合网 亚洲性久久久影院 久久99国产精一区二区三区! 91精品国产欧美一区二区 欧美日韩亚洲区久久综合 日韩精品一二三区 久久久夜色精品国产噜噜 国产在线精品福利91香蕉 久久久久久久亚洲精品 97se色综合一区二区二区 91国语精品自产拍在线观看性色 91久久国产综合精品女同我 日韩中文字幕a 国产成人亚洲日本精品 久久国产精品-国产精品 久久国产经典视频 久久国产精品伦理 亚洲第一页在线视频 国产精品久久久久三级 日韩毛片网 久久免费高清视频 麻豆国产在线观看一区二区 91麻豆国产福利在线观看 国产成人精品男人的天堂538 一区二区三区中文字幕 免费在线视频一区 欧美日韩国产成人精品 国产综合网站 国产资源免费观看 亚洲精品亚洲人成在线播放 精品久久久久久中文字幕专区 亚洲人成人毛片无遮挡 国产一起色一起爱 国产香蕉精品视频在 九九热免费观看 日韩亚洲欧美一区 九九热精品在线观看 精品久久久久久中文字幕专区 亚洲欧美自拍偷拍 国产精品每日更新 久久久久国产一级毛片高清板 久久天天躁狠狠躁夜夜中文字幕 久久精品片 日韩在线毛片 国产成人精品本亚洲 国产成人精品一区二区三区 九九热在线观看 国产r级在线观看 国产欧美日韩精品高清二区综合区 韩国电影一区二区 国产精品毛片va一区二区三区 五月婷婷伊人网 久久一区二区三区免费 一本色道久久综合狠狠躁篇 亚洲综合色站 国产尤物在线观看 亚洲一区亚洲二区 免费在线视频一区 欧洲精品视频在线观看 日韩中文字幕a 中文字幕日本在线mv视频精品 91精品在线免费视频 精品国产免费人成在线观看 精品a级片 中文字幕日本在线mv视频精品 日韩在线精品视频 婷婷丁香色 91精品国产高清久久久久 国产成人精品日本亚洲直接 五月综合视频 欧美日韩在线亚洲国产人 精液呈暗黄色 亚洲乱码一区 久久精品中文字幕不卡一二区 亚洲天堂精品在线 激情婷婷综合 国产免费久久精品久久久 国产精品亚洲二区在线 久久免费播放视频 五月婷婷丁香综合 在线亚洲欧美日韩 久久免费精品高清麻豆 精品久久久久久中文字幕 亚洲一区网站 国产精品福利社 日韩中文字幕免费 亚洲综合丝袜 91精品在线播放 国产精品18 亚洲日日夜夜 伊人久久大香线蕉综合影 亚洲精品中文字幕乱码影院 亚洲一区二区黄色 亚洲第一页在线视频 一区二区在线观看视频 国产成人福利精品视频 亚洲高清二区 国内成人免费视频 精品亚洲性xxx久久久 国产精品合集一区二区三区 97av免费视频 国产一起色一起爱 国产区久久 国产资源免费观看 99精品视频免费 国产成人一级 国产精品九九免费视频 欧美91精品久久久久网免费 99热国产免费 久久精品色 98精品国产综合久久 久久精品播放 中文字幕视频免费 国产欧美日韩一区二区三区在线 精品久久蜜桃 国产小视频精品 一本色道久久综合狠狠躁篇 91在线免费观看 亚洲精品区 伊人成综合网 伊人热人久久中文字幕 伊人黄色片 99国产精品热久久久久久夜夜嗨 久久免费精品视频 亚洲一区二区三区高清不卡 久久久久国产一级毛片高清板 国产片一区二区三区 久久狠狠干 99久久婷婷国产综合精品电影 国产99区 国产精品成人久久久久 久久狠狠干 青青国产在线观看 亚洲高清国产拍精品影院 国产精品一区二区av 九九热在线免费视频 伊人久久国产 国产精品久久久久久久久久一区 在线观看免费视频一区 国产精品自在在线午夜区app 国产精品综合色区在线观看 国产毛片久久久久久国产毛片 97国产免费全部免费观看 国产精品每日更新 国产尤物视频在线 九九视频这里只有精品99 一本一道久久a久久精品综合 久久综合给会久久狠狠狠 国产成人精品男人的天堂538 欧美一区二区高清 毛片一区二区三区 国产欧美日韩在线观看一区二区三区 在线国产二区 欧美不卡网 91在线精品中文字幕 在线国产福利 国内精品91久久久久 91亚洲福利 日韩欧美国产中文字幕 91久久精品国产性色也91久久 亚洲性久久久影院 欧美精品1区 国产热re99久久6国产精品 九九热免费观看 国产精品欧美日韩 久久久久国产一级毛片高清板 久久国产经典视频 日韩欧美亚洲国产一区二区三区 欧美亚洲综合另类在线观看 国产精品自在在线午夜区app 97中文字幕在线观看 视频一二三区 精品国产一区在线观看 国产欧美日韩在线一区二区不卡 欧美一区二三区 伊人成人在线观看 国内精品91久久久久 97在线亚洲 国产在线不卡一区 久久久全免费全集一级全黄片 国产精品v欧美精品∨日韩 亚洲毛片网站 在线不卡一区二区 99re热在线视频 久久激情网 国产毛片一区二区三区精品 久久亚洲综合色 中文字幕视频免费 国产视频亚洲 婷婷伊人久久 国产一区二区免费播放 久久99国产精品成人欧美 99国产在线视频 国产成人免费视频精品一区二区 国产不卡一区二区三区免费视 国产码欧美日韩高清综合一区 久久精品国产主播一区二区 国产一区电影 久久精品国产夜色 国产精品国产三级国产 日韩一区二区三区在线 久久97久久97精品免视看 久久国产免费一区二区三区 伊人久久大香线蕉综合电影网 99re6久精品国产首页 久久激情网 亚洲成人高清在线 国产精品网址 国产成人精品男人的天堂538 香蕉国产综合久久猫咪 国产专区中文字幕 91麻豆精品国产高清在线 久久国产经典视频 国产精品成人va在线观看 国产精品爱啪在线线免费观看 日本精品久久久久久久久免费 亚洲综合一区二区三区 久久五月网 精品国产网红福利在线观看 久久综合亚洲伊人色 亚洲国产精品久久久久久网站 在线日韩国产 99国产精品热久久久久久夜夜嗨 国产综合精品在线 国产区福利 精品亚洲综合久久中文字幕 国产制服丝袜在线 毛片在线播放网站 在线观看免费视频一区 国产精品久久久精品三级 亚洲国产电影在线观看 最新日韩欧美不卡一二三区 狠狠综合久久综合鬼色 日本精品1在线区 国产日韩一区二区三区在线播放 欧美日韩精品在线播放 亚洲欧美日韩国产一区二区三区精品 久久综合久久网 婷婷六月激情在线综合激情 亚洲乱码一区 国产专区91 97av视频在线观看 精品久久久久久中文字幕 久久五月视频 国产成人福利精品视频 国产精品网址 中文字幕视频在线 精品一区二区三区免费视频 伊人手机在线视频 亚洲精品中文字幕乱码 国产在线视频www色 色噜噜国产精品视频一区二区 精品亚洲成a人在线观看 国产香蕉尹人综合在线 成人免费一区二区三区在线观看 国产不卡一区二区三区免费视 欧美精品久久天天躁 国产专区中文字幕 久久精品国产免费中文 久久精品国产免费一区 久久无码精品一区二区三区 国产欧美另类久久久精品免费 欧美精品久久天天躁 亚洲精品在线视频 国产视频91在线 91精品福利一区二区三区野战 日韩中文字幕免费 国产精品99一区二区三区 欧美成人高清性色生活 国产精品系列在线观看 亚洲国产福利精品一区二区 国产成人在线小视频 国产精品久久久久免费 99re热在线视频 久久久久久久综合 一区二区国产在线播放 成人国产在线视频 亚洲精品乱码久久久久 欧美日韩一区二区综合 精品久久久久免费极品大片 中文字幕视频二区 激情粉嫩精品国产尤物 国产成人精品一区二区视频 久久精品中文字幕首页 亚洲高清在线 国产精品亚洲一区二区三区 伊人久久艹 中文在线亚洲 国产精品一区二区在线播放 国产精品九九免费视频 亚洲二区在线播放 亚洲狠狠婷婷综合久久久久网站 亚洲欧美日韩网站 日韩成人精品 亚洲国产一区二区三区青草影视 91精品国产福利在线观看 国产精品久久久久久久久99热 国产一区二区精品尤物 久碰香蕉精品视频在线观看 亚洲日日夜夜 在线不卡一区二区 国产午夜亚洲精品 九九热在线视频观看这里只有精品 伊人手机在线视频 91免费国产精品 日韩欧美中字 91精品国产91久久久久 国产全黄三级播放 视频一区二区三区免费观看 国产开裆丝袜高跟在线观看 国产成人欧美 激情综合丝袜美女一区二区 国产成人亚洲综合无 欧美精品一区二区三区免费观看 欧美亚洲国产日韩 日韩亚州 国产欧美日韩精品高清二区综合区 亚洲午夜国产片在线观看 精品久久久久久中文字幕 欧美精品1区 久久伊人久久亚洲综合 亚洲欧美日韩精品 国产成人精品久久亚洲高清不卡 久久福利影视 国产精品99精品久久免费 久久久久免费精品视频 国产日产亚洲精品 亚洲国产午夜电影在线入口 精品无码一区在线观看 午夜国产精品视频 亚洲一级片免费 伊人久久大香线蕉综合影 国产精品久久影院 久碰香蕉精品视频在线观看 www.欧美精品 在线小视频国产 亚洲国产天堂久久综合图区 欧美一区二区三区不卡 日韩美女福利视频 九九精品免视频国产成人 不卡国产00高中生在线视频 亚洲第一页在线视频 欧美日韩在线播放成人 99re视频这里只有精品 国产精品91在线 精品乱码一区二区三区在线 国产区久久 91麻豆精品国产自产在线观看一区 日韩精品成人在线 九九热在线观看 国产精品久久不卡日韩美女 欧美一区二区三区综合色视频 欧美精品免费一区欧美久久优播 国产精品网址 国产专区中文字幕 国产精品欧美亚洲韩国日本久久 日韩美香港a一级毛片 久久精品123 欧美一区二区三区免费看 99r在线视频 亚洲精品国产字幕久久vr 国产综合激情在线亚洲第一页 91免费国产精品 日韩免费小视频 亚洲国产精品综合一区在线 国产亚洲第一伦理第一区 在线亚洲精品 国产精品一区二区制服丝袜 国产在线成人精品 九九精品免视频国产成人 亚洲国产网 欧美日韩亚洲一区二区三区在线观看 在线亚洲精品 欧美一区二区三区高清视频 国产成人精品男人的天堂538 欧美日韩在线观看区一二 亚洲欧美一区二区久久 久久精品中文字幕首页 日本高清www午夜视频 久久精品国产免费 久久999精品 亚洲国产精品欧美综合 88国产精品视频一区二区三区 91久久偷偷做嫩草影院免费看 国产精品夜色视频一区二区 欧美日韩导航 国产成人啪精品午夜在线播放 一区二区视频在线免费观看 99久久精品国产自免费 精液呈暗黄色 久久99国产精品 日本精品久久久久久久久免费 精品国产97在线观看 99re视频这里只有精品 国产视频91在线 999av视频 亚洲美女视频一区二区三区 久久97久久97精品免视看 亚洲国产成人久久三区 99久久亚洲国产高清观看 日韩毛片在线视频 综合激情在线 91福利一区二区在线观看 一区二区视频在线免费观看 激情粉嫩精品国产尤物 国产成人精品曰本亚洲78 国产成人精品本亚洲 国产精品成人免费视频 国产成人啪精品视频免费软件 久久精品国产亚洲妲己影院 国产精品成人久久久久久久 久久大香线蕉综合爱 欧美一区二区三区高清视频 99热国产免费 在线观看欧美国产 91精品视频在线播放 国产精品福利社 欧美精品一区二区三区免费观看 国产一区二区免费视频 国产午夜精品一区二区 精品视频在线观看97 91精品福利久久久 国产一区福利 国产综合激情在线亚洲第一页 国产精品久久久久久久久久久不卡 九色国产 在线日韩国产 黄网在线观看 亚洲一区小说区中文字幕 中文字幕丝袜 日本二区在线观看 日本国产一区在线观看 欧美日韩一区二区三区久久 欧美精品亚洲精品日韩专 国产日产亚洲精品 久久综合九色综合欧美播 亚洲国产欧美无圣光一区 欧美视频区 亚洲乱码视频在线观看 久久无码精品一区二区三区 九九热精品免费视频 久久99精品久久久久久牛牛影视 国产精品成久久久久三级 国产一区福利 午夜国产精品视频 日本二区在线观看 99久久网站 国产亚洲天堂 精品国产一区二区三区不卡 亚洲国产日韩在线一区 国产成人综合在线观看网站 久久免费高清视频 欧美在线导航 午夜精品久久久久久99热7777 欧美久久综合网 国产小视频精品 国产尤物在线观看 亚洲国产精品综合一区在线 欧美一区二区三区不卡视频 欧美黑人欧美精品刺激 日本福利在线观看 久久国产偷 国产手机精品一区二区 国产热re99久久6国产精品 国产高清啪啪 欧美亚洲国产成人高清在线 国产在线第三页 亚洲综合一区二区三区 99r在线视频 99精品久久久久久久婷婷 国产精品乱码免费一区二区 国产在线精品福利91香蕉 国产尤物视频在线 五月婷婷亚洲 中文字幕久久综合伊人 亚洲精品一级毛片 99国产精品电影 在线视频第一页 久久99国产精品成人欧美 国产白白视频在线观看2 成人精品一区二区www 亚洲成人网在线观看 麻豆91在线视频 色综合合久久天天综合绕视看 久久精品国产免费高清 国产不卡一区二区三区免费视 欧美国产中文 99精品欧美 九九在线精品 国产中文字幕在线免费观看 国产一区中文字幕在线观看 国产成人一级 国产精品一区二区制服丝袜 国产一起色一起爱 亚洲精品成人在线 亚洲欧美精品在线 国产欧美自拍视频 99精品久久久久久久婷婷 久99视频 国产热re99久久6国产精品 视频一区亚洲 国产精品视频分类 国产精品成在线观看 99re6久精品国产首页 亚洲在成人网在线看 亚洲国产日韩在线一区 久久国产三级 日韩国产欧美 欧美在线一区二区三区 国产精品美女一级在线观看 成人午夜免费福利视频 亚洲天堂精品在线 91精品国产手机 欧美日韩视频在线播放 狠狠综合久久综合鬼色 九一色视频 青青视频国产 亚洲欧美自拍一区 中文字幕天天躁日日躁狠狠躁97 日韩免费大片 996热视频 伊人成综合网 亚洲天堂欧美 日韩精品亚洲人成在线观看 久久综合给会久久狠狠狠 日韩精品亚洲人成在线观看 日韩国产欧美 亚洲成aⅴ人片在线影院八 亚洲精品1区 99久久精品免费 国产精品高清在线观看 国产精品久久久免费视频 在线亚洲欧美日韩 91在线看视频 国产精品96久久久久久久 欧美日韩国产成人精品 91在线亚洲 热久久亚洲 国产精品美女免费视频观看 日韩在线毛片 亚洲永久免费视频 九九免费在线视频 亚洲一区网站 日本高清二区视频久二区 精品国产美女福利在线 伊人久久艹 国产精品久久久久三级 欧美成人精品第一区二区三区 99久久精品国产自免费 在线观看日韩一区 国产中文字幕一区 成人免费午夜视频 欧美日韩另类在线 久久99国产精品成人欧美 色婷婷中文网 久久天天躁夜夜躁狠狠躁2020 欧美成人伊人久久综合网 国产精品福利资源在线 国产伦精品一区二区三区高清 国产精品亚洲综合色区韩国 亚洲一区欧美日韩 色综合视频 国语自产精品视频在线区 国产高清a 成人国内精品久久久久影 国产在线精品香蕉综合网一区 国产不卡在线看 国产成人精品精品欧美 国产欧美日韩综合精品一区二区三区 韩国电影一区二区 国产在线视频www色 91中文字幕在线一区 国产人成午夜免视频网站 亚洲综合一区二区三区 色综合视频一区二区观看 久久五月网 九九热精品在线观看 国产一区二区三区国产精品 99久热re在线精品996热视频 亚洲国产网 在线视频亚洲一区 日韩字幕一中文在线综合 国产高清一级毛片在线不卡 精品国产色在线 国产高清视频一区二区 精品日本久久久久久久久久 亚洲国产午夜精品乱码 成人免费国产gav视频在线 日韩欧美一区二区在线观看 欧美曰批人成在线观看 韩国电影一区二区 99re这里只有精品6 日韩精品一区二区三区视频 99re6久精品国产首页 亚洲欧美一区二区三区导航 欧美色图一区二区三区 午夜精品视频在线观看 欧美激情在线观看一区二区三区 亚洲热在线 成人国产精品一区二区网站 亚洲一级毛片在线播放 亚洲一区小说区中文字幕 亚洲午夜久久久久影院 国产自产v一区二区三区c 国产精品视频免费 久久调教视频 国产成人91激情在线播放 国产精品欧美亚洲韩国日本久久 久久亚洲日本不卡一区二区 91中文字幕网 成人国产在线视频 国产视频91在线 欧美成人精品第一区二区三区 国产精品福利在线 久久综合九色综合精品 欧美一区二区三区精品 久久国产综合尤物免费观看 久久99青青久久99久久 日韩精品免费 久久国产精品999 91亚洲视频在线观看 国产精品igao视频 色综合区 在线亚洲欧国产精品专区 国产一区二区三区在线观看视频 亚洲精品成人在线 一区二区国产在线播放 中文在线亚洲 亚洲精品第一国产综合野 国产一区二区精品久久 一区二区三区四区精品视频 99热精品久久 中文字幕视频二区 国产成人精品男人的天堂538 99精品影视 美女福利视频一区二区 久久午夜夜伦伦鲁鲁片 综合久久久久久久综合网 国产精品国产欧美综合一区 国产99视频在线观看 国产亚洲女在线精品 婷婷影院在线综合免费视频 国产亚洲3p一区二区三区 91成人爽a毛片一区二区 亚洲一区二区高清 国产欧美亚洲精品第二区首页 欧美日韩导航 亚洲高清二区 欧美激情观看一区二区久久 日韩毛片在线播放 亚洲欧美日韩高清中文在线 亚洲日本在线播放 国产精品一区二区制服丝袜 精品国产一区二区三区不卡 国产不卡在线看 国产欧美网站 四虎永久在线观看视频精品 国产黄色片在线观看 夜夜综合 一本色道久久综合狠狠躁篇 欧美亚洲综合另类在线观看 国产91在线看 伊人久久国产 欧美一区二区在线观看免费网站 国产精品久久久久三级 久久福利 日韩中文字幕a 亚洲午夜久久久久影院 91在线高清视频 国产亚洲一区二区三区啪 久久人精品 国产精品亚洲午夜一区二区三区 综合久久久久久 久久伊人一区二区三区四区 国产综合久久久久久 日韩一区精品视频在线看 国产精品日韩欧美制服 日本精品1在线区 99re视频 无码av免费一区二区三区试看 国产视频1区 日韩欧美中文字幕一区 日本高清中文字幕一区二区三区a 亚洲国产欧美无圣光一区 国产在线视频一区二区三区 欧美国产第一页 在线亚洲欧美日韩 日韩中文字幕第一页 在线不卡一区二区 伊人久久青青 国产精品一区二区在线播放 www.五月婷婷 麻豆久久婷婷国产综合五月 亚洲精品区 久久国产欧美另类久久久 99在线视频免费 伊人久久中文字幕久久cm 久久精品成人免费看 久久这里只有精品首页 88国产精品视频一区二区三区 中文字幕日本在线mv视频精品 国产在线精品成人一区二区三区 伊人精品线视天天综合 亚洲一区二区黄色 国产尤物视频在线 亚洲精品99久久久久中文字幕 国产一区二区三区免费观看 伊人久久大香线蕉综合电影网 国产成人精品区在线观看 日本精品一区二区三区视频 日韩高清在线二区 久久免费播放视频 一区二区成人国产精品 国产精品免费精品自在线观看 亚洲精品视频二区 麻豆国产精品有码在线观看 精品日本一区二区 亚洲欧洲久久 久久中文字幕综合婷婷 中文字幕视频在线 国产成人精品综合在线观看 91精品国产91久久久久福利 精液呈暗黄色 香蕉国产综合久久猫咪 国产专区精品 亚洲精品无码不卡 国产永久视频 亚洲成a人片在线播放观看国产 一区二区国产在线播放 亚洲一区二区黄色 欧美日韩在线观看视频 亚洲精品另类 久久国产综合尤物免费观看 国产一区二区三区国产精品 高清视频一区 国产精品igao视频 国产精品资源在线 久久综合精品国产一区二区三区 www.五月婷婷 精品色综合 99热国产免费 麻豆福利影院 亚洲伊人久久大香线蕉苏妲己 久久电影院久久国产 久久精品伊人 在线日韩理论午夜中文电影 亚洲国产欧洲综合997久久 伊人国产精品 久草国产精品 欧美一区精品二区三区 亚洲成人高清在线 91免费国产精品 日韩精品福利在线 国产一线在线观看 国产不卡在线看 久久99青青久久99久久 亚洲精品亚洲人成在线播放 99久久免费看国产精品 国产日本在线观看 青草国产在线视频 麻豆久久婷婷国产综合五月 国产中文字幕一区 91久久精品国产性色也91久久 国产一区a 国产欧美日韩成人 国产亚洲女在线精品 一区二区美女 中文字幕在线2021一区 在线小视频国产 久久这里只有精品首页 国产在线第三页 欧美日韩中文字幕 在线亚洲+欧美+日本专区 精品国产一区二区三区不卡 久久这里精品 欧美在线va在线播放 精液呈暗黄色 91精品国产手机 91在线免费播放 欧美视频亚洲色图 欧美国产日韩精品 日韩高清不卡在线 精品视频免费观看 欧美日韩一区二区三区四区 国产欧美亚洲精品第二区首页 亚洲韩精品欧美一区二区三区 国产精品视频免费 在线精品小视频 久久午夜夜伦伦鲁鲁片 国产无套在线播放 久热这里只精品99re8久 欧美久久久久 久久香蕉国产线看观看精品蕉 国产成人精品男人的天堂538 亚洲人成网站色7799在线观看 日韩在线第二页 一本色道久久综合狠狠躁篇 国产一区二区三区不卡在线观看 亚洲乱码在线 在线观看欧美国产 久久福利青草精品资源站免费 国产玖玖在线观看 在线亚洲精品 亚洲成aⅴ人在线观看 精品91在线 欧美一区二三区 日韩中文字幕视频在线 日本成人一区二区 日韩免费专区 国内精品在线观看视频 久久国产综合尤物免费观看 国产精品系列在线观看 一本一道久久a久久精品综合 亚洲免费播放 久久精品国产免费 久久人精品 亚洲毛片网站 亚洲成a人一区二区三区 韩国福利一区二区三区高清视频 亚洲精品天堂在线 一区二区三区中文字幕 亚洲国产色婷婷精品综合在线观看 亚洲国产成人久久笫一页 999国产视频 国产精品香港三级在线电影 欧美日韩一区二区三区四区 日韩国产欧美 国产精品99一区二区三区 午夜国产精品理论片久久影院 亚洲精品中文字幕麻豆 亚洲国产高清视频 久久免费手机视频 日韩a在线观看 五月婷婷亚洲 亚洲精品中文字幕麻豆 中文字幕丝袜 www国产精品 亚洲天堂精品在线 亚洲乱码一区 国产日韩欧美三级 久久999精品 伊人热人久久中文字幕 久热国产在线视频 国产欧美日韩在线观看一区二区三区 国产一二三区在线 日韩国产欧美 91精品国产91久久久久 亚洲一区小说区中文字幕 精品一区二区免费视频 国产精品视频免费 国产精品亚洲综合色区韩国 亚洲国产精品成人午夜在线观看 欧美国产日韩精品 中文字幕精品一区二区精品