[發明專利]有支座沉降時空間坐標監測識別松弛的支承索的遞進式方法

有效

• 申請號: 201010127947.0
• 申請日: 2010-03-17
• 公開（公告）號: CN101782946A
• 公開（公告）日: 2010-07-21
• 發明（設計）人: 韓玉林
• 申請（專利權）人: 東南大學
• 主分類號: G06F19/00
• 分類號: G06F19/00;G01N19/00;G01L5/04;G01B21/02
• 代理公司: 南京蘇高專利商標事務所(普通合伙) 32204
• 代理人: 柏尚春

1.一種有支座沉降時空間坐標監測識別松弛的支承索的遞進式方法，其特征是該方法包括：

a.設共有N根索，首先確定索的編號規則，按此規則將索結構中所有的索編號，該編號在后續步驟中將用于生成向量和矩陣；

b.確定指定的將被監測空間坐標的被測量點，給所有指定點編號；確定過每一測量點的將被監測的空間坐標分量，給所有被測量空間坐標分量編號；上述編號在后續步驟中將用于生成向量和矩陣；“結構的全部被監測的空間坐標數據”由上述所有被測量空間坐標分量組成；為方便起見，將“結構的被監測的空間坐標數據”稱為“被監測量”；測量點的數量不得小于索的數量；所有被測量空間坐標分量的數量之和不得小于索的數量；

c.利用包括索的無損檢測數據的能夠表達索的健康狀態的數據建立初始虛擬損傷向量dⁱ_o，其中i表示循環次數，后面i及上標i都表示循環次數，i＝1，2，3，......；第一次循環時dⁱ_o記為d¹_o；如果沒有索的無損檢測數據及其他能夠表達索的健康狀態的數據時，或者可以認為結構初始狀態為無松弛、無損傷狀態時，向量d¹_o的各元素數值取0；

d.在建立初始虛擬損傷向量d¹_o的同時，直接測量計算得到索結構的所有被監測量的初始數值，組成被監測量的初始數值向量C¹_o；

e.在建立初始虛擬損傷向量d¹_o和被監測量的初始數值向量C¹_o的同時，直接測量計算得到所有支承索的初始索力，組成初始索力向量F_o；同時，依據結構設計數據、竣工數據得到所有支承索的初始自由長度，組成初始自由長度向量l_o；同時，依據結構設計數據、竣工數據或實測得到索結構的初始幾何數據；同時，實測或根據結構設計、竣工資料得到所有索的彈性模量、密度、初始橫截面面積；

f.建立索結構的初始力學計算基準模型A_o，建立初始索結構支座坐標向量U_o，建立第一次循環開始時需要的索結構的力學計算基準模型A¹；依據索結構竣工之時的索結構的實測數據，該實測數據包括索結構形狀數據、索力數據、拉桿拉力數據、索結構支座坐標數據、索結構模態數據、所有索的彈性模量、密度、初始橫截面面積，以及包括索的無損檢測數據的能夠表達索的健康狀態的數據，依據設計圖和竣工圖，利用力學方法建立索結構的初始力學計算基準模型A_o；如果沒有索結構竣工之時的結構的實測數據，那么就在建立健康監測系統前對該索結構進行實測，同樣得到索結構的實測數據，根據此數據和索結構的設計圖、竣工圖，同樣利用力學方法建立索結構的初始力學計算基準模型A_o；不論用何種方法獲得A_o，基于A_o計算得到的索結構計算數據必須非常接近其實測數據，其間的差異不得大于5％；對應于A_o的索結構支座坐標數據組成初始索結構支座坐標向量U_o；A_o和U_o是不變的，只在第一次循環開始時建立；第i次循環開始時建立的索結構的力學計算基準模型記為Aⁱ，其中i表示循環次數；其中字母i除了明顯地表示步驟編號的地方外，字母i僅表示循環次數，即第i次循環；因此第一次循環開始時建立的索結構的力學計算基準模型記為A¹，其中A¹就等于A_o；為敘述方便，命名“索結構當前力學計算基準模型A^ti_o”，在每一次循環中A^ti_o根據需要會不斷更新，每一次循環開始時，A^ti_o等于Aⁱ；同樣為敘述方便，命名“索結構實測支座坐標向量U^ti”，在每一次循環中，不斷實測獲得索結構支座坐標當前數據，所有索結構支座坐標當前數據組成當前索結構實測支座坐標向量U^ti，向量U^ti的元素與向量U_o相同位置的元素表示相同支座的相同方向的坐標；為敘述方便起見，對于第i次循環，將上一次更新A^ti_o時的索結構支座坐標當前數據記為當前索結構支座坐標向量U^ti_o；第一次循環開始時，A^t1_o等于A¹，U^t1_o等于U_o；A¹對應的索的健康狀態由d¹_o描述；力學計算基準模型Aⁱ對應的索的健康狀態由dⁱ_o描述；

