[發明專利]支座角位移時基于應變監測的識別松弛索的遞進式方法

有效

• 申請號: 201110122643.X
• 申請日: 2011-05-13
• 公開（公告）號: CN102288432A
• 公開（公告）日: 2011-12-21
• 發明（設計）人: 韓玉林;張居鎖
• 申請（專利權）人: 東南大學
• 主分類號: G01M99/00
• 分類號: G01M99/00;G01B21/22;G01B21/32;G01L5/00
• 代理公司: 南京蘇高專利商標事務所(普通合伙) 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，在必要時更新索結構當前力學計算基準模型A^ti_o和當前索結構支座角坐標向量U^ti_o；

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

i.?實測得到索結構的所有支承索的當前索力，組成當前索力向量Fⁱ；同時，實測得到索結構的所有指定被監測量的當前實測數值，組成“被監測量的當前數值向量Cⁱ”。實測計算得到所有支承索的兩個支承端點的空間坐標，兩個支承端點的空間坐標在水平方向分量的差就是兩個支承端點水平距離。給本步及本步之前出現的所有向量的元素編號時，應使用同一編號規則，這樣可以保證本步及本步之前和之后出現的各向量的、編號相同的元素，表示同一被監測量的、對應于該元素所屬向量所定義的相關信息；

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

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

????????????????????????????????式1

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

?????????????????????????????式2

式2中j=1,2,3,……,N；

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

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

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

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

???????????????????式3

?????????????????????????????????????式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；?

???????????????????????????????式5

式5中Dⁱ_uj是第i次循環名義虛擬單位損傷向量Dⁱ_u的第j個元素，dⁱ_cj是第i次循環當前名義虛擬損傷向量d^?i_c的第j個元素，Bⁱ_j是第i次循環標識向量Bⁱ的第j個元素；式5中j=1,?2,?3,……,N；

r.?在索結構當前力學計算基準模型A^ti_o的基礎上，令索的健康狀況為dⁱ⁺¹_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，開始下一次循環。