CN103616217A  Progressive identifying method for problem cable and concentrated loads based on mixture monitoring in time of angular displacement  Google Patents
Progressive identifying method for problem cable and concentrated loads based on mixture monitoring in time of angular displacement Download PDFInfo
 Publication number
 CN103616217A CN103616217A CN201310662024.9A CN201310662024A CN103616217A CN 103616217 A CN103616217 A CN 103616217A CN 201310662024 A CN201310662024 A CN 201310662024A CN 103616217 A CN103616217 A CN 103616217A
 Authority
 CN
 China
 Prior art keywords
 cable
 temperature
 data
 vector
 initial
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Pending
Links
 238000006073 displacement reactions Methods 0.000 title claims abstract description 109
 239000000203 mixtures Substances 0.000 title abstract description 5
 230000000750 progressive Effects 0.000 title abstract 2
 238000004364 calculation method Methods 0.000 claims abstract description 74
 239000011159 matrix material Substances 0.000 claims abstract description 12
 238000011156 evaluation Methods 0.000 claims description 229
 238000005259 measurement Methods 0.000 claims description 149
 230000000875 corresponding Effects 0.000 claims description 134
 239000002965 rope Substances 0.000 claims description 121
 230000003862 health status Effects 0.000 claims description 79
 238000009826 distribution Methods 0.000 claims description 37
 230000001131 transforming Effects 0.000 claims description 31
 239000000463 material Substances 0.000 claims description 30
 238000000034 method Methods 0.000 claims description 15
 238000009659 nondestructive testing Methods 0.000 claims description 14
 241001269238 Data Species 0.000 claims description 11
 206010022114 Injury Diseases 0.000 claims description 8
 230000035852 Tmax Effects 0.000 claims description 8
 230000000694 effects Effects 0.000 claims description 7
 229910000975 Carbon steel Inorganic materials 0.000 claims description 5
 238000009413 insulation Methods 0.000 claims description 5
 238000004458 analytical method Methods 0.000 claims description 4
 229910052799 carbon Inorganic materials 0.000 claims description 4
 239000010962 carbon steel Substances 0.000 claims description 4
 239000002131 composite material Substances 0.000 claims description 4
 238000004519 manufacturing process Methods 0.000 claims description 4
 239000010902 straw Substances 0.000 claims description 4
 230000002123 temporal effect Effects 0.000 claims description 4
 238000000547 structure data Methods 0.000 claims description 2
 238000005457 optimization Methods 0.000 description 9
 238000004861 thermometry Methods 0.000 description 9
 238000007796 conventional method Methods 0.000 description 8
 230000001771 impaired Effects 0.000 description 8
 238000002910 structure generation Methods 0.000 description 4
 238000009529 body temperature measurement Methods 0.000 description 3
 238000005516 engineering process Methods 0.000 description 3
 VTYYLEPIZMXCLOUHFFFAOYSAL calcium carbonate Chemical compound data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwJyBoZWlnaHQ9JzMwMCcgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHBhdGggY2xhc3M9J2JvbmQtMCcgZD0nTSAyMDYuOTQ1LDEyMi4yODIgTCAxNzYuMTE0LDE0MC4wODMnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiNFODQyMzU7c3Ryb2tlLXdpZHRoOjIuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+CjxwYXRoIGNsYXNzPSdib25kLTAnIGQ9J00gMTc2LjExNCwxNDAuMDgzIEwgMTQ1LjI4MywxNTcuODgzJyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojM0I0MTQzO3N0cm9rZS13aWR0aDoyLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0xJyBkPSdNIDE0NS4yODMsMTU3Ljg4MyBMIDExNC40NTIsMTQwLjA4Mycgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6IzNCNDE0MztzdHJva2Utd2lkdGg6Mi4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMScgZD0nTSAxMTQuNDUyLDE0MC4wODMgTCA4My42MjA2LDEyMi4yODInIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiNFODQyMzU7c3Ryb2tlLXdpZHRoOjIuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+CjxwYXRoIGNsYXNzPSdib25kLTInIGQ9J00gMTM0LjI1MywxNTcuODgzIEwgMTM0LjI1MywxOTQuNjI1JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojM0I0MTQzO3N0cm9rZS13aWR0aDoyLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0yJyBkPSdNIDEzNC4yNTMsMTk0LjYyNSBMIDEzNC4yNTMsMjMxLjM2Nycgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6I0U4NDIzNTtzdHJva2Utd2lkdGg6Mi4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMicgZD0nTSAxNTYuMzEzLDE1Ny44ODMgTCAxNTYuMzEzLDE5NC42MjUnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiMzQjQxNDM7c3Ryb2tlLXdpZHRoOjIuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+CjxwYXRoIGNsYXNzPSdib25kLTInIGQ9J00gMTU2LjMxMywxOTQuNjI1IEwgMTU2LjMxMywyMzEuMzY3JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojRTg0MjM1O3N0cm9rZS13aWR0aDoyLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8dGV4dCB4PScxMzMuMjgzJyB5PSc0OS4xOTY4JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjQwcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojM0I0MTQzJyA+QzwvdGV4dD4KPHRleHQgeD0nMTYwLjg4MycgeT0nNDkuMTk2OCcgY2xhc3M9J2F0b20tMCcgc3R5bGU9J2ZvbnQtc2l6ZTo0MHB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6IzNCNDE0MycgPmE8L3RleHQ+Cjx0ZXh0IHg9JzE4Mi4xMzcnIHk9JzMzLjE5NjgnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MjZweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiMzQjQxNDMnID4yPC90ZXh0Pgo8dGV4dCB4PScxOTQuOTQ1JyB5PSczMy4xOTY4JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjI2cHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojM0I0MTQzJyA+KzwvdGV4dD4KPHRleHQgeD0nMjI4LjgwOCcgeT0nMTIyLjczMicgY2xhc3M9J2F0b20tMScgc3R5bGU9J2ZvbnQtc2l6ZTo0MHB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6I0U4NDIzNScgPk88L3RleHQ+Cjx0ZXh0IHg9JzI1Ni40MDgnIHk9JzEwNi43MzInIGNsYXNzPSdhdG9tLTEnIHN0eWxlPSdmb250LXNpemU6MjZweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID4tPC90ZXh0Pgo8dGV4dCB4PSczNy43NTgxJyB5PScxMjIuNzMyJyBjbGFzcz0nYXRvbS0zJyBzdHlsZT0nZm9udC1zaXplOjQwcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+TzwvdGV4dD4KPHRleHQgeD0nNjUuMzU4MScgeT0nMTA2LjczMicgY2xhc3M9J2F0b20tMycgc3R5bGU9J2ZvbnQtc2l6ZToyNnB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6I0U4NDIzNScgPi08L3RleHQ+Cjx0ZXh0IHg9JzEzMy4yODMnIHk9JzI4OC4xODYnIGNsYXNzPSdhdG9tLTQnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID5PPC90ZXh0Pgo8L3N2Zz4K data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NScgaGVpZ2h0PSc4NScgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHBhdGggY2xhc3M9J2JvbmQtMCcgZD0nTSA2Mi4zNzM1LDMxLjQyMzkgTCA1MS4zOTY2LDM3Ljc2MTQnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiNFODQyMzU7c3Ryb2tlLXdpZHRoOjEuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+CjxwYXRoIGNsYXNzPSdib25kLTAnIGQ9J00gNTEuMzk2NiwzNy43NjE0IEwgNDAuNDE5OCw0NC4wOTg5JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojM0I0MTQzO3N0cm9rZS13aWR0aDoxLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0xJyBkPSdNIDQwLjQxOTgsNDQuMDk4OSBMIDI5LjQ0MjksMzcuNzYxNCcgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6IzNCNDE0MztzdHJva2Utd2lkdGg6MS4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMScgZD0nTSAyOS40NDI5LDM3Ljc2MTQgTCAxOC40NjYxLDMxLjQyMzknIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiNFODQyMzU7c3Ryb2tlLXdpZHRoOjEuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+CjxwYXRoIGNsYXNzPSdib25kLTInIGQ9J00gMzcuMzc1LDQ0LjA5ODkgTCAzNy4zNzUsNTcuMTE1NCcgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6IzNCNDE0MztzdHJva2Utd2lkdGg6MS4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMicgZD0nTSAzNy4zNzUsNTcuMTE1NCBMIDM3LjM3NSw3MC4xMzE4JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojRTg0MjM1O3N0cm9rZS13aWR0aDoxLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0yJyBkPSdNIDQzLjQ2NDYsNDQuMDk4OSBMIDQzLjQ2NDYsNTcuMTE1NCcgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6IzNCNDE0MztzdHJva2Utd2lkdGg6MS4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMicgZD0nTSA0My40NjQ2LDU3LjExNTQgTCA0My40NjQ2LDcwLjEzMTgnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiNFODQyMzU7c3Ryb2tlLXdpZHRoOjEuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+Cjx0ZXh0IHg9JzM2Ljc2NicgeT0nMTQuNjY1OScgY2xhc3M9J2F0b20tMCcgc3R5bGU9J2ZvbnQtc2l6ZToxMnB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6IzNCNDE0MycgPkM8L3RleHQ+Cjx0ZXh0IHg9JzQ1LjE2OTcnIHk9JzE0LjY2NTknIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MTJweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiMzQjQxNDMnID5hPC90ZXh0Pgo8dGV4dCB4PSc1MS42NDExJyB5PSc5Ljc5NDIzJyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjhweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiMzQjQxNDMnID4yPC90ZXh0Pgo8dGV4dCB4PSc1NS41NDA5JyB5PSc5Ljc5NDIzJyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjhweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiMzQjQxNDMnID4rPC90ZXh0Pgo8dGV4dCB4PSc2My4xMzQ3JyB5PSczNC45NjQ1JyBjbGFzcz0nYXRvbS0xJyBzdHlsZT0nZm9udC1zaXplOjEycHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+TzwvdGV4dD4KPHRleHQgeD0nNzEuNTM4MycgeT0nMzAuMDkyOCcgY2xhc3M9J2F0b20tMScgc3R5bGU9J2ZvbnQtc2l6ZTo4cHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+LTwvdGV4dD4KPHRleHQgeD0nMTAuMzk3NCcgeT0nMzQuOTY0NScgY2xhc3M9J2F0b20tMycgc3R5bGU9J2ZvbnQtc2l6ZToxMnB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6I0U4NDIzNScgPk88L3RleHQ+Cjx0ZXh0IHg9JzE4LjgwMScgeT0nMzAuMDkyOCcgY2xhc3M9J2F0b20tMycgc3R5bGU9J2ZvbnQtc2l6ZTo4cHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+LTwvdGV4dD4KPHRleHQgeD0nMzYuNzY2JyB5PSc4MC42MzY0JyBjbGFzcz0nYXRvbS00JyBzdHlsZT0nZm9udC1zaXplOjEycHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+TzwvdGV4dD4KPC9zdmc+Cg== [Ca+2].[O]C([O])=O VTYYLEPIZMXCLOUHFFFAOYSAL 0.000 description 2
 238000004891 communication Methods 0.000 description 2
 238000007906 compression Methods 0.000 description 2
 238000001514 detection method Methods 0.000 description 2
 230000002040 relaxant effect Effects 0.000 description 2
 238000003860 storage Methods 0.000 description 2
 239000000725 suspension Substances 0.000 description 2
 230000036962 time dependent Effects 0.000 description 2
 230000000007 visual effect Effects 0.000 description 2
 229960003563 Calcium Carbonate Drugs 0.000 description 1
 206010043431 Thinking abnormal Diseases 0.000 description 1
 230000005540 biological transmission Effects 0.000 description 1
 230000015572 biosynthetic process Effects 0.000 description 1
 229910000019 calcium carbonate Inorganic materials 0.000 description 1
 238000005755 formation reaction Methods 0.000 description 1
 230000002068 genetic Effects 0.000 description 1
 230000005484 gravity Effects 0.000 description 1
 238000003384 imaging method Methods 0.000 description 1
 238000009434 installation Methods 0.000 description 1
 230000001537 neural Effects 0.000 description 1
 230000003287 optical Effects 0.000 description 1
 238000002360 preparation method Methods 0.000 description 1
 230000001681 protective Effects 0.000 description 1
 239000007787 solid Substances 0.000 description 1
 239000000700 tracer Substances 0.000 description 1
 230000017105 transposition Effects 0.000 description 1
 238000005303 weighing Methods 0.000 description 1