g.每一次循環開始時，令A^ti_o等于Aⁱ；實測獲得索結構支座坐標當前數據，所有索結構支座坐標當前數據組成當前索結構實測支座坐標向量U^ti，根據當前索結構實測支座坐標向量U^ti，按步驟g1、g2更新索結構當前力學計算基準模型A^ti_o和當前索結構支座坐標向量U^ti_o；

g1.實測得到當前索結構實測支座坐標向量U^ti后，比較U^ti和U^ti_o，如果U^ti等于U^ti_o，則不需要對A^ti_o進行更新；

g2.實測得到當前索結構實測支座坐標向量U^ti后，比較U^ti和U^ti_o，如果U^ti不等于U^ti_o，則需要對A^ti_o進行更新，更新方法是：先計算U^ti與U_o的差，U^ti與U_o的差就是當前索結構支座關于在建立A_o時的索結構支座的當前支座位移，用當前支座位移向量V表示支座位移，當前支座位移向量V中的元素與支座位移分量之間是一一對應關系，當前支座位移向量V中一個元素的數值對應于一個指定支座的一個指定方向的位移，其中支座位移在重力方向的分量就是支座沉降量；更新A^ti_o的方法是：在A_o的基礎上令索的健康狀況為索系統初始損傷向量dⁱ_o，再進一步對A_o中的索結構支座施加當前支座位移約束，當前支座位移約束的數值就取自當前支座位移向量V中對應元素的數值，對A_o中的索結構支座施加當前支座位移約束后，最終得到的就是更新的當前力學計算基準模型A^ti_o，更新A^ti_o的同時，U^ti_o所有元素數值也用U^ti所有元素數值代替，即更新了U^ti_o，這樣就得到了正確地對應于A^ti_o的U^ti_o；

h.在索結構當前力學計算基準模型A^ti_o的基礎上按步驟h1、h2、h3、h4進行力學計算，通過計算獲得索結構虛擬單位損傷被監測量數值變化矩陣ΔCⁱ和名義虛擬單位損傷向量Dⁱ_u；

h1.在第i次循環開始時，直接按步驟h2至步驟h4所列方法獲得索結構虛擬單位損傷被監測量變化矩陣ΔCⁱ和名義虛擬單位損傷向量Dⁱ_u；在非第i次循環開始的時刻，當步驟g中對A^ti_o進行更新后，直接按步驟h2至步驟h4所列方法獲得索結構虛擬單位損傷被監測量變化矩陣ΔCⁱ和名義虛擬單位損傷向量Dⁱ_u；在非第i次循環開始的時刻，如果在步驟g中沒有對A^ti_o進行更新，則在此處直接轉入步驟i進行后續工作；

h2.在索結構當前力學計算基準模型A^ti_o的基礎上進行若干次力學計算，計算次數數值上等于所有索的數量，有N根索就有N次計算，每一次計算假設索系統中只有一根索在原有虛擬損傷的基礎上再增加虛擬單位損傷，每一次計算中出現虛擬單位損傷的索不同于其它次計算中出現虛擬單位損傷的索，并且每一次假定有虛擬單位損傷的索的虛擬單位損傷值可以不同于其他索的虛擬單位損傷值，用“名義虛擬單位損傷向量Dⁱ_u”記錄所有索的假定的單位損傷，每一次計算得到所有被監測量的當前數值，每一次計算得到的所有被監測量的當前數值組成一個“被監測量的計算當前數值向量”；當假設第j根索有單位損傷時，可用Cⁱ_tj表示對應的“被監測量的當前計算數值向量Cⁱ_tj”；

h3.每一次計算得到的那個“被監測量的當前計算數值向量Cⁱ_tj”減去“被監測量的初始數值向量Cⁱ_o”得到一個向量，再將該向量的每一個元素都除以本次計算中假定的虛擬單位損傷值后得到一個“被監測量的數值變化向量”；有N根索就有N個“被監測量的數值變化向量”；

h4.由這N個“被監測量的數值變化向量”依次組成有N列的“虛擬單位損傷被監測量數值變化矩陣ΔCⁱ”；“虛擬單位損傷被監測量數值變化矩陣ΔCⁱ”的每一列對應于一個“被監測量的數值變化向量”；“虛擬單位損傷被監測量變化矩陣”的列的編號規則與當前名義虛擬損傷向量dⁱ_c和當前實際虛擬損傷向量dⁱ的元素編號規則相同；

i.實測得到索結構的所有支承索的當前索力，組成當前索力向量Fⁱ；同時，實測得到索結構的所有指定被監測量的當前實測數值，組成“被監測量的當前數值向量Cⁱ”；實測計算得到所有支承索的兩個支承端點的空間坐標，兩個支承端點的空間坐標在水平方向分量的差就是兩個支承端點水平距離；

j.定義待求的當前名義虛擬損傷向量dⁱ_c和當前實際虛擬損傷向量dⁱ；損傷向量dⁱ_o、dⁱ_c和dⁱ的元素個數等于索的數量，損傷向量的元素和索之間是一一對應關系，損傷向量的元素數值代表對應索的虛擬損傷程度或健康狀態；

