應變監測受損索載荷廣義位移遞進式識別方法
【專利摘要】應變監測受損索載荷廣義位移遞進式識別方法基于應變監測、通過監測索結構溫度和環境溫度來決定是否需要更新索結構的力學計算基準模型，得到計入索結構溫度和環境溫度的索結構的力學計算基準模型，在此模型的基礎上計算獲得單位損傷被監測量數值變化矩陣。依據被監測量當前數值向量同被監測量當前初始數值向量、單位損傷被監測量數值變化矩陣和待求的被評估對象當前名義損傷向量間存在的近似線性關系算出被評估對象當前名義損傷向量的非劣解，據此可以在有溫度變化時，識別出支座廣義位移、載荷變化量和受損索。
【專利說明】應變監測受損索載荷廣義位移遞進式識別方法
【技術領域】
[0001]斜拉橋、懸索橋、桁架結構等結構有一個共同點，就是它們有許多承受拉伸載荷的部件，如斜拉索、主纜、吊索、拉桿等等，該類結構的共同點是以索、纜或僅承受拉伸載荷的桿件為支承部件，為方便起見，本方法將該類結構表述為“索結構”，并將索結構的所有承載索、承載纜，及所有僅承受軸向拉伸或軸向壓縮載荷的桿件(又稱為二力桿件)，為方便起見統一稱為“索系統”，本方法中用“支承索”這一名詞指稱承載索、承載纜及僅承受軸向拉伸或軸向壓縮載荷的桿件，有時簡稱為“索”，所以在后面使用“索”這個字的時候，對桁架結構實際就是指二力桿件。在結構服役過程中，對支承索或索系統的健康狀態的正確識別關系到整個索結構的安全。在環境溫度發生變化時，索結構的溫度一般也會隨著發生變化，在索結構溫度發生變化時，索結構支座可能發生廣義位移，索結構承受的載荷也可能發生變化，同時索結構的健康狀態也可能在發生變化，在這種復雜條件下，本方法基于應變監測(本方法將被監測的應變稱為“被監測量”)來識別支座廣義位移、受損索和索結構承受的載荷的變化量，屬工程結構健康監測領域。
【背景技術】
[0002]剔除載荷變化、索結構支座廣義位移和結構溫度變化對索結構健康狀態識別結果的影響，從而準確地識別結構的健康狀態的變化，是目前迫切需要解決的問題；剔除載荷變化、索結構健康狀態變化和結構溫度變化對索結構支座廣義位移識別結果的影響，從而準確地識別索結構支座廣義位移，也是目前迫切需要解決的問題；同樣的，剔除結構溫度變化、索結構支座廣義位移和結構健康狀態變化對結構承受的載荷的變化量的識別結果的影響，對結構安全同樣具有重要意義?；诮Y構健康監測技術，本方法公開了一種解決這三個問題的有效方法。
【發明內容】
[0003]技術問題:本方法公開了一種方法，實現了已有方法不可能具備的三種功能，分別是，一、剔除支座廣義位移、載荷變化和結構溫度變化對索結構健康狀態識別結果的影響，從而準確地識別出支承索的健康狀態；二、本方法在識別出受損索的同時，還能同時識別出載荷的變化，即本方法能夠剔除支座廣義位移、結構溫度變化和支承索健康狀態變化的影響，實現載荷變化程度的正確識別；三、本方法還能夠剔除載荷變化、索結構健康狀態變化和結構溫度變化對索結構支座廣義位移識別結果的影響，從而準確地識別索結構支座廣義位移。
[0004]技術方案:在本方法中，用“支座空間坐標”指稱支座關于笛卡爾直角坐標系的X、Y、z軸的坐標，也可以說成是支座關于X、Y、Z軸的空間坐標，支座關于某一個軸的空間坐標的具體數值稱為支座關于該軸的空間坐標分量，本方法中也用支座的一個空間坐標分量表達支座關于某一個軸的空間坐標的具體數值；用“支座角坐標”指稱支座關于χ、、ζ軸的角坐標，支座關于某一個軸的角坐標的具體數值稱為支座關于該軸的角坐標分量，本方法中也用支座的一個角坐標分量表達支座關于某一個軸的角坐標的具體數值；用“支座廣義坐標”指稱支座角坐標和支座空間坐標全體，本方法中也用支座的一個廣義坐標分量表達支座關于一個軸的空間坐標或角坐標的具體數值；支座關于Χ、y、ζ軸的坐標的改變稱為支座線位移，也可以說支座空間坐標的改變稱為支座線位移，本方法中也用支座的一個線位移分量表達支座關于某一個軸的線位移的具體數值；支座關于χ、y、ζ軸的角坐標的改變稱為支座角位移，本方法中也用支座的一個角位移分量表達支座關于某一個軸的角位移的具體數值；支座廣義位移指稱支座線位移和支座角位移全體，本方法中也用支座的一個廣義位移分量表達支座關于某一個軸的線位移或角位移的具體數值；支座線位移也可稱為平移位移，支座沉降是支座線位移或平移位移在重力方向的分量。
[0005]物體、結構承受的外力可稱為載荷，載荷包括面載荷和體積載荷。面載荷又稱表面載荷，是作用于物體表面的載荷，包括集中載荷和分布載荷兩種。體積載荷是連續分布于物體內部各點的載荷，如物體的自重和慣性力。
[0006]集中載荷分為集中力和集中力偶兩種，在坐標系中，例如在笛卡爾直角坐標系中，一個集中力可以分解成三個分量，同樣的，一個集中力偶也可以分解成三個分量，如果載荷實際上是集中載荷，在本方法中將一個集中力分量或一個集中力偶分量稱為一個載荷，此時載荷的變化具體化為一個集中力分量或一個集中力偶分量的變化。
[0007]分布載荷分為線分布載荷和面分布載荷，分布載荷的描述至少包括分布載荷的作用區域和分布載荷的 大小，分布載荷的大小用分布集度來表達，分布集度用分布特征(例如均布、正弦函數等分布特征)和幅值來表達(例如兩個分布載荷都是均布，但其幅值不同，可以均布壓力為例來說明幅值的概念:同一個結構承受兩個不同的均布壓力，兩個分布載荷都是均布載荷，但一個分布載荷的幅值是IOMPa,另一個分布載荷的幅值是50MPa)。如果載荷實際上是分布載荷，本方法談論載荷的變化時，實際上是指分布載荷分布集度的幅值的改變，而分布載荷的作用區域和分布集度的分布特征是不變的。在坐標系中，一個分布載荷可以分解成若干個分量，如果這分布載荷的若干個分量的各自的分布集度的幅值發生變化，且變化的比率不全部相同，那么在本方法中把這若干個分布載荷的分量看成同樣數量的獨立的分布載荷，此時一個載荷就代表一個分布載荷的分量，也可以將其中分布集度的幅值變化比率相同的分量合成為一個分布載荷或稱為一個載荷。
[0008]體積載荷是連續分布于物體內部各點的載荷，如物體的自重和慣性力，體積載荷的描述至少包括體積載荷的作用區域和體積載荷的大小，體積載荷的大小用分布集度來表達，分布集度用分布特征(例如均布、線性函數等分布特征)和幅值來表達(例如兩個體積載荷都是均布，但其幅值不同，可以自重為例來說明幅值的概念:同一個結構的兩個部分的材料不同，故密度不同，所以雖然這兩個部分所受的體積載荷都是均布的，但一個部分所受的體積載荷的幅值可能是10kN/m3，另一個部分所受的體積載荷的幅值是50kN/m3)。如果載荷實際上是體積載荷，在本方法中實際處理的是體積載荷分布集度的幅值的改變，而體積載荷的作用區域和分布集度的分布特征是不變的，此時在本方法中提到載荷的改變時實際上是指體積載荷的分布集度的幅值的改變，此時，發生變化的載荷是指那些分布集度的幅值發生變化的體積載荷。在坐標系中，一個體積載荷可以分解成若干個分量(例如在笛卡爾直角坐標系中，體積載荷可以分解成關于坐標系的三個軸的分量，也就是說，在笛卡爾直角坐標系中體積載荷可以分解成三個分量)，如果這體積載荷的若干個分量的各自的分布集度的幅值發生變化，且變化的比率不全部相同，那么在本方法中把這若干個體積載荷的分量看成同樣數量的獨立的載荷，也可以將其中分布集度的幅值變化比率相同的體積載荷分量合成為一個體積載荷或稱為一個載荷。
[0009]當載荷具體化為集中載荷時，在本方法中，“載荷單位變化”實際上是指“集中載荷的單位變化”，類似的，“載荷變化”具體指“集中載荷的大小的變化”，“載荷變化量”具體指“集中載荷的大小的變化量”，“載荷變化程度”具體指“集中載荷的大小的變化程度”，“載荷的實際變化量”是指“集中載荷的大小的實際變化量”，“發生變化的載荷”是指“大小發生變化的集中載荷”，簡單地說，此時“某某載荷的某某變化”是指“某某集中載荷的大小的某某變化”。
[0010]當載荷具體化為分布載荷時，在本方法中，“載荷單位變化”實際上是指“分布載荷的分布集度的幅值的單位變化”，而分布載荷的分布特征是不變的，類似的，“載荷變化”具體指“分布載荷的分布集度的幅值的變化”，而分布載荷的分布特征是不變的，“載荷變化量”具體指“分布載荷的分布集度的幅值的變化量”，“載荷變化程度”具體指“分布載荷的分布集度的幅值的變化程度”，“載荷的實際變化量”具體指“分布載荷的分布集度的幅值的實際變化量”，“發生變化的載荷”是指“分布集度的幅值發生變化的分布載荷”，簡單地說，此時“某某載荷的某某變化”是指“某某分布載荷的分布集度的幅值的某某變化”，而所有分布載荷的作用區域和分布集度的分布特征是不變的。
[0011]當載荷具體化為體積載荷時，在本方法中，“載荷單位變化”實際上是指“體積載荷的分布集度的幅值的單位變化”，類似的，“載荷變化”是指“體積載荷的分布集度的幅值的變化”，“載荷變化量”是指“體積載荷的分布集度的幅值的變化量”，“載荷變化程度”是指“體積載荷的分布集度的幅值的變化程度”，“載荷的實際變化量”是指“體積載荷的分布集度的幅值的實際變化量”，“發生變化的載荷”是指“分布集度的幅值發生變化的體積載荷”，簡單地說，“某某載荷的某某變化”是指“某某體積載荷的分布集度的幅值的某某變化”，而所有體積載荷的作用區域和分布集度的分布特征是不變的。
[0012]本方法具體包括:
[0013]a.為敘述方便起見，本方法統一稱被評估的支座廣義位移分量、支承索和載荷為被評估對象，設被評估的支座廣義位移分量的數量、支承索的數量和載荷的數量之和為N，即被評估對象的數量為N ;確定被評估對象的編號規則，按此規則將索結構中所有的被評估對象編號，該編號在后續步驟中將用于生成向量和矩陣；本方法用變量k表示這一編號，k=l, 2，3，…，N ;確定指定的被監測點，被監測點即表征索結構應變信息的所有指定點，并給所有指定點編號；確定被監測點的被監測的應變方向，并給所有指定的被監測應變編號，“被監測應變編號”在后續步驟中將用于生成向量和矩陣，“索結構的全部被監測的應變數據”由上述所有被監測應變組成；本方法將“索結構的被監測的應變數據”簡稱為“被監測量”;所有被監測量的數量之和記為M，M不得小于N ;本方法中對同一個量實時監測的任何兩次測量之間的時間間隔不得大于30分鐘，測量記錄數據的時刻稱為實際記錄數據時刻；物體、結構承受的外力可稱為載荷，載荷包括面載荷和體積載荷；面載荷又稱表面載荷，是作用于物體表面的載荷，包括集中載荷和分布載荷兩種；體積載荷是連續分布于物體內部各點的載荷，包括物體的自重和慣性力在內；集中載荷分為集中力和集中力偶兩種，在包括笛卡爾直角坐標系在內的坐標系中，一個集中力可以分解成三個分量，同樣的，一個集中力偶也可以分解成三個分量，如果載荷實際上是集中載荷，在本方法中將一個集中力分量或一個集中力偶分量計為或統計為一個載荷，此時載荷的變化具體化為一個集中力分量或一個集中力偶分量的變化；分布載荷分為線分布載荷和面分布載荷，分布載荷的描述至少包括分布載荷的作用區域和分布載荷的大小，分布載荷的大小用分布集度來表達，分布集度用分布特征和幅值來表達；如果載荷實際上是分布載荷，本方法談論載荷的變化時，實際上是指分布載荷分布集度的幅值的改變，而所有分布載荷的作用區域和分布集度的分布特征是不變的；在包括笛卡爾直角坐標系在內的坐標系中，一個分布載荷可以分解成三個分量，如果這分布載荷的三個分量的各自的分布集度的幅值發生變化，且變化的比率不全部相同，那么在本方法中把這分布載荷的三個分量計為或統計為三個分布載荷，此時一個載荷就代表分布載荷的一個分量；體積載荷是連續分布于物體內部各點的載荷，體積載荷的描述至少包括體積載荷的作用區域和體積載荷的大小，體積載荷的大小用分布集度來表達，分布集度用分布特征和幅值來表達；如果載荷實際上是體積載荷，在本方法中實際處理的是體積載荷分布集度的幅值的改變，而所有體積載荷的作用區域和分布集度的分布特征是不變的，此時在本方法中提到載荷的改變時實際上是指體積載荷的分布集度的幅值的改變，此時，發生變化的載荷是指那些分布集度的幅值發生變化的體積載荷；在包括笛卡爾直角坐標系在內的坐標系中，一個體積載荷可以分解成三個分量，如果這體積載荷的三個分量的各自的分布集度的幅值發生變化，且變化的比率不全部相同，那么在本方法中把這體積載荷的三個分量計為或統計為三個分布載荷；
[0014]b.本方法定義“本方法的索結構的溫度測量計算方法”按步驟bl至b3進行；
[0015]bl:查詢或實測得到索結構組成材料及索結構所處環境的隨溫度變化的傳熱學參數，利用索結構的設計圖、竣工圖和索結構的幾何實測數據，利用這些數據和參數建立索結構的傳熱學計算模型；查詢索結構所在地不少于2年的近年來的氣象資料，統計得到這段時間內的陰天數量記為T個陰天，在本方法中將白天不能見到太陽的一整日稱為陰天，統計得到T個陰天中每一個陰天的O時至次日日出時刻后30分鐘之間的最高氣溫與最低氣溫，日出時刻是指根據地球自轉和公轉規律確定的氣象學上的日出時刻，不表示當天一定可以看見太陽，可以查詢資料或通過常規氣象學計算得到所需的每一日的日出時刻，每一個陰天的O時至次日日出時刻后30分鐘之間的最高氣溫減去最低氣溫稱為該陰天的日氣溫的最大溫差，有T個陰天， 就有T個陰天的日氣溫的最大溫差，取T個陰天的日氣溫的最大溫差中的最大值為參考日溫差，參考日溫差記為Λ?；;查詢索結構所在地和所在海拔區間不少于2年的近年來的氣象資料或實測得到索結構所處環境的溫度隨時間和海拔高度的變化數據和變化規律，計算得到索結構所在地和所在海拔區間不少于2年的近年來的索結構所處環境的溫度關于海拔高度的最大變化率ATh,為方便敘述取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個索結構沿厚度溫度計算數據” 為B和E的乘積，本方法中在每一個選取的海拔高度處共有BE個“相同海拔高度索結構沿厚度的溫度分布數據”;在索結構所在地按照氣象學測量氣溫要求選取一個位置，將在此位置實測得到符合氣象學測量氣溫要求的索結構所在環境的氣溫；在索結構所在地的空曠無遮擋處選取一個位置，該位置應當在全年的每一日都能得到該地所能得到的該日的最充分的日照，在該位置安放一塊碳鋼材質的平板，稱為參考平板，參考平板與地面不可接觸，參考平板離地面距離不小于1.5米，該參考平板的一面向陽，稱為向陽面，參考平板的向陽面是粗糙的和深色的，參考平板的向陽面應當在全年的每一日都能得到一塊平板在該地所能得到的該日的最充分的日照，參考平板的非向陽面覆有保溫材料，將實時監測得到參考平板的向陽面的溫度；
[0016]b2:實時監測得到上述R個索結構表面點的R個索結構表面溫度實測數據，同時實時監測得到前面定義的索結構沿厚度的溫度分布數據，同時實時監測得到符合氣象學測量氣溫要求的索結構所在環境的氣溫數據；通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的索結構所在環境的氣溫實測數據序列，索結構所在環境的氣溫實測數據序列由當日日出時刻到次日日出時刻后30分鐘之間的索結構所在環境的氣溫實測數據按照時間先后順序排列，找到索結構所在環境的氣溫實測數據序列中的最高溫度和最低溫度，用索結構所在環境的氣溫實測數據序列中的最高溫度減去最低溫度得到索結構所在環境的當日日出時刻到次日日出時刻后30分鐘之間的最大溫差，稱為環境最大溫差，記為Λ Traiax ;由索結構所在環境的氣溫實測數據序列通過常規數學計算得到索結構所在環境的氣溫關于時間的變化率，該變化率也隨著時間變化；通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的參考平板的向陽面的溫度的實測數據序列，參考平板的向陽面的溫度的實測數據序列由當日日出時刻到次日日出時刻后30分鐘之間的參考平板的向陽面的溫度的實測數據按照時間先后順序排列，找到參考平板的向陽面的溫度的實測數據序列中的最高溫度和最低溫度，用參考平板的向陽面的溫度的實測數據序列中的最高溫度減去最低溫度得到參考平板的向陽面的溫度的當日日出時刻到次日日出時刻后30分鐘之間的最大溫差，稱為參考平板最大溫差，記為ATpmax ;通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的所有R個索結構表面點的索結構表面溫度實測數據序列，有R個索結構表面點就有R個索結構表面溫度實測數據序列，每一個索結構表面溫度實測數據序列由一個索結構表面點的當日日出時刻到次日日出時刻后30分鐘之間的索結構表面溫度實測數據按照時間先后順序排列，找到每一個索結構表面溫度實測數據序列中的最高溫度和最低溫度，用每一個索結構表面溫度實測數據序列中的最高溫度減去最低溫度得到每一個索結構表面點的溫度的當日日出時刻到次日日出時刻后30分鐘之間的最大溫差，有R個索結構表面點就有R個當日日出時刻到次日日出時刻后30分鐘之間的最大溫差數值，其中的最大值稱為索結構表面最大溫差，記為ATsmax ;由每一索結構表面溫度實測數據序列通過常規數學計算得到每一個索結構表面點的溫度關于時間的變化率，每一個索結構表面點的溫度關于時間的變化率也隨著時間變化；通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的、在同一時刻、HBE個“索結構沿厚度的溫度分布數據”后，計算在每一個選取的海拔高度處共計BE個“相同海拔高度索結構沿厚度的溫度分布數據”中的最高溫度與最低溫度的差值，這個差值的絕對值稱為“相同海拔高度處索結構厚度方向最大溫差”，選取了 H個不同的海拔高度就有H個“相同海拔高度處索結構厚度方向最大溫差”，稱這H個“相同海拔高度處索結構厚度方向最大溫差”中的最大值為“索結構厚度方向最大溫差”，記為Λ Ttmax ；
[0017]b3:測量計算獲得索結構穩態溫度數據；首先，確定獲得索結構穩態溫度數據的時刻，與決定獲得索結構穩態溫度數據的時刻相關的條件有六項，第一項條件是獲得索結構穩態溫度數據的時刻介于當日日落時刻到次日日出時刻后30分鐘之間，日落時刻是指根據地球自轉和公轉規律確定的氣象學上的日落時刻，可以查詢資料或通過常規氣象學計算得到所需的每一日的日落時刻；第二項條件的a條件是在當日日出時刻到次日日出時刻后30分鐘之間的這段時間內，參考平板最大溫差Λ Tpmax和索結構表面最大溫差八1；_都不大于5攝氏度；第二項條件的b條件是在當日日出時刻到次日日出時刻后30分鐘之間的這段時間內，在前面測量計算得到的環境最大誤差不大于參考日溫差Λ?；，且參考平板最大溫差Λ Tpmax減去2攝氏度后不大于Δ Temax,且索結構表面最大溫差Δ 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除以Λ Th得到的數值，為方便敘述取Λ Th的單位為。C/m，為方便敘述取Ah的單位為m ;“R個索結構表面點”沿海拔高度的兩兩相鄰索結構表面點的定義是指只考慮海拔高度時，在“R個索結構表面點”中不存在一個索結構表面點，該索結構表面點的海拔高度數值介于兩兩相鄰索結構表面點的海拔高度數值之間；第三個條件是查詢或按氣象學常規計算得到索結構所在地和所在海拔區間的日照規律，再根據索結構的幾何特征及方位數據，在索結構上找到全年受日照時間最充分的那些表面點的位置，“R個索結構表面點”中至少有一個索結構表面點是索結構上全年受日照時間最充分的那些表面點中的一個占.[0018]c.按照“本方法的索結構的溫度測量計算方法”直接測量計算得到初始狀態下的索結構穩態溫度數據，初始狀態下的索結構穩態溫度數據稱為初始索結構穩態溫度數據，記為“初始索結構穩態溫度數據向量T。”;實測或查資料得到索結構所使用的各種材料的隨溫度變化的物理和力學性能參數；在實測得到T。的同時，也就是在獲得初始索結構穩態溫度數據向量T。的時刻的同一時刻，直接測量計算得到初始索結構的實測數據，初始索結構的實測數據是包括索結構集中載荷測量數據、索結構分布載荷測量數據、索結構體積載荷測量數據、索結構支座初始廣義位移測量數據、所有被監測量的初始數值、所有支承索的初始索力數據、初始索結構模態數據、初始索結構應變數據、初始索結構幾何數據、初始索結構支座廣義坐標數據、初始索結構角度數據、初始索結構空間坐標數據在內的實測數據，在得到初始索結構的實測數據的同時，測量計算得到包括支承索的無損檢測數據在內的能夠表達支承索的健康狀態的數據，此時的能夠表達支承索的健康狀態的數據稱為支承索初始健康狀態數據；所有被監測量的初始數值組成被監測量初始數值向量C。，被監測量初始數值向量C。的編號規則與M個被監測量的編號規則相同；利用能表達支承索初始健康狀態數據、索結構支座初始廣義位移測量數據以及索結構載荷測量數據建立被評估對象初始損傷向量d。，向量d。表示用初始力學計算基準模型A。表示的索結構的被評估對象的初始健康狀態；被評估對象初始損傷向量d。的元素個數等于N，d。的元素與被評估對象是一一對應關系，向量d。的元素的編號規則與被評估對象的編號規則相同；如果d。的某一個元素對應的被評估對象是索系統中的一根支承索，那么d。的該元素的數值代表對應支承索的初始損傷程度，若該元素的數值為O，表示該元素所對應的支承索是完好的，沒有損傷的，若其數值為100%，則表示該元素所對應的支承索已經完全喪失承載能力，若其數值介于O和100%之間，則表示該支承索喪失了相應比例的承載能力；如果d。的某一個元素對應的被評估對象是某一個支座的某一個廣義位移分量，那么d。的該元素的數值代表這個支座的該廣義位移分量的初始數值；如果d。的某一個元素對應的被評估對象是某一個載荷，本方法中取d。的該元素數值為0，代表這個載荷的變化的初始數值為O ;如果沒有索結構支座初始廣義位移測量數據或者可以認為索結構支座初始廣義位移為O時，向量d。中與索結構支座廣義位移相關的各元素數值取O ;如果沒有支承索的無損檢測數據及其他能夠表達支承索的健康狀態的數據時，或者可以認為結構初始狀態為無損傷無松弛狀態時，向量d。中與支承索相關的各元素數值取O ;初始索結構支座廣義坐標數據指索結構設計狀態下的支座廣義坐標數據，索結構支座初始廣義位移測量數據指在建立初始力學計算基準模型A。時，索結構支座相對于索結構設計狀態下的支座所發生的廣義位移；
[0019]d.根據索結構的設計圖、竣工圖和初始索結構的實測數據、支承索初始健康狀態數據、索結構支座初始廣義位移測量數據、索結構集中載荷測量數據、索結構分布載荷測量數據、索結構體積載荷測量數據、索結構所使用的各種材料的隨溫度變化的物理和力學性能參數、初始索結構穩態溫度數據向量T。和和前面步驟得到的所有的索結構數據，建立計入“索結構穩態溫度數據” 的索結構的初始力學計算基準模型A。，基于A。計算得到的索結構計算數據必須非常接近其實測數據，其間的差異不得大于5% ;對應于A。的“索結構穩態溫度數據”就是“初始索結構穩態溫度數據向量T。” ;對應于A。的被評估對象健康狀態用被評估對象初始損傷向量d。表示；對應于A。的所有被監測量的初始數值用被監測量初始數值向量C。表示；T。和d。是A。的參數，由A。的力學計算結果得到的所有被監測量的初始數值與C0表示的所有被監測量的初始數值相同，因此也可以說C。由A。的力學計算結果組成，在本方法中K、C。、d。和T。是不變的；
[0020]e.在本方法中，字母i除了明顯地表示步驟編號的地方外，字母i僅表示循環次數，即第i次循環；第i次循環開始時需要建立的或已建立的索結構的當前初始力學計算基準模型記為當前初始力學計算基準模型A^A。和Ai。計入了溫度參數，可以計算溫度變化對索結構的力學性能影響；第i次循環開始時，對應于Aitj的“索結構穩態溫度數據”用當前初始索結構穩態溫度數據向量Titj表示，向量Titj的定義方式與向量T。的定義方式相同，Τ1。的元素與T。的元素一一對應；第i次循環開始時需要的被評估對象當前初始損傷向量記為Cli0, Cli0表示該次循環開始時索結構Ai0的被評估對象的健康狀態，Cli0的定義方式與d。的定義方式相同，dt的元素與d。的元素一一對應；第i次循環開始時，所有被監測量的初始值，用被監測量當前初始數值向量C1。表示，向量C1。的定義方式與向量C。的定義方式相同，Ci0
的元素與C。的元素--對應，被監測量當前初始數值向量C1。表示對應于A1。的所有被監測
量的具體數值Ji。和d、是Ai0的特性參數，Ci0由Aitj的力學計算結果組成；第一次循環開始時，Ai0記為A1。，建立A1。的方法為使A1。等于A。；第一次循環開始時，Ti0記為T1。，建立T1。的方法為使T1。等于T。；第一次循環開始時，Cli0記為d1。，建立d1。的方法為使d1。等于d。；第一次循環開始時，Ci0記為C1。，建立C1。的方法為使C1。等于C。；
[0021]f.從這里進入由第f步到第q步的循環；在結構服役過程中，按照“本方法的索結構的溫度測量計算方法”不斷實測計算獲得索結構穩態溫度數據的當前數據，所有“索結構穩態溫度數據”的當前數據組成當前索結構穩態溫度數據向量Ti,向量Ti的定義方式與向量T。的定義方式相同，Ti的元素與T。的元素一一對應；在實測得到向量Ti的同時，實測得到在獲得當前索結構穩態溫度數據向量Ti的時刻的同一時刻的索結構中所有被監測量的當前值，所有這些數值組成被監測量當前數值向量Ci,向量Ci的定義方式與向量C。的定義方式相同，C1的兀素與C。的兀素對應，表不相同被監測量在不同時刻的數值；
[0022]g.根據當前索結構穩態溫度數據向量Ti,按照步驟gl至g3更新當前初始力學計算基準模型A1。、被監測量當前初始數值向量C1。和當前初始索結構穩態溫度數據向量f。，而被評估對象當前初始損傷向量f。保持不變；
[0023]gl.比較Ti與T:如果Ti等于f。，則A1。、Ci0和Titj保持不變；否則需要按下列步驟對Ai。和?進行更新；
[0024]g2.計算Ti與T。的差`，Ti與T。的差就是當前索結構穩態溫度數據關于初始索結構穩態溫度數據的變化，Ti與T。的差用穩態溫度變化向量S表示，S等于Ti減去T。，S表示索結構穩態溫度數據的變化；
[0025]g3.對A。中的索結構施加溫度變化，施加的溫度變化的數值就取自穩態溫度變化向量S，對A。中的索結構施加的溫度變化后得到更新的當前初始力學計算基準模型Ai。，更新Aitj的同時,Ti。所有元素數值也用Ti的所有元素數值對應代替，即更新了 Τ1。，這樣就得到了正確地對應于Aitj的Titj ;此時(保持不變；當更新Aitj后，Aitj的索的健康狀況用被評估對象當前初始損傷向量f。表示，Aitj的索結構穩態溫度用當前索結構穩態溫度數據向量Ti表示，更新C1。的方法是:當更新(后，通過力學計算得到Aitj中所有被監測量的、當前的具體數值，這些具體數值組成C1。；
[0026]h.在當前初始力學計算基準模型Aitj的基礎上，按照步驟hi至步驟h4進行若干次力學計算，通過計算建立單位損傷被監測量數值變化矩陣ACi和被評估對象單位變化向
量 Diu;
[0027]h1.在第i次循環開始時，直接按步驟h2至步驟h4所列方法獲得ACi和Diu;在其它時刻，當在步驟g中對Ai0進行更新后，必須按步驟h2至步驟h4所列方法重新獲得Λ Ci和Diu,如果在步驟g中沒有對Ai0進行更新，則在此處直接轉入步驟i進行后續工作；
[0028]h2.在當前初始力學計算基準模型Aitj的基礎上進行若干次力學計算，計算次數數值上等于所有被評估對象的數量N，有N個評估對象就有N次計算；依據被評估對象的編號規則，依次進行計算；每一次計算假設只有一個被評估對象在原有損傷或載荷的基礎上再增加單位損傷或載荷單位變化，具體的，如果該被評估對象是索系統中的一根支承索，那么就假設該支承索再增加單位損傷，如果該被評估對象是一個支座的一個方向的廣義位移分量，就假設該支座在該位移方向再增加單位廣義位移，如果該被評估對象是一個載荷，就假設該載荷再增加載荷單位變化，用Diuk記錄這一增加的單位損傷或單位廣義位移或載荷單位變化，其中k表示增加單位損傷或單位廣義位移或載荷單位變化的被評估對象的編號，Diuk是被評估對象單位變化向量Diu的一個元素，被評估對象單位變化向量Diu的元素的編號規則與向量d。的元素的編號規則相同；每一次計算中再增加單位損傷或單位廣義位移或載荷單位變化的被評估對象不同于其它次計算中再增加單位損傷或單位廣義位移或載荷單位變化的被評估對象，每一次計算都利用力學方法計算索結構的所有被監測量的當前計算值，每一次計算得到的所有被監測量的當前計算值組成一個被監測量計算當前向量；當假設第k個被評估對象再增加單位損傷或單位廣義位移或載荷單位變化時，用Citk表示對應的“被監測量計算當前向量”;在本步驟中給各向量的元素編號時，應同本方法中其它向量使用同一編號規則，以保證本步驟中各向量中的任意一個元素，同其它向量中的、編號相同的兀素，表達了同一被監測量或同一對象的相關信息A1 tk的定義方式與向量C。的定義方式相同，C1tk的兀素與C。的兀素對應；
[0029]h3.每一次計算得到的向量Citk減去向量C1。得到一個向量，再將該向量的每一個元素都除以本次計算所假設的單位損傷或單位廣義位移或載荷單位變化數值后得到一個“被監測量的數值變化向量δ CiJ;有N個被評估對象就有N個“被監測量的數值變化向量”；
[0030]h4.由這N個“被監測量的數值變化向量”按照N個被評估對象的編號規則，依次組成有N列的“單位損傷被監測量數值變化矩陣單位損傷被監測量數值變化矩陣ACi的每一列對應于一個被監測量單位變化向量；單位損傷被監測量數值變化矩陣ACi的每一行對應于同一個被監測量在不同被評估對象增加單位損傷或單位廣義位移或載荷單位變化時的不同的單位變化幅度；單位損傷被監測量數值變化矩陣ACi的列的編號規則與向量
d。的元素的編號規則相同，單位損傷被監測量數值變化矩陣ACi的行的編號規則與M個被監測量的編號規則相同；
[0031]1.定義當前名義損傷向量d、和當前實際損傷向量ClSdic^Pdi的元素個數等于被評估對象的數量，Clic和Cli的元素和被評估對象之間是一一對應關系，f。的元素數值代表對應被評估對象的名義損傷程度或名義廣義位移或名義載荷變化量，Clic和Cli與被評估對象初始損傷向量d。的元素編號規則相同，Clic的元素、Cli的元素與d。的元素是一一對應關系；
[0032]j.依據被監測量當前數值向量Ci同“被監測量當前初始數值向量 <”、“單位損傷被監測量數值變化矩陣AC1”和“當前名義損傷向量屮。”間存在的近似線性關系，該近似線性關系可表達為式1，式I中除(外的其它量均為已知，求解式I就可以算出當前名義損傷向量d1。；
[0033]Ci =C0+AC md'c式 I
[0034]k.利用式2表達的當前實際損傷向量Cli的第k個元素Clik同被評估對象當前初始損傷向量必的第k個元素d、和當前名義損傷向量 < 的第k個元素Clidt間的關系，計算得到當前實際損傷向量Cli的所有元素；
[0035]
【權利要求】
1.應變監測受損索載荷廣義位移遞進式識別方法，其特征在于所述方法包括: a.為敘述方便起見，本方法統一稱被評估的支座廣義位移分量、支承索和載荷為被評估對象，設被評估的支座廣義位移分量的數量、支承索的數量和載荷的數量之和為N，即被評估對象的數量為N ;確定被評估對象的編號規則，按此規則將索結構中所有的被評估對象編號，該編號在后續步驟中將用于生成向量和矩陣；本方法用變量k表示這一編號，k=l, 2，3，…，N ;確定指定的被監測點，被監測點即表征索結構應變信息的所有指定點，并給所有指定點編號；確定被監測點的被監測的應變方向，并給所有指定的被監測應變編號，“被監測應變編號”在后續步驟中將用于生成向量和矩陣，“索結構的全部被監測的應變數據”由上述所有被監測應變組成；本方法將“索結構的被監測的應變數據”簡稱為“被監測量”;所有被監測量的數量之和記為M，M不得小于N ;本方法中對同一個量實時監測的任何兩次測量之間的時間間隔不得大于30分鐘，測量記錄數據的時刻稱為實際記錄數據時刻；物體、結構承受的外力可稱為載荷，載荷包括面載荷和體積載荷；面載荷又稱表面載荷，是作用于物體表面的載荷，包括集中載荷和分布載荷兩種；體積載荷是連續分布于物體內部各點的載荷，包括物體的自重和慣性力在內；集中載荷分為集中力和集中力偶兩種，在包括笛卡爾直角坐標系在內的坐標系中，一個集中力可以分解成三個分量，同樣的，一個集中力偶也可以分解成三個分量，如果載荷實際上是集中載荷，在本方法中將一個集中力分量或一個集中力偶分量計為或統計為一個載荷，此時載荷的變化具體化為一個集中力分量或一個集中力偶分量的變化；分布載荷分為線分布載荷和面分布載荷，分布載荷的描述至少包括分布載荷的作用區域和分布載荷的大小，分布載荷的大小用分布集度來表達，分布集度用分布特征和幅值來表達；如果載荷實際上是分布載荷，本方法談論載荷的變化時，實際上是指分布載荷分布集度的幅值的改變，而所有分布載荷的作用區域和分布集度的分布特征是不變的；在包括笛卡爾直角坐標系在內的坐標系中，一個分布載荷可以分解成三個分量，如果這分布載荷的三個分量的各自的分布集度的幅值發生變化，且變化的比率不全部相同，那么在本方法中把這分布載荷的三個分量計為或統計為三個分布載荷，此時一個載荷就代表分布載荷的一個分量；體積載荷是連續分布于物體內部各點的載荷，體積載荷的描述至少包括體積載荷的作用區域和體積載荷的大小，體積載荷的大小用分布集度來表達，分布集度用分布特征和幅值來表達；如果載荷實際上是體積載荷，在本方法中實際處理的是體積載荷分布集度的幅值的·改變，而所有體積載荷的作用區域和分布集度的分布特征是不變的，此時在本方法中提到載荷的改變時實際上是指體積載荷的分布集度的幅值的改變，此時，發生變化的載荷是指那些分布集度的幅值發生變化的體積載荷；在包括笛卡爾直角坐標系在內的坐標系中，一個體積載荷可以分解成三個分量，如果這體積載荷的三個分量的各自的分布集度的幅值發生變化，且變化的比率不全部相同，那么在本方法中把這體積載荷的三個分量計為或統計為三個分布載荷； b.本方法定義“本方法的索結構的溫度測量計算方法”按步驟bl至b3進行； bl:查詢或實測得到索結構組成材料及索結構所處環境的隨溫度變化的傳熱學參數，利用索結構的設計圖、竣工圖和索結構的幾何實測數據，利用這些數據和參數建立索結構的傳熱學計算模型；查詢索結構所在地不少于2年的近年來的氣象資料，統計得到這段時間內的陰天數量記為T個陰天，在本方法中將白天不能見到太陽的一整日稱為陰天，統計得到T個陰天中每一個陰天的O時至次日日出時刻后30分鐘之間的最高氣溫與最低氣溫，日出時刻是指根據地球自轉和公轉規律確定的氣象學上的日出時刻，不表示當天一定可以看見太陽，可以查詢資料或通過常規氣象學計算得到所需的每一日的日出時刻，每一個陰天的O時至次日日出時刻后30分鐘之間的最高氣溫減去最低氣溫稱為該陰天的日氣溫的最大溫差，有T個陰天，就有T個陰天的日氣溫的最大溫差，取T個陰天的日氣溫的最大溫差中的最大值為參考日溫差，參考日溫差記為△ I；;查詢索結構所在地和所在海拔區間不少于2年的近年來的氣象資料或實測得到索結構所處環境的溫度隨時間和海拔高度的變化數據和變化規律，計算得到索結構所在地和所在海拔區間不少于2年的近年來的索結構所處環境的溫度關于海拔高度的最大變化率ATh,為方便敘述取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分鐘之間的最大溫差，稱為環境最大溫差，記為ATemax ；由索結構所在環境的氣溫實測數據序列通過常規數學計算得到索結構所在環境的氣溫關于時間的變化率，該變化率也隨著時間變化；通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的參考平板的向陽面的溫度的實測數據序列，參考平板的向陽面的溫度的實測數據序列由當日日出時刻到次日日出時刻后30分鐘之間的參考平板的向陽面的溫度的實測數據按照時間先后順序排列，找到參考平板的向陽面的溫度的實測數據序列中的最高溫度和最低溫度，用參考平板的向陽面的溫度的實測數據序列中的最高溫度減去最低溫度得到參考平板的向陽面的溫度的當日日出時刻到次日日出時刻后30分鐘之間的最大溫差，稱為參考平板最大溫差，記為ATpmax ;通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的所有R個索結構表面點的索結構表面溫度實測數據序列，有R個索結構表面點就有R個索結構表面溫度實測數據序列，每一個索結構表面溫度實測數據序列由一個索結構表面點的當日日出時刻到次日日出時刻后30分鐘之間的索結構表面溫度實測數據按照時間先后順序排列，找到每一個索結構表面溫度實測數據序列中的最高溫度和最低溫度，用每一個索結構表面溫度實測數據序列中的最高溫度減去最低溫度得到每一個索結構表面點的溫度的當日日出時刻到次日日出時刻后30分鐘之間的最大溫差，有R個索結構表面點就有R個當日日出時刻到次日日出時刻后30分鐘之間的最大溫差數值，其中的最大值稱為索結構表面最大溫差，記為;由每一索結構表面溫度實測數據序列通過常規數學計算得到每一個索結構表面點的溫度關于時間的變化率，每一個索結構表面點的溫度關于時間的變化率也隨著時間變化；通過實時監測得到當日日出時刻到次日日出時刻后30分鐘之間的、在同一時刻、HBE個“索結構沿厚度的溫度分布數據”后，計算在每一個選取的海拔高度處共計BE個“相同海拔高度索結構沿厚度的溫度分布數據”中的最高溫度與最低溫度的差值，這個差值的絕對值稱為“相同海拔高度處索結構厚度方向最大溫差”，選取了 H個不同的海拔高度就有H個“相同`海拔高度處索結構厚度方向最大溫差”，稱這H個“相同海拔高度處索結構厚度方向最大溫差”中的最大值為“索結構厚度方向最大溫差”，記為Δ Ttmax ; b3:測量計算獲得索結構穩態溫度數據；首先，確定獲得索結構穩態溫度數據的時刻，與決定獲得索結構穩態溫度數據的時刻相關的條件有六項，第一項條件是獲得索結構穩態溫度數據的時刻介于當日日落時刻到次日日出時刻后30分鐘之間，日落時刻是指根據地球自轉和公轉規律確定的氣象學上的日落時刻，可以查詢資料或通過常規氣象學計算得到所需的每一日的日落時刻；第二項條件的a條件是在當日日出時刻到次日日出時刻后30分鐘之間的這段時間內，參考平板最大溫差Λ Tpmax和索結構表面最大溫差ATsmax都不大于5攝氏度；第二項條件的b條件是在當日日出時刻到次日日出時刻后30分鐘之間的這段時間內，在前面測量計算得到的環境最大誤差Λ Traiax不大于參考日溫差Λ ?；，且參考平板最大溫差ATpmax減去2攝氏度后不大于Λ Temax，且索結構表面最大溫差ATsmax不大于ATpmax;只需滿足第二項的a條件和b條件中的一項就稱為滿足第二項條件；第三項條件是在獲得索結構穩態溫度數據的時刻，索結構所在環境的氣溫關于時間的變化率的絕對值不大于每小時0.1攝氏度；第四項條件是在獲得索結構穩態溫度數據的時刻，R個索結構表面點中的每一個索結構表面點的溫度關于時間的變化率的絕對值不大于每小時0.1攝氏度；第五項條件是在獲得索結構穩態溫度數據的時刻，R個索結構表面點中的每一個索結構表面點的索結構表面溫度實測數據為當日日出時刻到次日日出時刻后30分鐘之間的極小值；第六項條件是在獲得索結構穩態溫度數據的時刻，“索結構厚度方向最大溫差” Λ Ttmax不大于I攝氏度；本方法利用上述六項條件，將下列三種時刻中的任意一種稱為“獲得索結構穩態溫度數據的數學時刻”，第一種時刻是滿足上述“與決定獲得索結構穩態溫度數據的時刻相關的條件”中的第一項至第五項條件的時刻，第二種時刻是僅僅滿足上述“與決定獲得索結構穩態溫度數據的時刻相關的條件”中的第六項條件的時刻，第三種時刻是同時滿足上述“與決定獲得索結構穩態溫度數據的時刻相關的條件”中的第一項至第六項條件的時刻；當獲得索結構穩態溫度數據的數學時刻就是本方法中實際記錄數據時刻中的一個時，獲得索結構穩態溫度數據的時刻就是獲得索結構穩態溫度數據的數學時刻；如果獲得索結構穩態溫度數據的數學時刻不是本方法中實際記錄數據時刻中的任一個時刻，則取本方法最接近于獲得索結構穩態溫度數據的數學時刻的那個實際記錄數據的時刻為獲得索結構穩態溫度數據的時刻；本方法將使用在獲得索結構穩態溫度數據的時刻測量記錄的量進行索結構相關健康監測分析；本方法近似認為獲得索結構穩態溫度數據的時刻的索結構溫度場處于穩態，即此時刻的索結構溫度不隨時間變化，此時刻就是本方法的“獲得索結構穩態溫度數據的時刻”;然后，根據索結構傳熱特性，利用獲得索結構穩態溫度數據的時刻的“R個索結構表面溫度實測數據”和“ΗΒΕ個索結構沿厚度溫度實測數據”，利用索結構的傳熱學計算模型，通過常規傳熱計算得到在獲得索結構穩態溫度數據的時刻的索結構的溫度分布，此時索結構的溫度場按穩態進行計算，計算得到的在獲得索結構穩態溫度數據的時刻的索結構的溫度分布數據包括索結構上R個索結構表面點的計算溫度，R個索結構表面點的計算溫度稱為R個索結構穩態表面溫度計算數據，還包括索結構在前面選定的HBE個“測量索結構沿厚度的溫度分布數據的點”的計算溫度，HBE個“測量索結構沿厚度的溫度分布數據的點”的計算溫度稱為“ΗΒΕ個索結構沿厚度溫度計算數據”，當R個索結構表面溫度實測數據與R個索結構穩態表面溫度計算數據對應相等時，且“ΗΒΕ個索結構沿厚度溫度實測數據”與“ΗΒΕ個索結構沿厚度溫度計·算數據”對應相等時，計算得到的在獲得索結構穩態溫度數據的時刻的索結構的溫度分布數據在本方法中稱為“索結構穩態溫度數據”，此時的“R個索結構表面溫度實測數據”稱為“R個索結構穩態表面溫度實測數據”，“ΗΒΕ個索結構沿厚度溫度實測數據”稱為“ΗΒΕ個索結構沿厚度穩態溫度實測數據”;在索結構的表面上取“R個索結構表面點”時，“R個索結構表面點”的數量與分布必須滿足三個條件，第一個條件是當索結構溫度場處于穩態時，當索結構表面上任意一點的溫度是通過“R個索結構表面點”中與索結構表面上該任意點相鄰的點的實測溫度線性插值得到時，線性插值得到的索結構表面上該任意點的溫度與索結構表面上該任意點的實際溫度的誤差不大于5% ;索結構表面包括支承索表面；第二個條件是“R個索結構表面點”中在同一海拔高度的點的數量不小于4，且“R個索結構表面點”中在同一海拔高度的點沿著索結構表面均布；“R個索結構表面點”沿海拔高度的所有兩兩相鄰索結構表面點的海拔高度之差的絕對值中的最大值Λ h不大于0.2°C除以Λ Th得到的數值，為方便敘述取Λ Th的單位為。C/m，為方便敘述取Ah的單位為m ;“R個索結構表面點”沿海拔高度的兩兩相鄰索結構表面點的定義是指只考慮海拔高度時，在“R個索結構表面點”中不存在一個索結構表面點，該索結構表面點的海拔高度數值介于兩兩相鄰索結構表面點的海拔高度數值之間；第三個條件是查詢或按氣象學常規計算得到索結構所在地和所在海拔區間的日照規律，再根據索結構的幾何特征及方位數據，在索結構上找到全年受日照時間最充分的那些表面點的位置，“R個索結構表面點”中至少有一個索結構表面點是索結構上全年受日照時間最充分的那些表面點中的一個點； c.按照“本方法的索結構的溫度測量計算方法”直接測量計算得到初始狀態下的索結構穩態溫度數據，初始狀態下的索結構穩態溫度數據稱為初始索結構穩態溫度數據，記為“初始索結構穩態溫度數據向量T?！?；實測或查資料得到索結構所使用的各種材料的隨溫度變化的物理和力學性能參數；在實測得到T。的同時，也就是在獲得初始索結構穩態溫度數據向量T。的時刻的同一時刻，直接測量計算得到初始索結構的實測數據，初始索結構的實測數據是包括索結構集中載荷測量數據、索結構分布載荷測量數據、索結構體積載荷測量數據、索結構支座初始廣義位移測量數據、所有被監測量的初始數值、所有支承索的初始索力數據、初始索結構模態數據、初始索結構應變數據、初始索結構幾何數據、初始索結構支座廣義坐標數據、初始索結構角度數據、初始索結構空間坐標數據在內的實測數據，在得到初始索結構的實測數據的同時，測量計算得到包括支承索的無損檢測數據在內的能夠表達支承索的健康狀態的數據，此時的能夠表達支承索的健康狀態的數據稱為支承索初始健康狀態數據；所有被監測量的初始數值組成被監測量初始數值向量C。，被監測量初始數值向量C。的編號規則與M個被監測量的編號規則相同；利用能表達支承索初始健康狀態數據、索結構支座初始廣義位移測量數據以及索結構載荷測量數據建立被評估對象初始損傷向量d。，向量d。表示用初始力學計算基準模型A。表示的索結構的被評估對象的初始健康狀態；被評估對象初始損傷向量d。的元素個數等于N，d。的元素與被評估對象是一一對應關系，向量d。的元素的編號規則與被評估對象的編號規則相同；如果d。的某一個元素對應的被評估對象是索系統中的一根支承索，那么d。的該元素的數值代表對應支承索的初始損傷程度，若該元素的數值為O，表示該元素所對應的支承索是完好的，沒有損傷的，若其數值為100%，則表示該元素所對應的支承索已經完全喪失承載能力，若其數值介于O和100%之間，則表示該支承索喪失了相應比例的承載能力；如果d。的某一個元素對應的被評估對象是某一個支座的某一個廣義位移分量`，那么d。的該元素的數值代表這個支座的該廣義位移分量的初始數值；如果d。的某一個元素對應的被評估對象是某一個載荷，本方法中取d。的該元素數值為0，代表這個載荷的變化的初始數值為O ;如果沒有索結構支座初始廣義位移測量數據或者可以認為索結構支座初始廣義位移為O時，向量d。中與索結構支座廣義位移相關的各元素數值取O ;如果沒有支承索的無損檢測數據及其他能夠表達支承索的健康狀態的數據時，或者可以認為結構初始狀態為無損傷無松弛狀態時，向量d。中與支承索相關的各元素數值取O ;初始索結構支座廣義坐標數據指索結構設計狀態下的支座廣義坐標數據，索結構支座初始廣義位移測量數據指在建立初始力學計算基準模型K時，索結構支座相對于索結構設計狀態下的支座所發生的廣義位移； d.根據索結構的設計圖、竣工圖和初始索結構的實測數據、支承索初始健康狀態數據、索結構支座初始廣義位移測量數據、索結構集中載荷測量數據、索結構分布載荷測量數據、索結構體積載荷測量數據、索結構所使用的各種材料的隨溫度變化的物理和力學性能參數、初始索結構穩態溫度數據向量T。和和前面步驟得到的所有的索結構數據，建立計入“索結構穩態溫度數據”的索結構的初始力學計算基準模型A。，基于A。計算得到的索結構計算數據必須非常接近其實測數據，其間的差異不得大于5% ;對應于A。的“索結構穩態溫度數據”就是“初始索結構穩態溫度數據向量T。” ;對應于A。的被評估對象健康狀態用被評估對象初始損傷向量d。表示；對應于A。的所有被監測量的初始數值用被監測量初始數值向量C。表示；T。和d。是A。的參數，由A。的力學計算結果得到的所有被監測量的初始數值與C。表示的所有被監測量的初始數值相同，因此也可以說C。由A。的力學計算結果組成，在本方法中A0、C0、d0和T0是不變的； e.在本方法中，字母i除了明顯地表示步驟編號的地方外，字母i僅表示循環次數，即第i次循環；第i次循環開始時需要建立的或已建立的索結構的當前初始力學計算基準模型記為當前初始力學計算基準模型A^A。和Ai。計入了溫度參數，可以計算溫度變化對索結構的力學性能影響；第i次循環開始時，對應于Aitj的“索結構穩態溫度數據”用當前初始索結構穩態溫度數據向量Titj表示，向量Titj的定義方式與向量T。的定義方式相同，Titj的元素與T。的元素一一對應；第i次循環開始時需要的被評估對象當前初始損傷向量記為(^。，屮。表示該次循環開始時索結構Ai0的被評估對象的健康狀態，Cli0的定義方式與d。的定義方式相同，dt的元素與d。的元素一一對應；第i次循環開始時，所有被監測量的初始值，用被監測量當前初始數值向量C1。表示，向量C1。的定義方式與向量C。的定義方式相同，C1。的元素與C。的元素一一對應，被監測量當前初始數值向量Ci。表示對應于Aitj的所有被監測量的具體數值；1"。和f。是Aitj的特性參數，C1。由Ai0的力學計算結果組成；第一次循環開始時，Aitj記為A1。，建立A1。的方法為使A1。等于A。；第一次循環開始時，Titj記為T1。，建立T1。的方法為使T1。等于T。；第一次循環開始時，Cli0記為d1。，建立d1。的方法為使d1。等于d。；第一次循環開始時，Ci0記為C10,建立C1。的方法為使C1。等于C0 ； f.從這里進入由第f步到第q步的循環；在結構服役過程中，按照“本方法的索結構的溫度測量計算方法”不斷實測計`算獲得索結構穩態溫度數據的當前數據，所有“索結構穩態溫度數據”的當前數據組成當前索結構穩態溫度數據向量Ti,向量Ti的定義方式與向量T。的定義方式相同，Ti的元素與T。的元素一一對應；在實測得到向量Ti的同時，實測得到在獲得當前索結構穩態溫度數據向量Ti的時刻的同一時刻的索結構中所有被監測量的當前值，所有這些數值組成被監測量當前數值向量Ci,向量Ci的定義方式與向量C。的定義方式相同，C1的兀素與C。的兀素對應，表不相同被監測量在不同時刻的數值； g.根據當前索結構穩態溫度數據向量T1,按照步驟gl至g3更新當前初始力學計算基準模型A1。、被監測量當前初始數值向量C1。和當前初始索結構穩態溫度數據向量Τ1。，而被評估對象當前初始損傷向量f。保持不變； gl.比較Ti與f。，如果Ti等于f。，則A1。、Ci0和Titj保持不變；否則需要按下列步驟對Ai0和Ti0進行更新； g2.計算Ti與T。的差，Ti與T。的差就是當前索結構穩態溫度數據關于初始索結構穩態溫度數據的變化，Ti與T。的差用穩態溫度變化向量S表示，S等于Ti減去T。，S表示索結構穩態溫度數據的變化； g3.對A。中的索結構施加溫度變化，施加的溫度變化的數值就取自穩態溫度變化向量S，對A。中的索結構施加的溫度變化后得到更新的當前初始力學計算基準模型A1。，更新A1。的同時，Titj所有元素數值也用Ti的所有元素數值對應代替，即更新了 Τ1。，這樣就得到了正確地對應于Aitj的Titj ;此時(保持不變；當更新Aitj后，Aitj的索的健康狀況用被評估對象當前初始損傷向量(表示，Aitj的索結構穩態溫度用當前索結構穩態溫度數據向量Ti表示，更新C1。的方法是:當更新(后，通過力學計算得到Aitj中所有被監測量的、當前的具體數值，這些具體數值組成C、；h.在當前初始力學計算基準模型Ai。的基礎上，按照步驟hi至步驟h4進行若干次力學計算，通過計算建立單位損傷被監測量數值變化矩陣Λ Ci和被評估對象單位變化向量Diu ；h1.在第i次循環開始時，直接按步驟h2至步驟h4所列方法獲得Λ Ci和Diu ;在其它時刻，當在步驟g中對Ai0進行更新后，必須按步驟h2至步驟h4所列方法重新獲得Λ Ci和Diu,如果在步驟g中沒有對Aitj進行更新，則在此處直接轉入步驟i進行后續工作； h2.在當前初始力學計算基準模型Aitj的基礎上進行若干次力學計算，計算次數數值上等于所有被評估對象的數量N，有N個評估對象就有N次計算；依據被評估對象的編號規則，依次進行計算；每一次計算假設只有一個被評估對象在原有損傷或載荷的基礎上再增加單位損傷或載荷單位變化，具體的，如果該被評估對象是索系統中的一根支承索，那么就假設該支承索再增加單位損傷，如果該被評估對象是一個支座的一個方向的廣義位移分量，就假設該支座在該位移方向再增加單位廣義位移，如果該被評估對象是一個載荷，就假設該載荷再增加載荷單位變化，用Diuk記錄這一增加的單位損傷或單位廣義位移或載荷單位變化，其中k表示增加單位損傷或單位廣義位移或載荷單位變化的被評估對象的編號，Diuk是被評估對象單位變化向量Diu的一個元素，被評估對象單位變化向量Diu的元素的編號規則與向量d。的元素的編號規則相同；每一次計算中再增加單位損傷或單位廣義位移或載荷單位變化的被評估對象不同于其它次計算中再增加單位損傷或單位廣義位移或載荷單位變化的被評估對象，每一次計算都利用力學方法計算索結構的所有被監測量的當前計算值，每一次計算得到的所有被監測量的當前計算值組成一個被監測量計算當前向量；當假設第k個被評估對象再增加單位損傷或單位廣義位移或載荷單位變化時，用Citk表示對應的“被監測量計算當前向量”;在本步`驟中給各向量的元素編號時，應同本方法中其它向量使用同一編號規則，以保證本步驟中各向量中的任意一個元素，同其它向量中的、編號相同的兀素，表達了同一被監測量或同一對象的相關信息A1 tk的定義方式與向量C。的定義方式相同，C1tk的兀素與C。的兀素對應； h3.每一次計算得到的向量Citk減去向量C1。得到一個向量，再將該向量的每一個元素都除以本次計算所假設的單位損傷或單位廣義位移或載荷單位變化數值后得到一個“被監測量的數值變化向量δ CY ;有N個被評估對象就有N個“被監測量的數值變化向量”；h4.由這N個“被監測量的數值變化向量”按照N個被評估對象的編號規則，依次組成有N列的“單位損傷被監測量數值變化矩陣△ Ci”;單位損傷被監測量數值變化矩陣八^勺每一列對應于一個被監測量單位變化向量；單位損傷被監測量數值變化矩陣ACi的每一行對應于同一個被監測量在不同被評估對象增加單位損傷或單位廣義位移或載荷單位變化時的不同的單位變化幅度；單位損傷被監測量數值變化矩陣ACi的列的編號規則與向量d。的元素的編號規則相同，單位損傷被監測量數值變化矩陣ACi的行的編號規則與M個被監測量的編號規則相同；i.定義當前名義損傷向量(和當前實際損傷向量dSdt和Cli的元素個數等于被評估對象的數量，Clic和Cli的元素和被評估對象之間是--對應關系，f。的元素數值代表對應被評估對象的名義損傷程度或名義廣義位移或名義載荷變化量，Cli。和Cli與被評估對象初始損傷向量d。的元素編號規則相同，Clic的元素、Cli的元素與d。的元素是一一對應關系；j.依據被監測量當前數值向量Ci同“被監測量當前初始數值向量C1?！?、“單位損傷被監測量數值變化矩陣AC1”和“當前名義損傷向量屮?！遍g存在的近似線性關系，該近似線性關系可表達為式1，式I中除f。外的其它量均為已知，求解式I就可以算出當前名義損傷向Sdi,;
【文檔編號】G01L5/00GK103852287SQ201410085857
【公開日】2014年6月11日 申請日期:2014年3月10日 優先權日:2014年3月10日
【發明者】韓玉林, 韓佳邑 申請人:東南大學
專業數控銑床產品供應商
20年生產經驗，實力工廠，精工品質，質量有保障