Abstract
The invention discloses a progressive indentifying method for a problem cable and concentrated loads based on mixture monitoring in time of angular displacement. The identifying method determines whether the mechanical calculation reference model needs to be updated or not based on the mixture monitoring through monitoring the support angular displacement, cable structure temperature and environmental temperature, so that a mechanical calculation reference mode, considering the support angular displacement, the cable structure temperature and the environmental temperature, of a cable structure is obtained. On the basis of the model, a unit damage monitored parameter numerical value matrix is obtained through calculation. According to the approximate linear relation existing between the current numerical value vector of monitored coordinates and the current initial numerical value vector of the monitored coordinates, the unit damage monitored coordinate numerical value change matrix and a current evaluated object nominal damage vector to be solved, the noninferior solution of the current evaluated object nominal damage vector is calculated, and therefore when the support angular displacement and temperature changes happen, influences of interference factors can be removed, and the variable quantity of concentrated loads and the problem cable can be identified.
Description
Technical field
Cablestayed bridge, suspension bridge, the structures such as trussframe structure have a common ground, be exactly that they have many parts that bear tensile load, as suspension cable, main pushtowing rope, hoist cable, pull bar etc., the common ground of this class formation is with rope, cable or the rod member that only bears tensile load are support unit, for simplicity, this method is " Cable Structure " by such structure representation, and by all ropeway carryingropes of Cable Structure, carrying cable, and all rod members (being called again two power rod members) that only bear axial tension or axial compression load, unification is called " cable system " for simplicity, in this method, with " support cable " this noun, censure ropeway carryingrope, carrying cable and only bear the rod member of axial tension or axial compression load, sometimes referred to as " rope ", so when using " rope " this word in the back, trussframe structure reality is just referred to two power rod members.Impaired and the lax pair Cable Structure of support cable is safely a significant threat, and this method is referred to as damaged cable and slack line the support cable of unsoundness problem, referred to as problem rope.In structure military service process, the correct identification of the health status of support cable or cable system is related to the safety of whole Cable Structure.When environment temperature changes, the temperature of Cable Structure generally also can be along with changing, when Cable Structure temperature changes, may there is angular displacement in Cable Structure bearing, the centrepoint load that Cable Structure is born also may change, the health status of Cable Structure also may change simultaneously, at this complex condition, (this method is by judging the health status of Cable Structure to the hybrid monitoring of the variation of the measurable parameter of the aforementioned dissimilar Cable Structure of this section based on hybrid monitoring for this method, this method is referred to as " monitored amount " by all monitored Cable Structure characteristic parameters, because monitored amount is now mixed and formed by the dissimilar measurable parameter of Cable Structure, this method claims that this is hybrid monitoring) carry out the variable quantity of the centrepoint load that identification problem rope and Cable Structure bear, belong to engineering structure health monitoring field.
Background technology
Reject load change, Cable Structure angular displacement of support and structure temperature and change the impact on Cable Structure health status recognition result, thereby the variation of the health status of recognition structure is exactly current problem in the urgent need to address; Same, the impact of the recognition result of the variable quantity of the centrepoint load that the variation of rejecting structure temperature, Cable Structure angular displacement of support and structural health conditions variation are born structure, significant equally to structural safety, this method discloses a kind of effective ways that solve these two problems.
Support cable is impaired, lax pair Cable Structure is safely a significant threat, and the problem rope of identifying based on structural health monitoring technology in the cable system of Cable Structure is a kind of method that has potentiality.
When changing appears in the centrepoint load of bearing when Cable Structure, or Cable Structure angular displacement of support, or the temperature of Cable Structure is when change, for example, or the health status of cable system is while changing (damaging), or when four kinds of situations occur simultaneously, can cause the variation of the measurable parameter of Cable Structure, for example can cause the variation of Suo Li, can affect distortion or the strain of Cable Structure, can affect shape or the volume coordinate of Cable Structure, can cause variation (for example variation of the angle coordinate of the straight line of any this point of mistake in the section of body structure surface any point of angle coordinate of any imaginary line of the every bit of Cable Structure, or the variation of the angle coordinate of the normal of body structure surface any point), all these change the health status information that has all comprised cable system, also the variable quantity information that has comprised centrepoint load, that is to say the variable quantity that can utilize the measurable parameter of Cable Structure to identify damaged cable and centrepoint load.
When bearing has angular displacement, current published technology, in method, some only can be when other all conditions be constant the variation of (load of only only having structure to bear changes) recognition structure bearing load, the variation of some recognition structure health status of only can (only only having structural health conditions to change) when other all conditions is constant, the variation of some only can (only only have structure temperature and structural health conditions to change) when other all conditions is constant recognition structure (environment) temperature and structural health conditions, also do not have at present a kind of disclosed, effective method is recognition structure bearing load simultaneously, the variation of structure (environment) temperature and structural health conditions, when the load of bearing in structure in other words and structure (environment) temperature changes simultaneously, also there is no the variation that effective method can recognition structure health status, and the load that structure is born and structure (environment) temperature usually changes, so during the load of how to bear in structure and structure (environment) temperature variation, reject load change and structure temperature and change the impact on Cable Structure health status recognition result, thereby the variation of the health status of recognition structure exactly, it is current problem in the urgent need to address, this method discloses a kind of method, when bearing has angular displacement, when the centrepoint load that can bear in Cable Structure and structure (environment) temperature changes, reject angular displacement of support, load change and structure temperature change the impact on Cable Structure health status recognition result, based on monitored amount, monitor identification problem rope, the safety of Cable Structure is had to important value.
Same, in current disclosed method, thereby also do not occur rejecting the correct knowledge method for distinguishing of realizing centrepoint load intensity of variation of angular displacement of support, structure temperature variation and the impact of support cable health status, and concerning structure, the identification of load change is also very important.This method, when identifying problem rope, can also identify the variation of centrepoint load simultaneously, and this method can be rejected angular displacement of support, structure temperature changes and the impact of support cable health status variation, realizes the correct identification of centrepoint load intensity of variation.
That is to say, this method has realized existing method can not possess function.
Summary of the invention
Technical matters: this method discloses a kind of method, two kinds of functions that existing method can not possess have been realized, be respectively, one, when bearing has angular displacement, during the centrepoint load of bearing in structure and structure (environment) temperature variation, can reject angular displacement of support, centrepoint load variation and structure temperature and change the impact on Cable Structure health status recognition result, thereby identify exactly the health status of support cable; Two, this method, when identifying problem rope, can also identify the variation of centrepoint load simultaneously, and this method can be rejected angular displacement of support, structure temperature changes and the impact of support cable health status variation, realizes the correct identification of centrepoint load intensity of variation.
Technical scheme: this method is comprised of three parts.Respectively: one, " the temperature survey calculating method of the Cable Structure of this method "; Two, set up method, the structural health conditions appraisal procedure based on knowledge base (containing parameter) and the monitored amount of actual measurement of the required knowledge base of cable structure health monitoring system and parameter; Three, the software and hardware part of health monitoring systems.
In the method, with " bearing volume coordinate ", censure bearing about the coordinate of the X, Y, Z axis of Descartes's rectangular coordinate system, also can be said to is that bearing is about the volume coordinate of X, Y, Z axis, bearing is called bearing about the volume coordinate component of this axle about the concrete numerical value of the volume coordinate of some axles, in this method, also with a volume coordinate component of bearing, expresses bearing about the concrete numerical value of the volume coordinate of some axles; With " bearing angular coordinate ", censure bearing about the angular coordinate of X, Y, Z axis, bearing is called bearing about the angular coordinate component of this axle about the concrete numerical value of the angular coordinate of some axles, in this method, also with an angular coordinate component of bearing, expresses bearing about the concrete numerical value of the angular coordinate of some axles; All by " bearing generalized coordinate " denotion bearing angular coordinate and bearing volume coordinate, in this method, also with a generalized coordinate component of bearing, express bearing about the concrete numerical value of volume coordinate or the angular coordinate of an axle; Bearing is called support wire displacement about the change of the coordinate of X, Y, Z axis, also can say that the change of bearing volume coordinate is called support wire displacement, in this method, also with a translational component of bearing, expresses bearing about the concrete numerical value of the displacement of the lines of some axles; Bearing is called angular displacement of support about the change of the angular coordinate of X, Y, Z axis, in this method, also with an angular displacement component of bearing, expresses bearing about the concrete numerical value of the angular displacement of some axles; Generalized displacement of support denotion support wire displacement and angular displacement of support are all, in this method, also with a generalized displacement component of bearing, express bearing about the displacement of the lines of some axles or the concrete numerical value of angular displacement; Support wire displacement also can be described as translational displacement, and support settlement is that support wire displacement or translational displacement are at the component of gravity direction.
First confirm the quantity of the centrepoint load that may change that Cable Structure is born.The feature of the centrepoint load of bearing according to Cable Structure, confirm wherein " centrepoint load likely changing ", or all centrepoint load is considered as " centrepoint load likely changing ", establishes total JZW the centrepoint load that may change.
Centrepoint load is divided into two kinds of concentrated force and concentrated couples, in coordinate system, for example, in Descartes's rectangular coordinate system, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, in the method a concentrated force component or a concentrated couple component is called to a centrepoint load.
If the quantity sum of the quantity of the support cable of Cable Structure and JZW " centrepoint load likely changing " is N.For sake of convenience, it is " evaluation object " that this method unitedly calls evaluated support cable and " centrepoint load likely changing ", total N evaluation object.Give evaluation object serial number, this numbering will be for generating vector sum matrix in subsequent step.
Monitored multiclass parameter can comprise: Suo Li, strain, angle and volume coordinate, be described below respectively:
If total Q root support cable in cable system, the monitored rope force data of Cable Structure is by M in Cable Structure
_{1}the M of individual appointment rope
_{1}individual rope force data is described, and the variation of Cable Structure Suo Li is exactly the variation of the Suo Li of all appointment ropes.Each total M
_{1}individual cable force measurement value or calculated value characterize the rope force information of Cable Structure.M
_{1}be one and be not less than 0 integer.
The monitored strain data of Cable Structure can be by K in Cable Structure
_{2}l individual specified point and each specified point
_{2}the strain of individual assigned direction is described, and the variation of Cable Structure strain data is exactly K
_{2}the variation of all tested strains of individual specified point.Each total M
_{2}(M
_{2}=K
_{2}* L
_{2}) individual strain measurement value or calculated value characterize Cable Structure strain.M
_{2}be one and be not less than 0 integer.
The monitored angledata of Cable Structure is by K in Cable Structure
_{3}l individual specified point, that cross each specified point
_{3}h individual appointment straight line, each appointment straight line
_{3}individual angle coordinate component is described, and the variation of Cable Structure angle is exactly the variation of angle coordinate components appointment straight lines all specified points, all, all appointments.Each total M
_{3}(M
_{3}=K
_{3}* L
_{3}* H
_{3}) individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure.M
_{3}be one and be not less than 0 integer.
The monitored shape data of Cable Structure is by K in Cable Structure
_{4}l individual specified point and each specified point
_{4}the volume coordinate of individual assigned direction is described, and the variation of Cable Structure shape data is exactly K
_{4}the variation of all coordinate components of individual specified point.Each total M
_{4}(M
_{4}=K
_{4}* L
_{4}) individual measurement of coordinates value or calculated value characterize Cable Structure shape.M
_{4}be one and be not less than 0 integer.
Comprehensive abovementioned monitored amount, whole Cable Structure has M(M=M
_{1}+ M
_{2}+ M
_{3}+ M
_{4}) individual monitored amount, definition parameter K(K=M
_{1}+ K
_{2}+ K
_{3}+ K
_{4}), K and M must not be less than N.
For simplicity, in the method by " monitored all parameters of Cable Structure " referred to as " monitored amount ".Give M monitored amount serial number, this numbering will be for generating vector sum matrix in subsequent step.This method is with representing this numbering with variable j, j=1, and 2,3 ..., M.
The first of this method: " the temperature survey calculating method of the Cable Structure of this method ".
First determine " the temperature survey calculating method of the Cable Structure of this method ".Because the temperature of Cable Structure may change, for example the temperature of the different parts of Cable Structure is to change along with the variation of intensity of sunshine, along with the variation of environment temperature changes, the surface of Cable Structure and inner temperature may be time dependent sometimes, the surface of Cable Structure may be different from inner temperature, the surface of Cable Structure and inner temperature difference are time dependent, the Mechanics Calculation of Cable Structure when this just makes to consider temperature conditions and monitoring very complex, for simplification problem, reduce calculated amount and reduce and measure cost, especially in order to improve computational accuracy, this method proposes " the temperature survey calculating method of the Cable Structure of this method ", specific as follows:
The first step, inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of Cable Structure composition material and Cable Structure environment of living in, utilize the geometry measured data of design drawing, asconstructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure.Inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, statistics obtains interior during this period of time cloudy quantity and is designated as T cloudy day, statistics obtains 0 highest temperature and the lowest temperature between latter 30 minutes of the moment of sunrise next day at each cloudy day in T cloudy day, sunrise on the meteorology that sunrise refers to base area revolutions constantly and the rule that revolves round the sun is definite constantly, the sunrise that can inquire about data or calculate each required day by conventional meteorology constantly, each cloudy day 0 up to next day sunrise constantly the highest temperature between latter 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, there is T cloudy day, the maximum temperature difference that just has the daily temperature at T cloudy day, get maximal value in the maximum temperature difference of daily temperature at T cloudy day for reference to temperature difference per day, with reference to temperature difference per day, be designated as Δ T
_{r}.Between inquiry Cable Structure location and Altitude Region, place, be no less than temperature that the meteorological data in recent years of 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and Altitude Region, place about the maximum rate of change Δ T of sea level elevation
_{h}, for convenience of narration, get Δ T
_{h}unit be ℃/m.On the surface of Cable Structure, get " R Cable Structure surface point ", after will by actual measurement, obtain the temperature of this R Cable Structure surface point, claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, by Calculation of Heat Transfer, obtain the temperature of this R Cable Structure surface point, just claim that the temperature data calculating is " R Cable Structure surface temperature computational data ".While getting " R Cable Structure surface point " on the surface of Cable Structure, the quantity of " R Cable Structure surface point " is narrated in the back with the condition that distribution must be satisfied.From the residing minimum height above sea level of Cable Structure to the highest height above sea level, in Cable Structure, uniform choosing is no less than three different sea level elevations, the sea level elevation place choosing at each, at the intersection place on surface level and Cable Structure surface, at least choose two points, outer normal from selected point straw line body structure surface, all outer normal directions of choosing are called " measuring Cable Structure along the direction of the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in the measurement Cable Structure of choosing along comprising the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure in the direction of the Temperature Distribution of wall thickness, direction uniform choosing in Cable Structure along each measurement Cable Structure along the Temperature Distribution of wall thickness is no less than three points, measure all temperature that are selected a little, the temperature recording is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure Cable Structure along the direction of the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if chosen H different sea level elevation, at each sea level elevation place, choose B and measured Cable Structure along the direction of the Temperature Distribution of wall thickness, along each, measure Cable Structure and in Cable Structure, chosen E point along the direction of the Temperature Distribution of wall thickness, wherein H and E are not less than 3, B is not less than 2, if HBE is the product of H and B and E, corresponding total HBE " measuring Cable Structure along the point of the temperature profile data of thickness ", after will obtain the temperature that this HBE " measures Cable Structure along the point of the temperature profile data of thickness " by actual measurement, claim that the temperature data that actual measurement obtains is " HBE Cable Structure is along thickness temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, by Calculation of Heat Transfer, obtain this HBE and measure Cable Structure along the temperature of the point of the temperature profile data of thickness, just claim that the temperature data calculating is " HBE Cable Structure is along thickness temperature computation data ", if BE is the product of B and E, total BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in sea level elevation place of choosing at each in this method.In Cable Structure location, according to meteorology, measure temperature and require to choose a position, will obtain meeting the temperature that meteorology is measured the Cable Structure place environment of temperature requirement in this position actual measurement, in the onsite spaciousness of Cable Structure, unobstructed place chooses a position, this position should can obtain in each day of the whole year this ground the most sufficient sunshine of getable this day, flat board at a carbon steel material of this position of sound production, be called reference plate, the one side of this reference plate on the sunny side, be called sunny slope, the sunny slope of reference plate is coarse and dark color, the sunny slope of reference plate should can obtain in each day of the whole year one flat plate on this ground the most sufficient sunshine of getable this day, the nonsunny slope of reference plate is covered with insulation material, RealTime Monitoring is obtained to the temperature of the sunny slope of reference plate.In this method, to the time interval between any twice measurement of same amount RealTime Monitoring, must not be greater than 30 minutes, the moment of survey record data is called physical record data constantly.
Second step, RealTime Monitoring obtains R Cable Structure surface temperature measured data of abovementioned R Cable Structure surface point, RealTime Monitoring obtains previously defined Cable Structure along the temperature profile data of thickness simultaneously, and RealTime Monitoring obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement simultaneously, by RealTime Monitoring, obtain being carved at sunrise the same day sunrise next day temperature measured data sequence of the Cable Structure place environment between latter 30 minutes constantly, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data that was carved at sunrise the Cable Structure place environment between latter 30 minutes of the moment of sunrise next day the same day, find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment, by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, deduct and be carved at sunrise the sunrise next day maximum temperature difference between latter 30 minutes constantly on same day that minimum temperature obtains Cable Structure place environment, be designated as Δ T
_{emax}, by the temperature measured data sequence of Cable Structure place environment, by conventional mathematical computations, obtain the temperature of Cable Structure place environment about the rate of change of time, this rate of change is also along with the time changes, by RealTime Monitoring, obtain being carved at sunrise the same day sunrise next day measured data sequence of the temperature of the sunny slope of the reference plate between latter 30 minutes constantly, the measured data sequence of the temperature of the sunny slope of reference plate by be carved at sunrise the same day next day sunrise constantly the measured data of the temperature of the sunny slope of the reference plate between latter 30 minutes according to time order and function order, arrange, find maximum temperature and minimum temperature in the measured data sequence of temperature of sunny slope of reference plate, by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, deduct and be carved at sunrise the sunrise next day maximum temperature difference between latter 30 minutes constantly on same day of temperature that minimum temperature obtains the sunny slope of reference plate, be designated as Δ T
_{pmax}, by RealTime Monitoring, obtain being carved at sunrise the same day sunrise next day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between latter 30 minutes constantly, there is R Cable Structure surface point just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence by be carved at sunrise on same day of a Cable Structure surface point sunrise next day constantly the Cable Structure surface temperature measured data between latter 30 minutes according to time order and function order, arrange, find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence, by the maximum temperature in each Cable Structure surface temperature measured data sequence, deduct and be carved at sunrise the sunrise next day maximum temperature difference between latter 30 minutes constantly on same day that minimum temperature obtains the temperature of each Cable Structure surface point, there is R Cable Structure surface point just to have and be carved at sunrise the sunrise next day maximum temperature difference numerical value between latter 30 minutes constantly R the same day, maximal value is wherein designated as Δ T
_{smax}, by each Cable Structure surface temperature measured data sequence, by conventional mathematical computations, obtain the temperature of each Cable Structure surface point about the rate of change of time, the temperature of each Cable Structure surface point about the rate of change of time also along with the time changes.By RealTime Monitoring, obtain being carved at sunrise the same day between latter 30 minutes of the moment of sunrise next day, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", the sea level elevation place that calculating is chosen at each amounts to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and the difference of minimum temperature, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", chosen H different sea level elevation and just had H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T
_{tmax}.
The 3rd step, measures and calculates acquisition Cable Structure steady temperature data, first, determine the moment that obtains Cable Structure steady temperature data, the condition relevant to the moment that determines acquisition Cable Structure steady temperature data has six, first condition is the moment that obtains Cable Structure steady temperature data to be carved at sunset sunrise next day constantly between latter 30 minutes between the same day, sunset refers to sunset on base area revolutions and the definite meteorology of revolution rule constantly constantly, and the sunset that can inquire about data or calculate each required day by conventional meteorology constantly, the a condition of second condition be the same day be carved at sunrise next day sunrise constantly between latter 30 minutes during this period of time in, reference plate maximum temperature difference Δ T
_{pmax}with Cable Structure surface maximum temperature difference Δ T
_{smax}all be not more than 5 degrees Celsius, the b condition of second condition be the same day be carved at sunrise next day sunrise constantly between latter 30 minutes during this period of time in, above, measure the environment maximum error Δ T calculating
_{emax}be not more than with reference to temperature difference per day Δ T
_{r}, and reference plate maximum temperature difference Δ T
_{pmax}after deducting 2 degrees Celsius, be not more than Δ T
_{emax}, and Cable Structure surface maximum temperature difference Δ T
_{smax}be not more than Δ T
_{pmax}, only needing to meet in a condition of second and b condition one is just called and meets second condition, the 3rd condition is that the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, the 4th condition is that the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, the 5th condition is in the moment that obtains Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is to be carved at sunrise the sunrise next day minimal value between latter 30 minutes constantly the same day, the 6th condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T
_{tmax}be not more than 1 degree Celsius, this method is utilized abovementioned six conditions, any one in following three kinds of moment is called to " obtaining the mathematics of Cable Structure steady temperature data constantly ", the first is first moment to the 5th condition meeting in abovementioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data " constantly, the second is the moment that only meets the 6th condition in abovementioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data " constantly, the third is first moment to the 6th condition simultaneously meeting in abovementioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data " constantly, when obtaining the mathematics of Cable Structure steady temperature data, be exactly in this method during in constantly one of physical record data constantly, the moment that obtains Cable Structure steady temperature data be exactly obtain Cable Structure steady temperature data mathematics constantly, if obtain the mathematics of Cable Structure steady temperature data and be not constantly any in constantly of physical record data in this method constantly, get this method close to moment of the mathematics that obtains Cable Structure steady temperature data those physical record data constantly for obtaining the moment of Cable Structure steady temperature data, this method is carried out the relevant health monitoring analysis of Cable Structure by using in the amount that obtains the moment survey record of Cable Structure steady temperature data, this method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, i.e. this Cable Structure temperature constantly temporal evolution not, and this is exactly " obtaining the moment of Cable Structure steady temperature data " of this method constantly, then, according to Cable Structure heat transfer characteristic, utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment that obtains Cable Structure steady temperature data, utilize the thermal conduction study computation model of Cable Structure, by conventional Calculation of Heat Transfer, obtain in the Temperature Distribution of Cable Structure that obtains the moment of Cable Structure steady temperature data, now calculate by stable state in the temperature field of Cable Structure, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data, also comprise that Cable Structure is in the accounting temperature of selected HBE " measuring Cable Structure along the point of the temperature profile data of thickness " above, the accounting temperature of HBE " measuring Cable Structure along the point of the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", when R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure stable state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", while getting " R Cable Structure surface point " on the surface of Cable Structure, the quantity of " R Cable Structure surface point " and necessary three conditions that meet that distribute, first condition is when Cable Structure temperature field is during in stable state, on Cable Structure surface the temperature of any point be by " R Cable Structure surface point " with Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point while obtaining, on the Cable Structure surface that linear interpolation obtains, on the temperature of this arbitrfary point and Cable Structure surface, the error of the actual temperature of this arbitrfary point is not more than 5%, Cable Structure surface comprises support cable surface, second condition is that in " R Cable Structure surface point ", the quantity at the point of same sea level elevation is not less than 4, and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 ℃ divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation
_{h}the numerical value obtaining, gets Δ T for convenience of narration
_{h}unit be ℃/m that the unit of getting Δ h for convenience of narration is m, " R Cable Structure surface point " refers to while only considering sea level elevation along the definition of adjacent Cable Structure surface point between two of sea level elevation, in " R Cable Structure surface point ", do not have a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two, the 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to the geometric properties of Cable Structure and bearing data, in Cable Structure, find and be subject to the sunshineduration position of those surface points the most fully the whole year, in " R Cable Structure surface point ", having a Cable Structure surface point at least is the annual point being subject in the most sufficient those surface points of sunshineduration in Cable Structure.
The second portion of this method: method, the structural health conditions appraisal procedure based on knowledge base (containing parameter) and the monitored amount of actual measurement of setting up the required knowledge base of cable structure health monitoring system and parameter.Can carry out successively as follows, to obtain the health status assessment of evaluation object more accurately.
The first step: set up initial mechanical calculating benchmark model A
_{o}in Cable Structure completion, or before setting up health monitoring systems, according to " the temperature survey calculating method of the Cable Structure of this method " measurement, calculating " Cable Structure steady temperature data " (can measure by conventional thermometry, for example use thermal resistance to measure), " Cable Structure steady temperature data " now use vector T
_{o}represent, be called initial Cable Structure steady temperature data vector T
_{o}.In actual measurement, obtain T
_{o}time, namely obtaining initial Cable Structure steady temperature data vector T
_{o}the synchronization in the moment, use conventional method directly to measure the initial number of all monitored amounts that calculate Cable Structure.At Actual measurement, obtain initial Cable Structure steady temperature data vector T
_{o}time, use conventional method (consult reference materials or survey) to obtain temperature variant physical parameter (for example thermal expansivity) and the mechanical property parameters (for example elastic modulus, Poisson ratio) of the various materials that Cable Structure used; At Actual measurement, obtain initial Cable Structure steady temperature data vector T
_{o}time, namely obtaining initial Cable Structure steady temperature data vector T
_{o}the synchronization in the moment, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.First the Actual measurement data of Cable Structure are the data of the health status that can express rope that comprise the Nondestructive Testing Data of support cable, and the Actual measurement data of Cable Structure still comprise the measured data of the initial geometric data of Cable Structure, rope force data, drawbar pull data, initial Cable Structure bearing generalized coordinate data (comprising initial Cable Structure bearing spatial data and initial Cable Structure bearing angular data), Cable Structure modal data, structural strain data, structure angle measurement data, structure space measurement of coordinates data, load data.The initial geometric data of Cable Structure can be the spatial data that the spatial data of the end points of all ropes adds a series of point in structure, and object is to determine according to these coordinate datas the geometric properties of Cable Structure.For cablestayed bridge, the spatial data that initial geometric data can be the end points of all ropes adds the spatial data of some points on bridge two ends, socalled bridge type data that Here it is.The variable quantity of " centrepoint load likely changing " is being set up initial mechanical calculating benchmark model A
_{o}time be all 0, the variable quantity that is to say " centrepoint load likely changing " that identifies is below with respect to setting up initial mechanical calculating benchmark model A
_{o}time the structure corresponding centrepoint load of bearing variable quantity.The data of health status and the variable quantity data of " centrepoint load likely changing " of utilizing the Nondestructive Testing Data etc. of support cable can express support cable are set up evaluation object initial damage vector d
_{o}(as the formula (1)), use d
_{o}represent that Cable Structure is (with initial mechanical calculating benchmark model A
_{o}the initial health of evaluation object expression).If while there is no the Nondestructive Testing Data of support cable and the data of other health status that can express support cable, or can think that structure original state is not damaged during without relaxed state, vectorial d
_{o}in each element numerical value relevant to support cable get 0.Vector d
_{o}in each element numerical value relevant to the variable quantity of centrepoint load get 0.Utilize the design drawing, asconstructed drawing of Cable Structure and the initial measured data of Cable Structure, temperature variant physical and mechanical properties parameter and the initial Cable Structure steady temperature data vector T of the various materials that the Nondestructive Testing Data of support cable, Cable Structure are used
_{o}, utilize mechanics method (for example finite element method) to count " Cable Structure steady temperature data " and set up initial mechanical calculating benchmark model A
_{o}.
Corresponding to A
_{o}cable Structure bearing angular data form initial Cable Structure bearing angular coordinate vector U
_{o}.
d
_{o}＝[d
_{o1}?d
_{o2}?·?·?·?d
_{ok}?·?·?·?d
_{oN}]
^{T}?????(1)
D in formula (1)
_{ok}(k=1,2,3 ...., N) represent initial mechanical calculating benchmark model A
_{o}in the original state of k evaluation object, if this evaluation object is the rope (or pull bar) in cable system, so d
_{ok}represent its initial damage, d
_{ok}be to represent not damaged at 0 o'clock, while being 100%, represent that this rope thoroughly loses loadbearing capacity, in the time of between 0 and 100%, represent to lose the loadbearing capacity of corresponding proportion; If this evaluation object is " centrepoint load that may change ", so a d
_{ok}represent its initial value, d
_{ok}be 0, the variable quantity that is to say " centrepoint load likely changing " that identifies is below with respect to setting up initial mechanical calculating benchmark model A
_{o}time the structure corresponding centrepoint load of bearing variable quantity.In formula, subscript T represents vectorial transposition (same afterwards).
In actual measurement, obtain T
_{o}time, namely at the synchronization that obtains the moment of Cable Structure steady temperature data, use conventional method directly to measure the initial value of all monitored amounts of the Cable Structure calculating, form monitored amount initial value vector C
_{o}(seeing formula (2)).Requirement is obtaining A
_{o}time obtain C
_{o}, monitored amount initial value vector C
_{o}expression is corresponding to A
_{o}the concrete numerical value of " monitored amount ".Because of subject to the foregoing, the calculating benchmark model based on Cable Structure calculates the monitored amount of gained reliably close to the measured data of initial monitored amount, in narration below, will represent this calculated value and measured value with prosign.
C
_{o}＝[C
_{o1}?C
_{o2}?·?·?·?C
_{oj}?·?·?·?C
_{oM}]
^{T}?????(2)
C in formula (2)
_{oj}(j=1,2,3 ...., M) be the original bulk of j monitored amount in Cable Structure, this component according to coding rule corresponding to specific j monitored amount.Vector C
_{o}be to arrange and form according to a definite sequence by the monitored amount of M, this put in order and there is no specific (special) requirements, only require below also array data in this order of all associated vector.
No matter which kind of method to obtain initial mechanical calculating benchmark model A by
_{o}, counting " Cable Structure steady temperature data " (is initial Cable Structure steady temperature data vector T
_{o}), based on A
_{o}the Cable Structure computational data calculating must approach its measured data very much, and error generally must not be greater than 5%.Like this can utility A
_{o}calculate Suo Li computational data, strain computational data, Cable Structure shape computational data and displacement computational data under the analog case of gained, Cable Structure angledata, Cable Structure spatial data etc., the measured data when approaching reliably institute's analog case and truly occurring.Model A
_{o}evaluation object initial damage vector d for the health status of middle evaluation object
_{o}represent model A
_{o}the vectorial U of middle bearing angular coordinate
_{o}represent initial Cable Structure steady temperature data vector T for Cable Structure steady temperature data
_{o}represent.Due to based on A
_{o}the evaluation that calculates all monitored amounts approaches the initial value (actual measurement obtains) of all monitored amounts very much, so also can be used in A
_{o}basis on, carry out Mechanics Calculation obtains, A
_{o}the evaluation of each monitored amount form monitored amount initial value vector C
_{o}.U
_{o}, T
_{o}and d
_{o}a
_{o}parameter, also can say C
_{o}by A
_{o}mechanics Calculation result form.
Second step: circulation starts.When circulation starts each time, the current initial damage vector of the evaluation object d while first needing to set up or set up this circulation beginning
^{i} _{o}(i=1,2,3 ...), set up the current initial mechanical calculating benchmark model A of Cable Structure
^{i} _{o}(finite element benchmark model for example, A in circulation each time
^{i} _{o}constantly update), A
^{i} _{o}" current initial Cable Structure steady temperature data vector T for Temperature Distribution
^{i} _{o}" express.Letter i, except representing that significantly, the place of number of steps, alphabetical i only represents cycle index in the method, circulates for the i time.A
_{o}and A
^{i} _{o}count temperature parameter, can accounting temperature change the Effect on Mechanical Properties to Cable Structure.
The current initial damage vector of evaluation object that the i time circulation needs while starting is designated as d
^{i} _{o}(as the formula (3)), use d
^{i} _{o}while representing this circulation beginning, Cable Structure is (with current initial mechanical calculating benchmark model A
^{i} _{o}the health status of evaluation object expression).
D in formula (3)
^{i} _{ok}(i=1,2,3, K=1,2,3 ...., while N) representing that the i time circulation starts, current initial mechanical calculating benchmark model A
^{i} _{o}in the original state of k evaluation object, if this evaluation object is the rope (or pull bar) in cable system, so d
^{i} _{ok}represent its initial damage, d
^{i} _{ok}be to represent not damaged at 0 o'clock, while being 100%, represent that this rope thoroughly loses loadbearing capacity, in the time of between 0 and 100%, represent to lose the loadbearing capacity of corresponding proportion, if this evaluation object is " centrepoint load that may change ", so a d
^{i} _{ok}represent that it is with respect to setting up initial mechanical calculating benchmark model A
_{o}time the structure corresponding centrepoint load of bearing variable quantity.
Current initial mechanical calculating benchmark model A corresponding to Cable Structure
^{i} _{o}cable Structure bearing angular data form current initial Cable Structure bearing angular coordinate vector U
^{i} _{o}, at initial time, namely set up for the first time the current initial mechanical calculating benchmark model A of Cable Structure
^{i} _{o}time, U
^{i} _{o}just equal U
_{o}.
Set up and renewal d
^{i} _{o}method as follows:
When circulation starts for the first time, set up the current initial damage vector of evaluation object and (according to formula (3), be designated as d
^{1} _{o}) time, d
^{1} _{o}just equal d
_{o}.I (i=2,3,4,5,6 ...) the current initial damage vector of the evaluation object d that needs while starting of inferior circulation
^{i} _{o}, be front once (the i1 time, i=2,3,4,5,6 ...) circulation finishes that front calculating obtains, concrete grammar is described below.
I (i=1,2,3,4,5,6 ...) the Mechanics Calculation benchmark model of the Mechanics Calculation benchmark model that inferior circulation need to be set up while starting or the Cable Structure of having set up is designated as current initial mechanical calculating benchmark model A
^{i} _{o}.Corresponding to A
^{i} _{o}" Cable Structure steady temperature data " use vector T
^{i} _{o}represent, be called current initial Cable Structure steady temperature data vector T
^{i} _{o}.Vector T
^{i} _{o}definition mode and vector T
_{o}definition mode identical, when circulation starts each time, must set up or set up current initial Cable Structure steady temperature data vector T
^{i} _{o}.
Set up, upgrade A
^{i} _{o}, U
^{i} _{o}and T
^{i} _{o}method as follows:
The Mechanics Calculation benchmark model of the Cable Structure of setting up when circulation starts is for the first time designated as A
^{1} _{o}, A
^{1} _{o}equal A
_{o}, T
^{1} _{o}equal T
_{o}.A in circulation each time
^{i} _{o}, U
^{i} _{o}and T
^{i} _{o}constantly update, concrete grammar is described below; When circulation end each time, upgrade A
^{i} _{o}, U
^{i} _{o}and T
^{i} _{o}the Mechanics Calculation benchmark model of required Cable Structure while obtaining next time circulating beginning, concrete grammar is described below.
" the current initial value vector of monitored amount C for this method
^{i} _{o}" (i=1,2,3 ...) represent the initial value (referring to formula (4)) of the monitored amount of all appointments when the i time (i=1,2,3,4,5,6 ...) circulation starts, C
^{i} _{o}also can be called " the i time current initial value of the monitored amount of circulation vector ".
C in formula (2)
^{i} _{oj}(i=1,2,3, J=1,2,3 ...., j monitored amount while M) being the i time circulation beginning, in Cable Structure.Vector C
^{i} _{o}be to arrange and form according to a definite sequence by the monitored amount of previously defined M, this put in order and there is no specific (special) requirements, only require below also array data in this order of all associated vector.
Setting up model A
^{i} _{o}time set up " the current initial value vector of monitored amount C
^{i} _{o}", the current initial value vector of monitored amount C
^{i} _{o}expression is corresponding to A
^{i} _{o}the concrete numerical value of all monitored amounts, C
^{i} _{o}element and C
_{o}element corresponding one by one, represent respectively all monitored amounts in Cable Structure in A
^{i} _{o}and A
_{o}concrete numerical value during two states.
Set up and renewal C
^{i} _{o}concrete grammar as follows:
When circulation starts for the first time, C
^{1} _{o}(i=1, C
^{i} _{o}be embodied as C
^{1} _{o}) equal C
_{o}; I (i=2,3,4,5,6 ...) the i time circulation " vectorial C of the current initial value of monitored amount of needing while starting of inferior circulation
^{i} _{o}", be front once (the i1 time, i=2,3,4,5,6 ...) circulation finishes that front calculating obtains, concrete grammar is described below.The i time (i=1,2,3,4,5,6 ...) in circulation, " the current initial value vector of monitored amount C
^{i} _{o}" constantly update, concrete grammar is described below.Due to according to model A
^{i} _{o}calculate the initial value of the monitored amount of gained reliably close to corresponding measured value, in narration below, will represent this calculated value composition of vector and measured value composition of vector with prosign.
T
^{i} _{o}, U
^{i} _{o}and d
^{i} _{o}a
^{i} _{o}characterisitic parameter, C
^{i} _{o}a
^{i} _{o}mechanics Calculation result form.
The 3rd step: in Cable Structure military service process, in circulation each time, in other words in i (i=1,2,3,4,5,6 ...) in inferior circulation, at known A
^{i} _{o}, T
^{i} _{o}, U
^{i} _{o}, C
^{i} _{o}and d
^{i} _{o}after, the current data that obtains " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement, the current data of all " Cable Structure steady temperature data " forms " current cable structure steady temperature data vector T
^{i}", vector T
^{i}definition mode and vector T
_{o}definition mode identical; In actual measurement vector T
^{i}time, namely obtaining current cable structure steady temperature data vector T
^{i}the synchronization in the moment, actual measurement obtains the currency of all monitored amounts in Cable Structure, all these numerical value form monitored amount current value vector C
^{i}.C
^{i}element and C
_{o}element corresponding one by one, represent that identical monitored amount is at numerical value in the same time not.Obtaining vector T
^{i}time, actual measurement obtains Cable Structure bearing angular coordinate current data, and all Cable Structure bearing angular coordinate current datas form current cable structure actual measurement bearing angular coordinate vector U
^{i}.
Obtaining vector T
^{i}after, according to following concrete grammar, upgrade A
^{i} _{o}, T
^{i} _{o}, U
^{i} _{o}, C
^{i} _{o}and d
^{i} _{o}:
Compare respectively T
^{i}and T
^{i} _{o}, U
^{i}and U
^{i} _{o}if, T
^{i}equal T
^{i} _{o}and U
^{i}equal U
^{i} _{o}, do not need A
^{i} _{o}upgrade, otherwise need to be to A
^{i} _{o}, U
^{i} _{o}and T
^{i} _{o}upgrade, update method is: the first step is calculated U
^{i}with U
_{o}poor, U
^{i}with U
_{o}difference be exactly that Cable Structure bearing is about the angular displacement of support of initial position, with angular displacement of support vector, V represents angular displacement of support, between element in angular displacement of support vector V and angular displacement of support component, be onetoone relationship, in angular displacement of support vector V, the numerical value of an element is corresponding to the angular displacement of an assigned direction of an appointment bearing; Second step calculates T
^{i}with T
_{o}poor, T
^{i}with T
_{o}difference be exactly that current cable structure steady temperature data are about the variation of initial Cable Structure steady temperature data, T
^{i}with T
_{o}poor with steady temperature change vector S, represent, S equals T
^{i}deduct T
_{o}, S represents the variation of Cable Structure steady temperature data; The 3rd step is first to A
_{o}in Cable Structure bearing apply angular displacement of support constraint, the numerical value of angular displacement of support constraint is just taken from the numerical value of corresponding element in angular displacement of support vector V, then to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation applying is just taken from steady temperature change vector S, to A
_{o}middle Cable Structure bearing applies angular displacement of support constraint and to A
_{o}in Cable Structure apply and obtain the current initial mechanical calculating benchmark model A that upgrades after temperature variation
^{i} _{o}, upgrade A
^{i} _{o}time, U
^{i} _{o}all elements numerical value is also used U
^{i}all elements numerical value is corresponding to be replaced, and has upgraded U
^{i} _{o}, T
^{i} _{o}all elements numerical value is also used T
^{i}corresponding replacement of all elements numerical value, upgraded T
^{i} _{o}, so just obtained correctly corresponding to A
^{i} _{o}u
^{i} _{o}and T
^{i} _{o}; D now
^{i} _{o}remain unchanged.When upgrading A
^{i} _{o}after, A
^{i} _{o}the current initial damage of the cable system vector d for health status of rope
^{i} _{o}represent A
^{i} _{o}current cable structure steady temperature data vector T for Cable Structure steady temperature
^{i}represent A
^{i} _{o}current initial Cable Structure bearing angular coordinate vector U for bearing angular coordinate
^{i} _{o}represent, by Mechanics Calculation, obtain A
^{i} _{o}in concrete numerical value all monitored amounts, current, with these concrete numerical value, replace C
^{i} _{o}the element of middle correspondence, has so just realized the current initial value vector of monitored amount C
^{i} _{o}renewal.
The 4th step: circulation time must first be set up " unit damage monitored numerical quantity transformation matrices " and " evaluation object unit change vector " each time, and " unit damage monitored numerical quantity transformation matrices " that the i time circulation set up is designated as Δ C
^{i}(i=1,2,3 ...)." evaluation object unit change vector " that the i time circulation set up is designated as D
^{i} _{u}.Δ C in circulation each time
^{i}and D
^{i} _{u}need to according to circumstances constantly update, upgrade current initial mechanical calculating benchmark model A
^{i} _{o}, current initial Cable Structure steady temperature data vector T
^{i} _{o}with the current initial value vector of monitored amount C
^{i} _{o}after, upgrade unit damage monitored numerical quantity transformation matrices Δ C
^{i}with evaluation object unit change vector D
^{i} _{u}.
When circulation starts each time, first set up in the steps below unit damage monitored numerical quantity transformation matrices Δ Ci and evaluation object unit change vector D
^{i} _{u}; If upgraded A in the 3rd step
^{i} _{o}, in this step, must reestablish so (upgrading) unit damage monitored numerical quantity transformation matrices Δ C
^{i}with evaluation object unit change vector D
^{i} _{u}; If do not upgrade A in the 3rd step
^{i} _{o}, in this step, needn't reestablish so unit damage monitored numerical quantity transformation matrices Δ C
^{i}with evaluation object unit change vector D
^{i} _{u}; Set up and reestablish (upgrading) Δ C
^{i}and D
^{i} _{u}detailed process identical, be listed as follows:
Current initial mechanical calculating benchmark model A in Cable Structure
^{i} _{o}basis on carry out several times calculating, on calculation times numerical value, equal the quantity of all ropes.Calculating each time hypothesis only has an evaluation object (original centrepoint load variable quantity can be 0 at original damage or centrepoint load variable quantity, also can not be 0) basis on increase again unit damage or centrepoint load unit change, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable increases unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage) again, if this evaluation object is a centrepoint load, just suppose that this centrepoint load is at vectorial d
^{i} _{o}(if this centrepoint load is couple, centrepoint load unit change can be got 1kNm, 2kNm, 3kNm etc. for unit change on the basis of the existing variable quantity of this centrepoint load representing, to increase centrepoint load unit change again; If this centrepoint load is concentrated force, centrepoint load unit change can be got 1kN, 2kN, 3kN etc. for unit change).For convenience of calculating, while set increasing unit damage or centrepoint load unit change in circulation each time, can be all structural health conditions when this circulation started as being completely healthy, and set on this basis unit damage or centrepoint load unit change (in subsequent step, damage numerical value that calculate, evaluation object or centrepoint load variable quantitybe called name damage d
^{i} _{c}(i=1,2,3 ...), all when this circulation started, by the health status of evaluation object as being completely healthy speech, the formula therefore must foundation hereinafter providing damages the name calculating to be converted into true damage).The evaluation object that occurs unit damage or centrepoint load unit change in calculating each time with once circulation is different from the evaluation object that occurs unit damage or centrepoint load unit change in other calculating, and suppose each time the unit damage value of the evaluation object that has unit damage or centrepoint load unit change or unit damage value or the centrepoint load unit change numerical value that centrepoint load unit change numerical value can be different from other evaluation objects, with " evaluation object unit change vector D
^{i} _{u}" (as the formula (5)) record unit damage or the centrepoint load unit change of supposition of all evaluation objects in each circulation, circulation time is designated as D for the first time
^{1} _{u}calculate each time all utilize mechanics method (for example finite element method) calculate Cable Structure, at the current calculated value of the M of appointment monitored amount above, the current calculated value that calculates each time gained M monitored amount forms one " monitored amount calculation current vector ", and (when k evaluation object of hypothesis has unit damage, available formula (6) represents the monitored amount calculation current vector C of M monitored amount of all appointments
^{i} _{tk}); The monitored amount calculation current vector calculating each time deducts the current initial value vector of monitored amount C
^{i} _{o}, gained vector is exactly that " the numerical value change vector of monitored amount " of (take the position of the rope that has unit damage or numbering etc. are mark) (when k evaluation object has unit damage, used δ C under this condition
^{i} _{k}the numerical value change vector that represents monitored amount, δ C
^{i} _{k}definition see formula (7), formula (8) and formula (9), formula (7) deducts after formula (4) again divided by vectorial D for formula (6)
^{i} _{u}k element D
^{i} _{uk}gained), the numerical value change of monitored amount vector δ C
^{i} _{k}the supposition owing to calculating of each element representation have the unit damage of that evaluation object (for example k evaluation object) of unit damage or centrepoint load unit change or centrepoint load unit change (D for example
^{i} _{uk}), and the numerical value change amount of the corresponding monitored amount of this element causing is with respect to unit damage or the centrepoint load unit change numerical value D of supposition
^{i} _{uk}rate of change; There is N evaluation object just to have N " the numerical value change vector of monitored amount ", the numerical value change vector of each monitored amount has M element, by this N " the numerical value change vector of monitored amount ", forms successively " the unit damage monitored numerical quantity transformation matrices Δ C that has M * N element
^{i}" (the capable N row of M), each vectorial δ C
^{i} _{k}(k=1,2,3 ...., N) be matrix Δ C
^{i}row, Δ C
^{i}definition as the formula (10).
Evaluation object unit change vector D in formula (5)
^{i} _{u}element D
^{i} _{uk}(i=1,2,3, K=1,2,3 ...., N) represent unit damage or the centrepoint load unit change numerical value of k evaluation object of supposition in the i time circulation, vectorial D
^{i} _{u}in the numerical value of each element can be the same or different.
Elements C in formula (6)
^{i} _{tkj}(i=1,2,3 ...; K=1,2,3 ...., N; J=1,2,3 ...., M) represent that the i time circulation is while having unit damage or centrepoint load unit change due to k evaluation object, according to the calculating current value of the monitored amount of corresponding j the appointment of coding rule.
The subscript i(i=1 of each amount in formula (7), 2,3 ...) represent the i time circulation, subscript k(k=1,2,3 ...., N) represent unit damage or the centrepoint load unit change that k evaluation object increases, D in formula
^{i} _{uk}vectorial D
^{i} _{u}in k element.Vector δ C
^{i} _{k}definition suc as formula shown in (7) and formula (8), δ C
^{i} _{k}j(j=1,2,3 ...., M) individual element δ C
^{i}k
_{j}(definition as the formula (9)) represents, in the i time circulation, to set up matrix Δ C
^{i}time, suppose and when k evaluation object has unit damage or centrepoint load unit change, calculate the change amount of gained j monitored amount with respect to unit damage or the centrepoint load unit change D of supposition
^{i} _{uk}rate of change.
Vectorial δ C in formula (10)
^{i} _{k}(i=1,2,3 ....,, k=1,2,3 ...., N) represent in the i time circulation, because k evaluation object increases unit damage or centrepoint load unit change D
^{i} _{uk}cause, the relative value of all monitored amounts changes.Matrix Δ C
^{i}the coding rule and vectorial d above of row (subscript k)
^{i} _{o}the coding rule of subscript k of element identical.
The 5th step: the current health status of identification Cable Structure.Detailed process is as follows.
I(i=1,2,3 ...) in inferior circulation, utilize " the monitored amount current value vector C obtaining in second step actual measurement
^{i}" " the current initial value of monitored amount vector C together
^{i} _{o}", " unit damage monitored numerical quantity transformation matrices Δ C
^{i}" and " the vectorial d of current name damage
^{i} _{c}" between linear approximate relationship, shown in (11) or formula (12).
Monitored amount current value vector C in formula (11) and formula (12)
^{i}definition be similar to the current initial value of monitored amount vector C
^{i} _{o}definition, see formula (13); The vectorial d of the current name damage of evaluation object
^{i} _{c}definition see formula (14).
Elements C in formula (13)
^{i} _{j}(i=1,2,3 ....; J=1,2,3 ...., M) be the i time circulation time Cable Structure, according to the current value of the monitored amount of the corresponding j of being numbered of coding rule.
D in formula (14)
^{i} _{ck}(i=1,2,3 ....; K=1,2,3 ...., N) be current name damage or the current nominal centrepoint load changing value of k evaluation object in the i time circulation, vectorial d
^{i} _{c}coding rule and the matrix Δ C of subscript k of element
^{i}the coding rule of row identical.
When support cable actual damage is not too large, because Cable Structure material is still in the linear elasticity stage, the distortion of Cable Structure is also less, and the represented a kind of like this linear relationship of formula (11) or formula (12) is less with the error of actual conditions, and error can be used error vector e
^{i}(formula (15)) definition, the error of linear relationship shown in expression (11) or formula (12).
In formula (15), abs () is the function that takes absolute value, and each vectorial element of trying to achieve in bracket is taken absolute value.
Owing to there is certain error in formula (11) or the represented linear relationship of formula (12), therefore can not be simply according to formula (11) or formula (12) and " the vectorial C of monitored amount current value
^{i}" come direct solution to obtain the vectorial d of the current name damage of evaluation object
^{i} _{c}.If done like this, the vectorial d of the current name damage of evaluation object obtaining
^{i} _{c}in element even there will be larger negative value, namely negative damage, this is obviously irrational.Therefore obtain the vectorial d of the current name damage of evaluation object
^{i} _{c}acceptable solution (with reasonable error, but can be more exactly from cable system, determine damaged cable position and degree of injury thereof, also can determine more exactly centrepoint load variation numerical value) become a rational solution, available formula (16) is expressed this method.
In formula (16), abs () is the function that takes absolute value, vectorial g
^{i}description departs from the legitimate skew of ideal linearity relation (formula (11) or formula (12)), by formula (17), is defined.
G in formula (17)
^{i} _{j}(i=1,2,3 ....; J=1,2,3 ...., M) maximum allowable offset that departs from the ideal linearity relation shown in formula (11) or formula (12) in the i time circulation has been described.Vector g
^{i}can be according to the error vector e of formula (15) definition
^{i}tentative calculation is selected.
At the current initial value vector of monitored amount C
^{i} _{o}, unit damage monitored numerical quantity transformation matrices Δ C
^{i}with monitored amount current value vector C
^{i}when known, can utilize suitable algorithm (for example multiobjective optimization algorithm) to solve formula (16), obtain the vectorial d of the current name damage of evaluation object
^{i} _{c}acceptable solution, current actual damage vector d
^{i}the element of (formula (18) is shown in definition) can calculate according to formula (19), thereby can be by d
^{i}determine the health status of evaluation object.
D in formula (18)
^{i} _{k}(i=1,2,3, K=1,2,3 ...., N) representing the current actual health status of k evaluation object in the i time circulation, formula (19) is shown in its definition, if this evaluation object is the support cable (or pull bar) in cable system, so d
^{i} _{k}represent its current actual damage, d
^{i} _{k}be to represent not damaged at 0 o'clock, while being 100%, represent that this support cable thoroughly loses loadbearing capacity, in the time of between 0 and 100%, represent to lose the loadbearing capacity of corresponding proportion; If this evaluation object is centrepoint load, so a d
^{i} _{k}the current actual change numerical value that represents the centrepoint load that it is corresponding, vectorial d
^{i}the coding rule of element and formula (1) in vectorial d
_{o}the coding rule of element identical.
D in formula (19)
^{i} _{ok}(i=1,2,3,4, K=1,2,3 ...., N) be the current initial damage vector of evaluation object d
^{i} _{o}k element, d
^{i} _{ck}the vectorial d of the current name damage of evaluation object
^{i} _{c}k element.
Narration has below obtained the current actual damage vector of evaluation object d
^{i}after, the position of how to confirm slack line and relax level.
Total M in front known Cable Structure
_{1}root support cable, Cable Structure rope force data is by M
_{1}the Suo Li of root support cable describes.Available " initial rope force vector F
_{o}" represent the initial Suo Li (formula (20) is shown in definition) of all support cables in Cable Structure.
F in formula (20)
_{o}(h=1,2,3 ...., M
_{1}) be the initial Suo Li of h root support cable in Cable Structure, this element is the Suo Li corresponding to appointment support cable according to coding rule.Vector F
_{o}it is constant.In actual measurement, obtain T
_{o}synchronization, use conventional method directly to measure the rope force data that calculates all support cables, all these rope force datas form initial rope force vector F
_{o}.Setting up the initial mechanical calculating benchmark model A of Cable Structure
_{o}time in fact used vectorial F
_{o}.
By the current actual damage vector of evaluation object d
^{i}in the M relevant to support cable
_{1}individual element takes out, and forms the current actual damage vector of support cable d
^{ci}, the current actual damage vector of support cable d
^{ci}the coding rule and initial rope force vector F of element
_{o}the coding rule of element identical.The current actual damage vector of support cable d
^{ci}h element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., M
_{1}; The current actual damage vector of support cable d
^{ci}middle numerical value be not 0 element corresponding to the support cable of unsoundness problem, the support cable of these unsoundness problems is carried out to NonDestructive Testing, through NonDestructive Testing, find out that this support cable is not after damage, this element numerical value (is used d so
^{ci} _{h}expression) represent this support cable and d
^{ci} _{h}relaxing of impairment value mechanics equivalence, has just determined slack line thus, and the computing method of concrete slack illustrate below.
In this method, use " current cable force vector F
^{i}" represent the current cable power (formula (21) is shown in definition) of all support cables in Cable Structure that actual measurement obtains.
F in formula (21)
^{i} _{h}(h=1,2,3 ...., M
_{1}) be the current cable power of h root support cable in Cable Structure.In actual measurement, obtain current cable structure steady temperature data vector T
^{i}synchronization, actual measurement obtains the rope force data of all support cables in Cable Structure, all these rope force datas form current cable force vector F
^{i}.Vector F
^{i}element and vectorial F
_{o}the coding rule of element identical.According to narration above, vector T
^{i} _{o}equal vector T
^{i}.
In this method, under support cable original state, at the initial Cable Structure steady temperature data vector T for steady temperature data of Cable Structure
_{o}during expression, and support cable is when free state (free state refers to that Suo Li is 0, rear with), and the length of support cable is called initial drift, with " initial drift vector l
_{o}" represent the initial drift (formula (22) is shown in definition) of all support cables in Cable Structure.According to " the temperature survey calculating method of the Cable Structure of this method ", pass through vector T
_{o}can determine and obtain vector T
_{o}the Temperature Distribution of all support cables constantly.
L in formula (22)
_{oh}(h=1,2,3 ...., M
_{1}) be the initial drift of h root support cable in Cable Structure.Vector l
_{o}be constant, after determining when starting, just no longer change.
Similarly, under support cable original state, at the initial Cable Structure steady temperature data vector T for steady temperature data of Cable Structure
_{o}during expression, and support cable is when free state, and the crosssectional area of support cable is called initial free crosssectional area, with " initial free crosssectional area vector A
_{o}" representing the initial free crosssectional area (formula (23) is shown in definition) of all support cables in Cable Structure, the weight of the unit length of support cable is called the weight of initial free unit length, with " the weight vector ω of initial free unit length
_{o}" represent the weight (formula (24) is shown in definition) of the initial free unit length of all support cables in Cable Structure.
A in formula (23)
_{oh}(h=1,2,3 ...., M
_{1}) be the initial free crosssectional area of h root support cable in Cable Structure.Vector A
_{o}be constant, after determining when starting, just no longer change.
ω in formula (24)
_{oh}(h=1,2,3 ...., M1) be the weight of the free unit length of initial freedom of h root support cable in Cable Structure.Vector ω
_{o}be constant, after determining when starting, just no longer change.
In this method, at the current initial Cable Structure steady temperature data vector T for steady temperature data of Cable Structure
^{i} _{o}during expression, with " current initial drift vector l
^{i} _{o}" represent all support cables in Cable Structure current initial drift (formula (25) is shown in definition, refers to that hypothesis supporting cable force is at 0 o'clock, after having considered that thermal expansivity and temperature variation are on the impact of support cable drift, initial drift vector l
_{o}with initial Cable Structure steady temperature data vector T
_{o}the support cable representing is at current initial Cable Structure steady temperature data vector T for temperature
^{i} _{o}support cable drift during expression).According to " the temperature survey calculating method of the Cable Structure of this method ", pass through vector T
^{i} _{o}can determine and obtain vector T
^{i} _{o}the Temperature Distribution of all support cables constantly.
L in formula (25)
^{i} _{oh}(h=1,2,3 ...., M
_{1}) be the current initial Cable Structure steady temperature data vector T for steady temperature data in Cable Structure
^{i} _{o}during expression, the current initial drift of h root support cable in Cable Structure, can utilize thermal expansivity, the l of support cable
_{oh}, T
_{o}and T
^{i} _{o}by conventional physical computing, obtain l
^{i} _{oh}, now impossible this support cable is impaired and lax.
Vector d
^{ci}element, the element of vectorial F, vectorial l
_{o}element, vectorial l
^{i} _{o}element, vectorial A
_{o}element, vectorial ω
_{o}element and vectorial F
_{o}the coding rule of element identical, the different information of the same support cable of element representation of the identical numbering of these vectors.
In this method, at the current initial Cable Structure steady temperature data vector T for steady temperature data of Cable Structure
^{i} _{o}during expression, with " current drift vector l
^{i}" represent the current drift (formula (26) is shown in definition, and now support cable may be intact, may be also impaired, also may relax) of all support cables in Cable Structure.
L in formula (26)
^{i} _{h}(h=1,2,3 ...., M
_{1}) be the current drift of h root support cable in Cable Structure.
In this method, with " drift changes vectorial Δ l
^{i}" (or claiming support cable current relax level vector) represent the change amount (formula (27) and formula (28) are shown in definition) of the drift of all support cables in Cable Structure.
Δ l in formula (27)
^{i} _{h}(h=1,2,3 ...., M
_{1}) be the change amount of the drift of h root support cable in current cable structure, formula (28), Δ l are shown in its definition
^{i} _{h}be not that 0 rope is slack line, Δ l
^{i} _{h}the numerical value slack that is rope, and represent the current relax level of cable system h root support cable, be also the long adjustment amount of rope of this rope while adjusting Suo Li.
By slack line is carried out to mechanics equivalence with damaged cable, carry out the relax level of slack line and identify in the method, mechanics equivalent condition is:
One, the mechanics parameters without lax initial drift, geometrical property parameter and material during with not damaged of the rope of two equivalences is identical;
Two, after lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after overall length identical.
While meeting abovementioned two mechanics equivalent conditions, the mechanics function of two such support cables in structure is exactly identical, if replaced after slack line with equivalent damaged cable, Cable Structure any variation can not occur, and vice versa.
Obtained the current actual damage vector of support cable d
^{ci}after, d
^{ci}h element d
^{ci} _{h}(h=1,2,3 ...., M
_{1}) represent the actual damage value of h root support cable, although by d
^{ci} _{h}be called the actual damage value of h root rope or the actual damage degree of h root rope, but due to h root support cable may be impaired may be also lax, so d
^{ci}h element d
^{ci} _{h}the actual damage value of the h root support cable representing is actually the actual equivalent damage value of h root support cable, when h root support cable is actually impaired, and d
^{ci} _{h}the actual damage value of the h root support cable just representing, when h root support cable is actually lax, d
^{ci} _{h}h root support cable and lax equivalent actual damage value with regard to representing, for sake of convenience, claims d in the method
^{ci} _{h}be to represent h root support cable not damaged at 0 o'clock, while being 100%, represent that this rope thoroughly loses loadbearing capacity, in the time of between 0 and 100%, represent that h root support cable loses the loadbearing capacity of corresponding proportion, by the current actual damage vector of support cable d
^{ci}just can identify the support cable that health status goes wrong, but in the support cable that these health status go wrong, some is impaired, some is to have relaxed, if h support cable is actually (its current relax level Δ l that relaxed
^{i} _{h}the current relax level Δ l of h support cable definition), relaxing so
^{i} _{h}(Δ l
^{i} _{h}definition see formula (27)) with the current actual damage degree d of equivalent damaged cable
^{ci} _{h}between relation by aforementioned two mechanics equivalent conditions, determined.Δ l
^{i} _{h}same d
^{ci} _{h}between physical relationship can adopt accomplished in many ways, for example can directly according to aforementioned equivalent condition, determine (referring to formula (29)), also can adopt based on Ernst equivalent elastic modulus and replace the E in formula (29) to revise rear determine (referring to formula (30)), also can adopt other methods such as trial and error procedure based on finite element method to determine.
E in formula (29) and formula (30)
^{i} _{h}the current initial Cable Structure steady temperature data vector T for steady temperature data in Cable Structure
^{i} _{o}during expression, the elastic modulus of h root support cable, A
^{i} _{h}the current initial Cable Structure steady temperature data vector T for steady temperature data in Cable Structure
^{i} _{o}during expression, the crosssectional area of h root support cable, F
^{i} _{h}the current initial Cable Structure steady temperature data vector T for steady temperature data in Cable Structure
^{i} _{o}during expression, the current cable power of h root support cable, d
^{ci} _{h}the current actual damage degree of h root support cable, ω
^{i} _{h}the current initial Cable Structure steady temperature data vector T for steady temperature data in Cable Structure
^{i} _{o}during expression, the weight of the unit length of h root support cable, l
^{i} _{xh}the current initial Cable Structure steady temperature data vector T for steady temperature data in Cable Structure
^{i} _{o}during expression, the horizontal range of two supporting end points of h root support cable, l
^{i} _{xh}current support cable two supporting end points horizontal range vector l
^{i} _{x}an element, current support cable two supporting end points horizontal range vector l
^{i} _{x}the coding rule and initial drift vector l of element
_{o}the coding rule of element identical, E
^{i} _{h}can obtain according to the characteristic material data of looking into or survey h root support cable A
^{i} _{h}and ω
^{i} _{h}can be according to the thermal expansivity of h root support cable, A
_{oh}, ω
_{oh}, F
^{i} _{h}, T
_{o}and T
^{i} _{o}by conventional physics and Mechanics Calculation, obtain.Item in formula (30) in [] is the Ernst equivalent elastic modulus of this support cable, by formula (29) or formula (30), can just can determine the current relax level vector of support cable Δ l
^{i}.Formula (30) is the correction to formula (29).
In a word, when bearing has angular displacement, this method has so far realized two kinds of functions that existing method can not possess, be respectively, during the centrepoint load of one, bearing in structure and structure (environment) temperature variation, can reject angular displacement of support, centrepoint load variation and structure temperature and change the impact on Cable Structure health status recognition result, thereby identify exactly the structure health monitoring method of problem rope; Two, this method, when identifying problem rope, can also identify the variation of centrepoint load simultaneously, and this method can be rejected angular displacement of support, structure temperature changes and the impact of support cable health status variation, realizes the correct identification of centrepoint load intensity of variation.
The 6th step: judge whether to finish this (the i time) circulation, if so, complete this circulation and finish front tailing in work, for next time (the i+1 time, i=1,2,3,4 ...) circulation preparation Mechanics Calculation benchmark model and necessary vector.Detailed process is as follows:
In this (the i time) circulation, try to achieve the vectorial d of current name damage
^{i} _{c}after, first, according to formula (31), set up mark vector B
^{i}, formula (32) has provided mark vector B
^{i}the definition of k element; If mark vector B
^{i}element be 0 entirely, get back to second step and proceed the health monitoring of Cable Structure and calculating; If mark vector B
^{i}element be not 0 entirely, complete after subsequent step, enter next time circulation.
Socalled subsequent step is: first, according to formula (33) calculate next time (the i+1 time, i=1,2,3,4 ...) circulate initial damage vector d required
^{i+1} _{o}each element d
^{i+1} _{ok}; The second, at Mechanics Calculation benchmark model A
_{o}basis on, make A
_{o}in the health status of evaluation object be d
^{i+1} _{o}rather than be d
_{o}after, more further to A
_{o}in Cable Structure apply temperature variation (as previously mentioned, the numerical value of the temperature variation applying just taken from steady temperature change vector S, and steady temperature change vector S equals T
^{i}deduct T
_{o}), so just obtained next time (the i+1 time, i=1,2,3,4 ...) required current initial mechanical calculating benchmark mould A circulates
^{i+1} _{o}, next time (the i+1 time, i=1,2,3,4 ...) required current initial Cable Structure steady temperature data vector T circulates
^{i+1} _{o}equal T
^{i} _{o}, more further to A
_{o}in Cable Structure apply angular displacement of support constraint (as previously mentioned, the numerical value of the angular displacement of support constraint applying just taken from angular displacement of support vector V, and angular displacement of support vector V equals U
^{i}deduct U
_{o}), so just obtained next time (the i+1 time, i=1,2,3,4 ...) required current initial mechanical calculating benchmark mould A circulates
^{i+1} _{o}, next time (the i+1 time, i=1,2,3,4 ...) required current initial Cable Structure bearing angular coordinate vector U circulates
^{i+1} _{o}equal U
^{i} _{o}, to A
^{i+1} _{o}carrying out Mechanics Calculation obtains corresponding to A
^{i+1} _{o}concrete numerical value all monitored amounts, current, these concrete numerical value form next time (the i+1 time, i=1,2,3,4 ...) the current initial value vector C of the required monitored amount that circulates
^{i+1} _{o}.
Mark vector B in formula (31)
^{i}subscript i represent the i time circulation, its element B
^{i} _{k}(k=1,2,3 ..., subscript k N) represents the health status feature of k evaluation object, can only get 0 and 1 two amount, concrete value rule is shown in formula (31).
Element B in formula (32)
^{i} _{k}k the element of mark vector Bi, D
^{i} _{uk}evaluation object unit change vector D
^{i} _{u}k element (seeing formula (5)), d
^{i} _{ck}the vectorial d of the current name damage of evaluation object
^{i} _{c}k element (seeing formula (14)), they all represent the relevant information of k evaluation object.
D in formula (33)
^{i} _{uk}evaluation object unit change vector D
^{i} _{u}k element (seeing formula (5)), d
^{i} _{ok}the current initial damage vector of evaluation object d
^{i} _{o}k element (seeing formula (3)).
The third part of this method: the software and hardware part of health monitoring systems.
Hardware components comprises monitoring system (the space coordinate monitoring system that comprises the supporting end points of monitored amount monitoring system, temperature monitoring system, Cable Structure bearing angular coordinate monitoring system, support cable cable force monitoring system, support cable), signal picker and computing machine etc.Require RealTime Monitoring to obtain the measured data of volume coordinate of the supporting end points of temperature required, Cable Structure bearing angular coordinate, supporting cable force and support cable, require each monitored amount of RealTime Monitoring simultaneously.
Software section should complete the process that this method sets, complete needed in this method, can be by functions such as computer implemented monitoring, record, control, storage, calculating, notice, warnings.
This method specifically comprises:
A. for sake of convenience, it is evaluation object that this method unitedly calls evaluated support cable and centrepoint load, establishes the quantity of evaluated support cable and the quantity sum of centrepoint load is N, and the quantity of evaluation object is N; Determine the coding rule of evaluation object, by this rule, by evaluation object numberings all in Cable Structure, this numbering will be for generating vector sum matrix in subsequent step; This method represents this numbering with variable k, k=1, and 2,3 ..., N; Appointment by the support cable of monitored Suo Li while determining hybrid monitoring, if total Q root support cable in cable system, the monitored rope force data of Cable Structure specifies M1 rope force data of support cable to describe by M1 in Cable Structure, and the variation of Cable Structure Suo Li is exactly the variation of the Suo Li of all appointment support cables; Each total M
_{1}individual cable force measurement value or calculated value characterize the rope force information of Cable Structure; M
_{1}be one and be not less than 0 integer that is not more than Q; While determining hybrid monitoring appointment by the measured point of monitored strain, the monitored strain data of Cable Structure can be by K in Cable Structure
_{2}l individual specified point and each specified point
_{2}the strain of individual assigned direction is described, and the variation of Cable Structure strain data is exactly K
_{2}the variation of all tested strains of individual specified point; Each total M
_{2}individual strain measurement value or calculated value characterize Cable Structure strain, M
_{2}for K
_{2}and L
_{2}longpending; M
_{2}to be not less than 0 integer; While determining hybrid monitoring appointment by the measured point of monitored angle, the monitored angledata of Cable Structure is by K in Cable Structure
_{3}l individual specified point, that cross each specified point
_{3}h individual appointment straight line, each appointment straight line
_{3}individual angle coordinate component is described, and the variation of Cable Structure angle is exactly the variation of angle coordinate components appointment straight lines all specified points, all, all appointments; Each total M
_{3}individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M
_{3}for K
_{3}, L
_{3}and H
_{3}longpending; M
_{3}be one and be not less than 0 integer; While determining hybrid monitoring appointment by monitored shape data, the monitored shape data of Cable Structure is by K in Cable Structure
_{4}l individual specified point and each specified point
_{4}the volume coordinate of individual assigned direction is described, and the variation of Cable Structure shape data is exactly K
_{4}the variation of all coordinate components of individual specified point; Each total M
_{4}individual measurement of coordinates value or calculated value characterize Cable Structure shape, M
_{4}for K
_{4}and L
_{4}longpending; M
_{4}be one and be not less than 0 integer; The monitored amount of comprehensive abovementioned hybrid monitoring, total M the monitored amount of whole Cable Structure, M is M
_{1}, M
_{2}, M
_{3}and M
_{4}sum, definition parameter K, K is M
_{1}, K
_{2}, K
_{3}and K
_{4}sum, K and M must not be less than the quantity N of evaluation object; For simplicity, in the method this is walked to listed M monitored amount referred to as " monitored amount "; In this method, to the time interval between any twice measurement of same amount RealTime Monitoring, must not be greater than 30 minutes, the moment of survey record data is called physical record data constantly;
B. this method definition " the temperature survey calculating method of the Cable Structure of this method " is undertaken by step b1 to b3;
B1: inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of Cable Structure composition material and Cable Structure environment of living in, utilize the geometry measured data of design drawing, asconstructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure, inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, statistics obtains interior during this period of time cloudy quantity and is designated as T cloudy day, in the method can not be seen to one of the sun daytime and be called all day the cloudy day, statistics obtains 0 highest temperature and the lowest temperature between latter 30 minutes of the moment of sunrise next day at each cloudy day in T cloudy day, sunrise on the meteorology that sunrise refers to base area revolutions constantly and the rule that revolves round the sun is definite constantly, do not represent necessarily can see the sun same day, the sunrise that can inquire about data or calculate each required day by conventional meteorology constantly, each cloudy day 0 up to next day sunrise constantly the highest temperature between latter 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, there is T cloudy day, the maximum temperature difference that just has the daily temperature at T cloudy day, get maximal value in the maximum temperature difference of daily temperature at T cloudy day for reference to temperature difference per day, with reference to temperature difference per day, be designated as Δ T
_{r}, between inquiry Cable Structure location and Altitude Region, place, be no less than temperature that the meteorological data in recent years of 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and Altitude Region, place about the maximum rate of change Δ T of sea level elevation
_{h}, for convenience of narration, get Δ T
_{h}unit be ℃/m, on the surface of Cable Structure, get " R Cable Structure surface point ", get the Specific Principles of " R Cable Structure surface point " narrates in step b3, after will by actual measurement, obtain the temperature of this R Cable Structure surface point, claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, by Calculation of Heat Transfer, obtain the temperature of this R Cable Structure surface point, just claim that the temperature data calculating is " R Cable Structure surface temperature computational data ", from the residing minimum height above sea level of Cable Structure to the highest height above sea level, in Cable Structure, uniform choosing is no less than three different sea level elevations, the sea level elevation place choosing at each, at the intersection place on surface level and Cable Structure surface, at least choose two points, outer normal from selected point straw line body structure surface, all outer normal directions of choosing are called " measuring Cable Structure along the direction of the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in the measurement Cable Structure of choosing along comprising the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure in the direction of the Temperature Distribution of wall thickness, direction uniform choosing in Cable Structure along each measurement Cable Structure along the Temperature Distribution of wall thickness is no less than three points, measure all temperature that are selected a little, the temperature recording is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure Cable Structure along the direction of the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if chosen H different sea level elevation, at each sea level elevation place, choose B and measured Cable Structure along the direction of the Temperature Distribution of wall thickness, along each, measure Cable Structure and in Cable Structure, chosen E point along the direction of the Temperature Distribution of wall thickness, wherein H and E are not less than 3, B is not less than 2, if HBE is the product of H and B and E, corresponding total HBE " measuring Cable Structure along the point of the temperature profile data of thickness ", after will obtain the temperature that this HBE " measures Cable Structure along the point of the temperature profile data of thickness " by actual measurement, claim that the temperature data that actual measurement obtains is " HBE Cable Structure is along thickness temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, by Calculation of Heat Transfer, obtain this HBE and measure Cable Structure along the temperature of the point of the temperature profile data of thickness, just claim that the temperature data calculating is " HBE Cable Structure is along thickness temperature computation data ", if BE is the product of B and E, total BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in sea level elevation place of choosing at each in this method, in Cable Structure location, according to meteorology, measure temperature and require to choose a position, will obtain meeting the temperature that meteorology is measured the Cable Structure place environment of temperature requirement in this position actual measurement, in the onsite spaciousness of Cable Structure, unobstructed place chooses a position, this position should can obtain in each day of the whole year this ground the most sufficient sunshine of getable this day, flat board at a carbon steel material of this position of sound production, be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, the one side of this reference plate on the sunny side, be called sunny slope, the sunny slope of reference plate is coarse and dark color, the sunny slope of reference plate should can obtain in each day of the whole year one flat plate on this ground the most sufficient sunshine of getable this day, the nonsunny slope of reference plate is covered with insulation material, RealTime Monitoring is obtained to the temperature of the sunny slope of reference plate,
B2: RealTime Monitoring obtains R Cable Structure surface temperature measured data of abovementioned R Cable Structure surface point, RealTime Monitoring obtains previously defined Cable Structure along the temperature profile data of thickness simultaneously, and RealTime Monitoring obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement simultaneously, by RealTime Monitoring, obtain being carved at sunrise the same day sunrise next day temperature measured data sequence of the Cable Structure place environment between latter 30 minutes constantly, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data that was carved at sunrise the Cable Structure place environment between latter 30 minutes of the moment of sunrise next day the same day, find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment, by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, deduct and be carved at sunrise the sunrise next day maximum temperature difference between latter 30 minutes constantly on same day that minimum temperature obtains Cable Structure place environment, be called environment maximum temperature difference, be designated as Δ T
_{emax}, by the temperature measured data sequence of Cable Structure place environment, by conventional mathematical computations, obtain the temperature of Cable Structure place environment about the rate of change of time, this rate of change is also along with the time changes, by RealTime Monitoring, obtain being carved at sunrise the same day sunrise next day measured data sequence of the temperature of the sunny slope of the reference plate between latter 30 minutes constantly, the measured data sequence of the temperature of the sunny slope of reference plate by be carved at sunrise the same day next day sunrise constantly the measured data of the temperature of the sunny slope of the reference plate between latter 30 minutes according to time order and function order, arrange, find maximum temperature and minimum temperature in the measured data sequence of temperature of sunny slope of reference plate, by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, deduct and be carved at sunrise the sunrise next day maximum temperature difference between latter 30 minutes constantly on same day of temperature that minimum temperature obtains the sunny slope of reference plate, be called reference plate maximum temperature difference, be designated as Δ T
_{pmax}, by RealTime Monitoring, obtain being carved at sunrise the same day sunrise next day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between latter 30 minutes constantly, there is R Cable Structure surface point just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence by be carved at sunrise on same day of a Cable Structure surface point sunrise next day constantly the Cable Structure surface temperature measured data between latter 30 minutes according to time order and function order, arrange, find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence, by the maximum temperature in each Cable Structure surface temperature measured data sequence, deduct and be carved at sunrise the sunrise next day maximum temperature difference between latter 30 minutes constantly on same day that minimum temperature obtains the temperature of each Cable Structure surface point, there is R Cable Structure surface point just to have and be carved at sunrise the sunrise next day maximum temperature difference numerical value between latter 30 minutes constantly R the same day, maximal value is wherein called Cable Structure surface maximum temperature difference, be designated as Δ T
_{smax}, by each Cable Structure surface temperature measured data sequence, by conventional mathematical computations, obtain the temperature of each Cable Structure surface point about the rate of change of time, the temperature of each Cable Structure surface point about the rate of change of time also along with the time changes, by RealTime Monitoring, obtain being carved at sunrise the same day between latter 30 minutes of the moment of sunrise next day, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", the sea level elevation place that calculating is chosen at each amounts to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and the difference of minimum temperature, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", chosen H different sea level elevation and just had H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T
_{tmax},
B3: measure and calculate acquisition Cable Structure steady temperature data, first, determine the moment that obtains Cable Structure steady temperature data, the condition relevant to the moment that determines acquisition Cable Structure steady temperature data has six, first condition is the moment that obtains Cable Structure steady temperature data to be carved at sunset sunrise next day constantly between latter 30 minutes between the same day, sunset refers to sunset on base area revolutions and the definite meteorology of revolution rule constantly constantly, and the sunset that can inquire about data or calculate each required day by conventional meteorology constantly, the a condition of second condition be the same day be carved at sunrise next day sunrise constantly between latter 30 minutes during this period of time in, reference plate maximum temperature difference Δ T
_{pmax}with Cable Structure surface maximum temperature difference Δ T
_{smax}all be not more than 5 degrees Celsius, the b condition of second condition be the same day be carved at sunrise next day sunrise constantly between latter 30 minutes during this period of time in, above, measure the environment maximum error Δ T calculating
_{emax}be not more than with reference to temperature difference per day Δ T
_{r}, and reference plate maximum temperature difference Δ T
_{pmax}after deducting 2 degrees Celsius, be not more than Δ T
_{emax}, and Cable Structure surface maximum temperature difference Δ T
_{smax}be not more than Δ T
_{pmax}, only needing to meet in a condition of second and b condition one is just called and meets second condition, the 3rd condition is that the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, the 4th condition is that the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, the 5th condition is in the moment that obtains Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is to be carved at sunrise the sunrise next day minimal value between latter 30 minutes constantly the same day, the 6th condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T
_{tmax}be not more than 1 degree Celsius, this method is utilized abovementioned six conditions, any one in following three kinds of moment is called to " obtaining the mathematics of Cable Structure steady temperature data constantly ", the first is first moment to the 5th condition meeting in abovementioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data " constantly, the second is the moment that only meets the 6th condition in abovementioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data " constantly, the third is first moment to the 6th condition simultaneously meeting in abovementioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data " constantly, when obtaining the mathematics of Cable Structure steady temperature data, be exactly in this method during in constantly one of physical record data constantly, the moment that obtains Cable Structure steady temperature data be exactly obtain Cable Structure steady temperature data mathematics constantly, if obtain the mathematics of Cable Structure steady temperature data and be not constantly any in constantly of physical record data in this method constantly, get this method close to moment of the mathematics that obtains Cable Structure steady temperature data those physical record data constantly for obtaining the moment of Cable Structure steady temperature data, this method is carried out the relevant health monitoring analysis of Cable Structure by using in the amount that obtains the moment survey record of Cable Structure steady temperature data, this method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, i.e. this Cable Structure temperature constantly temporal evolution not, and this is exactly " obtaining the moment of Cable Structure steady temperature data " of this method constantly, then, according to Cable Structure heat transfer characteristic, utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment that obtains Cable Structure steady temperature data, utilize the thermal conduction study computation model of Cable Structure, by conventional Calculation of Heat Transfer, obtain in the Temperature Distribution of Cable Structure that obtains the moment of Cable Structure steady temperature data, now calculate by stable state in the temperature field of Cable Structure, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data, also comprise that Cable Structure is in the accounting temperature of selected HBE " measuring Cable Structure along the point of the temperature profile data of thickness " above, the accounting temperature of HBE " measuring Cable Structure along the point of the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", when R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure stable state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", while getting " R Cable Structure surface point " on the surface of Cable Structure, the quantity of " R Cable Structure surface point " and necessary three conditions that meet that distribute, first condition is when Cable Structure temperature field is during in stable state, on Cable Structure surface the temperature of any point be by " R Cable Structure surface point " with Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point while obtaining, on the Cable Structure surface that linear interpolation obtains, on the temperature of this arbitrfary point and Cable Structure surface, the error of the actual temperature of this arbitrfary point is not more than 5%, Cable Structure surface comprises support cable surface, second condition is that in " R Cable Structure surface point ", the quantity at the point of same sea level elevation is not less than 4, and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 ℃ divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation
_{h}the numerical value obtaining, gets Δ T for convenience of narration
_{h}unit be ℃/m that the unit of getting Δ h for convenience of narration is m, " R Cable Structure surface point " refers to while only considering sea level elevation along the definition of adjacent Cable Structure surface point between two of sea level elevation, in " R Cable Structure surface point ", do not have a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two, the 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to the geometric properties of Cable Structure and bearing data, in Cable Structure, find and be subject to the sunshineduration position of those surface points the most fully the whole year, in " R Cable Structure surface point ", having a Cable Structure surface point at least is the annual point being subject in the most sufficient those surface points of sunshineduration in Cable Structure,
C. according to " the temperature survey calculating method of the Cable Structure of this method ", directly measure and calculate the Cable Structure steady temperature data under original state, Cable Structure steady temperature data under original state are called initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T
_{o}", actual measurement or consult reference materials and obtain the temperature variant physical and mechanical properties parameter of the various materials that Cable Structure used, in actual measurement, obtain initial Cable Structure steady temperature data vector T
_{o}synchronization, directly measure the initial Suo Li that calculates all support cables, form initial rope force vector F
_{o}, according to the data that comprise Cable Structure design data, completion data obtain that all support cables are in free state that Suo Li is the length of 0 o'clock, the weight of crosssectional area during in free state and the unit length during in free state, and the temperature of all support cables while obtaining these three kinds of data, utilize on this basis temperature variant physical function parameter and the mechanical property parameters of all support cables, according to conventional physical computing, obtain all support cables at initial Cable Structure steady temperature data vector T
_{o}suo Li under condition is that the length of 0 o'clock all support cable, crosssectional area and the Suo Li that Suo Li is 0 o'clock all support cable are the weight of the unit length of 0 o'clock all support cable, form successively the initial drift vector of support cable, the weight vector of the initial free unit length of initial free crosssectional area vector sum, the coding rule and initial rope force vector F of the element of the initial drift vector of support cable, the weight vector of the initial free unit length of initial free crosssectional area vector sum
_{o}the coding rule of element identical, in actual measurement, obtain T
_{o}time, namely obtaining initial Cable Structure steady temperature data vector T
_{o}the synchronization in the moment, directly measure the measured data that calculates initial Cable Structure, the measured data of initial Cable Structure is to comprise Cable Structure centrepoint load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the initial value of all monitored amounts, the initial rope force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing generalized coordinate data, initial Cable Structure angledata, initial Cable Structure spatial data is in interior measured data, initial Cable Structure bearing generalized coordinate data comprise initial Cable Structure bearing spatial data and initial Cable Structure bearing angular data, when obtaining the measured data of initial Cable Structure, measurement calculates the data of the health status that can express support cable of the Nondestructive Testing Data that comprises support cable, and the data of the health status that can express support cable are now called support cable initial health data, the initial value of all monitored amounts forms monitored amount initial value vector C
_{o}, monitored amount initial value vector C
_{o}the coding rule of coding rule and M monitored amount identical, utilize support cable initial health data and Cable Structure centrepoint load measurement data to set up evaluation object initial damage vector d
_{o}, vectorial d
_{o}represent with initial mechanical calculating benchmark model A
_{o}the initial health of the evaluation object of the Cable Structure representing, evaluation object initial damage vector d
_{o}element number equal N, d
_{o}element and evaluation object be onetoone relationship, vectorial d
_{o}the coding rule of element identical with the coding rule of evaluation object, if d
_{o}evaluation object corresponding to some elements be support cable, so a d in cable system
_{o}the numerical value of this element represent the initial damage degree of corresponding support cable, if the numerical value of this element is 0, represent that the corresponding support cable of this element is intact, do not damage, if its numerical value is 100%, represent that the corresponding support cable of this element has completely lost loadbearing capacity, if its numerical value between 0 and 100%, represents this support cable, lost the loadbearing capacity of corresponding proportion, if d
_{o}evaluation object corresponding to some elements be some centrepoint load, in this method, get d
_{o}this element numerical value be 0, the initial value that represents the variation of this centrepoint load is 0, if while there is no the Nondestructive Testing Data of support cable and the data of other health status that can express support cable, or can think that structure original state is not damaged during without relaxed state, vectorial d
_{o}in each element numerical value relevant to support cable get 0, initial Cable Structure bearing angular data forms initial Cable Structure bearing angular coordinate vector U
_{o},
The temperature variant physical and mechanical properties parameter of the various materials that d. use according to measured data, support cable initial health data, Cable Structure centrepoint load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the Cable Structure of the design drawing of Cable Structure, asconstructed drawing and initial Cable Structure, initial Cable Structure bearing angular coordinate vector U
_{o}, initial Cable Structure steady temperature data vector T
_{o}with all Cable Structure data that preceding step obtains, set up the initial mechanical calculating benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data "
_{o}, based on A
_{o}the Cable Structure computational data calculating must approach its measured data very much, and difference therebetween must not be greater than 5%; Corresponding to A
_{o}" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T
_{o}"; Corresponding to A
_{o}cable Structure bearing angular data be exactly initial Cable Structure bearing angular coordinate vector U
_{o}; Corresponding to A
_{o}evaluation object initial damage vector d for evaluation object health status
_{o}represent; Corresponding to A
_{o}monitored amount initial value vector C for the initial value of all monitored amounts
_{o}represent; U
_{o}, T
_{o}and d
_{o}a
_{o}parameter, by A
_{o}initial value and the C of all monitored amounts of obtaining of Mechanics Calculation result
_{o}the initial value of all monitored amounts that represent is identical, therefore also can say C
_{o}by A
_{o}mechanics Calculation result form, A in the method
_{o}, C
_{o}, d
_{o}, U
_{o}and T
_{o}constant;
E. in the method, alphabetical i, except representing that significantly, the place of number of steps, alphabetical i only represents cycle index, circulates for the i time; When the i time circulation starts, the current initial mechanical calculating benchmark model of Cable Structure that need to set up or that set up is designated as current initial mechanical calculating benchmark model A
^{i} _{o}, A
_{o}and A
^{i} _{o}count temperature parameter, can accounting temperature change the Effect on Mechanical Properties to Cable Structure; When the i time circulation starts, corresponding to A
^{i} _{o}" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T
^{i} _{o}represent vector T i
_{o}definition mode and vector T
_{o}definition mode identical, T
^{i} _{o}element and T
_{o}element corresponding one by one; When the i time circulation starts, corresponding to A
^{i} _{o}" Cable Structure bearing angular data " with current initial Cable Structure bearing angular coordinate vector U
^{i} _{o}represent vectorial U
^{i} _{o}definition mode and vectorial U
_{o}definition mode identical, U
^{i} _{o}element and U
_{o}element corresponding one by one; The current initial damage vector of evaluation object that the i time circulation needs while starting is designated as d
^{i} _{o}, d
^{i} _{o}cable Structure A while representing this circulation beginning
^{i} _{o}the health status of evaluation object, d
^{i} _{o}definition mode and d
_{o}definition mode identical, d
^{i} _{o}element and d
_{o}element corresponding one by one; When the i time circulation starts, the initial value of all monitored amounts, with the current initial value vector of monitored amount C
^{i} _{o}represent vectorial C
^{i} _{o}definition mode and vectorial C
_{o}definition mode identical, C
^{i} _{o}element and C
_{o}element corresponding one by one, the current initial value vector of monitored amount C
^{i} _{o}expression is corresponding to A
^{i} _{o}the concrete numerical value of all monitored amounts; U
^{i} _{o}, T
^{i} _{o}and d
^{i} _{o}a
^{i} _{o}characterisitic parameter, C
^{i} _{o}by A
^{i} _{o}mechanics Calculation result form; When circulation starts for the first time, A
^{i} _{o}be designated as A
^{1} _{o}, set up A
^{1} _{o}method for making A
^{1} _{o}equal A
_{o}; When circulation starts for the first time, T
^{i} _{o}be designated as T
^{1} _{o}, set up T
^{1} _{o}method for making T
^{1} _{o}equal T
_{o}; When circulation starts for the first time, U
^{i} _{o}be designated as U
^{1} _{o}, set up U
^{1} _{o}method for making U
^{1} _{o}equal U
_{o}; When circulation starts for the first time, d
^{i} _{o}be designated as d
^{1} _{o}, set up d
^{1} _{o}method for making d
^{1} _{o}equal d
_{o}; When circulation starts for the first time, C
^{i} _{o}be designated as C
^{1} _{o}, set up C
^{1} _{o}method for making C
^{1} _{o}equal C
_{o};
F. from entering the circulation that is walked q step by f here; In structure military service process, according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement, obtain the current data of Cable Structure steady temperature data, the current data of all " Cable Structure steady temperature data " forms current cable structure steady temperature data vector Ti, the definition mode of vector T i and vector T
_{o}definition mode identical, T
^{i}element and T
_{o}element corresponding one by one; In actual measurement, obtain current cable structure steady temperature data vector T
^{i}synchronization, actual measurement obtains Cable Structure bearing angular coordinate current data, all Cable Structure bearing angular coordinate current datas form current cable structures actual measurement bearing angular coordinate vector U
^{i}, vectorial U
^{i}definition mode and vectorial U
_{o}definition mode identical; In actual measurement, obtain vector T
^{i}time, actual measurement obtains obtaining current cable structure steady temperature data vector T
^{i}the Cable Structure of synchronization in the moment in the currency of all monitored amounts, all these numerical value form monitored amount current value vector C
^{i}, vectorial C
^{i}definition mode and vectorial C
_{o}definition mode identical, C
^{i}element and C
_{o}element corresponding one by one, represent that identical monitored amount is at numerical value in the same time not; In actual measurement, obtain the synchronization of current cable structure steady temperature data vector Ti, actual measurement obtains all M in Cable Structure
_{1}the rope force data of root support cable, all these rope force datas form current cable force vector F
^{i}, vectorial F
^{i}element and vectorial F
_{o}the coding rule of element identical; In actual measurement, obtain current cable structure steady temperature data vector T
^{i}synchronization, Actual measurement obtains all M
_{1}the volume coordinate of two supporting end points of root support cable, the volume coordinate of two the supporting end points in the horizontal direction difference of component is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables form current support cable two supporting end points horizontal range vectors, the coding rule and initial rope force vector F of the element of current support cable two supporting end points horizontal range vectors
_{o}the coding rule of element identical;
G. according to current cable structure actual measurement bearing angular coordinate vector U
^{i}with current cable structure steady temperature data vector T
^{i}, according to step g 1 to g3, upgrade current initial mechanical calculating benchmark model A
^{i} _{o}, the current initial value of monitored amount vector C
^{i} _{o}, current initial Cable Structure steady temperature data vector T
^{i} _{o}with current initial Cable Structure bearing angular coordinate vector U
^{i} _{o}, and the current initial damage vector of evaluation object d
^{i} _{o}remain unchanged;
G1. compare respectively T
^{i}and T
^{i} _{o}, U
^{i}and U
^{i} _{o}if, T
^{i}equal T
^{i} _{o}and U
^{i}equal U
^{i} _{o}, do not need A
^{i} _{o}upgrade, otherwise need to follow these steps to A
^{i} _{o}, U
^{i} _{o}and T
^{i} _{o}upgrade;
G2. calculate U
^{i}with U
_{o}poor, U
^{i}with U
_{o}difference be exactly Cable Structure bearing about the angular displacement of support of initial position, with angular displacement of support vector V, represent angular displacement of support, V equals U
^{i}deduct U
_{o}; Calculate T
^{i}with T
_{o}poor, T
^{i}with T
_{o}difference be exactly that current cable structure steady temperature data are about the variation of initial Cable Structure steady temperature data, T
^{i}with T
_{o}poor with steady temperature change vector S, represent, S equals T
^{i}deduct T
_{o}, S represents the variation of Cable Structure steady temperature data;
G3. first to A
_{o}in Cable Structure bearing apply angular displacement of support constraint, the numerical value of angular displacement of support constraint is just taken from the numerical value of corresponding element in angular displacement of support vector V, then to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation applying is just taken from steady temperature change vector S, to A
_{o}middle Cable Structure bearing applies angular displacement of support constraint and Cable Structure is applied and obtains the current initial mechanical calculating benchmark model A that upgrades after temperature variation
^{i} _{o}, upgrade A
^{i} _{o}time, U
^{i} _{o}all elements numerical value is also used U
^{i}all elements numerical value is corresponding to be replaced, and has upgraded U
^{i} _{o}, T
^{i} _{o}all elements numerical value is also used T
^{i}corresponding replacement of all elements numerical value, upgraded T
^{i} _{o}, so just obtained correctly corresponding to A
^{i} _{o}u
^{i} _{o}and T
^{i} _{o}, d now
^{i} _{o}remain unchanged; When upgrading A
^{i} _{o}after, A
^{i} _{o}the current initial damage of the evaluation object vector d for health status of rope
^{i} _{o}represent A
^{i} _{o}current cable structure steady temperature data vector T for Cable Structure steady temperature
^{i} _{o}represent A
^{i} _{o}current initial Cable Structure bearing angular coordinate vector U for bearing angular coordinate
^{i} _{o}represent; Upgrade C
^{i} _{o}method be: when upgrading A
^{i} _{o}after, by Mechanics Calculation, obtain A
^{i} _{o}in concrete numerical value all monitored amounts, current, these concrete numerical value form C
^{i} _{o};
H. at current initial mechanical calculating benchmark model A
^{i} _{o}basis on, according to step h1, carry out several times Mechanics Calculation to step h4, by calculating, set up unit damage monitored numerical quantity transformation matrices Δ C
^{i}with evaluation object unit change vector D
^{i} _{u};
H1. when the i time circulation starts, directly press step h2 to the listed method acquisition of step h4 Δ C
^{i}and D
^{i} _{u}; At other constantly, when in step g to A
^{i} _{o}after upgrading, must regain Δ C to the listed method of step h4 by step h2
^{i}and D
^{i} _{u}if, in step g not to A
^{i} _{o}upgrade, directly proceed to herein step I and carry out followup work;
H2. at current initial mechanical calculating benchmark model A
^{i} _{o}basis on carry out several times Mechanics Calculation, on calculation times numerical value, equal the quantity N of all evaluation objects, have N evaluation object just to have N calculating; Coding rule according to evaluation object, calculates successively; Calculating each time hypothesis only has an evaluation object on the basis of original damage or centrepoint load, to increase unit damage or centrepoint load unit change again, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable increases unit damage again, if this evaluation object is a centrepoint load, just suppose that this centrepoint load increases centrepoint load unit change again, uses D
^{i} _{uk}record unit damage or the centrepoint load unit change of this increase, wherein k represents to increase the numbering of the evaluation object of unit damage or centrepoint load unit change, D
^{i} _{uk}evaluation object unit change vector D
^{i} _{u}an element, evaluation object unit change vector D
^{i} _{u}coding rule and the vectorial d of element
_{o}the coding rule of element identical; The evaluation object that increases again unit damage or centrepoint load unit change in calculating is each time different from the evaluation object that increases again unit damage or centrepoint load unit change in other calculating, calculate each time the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure, the current calculated value of all monitored amounts that calculate each time forms a monitored amount calculation current vector; When k evaluation object of hypothesis increases unit damage or centrepoint load unit change again, use C
^{i} _{tk}represent corresponding " monitored amount calculation current vector "; While giving in this step each vectorial element numbering, should use same coding rule with other vector in this method, to guarantee any one element in each vector in this step, with in other vector, number identical element, expressed the relevant information of same monitored amount or same target; C
^{i} _{tk}definition mode and vectorial C
_{o}definition mode identical, C
^{i} _{tk}element and C
_{o}element corresponding one by one;
H3. the vectorial C calculating each time
^{i} _{tk}deduct vectorial C
^{i} _{o}obtain a vector, then obtain " numerical value change vector δ a C for monitored amount after each element of this vector is calculated to the unit damage suppose or centrepoint load unit change numerical value divided by this
^{i} _{k}"; There is N evaluation object just to have N " the numerical value change vector of monitored amount ";
H4. by this N " the numerical value change vector of monitored amount ", according to the coding rule of N evaluation object, form successively " the unit damage monitored numerical quantity transformation matrices Δ C that has N row
^{i}"; Unit damage monitored numerical quantity transformation matrices Δ C
^{i}each row corresponding to a monitored amount unit change vector; Unit damage monitored numerical quantity transformation matrices Δ C
^{i}every a line corresponding to same monitored amount the different unit change amplitude when different evaluation objects increase unit damage or centrepoint load unit change; Unit damage monitored numerical quantity transformation matrices Δ C
^{i}coding rule and the vectorial d of row
_{o}the coding rule of element identical, unit damage monitored numerical quantity transformation matrices Δ C
^{i}the coding rule of coding rule and M monitored amount of row identical;
I. define the vectorial d of current name damage
^{i} _{c}with current actual damage vector d
^{i}, d
^{i} _{c}and d
^{i}element number equal the quantity of evaluation object, d
^{i} _{c}and d
^{i}element and evaluation object between be onetoone relationship, d
^{i} _{c}element numerical value represent nominal degree of injury or the nominal centrepoint load variable quantity of corresponding evaluation object, d
^{i} _{c}and d
^{i}with evaluation object initial damage vector d
_{o}element coding rule identical, d
^{i} _{c}element, d
^{i}element and d
_{o}element be onetoone relationship;
J. according to monitored amount current value vector C
^{i}with " the current initial value vector of monitored amount C
^{i} _{o}", " unit damage monitored numerical quantity transformation matrices Δ C
^{i}" and " the vectorial d of current name damage
^{i} _{c}" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, in formula 1 except d
^{i} _{c}other outer amount is known, solves formula 1 and just can calculate the vectorial d of current name damage
^{i} _{c};
K. the current actual damage vector d that utilizes formula 2 to express
^{i}k element d
^{i} _{k}with the current initial damage vector of evaluation object d
^{i} _{o}k element d
^{i} _{ok}with the vectorial d of current name damage
^{i} _{c}k element d
^{i} _{ck}between relation, calculate all elements of current actual damage vector di;
K=1 in formula 2,2,3 ..., N; Vector d
^{i}the coding rule of element and formula (1) in vectorial d
_{o}the coding rule of element identical; d
^{i} _{k}the current actual health status that represents k evaluation object in the i time circulation, if this evaluation object is support cable, so a d in cable system
^{i} _{k}the order of severity that represents its current health problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, d
^{i} _{k}numerical response the degree of lax or damage of this support cable; If this evaluation object is centrepoint load, so a d
^{i} _{k}the actual change amount that represents this centrepoint load; If this evaluation object is centrepoint load, so a d
^{i} _{k}the actual change amount that represents this centrepoint load; By the current actual damage vector of evaluation object d
^{i}in with M
_{1}the M that root support cable is relevant
_{1}individual element takes out, and forms the current actual damage vector of support cable d
^{ci}, the current actual damage vector of support cable d
^{ci}the coding rule and initial rope force vector F of element
_{o}the coding rule of element identical; The current actual damage vector of support cable d
^{ci}h element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., M
_{1}; The current actual damage vector of support cable d
^{ci}middle numerical value be not 0 element corresponding to the support cable of unsoundness problem, from the support cable of these unsoundness problems, identify damaged cable, remaining is exactly slack line; The current actual damage vector of the support cable corresponding with damaged cable d
^{ci}in element numerical expression be the current actual damage of this damaged cable, element numerical value represents while being 100% that this support cable thoroughly loses loadbearing capacity, represents to lose the loadbearing capacity of corresponding proportion in the time of between 0 and 100%; Utilization is at current cable structure steady temperature data vector T
^{i}slack line under condition, that identify in the 1st step and with the current actual damage of support cable vector d
^{ci}these slack lines of expressing, with the current actual equivalent damage degree of its relax level mechanics equivalence, utilize in f step, obtain at current cable structure steady temperature data vector T
^{i}current cable force vector F under condition
^{i}with current support cable two supporting end points horizontal ranges vectors, utilize in c step, obtain at initial Cable Structure steady temperature data vector T
_{o}the initial drift vector of the support cable under condition, the weight vector of the initial free unit length of initial free crosssectional area vector sum, initial rope force vector F
_{o}, utilize current cable structure steady temperature data vector T
^{i}the current steady temperature data of support cable that represent, utilize in c step, obtain at initial Cable Structure steady temperature data vector T
_{o}the support cable initial steady state temperature data representing, the temperature variant physical and mechanical properties parameter of the various materials that the Cable Structure that utilization obtains in c step is used, count the impact of temperature variation on support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanics equivalence calculate slack line, with the relax level of current actual equivalent damage degree equivalence, mechanics equivalent condition is: one, the mechanics parameters without lax initial drift, geometrical property parameter, density and material during with not damaged of the rope of two equivalences is identical; Two, after lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after overall length identical; While meeting abovementioned two mechanics equivalent conditions, the mechanics function of two such support cables in Cable Structure is exactly identical, if replaced after damaged cable with equivalent slack line, Cable Structure any variation can not occur, and vice versa; According to aforementioned mechanics equivalent condition, try to achieve the relax level that those are judged as slack line, relax level is exactly the change amount of support cable drift, has namely determined the long adjustment amount of rope of the support cable that those need adjust Suo Li; Lax identification and the damage identification of support cable have so just been realized; During calculating, institute's demand power is provided by current cable force vector Fi corresponding element; This method is referred to as damaged cable and slack line the support cable of unsoundness problem, referred to as problem rope, so far this method has realized and has rejected problem rope identification impact, Cable Structure that angular displacement of support, centrepoint load variation and structure temperature change, and has realized simultaneously and has rejected angular displacement of support, structure temperature variation and identification support cable health status variable effect, centrepoint load variable quantity;
1. try to achieve the vectorial d of current name damage
^{i} _{c}after, according to formula 3, set up mark vector B
^{i}, formula 4 has provided mark vector B
^{i}the definition of k element;
Element B in formula 4
^{i} _{k}mark vector B
^{i}k element, D
^{i} _{uk}evaluation object unit change vector D
^{i} _{u}k element, d
^{i} _{ck}the vectorial d of the current name damage of evaluation object
^{i} _{c}k element, they all represent the relevant information of k evaluation object, k=1 in formula 4,2,3 ..., N;
If mark vector B m.
^{i}element be 0 entirely, get back to step f and continue this circulation; If mark vector B
^{i}element be not 0 entirely, enter next step, be step n;
N. according to formula 5 calculate next time, i.e. the i+1 time current initial damage vector of the required evaluation object of circulation d
^{i+1} _{o}each element;
D in formula 5
^{i+1} _{ok}the current initial damage vector of the required evaluation object d that next time, circulates for the i+1 time
^{i+1} _{o}k element, d
^{i} _{ok}the current initial damage vector of the evaluation object d that is this, circulates for the i time
^{i} _{o}k element, D
^{i} _{uk}the evaluation object unit change vector D of the i time circulation
^{i} _{u}k element, B
^{i} _{k}the mark vector B of the i time circulation
^{i}k element, k=1 in formula 5,2,3 ..., N;
O. at initial mechanical calculating benchmark model A
_{o}basis on, first to A
_{o}in Cable Structure bearing apply angular displacement of support constraint, the numerical value of angular displacement of support constraint is just taken from the numerical value of corresponding element in angular displacement of support vector V, then to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation applying is just taken from steady temperature change vector S, then to make the health status of rope be d
^{i+1} _{o}after obtain be exactly next time, i.e. the i+1 time required Mechanics Calculation benchmark model A of circulation
^{i+1}; Obtain A
^{i+1}after, by Mechanics Calculation, obtain A
^{i+1}in concrete numerical value all monitored amounts, current, these concrete numerical value form next time, the vectorial C of the current initial value of required monitored amount that circulates for the i+1 time
^{i+1} _{o};
P. take off once, i.e. the i+1 time required current initial Cable Structure steady temperature data vector T of circulation
^{i+1}o equals the current initial Cable Structure steady temperature data vector T of the i time circulation
^{i} _{o}; The required current initial Cable Structure bearing angular coordinate vector U of the i+1 time circulation next time, i.e.
^{i+1} _{o}equal the current initial Cable Structure bearing angular coordinate vector U of the i time circulation
^{i} _{o};
Q. get back to step f, start circulation next time.
Beneficial effect: in current published correlation technique, some only can when other all conditions is constant, (load of only only having structure to bear changes, and structural health conditions etc. are all constant) variation of recognition structure bearing load, some only can (only only have structural health conditions to change when other all conditions is constant, and the load that structure is born etc. are constant) variation of recognition structure health status, some only can (only only have structure temperature and structural health conditions to change when other all conditions is constant, and the load that structure is born is constant) variation of recognition structure health status, when the load of bearing in structure, structure (environment) temperature and structural health conditions change simultaneously, when Cable Structure generation angular displacement of support, also there is no at present a kind of disclosed, effective method load that recognition structure bears simultaneously and the variation of structural health conditions, when structure bearing load and structure (environment) temperature change simultaneously, the variation and the structure temperature that also do not have effective method can reject Cable Structure angular displacement of support, structure bearing load change the impact on structural health conditions recognition result in other words, changing an angle sees, in current disclosed method, thereby also there is not rejecting the correct knowledge method for distinguishing of realizing centrepoint load intensity of variation of Cable Structure angular displacement of support, structure temperature variation and the impact of support cable health status, and concerning structure, the identification of load change is also very important, with existing method, compare, this method can be when Cable Structure generation angular displacement of support, when the centrepoint load of bearing in structure and structure temperature change, can reject Cable Structure angular displacement of support, variation and structure temperature that structure is born centrepoint load change the impact on structural health conditions recognition result, can identify very exactly problem rope, solved monitoring structural health conditions field problem in the urgent need to address, otherwise, if can not reject Cable Structure angular displacement of support, structure temperature changes and the impact of the variation of the centrepoint load that structure is born, identification problem rope exactly just, moreover, this method is when identifying damaged cable, can also identify the variation of centrepoint load simultaneously, be that this method can be rejected Cable Structure angular displacement of support, structure temperature changes and the impact of support cable health status variation, realize the correct identification of centrepoint load intensity of variation, otherwise, if can not reject the impact that Cable Structure angular displacement of support, structure temperature variation and support cable health status change, just can not identify exactly the intensity of variation of centrepoint load.That is to say, this method has realized two kinds of functions that existing method can not possess, respectively: one, when Cable Structure generation angular displacement of support, during the centrepoint load of bearing in structure and structure (environment) temperature variation, can reject Cable Structure angular displacement of support, centrepoint load variation and structure temperature and change the impact on Cable Structure health status recognition result, thereby identify exactly the structure health monitoring method of problem rope; Two, this method is when identifying problem rope, can also identify the variation of centrepoint load simultaneously, be that this method can be rejected Cable Structure angular displacement of support, structure temperature changes and the impact of support cable health status variation, realize the correct identification of centrepoint load intensity of variation.
Embodiment
For cable structure health monitoring problem, this method has realized two kinds of functions that existing method can not possess, respectively: one, when Cable Structure generation angular displacement of support, during the centrepoint load of bearing in structure and structure (environment) temperature variation, can reject Cable Structure angular displacement of support, centrepoint load variation and structure temperature and change the impact on Cable Structure health status recognition result, thereby identify exactly the structure health monitoring method of problem rope (support cable of unsoundness problem); Two, this method is when identifying problem rope, can also identify the variation of centrepoint load simultaneously, be that this method can be rejected Cable Structure angular displacement of support, structure temperature changes and the impact of support cable health status variation, realize the correct identification of centrepoint load intensity of variation.The following describes of the embodiment of this method is in fact only exemplary, and object is never to limit application or the use of this method.
This method adopts a kind of algorithm, and this algorithm is for the variation of identification problem rope and centrepoint load.During concrete enforcement, the following step is a kind of in the various steps that can take.
The first step: the quantity of first confirming the centrepoint load that may change that Cable Structure is born.The feature of the centrepoint load of bearing according to Cable Structure, confirm wherein " centrepoint load likely changing ", or all centrepoint load is considered as " centrepoint load likely changing ", establishes total JZW the centrepoint load that may change.
Centrepoint load is divided into two kinds of concentrated force and concentrated couples, in coordinate system, for example, in Descartes's rectangular coordinate system, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, in the method a concentrated force component or a concentrated couple component is called to a centrepoint load.
If the quantity sum of the quantity of the support cable of Cable Structure and JZW " centrepoint load likely changing " is N.For sake of convenience, it is " evaluation object " that this method unitedly calls evaluated support cable and " centrepoint load likely changing ", total N evaluation object.Give evaluation object serial number, this numbering will be for generating vector sum matrix in subsequent step.
Monitored multiclass parameter can comprise: Suo Li, strain, angle and volume coordinate, be described below respectively:
If total Q root support cable in cable system, the monitored rope force data of Cable Structure is by M in Cable Structure
_{1}the M of individual appointment rope
_{1}individual rope force data is described, and the variation of Cable Structure Suo Li is exactly the variation of the Suo Li of all appointment ropes.Each total M
_{1}individual cable force measurement value or calculated value characterize the rope force information of Cable Structure.M
_{1}be one and be not less than 0 integer.
The monitored strain data of Cable Structure can be by K in Cable Structure
_{2}l individual specified point and each specified point
_{2}the strain of individual assigned direction is described, and the variation of Cable Structure strain data is exactly K
_{2}the variation of all tested strains of individual specified point.Each total M
_{2}(M
_{2}=K
_{2}* L
_{2}) individual strain measurement value or calculated value characterize Cable Structure strain.M
_{2}be one and be not less than 0 integer.
The monitored angledata of Cable Structure is by K in Cable Structure
_{3}l individual specified point, that cross each specified point
_{3}h individual appointment straight line, each appointment straight line
_{3}individual angle coordinate component is described, and the variation of Cable Structure angle is exactly the variation of angle coordinate components appointment straight lines all specified points, all, all appointments.Each total M
_{3}(M
_{3}=K
_{3}* L
_{3}* H
_{3}) individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure.M
_{3}be one and be not less than 0 integer.
The monitored shape data of Cable Structure is by K in Cable Structure
_{4}l individual specified point and each specified point
_{4}the volume coordinate of individual assigned direction is described, and the variation of Cable Structure shape data is exactly K
_{4}the variation of all coordinate components of individual specified point.Each total M
_{4}(M
_{4}=K
_{4}* L
_{4}) individual measurement of coordinates value or calculated value characterize Cable Structure shape.M
_{4}be one and be not less than 0 integer.
Comprehensive abovementioned monitored amount, whole Cable Structure has M(M=M
_{1}+ M
_{2}+ M
_{3}+ M
_{4}) individual monitored amount, definition parameter K(K=M
_{1}+ K
_{2}+ K
_{3}+ K
_{4}), K and M must not be less than N.
For simplicity, in the method by " monitored all parameters of Cable Structure " referred to as " monitored amount ".Give M monitored amount serial number, this numbering will be for generating vector sum matrix in subsequent step.This method is with representing this numbering with variable j, j=1, and 2,3 ..., M.
Determine " the temperature survey calculating method of the Cable Structure of this method ", the method concrete steps are as follows:
A step: inquiry or actual measurement (can be measured by conventional thermometry, for example use thermal resistance to measure) obtain the temperature variant thermal conduction study parameter of Cable Structure composition material and Cable Structure environment of living in, utilize the geometry measured data of design drawing, asconstructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure (for example finite element model).Inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, statistics obtains interior during this period of time cloudy quantity and is designated as T cloudy day, statistics obtains 0 highest temperature and the lowest temperature between latter 30 minutes of the moment of sunrise next day at each cloudy day in T cloudy day, sunrise on the meteorology that sunrise refers to base area revolutions constantly and the rule that revolves round the sun is definite constantly, the sunrise that can inquire about data or calculate each required day by conventional meteorology constantly, each cloudy day 0 up to next day sunrise constantly the highest temperature between latter 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, there is T cloudy day, the maximum temperature difference that just has the daily temperature at T cloudy day, get maximal value in the maximum temperature difference of daily temperature at T cloudy day for reference to temperature difference per day, with reference to temperature difference per day, be designated as Δ T
_{r}, between inquiry Cable Structure location and Altitude Region, place, be no less than temperature that the meteorological data in recent years of 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and Altitude Region, place about the maximum rate of change Δ T of sea level elevation
_{h}, for convenience of narration, get Δ T
_{h}unit be ℃/m, on the surface of Cable Structure, get " R Cable Structure surface point ", get the Specific Principles of " R Cable Structure surface point " narrates in step b3, after will by actual observation record, obtain the temperature of this R Cable Structure surface point, claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, by Calculation of Heat Transfer, obtain the temperature of this R Cable Structure surface point, just claim that the temperature data calculating is " R Cable Structure surface temperature computational data ".From the residing minimum height above sea level of Cable Structure to the highest height above sea level, in Cable Structure, uniform choosing is no less than three different sea level elevations, if for example the sea level elevation of Cable Structure is between 0m to 200m, can choose height above sea level 0m so, 50m, 100m and height above sea level 200m, intersect with imaginary surface level and Cable Structure surface at the sea level elevation place choosing at each, obtain intersection, the crossing cross surface that obtains of surface level and Cable Structure, intersection is the outer edge line of cross surface, at the intersection place on surface level and Cable Structure surface, choose 6 points, outer normal from selected point straw line body structure surface, all outer normal directions of choosing are called " measuring Cable Structure along the direction of the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness.In the measurement Cable Structure of choosing along in 6 directions of the Temperature Distribution of wall thickness, first according to the meteorological data throughout the year in region, Cable Structure position and the physical dimension of Cable Structure, volume coordinate, the sunny slope of definite Cable Structure such as Cable Structure surrounding environment and in the shade, the sunny slope of Cable Structure and in the shade face are surperficial parts for Cable Structure, the sea level elevation place choosing at each, aforementioned intersection respectively has one section in sunny slope and in the shade, two sections of these of intersection respectively have a mid point, cross these two mid points and get the outer normal of Cable Structure, this method is called the sunny slope outer normal of Cable Structure and in the shade outer normal of Cable Structure by these two outer normals, this method is called the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure by these two outer normal directions, the outer normal of obvious sunny slope and the outer normal of in the shade all intersect with aforementioned intersection, also just there are two intersection points, these two intersection points are divided into two line segments by intersection, on two line segments, get respectively 2 points, totally 4 points, taken point is divided into equal in length 3 sections by each line segment in two line segments of intersection, at these 4 some places, get the outer normal on Cable Structure surface, at each selected sea level elevation place, just chosen altogether like this outer normal on 6 Cable Structure surfaces, the direction of 6 outer normals is exactly " measuring Cable Structure along the direction of the Temperature Distribution of wall thickness ".There are two intersection points on the surface of each " measures Cable Structure along the direction of the Temperature Distribution of wall thickness " line and Cable Structure, if Cable Structure is hollow, one, these two intersection points are on Cable Structure outside surface, another is on inside surface, if Cable Structure is solid, these two intersection points are all on Cable Structure outside surface, connect these two intersection points and obtain a straightline segment, on straightline segment, choose again three points, these three these straightline segments of naming a person for a particular job are divided into four sections, three points measuring that Cable Structure chooses at this and two end points of straightline segment, the temperature that amounts to 5 points, concrete can first hole in Cable Structure, how temperature sensor is embedded in to this 5 some places, the temperature recording is called this place " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure Cable Structure along the direction of the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ".If chosen H different sea level elevation, at each sea level elevation place, choose B and measured Cable Structure along the direction of the Temperature Distribution of wall thickness, along each, measure Cable Structure and in Cable Structure, chosen E point along the direction of the Temperature Distribution of wall thickness, wherein H and E are not less than 3, B is not less than 2, if HBE is the product of H and B and E, corresponding total HBE " measuring Cable Structure along the point of the temperature profile data of thickness ", after will obtain the temperature that this HBE " measures Cable Structure along the point of the temperature profile data of thickness " by actual measurement, claim that the temperature data that actual measurement obtains is " HBE Cable Structure is along thickness temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, by Calculation of Heat Transfer, obtain this HBE and measure Cable Structure along the temperature of the point of the temperature profile data of thickness, just claim that the temperature data calculating is " HBE Cable Structure is along thickness temperature computation data ", if BE is the product of B and E, total BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in sea level elevation place of choosing at each in this method.In Cable Structure location, according to meteorology, measure temperature and require to choose a position, will in this position actual observation record, obtain meeting the temperature that meteorology is measured the Cable Structure place environment of temperature requirement, in the onsite spaciousness of Cable Structure, unobstructed place chooses a position, this position should can obtain in each day of the whole year this ground the most sufficient sunshine of getable this day (as long as had sunrise the same day, this position just should be arrived by solar radiation), for example, flat board (square that for example the wide 3mm of 30cm is thick is dull and stereotyped) at carbon steel material of this position of sound production (No. 45 carbon steels), be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, reference plate can be placed in and meet the top that meteorology temperature is measured the wooden thermometer screen requiring, the one side of this reference plate on the sunny side, (be for example called sunny slope, in the time of on the Northern Hemisphere, sunny slope faces up towards south, full daytime is all by sunshine, sunny slope should have the suitable gradient to make snow can not accumulate or clear up sunny slope after snow), the sunny slope of reference plate is coarse and (being conducive to accept solar irradiation) dark color, the sunny slope of reference plate should can obtain in each day of the whole year one flat plate on this ground the most sufficient sunshine of getable this day, the nonsunny slope of reference plate is covered with insulation material (for example thick calcium carbonate insulation material of 5mm), RealTime Monitoring record is obtained to the temperature of the sunny slope of reference plate.
B step, RealTime Monitoring (can be measured by conventional thermometry, for example use thermal resistance to measure, for example, every temperature data of 10 minutes survey records) record obtains R Cable Structure surface temperature measured data of abovementioned R Cable Structure surface point, RealTime Monitoring (can be measured by conventional thermometry simultaneously, for example use thermal resistance to measure, for example, every temperature data of 10 minutes survey records) obtain previously defined Cable Structure along the temperature profile data of thickness, RealTime Monitoring (can be measured by conventional thermometry simultaneously, for example in the wooden thermometer screen that meets meteorology temperature measurement requirement, lay thermal resistance and measure temperature, for example, every temperature data of 10 minutes survey records) record obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement, by RealTime Monitoring, (can measure by conventional thermometry, for example in the wooden thermometer screen that meets meteorology temperature measurement requirement, lay thermal resistance and measure temperature, for example, every temperature data of 10 minutes survey records) record obtains being carved at sunrise the sunrise next day temperature measured data sequence of the Cable Structure place environment between latter 30 minutes constantly the same day, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data that was carved at sunrise the Cable Structure place environment between latter 30 minutes of the moment of sunrise next day the same day, find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment, by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, deduct and be carved at sunrise the sunrise next day maximum temperature difference between latter 30 minutes constantly on same day that minimum temperature obtains Cable Structure place environment, be designated as Δ T
_{emax}, temperature measured data sequence by Cable Structure place environment (for example first carries out curve fitting to the temperature measured data sequence of Cable Structure place environment by conventional mathematical computations, then by asking curve to the derivative of time or by ask the rate of change of each point corresponding to survey record data time to the time on curve by numerical method) obtain the temperature of Cable Structure place environment about the rate of change of time, this rate of change is also along with the time changes, by RealTime Monitoring, (can measure by conventional thermometry, for example use the temperature of the dull and stereotyped sunny slope of thermal resistance witness mark, for example, every temperature data of 10 minutes survey records) obtain being carved at sunrise the same day sunrise next day measured data sequence of the temperature of the sunny slope of the reference plate between latter 30 minutes constantly, the measured data sequence of the temperature of the sunny slope of reference plate by be carved at sunrise the same day next day sunrise constantly the measured data of the temperature of the sunny slope of the reference plate between latter 30 minutes according to time order and function order, arrange, find maximum temperature and minimum temperature in the measured data sequence of temperature of sunny slope of reference plate, by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, deduct and be carved at sunrise the sunrise next day maximum temperature difference between latter 30 minutes constantly on same day of temperature that minimum temperature obtains the sunny slope of reference plate, be designated as Δ T
_{pmax}, by RealTime Monitoring, (can measure by conventional thermometry, for example use thermal resistance to measure Cable Structure surface point, for example, every temperature data of 10 minutes survey records) record obtains being carved at sunrise the sunrise next day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between latter 30 minutes constantly the same day, there is R Cable Structure surface point just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence by be carved at sunrise on same day of a Cable Structure surface point sunrise next day constantly the Cable Structure surface temperature measured data between latter 30 minutes according to time order and function order, arrange, find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence, by the maximum temperature in each Cable Structure surface temperature measured data sequence, deduct and be carved at sunrise the sunrise next day maximum temperature difference between latter 30 minutes constantly on same day that minimum temperature obtains the temperature of each Cable Structure surface point, there is R Cable Structure surface point just to have and be carved at sunrise the sunrise next day maximum temperature difference numerical value between latter 30 minutes constantly R the same day, maximal value is wherein designated as Δ T
_{smax}, by each Cable Structure surface temperature measured data sequence, by conventional mathematical computations, (for example first each Cable Structure surface temperature measured data sequence is carried out curve fitting, then by asking curve to the derivative of time or by ask the rate of change of each point corresponding to survey record data time to the time on curve by numerical method) obtain the temperature of each Cable Structure surface point about the rate of change of time, the temperature of each Cable Structure surface point about the rate of change of time also along with the time changes.By RealTime Monitoring, obtain being carved at sunrise the same day between latter 30 minutes of the moment of sunrise next day, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", the sea level elevation place that calculating is chosen at each amounts to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and the difference of minimum temperature, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", chosen H different sea level elevation and just had H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T
_{tmax}.
C step, measures and calculates acquisition Cable Structure steady temperature data; First, determine the moment that obtains Cable Structure steady temperature data, the condition relevant to the moment that determines acquisition Cable Structure steady temperature data has six, first condition is the moment that obtains Cable Structure steady temperature data to be carved at sunset sunrise next day constantly between latter 30 minutes between the same day, sunset refers to sunset on base area revolutions and the definite meteorology of revolution rule constantly constantly, and the sunset that can inquire about data or calculate each required day by conventional meteorology constantly; The a condition of second condition be the same day be carved at sunrise next day sunrise constantly between latter 30 minutes during this period of time in, Δ T
_{pmax}with Δ T
_{smax}all be not more than 5 degrees Celsius; Second must be satisfied b condition be the same day be carved at sunrise next day sunrise constantly between latter 30 minutes during this period of time in, above, measure the Δ T calculating
_{emax}be not more than with reference to temperature difference per day Δ T
_{r}, and above, measure the Δ T calculating
_{pmax}deduct 2 degrees Celsius and be not more than Δ T
_{emax}, and above, measure the Δ T calculating
_{smax}be not more than Δ T
_{pmax}; Only needing to meet in a condition of second and b condition one is just called and meets second condition; The 3rd condition is that the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data; The 4th condition is that the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data; The 5th condition is in the moment that obtains Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is to be carved at sunrise the sunrise next day minimal value between latter 30 minutes constantly the same day; The 6th condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T
_{tmax}be not more than 1 degree Celsius.This method is utilized abovementioned six conditions, any one in following three kinds of moment is called to " obtaining the mathematics of Cable Structure steady temperature data constantly ", the first is first moment to the 5th condition meeting in abovementioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data " constantly, the second is the moment that only meets the 6th condition in abovementioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data " constantly, the third is first moment to the 6th condition simultaneously meeting in abovementioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data " constantly, when obtaining the mathematics of Cable Structure steady temperature data, be exactly in this method during in constantly one of physical record data constantly, the moment that obtains Cable Structure steady temperature data be exactly obtain Cable Structure steady temperature data mathematics constantly, if obtain the mathematics of Cable Structure steady temperature data and be not constantly any in constantly of physical record data in this method constantly, get this method close to moment of the mathematics that obtains Cable Structure steady temperature data those physical record data constantly for obtaining the moment of Cable Structure steady temperature data, this method is carried out the relevant health monitoring analysis of Cable Structure by using in the amount that obtains the moment survey record of Cable Structure steady temperature data, this method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, i.e. this Cable Structure temperature constantly temporal evolution not, and this is exactly the moment of the acquisition Cable Structure steady temperature data of this method constantly, then, according to Cable Structure heat transfer characteristic, utilize to obtain R Cable Structure surface temperature measured data and " HBE Cable Structure is along thickness temperature measured data " in the moment of Cable Structure steady temperature data, utilize the thermal conduction study computation model (for example finite element model) of Cable Structure, for example, by conventional Calculation of Heat Transfer (finite element method), obtain in the Temperature Distribution of Cable Structure that obtains the moment of Cable Structure steady temperature data, now calculate by stable state in the temperature field of Cable Structure, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data, also comprise that Cable Structure is in the accounting temperature of selected HBE " measuring Cable Structure along the point of the temperature profile data of thickness " above, the accounting temperature of HBE " measuring Cable Structure along the point of the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", when R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure stable state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ".While getting " R Cable Structure surface point " on the surface of Cable Structure, the quantity of " R Cable Structure surface point " and necessary three conditions that meet that distribute, first condition is when Cable Structure temperature field is during in stable state, on Cable Structure surface the temperature of any point be by " R Cable Structure surface point " with Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point while obtaining, on the Cable Structure surface that linear interpolation obtains, on the temperature of this arbitrfary point and Cable Structure surface, the error of the actual temperature of this arbitrfary point is not more than 5%; Cable Structure surface comprises support cable surface; Second condition is that in " R Cable Structure surface point ", the quantity at the point of same sea level elevation is not less than 4, and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point "; " R Cable Structure surface point " is not more than 0.2 ℃ divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation
_{h}the numerical value obtaining, gets Δ T for convenience of narration
_{h}unit be ℃/m that the unit of getting Δ h for convenience of narration is m; " R Cable Structure surface point " refers to while only considering sea level elevation along the definition of adjacent Cable Structure surface point between two of sea level elevation, in " R Cable Structure surface point ", do not have a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two; The 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to the geometric properties of Cable Structure and bearing data, in Cable Structure, find and be subject to the sunshineduration position of those surface points the most fully the whole year, in " R Cable Structure surface point ", having a Cable Structure surface point at least is the annual point being subject in the most sufficient those surface points of sunshineduration in Cable Structure.
Second step: set up initial mechanical calculating benchmark model A
_{o}.
In Cable Structure completion, or before setting up health monitoring systems, according to " the temperature survey calculating method of the Cable Structure of this method " measurement, calculating " Cable Structure steady temperature data " (can measure by conventional thermometry, for example use thermal resistance to measure), " Cable Structure steady temperature data " now use vector T
_{o}represent, be called initial Cable Structure steady temperature data vector T
_{o}.In actual measurement, obtain T
_{o}time, namely at the synchronization that obtains the moment of initial Cable Structure steady temperature data vector, use conventional method directly to measure the initial value of all monitored amounts that calculate Cable Structure, form monitored amount initial value vector C
_{o}.
Can be specifically in this method according to following method at the synchronization that obtains the moment of soandso (such as initial or current etc.) Cable Structure steady temperature data vector, use soandso method measurement to calculate the data of the monitored amount of soandso measured amount (for example all monitored amount of Cable Structure): at the survey record temperature (temperature that comprises Cable Structure place environment, the temperature of the sunny slope of reference plate and Cable Structure surface temperature) time, for example, every temperature of 10 minutes survey records, so simultaneously equally also every 10 minutes the monitored amount of soandso measured amount of survey record (for example all monitored amount of Cable Structure) data.Once determine the moment that obtains Cable Structure steady temperature data, for example, be just called at the synchronization that obtains the moment of Cable Structure steady temperature data with the data of the monitored amount of soandso measured amount (all monitored amount of Cable Structure) that obtain the moment synchronization of Cable Structure steady temperature data so, use soandso method to measure the data of the monitored amount of soandso measured amount that computing method obtain.
Use conventional method (consult reference materials or survey) to obtain temperature variant physical parameter (for example thermal expansivity) and the mechanical property parameters (for example elastic modulus, Poisson ratio) of the various materials that Cable Structure used.
In actual measurement, obtain initial Cable Structure steady temperature data vector T
_{o}synchronization, directly measure the initial Suo Li that calculates all support cables, form initial rope force vector F
_{o}; According to Cable Structure design data, completion data obtain that all support cables are in free state that Suo Li is the length of 0 o'clock, the weight of crosssectional area during in free state and the unit length during in free state, and the temperature of all support cables while obtaining these three kinds of data, utilize on this basis temperature variant physical function parameter and the mechanical property parameters of all support cables, according to conventional physical computing, obtain all support cables at initial Cable Structure steady temperature data vector T
_{o}suo Li under condition is that the length of 0 o'clock all support cable, crosssectional area and the Suo Li that Suo Li is 0 o'clock all support cable are the weight of the unit length of 0 o'clock all support cable, forms successively the initial drift vector l of support cable
_{o}, initial free crosssectional area vector A
_{o}weight vector ω with initial free unit length
_{o}, the initial drift vector l of support cable
_{o}, initial free crosssectional area vector A
_{o}weight vector ω with initial free unit length
_{o}the coding rule and initial rope force vector F of element
_{o}the coding rule of element identical.
At Actual measurement, obtain initial Cable Structure steady temperature data vector T
_{o}time, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Actual measurement data of Cable Structure comprise that Nondestructive Testing Data of support cable etc. can express the data of the health status of rope, the initial geometric data of Cable Structure, rope force data, drawbar pull data, initial Cable Structure bearing generalized coordinate data (comprise that bearing is about Descartes's rectangular coordinate system X, Y, the volume coordinate of Z axis and angular coordinate are initial Cable Structure bearing spatial data and initial Cable Structure bearing angular data), Cable Structure centrepoint load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure modal data, structural strain data, structure angle measurement data, the measured datas such as structure space measurement of coordinates data.Initial Cable Structure bearing angular data forms initial Cable Structure bearing angular coordinate vector U
_{o}.The initial geometric data of Cable Structure can be the spatial data that the spatial data of the end points of all ropes adds a series of point in structure, and object is to determine according to these coordinate datas the geometric properties of Cable Structure.For cablestayed bridge, the spatial data that initial geometric data can be the end points of all ropes adds the spatial data of some points on bridge two ends, socalled bridge type data that Here it is.Data and the Cable Structure centrepoint load measurement data of utilizing the Nondestructive Testing Data etc. of support cable can express the health status of support cable are set up evaluation object initial damage vector d
_{o}, use d
_{o}represent that Cable Structure is (with initial mechanical calculating benchmark model A
_{o}the initial health of evaluation object expression).If while there is no the Nondestructive Testing Data of support cable and the data of other health status that can express support cable, or can think that structure original state is not damaged during without relaxed state, vectorial d
_{o}in each element numerical value relevant to support cable get 0, if d
_{o}evaluation object corresponding to some elements be some centrepoint load, in this method, get d
_{o}this element numerical value be 0, the initial value that represents the variation of this centrepoint load is 0.Utilize the design drawing, asconstructed drawing of Cable Structure and the initial measured data of Cable Structure, temperature variant physical and mechanical properties parameter, the initial Cable Structure bearing angular coordinate vector U of the various materials that the Nondestructive Testing Data of support cable, Cable Structure are used
_{o}with initial Cable Structure steady temperature data vector T
_{o}, utilize mechanics method (for example finite element method) to count " Cable Structure steady temperature data " and set up initial mechanical calculating benchmark model A
_{o}.
No matter which kind of method to obtain initial mechanical calculating benchmark model A by
_{o}, counting " Cable Structure steady temperature data " (is initial Cable Structure steady temperature data vector T
_{o}), based on A
_{o}the Cable Structure computational data calculating must approach its measured data very much, and error generally must not be greater than 5%.Like this can utility A
_{o}calculate Suo Li computational data, strain computational data, Cable Structure shape computational data and displacement computational data under the analog case of gained, Cable Structure angledata, Cable Structure spatial data etc., the measured data when approaching reliably institute's analog case and truly occurring.Model A
_{o}evaluation object initial damage vector d for the health status of middle support cable
_{o}represent initial Cable Structure steady temperature data vector T for Cable Structure steady temperature data
_{o}represent.Due to based on A
_{o}the evaluation that calculates all monitored amounts approaches the initial value (actual measurement obtains) of all monitored amounts very much, so also can be used in A
_{o}basis on, carry out Mechanics Calculation obtains, A
_{o}the evaluation of each monitored amount form monitored amount initial value vector C
_{o}.Corresponding to A
_{o}" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T
_{o}"; Corresponding to A
_{o}evaluation object initial damage vector d for evaluation object health status
_{o}represent; Corresponding to A
_{o}monitored amount initial value vector C for the initial value of all monitored amounts
_{o}represent.Corresponding to A
_{o}initial Cable Structure bearing angular coordinate vector U for Cable Structure bearing angular data
_{o}represent; T
_{o}, U
_{o}and d
_{o}a
_{o}parameter, C
_{o}by A
_{o}mechanics Calculation result form.
The 3rd step: in the method, alphabetical i, except representing that significantly, the place of number of steps, alphabetical i only represents cycle index, circulates for the i time; When the i time circulation starts, the current initial mechanical calculating benchmark model of Cable Structure that need to set up or that set up is designated as current initial mechanical calculating benchmark model A
^{i} _{o}, A
_{o}and A
^{i} _{o}count temperature parameter, can accounting temperature change the Effect on Mechanical Properties to Cable Structure; When the i time circulation starts, corresponding to A
^{i} _{o}" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T
^{i} _{o}represent vector T
^{i} _{o}definition mode and vector T
_{o}definition mode identical, T
^{i} _{o}element and T
_{o}element corresponding one by one; That the i time circulation needs while starting, corresponding to the current initial mechanical calculating benchmark model A of Cable Structure
^{i} _{o}cable Structure bearing angular data form current initial Cable Structure bearing angular coordinate vector U
^{i} _{o}, set up for the first time the current initial mechanical calculating benchmark model A of Cable Structure
^{i} _{o}time, U
^{i} _{o}just equal U
_{o}.The current initial damage vector of evaluation object that the i time circulation needs while starting is designated as d
^{i} _{o}, d
^{i} _{o}cable Structure A while representing this circulation beginning
^{i} _{o}the health status of evaluation object, d
^{i} _{o}definition mode and d
_{o}definition mode identical, d
^{i} _{o}element and d
_{o}element corresponding one by one; When the i time circulation starts, the initial value of all monitored amounts, with the current initial value vector of monitored amount C
^{i} _{o}represent vectorial C
^{i} _{o}definition mode and vectorial C
_{o}definition mode identical, C
^{i} _{o}element and C
_{o}element corresponding one by one, the current initial value vector of monitored amount C
^{i} _{o}expression is corresponding to A
^{i} _{o}the concrete numerical value of all monitored amounts; T
^{i} _{o}and d
^{i} _{o}a
^{i} _{o}characterisitic parameter; C
^{i} _{o}by A
^{i} _{o}mechanics Calculation result form; When circulation starts for the first time, A
^{i} _{o}be designated as A
^{1} _{o}, set up A
^{1} _{o}method for making A
^{1} _{o}equal A
_{o}; When circulation starts for the first time, T
^{i} _{o}be designated as T
^{1} _{o}, set up T
^{1} _{o}method for making T
^{1} _{o}equal T
_{o}; When circulation starts for the first time, U
^{i} _{o}be designated as U
^{1} _{o}, set up U
^{1} _{o}method for making U
^{1} _{o}equal U
_{o}; When circulation starts for the first time, d
^{i} _{o}be designated as d
^{1} _{o}, set up d
^{1} _{o}method for making d
^{1} _{o}equal d
_{o}; When circulation starts for the first time, C
^{i} _{o}be designated as C
^{1} _{o}, set up C
^{1} _{o}method for making C
^{1} _{o}equal C
_{o}.
The 4th step: the hardware components of pass line structural healthy monitoring system.Hardware components at least comprises: monitored amount monitoring system is (for example, containing measurement of angle subsystem, cable force measurement subsystem, strain measurement subsystem, volume coordinate is measured subsystem, signal conditioner etc.), Cable Structure bearing angular coordinate monitoring system is (containing angle measuring sensor, signal conditioner etc.), Cable Structure temperature monitoring system is (containing temperature sensor, signal conditioner etc.) and Cable Structure ambient temperature measurement system (containing temperature sensor, signal conditioner etc.), support cable cable force monitoring system, the space coordinate monitoring system of the supporting end points of support cable, signal (data) collector, computing machine and the panalarm of communicating by letter.Each bearing angular coordinate of each monitored amount, Cable Structure, each temperature, the Suo Li of each root support cable, the volume coordinate of the supporting end points of each root support cable must arrive by monitored system monitoring, and monitoring system is transferred to signal (data) collector by the signal monitoring; Signal is delivered to computing machine through signal picker; Computing machine is responsible for the health monitoring software of the evaluation object of operation Cable Structure, comprises the signal that the transmission of tracer signal collector comes; When monitoring evaluation object health status and change, computer control communication panalarm to monitor staff, owner and (or) personnel of appointment report to the police.
The 5th step: establishment the system software of installation and operation this method on computers, this software will complete the functions (being all work that can complete with computing machine in this specific implementation method) such as monitoring that this method required by task wants, record, control, storage, calculating, notice, warning.
The 6th step: step starts circulation running thus, in structure military service process, according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement, obtain the current data of Cable Structure steady temperature data, the current data of all " Cable Structure steady temperature data " forms current cable structure steady temperature data vector T
^{i}, vector T
^{i}definition mode and vector T
_{o}definition mode identical, T
^{i}element and T
_{o}element corresponding one by one; In actual measurement vector T
^{i}time, namely obtaining current cable structure steady temperature data vector T
^{i}the synchronization in the moment, actual measurement obtains the currency of all monitored amounts in Cable Structure, all these numerical value form monitored amount current value vector C
^{i}, vectorial C
^{i}definition mode and vectorial C
_{o}definition mode identical, C
^{i}element and C
_{o}element corresponding one by one, represent that identical monitored amount is at numerical value in the same time not.
In actual measurement, obtain current cable structure steady temperature data vector T
^{i}time, actual measurement obtains Cable Structure bearing angular coordinate current data, and all data form current cable structure actual measurement bearing angular coordinate vector U
^{i}.
In actual measurement, obtain current cable structure steady temperature data vector T
^{i}synchronization, actual measurement obtains all M in Cable Structure
_{1}the rope force data of root support cable, all these rope force datas form current cable force vector F
^{i}, the element of vectorial Fi and vectorial F
_{o}the coding rule of element identical; In actual measurement, obtain current cable structure steady temperature data vector T
^{i}synchronization, Actual measurement obtains all M
_{1}the volume coordinate of two supporting end points of root support cable, the volume coordinate of two the supporting end points in the horizontal direction difference of component is exactly two supporting end points horizontal ranges, all M
_{1}two supporting end points horizontal range data of root support cable form current support cable two supporting end points horizontal range vector l
^{i} _{x}, current support cable two supporting end points horizontal range vector l
^{i} _{x}the coding rule and initial rope force vector F of element
_{o}the coding rule of element identical.
The 7th step: obtaining current cable structure actual measurement bearing angular coordinate vector U
^{i}with current cable structure steady temperature data vector T
^{i}after, compare respectively U
^{i}and U
^{i} _{o}, T
^{i}and T
^{i} _{o}if, U
^{i}equal U
^{i} _{o}and T
^{i}equal T
^{i} _{o}, do not need A
^{i} _{o}, U
^{i} _{o}and T
^{i} _{o}upgrade, otherwise need to be to current initial mechanical calculating benchmark model A
^{i} _{o}, current initial Cable Structure bearing angular coordinate vector U
^{i} _{o}, current initial Cable Structure steady temperature data vector T
^{i} _{o}with the current initial value vector of monitored amount C
^{i} _{o}upgrade, and the current initial damage vector of evaluation object d
^{i} _{o}remain unchanged, update method follows these steps to a and carries out to step c:
A. calculate U
^{i}with U
_{o}poor, U
^{i}with U
_{o}difference be exactly Cable Structure bearing about the angular displacement of support of initial position, with angular displacement of support vector V, represent angular displacement of support, V equals U
^{i}deduct U
_{o}, between the element in angular displacement of support vector V and angular displacement of support component, be onetoone relationship, in angular displacement of support vector V, the numerical value of an element is corresponding to the angular displacement of an assigned direction of an appointment bearing.
B. calculate T
^{i}with T
_{o}poor, T
^{i}with T
_{o}difference be exactly that current cable structure steady temperature data are about the variation of initial Cable Structure steady temperature data, T
^{i}with T
_{o}poor with steady temperature change vector S, represent, S equals T
^{i}deduct T
_{o}, S represents the variation of Cable Structure steady temperature data.
C. first to A
_{o}in Cable Structure bearing apply angular displacement of support constraint, the numerical value of angular displacement of support constraint is just taken from the numerical value of corresponding element in angular displacement of support vector V, then to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation applying is just taken from steady temperature change vector S, to A
_{o}in Cable Structure apply and obtain the current initial mechanical calculating benchmark model A that upgrades after temperature variation
^{i} _{o}, upgrade A
^{i} _{o}time, T
^{i} _{o}all elements numerical value is also used T
^{i}corresponding replacement of all elements numerical value, upgraded T
^{i} _{o}, so just obtained correctly corresponding to A
^{i} _{o}t
^{i} _{o}; D now
^{i} _{o}remain unchanged; When upgrading A
^{i} _{o}after, A
^{i} _{o}the current initial damage of the cable system vector d for health status of rope
^{i} _{o}represent A
^{i} _{o}current cable structure steady temperature data vector T for Cable Structure steady temperature
^{i}represent A
^{i} _{o}current initial Cable Structure bearing angular coordinate vector U for bearing angular coordinate
^{i} _{o}represent.Upgrade C
^{i} _{o}method be: when upgrading A
^{i} _{o}after, A
^{i} _{o}the current initial damage of the evaluation object vector d for health status of evaluation object
^{i} _{o}represent A
^{i} _{o}current cable structure steady temperature data vector T for Cable Structure steady temperature
^{i}represent A
^{i} _{o}current initial Cable Structure bearing angular coordinate vector U for bearing angular coordinate
^{i} _{o}represent, upgrade C
^{i} _{o}method be: when upgrading A
^{i} _{o}after, by Mechanics Calculation, obtain A
^{i} _{o}in concrete numerical value all monitored amounts, current, these concrete numerical value form C
^{i} _{o};
The 8th step: at current initial mechanical calculating benchmark model A
^{i} _{o}basis on, according to step a, to steps d, carry out several times Mechanics Calculation, by calculating, set up unit damage monitored numerical quantity transformation matrices Δ C
^{i}with evaluation object unit change vector D
^{i} _{u}.
A. when the i time circulation starts, directly press step b to the listed method acquisition of steps d Δ C
^{i}and D
^{i} _{u}; At other constantly, when in the 7th step to A
^{i} _{o}after upgrading, must regain Δ C to the listed method of steps d by step b
^{i}and D
^{i} _{u}if, in the 7th step not to A
^{i} _{o}upgrade, directly proceed to herein the 9th step and carry out followup work.
B. at current initial mechanical calculating benchmark model A
^{i} _{o}basis on carry out several times Mechanics Calculation, vectorial d
^{i} _{o}represent A
^{i} _{o}the health status of evaluation object, on calculation times numerical value, equal the quantity N of all evaluation objects, have N evaluation object just to have N calculating; Calculate each time hypothesis and only have an evaluation object at vectorial d
^{i} _{o}on the basis of the health status of the evaluation object representing, there is unit damage or centrepoint load unit change, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable is at vectorial d
^{i} _{o}on the basis of the existing damage of this support cable representing, there is again unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage), if this evaluation object is a centrepoint load, just suppose that this centrepoint load is at vectorial d
^{i} _{o}(if this centrepoint load is couple, centrepoint load unit change can be got 1kNm, 2kNm, 3kNm etc. for unit change on the basis of the existing variable quantity of this centrepoint load representing, to increase centrepoint load unit change again; If this centrepoint load is concentrated force, centrepoint load unit change can be got 1kN, 2kN, 3kN etc. for unit change), use D
^{i} _{uk}record this unit damage or centrepoint load unit change, wherein k represents to occur the numbering of the evaluation object of unit damage or centrepoint load unit change, D
^{i} _{uk}evaluation object unit change vector D
^{i} _{u}an element, evaluation object unit change vector D
^{i} _{u}coding rule and the vectorial d of element
_{o}the coding rule of element identical; The evaluation object that occurs unit damage or centrepoint load unit change in calculating is each time different from the evaluation object that occurs unit damage or centrepoint load unit change in other calculating, calculate each time the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure, the current calculated value of all monitored amounts that calculate each time forms a monitored amount calculation current vector; When k evaluation object of hypothesis has unit damage or centrepoint load unit change, available C
^{i} _{tk}represent corresponding " monitored amount calculation current vector "; While giving in this step each vectorial element numbering, should use same coding rule with other vector in this method, to guarantee any one element in each vector in this step, with in other vector, number identical element, expressed the relevant information of same monitored amount or same target; C
^{i} _{tk}definition mode and vectorial C
_{o}definition mode identical, C
^{i} _{tk}element and C
_{o}element corresponding one by one.
C. the vectorial C calculating each time
^{i} _{tk}deduct vectorial C
^{i} _{o}obtain a vector, then unit damage or the centrepoint load unit change numerical value D of supposition during each element of this vector is calculated divided by this
^{i} _{uk}after obtain " the numerical value change of a monitored amount vector δ C
^{i} _{k}"; There is N evaluation object just to have N " the numerical value change vector of monitored amount ".
D. by this N " the numerical value change vector of monitored amount ", according to the coding rule of N evaluation object, form successively " the unit damage monitored numerical quantity transformation matrices Δ C that has N row
^{i}"; Unit damage monitored numerical quantity transformation matrices Δ C
^{i}each row corresponding to a monitored amount unit change vector; Unit damage monitored numerical quantity transformation matrices Δ C
^{i}every a line corresponding to same monitored amount the different unit change amplitude when different evaluation objects increase unit damage or centrepoint load unit change; Unit damage monitored numerical quantity transformation matrices Δ C
^{i}coding rule and the vectorial d of row
_{o}the coding rule of element identical, unit damage monitored numerical quantity transformation matrices Δ C
^{i}the coding rule of coding rule and M monitored amount of row identical.
The 9th step: set up linear relationship error vector e
^{i}with vectorial g
^{i}.Utilize data (" the current initial value vector of monitored amount C above
^{i} _{o}", " unit damage monitored numerical quantity transformation matrices Δ C
^{i}"); when the 8th step is calculated each time; when only having the increase unit damage or centrepoint load unit change of an evaluation object in calculating each time hypothesis evaluation object; when hypothesis k(k=1,2,3; ..., when N) individual evaluation object increases unit damage or centrepoint load unit change, calculate each time and form a damage vector, use d
^{i} _{tk}represent this damage vector, corresponding monitored amount calculation current vector is C
^{i} _{tk}(referring to the 8th step), damages vectorial d
^{i} _{tk}element number equal the quantity of evaluation object, vectorial d
^{i} _{tk}all elements in only have the numerical value of an element to get to calculate each time in hypothesis increase unit damage or the centrepoint load unit change value of the evaluation object of unit damage or centrepoint load unit change, d
^{i} _{tk}the numerical value of other element get 0, that is not the numbering of 0 element and corresponding relation that supposition increases the evaluation object of unit damage or centrepoint load unit change, with the element of the same numbering of other vectors, with the corresponding relation of this evaluation object, is identical; d
^{i} _{tk}with evaluation object initial damage vector d
_{o}element coding rule identical, d
^{i} _{tk}element and d
_{o}element be onetoone relationship.By C
^{i} _{tk}, C
^{i} _{o}, Δ C
^{i}, d
^{i} _{tk}bring formula (34) into, obtain a linear relationship error vector e
^{i} _{k}, calculate each time a linear relationship error vector e
^{i} _{k}; e
^{i} _{k}subscript k represent k(k=1,2,3 ..., N) individual evaluation object increases unit damage or centrepoint load unit change.There is N evaluation object just to have N calculating, just have N linear relationship error vector e
^{i} _{k}, by this N linear relationship error vector e
^{i} _{k}after addition, obtain a vector, the new vector that each element of this vector is obtained after divided by N is exactly final linear relationship error vector e
^{i}.Vector g
^{i}equal final error vector e
^{i}.By vectorial g
^{i}be kept on the hard disc of computer of operation health monitoring systems software, for health monitoring systems software application.
The tenth step: define the vectorial d of current name damage
^{i} _{c}with current actual damage vector d
^{i}, d
^{i} _{c}and d
^{i}element number equal the quantity of evaluation object, d
^{i} _{c}and d
^{i}element and evaluation object between be onetoone relationship, d
^{i} _{c}and d
^{i}element numerical value represent degree of injury or the centrepoint load intensity of variation of corresponding evaluation object, d
^{i} _{c}and d
^{i}with evaluation object initial damage vector d
_{o}element coding rule identical, d
^{i} _{c}element, d
^{i}element and d
_{o}element be onetoone relationship.
The 11 step: according to monitored amount current value vector C
^{i}with " the current initial value vector of monitored amount C
^{i} _{o}", " unit damage monitored numerical quantity transformation matrices Δ C
^{i}" and " the vectorial d of current name damage
^{i} _{c}" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula (11), according to multiobjective optimization algorithm, calculates the vectorial d of current name damage
^{i} _{c}noninferior solution, namely with reasonable error but can be more exactly determine the position of damaged cable and the solution of nominal degree of injury thereof from all ropes.
The multiobjective optimization algorithm that can adopt has a variety of, for example: the multipleobjection optimization based on genetic algorithm, the multipleobjection optimization based on artificial neural network, the multiobjective optimization algorithm based on population, the multipleobjection optimization based on ant group algorithm, leash law (Constrain Method), weighted method (Weighted SUm Method), Objective Programming (Goal Attainment Method) etc.Because various multiobjective optimization algorithms are all conventional algorithms, can realize easily, this implementation step only be take Objective Programming and as example provides, is solved the vectorial d of current name damage
^{i} _{c}process, the specific implementation process of other algorithm can realize in a similar fashion according to the requirement of its specific algorithm.
According to Objective Programming, formula (11) can transform the multiobjective optimization question shown in an accepted way of doing sth (35) and formula (36), and in formula (35), γ is a real number, and R is real number field, and area of space Ω has limited vectorial d
^{i} _{c}span (the present embodiment requirements vector d of each element
^{i} _{c}each element be not less than 0, be not more than 1).The meaning of formula (35) is to find a minimum real number γ, and formula (36) is met.G (d in formula (36)
^{i} _{c}) by formula (36) definition, G (d in the product representation formula (36) of weighing vector W and γ in formula (36)
^{i} _{c}) and vectorial g
^{i}between the deviation that allows, g
^{i}definition referring to formula (17), its value calculates in the 9th step.During actual computation vector W can with vectorial g
^{i}identical.The concrete programming of Objective Programming realizes has had universal program directly to adopt.Use Objective Programming just can damage vectorial d in the hope of current name
^{i} _{c}.
The 12 step: according to the current actual damage vector of cable system d
^{i}definition (seeing formula (18)) and the definition (seeing formula (19)) of its element calculate current actual damage vector d
^{i}each element, thereby can be by d
^{i}determine the health status of evaluation object.Current actual damage vector d
^{i}k element d
^{i} _{k}the current actual health status that represents k evaluation object in the i time circulation.
D
^{i} _{k}the current actual health status that represents k evaluation object in the i time circulation, if this evaluation object is support cable, so a d in cable system
^{i} _{k}represent its current actual damage, d
^{i} _{k}be to represent that its corresponding support cable is without health problem at 0 o'clock, d
^{i} _{k}numerical value is not the support cable that represents that its corresponding support cable is unsoundness problem at 0 o'clock, and the support cable of unsoundness problem may be slack line, also may be damaged cable, its numerical response the degree of lax or damage.
D
^{i} _{k}the current actual health status that represents k evaluation object in the i time circulation, if this evaluation object is centrepoint load, so a d
^{i} _{k}represent that its current actual centrepoint load changes numerical value.
By the current actual damage vector of evaluation object d
^{i}in the M relevant to support cable
_{1}individual element takes out, and forms the current actual damage vector of support cable d
^{ci}, the current actual damage vector of support cable d
^{ci}the coding rule and initial rope force vector F of element
_{o}the coding rule of element identical.The current actual damage vector of support cable d
^{ci}h element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., M
_{1}; The current actual damage vector of support cable d
^{ci}middle numerical value be not 0 element corresponding to the support cable of unsoundness problem, from the support cable of these unsoundness problems, identify damaged cable, remaining is exactly slack line.Mirror method for distinguishing is varied; can be by removing the protective seam of the support cable of unsoundness problem; support cable is carried out to visual discriminating; or carry out visual discriminating by optical imaging apparatus; also can to whether support cable is impaired, differentiate by lossless detection method, UT (Ultrasonic Testing) is exactly a kind of now widely used lossless detection method.After differentiating, those do not find that support cable damage, unsoundness problem is exactly that lax rope has occurred, and need exactly to adjust the rope of Suo Li, are exactly slack line, and these ropes that need adjust Suo Li are at the current actual damage vector of support cable d
^{ci}in corresponding element numerical value (for example one of them element can be used d
^{ci} _{h}expression) degree of injury of the relax level mechanics equivalence of expression and these support cables, has just determined slack line thus, and the computing method of concrete slack illustrate below.The current actual damage vector of the support cable corresponding with damaged cable d
^{ci}in element numerical expression be the current actual damage of this damaged cable, element numerical value represents while being 100% that this support cable thoroughly loses loadbearing capacity, represents to lose the loadbearing capacity of corresponding proportion in the time of between 0 and 100%; According to the current actual damage vector of support cable d
^{ci}, from the support cable of unsoundness problem, identify after slack line, remaining is exactly damaged cable, and these damaged cables are at the current actual damage vector of support cable d
^{ci}the numerical value of the element of middle correspondence just represents its degree of injury, the numerical value of corresponding element represents while being 100% that this support cable thoroughly loses loadbearing capacity, the loadbearing capacity that represents this support cable forfeiture corresponding proportion in the time of between 0 and 100%, has so far just identified damaged cable and degree of injury thereof.
Utilize the current actual damage vector of support cable d
^{ci}slack line and the current actual equivalent damage degree its relax level mechanics equivalence of expressing, utilize in the 6th step, obtain at current cable structure steady temperature data vector T
^{i}current cable force vector F under condition
^{i}with current support cable two supporting end points horizontal range vector l
^{i} _{x}, utilize at second step, obtain at initial Cable Structure steady temperature data vector T
_{o}the initial drift vector l of the support cable under condition
_{o}, initial free crosssectional area vector A
_{o}weight vector ω with initial free unit length
_{o}, utilize current cable structure steady temperature data vector T
^{i}the current steady temperature data of support cable that represent, utilize at second step, obtain at initial Cable Structure steady temperature data vector T
_{o}the support cable initial steady state temperature data representing, the temperature variant physical and mechanical properties parameter of the various materials that the Cable Structure that utilization obtains at second step is used, count the impact of temperature variation on support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanics equivalence calculate slack line, with the relax level of current actual equivalent damage degree equivalence, mechanics equivalent condition is: one, the mechanics parameters without lax initial drift, geometrical property parameter, density and material during with not damaged of the rope of two equivalences is identical; Two, after lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after overall length identical.While meeting abovementioned two mechanics equivalent conditions, the mechanics function of two such support cables in Cable Structure is exactly identical, if replaced after damaged cable with equivalent slack line, Cable Structure any variation can not occur, and vice versa.According to aforementioned mechanics equivalent condition, try to achieve the relax level that those are judged as slack line, relax level is exactly the change amount of support cable drift, has namely determined the long adjustment amount of rope of the support cable that those need adjust Suo Li.Particularly can be in the hope of the relax level (being the long adjustment amount of rope) of these ropes according to formula (29) or formula (30).So just realized the lax identification of support cable.So far damaged cable and slack line have just all been identified.
So according to the current actual damage vector of evaluation object d
^{i}can problem identificatioin rope and health degree thereof, can define which centrepoint load variation and numerical value thereof have occurred simultaneously.
So far, can say that this method has realized two kinds of functions that existing method can not possess, be respectively, when the centrepoint load of one, bearing when bearing has angular displacement, in structure and structure (environment) temperature changes simultaneously, can reject angular displacement of support, centrepoint load variation and structure temperature and change the impact on Cable Structure health status recognition result, thereby identify exactly the structure health monitoring method of problem rope; Two, this method, when identifying problem rope, can also identify the variation of centrepoint load simultaneously, and this method can be rejected angular displacement of support, structure temperature changes and the impact of support cable health status variation, realizes the correct identification of centrepoint load intensity of variation.
The 13 step: the computing machine in health monitoring systems regularly generates cable system health condition form automatically or by personnel's operational health monitoring system.
The 14 step: under specified requirements, computing machine automatic operation in health monitoring systems communication panalarm to monitor staff, owner and (or) personnel of appointment report to the police.
The 15 step: set up mark vector B according to formula (31)
^{i}, formula (32) has provided mark vector B
^{i}the definition of k element; If mark vector B
^{i}element be 0 entirely, get back to the 6th step and proceed the health monitoring of cable system and calculating; If mark vector B
^{i}element be not 0 entirely, complete after subsequent step, enter next time circulation.
The 16 step: according to formula (33) calculate next time (the i+1 time, i=1,2,3,4 ...) circulate initial damage vector d required
^{i+1} _{o}each element d
^{i+1} _{ok}(k=1,2,3 ..., N); The second, at initial mechanical calculating benchmark model A
_{o}basis on, first to A
_{o}in Cable Structure bearing apply angular displacement of support constraint, the numerical value of angular displacement of support constraint is just taken from the numerical value of corresponding element in angular displacement of support vector V, then to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation applying is just taken from steady temperature change vector S, then to make the health status of rope be d
^{i+1} _{o}after obtain be exactly next time, the i+1 time (i=1,2,3,4 ...) circulate Mechanics Calculation benchmark model A required
^{i+1}; Next time (the i+1 time, i=1,2,3,4 ...) required current initial Cable Structure steady temperature data vector T circulates
^{i+1} _{o}equal T
^{i} _{o}, next time (the i+1 time, i=1,2,3,4 ...) required current initial Cable Structure bearing angular coordinate vector U circulates
^{i+1} _{o}equal U
^{i} _{o}.Obtain A
^{i+1}, d
^{i+1} _{o}, U
^{i+1} _{o}and T
^{i+1} _{o}after, by Mechanics Calculation, obtain A
^{i+1}in concrete numerical value all monitored amounts, current, these concrete numerical value form next time, the vectorial C of the current initial value of required monitored amount that circulates for the i+1 time
^{i+1} _{o}.
The 17 step: get back to the 6th step, start the circulation by the 6th step to the 17 steps.
Claims (1)
1. the laddering recognition methods of problem rope centrepoint load of hybrid monitoring during angular displacement, is characterized in that described method comprises:
A. for sake of convenience, it is evaluation object that this method unitedly calls evaluated support cable and centrepoint load, establishes the quantity of evaluated support cable and the quantity sum of centrepoint load is N, and the quantity of evaluation object is N; Determine the coding rule of evaluation object, by this rule, by evaluation object numberings all in Cable Structure, this numbering will be for generating vector sum matrix in subsequent step; This method represents this numbering with variable k, k=1, and 2,3 ..., N; While determining hybrid monitoring appointment by the support cable of monitored Suo Li, establish in cable system total Q root support cable, the monitored rope force data of Cable Structure is by M in Cable Structure
_{1}the M of individual appointment support cable
_{1}individual rope force data is described, and the variation of Cable Structure Suo Li is exactly the variation of the Suo Li of all appointment support cables; Each total M
_{1}individual cable force measurement value or calculated value characterize the rope force information of Cable Structure; M
_{1}be one and be not less than 0 integer that is not more than Q; While determining hybrid monitoring appointment by the measured point of monitored strain, the monitored strain data of Cable Structure can be by K in Cable Structure
_{2}l individual specified point and each specified point
_{2}the strain of individual assigned direction is described, and the variation of Cable Structure strain data is exactly K
_{2}the variation of all tested strains of individual specified point; Each total M
_{2}individual strain measurement value or calculated value characterize Cable Structure strain, M
_{2}for K
_{2}and L
_{2}longpending; M
_{2}to be not less than 0 integer; While determining hybrid monitoring appointment by the measured point of monitored angle, the monitored angledata of Cable Structure is by K in Cable Structure
_{3}l individual specified point, that cross each specified point
_{3}h individual appointment straight line, each appointment straight line
_{3}individual angle coordinate component is described, and the variation of Cable Structure angle is exactly the variation of angle coordinate components appointment straight lines all specified points, all, all appointments; Each total M
_{3}individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M
_{3}for K
_{3}, L
_{3}and H
_{3}longpending; M
_{3}be one and be not less than 0 integer; While determining hybrid monitoring appointment by monitored shape data, the monitored shape data of Cable Structure is by K in Cable Structure
_{4}l individual specified point and each specified point
_{4}the volume coordinate of individual assigned direction is described, and the variation of Cable Structure shape data is exactly K
_{4}the variation of all coordinate components of individual specified point; Each total M
_{4}individual measurement of coordinates value or calculated value characterize Cable Structure shape, M
_{4}for K
_{4}and L
_{4}longpending; M
_{4}be one and be not less than 0 integer; The monitored amount of comprehensive abovementioned hybrid monitoring, total M the monitored amount of whole Cable Structure, M is M
_{1}, M
_{2}, M
_{3}and M
_{4}sum, definition parameter K, K is M
_{1}, K
_{2}, K
_{3}and K
_{4}sum, K and M must not be less than the quantity N of evaluation object; For simplicity, in the method this is walked to listed M monitored amount referred to as " monitored amount "; In this method, to the time interval between any twice measurement of same amount RealTime Monitoring, must not be greater than 30 minutes, the moment of survey record data is called physical record data constantly;
B. this method definition " the temperature survey calculating method of the Cable Structure of this method " is undertaken by step b1 to b3;
B1: inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of Cable Structure composition material and Cable Structure environment of living in, utilize the geometry measured data of design drawing, asconstructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure, inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, statistics obtains interior during this period of time cloudy quantity and is designated as T cloudy day, in the method can not be seen to one of the sun daytime and be called all day the cloudy day, statistics obtains 0 highest temperature and the lowest temperature between latter 30 minutes of the moment of sunrise next day at each cloudy day in T cloudy day, sunrise on the meteorology that sunrise refers to base area revolutions constantly and the rule that revolves round the sun is definite constantly, do not represent necessarily can see the sun same day, the sunrise that can inquire about data or calculate each required day by conventional meteorology constantly, each cloudy day 0 up to next day sunrise constantly the highest temperature between latter 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, there is T cloudy day, the maximum temperature difference that just has the daily temperature at T cloudy day, get maximal value in the maximum temperature difference of daily temperature at T cloudy day for reference to temperature difference per day, with reference to temperature difference per day, be designated as Δ T
_{r}, between inquiry Cable Structure location and Altitude Region, place, be no less than temperature that the meteorological data in recent years of 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and Altitude Region, place about the maximum rate of change Δ T of sea level elevation
_{h}, for convenience of narration, get Δ T
_{h}unit be ℃/m, on the surface of Cable Structure, get " R Cable Structure surface point ", get the Specific Principles of " R Cable Structure surface point " narrates in step b3, after will by actual measurement, obtain the temperature of this R Cable Structure surface point, claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, by Calculation of Heat Transfer, obtain the temperature of this R Cable Structure surface point, just claim that the temperature data calculating is " R Cable Structure surface temperature computational data ", from the residing minimum height above sea level of Cable Structure to the highest height above sea level, in Cable Structure, uniform choosing is no less than three different sea level elevations, the sea level elevation place choosing at each, at the intersection place on surface level and Cable Structure surface, at least choose two points, outer normal from selected point straw line body structure surface, all outer normal directions of choosing are called " measuring Cable Structure along the direction of the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in the measurement Cable Structure of choosing along comprising the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure in the direction of the Temperature Distribution of wall thickness, direction uniform choosing in Cable Structure along each measurement Cable Structure along the Temperature Distribution of wall thickness is no less than three points, measure all temperature that are selected a little, the temperature recording is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure Cable Structure along the direction of the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if chosen H different sea level elevation, at each sea level elevation place, choose B and measured Cable Structure along the direction of the Temperature Distribution of wall thickness, along each, measure Cable Structure and in Cable Structure, chosen E point along the direction of the Temperature Distribution of wall thickness, wherein H and E are not less than 3, B is not less than 2, if HBE is the product of H and B and E, corresponding total HBE " measuring Cable Structure along the point of the temperature profile data of thickness ", after will obtain the temperature that this HBE " measures Cable Structure along the point of the temperature profile data of thickness " by actual measurement, claim that the temperature data that actual measurement obtains is " HBE Cable Structure is along thickness temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, by Calculation of Heat Transfer, obtain this HBE and measure Cable Structure along the temperature of the point of the temperature profile data of thickness, just claim that the temperature data calculating is " HBE Cable Structure is along thickness temperature computation data ", if BE is the product of B and E, total BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in sea level elevation place of choosing at each in this method, in Cable Structure location, according to meteorology, measure temperature and require to choose a position, will obtain meeting the temperature that meteorology is measured the Cable Structure place environment of temperature requirement in this position actual measurement, in the onsite spaciousness of Cable Structure, unobstructed place chooses a position, this position should can obtain in each day of the whole year this ground the most sufficient sunshine of getable this day, flat board at a carbon steel material of this position of sound production, be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, the one side of this reference plate on the sunny side, be called sunny slope, the sunny slope of reference plate is coarse and dark color, the sunny slope of reference plate should can obtain in each day of the whole year one flat plate on this ground the most sufficient sunshine of getable this day, the nonsunny slope of reference plate is covered with insulation material, RealTime Monitoring is obtained to the temperature of the sunny slope of reference plate,
B2: RealTime Monitoring obtains R Cable Structure surface temperature measured data of abovementioned R Cable Structure surface point, RealTime Monitoring obtains previously defined Cable Structure along the temperature profile data of thickness simultaneously, and RealTime Monitoring obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement simultaneously, by RealTime Monitoring, obtain being carved at sunrise the same day sunrise next day temperature measured data sequence of the Cable Structure place environment between latter 30 minutes constantly, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data that was carved at sunrise the Cable Structure place environment between latter 30 minutes of the moment of sunrise next day the same day, find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment, by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, deduct and be carved at sunrise the sunrise next day maximum temperature difference between latter 30 minutes constantly on same day that minimum temperature obtains Cable Structure place environment, be called environment maximum temperature difference, be designated as Δ T
_{emax}, by the temperature measured data sequence of Cable Structure place environment, by conventional mathematical computations, obtain the temperature of Cable Structure place environment about the rate of change of time, this rate of change is also along with the time changes, by RealTime Monitoring, obtain being carved at sunrise the same day sunrise next day measured data sequence of the temperature of the sunny slope of the reference plate between latter 30 minutes constantly, the measured data sequence of the temperature of the sunny slope of reference plate by be carved at sunrise the same day next day sunrise constantly the measured data of the temperature of the sunny slope of the reference plate between latter 30 minutes according to time order and function order, arrange, find maximum temperature and minimum temperature in the measured data sequence of temperature of sunny slope of reference plate, by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, deduct and be carved at sunrise the sunrise next day maximum temperature difference between latter 30 minutes constantly on same day of temperature that minimum temperature obtains the sunny slope of reference plate, be called reference plate maximum temperature difference, be designated as Δ T
_{pmax}, by RealTime Monitoring, obtain being carved at sunrise the same day sunrise next day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between latter 30 minutes constantly, there is R Cable Structure surface point just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence by be carved at sunrise on same day of a Cable Structure surface point sunrise next day constantly the Cable Structure surface temperature measured data between latter 30 minutes according to time order and function order, arrange, find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence, by the maximum temperature in each Cable Structure surface temperature measured data sequence, deduct and be carved at sunrise the sunrise next day maximum temperature difference between latter 30 minutes constantly on same day that minimum temperature obtains the temperature of each Cable Structure surface point, there is R Cable Structure surface point just to have and be carved at sunrise the sunrise next day maximum temperature difference numerical value between latter 30 minutes constantly R the same day, maximal value is wherein called Cable Structure surface maximum temperature difference, be designated as Δ T
_{smax}, by each Cable Structure surface temperature measured data sequence, by conventional mathematical computations, obtain the temperature of each Cable Structure surface point about the rate of change of time, the temperature of each Cable Structure surface point about the rate of change of time also along with the time changes, by RealTime Monitoring, obtain being carved at sunrise the same day between latter 30 minutes of the moment of sunrise next day, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", the sea level elevation place that calculating is chosen at each amounts to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and the difference of minimum temperature, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", chosen H different sea level elevation and just had H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T
_{tmax},
B3: measure and calculate acquisition Cable Structure steady temperature data, first, determine the moment that obtains Cable Structure steady temperature data, the condition relevant to the moment that determines acquisition Cable Structure steady temperature data has six, first condition is the moment that obtains Cable Structure steady temperature data to be carved at sunset sunrise next day constantly between latter 30 minutes between the same day, sunset refers to sunset on base area revolutions and the definite meteorology of revolution rule constantly constantly, and the sunset that can inquire about data or calculate each required day by conventional meteorology constantly, the a condition of second condition be the same day be carved at sunrise next day sunrise constantly between latter 30 minutes during this period of time in, reference plate maximum temperature difference Δ T
_{pmax}with Cable Structure surface maximum temperature difference Δ T
_{smax}all be not more than 5 degrees Celsius, the b condition of second condition be the same day be carved at sunrise next day sunrise constantly between latter 30 minutes during this period of time in, above, measure the environment maximum error Δ T calculating
_{emax}be not more than with reference to temperature difference per day Δ T
_{r}, and reference plate maximum temperature difference Δ T
_{pmax}after deducting 2 degrees Celsius, be not more than Δ T
_{emax}, and Cable Structure surface maximum temperature difference Δ T
_{smax}be not more than Δ T
_{pmax}, only needing to meet in a condition of second and b condition one is just called and meets second condition, the 3rd condition is that the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, the 4th condition is that the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, the 5th condition is in the moment that obtains Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is to be carved at sunrise the sunrise next day minimal value between latter 30 minutes constantly the same day, the 6th condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T
_{tmax}be not more than 1 degree Celsius, this method is utilized abovementioned six conditions, any one in following three kinds of moment is called to " obtaining the mathematics of Cable Structure steady temperature data constantly ", the first is first moment to the 5th condition meeting in abovementioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data " constantly, the second is the moment that only meets the 6th condition in abovementioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data " constantly, the third is first moment to the 6th condition simultaneously meeting in abovementioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data " constantly, when obtaining the mathematics of Cable Structure steady temperature data, be exactly in this method during in constantly one of physical record data constantly, the moment that obtains Cable Structure steady temperature data be exactly obtain Cable Structure steady temperature data mathematics constantly, if obtain the mathematics of Cable Structure steady temperature data and be not constantly any in constantly of physical record data in this method constantly, get this method close to moment of the mathematics that obtains Cable Structure steady temperature data those physical record data constantly for obtaining the moment of Cable Structure steady temperature data, this method is carried out the relevant health monitoring analysis of Cable Structure by using in the amount that obtains the moment survey record of Cable Structure steady temperature data, this method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, i.e. this Cable Structure temperature constantly temporal evolution not, and this is exactly " obtaining the moment of Cable Structure steady temperature data " of this method constantly, then, according to Cable Structure heat transfer characteristic, utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment that obtains Cable Structure steady temperature data, utilize the thermal conduction study computation model of Cable Structure, by conventional Calculation of Heat Transfer, obtain in the Temperature Distribution of Cable Structure that obtains the moment of Cable Structure steady temperature data, now calculate by stable state in the temperature field of Cable Structure, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data, also comprise that Cable Structure is in the accounting temperature of selected HBE " measuring Cable Structure along the point of the temperature profile data of thickness " above, the accounting temperature of HBE " measuring Cable Structure along the point of the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", when R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure stable state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", while getting " R Cable Structure surface point " on the surface of Cable Structure, the quantity of " R Cable Structure surface point " and necessary three conditions that meet that distribute, first condition is when Cable Structure temperature field is during in stable state, on Cable Structure surface the temperature of any point be by " R Cable Structure surface point " with Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point while obtaining, on the Cable Structure surface that linear interpolation obtains, on the temperature of this arbitrfary point and Cable Structure surface, the error of the actual temperature of this arbitrfary point is not more than 5%, Cable Structure surface comprises support cable surface, second condition is that in " R Cable Structure surface point ", the quantity at the point of same sea level elevation is not less than 4, and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 ℃ divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation
_{h}the numerical value obtaining, gets Δ T for convenience of narration
_{h}unit be ℃/m that the unit of getting Δ h for convenience of narration is m, " R Cable Structure surface point " refers to while only considering sea level elevation along the definition of adjacent Cable Structure surface point between two of sea level elevation, in " R Cable Structure surface point ", do not have a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two, the 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to the geometric properties of Cable Structure and bearing data, in Cable Structure, find and be subject to the sunshineduration position of those surface points the most fully the whole year, in " R Cable Structure surface point ", having a Cable Structure surface point at least is the annual point being subject in the most sufficient those surface points of sunshineduration in Cable Structure,
C. according to " the temperature survey calculating method of the Cable Structure of this method ", directly measure and calculate the Cable Structure steady temperature data under original state, Cable Structure steady temperature data under original state are called initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T
_{o}", actual measurement or consult reference materials and obtain the temperature variant physical and mechanical properties parameter of the various materials that Cable Structure used, in actual measurement, obtain initial Cable Structure steady temperature data vector T
_{o}synchronization, directly measure the initial Suo Li that calculates all support cables, form initial rope force vector F
_{o}, according to the data that comprise Cable Structure design data, completion data obtain that all support cables are in free state that Suo Li is the length of 0 o'clock, the weight of crosssectional area during in free state and the unit length during in free state, and the temperature of all support cables while obtaining these three kinds of data, utilize on this basis temperature variant physical function parameter and the mechanical property parameters of all support cables, according to conventional physical computing, obtain all support cables at initial Cable Structure steady temperature data vector T
_{o}suo Li under condition is that the length of 0 o'clock all support cable, crosssectional area and the Suo Li that Suo Li is 0 o'clock all support cable are the weight of the unit length of 0 o'clock all support cable, form successively the initial drift vector of support cable, the weight vector of the initial free unit length of initial free crosssectional area vector sum, the coding rule and initial rope force vector F of the element of the initial drift vector of support cable, the weight vector of the initial free unit length of initial free crosssectional area vector sum
_{o}the coding rule of element identical, in actual measurement, obtain T
_{o}time, namely obtaining initial Cable Structure steady temperature data vector T
_{o}the synchronization in the moment, directly measure the measured data that calculates initial Cable Structure, the measured data of initial Cable Structure is to comprise Cable Structure centrepoint load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the initial value of all monitored amounts, the initial rope force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing generalized coordinate data, initial Cable Structure angledata, initial Cable Structure spatial data is in interior measured data, initial Cable Structure bearing generalized coordinate data comprise initial Cable Structure bearing spatial data and initial Cable Structure bearing angular data, when obtaining the measured data of initial Cable Structure, measurement calculates the data of the health status that can express support cable of the Nondestructive Testing Data that comprises support cable, and the data of the health status that can express support cable are now called support cable initial health data, the initial value of all monitored amounts forms monitored amount initial value vector C
_{o}, monitored amount initial value vector C
_{o}the coding rule of coding rule and M monitored amount identical, utilize support cable initial health data and Cable Structure centrepoint load measurement data to set up evaluation object initial damage vector d
_{o}, vectorial d
_{o}represent with initial mechanical calculating benchmark model A
_{o}the initial health of the evaluation object of the Cable Structure representing, evaluation object initial damage vector d
_{o}element number equal N, d
_{o}element and evaluation object be onetoone relationship, vectorial d
_{o}the coding rule of element identical with the coding rule of evaluation object, if d
_{o}evaluation object corresponding to some elements be support cable, so a d in cable system
_{o}the numerical value of this element represent the initial damage degree of corresponding support cable, if the numerical value of this element is 0, represent that the corresponding support cable of this element is intact, do not damage, if its numerical value is 100%, represent that the corresponding support cable of this element has completely lost loadbearing capacity, if its numerical value between 0 and 100%, represents this support cable, lost the loadbearing capacity of corresponding proportion, if d
_{o}evaluation object corresponding to some elements be some centrepoint load, in this method, get d
_{o}this element numerical value be 0, the initial value that represents the variation of this centrepoint load is 0, if while there is no the Nondestructive Testing Data of support cable and the data of other health status that can express support cable, or can think that structure original state is not damaged during without relaxed state, vectorial d
_{o}in each element numerical value relevant to support cable get 0, initial Cable Structure bearing angular data forms initial Cable Structure bearing angular coordinate vector U
_{o},
The temperature variant physical and mechanical properties parameter of the various materials that d. use according to measured data, support cable initial health data, Cable Structure centrepoint load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the Cable Structure of the design drawing of Cable Structure, asconstructed drawing and initial Cable Structure, initial Cable Structure bearing angular coordinate vector U
_{o}, initial Cable Structure steady temperature data vector T
_{o}with all Cable Structure data that preceding step obtains, set up the initial mechanical calculating benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data "
_{o}, based on A
_{o}the Cable Structure computational data calculating must approach its measured data very much, and difference therebetween must not be greater than 5%; Corresponding to A
_{o}" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T
_{o}"; Corresponding to A
_{o}cable Structure bearing angular data be exactly initial Cable Structure bearing angular coordinate vector U
_{o}; Corresponding to A
_{o}evaluation object initial damage vector d for evaluation object health status
_{o}represent; Corresponding to A
_{o}monitored amount initial value vector C for the initial value of all monitored amounts
_{o}represent; U
_{o}, T
_{o}and d
_{o}a
_{o}parameter, by A
_{o}initial value and the C of all monitored amounts of obtaining of Mechanics Calculation result
_{o}the initial value of all monitored amounts that represent is identical, therefore also can say C
_{o}by A
_{o}mechanics Calculation result form, A in the method
_{o}, C
_{o}, d
_{o}, U
_{o}and T
_{o}constant;
E. in the method, alphabetical i, except representing that significantly, the place of number of steps, alphabetical i only represents cycle index, circulates for the i time; When the i time circulation starts, the current initial mechanical calculating benchmark model of Cable Structure that need to set up or that set up is designated as current initial mechanical calculating benchmark model A
^{i} _{o}, A
_{o}and A
^{i} _{o}count temperature parameter, can accounting temperature change the Effect on Mechanical Properties to Cable Structure; When the i time circulation starts, corresponding to A
^{i} _{o}" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T
^{i} _{o}represent vector T
^{i} _{o}definition mode and vector T
_{o}definition mode identical, T
^{i} _{o}element and T
_{o}element corresponding one by one; When the i time circulation starts, corresponding to A
^{i} _{o}" Cable Structure bearing angular data " with current initial Cable Structure bearing angular coordinate vector U
^{i} _{o}represent vectorial U
^{i} _{o}definition mode and vectorial U
_{o}definition mode identical, Ui
_{o}element and U
_{o}element corresponding one by one; The current initial damage vector of evaluation object that the i time circulation needs while starting is designated as d
^{i} _{o}, d
^{i} _{o}cable Structure A while representing this circulation beginning
^{i} _{o}the health status of evaluation object, d
^{i} _{o}definition mode and d
_{o}definition mode identical, d
^{i} _{o}element and d
_{o}element corresponding one by one; When the i time circulation starts, the initial value of all monitored amounts, with the current initial value vector of monitored amount C
^{i} _{o}represent vectorial C
^{i} _{o}definition mode and vectorial C
_{o}definition mode identical, C
^{i} _{o}element and C
_{o}element corresponding one by one, the current initial value vector of monitored amount C
^{i} _{o}expression is corresponding to A
^{i} _{o}the concrete numerical value of all monitored amounts; U
^{i} _{o}, T
^{i} _{o}and d
^{i} _{o}a
^{i} _{o}characterisitic parameter, C
^{i} _{o}by A
^{i} _{o}mechanics Calculation result form; When circulation starts for the first time, A
^{i} _{o}be designated as A
^{1} _{o}, set up A
^{1} _{o}method for making A
^{1} _{o}equal A
_{o}; When circulation starts for the first time, T
^{i} _{o}be designated as T
^{1} _{o}, set up T
^{1} _{o}method for making T
^{1} _{o}equal T
_{o}; When circulation starts for the first time, U
^{i} _{o}be designated as U
^{1} _{o}, set up U
^{1} _{o}method for making U
^{1} _{o}equal U
_{o}; When circulation starts for the first time, d
^{i} _{o}be designated as d
^{1} _{o}, set up d
^{1} _{o}method for making d
^{1} _{o}equal d
_{o}; When circulation starts for the first time, C
^{i} _{o}be designated as C
^{1} _{o}, set up C
^{1} _{o}method for making C
^{1} _{o}equal C
_{o};
F. from entering the circulation that is walked q step by f here; In structure military service process, according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement, obtain the current data of Cable Structure steady temperature data, the current data of all " Cable Structure steady temperature data " forms current cable structure steady temperature data vector T
^{i}, vector T
^{i}definition mode and vector T
_{o}definition mode identical, T
^{i}element and T
_{o}element corresponding one by one; In actual measurement, obtain current cable structure steady temperature data vector T
^{i}synchronization, actual measurement obtains Cable Structure bearing angular coordinate current data, all Cable Structure bearing angular coordinate current datas form current cable structures actual measurement bearing angular coordinate vector U
^{i}, vectorial U
^{i}definition mode and vectorial U
_{o}definition mode identical; In actual measurement, obtain vector T
^{i}time, actual measurement obtains obtaining current cable structure steady temperature data vector T
^{i}the Cable Structure of synchronization in the moment in the currency of all monitored amounts, all these numerical value form monitored amount current value vector C
^{i}, vectorial C
^{i}definition mode and vectorial C
_{o}definition mode identical, C
^{i}element and C
_{o}element corresponding one by one, represent that identical monitored amount is at numerical value in the same time not; In actual measurement, obtain the synchronization of current cable structure steady temperature data vector Ti, actual measurement obtains all M in Cable Structure
_{1}the rope force data of root support cable, all these rope force datas form current cable force vector F
^{i}, vectorial F
^{i}element and vectorial F
_{o}the coding rule of element identical; In actual measurement, obtain current cable structure steady temperature data vector T
^{i}synchronization, Actual measurement obtains all M
_{1}the volume coordinate of two supporting end points of root support cable, the volume coordinate of two the supporting end points in the horizontal direction difference of component is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables form current support cable two supporting end points horizontal range vectors, the coding rule and initial rope force vector F of the element of current support cable two supporting end points horizontal range vectors
_{o}the coding rule of element identical;
G. according to current cable structure actual measurement bearing angular coordinate vector U
^{i}with current cable structure steady temperature data vector T
^{i}, according to step g 1 to g3, upgrade current initial mechanical calculating benchmark model A
^{io}, the current initial value of monitored amount vector C
^{i} _{o}, current initial Cable Structure steady temperature data vector T
^{i} _{o}with current initial Cable Structure bearing angular coordinate vector U
^{i} _{o}, and the current initial damage vector of evaluation object d
^{i} _{o}remain unchanged;
G1. compare respectively T
^{i}and T
^{i} _{o}, U
^{i}and U
^{i} _{o}if, T
^{i}equal T
^{i} _{o}and U
^{i}equal U
^{i} _{o}, do not need A
^{i} _{o}upgrade, otherwise need to follow these steps to A
^{i} _{o}, U
^{i} _{o}and T
^{i} _{o}upgrade;
G2. calculate U
^{i}with U
_{o}poor, U
^{i}with U
_{o}difference be exactly Cable Structure bearing about the angular displacement of support of initial position, with angular displacement of support vector V, represent angular displacement of support, V equals U
^{i}deduct U
_{o}; Calculate T
^{i}with T
_{o}poor, T
^{i}with T
_{o}difference be exactly that current cable structure steady temperature data are about the variation of initial Cable Structure steady temperature data, T
^{i}with T
_{o}poor with steady temperature change vector S, represent, S equals T
^{i}deduct T
_{o}, S represents the variation of Cable Structure steady temperature data;
G3. first to A
_{o}in Cable Structure bearing apply angular displacement of support constraint, the numerical value of angular displacement of support constraint is just taken from the numerical value of corresponding element in angular displacement of support vector V, then to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation applying is just taken from steady temperature change vector S, to A
_{o}middle Cable Structure bearing applies angular displacement of support constraint and Cable Structure is applied and obtains the current initial mechanical calculating benchmark model A that upgrades after temperature variation
^{i} _{o}, upgrade A
^{i} _{o}time, U
^{i} _{o}all elements numerical value is also used U
^{i}all elements numerical value is corresponding to be replaced, and has upgraded U
^{i} _{o}, T
^{i} _{o}all elements numerical value is also used T
^{i}corresponding replacement of all elements numerical value, upgraded T
^{i} _{o}, so just obtained correctly corresponding to A
^{i} _{o}u
^{i} _{o}and T
^{i} _{o}, d now
^{i} _{o}remain unchanged; When upgrading A
^{i} _{o}after, A
^{i} _{o}the current initial damage of the evaluation object vector d for health status of rope
^{i} _{o}represent A
^{i} _{o}current cable structure steady temperature data vector T for Cable Structure steady temperature
^{i} _{o}represent A
^{i} _{o}current initial Cable Structure bearing angular coordinate vector U for bearing angular coordinate
^{i} _{o}represent; Upgrade C
^{i} _{o}method be: when upgrading A
^{i} _{o}after, by Mechanics Calculation, obtain A
^{i} _{o}in concrete numerical value all monitored amounts, current, these concrete numerical value form C
^{i} _{o};
H. at current initial mechanical calculating benchmark model A
^{i} _{o}basis on, according to step h1, carry out several times Mechanics Calculation to step h4, by calculating, set up unit damage monitored numerical quantity transformation matrices Δ C
^{i}with evaluation object unit change vector D
^{i} _{u};
H1. when the i time circulation starts, directly press step h2 to the listed method acquisition of step h4 Δ C
^{i}and D
^{i} _{u}; At other constantly, when in step g to A
^{i} _{o}after upgrading, must regain Δ C to the listed method of step h4 by step h2
^{i}and D
^{i} _{u}if, in step g not to A
^{i} _{o}upgrade, directly proceed to herein step I and carry out followup work;
H2. at current initial mechanical calculating benchmark model A
^{i} _{o}basis on carry out several times Mechanics Calculation, on calculation times numerical value, equal the quantity N of all evaluation objects, have N evaluation object just to have N calculating; Coding rule according to evaluation object, calculates successively; Calculating each time hypothesis only has an evaluation object on the basis of original damage or centrepoint load, to increase unit damage or centrepoint load unit change again, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable increases unit damage again, if this evaluation object is a centrepoint load, just suppose that this centrepoint load increases centrepoint load unit change again, uses D
^{i} _{uk}record unit damage or the centrepoint load unit change of this increase, wherein k represents to increase the numbering of the evaluation object of unit damage or centrepoint load unit change, D
^{i} _{uk}evaluation object unit change vector D
^{i} _{u}an element, evaluation object unit change vector D
^{i} _{u}coding rule and the vectorial d of element
_{o}the coding rule of element identical; The evaluation object that increases again unit damage or centrepoint load unit change in calculating is each time different from the evaluation object that increases again unit damage or centrepoint load unit change in other calculating, calculate each time the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure, the current calculated value of all monitored amounts that calculate each time forms a monitored amount calculation current vector; When k evaluation object of hypothesis increases unit damage or centrepoint load unit change again, use C
^{i} _{tk}represent corresponding " monitored amount calculation current vector "; While giving in this step each vectorial element numbering, should use same coding rule with other vector in this method, to guarantee any one element in each vector in this step, with in other vector, number identical element, expressed the relevant information of same monitored amount or same target; C
^{i} _{tk}definition mode and vectorial C
_{o}definition mode identical, C
^{i} _{tk}element and C
_{o}element corresponding one by one;
H3. the vectorial C calculating each time
^{i} _{tk}deduct vectorial C
^{i} _{o}obtain a vector, then obtain " numerical value change vector δ a C for monitored amount after each element of this vector is calculated to the unit damage suppose or centrepoint load unit change numerical value divided by this
^{i} _{k}"; There is N evaluation object just to have N " the numerical value change vector of monitored amount ";
H4. by this N " the numerical value change vector of monitored amount ", according to the coding rule of N evaluation object, form successively " the unit damage monitored numerical quantity transformation matrices Δ C that has N row
^{i}"; Unit damage monitored numerical quantity transformation matrices Δ C
^{i}each row corresponding to a monitored amount unit change vector; Unit damage monitored numerical quantity transformation matrices Δ C
^{i}every a line corresponding to same monitored amount the different unit change amplitude when different evaluation objects increase unit damage or centrepoint load unit change; Unit damage monitored numerical quantity transformation matrices Δ C
^{i}coding rule and the vectorial d of row
_{o}the coding rule of element identical, unit damage monitored numerical quantity transformation matrices Δ C
^{i}the coding rule of coding rule and M monitored amount of row identical;
I. define the vectorial d of current name damage
^{i} _{c}with current actual damage vector d
^{i}, d
^{i} _{c}and d
^{i}element number equal the quantity of evaluation object, d
^{i} _{c}and d
^{i}element and evaluation object between be onetoone relationship, d
^{i} _{c}element numerical value represent nominal degree of injury or the nominal centrepoint load variable quantity of corresponding evaluation object, d
^{i} _{c}and d
^{i}with evaluation object initial damage vector d
_{o}element coding rule identical, d
^{i} _{c}element, d
^{i}element and d
_{o}element be onetoone relationship;
J. according to monitored amount current value vector C
^{i}with " the current initial value vector of monitored amount C
^{i} _{o}", " unit damage monitored numerical quantity transformation matrices Δ C
^{i}" and " the vectorial d of current name damage
^{i} _{c}" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, in formula 1 except d
^{i} _{c}other outer amount is known, solves formula 1 and just can calculate the vectorial d of current name damage
^{i} _{c};
K. the current actual damage vector d that utilizes formula 2 to express
^{i}k element d
^{i} _{k}with the current initial damage vector of evaluation object d
^{i} _{o}k element d
^{i} _{ok}with the vectorial d of current name damage
^{i} _{c}k element d
^{i} _{ck}between relation, calculate current actual damage vector d
^{i}all elements;
K=1 in formula 2,2,3 ..., N; Vector d
^{i}the coding rule of element and formula (1) in vectorial d
_{o}the coding rule of element identical; d
^{i} _{k}the current actual health status that represents k evaluation object in the i time circulation, if this evaluation object is support cable, so a d in cable system
^{i} _{k}the order of severity that represents its current health problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, d
^{i} _{k}numerical response the degree of lax or damage of this support cable; If this evaluation object is centrepoint load, so a d
^{i} _{k}the actual change amount that represents this centrepoint load; If this evaluation object is centrepoint load, so a d
^{i} _{k}the actual change amount that represents this centrepoint load; By the current actual damage vector of evaluation object d
^{i}in with M
_{1}the M that root support cable is relevant
_{1}individual element takes out, and forms the current actual damage vector of support cable d
^{ci}, the current actual damage vector of support cable d
^{ci}the coding rule and initial rope force vector F of element
_{o}the coding rule of element identical; The current actual damage vector of support cable d
^{ci}h element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., M
_{1}; The current actual damage vector of support cable d
^{ci}middle numerical value be not 0 element corresponding to the support cable of unsoundness problem, from the support cable of these unsoundness problems, identify damaged cable, remaining is exactly slack line; The current actual damage vector of the support cable corresponding with damaged cable d
^{ci}in element numerical expression be the current actual damage of this damaged cable, element numerical value represents while being 100% that this support cable thoroughly loses loadbearing capacity, represents to lose the loadbearing capacity of corresponding proportion in the time of between 0 and 100%; Utilization slack line under current cable structure steady temperature data vector Ti condition, that identify in the 1st step and with the current actual damage of support cable vector d
^{ci}these slack lines of expressing, with the current actual equivalent damage degree of its relax level mechanics equivalence, utilize in f step, obtain at current cable structure steady temperature data vector T
^{i}current cable force vector F under condition
^{i}with current support cable two supporting end points horizontal ranges vectors, utilize in c step, obtain at initial Cable Structure steady temperature data vector T
_{o}the initial drift vector of the support cable under condition, the weight vector of the initial free unit length of initial free crosssectional area vector sum, initial rope force vector F
_{o}, utilize current cable structure steady temperature data vector T
^{i}the current steady temperature data of support cable that represent, utilize in c step, obtain at initial Cable Structure steady temperature data vector T
_{o}the support cable initial steady state temperature data representing, the temperature variant physical and mechanical properties parameter of the various materials that the Cable Structure that utilization obtains in c step is used, count the impact of temperature variation on support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanics equivalence calculate slack line, with the relax level of current actual equivalent damage degree equivalence, mechanics equivalent condition is: one, the mechanics parameters without lax initial drift, geometrical property parameter, density and material during with not damaged of the rope of two equivalences is identical; Two, after lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after overall length identical; While meeting abovementioned two mechanics equivalent conditions, the mechanics function of two such support cables in Cable Structure is exactly identical, if replaced after damaged cable with equivalent slack line, Cable Structure any variation can not occur, and vice versa; According to aforementioned mechanics equivalent condition, try to achieve the relax level that those are judged as slack line, relax level is exactly the change amount of support cable drift, has namely determined the long adjustment amount of rope of the support cable that those need adjust Suo Li; Lax identification and the damage identification of support cable have so just been realized; During calculating, institute's demand power is by current cable force vector F
^{i}corresponding element provides; This method is referred to as damaged cable and slack line the support cable of unsoundness problem, referred to as problem rope, so far this method has realized and has rejected problem rope identification impact, Cable Structure that angular displacement of support, centrepoint load variation and structure temperature change, and has realized simultaneously and has rejected angular displacement of support, structure temperature variation and identification support cable health status variable effect, centrepoint load variable quantity;
1. try to achieve the vectorial d of current name damage
^{i} _{c}after, according to formula 3, set up mark vector B
^{i}, formula 4 has provided mark vector B
^{i}the definition of k element;
Element B in formula 4
^{i} _{k}mark vector B
^{i}k element, D
^{i} _{uk}evaluation object unit change vector D
^{i} _{u}k element, d
^{i} _{ck}the vectorial d of the current name damage of evaluation object
^{i} _{c}k element, they all represent the relevant information of k evaluation object, k=1 in formula 4,2,3 ..., N;
If mark vector B m.
^{i}element be 0 entirely, get back to step f and continue this circulation; If mark vector B
^{i}element be not 0 entirely, enter next step, be step n;
N. according to formula 5 calculate next time, i.e. the i+1 time current initial damage vector of the required evaluation object of circulation d
^{i+1} _{o}each element;
D in formula 5
^{i+1} _{ok}the current initial damage vector of the required evaluation object d that next time, circulates for the i+1 time
^{i+1} _{o}k element, d
^{i} _{ok}the current initial damage vector of the evaluation object d that is this, circulates for the i time
^{i} _{o}k element, D
^{i} _{uk}the evaluation object unit change vector D of the i time circulation
^{i} _{u}k element, B
^{i} _{k}the mark vector B of the i time circulation
^{i}k element, k=1 in formula 5,2,3 ..., N;
O. at initial mechanical calculating benchmark model A
_{o}basis on, first to A
_{o}in Cable Structure bearing apply angular displacement of support constraint, the numerical value of angular displacement of support constraint is just taken from the numerical value of corresponding element in angular displacement of support vector V, then to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation applying is just taken from steady temperature change vector S, then to make the health status of rope be d
^{i+1} _{o}after obtain be exactly next time, i.e. the i+1 time required Mechanics Calculation benchmark model A of circulation
^{i+1}; Obtain A
^{i+1}after, by Mechanics Calculation, obtain A
^{i+1}in concrete numerical value all monitored amounts, current, these concrete numerical value form next time, the vectorial C of the current initial value of required monitored amount that circulates for the i+1 time
^{i+1} _{o};
P. take off once, i.e. the i+1 time required current initial Cable Structure steady temperature data vector T of circulation
^{i+1} _{o}equal the current initial Cable Structure steady temperature data vector T of the i time circulation
^{i} _{o}; The required current initial Cable Structure bearing angular coordinate vector U of the i+1 time circulation next time, i.e.
^{i+1} _{o}equal the current initial Cable Structure bearing angular coordinate vector U of the i time circulation
^{i} _{o};
Q. get back to step f, start circulation next time.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201310662024.9A CN103616217A (en)  20131209  20131209  Progressive identifying method for problem cable and concentrated loads based on mixture monitoring in time of angular displacement 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201310662024.9A CN103616217A (en)  20131209  20131209  Progressive identifying method for problem cable and concentrated loads based on mixture monitoring in time of angular displacement 
Publications (1)
Publication Number  Publication Date 

CN103616217A true CN103616217A (en)  20140305 
Family
ID=50166923
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201310662024.9A Pending CN103616217A (en)  20131209  20131209  Progressive identifying method for problem cable and concentrated loads based on mixture monitoring in time of angular displacement 
Country Status (1)
Country  Link 

CN (1)  CN103616217A (en) 
Cited By (6)
Publication number  Priority date  Publication date  Assignee  Title 

CN103852317A (en) *  20140310  20140611  东南大学  Angular displacement mixedmonitoring load progressive identification method for problematic cable 
CN103852311A (en) *  20140310  20140611  东南大学  Linear displacement hybrid monitoring problem cable load progressive identification method 
CN103868713A (en) *  20140310  20140618  东南大学  Hybrid monitoring damaged cable load angular displacement identification method 
CN103868729A (en) *  20140310  20140618  东南大学  Progressive type hybrid monitoring damaged cable load angular displacement identification method 
CN105021419A (en) *  20150723  20151104  东南大学  Simple damaged cable load progressive identification method based on angular displacement space coordinatemonitoring 
CN105067344A (en) *  20150723  20151118  东南大学  Method for progressively recognizing load angular displacement of damaged cable based on streamlined hybrid monitoring process 
Citations (4)
Publication number  Priority date  Publication date  Assignee  Title 

EP0090730A2 (en) *  19820330  19831005  S.T.A.S. SOCIETE TECHNIQUE D'ACCESSOIRES SPECIALISES Société Anonyme  Fatigue monitor 
JPS58169002A (en) *  19820331  19831005  Mitsubishi Heavy Ind Ltd  Attaching method of welding type strain gage 
CN102297777A (en) *  20110513  20111228  东南大学  Progressive health monitoring method for cable system based on hybrid monitoring during supportingbase angular displacement 
CN102706658A (en) *  20120530  20121003  东南大学  Defective cable and support angular displacement progressive identification method based on mixed monitoring of temperature change 

2013
 20131209 CN CN201310662024.9A patent/CN103616217A/en active Pending
Patent Citations (5)
Publication number  Priority date  Publication date  Assignee  Title 

EP0090730A2 (en) *  19820330  19831005  S.T.A.S. SOCIETE TECHNIQUE D'ACCESSOIRES SPECIALISES Société Anonyme  Fatigue monitor 
EP0090730B1 (en) *  19820330  19860528  S.T.A.S. SOCIETE TECHNIQUE D'ACCESSOIRES SPECIALISES Société Anonyme  Fatigue monitor 
JPS58169002A (en) *  19820331  19831005  Mitsubishi Heavy Ind Ltd  Attaching method of welding type strain gage 
CN102297777A (en) *  20110513  20111228  东南大学  Progressive health monitoring method for cable system based on hybrid monitoring during supportingbase angular displacement 
CN102706658A (en) *  20120530  20121003  东南大学  Defective cable and support angular displacement progressive identification method based on mixed monitoring of temperature change 
NonPatent Citations (5)
Title 

孙景惠等: "缆索吊车悬索系统运动载荷的动力响应", 《林业机械》 * 
岳丽娜: "大跨悬索桥安全监测方法及体系研究与应用", 《中国博士学位论文全文数据库》 * 
徐藏等: "悬索式管桥应变影响因素研究", 《石油化工设备》 * 
汝继星等: "斜拉索管桥应变影响因素的试验分析", 《科学技术与工程》 * 
赵翔: "拉索损伤对斜拉桥结构性能影响的研究", 《中国博士学位论文全文数据库》 * 
Cited By (6)
Publication number  Priority date  Publication date  Assignee  Title 

CN103852317A (en) *  20140310  20140611  东南大学  Angular displacement mixedmonitoring load progressive identification method for problematic cable 
CN103852311A (en) *  20140310  20140611  东南大学  Linear displacement hybrid monitoring problem cable load progressive identification method 
CN103868713A (en) *  20140310  20140618  东南大学  Hybrid monitoring damaged cable load angular displacement identification method 
CN103868729A (en) *  20140310  20140618  东南大学  Progressive type hybrid monitoring damaged cable load angular displacement identification method 
CN105021419A (en) *  20150723  20151104  东南大学  Simple damaged cable load progressive identification method based on angular displacement space coordinatemonitoring 
CN105067344A (en) *  20150723  20151118  东南大学  Method for progressively recognizing load angular displacement of damaged cable based on streamlined hybrid monitoring process 
Similar Documents
Publication  Publication Date  Title 

CN103616222A (en)  Hybridmonitoring linear displacement progressive identification method for defective cable and concentrated loads  
CN103604549A (en)  Linear displacement progressive type identification method for strain monitoring of damaged cable and intensive load  
CN103604642A (en)  Problematic cable/concentrated load progressive recognition method on basis of generalized displacement cable power monitoring  
CN103604644A (en)  Problematic cable/concentrated load progressive recognition method on basis of linear displacement strain monitoring  
CN103604628A (en)  Problematic cable/concentrated load progressive recognition method on basis of linear displacement angle monitoring  
CN103604649A (en)  Problematic cable/concentrated load/linear displacement progressive recognition method on basis of strain monitoring  
CN103604623A (en)  Progressive type damaged rope and concentrated load identification method based on angular displacement strain monitoring  
CN103604634A (en)  Linear displacement progressive type identification method for space coordinate monitoring of damaged cable and intensive load  
CN103616217A (en)  Progressive identifying method for problem cable and concentrated loads based on mixture monitoring in time of angular displacement  
CN103604640A (en)  Method for progressive recognition of concentrated loads of damaged cables based on generalized displacement angle monitoring  
CN103616231A (en)  Generalized displacement strain monitoring progressive identification method for damaged cable and concentrated loads  
CN103616242A (en)  Method for recognizing concentrated loads of damaged cable progressively under condition of generalized displacement based on space coordinate monitoring  
CN103604627A (en)  Problematic cable/concentrated load/generalized displacement progressive recognition method on basis of hybrid monitoring  
CN103616200A (en)  Defective cable concentrated load progressive type recognition method based on spatial coordinate monitoring in process of linear displacement  
CN103616248A (en)  Progressive identifying method for problem cable concentrated loads and angular displacement based on angle monitoring  
CN103616243A (en)  Strainmonitoring damaged cable centralized load generalized displacement progressive recognition method  
CN103604632A (en)  Problematic cable/concentrated load progressive recognition method on basis of generalized displacement angle monitoring  
CN103604654A (en)  Problematic cable/concentrated load/generalized displacement progressive recognition method on basis of space coordinate monitoring  
CN103604554A (en)  Problematic cable/concentrated load/generalized displacement progressive recognition method on basis of cable power monitoring  
CN103604631A (en)  Problematic cable/concentrated load progressive recognition method on basis of linear displacement cable power monitoring  
CN103616241A (en)  Generalized displacement strain monitoring progressive identification method for defective cable and concentrated loads  
CN103616197A (en)  Hybridmonitoring progressive identification method for defective cable and concentrated loads  
CN103616114A (en)  Progressive identifying method for problem cable and concentrated loads based on cable force monitoring in time of angular displacement  
CN103604652A (en)  Problematic cable/concentrated load/angular displacement progressive recognition method on basis of strain monitoring  
CN103604639A (en)  Progressive type damaged rope and concentrated load identification method based on hybrid angular displacement monitoring 
Legal Events
Date  Code  Title  Description 

PB01  Publication  
PB01  Publication  
SE01  Entry into force of request for substantive examination  
C10  Entry into substantive examination  
RJ01  Rejection of invention patent application after publication 
Application publication date: 20140305 

RJ01  Rejection of invention patent application after publication 