k.依據“被監測量的當前數值向量Cⁱ”同“被監測量的初始數值向量Cⁱ_o”、“虛擬單位損傷被監測量數值變化矩陣ΔCⁱ”和“當前名義虛擬損傷向量dⁱ_c”間存在的近似線性關系，該近似線性關系可表達為式1，式1中除dⁱ_c外的其它量均為已知，求解式1就可以算出當前名義虛擬損傷向量dⁱ_c；

Ci=Coi+ΔCi·dci]]>式1

l.利用式2表達的當前實際虛擬損傷向量dⁱ的元素dⁱ_j同初始虛擬損傷向量dⁱ_o的元素dⁱ_oj和當前名義虛擬損傷向量dⁱ_c的元素dⁱ_cj間的關系，計算得到當前實際虛擬損傷向量dⁱ的所有元素；

dij=1-(1-doji)(1-dcji)]]>式2

式2中j＝1，2，3，……，N；

由于當前實際虛擬損傷向量dⁱ的元素數值代表對應索的當前實際虛擬損傷程度，即實際松弛程度或實際損傷程度，當前實際虛擬損傷向量dⁱ中數值不為0的元素對應的支承索就是有問題的支承索，有問題的支承索可能是松弛索、也可能是受損索，其數值反應了松弛或損傷的程度；

m.從第l步中識別出的有問題的支承索中鑒別出受損索，剩下的就是松弛索；

n.利用在第l步獲得的當前實際虛擬損傷向量dⁱ得到松弛索的當前實際虛擬損傷程度，利用在第i步獲得的當前索力向量Fⁱ，利用在第i步獲得的所有支承索的兩個支承端點的空間坐標，利用在第e步獲得的初始自由長度向量l_o，利用在第e步獲得的所有索的彈性模量、密度、初始橫截面面積數據，通過將松弛索同受損索進行力學等效來計算松弛索的、與當前實際虛擬損傷程度等效的松弛程度，等效的力學條件是：一、兩等效的索的無松弛和無損傷時的初始自由長度、幾何特性參數、密度及材料的力學特性參數相同；二、松弛或損傷后，兩等效的松弛索和損傷索的索力和變形后的總長相同；滿足上述兩個等效條件時，這樣的兩根支承索在結構中的力學功能就是完全相同的，即如果用等效的松弛索代替受損索后，索結構不會發生任何變化，反之亦然；依據前述力學等效條件求得那些被判定為松弛索的松弛程度，松弛程度就是支承索自由長度的改變量，也就是確定了那些需調整索力的支承索的索長調整量；這樣就實現了支承索的松弛識別和損傷識別，計算時所需索力由當前索力向量Fⁱ對應元素給出；

o.在求得當前名義虛擬損傷向量dⁱ_c后，按照式3建立標識向量Bⁱ，式4給出了標識向量Bⁱ的第j個元素的定義；

Bi=B1iB2i...Bij...BNiT]]>式3

Bij=0,ifdcji<Duji1,ifdcji≥Duji]]>式4

式3、式4中元素Bⁱ_j是標識向量Bⁱ的第j個元素，Dⁱ_uj是名義虛擬單位損傷向量Dⁱ_u的第j個元素，dⁱ_cj是當前名義虛擬損傷向量dⁱ_c的第j個元素，它們都表示第j根索的相關信息，式4中j＝1，2，3，……，N；

p.如果標識向量Bⁱ的元素全為0，則回到第g步繼續本次循環；如果標識向量Bⁱ的元素不全為0，則進入下一步、即第q步；

q.根據式5計算得到下一次、即第i+1次循環所需的初始虛擬損傷向量dⁱ⁺¹_o的每一個元素dⁱ⁺¹_oj；

doji+1=1-(1-doji)(1-DujiBij)]]>式5

式5中Dⁱ_uj是第i次循環名義虛擬單位損傷向量Dⁱ_u的第j個元素，dⁱ_cj是第i次循環

當前名義虛擬損傷向量dⁱ_c的第j個元素，Bⁱ_j是第i次循環標識向量Bⁱ的第j個元素；式5中j＝1，2，3，……，N；

r.在索結構當前力學計算基準模型A^ti_o的基礎上，令索的健康狀況為aⁱ⁺¹_o后更新得到下一次、即第i+1次循環所需的力學計算基準模型Aⁱ⁺¹；

s.通過對力學計算基準模型Aⁱ⁺¹的計算得到對應于模型Aⁱ⁺¹的結構的所有被監測量的數值，這些數值組成下一次、即第i+1次循環所需的被監測量的初始數值向量Cⁱ⁺¹_o；

t.建立下一次、即第i+1次循環所需的索結構當前力學計算基準模型A^ti+1_o，即取A^ti+1_o等于Aⁱ⁺¹；

u.建立下一次、即第i+1次循環所需的當前索結構支座坐標向量U^ti+1_o，即取U^ti+1_o等于U^ti_o；

v.回到步驟g，開始下一次循環。