CN102735460B  Progressive generalized displacement identification method for problem cables and support based on temperature change angle monitoring  Google Patents
Progressive generalized displacement identification method for problem cables and support based on temperature change angle monitoring Download PDFInfo
 Publication number
 CN102735460B CN102735460B CN201210171652.2A CN201210171652A CN102735460B CN 102735460 B CN102735460 B CN 102735460B CN 201210171652 A CN201210171652 A CN 201210171652A CN 102735460 B CN102735460 B CN 102735460B
 Authority
 CN
 China
 Prior art keywords
 cable
 temperature
 data
 vector
 initial
 Prior art date
Links
 238000006073 displacement reaction Methods 0.000 title claims abstract description 190
 230000000750 progressive Effects 0.000 title abstract 3
 238000004364 calculation methods Methods 0.000 claims abstract description 163
 239000011159 matrix materials Substances 0.000 claims abstract description 43
 239000002965 ropes Substances 0.000 claims description 77
 230000000875 corresponding Effects 0.000 claims description 71
 230000003862 health status Effects 0.000 claims description 50
 238000009826 distribution Methods 0.000 claims description 43
 239000000203 mixtures Substances 0.000 claims description 34
 239000000463 materials Substances 0.000 claims description 30
 238000009659 nondestructive testing Methods 0.000 claims description 18
 241001081830 Degeneriaceae Species 0.000 claims description 12
 238000000034 methods Methods 0.000 claims description 11
 206010022114 Injuries Diseases 0.000 claims description 10
 230000004048 modification Effects 0.000 claims description 9
 238000006011 modification reactions Methods 0.000 claims description 9
 230000035852 Tmax Effects 0.000 claims description 8
 230000000694 effects Effects 0.000 claims description 7
 229910052799 carbon Inorganic materials 0.000 claims description 6
 229910000975 Carbon steel Inorganic materials 0.000 claims description 5
 241001269238 Data Species 0.000 claims description 5
 238000009413 insulation Methods 0.000 claims description 5
 238000004458 analytical methods Methods 0.000 claims description 4
 239000010962 carbon steel Substances 0.000 claims description 4
 239000002131 composite materials Substances 0.000 claims description 4
 238000004519 manufacturing process Methods 0.000 claims description 4
 238000005259 measurements Methods 0.000 claims description 4
 239000010902 straw Substances 0.000 claims description 4
 230000001066 destructive Effects 0.000 claims description 2
 238000000547 structure data Methods 0.000 claims description 2
 238000005457 optimization Methods 0.000 description 9
 238000007796 conventional methods Methods 0.000 description 8
 230000001771 impaired Effects 0.000 description 7
 238000009529 body temperature measurement Methods 0.000 description 3
 238000004891 communication Methods 0.000 description 3
 239000000725 suspensions Substances 0.000 description 3
 281000099506 Members Only companies 0.000 description 2
 230000005540 biological transmission Effects 0.000 description 2
 230000015572 biosynthetic process Effects 0.000 description 2
 238000001514 detection method Methods 0.000 description 2
 238000005755 formation reactions Methods 0.000 description 2
 238000003860 storage Methods 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
 241001251094 Formica Species 0.000 description 1
 206010043431 Thinking abnormal Diseases 0.000 description 1
 230000001133 acceleration Effects 0.000 description 1
 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 1
 229910000019 calcium carbonate Inorganic materials 0.000 description 1
 239000008105 calcium carbonate Substances 0.000 description 1
 238000005260 corrosion Methods 0.000 description 1
 238000005516 engineering processes 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
 230000001681 protective Effects 0.000 description 1
 230000002040 relaxant effect Effects 0.000 description 1
 239000007787 solids 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 generalized displacement identification method for problem cables and a support based on temperature change angle monitoring. The progressive generalized displacement identification method can be used for deciding whether to update a mechanical calculation reference model of a cable structure through monitoring a cable structure temperature, an environment temperature and a heath state of an evaluated object, thus obtaining the mechanical calculation reference model of the cable structure including the cable structure temperature and the environment temperature; and a unit change matrix of a unit damage monitored quantity is obtained through calculation on the basis of the mechanical calculation reference model. A noninferior solution of a current nominal damage vector of the evaluated object is figured out according to an approximate linear relationship existing among a current value vector of a monitored quantity, a current initial value vector of the monitored quantity, the unit change matrix of the unit damage monitored quantity, a unit damage or unit generalized displacement vector and a current nominal damage vector of the tobesolved evaluated object, therefore, generalized displacements of the damaged cable, the slack cable and the support can be rapidly identified when temperature changes exist.
Description
Technical field
The structures such as cablestayed bridge, suspension bridge, trussframe structure have a common ground, be exactly that they have many parts bearing 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 only bears the rod member of tensile load for support unit, and such structure representation is " Cable Structure " by this method for simplicity.Along with the change of environment temperature, the temperature of Cable Structure also can change, and when Cable Structure temperature changes, there being generalized displacement of support, (such as generalized displacement of support refers to that bearing is along the angular displacement around X, Y, Z axis of the displacement of the lines of X, Y, Z axis and bearing; Corresponding to generalized displacement of support, bearing generalized coordinate refers to that bearing is about the coordinate of X, Y, Z axis and the bearing angular coordinate about X, Y, Z axis) time, this method identifies damaged cable, slack line (just referring to the impaired or lax rod member only bearing tensile load to trussframe structure) and generalized displacement of support in the supporting system of Cable Structure based on hybrid monitoring, belong to engineering structure health monitoring field.The supporting system of Cable Structure refers to all ropeway carryingropes and all rod members only bearing tensile load play supporting role, for simplicity, whole support units of this class formation are collectively referred to as " cable system " by this patent, but in fact cable system not only refers to support cable, also comprise the rod member only bearing tensile load, in this method, censure all ropeway carryingropes and all rod members only bearing tensile load play supporting role with " support cable " this noun.Damaged cable and slack line are called the support cable of unsoundness problem, referred to as problem cable by this method.
Background technology
Cable Structure (particularly largescale Cable Structure, such as largescale cablestayed bridge, suspension bridge) its support cable there will be damage, the structural health problem such as lax after long service, its bearing there will be the structural health problems such as generalized displacement, these structural health problems all will cause the change of structural internal force, harmful effect is caused to the safety of structure, the inefficacy of structure will be caused time serious, therefore identify accurately and timely and need slack line, damaged cable and generalized displacement of support to be very important.
The health status of support cable system changes and after bearing generation generalized displacement, the change of the measurable parameter of structure can be caused, such as can cause the change of Suo Li, distortion or the strain of Cable Structure can be affected, shape or the volume coordinate of Cable Structure can be affected, the change of the angle coordinate of any imaginary line of the every bit of the Cable Structure (change of the angle coordinate of the straight line of any this point of mistake in the section of such as body structure surface any point can be caused, or the change of the angle coordinate of the normal of body structure surface any point), these all changes all contain the health status information of cable system, in fact the change of these measurable parameters contains the health status information of cable system, contain generalized displacement of support information, that is the measurable parameter of structure can be utilized to identify generalized displacement of support and damaged cable, therefore the health status of structure can be judged by the hybrid monitoring of the change of the characteristic parameter to these dissimilar structures, all monitored structure characteristic parameters are referred to as " monitored amount " by this method, because now monitored amount is made up of the dissimilar measurable parameter mixing of structure, this method claims this to be hybrid monitoring, that is hybrid monitoring can be utilized to identify damaged cable and generalized displacement of support.Monitored amount is except being subject to the impact of Cable Structure health status; also can be subject to the impact of Cable Structure temperature variation (usually can occur); under the condition that Cable Structure temperature changes; if can based on the identification of the support cable realized the monitoring of monitored amount unsoundness problem and generalized displacement of support; to the safety of Cable Structure, there is important value, also do not have a kind of disclosed, effective health monitoring systems and method to solve this problem at present.
Summary of the invention
Technical matters: the object of this method is when Cable Structure has temperature variation, disclose a kind of based on hybrid monitoring, the health monitor method that can identify damaged cable in Cable Structure, slack line and generalized displacement of support rationally and effectively.
According to the reason that the Suo Li of support cable changes, the Suo Li of support cable change can be divided into four kinds of situations: one is that support cable receives damage, and localized cracks and corrosion etc. have appearred in such as support cable, two are support cable and not damaged, but Suo Li also there occurs change, occur that the one of the main reasons of this change is that Suo Changdu under support cable free state (now rope tensility also claims Suo Li to be 0) (is called drift, this method specially refers to that support cable two supports the drift of that section of rope between end points) there occurs change, one of object of this method will identify the support cable that drift there occurs change exactly, and identify the knots modification of their drift, this knots modification is that the cable force adjustment of this rope provides direct basis, the reason that support cable drift changes is not single, conveniently, the support cable that drift changes by this method is referred to as slack line, three is Suo Li changes that Cable Structure bearing generation generalized displacement causes, four is Suo Li changes that Cable Structure temperature variation causes.
Technical scheme: this method is made up of three parts.Respectively: one, " the temperature survey calculating method of the Cable Structure of this method "; Two, the Cable Structure health state evaluation method of the method for knowledge base needed for cable structure health monitoring system and parameter, knowledge based storehouse (containing parameter) and the monitored amount of actual measurement is set up; Three, the software and hardware part of health monitoring systems.
If the quantity sum of the generalized displacement of support component of the quantity of the support cable of Cable Structure and Cable Structure is N.For sake of convenience, this method unitedly calls evaluated support cable and generalized displacement of support to be " evaluation object ", total N number of evaluation object.To evaluation object serial number, this numbering will be used 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 M in cable system
_{1}root support cable, the monitored rope force data of structure specifies Q rope force data of support cable to describe by Q in structure, and the change of structure Suo Li is exactly all Q change of specifying the Suo Li of support cable.Each total Q cable force measurement value or calculated value carry out the rope force information of characterisation of structures.Q is one and is not less than 0, is not more than M
_{1}integer.
The monitored strain data of structure can by K in structure
_{2}the L of individual specified point and each specified point
_{2}the strain of individual assigned direction describes, and the change of structural strain data is exactly K
_{2}the change of all tested strain of individual specified point.Each total M
_{2}(M
_{2}=K
_{2}× L
_{2}) individual strain measurement value or calculated value carry out characterisation of structures strain.M
_{2}it is an integer being not less than 0.
The monitored angledata of structure is by K in structure
_{3}individual specified point, cross the L of each specified point
_{3}the H of individual appointment straight line, each appointment straight line
_{3}individual angle coordinate component describes, and the change of structural point is exactly the change of all specified points, all appointment straight line, all angle coordinate component of specifying.Each total M
_{3}(M
_{3}=K
_{3}× L
_{3}× H
_{3}) individual angle coordinate component measurement value or calculated value carry out the angle information of characterisation of structures.M
_{3}it is an integer being not less than 0.
The monitored shape data of structure is by K in structure
_{4}the L of individual specified point and each specified point
_{4}the volume coordinate of individual assigned direction describes, and the change of planform data is exactly K
_{4}the change of all coordinate components of individual specified point.Each total M
_{4}(M
_{4}=K
_{4}× L
_{4}) individual coordinates measurements or calculated value carry out characterisation of structures shape.M
_{4}it is an integer being not less than 0.
Comprehensive abovementioned monitored amount, total has M (M=Q+M
_{2}+ M
_{3}+ M
_{4}) individual monitored amount, definition parameter K (K=Q+K
_{2}+ K
_{3}+ K
_{4}), K and M must not be less than N.
Comprehensive abovementioned monitored amount, whole Cable Structure has M monitored amount, and M must not be less than the quantity N of evaluation object.
For simplicity, in the method by " monitored all parameters of Cable Structure " referred to as " monitored amount ".To M monitored amount serial number, this numbering will be used for generating vector sum matrix in subsequent step.This method represents this numbering, j=1,2,3 with variable j ..., M.
The Part I 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 ".Temperature due to Cable Structure may be change, the temperature of the different parts of such as Cable Structure changes along with the change of intensity of sunshine, along with the change of environment temperature changes, the surface of Cable Structure may be time dependent with inner temperature sometimes, the surface of Cable Structure may be different from inner temperature, the surface of Cable Structure is time dependent with inner temperature difference, this just makes the Mechanics Calculation of Cable Structure when considering temperature conditions and monitoring quite complicated, for simplifying 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 environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, asconstructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model of Cable Structure.Inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day
_{r}.Be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation
_{h}, get Δ T for convenience of describing
_{h}unit be DEG C/m.The surface of Cable Structure is got " R Cable Structure surface point ", the temperature of this R Cable Structure surface point will be obtained below by actual measurement, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ".When the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " describes later with the condition that must meet that distributes.From the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, at the sea level elevation place that each is chosen, two points are at least chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along 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 in the shade the outer normal direction of the measurement Cable Structure chosen along the sunny slope outer normal direction and Cable Structure that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness, three points are no less than along each measurement Cable Structure along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, especially, along each, Cable Structure is measured for support cable and only gets a point along the direction of the Temperature Distribution of wall thickness, namely the temperature of the surface point of support cable is only measured, measure all temperature be selected a little, the temperature recorded 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 the direction of Cable Structure along the Temperature Distribution of wall thickness " to 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 have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, especially, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", the number temperature profile data at sea level elevation place " identical sea level elevation Cable Structure is along the temperature profile data of thickness " will chosen at each in this method ".Measure temperature in Cable Structure location according to meteorology to require to choose a position, obtain meeting the temperature that meteorology measures the Cable Structure place environment of temperature requirement by the actual measurement of this position, in the onsite spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable, at the flat board of this position of sound production one piece of carbon steel material, 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 with dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, the nonsunny slope of reference plate is covered with insulation material, RealTime Monitoring is obtained the temperature of the sunny slope of reference plate.Must not be greater than 30 minutes to the time interval between any twice measurement of same amount RealTime Monitoring in this method, the moment of survey record data is called the physical record data moment.
Second step, RealTime Monitoring obtains R Cable Structure surface temperature measured data of abovementioned R Cable Structure surface point, RealTime Monitoring obtains the temperature profile data of previously defined Cable Structure along thickness simultaneously, and RealTime Monitoring obtains meeting the temperature record that meteorology measures the Cable Structure place environment of temperature requirement simultaneously, the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be designated as Δ T
_{emax}, calculated the rate of change of temperature about the time of Cable Structure place environment by Conventional mathematical by the temperature measured data sequence of Cable Structure place environment, this rate of change is also along with time variations, the measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be designated as Δ T
_{pmax}, the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein designated as Δ T
_{smax}, calculated the rate of change of temperature about the time of each Cable Structure surface point by Conventional mathematical by each Cable Structure surface temperature measured data sequence, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations.Obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by RealTime Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and minimum temperature at the sea level elevation place that each is chosen, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T
_{tmax}.
3rd step, survey calculation obtains Cable Structure steady temperature data, first, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day 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 Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the environment maximum error Δ T that survey calculation obtains above
_{emax}be not more than with reference to temperature difference per day Δ T
_{r}, and reference plate maximum temperature difference Δ T
_{pmax}Δ T is not more than after deducting 2 degrees Celsius
_{emax}, and Cable Structure surface maximum temperature difference Δ T
_{smax}be not more than Δ T
_{pmax}, one of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition, Section 3 condition is when obtaining Cable Structure steady temperature data, and 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, Section 4 condition is when obtaining Cable Structure steady temperature data, and 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, Section 5 condition is when obtaining 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 be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise, Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T
_{tmax}be not more than 1 degree Celsius, this method utilizes abovementioned six conditions, any one in following three kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in abovementioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, the second moment is moment of the Section 6 condition only met in abovementioned " condition relevant to determining to obtain moment of Cable Structure steady temperature data ", simultaneously the third moment meets the moment of the Section 1 in abovementioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition, when the mathematics moment obtaining Cable Structure steady temperature data is exactly one in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data, if the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data, the amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method, this method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, 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 obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution of the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steadystate surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steadystate surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steadystate 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 ", when the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%, Cable Structure surface comprises support cable surface, second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", 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 DEG C 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 obtained, gets Δ T for convenience of describing
_{h}unit be DEG C/m, be m for convenience of describing the unit getting Δ h, " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not 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, 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 geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshineduration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshineduration the most fully those surface points in Cable Structure.
The Part II of this method: the Cable Structure health state evaluation method setting up the method for knowledge base needed for cable structure health monitoring system and parameter, knowledge based storehouse (containing parameter) and the monitored amount of actual measurement.Can carry out successively as follows, to obtain the health state evaluation of evaluation object more accurately.
The first step: set up initial mechanical Calculation Basis model A
_{o}.
Set up initial mechanical Calculation Basis model A
_{o}when Cable Structure is completed, or before setting up cable structure health monitoring system, obtain " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " survey calculation (to measure by ordinary temperature measuring method, thermal resistance is such as used to measure), " Cable Structure steady temperature data " now use vector T
_{o}represent, be called initial Cable Structure steady temperature data vector T
_{o}.T is obtained in actual measurement
_{o}while, namely at the initial Cable Structure steady temperature data vector T of acquisition
_{o}the synchronization in moment, use the direct survey calculation of conventional method to obtain the initial number of all monitored amount of Cable Structure.Initial Cable Structure steady temperature data vector T is obtained at Actual measurement
_{o}while, use conventional method (consult reference materials or survey) to obtain temperature variant physical parameter (such as thermal expansivity) and the mechanical property parameters (such as elastic modulus, Poisson ratio) of the various materials that Cable Structure uses; Initial Cable Structure steady temperature data vector T is obtained at Actual measurement
_{o}while, namely at the initial Cable Structure steady temperature data vector T of acquisition
_{o}the synchronization in moment, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Actual measurement data of Cable Structure comprise the measured data such as data, the initial geometric data of Cable Structure, rope force data, drawbar pull data, initial Cable Structure bearing generalized coordinate data, the initial generalized displacement measurement data of Cable Structure bearing, Cable Structure modal data, structural strain data, structural point measurement data, structure space measurement of coordinates data that the Nondestructive Testing Data of support cable etc. can express the health status of rope.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 the geometric properties according to these coordinate data determination Cable Structure.Cable Structure bearing initial generalized displacement measurement data refers to setting up initial mechanical Calculation Basis model A
_{o}time, the generalized displacement that Cable Structure bearing occurs relative to the bearing under Cable Structure design point.For cablestayed bridge, initial geometric data can be the spatial data that the spatial data of the end points of all ropes adds some points on bridge two ends, socalled bridge type data that Here it is.Utilize the Nondestructive Testing Data etc. of support cable can express the data of the health status of support cable and the initial generalized displacement measurement data of Cable Structure bearing sets up evaluation object initial damage vector d
_{o}(such as formula (1) Suo Shi), uses d
_{o}represent that Cable Structure is (with initial mechanical Calculation Basis model A
_{o}represent) the initial health of evaluation object.If there is no the Nondestructive Testing Data of support cable and other are when can express the data of the health status of support cable, or can think structure original state be not damaged without relaxed state time, vectorial d
_{o}in each element numerical value relevant to support cable get 0, if when there is no Cable Structure bearing initial generalized displacement measurement data or can think that the initial generalized displacement of Cable Structure bearing is 0, vectorial d
_{o}in each element numerical value relevant to Cable Structure generalized displacement of support get 0.The temperature variant physical and mechanical properties parameter of the various materials utilizing the measured data of the design drawing of Cable Structure, asconstructed drawing and initial Cable Structure, the Nondestructive Testing Data of support cable, Cable Structure bearing initial generalized displacement measurement data, Cable Structure to use and initial Cable Structure steady temperature data vector T
_{o}, utilize mechanics method (such as finite element method) to count " Cable Structure steady temperature data " and set up initial mechanical Calculation Basis model A
_{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 Calculation Basis model A
_{o}in the original state of a kth evaluation object, if this evaluation object is a rope (or pull bar) in cable system, so d
_{ok}represent its initial damage, d
_{ok}represent not damaged when being 0, when being 100%, represent that this rope thoroughly loses loadbearing capacity, represent the loadbearing capacity losing corresponding proportion time between 0 and 100%, if this evaluation object is generalized displacement component, so a d of a bearing
_{ok}represent its initial displacement numerical value, T represents the transposition (same afterwards) of vector.
T is obtained in actual measurement
_{o}while, namely at the synchronization in the moment of acquisition Cable Structure steady temperature data, the initial value of all monitored amount of the Cable Structure using the direct survey calculation of conventional method to obtain, forms monitored amount initial value vector C
_{o}(see formula (2)).Require at acquisition A
_{o}while obtain C
_{o}, monitored amount initial value vector C
_{o}represent and correspond to A
_{o}the concrete numerical value of " monitored amount ".Because of subject to the foregoing, the Calculation Basis model based on Cable Structure calculates the monitored amount of gained reliably close to the measured data of initial monitored amount, in describing 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 jth monitored amount in Cable Structure, this component corresponds to a specific jth monitored amount according to coding rule.Vector C
_{o}be to be arranged according to a definite sequence by the monitored amount of M to form, putting in order to this there is no particular/special requirement, only requires all associated vector also array data in this order below.
No matter which kind of method to obtain initial mechanical Calculation Basis model A by
_{o}, count " Cable Structure steady temperature data " (i.e. initial Cable Structure steady temperature data vector T
_{o}), based on A
_{o}the Cable Structure that calculates calculates data must closely its measured data, and error generally must not be greater than 5%.Like this can utility A
_{o}the Suo Li calculated under the analog case of gained calculates data, strain calculation data, Cable Structure shapometer count certificate and displacement meter counts certificate, Cable Structure angledata, Cable Structure spatial data etc., measured data when reliably truly occurring close to institute's analog case.Model A
_{o}the health status evaluation object initial damage vector d of middle evaluation object
_{o}represent, the initial Cable Structure steady temperature data vector T of Cable Structure steady temperature data
_{o}represent.Due to based on A
_{o}the initial value (actual measurement obtains) of the evaluation calculating all monitored amounts closely all monitored amounts, 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}.T
_{o}and d
_{o}a
_{o}parameter, alternatively C
_{o}by A
_{o}mechanics Calculation result composition.
Second step: circulation starts.When circulation starts each time, first need the evaluation object current initial damage vector d set up or set up when this circulation starts
^{i} _{o}(i=1,2,3 ...), set up the current initial mechanical Calculation Basis model A of Cable Structure
^{i} _{o}(such as finite element benchmark model, A in circulation each time
^{i} _{o}constantly update), A
^{i} _{o}temperature Distribution " current initial Cable Structure steady temperature data vector T
^{i} _{o}" express.Letter i is except representing the place of number of steps significantly, and alphabetical i only represents cycle index in the method, i.e. ith circulation.A
_{o}and A
^{i} _{o}count temperature parameter, can the Effect on Mechanical Properties of accounting temperature change to Cable Structure.
The current initial damage vector of evaluation object that ith circulation needs when starting is designated as d
^{i} _{o}(such as formula (3) Suo Shi), uses d
^{i} _{o}when representing that this circulation starts, Cable Structure is (with current initial mechanical Calculation Basis model A
^{i} _{o}represent) the health status of evaluation object.
D in formula (3)
^{i} _{ok}(i=1,2,3, K=1,2,3 ...., N) represent ith time circulation start time, current initial mechanical Calculation Basis model A
^{i} _{o}in the original state of a kth evaluation object, if this evaluation object is a rope (or pull bar) in cable system, so d
^{i} _{ok}represent its initial damage, d
^{i} _{ok}represent not damaged when being 0, when being 100%, represent that this rope thoroughly loses loadbearing capacity, represent the loadbearing capacity losing corresponding proportion time between 0 and 100%, if this evaluation object is generalized displacement component, so a d of a bearing
^{i} _{ok}represent its initial displacement numerical value.
Set up and upgrade d
^{i} _{o}method as follows:
When first time, circulation started, (foundation formula (3) is designated as d to set up the current initial damage vector of evaluation object
^{1} _{o}) time, d
^{1} _{o}just equal d
_{o}.Ith (i=2,3,4,5,6 ...) the secondary evaluation object current initial damage vector d needed when starting that circulates
^{i} _{o}, be front once (namely the ith1 time, i=2,3,4,5,6 ...) circulation terminate before calculate obtain, concrete grammar is described below.
Ith (i=1,2,3,4,5,6 ...) secondary circulation needs the Mechanics Calculation benchmark model set up or the Mechanics Calculation benchmark model of Cable Structure set up to be designated as current initial mechanical Calculation Basis model A when starting
^{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, must set up or set up and be called current initial Cable Structure steady temperature data vector T when circulation starts each time
^{i} _{o}.
Set up, upgrade A
^{i} _{o}and T
^{i} _{o}method as follows:
The Mechanics Calculation benchmark model of the Cable Structure set up when first time, circulation started is designated as A
^{1} _{o}, A
^{1} _{o}equal A
_{o}, T
^{1} _{o}equal T
_{o}.A in circulation each time
^{i} _{o}and T
^{i} _{o}be constantly update, concrete grammar is described below; At the end of circulation each time, upgrade A
^{i} _{o}and T
^{i} _{o}the Mechanics Calculation benchmark model of the Cable Structure required when starting that next time circulated, concrete grammar is described below.
This method " monitored amount current initial value vector C
^{i} _{o}" (i=1,2,3 ...) initial value (see formula (4)) of all monitored amounts of specifying when representing that ith time (i=1,2,3,4,5,6 ...) circulation starts, C
^{i} _{o}also can be called " the monitored amount that circulates for ith time current initial value vector ".
C in formula (2)
^{i} _{oj}(i=1,2,3, J=1,2,3 ...., M) be jth monitored amount when circulating beginning for ith time, in Cable Structure.Vector C
^{i} _{o}be to be arranged according to a definite sequence by the monitored amount of previously defined M to form, putting in order to this there is no particular/special requirement, only requires all associated vector also array data in this order below.
At Modling model A
^{i} _{o}while set up " monitored amount current initial value vector C
^{i} _{o}", monitored amount current initial value vector C
^{i} _{o}represent and correspond to A
^{i} _{o}the concrete numerical value of all monitored amount, C
^{i} _{o}element and C
_{o}element one_to_one corresponding, represent that all monitored amounts are in A in Cable Structure respectively
^{i} _{o}and A
_{o}concrete numerical value during two states.
Set up and upgrade C
^{i} _{o}concrete grammar as follows:
When first time, circulation started, C
^{1} _{o}(i=1, C
^{i} _{o}be embodied as C
^{1} _{o}) equal C
_{o}; Ith (i=2,3,4,5,6 ...) secondary ith circulation " the monitored amount current initial value vector C needed when starting that circulates
^{i} _{o}", be front once (namely the ith1 time, i=2,3,4,5,6 ...) circulation terminate before calculate obtain, concrete grammar is described below.Ith (i=1,2,3,4,5,6 ...) in circulation, " monitored amount current initial value vector C
^{i} _{o}" be constantly update, concrete grammar is described below.Due to according to model A
^{i} _{o}the initial value calculating gained monitored amount, reliably close to corresponding measured value, in describing below, will represent this calculated value composition of vector and measured value composition of vector with prosign.
T
^{i} _{o}and d
^{i} _{o}a
^{i} _{o}characterisitic parameter, C
^{i} _{o}a
^{i} _{o}at T
^{i} _{o}and d
^{i} _{o}mechanics Calculation result composition under condition.
3rd step: in Cable Structure military service process, in circulation each time, in other words in the ith (i=1,2,3,4,5,6 ...) in secondary circulation, at known A
^{i} _{o}, T
^{i} _{o}, C
^{i} _{o}and d
^{i} _{o}after, the current data of " Cable Structure steady temperature data " is obtained, current data composition " the current cable structure steady temperature data vector T of all " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement
^{i}", vector T
^{i}definition mode and vector T
_{o}definition mode identical; In actual measurement vector T
^{i}while, namely at acquisition current cable structure steady temperature data vector T
^{i}the synchronization in moment, actual measurement obtains the currency of all monitored amounts in Cable Structure, and all these numerical value form monitored amount current value vector C
^{i}.C
^{i}element and C
_{o}element one_to_one corresponding, represent that identical monitored amount is at not numerical value in the same time.
In acquisition vector T
^{i}after, upgrade A according to following concrete grammar
^{i} _{o}, T
^{i} _{o}, C
^{i} _{o}and d
^{i} _{o}:
Relatively T
^{i}and T
^{i} _{o}if, T
^{i}equal T
^{i} _{o}, then do not need A
^{i} _{o}upgrade, otherwise need A
^{i} _{o}and T
^{i} _{o}upgrade, update method is: the first step calculates T
^{i}with T
_{o}difference, T
^{i}with T
_{o}difference be exactly the changes of current cable structure steady temperature data about initial Cable Structure steady temperature data, T
^{i}with T
_{o}difference represent with steady temperature change vector S, S equals T
^{i}deduct T
_{o}, S represents the change of Cable Structure steady temperature data; Second step is to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A
_{o}in Cable Structure apply temperature variation after obtain upgrade current initial mechanical Calculation Basis model A
^{i} _{o}, upgrade A
^{i} _{o}while, T
^{i} _{o}all elements numerical value also uses T
^{i}all elements numerical value correspondence replace, namely have updated T
^{i} _{o}, so just obtain and correctly correspond to A
^{i} _{o}t
^{i} _{o}; Now d
^{i} _{o}remain unchanged.As renewal A
^{i} _{o}after, A
^{i} _{o}the health status evaluation object of evaluation object current initial damage vector d
^{i} _{o}represent, A
^{i} _{o}cable Structure steady temperature current cable structure steady temperature data vector T
^{i}represent, obtain A by Mechanics Calculation
^{i} _{o}in all monitored amounts, current concrete numerical value, replace C with these concrete numerical value
^{i} _{o}the element of middle correspondence, so just achieves monitored amount current initial value vector C
^{i} _{o}renewal.
4th step: circulation time must first be set up " unit damage monitored numerical quantity unit change matrix " and " unit damage or unit generalized displacement vector " each time, and " unit damage monitored numerical quantity unit change matrix " that ith circulation is set up is designated as Δ C
^{i}(i=1,2,3 ...)." unit damage or unit generalized displacement vector " that ith circulation is set up is designated as D
^{i} _{u}.Δ C in circulation each time
^{i}and D
^{i} _{u}need according to circumstances to constantly update, namely at the current initial mechanical Calculation Basis model A of renewal
^{i} _{o}, current initial Cable Structure steady temperature data vector T
^{i} _{o}with monitored amount current initial value vector C
^{i} _{o}after, upgrade unit damage monitored numerical quantity unit change matrix Δ C
^{i}with unit damage or unit generalized displacement vector D
^{i} _{u}.
First unit damage monitored numerical quantity unit change matrix Δ C is set up in the steps below when circulation starts each time
^{i}with unit damage or unit generalized displacement vector D
^{i} _{u}; If have updated A in the third step
^{i} _{o}, (namely upgrading) unit damage monitored numerical quantity unit change matrix Δ C so must be reestablished in this step
^{i}with unit damage or unit generalized displacement vector D
^{i} _{u}; If do not upgrade A in the third step
^{i} _{o}, unit damage monitored numerical quantity unit change matrix Δ C so need not be reestablished in this step
^{i}with unit damage or unit generalized displacement vector D
^{i} _{u}; Set up and reestablish (namely upgrading) Δ C
^{i}and D
^{i} _{u}detailed process identical, arrange as follows:
At the current initial mechanical Calculation Basis model A of Cable Structure
^{i} _{o}basis on carry out several times calculating, calculation times numerically equals the quantity of all ropes.Calculate hypothesis each time to only have an evaluation object (original damage or generalized displacement of support can be 0 at original damage or generalized displacement of support, also can not be 0) basis on increase unit damage or unit generalized displacement again, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable increases unit damage again and (such as gets 5%, 10%, 20% or 30% equivalent damage is unit damage), if this evaluation object is the generalized displacement component in a direction of a bearing, just suppose that this bearing increases unit generalized displacement again in this generalized displacement direction and (such as gets 1 millimeter, 2 millimeters, 3 millimeters of grades are unit displacement of the lines, get 100,000/radian, 2/100000ths radians, 3/100000ths radians etc. are unit angular displacement).For convenience of calculating, can be all structural health conditions when this circulation is started when in circulation, setting increases unit damage or unit generalized displacement each time as being completely healthy, and set on this basis unit damage or unit generalized displacement (in subsequent step, calculate, the damage numerical value of evaluation object or generalized displacementbe called nominal fatigue d
^{i} _{c}(i=1,2,3 ...), all relative to when this circulation is started, by the health status of evaluation object as being completely healthy speech, therefore must the formula that hereinafter provides of foundation the nominal fatigue calculated is converted into true damage).With occurring in the calculating each time once circulated that the evaluation object of unit damage or unit generalized displacement is different from during other time calculates the evaluation object occurring unit damage or unit generalized displacement, and suppose there is unit damage value or unit generalized displacement value that the unit damage value of the evaluation object of unit damage or unit generalized displacement or unit generalized displacement value can be different from other evaluation objects each time, with " unit damage or unit generalized displacement vector D
^{i} _{u}" (such as formula (5) Suo Shi) record the unit damage of supposition or the unit generalized displacement of all evaluation objects in each circulation, first time circulation time be designated as D
^{1} _{u}calculate the current calculated value all utilizing M that mechanics method (such as finite element method) calculates Cable Structure, that specified above monitored amount each time, the current calculated value calculating a gained M monitored amount each time forms one " monitored amount calculation current vector ", and (when supposing that a kth evaluation object has unit damage or unit generalized displacement, available formula (6) represents the monitored amount calculation current vector C of all M specified a monitored amount
^{i} _{tk}); The monitored amount calculation current vector calculated each time deducts monitored amount current initial value vector C
^{i} _{o}obtain a vector, obtain after again each element of this vector being calculated divided by this unit damage or unit generalized displacement numerical value supposed one " the unit change vector of monitored amount ", " the unit change vector of monitored amount " that be exactly (to have the position of the rope of unit damage or numbering etc. for mark) under this condition (when a kth evaluation object has unit damage or unit generalized displacement, uses δ C
^{i} _{k}represent the unit change vector of 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}a kth element D
^{i} _{uk}gained), the unit change vector δ C of monitored amount
^{i} _{k}each element representation owing to suppose there is unit damage or unit generalized displacement (the such as D of that evaluation object (such as a kth evaluation object) of unit damage or unit generalized displacement when calculating
^{i} _{uk}), and the numerical value knots modification of monitored amount corresponding to this element caused is relative to the unit damage of supposition or unit generalized displacement numerical value D
^{i} _{uk}rate of change; N number of evaluation object is had just to have N number of " the unit change vector of monitored amount ", the unit change vector of each monitored amount has M element, forms by this N number of " unit change vector of monitored amount " " the unit damage monitored numerical quantity unit change matrix Δ C having M × N number of element successively
^{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 such as formula shown in (10).
Unit damage or unit generalized displacement vector D in formula (5)
^{i} _{u}element D
^{i} _{uk}(i=1,2,3, K=1,2,3 ...., N) represent the unit damage of a kth evaluation object or the unit generalized displacement numerical value of supposition in ith 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 ith circulation due to a kth evaluation object have unit damage or a unit generalized displacement time, according to the calculating current value of the individual monitored amount of specifying of the jth corresponding to coding rule.
In formula (7) each amount subscript i (i=1,2,3 ...) and represent ith circulation, subscript k (k=1,2,3 ...., N) represent unit damage or the unit generalized displacement of a kth evaluation object increase, D in formula
^{i} _{uk}vectorial D
^{i} _{u}in a kth element.Vector δ C
^{i} _{k}definition such as formula shown in (7) and formula (8), δ C
^{i} _{k}jth (j=1,2,3 ...., M) individual element δ C
^{i} _{kj}(defining such as formula (9) Suo Shi) represents in ith circulation, sets up matrix Δ C
^{i}time, assuming that a kth evaluation object has the knots modification calculating a gained jth monitored amount when unit damage or unit generalized displacement relative to the unit damage supposed or unit generalized displacement D
^{i} _{uk}rate of change.
Vectorial δ C in formula (10)
^{i} _{k}(i=1,2,3 ...., k=1,2,3 ...., N) represent in ith circulation, because a kth evaluation object increases unit damage or unit generalized displacement D
^{i} _{uk}cause, the change of the relative value of all monitored amounts.Matrix Δ C
^{i}the coding rule of row (subscript k) and vectorial d above
^{i} _{o}the coding rule of subscript k of element identical.
5th step: the current health state identifying Cable Structure.Detailed process is as follows.
Ith (i=1,2,3 ...) in secondary circulation, utilize " the monitored amount current value vector C obtained in second step actual measurement
^{i}" " monitored amount current initial value vector C together
^{i} _{o}", " unit damage monitored numerical quantity unit change matrix Δ C
^{i}" and " current nominal fatigue vector d
^{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 monitored amount current initial value vector C
^{i} _{o}definition, see formula (13); Evaluation object current nominal fatigue vector d
^{i} _{c}definition see formula (14).
Elements C in formula (13)
^{i} _{j}(i=1,2,3 ....; J=1,2,3 ...., M) be ith circulation time Cable Structure, the current value that be numbered the monitored amount of j of foundation corresponding to coding rule.
D in formula (14)
^{i} _{ck}(i=1,2,3 ....; K=1,2,3 ...., N) be the current nominal fatigue of a kth evaluation object in ith circulation or current nominal generalized displacement value, vectorial d
^{i} _{c}the coding rule of subscript k of element and matrix Δ C
^{i}the coding rule of row identical.
When support cable actual damage is not too large, generalized displacement of support is less, because Cable Structure material is still in the linear elasticity stage, the distortion of Cable Structure is also less, formula (11) or a kind of like this linear relationship represented by formula (12) less with the error of actual conditions, error can use error vector e
^{i}(formula (15)) define, the error of expression (11) or the shown linear relationship of formula (12).
In formula (15), abs () is the function that takes absolute value, and takes absolute value to each element of the vector of trying to achieve in bracket.
There is certain error in the linear relationship represented by formula (11) or formula (12), therefore can not simply according to formula (11) or formula (12) and " monitored amount current value vector C
^{i}" carry out direct solution and obtain evaluation object current nominal fatigue vector d
^{i} _{c}.If this has been doned, the evaluation object obtained current nominal fatigue vector d
^{i} _{c}in element even there will be larger negative value, namely negative damage, this is obviously irrational.Therefore evaluation object current nominal fatigue vector d is obtained
^{i} _{c}acceptable solution (namely with reasonable error, but position and the degree of injury thereof of damaged cable can be determined more exactly from cable system, also can determine generalized displacement of support numerical value more exactly) become a rational solution, available formula (16) expresses 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)), is defined by formula (17).
G in formula (17)
^{i} _{j}(i=1,2,3 ....; J=1,2,3 ...., M) describe the maximum allowable offset departing from formula (11) or the ideal linearity relation shown in formula (12) in ith circulation.Vector g
^{i}the error vector e that can define according to formula (15)
^{i}tentative calculation is selected.
At monitored amount current initial value vector C
^{i} _{o}, unit damage monitored numerical quantity unit change matrix Δ C
^{i}with monitored amount current value vector C
^{i}time known, suitable algorithm (such as multiobjective optimization algorithm) can be utilized to solve formula (16), obtain evaluation object current nominal fatigue vector d
^{i} _{c}acceptable solution, evaluation object current actual damage vector d
^{i}the element of (formula (18) is shown in definition) can calculate according to formula (19), thus can 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) represent the current actual health status of a kth evaluation object in ith circulation, formula (19) is shown in its definition, if this evaluation object is a support cable (or pull bar) in cable system, so d
^{i} _{k}represent its current actual damage, d
^{i} _{k}represent not damaged when being 0, when being 100%, represent that this support cable thoroughly loses loadbearing capacity, represent the loadbearing capacity losing corresponding proportion time between 0 and 100%, if this evaluation object is generalized displacement component, so a d of a bearing
^{i} _{k}represent its current actual generalized displacement numerical value, 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 evaluation object current initial damage vector d
^{i} _{o}a kth element, d
^{i} _{ck}evaluation object current nominal fatigue vector d
^{i} _{c}a kth element.
Describe below and obtain evaluation object current actual damage vector d
^{i}after, how to determine position and the relax level of slack line.
By M total 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 cable force vector F
_{o}" represent the Initial cable force (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 cable force of h root support cable in Cable Structure, this element corresponds to the Suo Li specifying support cable according to coding rule.Vector F
_{o}it is constant.T is obtained in actual measurement
_{o}synchronization, use conventional method direct survey calculation to obtain the rope force data of all support cables, all these rope force datas composition Initial cable force vector F
_{o}.Setting up the initial mechanical Calculation Basis model A of Cable Structure
_{o}time in fact employ vectorial F
_{o}.
By current for evaluation object actual damage vector d
^{i}in the M relevant to support cable
_{1}individual element takes out, composition support cable current actual damage vector d
^{ic}, support cable current actual damage vector d
^{ic}the coding rule of element and Initial cable force vector F
_{o}the coding rule of element identical.Support cable current actual damage vector d
^{ic}h element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., M
_{1}; Current actual damage vector d
^{ic}middle numerical value be not 0 element correspond to the support cable of unsoundness problem, carry out NonDestructive Testing to the support cable of these unsoundness problems, after NonDestructive Testing finds out that this support cable is not damaged, so this element numerical value (uses d
^{ic} _{h}represent) represent this support cable and d
^{ic} _{h}relaxing of impairment value mechanic equivalent, just determine slack line thus, the computing method of concrete slack are described below.
Represent the current cable power (formula (21) is shown in definition) of all support cables in the Cable Structure of surveying and obtaining with " current cable force vector F " in this method.
F in formula (21)
_{h}(h=1,2,3 ...., M
_{1}) be the current cable power of h root support cable in Cable Structure.Current cable structure steady temperature data vector T is obtained in actual measurement
^{i}synchronization, actual measurement obtains the rope force data of all support cables in Cable Structure, all these rope force datas composition current cable force vector F.The element of vector F and vectorial F
_{o}the coding rule of element identical.According to describing 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 of steady temperature data of Cable Structure
_{o}during expression, and when support cable is in free state (free state refers to that Suo Li is 0, rear same), 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.Vector T is passed through according to " the temperature survey calculating method of the Cable Structure of this method "
_{o}can determine obtaining vector T
_{o}the Temperature Distribution of all support cables in moment.
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}being constant, after determining when starting, just no longer changing.
Similar, under support cable original state, at the initial Cable Structure steady temperature data vector T of steady temperature data of Cable Structure
_{o}during expression, and when support cable is in free state, 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}being constant, after determining when starting, just no longer changing.
ω in formula (24)
_{oh}(h=1,2,3 ...., M
_{1}) be the weight of the free unit length of initial freedom of h root support cable in Cable Structure.Vector ω
_{o}being constant, after determining when starting, just no longer changing.
In this method, at the current initial Cable Structure steady temperature data vector T of steady temperature data of Cable Structure
^{t} _{o}during expression, with " current initial drift vector l
^{t} _{o}" (formula (25) is shown in definition; when referring to that hypothesis supporting cable force is 0; consider thermal expansivity and temperature variation to after the impact of support cable drift, initial drift vector l to represent the current initial drift of all support cables in Cable Structure
_{o}with initial Cable Structure steady temperature data vector T
_{o}the support cable represented is at the current initial Cable Structure steady temperature data vector T of 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 obtaining vector T
^{i} _{o}the Temperature Distribution of all support cables in moment.
L in formula (25)
^{t} _{oh}(h=1,2,3 ...., M
_{1}) be the current initial Cable Structure steady temperature data vector T of steady temperature data in Cable Structure
^{i} _{o}during expression, the current initial drift of h root support cable in Cable Structure, can utilize the thermal expansivity of support cable, l
_{oh}, T
_{o}and T
^{i} _{o}l is calculated by Typical physical
^{t} _{oh}.
Element, the vectorial l of vector F
_{o}element, vectorial l
^{t} _{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 of steady temperature data of Cable Structure
^{i} _{o}during expression, represent the current drift (formula (26) is shown in definition, and now support cable may be intact, may be impaired, also may be lax) of all support cables in Cable Structure with " current drift vector l ".
L in formula (26)
_{h}(h=1,2,3 ...., M
_{1}) be the current drift of h root support cable in Cable Structure.
In this method, represent the knots modification (formula (27) and formula (28) are shown in definition) of the drift of all support cables in Cable Structure with " drift changes vectorial Δ l " (or claiming support cable current slack degree vector).
Δl＝[Δl
_{1}Δl
_{2}···Δl
_{j}···Δl
_{N}]
^{T}(27)
Δ l in formula (27)
_{h}(h=1,2,3 ...., M
_{1}) be the knots modification of the drift of h root support cable in current cable structure, its definition is shown in formula (28), Δ l
_{h}be not 0 rope be slack line, Δ l
_{h}numerical value be the slack of rope, and representing the current slack degree of cable system h root support cable, is also the long adjustment amount of rope of this rope during adjustment Suo Li.
In the method by slack line is carried out with damaged cable the relax level identification that mechanic equivalent carries out slack line, the mechanical condition of equivalence is:
One, the nothing of the rope of two equivalences is lax identical with the mechanics parameters of initial drift during not damaged, geometrical property parameter and material;
Two, after lax or damage, the slack line of two equivalences and the Suo Li of damage rope be out of shape after overall length identical.
When meeting abovementioned two equivalent conditions, such two support cables mechanics function is in the structure exactly identical, if after namely replacing slack line with the damaged cable of equivalence, any change can not occur Cable Structure, and vice versa.
Obtain support cable current actual damage vector d
^{ic}after, d
^{ic}h element d
^{ic} _{h}(h=1,2,3 ...., M
_{1}) representing the actual damage value of h root support cable, its definition is shown in formula (19), although by d
^{ic} _{h}be called the actual damage value of h root support cable or the actual damage degree of h root support cable, but also may be lax because h root Suo Keneng is impaired, so d
^{ic}h element d
^{ic} _{h}the actual damage value of the h root support cable represented is actually the actual equivalent damage value of h root support cable, when h root support cable is actually impaired, and d
^{ic} _{h}the actual damage value of the h root support cable just represented, when h root support cable is actually lax, d
^{ic} _{h}the h root support cable just represented with the actual damage value of lax equivalence, for sake of convenience, claim d in the method
^{ic} _{h}represent h root support cable not damaged when being 0, when being 100%, represent that this rope thoroughly loses loadbearing capacity, time between 0 and 100%, represent that h root support cable loses the loadbearing capacity of corresponding proportion, by evaluation object current actual damage vector d
^{ic}just can identify the support cable that health status goes wrong, but in the support cable that goes wrong of these health status, some is impaired, some relaxes, if h root support cable is actually relax (its current slack degree Δ l
_{j}definition), the current slack degree Δ l of so lax h root support cable
_{j}(Δ l
_{j}definition see formula (27)) with the current actual damage degree d of damaged cable of equivalence
^{ic} _{h}between relation determined by aforementioned two mechanic equivalent conditions.Δ l
_{h}same d
^{ic} _{h}between physical relationship can adopt accomplished in many ways, such as can directly determine (see formula (29)) according to aforementioned equivalent condition, also can adopt after replacing the E in formula (25) to revise based on Ernst equivalent elastic modulus and determine (see formula (30)), other method such as trial and error procedure based on finite element method also can be adopted to determine.
Formula (29) and the middle E of formula (30)
^{t} _{h}the current initial Cable Structure steady temperature data vector T of steady temperature data in Cable Structure
^{i} _{o}during expression, the elastic modulus of h root support cable, A
^{t} _{h}the current initial Cable Structure steady temperature data vector T of steady temperature data in Cable Structure
^{i} _{o}during expression, the crosssectional area of h root support cable, F
_{h}the current initial Cable Structure steady temperature data vector T of steady temperature data in Cable Structure
^{i} _{o}during expression, the current cable power of h root support cable, d
^{i} _{h}the current actual damage degree of h root support cable, ω
^{t} _{h}the current initial Cable Structure steady temperature data vector T of steady temperature data in Cable Structure
^{t} _{o}during expression, the weight of the unit length of h root support cable, l
^{t} _{xh}the current initial Cable Structure steady temperature data vector T of steady temperature data in Cable Structure
^{i} _{o}during expression, the horizontal range of two supporting end points of h root support cable, l
^{t} _{xh}that current support cable two supports end points horizontal range vector l
^{t} _{x}an element, current support cable two support end points horizontal range vector l
^{t} _{x}the coding rule of element and initial drift vector l
_{o}the coding rule of element identical, E
^{t} _{h}can obtain according to the characteristic material data looking into or survey h root support cable, A
^{t} _{h}and ω
^{t} _{h}can according to the thermal expansivity of h root support cable, A
_{oh}, ω
_{oh}, F
_{h}, T
_{o}and T
^{i} _{o}obtained by Typical physical and Mechanics Calculation.Item in formula (30) in [] is the Ernst equivalent elastic modulus of this support cable, just can determine support cable current slack degree vector Δ l by formula (29) or formula (30).Formula (30) is the correction to formula (29).
6th step: judge whether to terminate this (ith time) circulation, if so, then completes this tailing in work before terminating that circulates, for next time (namely the ith+1 time, i=1,2,3,4 ...) circulating prepares Mechanics Calculation benchmark model and necessary vector.Detailed process is as follows:
Current nominal fatigue vector d is tried to achieve in this (ith time) circulation
^{i} _{c}after, first, set up mark vector B according to formula (31)
^{i}, formula (32) gives mark vector B
^{i}the definition of a kth element; If mark vector B
^{i}element be 0 entirely, then get back to the 3rd step and proceed health monitoring to Cable Structure and calculating; If mark vector B
^{i}element be not 0, then after completing subsequent step entirely, enter and circulate next time.
Socalled subsequent step is: first, according to formula (33) calculate next time (namely the ith+1 time, i=1,2,3,4 ...) needed for circulation initial damage vector d
^{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}instead of 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 of applying just takes from steady temperature change vector S, and steady temperature change vector S equals T
^{i}deduct T
_{o}), so just obtain next time (namely the ith+1 time, i=1,2,3,4 ...) current initial mechanical Calculation Basis mould A needed for circulation
^{i+1} _{o}, next time (namely the ith+1 time, i=1,2,3,4 ...) current initial Cable Structure steady temperature data vector T needed for circulation
^{i+1} _{o}equal T
^{i} _{o}, to A
^{i+1} _{o}carry out Mechanics Calculation to obtain corresponding to A
^{i+1} _{o}all monitored amount, current concrete numerical value, these concrete numerical value compositions next time (namely the ith+1 time, i=1,2,3,4 ...) the current initial value vector C of monitored amount needed for circulation
^{i+1} _{o}.
Mark vector B in formula (31)
^{i}subscript i represent ith circulation, its element B
^{i} _{k}(k=1,2,3 ..., N) subscript k represent the health status feature of a kth evaluation object, can only get 0 and 1 two amount, concrete value rule is shown in formula (32).
Element B in formula (32)
^{i} _{k}mark vector B
^{i}a kth element, D
^{i} _{uk}unit damage or unit generalized displacement vector D
^{i} _{u}a kth element (see formula (3)), d
^{i} _{ck}evaluation object current nominal fatigue vector d
^{i} _{c}a kth element (see formula (14)), they all represent the relevant information of a kth evaluation object.
D in formula (33)
^{i} _{uk}unit damage or unit generalized displacement vector D
^{i} _{u}a kth element (see formula (5)), d
^{i} _{ok}evaluation object current initial damage vector d
^{i} _{o}a kth element (see formula (3)).
The Part III of this method: the software and hardware part of health monitoring systems.
Hardware components comprises monitoring system (comprising monitored amount monitoring system, temperature monitoring system), signal picker and computing machine etc.Require that RealTime Monitoring obtains temperature required measured data, require each monitored amount of RealTime Monitoring simultaneously.
Software section should complete the process set by this method, namely to complete in this method required, can by functions such as computer implemented monitoring, record, control, storage, calculating, notice, warnings.
This method specifically comprises:
A. for sake of convenience, this method unitedly calls evaluated support cable and generalized displacement of support component to be evaluation object, if the quantity sum of the quantity of evaluated support cable and generalized displacement of support component is N, namely the quantity of evaluation object is N; Determine the coding rule of evaluation object, evaluation objects all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step; This method variable k represents this numbering, k=1,2,3 ..., N; Specify when determining hybrid monitoring by the support cable of monitored Suo Li, if total M in cable system
_{1}root support cable, the monitored rope force data of Cable Structure specifies Q rope force data of support cable to describe by Q in Cable Structure, and the change of Cable Structure Suo Li is exactly the change of the Suo Li of all appointment support cables; Each total Q cable force measurement value or calculated value characterize the rope force information of Cable Structure; Q is one and is not less than 0, is not more than M
_{1}integer; Specify when determining hybrid monitoring by the measured point of monitored strain, the monitored strain data of Cable Structure can by K in Cable Structure
_{2}the L of individual specified point and each specified point
_{2}the strain of individual assigned direction describes, and the change of Cable Structure strain data is exactly K
_{2}the change of all tested strain 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}be be not less than 0 integer; Specify when determining hybrid monitoring by the measured point of monitored angle, the monitored angledata of Cable Structure is by K in Cable Structure
_{3}individual specified point, cross the L of each specified point
_{3}the H of individual appointment straight line, each appointment straight line
_{3}individual angle coordinate component describes, and the change of Cable Structure angle is exactly the change of all specified points, all appointment straight line, all angle coordinate component of specifying; 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}it is an integer being not less than 0; Specify when determining hybrid monitoring by monitored shape data, the monitored shape data of Cable Structure is by K in Cable Structure
_{4}the L of individual specified point and each specified point
_{4}the volume coordinate of individual assigned direction describes, and the change of Cable Structure shape data is exactly K
_{4}the change of all coordinate components of individual specified point; Each total M
_{4}individual coordinates measurements or calculated value characterize Cable Structure shape, M
_{4}for K
_{4}and L
_{4}longpending; M
_{4}it is an integer being not less than 0; The monitored amount of comprehensive abovementioned hybrid monitoring, whole Cable Structure has M monitored amount, and M is Q, M
_{2}, M
_{3}and M
_{4}sum, definition parameter K, K is Q, K
_{2}, K
_{3}and K
_{4}sum, K and M must not be less than the quantity N of evaluation object; Because M monitored amount is dissimilar, so this method is called " during temperature variation the damaged cable of hybrid monitoring and support translation identification method "; For simplicity, in the method by monitored all parameters of Cable Structure " during the hybrid monitoring " listed by this step referred to as " monitored amount "; The quantity sum of all monitored amounts is designated as M, and M must not be less than N; Must not be greater than 30 minutes to the time interval between any twice measurement of same amount RealTime Monitoring in this method, the moment of survey record data is called the physical record data moment;
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 environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, asconstructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model of Cable Structure, inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, in the method daytime can not be seen one of the sun and be called the cloudy day all day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, do not represent that the same day necessarily can see the sun, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day
_{r}, be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation
_{h}, get Δ T for convenience of describing
_{h}unit be DEG C/m, the surface of Cable Structure is got " R Cable Structure surface point ", the Specific Principles getting " R Cable Structure surface point " describes in step b3, the temperature of this R Cable Structure surface point will be obtained below by actual measurement, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ", from the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, at the sea level elevation place that each is chosen, two points are at least chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along 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 in the shade the outer normal direction of the measurement Cable Structure chosen along the sunny slope outer normal direction and Cable Structure that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness, three points are no less than along each measurement Cable Structure along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, especially, along each, Cable Structure is measured for support cable and only gets a point along the direction of the Temperature Distribution of wall thickness, namely the temperature of the surface point of support cable is only measured, measure all temperature be selected a little, the temperature recorded 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 the direction of Cable Structure along the Temperature Distribution of wall thickness " to 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 have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, especially, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", the number temperature profile data at sea level elevation place " identical sea level elevation Cable Structure is along the temperature profile data of thickness " will chosen at each in this method ", measure temperature in Cable Structure location according to meteorology to require to choose a position, obtain meeting the temperature that meteorology measures the Cable Structure place environment of temperature requirement by the actual measurement of this position, in the onsite spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable, at the flat board of this position of sound production one piece of carbon steel material, 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 with dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, the nonsunny slope of reference plate is covered with insulation material, RealTime Monitoring is obtained 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 the temperature profile data of previously defined Cable Structure along thickness simultaneously, and RealTime Monitoring obtains meeting the temperature record that meteorology measures the Cable Structure place environment of temperature requirement simultaneously, the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be called environment maximum temperature difference, be designated as Δ T
_{emax}, calculated the rate of change of temperature about the time of Cable Structure place environment by Conventional mathematical by the temperature measured data sequence of Cable Structure place environment, this rate of change is also along with time variations, the measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be called reference plate maximum temperature difference, be designated as Δ T
_{pmax}, the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein called Cable Structure surface maximum temperature difference, be designated as Δ T
_{smax}, calculated the rate of change of temperature about the time of each Cable Structure surface point by Conventional mathematical by each Cable Structure surface temperature measured data sequence, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations, obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by RealTime Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and minimum temperature at the sea level elevation place that each is chosen, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T
_{tmax},
B3: survey calculation obtains Cable Structure steady temperature data, first, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day 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 Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the environment maximum error Δ T that survey calculation obtains above
_{emax}be not more than with reference to temperature difference per day Δ T
_{r}, and reference plate maximum temperature difference Δ T
_{pmax}Δ T is not more than after deducting 2 degrees Celsius
_{emax}, and Cable Structure surface maximum temperature difference Δ T
_{smax}be not more than Δ T
_{pmax}, one of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition, Section 3 condition is when obtaining Cable Structure steady temperature data, and 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, Section 4 condition is when obtaining Cable Structure steady temperature data, and 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, Section 5 condition is when obtaining 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 be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise, Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T
_{tmax}be not more than 1 degree Celsius, this method utilizes abovementioned six conditions, any one in following three kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in abovementioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, the second moment is moment of the Section 6 condition only met in abovementioned " condition relevant to determining to obtain moment of Cable Structure steady temperature data ", simultaneously the third moment meets the moment of the Section 1 in abovementioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition, when the mathematics moment obtaining Cable Structure steady temperature data is exactly one in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data, if the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data, the amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method, this method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, 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 obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution of the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steadystate surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steadystate surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steadystate 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 ", when the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%, Cable Structure surface comprises support cable surface, second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", 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 DEG C 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 obtained, gets Δ T for convenience of describing
_{h}unit be DEG C/m, be m for convenience of describing the unit getting Δ h, " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not 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, 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 geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshineduration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshineduration the most fully those surface points in Cable Structure,
C. the Cable Structure steady temperature data under original state are obtained according to " the temperature survey calculating method of the Cable Structure of this method " direct survey calculation, 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}"; Survey or consult reference materials and obtain the temperature variant physical and mechanical properties parameter of the various materials that Cable Structure uses; Initial Cable Structure steady temperature data vector T is obtained in actual measurement
_{o}synchronization, direct survey calculation obtains the Initial cable force of all support cables, composition Initial cable force vector F
_{o}; Obtain the length of all support cables when free state and Suo Li are 0 according to Cable Structure design data, completion data, in free state time crosssectional area and in free state time the weight of unit length, and the temperature of all support cables when obtaining these three kinds of data, utilize temperature variant physical function parameter and the mechanical property parameters of all support cables on this basis, conveniently physical computing obtains all support cables at initial Cable Structure steady temperature data vector T
_{o}the weight of the unit length of all support cables when the crosssectional area of all support cables and Suo Li are 0 when the length of all support cables, Suo Li are 0 when Suo Li under condition is 0, form the initial drift vector of support cable, the weight vector of the initial free unit length of initial free crosssectional area vector sum successively, the coding rule 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 and Initial cable force vector F
_{o}the coding rule of element identical; T is obtained in actual measurement
_{o}while, namely at the initial Cable Structure steady temperature data vector T of acquisition
_{o}the synchronization in moment, direct survey calculation obtains the measured data of initial Cable Structure, and the measured data of initial Cable Structure comprises the Nondestructive Testing Data of the health status expressing support cable, Cable Structure bearing initial generalized displacement measurement data, the initial value of all monitored amounts, the Initial cable 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; 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; Utilization can express the Nondestructive Testing Data of the health status of support cable and the initial generalized displacement measurement data of Cable Structure bearing sets up evaluation object initial damage vector d
_{o}, vectorial d
_{o}represent with initial mechanical Calculation Basis model A
_{o}the initial health of the evaluation object of the Cable Structure represented; 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 support cable corresponding to this element is intact, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to this element completely loses loadbearing capacity, if its numerical value is between 0 and 100%, then represent that this support cable loses the loadbearing capacity of corresponding proportion; If d
_{o}evaluation object corresponding to some elements be some generalized displacement components of some bearings, so d
_{o}the numerical value of this element represent the initial value of this generalized displacement component of this bearing; If there is no the Nondestructive Testing Data of support cable and other are when can express the data of the health status of support cable, or can think Cable Structure original state be not damaged without relaxed state time, vectorial d
_{o}in each element numerical value relevant to support cable get 0, if when there is no Cable Structure bearing initial generalized displacement measurement data or can think that the initial generalized displacement of Cable Structure bearing is 0, vectorial d
_{o}in each element numerical value relevant to Cable Structure generalized displacement of support get 0; Initial Cable Structure bearing generalized coordinate data refer to the support coordinate data under Cable Structure design point, and Cable Structure bearing initial generalized displacement measurement data refers to setting up initial mechanical Calculation Basis model A
_{o}time, the generalized displacement that Cable Structure bearing occurs relative to the bearing under Cable Structure design point; Bearing generalized coordinate comprises line amount and angular amount two kinds;
Temperature variant physical and mechanical properties parameter, the initial Cable Structure steady temperature data vector T of the various materials d. used according to the measured data of the design drawing of Cable Structure, asconstructed drawing and initial Cable Structure, the Nondestructive Testing Data of support cable, Cable Structure bearing initial generalized displacement measurement data, Cable Structure
_{o}with all Cable Structure data obtained with preceding step, set up the initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data "
_{o}, based on A
_{o}the Cable Structure that calculates calculates data must closely its measured data, 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}evaluation object health status with evaluation object initial damage vector d
_{o}represent; Corresponding to A
_{o}the initial value monitored amount initial value vector C of all monitored amount
_{o}represent; T
_{o}and d
_{o}a
_{o}parameter, by A
_{o}the initial value of all monitored amount that obtains of Mechanics Calculation result and C
_{o}the initial value of all monitored amount represented is identical, therefore alternatively C
_{o}by A
_{o}mechanics Calculation result composition, A in the method
_{o}, C
_{o}, d
_{o}and T
_{o}constant;
E. in the method, alphabetical i is except representing the place of number of steps significantly, and alphabetical i only represents cycle index, i.e. ith circulation; Ith circulation needs the current initial mechanical Calculation Basis model of Cable Structure that is that set up or that set up to be designated as current initial mechanical Calculation Basis model A when starting
^{i} _{o}, A
_{o}and A
^{i} _{o}count temperature parameter, can the Effect on Mechanical Properties of accounting temperature change to Cable Structure; When ith 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 one_to_one corresponding; The current initial damage vector of evaluation object that ith circulation needs when starting is designated as d
^{i} _{o}, d
^{i} _{o}cable Structure A when representing that this circulation starts
^{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 one_to_one corresponding; When ith circulation starts, the initial value of all monitored amounts, with monitored amount current initial value vector 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 one_to_one corresponding, monitored amount current initial value vector C
^{i} _{o}represent and correspond to A
^{i} _{o}the concrete numerical value of all monitored amount; T
^{i} _{o}and d
^{i} _{o}a
^{i} _{o}characterisitic parameter, C
^{i} _{o}by A
^{i} _{o}mechanics Calculation result composition; When first time, circulation started, A
^{i} _{o}be designated as A
^{1} _{o}, set up A
^{1} _{o}method for making A
^{1} _{o}equal A
_{o}; When first time, circulation started, T
^{i} _{o}be designated as T
^{1} _{o}, set up T
^{1} _{o}method for making T
^{1} _{o}equal T
_{o}; When first time, circulation started, d
^{i} _{o}be designated as d
^{1} _{o}, set up d
^{1} _{o}method for making d
^{1} _{o}equal d
_{o}; When first time, circulation started, 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 being walked to s step by f here; In Cable Structure military service process, the current data of Cable Structure steady temperature data is obtained, the current data composition current cable structure steady temperature data vector T of all " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement
^{i}, vector T
^{i}definition mode and vector T
_{o}definition mode identical, T
^{i}element and T
_{o}element one_to_one corresponding; Vector T is obtained in actual measurement
^{i}while, actual measurement obtains at acquisition current cable structure steady temperature data vector T
^{i}moment synchronization Cable Structure 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 one_to_one corresponding, represent that identical monitored amount is at not numerical value in the same time; Current cable structure steady temperature data vector T is obtained in actual measurement
^{i}synchronization, actual measurement obtain all M in Cable Structure
_{l}the rope force data of root support cable, all these rope force datas composition current cable force vector F, the element of vectorial F and vectorial F
_{o}the coding rule of element identical; Current cable structure steady temperature data vector T is obtained in actual measurement
^{i}synchronization, Actual measurement obtains all M
_{1}the volume coordinate of two supporting end points of root support cable, the difference of the volume coordinate component in the horizontal direction of two supporting end points is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables form current support cable two and support end points horizontal range vector, and current support cable two supports coding rule and the Initial cable force vector F of the element of end points horizontal range vector
_{o}the coding rule of element identical;
G. according to current cable structure steady temperature data vector T
^{i}, upgrade current initial mechanical Calculation Basis model A according to step g 1 to g3
^{i} _{o}, monitored amount current initial value vector C
^{i} _{o}with current initial Cable Structure steady temperature data vector T
^{i} _{o}, and evaluation object current initial damage vector d
^{i} _{o}remain unchanged;
G1. T is compared
^{i}with T
^{i} _{o}if, T
^{i}equal T
^{i} _{o}, then A
^{i} _{o}, C
^{i} _{o}and T
^{i} _{o}remain unchanged; Otherwise need to follow these steps to A
^{i} _{o}and T
^{i} _{o}upgrade;
G2. T is calculated
^{i}with T
_{o}difference, T
^{i}with T
_{o}difference be exactly the changes of current cable structure steady temperature data about initial Cable Structure steady temperature data, T
^{i}with T
_{o}difference represent with steady temperature change vector S, S equals T
^{i}deduct T
_{o}, S represents the change of Cable Structure steady temperature data;
G3. to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A
_{o}in Cable Structure apply temperature variation after obtain upgrade current initial mechanical Calculation Basis model A
^{i} _{o}, upgrade A
^{i} _{o}while, T
^{i} _{o}all elements numerical value also uses T
^{i}all elements numerical value correspondence replace, namely have updated T
^{i} _{o}, so just obtain and correctly correspond to A
^{i} _{o}t
^{i} _{o}; Now d
^{i} _{o}remain unchanged; As renewal A
^{i} _{o}after, A
^{i} _{o}the health status evaluation object of rope current initial damage vector d
^{i} _{o}represent, A
^{i} _{o}cable Structure steady temperature current cable structure steady temperature data vector T
^{i}represent, upgrade C
^{i} _{o}method be: when renewal A
^{i} _{o}after, obtain A by Mechanics Calculation
^{i} _{o}in all monitored amounts, current concrete numerical value, these concrete numerical value composition C
^{i} _{o};
H. at current initial mechanical Calculation Basis model A
^{i} _{o}basis on, carry out several times Mechanics Calculation according to step h1 to step h4, by calculate set up unit damage monitored numerical quantity unit change matrix Δ C
^{i}with unit damage or unit generalized displacement vector D
^{i} _{u};
H1., when ith circulation starts, directly Δ C is obtained by method listed by step h2 to step h4
^{i}and D
^{i} _{u}; In other moment, when in step g to A
^{i} _{o}after upgrading, Δ C must be regained by method listed by step h2 to step h4
^{i}and D
^{i} _{u}if, not to A in step g
^{i} _{o}upgrade, then directly proceed to step I herein and carry out followup work;
H2. at the current initial mechanical Calculation Basis model A of Cable Structure
^{i} _{o}basis on carry out several times Mechanics Calculation, calculation times numerically equals the quantity N of all evaluation objects, has N number of evaluation object just to have N calculating; According to the coding rule of evaluation object, calculate successively; Calculating hypothesis each time only has an evaluation object to increase unit damage or unit generalized displacement again on the basis of original damage or generalized displacement, 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 the generalized displacement component in a direction of a bearing, just suppose that this bearing increases unit generalized displacement again at this sense of displacement, use D
^{i} _{uk}record unit damage or the unit generalized displacement of this increase, wherein k represents the numbering of the evaluation object increasing unit damage or unit generalized displacement, D
^{i} _{uk}unit damage or unit generalized displacement vector D
^{i} _{u}an element, unit damage or unit generalized displacement vector D
^{i} _{u}the coding rule of element and vectorial d
_{o}the coding rule of element identical; The evaluation object increasing unit damage or unit generalized displacement in calculating each time is again different from during other time calculates the evaluation object increasing unit damage or unit generalized displacement again, calculate the current calculated value all utilizing mechanics method to calculate all monitored amount of Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector; When supposing that a kth evaluation object increases unit damage or unit generalized displacement again, use C
^{i} _{tk}represent corresponding " monitored amount calculation current vector "; When giving the element number of each vector in this step, same coding rule should be used with other vector in this method, to ensure any one element in this step in each vector, with in other vector, number identical element, have 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 one_to_one corresponding;
H3. the vectorial C calculated each time
^{i} _{tk}deduct vectorial C
^{i} _{o}obtain a vector, then after each element of this vector is calculated divided by this unit damage or unit generalized displacement numerical value supposed, obtain " unit change vector δ a C for monitored amount
^{i} _{k}"; N number of evaluation object is had just to have N number of " the unit change vector of monitored amount ";
H4. by this is N number of " the unit change vector of monitored amount " according to the coding rule of N number of evaluation object, " the unit damage monitored numerical quantity unit change matrix Δ C having N to arrange is formed successively
^{i}"; Unit damage monitored numerical quantity unit change matrix Δ C
^{i}each row correspond to a monitored amount unit change vector; Unit damage monitored numerical quantity unit change matrix Δ C
^{i}every a line correspond to the different unit change amplitude of same monitored amount when different evaluation object increases unit damage or unit generalized displacement; Unit damage monitored numerical quantity unit change matrix Δ C
^{i}the coding rule of row and vectorial d
_{o}the coding rule of element identical, unit damage monitored numerical quantity unit change matrix Δ C
^{i}the coding rule of coding rule and M monitored amount of row identical;
I. current nominal fatigue vector d is defined
^{i} _{c}with current actual damage vector d
^{i}, d
^{i} _{c}and d
^{i}element number equal the quantity of support cable, 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 health status of corresponding evaluation object, d
^{i} _{c}and d
^{i}with evaluation object initial damage vector d
_{o}element number 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 " monitored amount current initial value vector C
^{i} _{o}", " unit damage monitored numerical quantity unit change matrix Δ C
^{i}" and " current nominal fatigue vector d
^{i} _{c}" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, except d in formula 1
^{i} _{c}other outer amount is known, solves formula 1 and just can calculate current nominal fatigue vector d
^{i} _{c};
K. the current actual damage vector d utilizing formula 2 to express
^{i}a kth element d
^{i} _{k}with evaluation object current initial damage vector d
^{i} _{o}a kth element d
^{i} _{ok}with current nominal fatigue vector d
^{i} _{c}a kth element d
^{i} _{ck}between relation, calculate current actual damage vector d
^{i}all elements;
K=1 in formula 2,2,3 ..., N; d
^{i} _{k}represent the current actual health status of a kth evaluation object in ith circulation, d
^{i} _{k}represent when being 0 that a kth evaluation object is without health problem, d
^{i} _{k}numerical value represents when not being 0 that a kth evaluation object is the evaluation object of unsoundness problem, if this evaluation object is support cable, so a d in cable system
^{i} _{k}represent the order of severity of its current health problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, d
^{i} _{k}the degree of the lax or damage of this support cable of numerical response, if this evaluation object is generalized displacement component, so a d of a bearing
^{i} _{k}represent its current actual generalized displacement numerical value; So according to evaluation object current actual damage vector d
^{i}which support cable unsoundness problem and order of severity thereof can be defined, define which bearing and there occurs generalized displacement and numerical value thereof;
L. by current for evaluation object actual damage vector d
^{i}in with M
_{1}the M that root support cable is relevant
_{1}individual element takes out, composition support cable current actual damage vector d
^{ca}, support cable current actual damage vector d
^{ca}the coding rule of element and Initial cable force vector F
_{o}the coding rule of element identical; Support cable current actual damage vector d
^{ca}h element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., M
_{1}; Support cable current actual damage vector d
^{ca}middle numerical value be not 0 element correspond to the support cable of unsoundness problem, from the support cable of these unsoundness problems, identify damaged cable, remaining is exactly slack line, support cable current actual damage vector d
^{ca}in correspond to slack line element numerical expression be the current actual equivalent damage degree with slack line relax level mechanic equivalent;
M. utilize at current cable structure steady temperature data vector T
^{i}slack line under condition, that identify in the 1st step and with support cable current actual damage vector d
^{ca}these slack lines of expressing, with the current actual equivalent damage degree of its relax level mechanic equivalent, utilize f step obtain at current cable structure steady temperature data vector T
^{i}current cable force vector F under condition and current support cable two support end points horizontal range vector, utilize 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, the vectorial F of Initial cable force
_{o}, utilize current cable structure steady temperature data vector T
^{t}represent support cable current steady state temperature data, utilize c step obtain at initial Cable Structure steady temperature data vector T
_{o}the support cable initial steady state temperature data represented, the temperature variant physical and mechanical properties parameter of the various materials utilizing the Cable Structure obtained in c step to use, count the impact of temperature variation on support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanic equivalent calculate slack line, with the relax level of current actual equivalent damage degree equivalence, the mechanical condition of equivalence is: one, two equivalences rope without lax identical with the mechanics parameters of initial drift during not damaged, geometrical property parameter, density and material; Two, after lax or damage, the slack line of two equivalences and the Suo Li of damage rope be out of shape after overall length identical; When meeting abovementioned two equivalent conditions, the such mechanics function of two support cables in Cable Structure is exactly identical, if after namely replacing damaged cable with the slack line of equivalence, any change can not occur Cable Structure, and vice versa; Try to achieve according to aforementioned mechanic equivalent condition the relax level that those are judged as slack line, relax level is exactly the knots modification of support cable drift, namely determines the long adjustment amount of those ropes that need adjust the support cable of Suo Li; So just achieve lax identification and the nondestructive tests of support cable; During calculating, institute's demand power is provided by current cable force vector F corresponding element;
N. current nominal fatigue vector d is tried to achieve
^{i} _{c}after, set up mark vector B according to formula 3
^{i}, formula 4 gives mark vector B
^{i}the definition of a kth element;
Element B in formula 4
^{i} _{k}mark vector B
^{i}a kth element, D
^{i} _{uk}unit damage or unit generalized displacement vector D
^{i} _{u}a kth element, d
^{i} _{ck}evaluation object current nominal fatigue vector d
^{i} _{c}a kth element, they all represent the relevant information of a kth evaluation object, k=1 in formula 4,2,3 ..., N;
If o. mark vector B
^{i}element be 0 entirely, then get back to step f continue this circulation; If mark vector B
^{i}element be not 0 entirely, then enter next step, i.e. step p;
P. calculate next time according to formula 5, evaluation object current initial damage vector d namely needed for the ith+1 time circulation
^{i+1} _{o}each element;
D in formula 5
^{i+1} _{ok}the evaluation object current initial damage vector d next time, namely needed for the ith+1 time circulation
^{i+1} _{o}a kth element, d
^{i} _{ok}this, i.e. the evaluation object of ith circulation current initial damage vector d
^{i} _{o}a kth element, D
^{i} _{uk}unit damage or the unit generalized displacement vector D of ith circulation
^{i} _{u}a kth element, B
^{i} _{k}the mark vector B of ith circulation
^{i}a kth element, k=1 in formula 5,2,3 ..., N;
Q. take off once, namely the ith+1 time circulation needed for current initial Cable Structure steady temperature data vector T
^{i+1} _{o}equal the current initial Cable Structure steady temperature data vector T of ith circulation
^{i} _{o};
R. at initial mechanical Calculation Basis model A
_{o}basis on, to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, then makes the health status of rope be d
^{i+1} _{o}after obtain be exactly next time, namely the ith+1 time circulation needed for Mechanics Calculation benchmark model A
^{i+1}; Obtain A
^{i+1}after, obtain A by Mechanics Calculation
^{i+1}in all monitored amounts, current concrete numerical value, these the monitored amounts of concrete numerical value composition next time, namely needed for the ith+1 time circulation current initial value vector C
^{i+1} _{o};
S. get back to step f, start to circulate next time.
Beneficial effect: when the temperature field of Cable Structure is subject to affecting of the factor such as sunshine and environment temperature, the temperature field of Cable Structure is constantly change, the change of temperature field of Cable Structure must affect the monitored amount of Cable Structure, only have and monitored amount is rejected could carry out rational cable structure health monitoring based on monitored amount by temperature profile effect part, and the temperature field measurement of Cable Structure and calculating are very complicated, this method discloses to comprise and is a kind ofly suitable for the simple of cable structure health monitoring, economical, feasible, the cable structure health monitoring method of efficient Cable Structure Calculation Method of Temperature Field, adopt this method when generalized displacement appears in Cable Structure bearing, when many ropes of Cable Structure are synchronously impaired or lax, and the temperature of Cable Structure along with time variations time, disclose a kind of employing linear segmented and approach nonlinear method, very monitor assessment can identify damaged cable, slack line and generalized displacement of support, the effective health monitoring of system and method disclosed in this method to Cable Structure is highly profitable.
Embodiment
When temperature variation, for damaged cable and the generalized displacement of support identification of Cable Structure, this method discloses a kind of system and method can monitoring the health status identifying each evaluation object in Cable Structure rationally and effectively.The following describes of embodiment of this method is in fact only exemplary, and object is never the application or the use that limit this method.
This method adopts a kind of algorithm, and this algorithm is for identifying damaged cable, slack line and generalized displacement of support.During concrete enforcement, the following step is the one in the various steps that can take.
The first step: set the quantity sum of the quantity of the support cable of Cable Structure and the generalized displacement of support component of Cable Structure as N.For sake of convenience, this method unitedly calls evaluated support cable and generalized displacement of support to be " evaluation object ", total N number of evaluation object.To evaluation object serial number, this numbering will be used for generating vector sum matrix in subsequent step, and this method variable k represents this numbering, k=1, and 2,3 ..., N.
If total M in cable system
_{1}root support cable, Cable Structure rope force data comprises this M
_{1}the Suo Li of root support cable, obvious M
_{1}be less than the quantity N of evaluation object.
Specify when determining hybrid monitoring by the support cable of monitored Suo Li, the monitored rope force data of Cable Structure specifies Q rope force data of support cable to describe by Q in Cable Structure, and the change of Cable Structure Suo Li is exactly the change of the Suo Li of all appointment support cables.Each total Q cable force measurement value or calculated value characterize the rope force information of Cable Structure.Q is one and is not less than 0, is not more than M
_{1}integer.When reality selectes the rope of monitored Suo Li, the rope that those Suo Li can be selected to be easy to measure is monitored rope.
Specify when determining hybrid monitoring by the measured point of monitored strain, the monitored strain data of Cable Structure can by K in Cable Structure
_{2}the L of individual specified point and each specified point
_{2}the strain of individual assigned direction describes, and the change of Cable Structure strain data is exactly K
_{2}the change of all tested strain 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}it is an integer being not less than 0.The measured point of monitored strain can be exactly a point near the fixed endpoint (being such as the stiff end of dragline on bridge of cablestayed bridge) of each root support cable by each, this specified point can also be a point near Cable Structure bearing, this point should not be generally stress concentration point, to avoid occurring excessive strain measurement value, the fixed endpoint of the rope of the monitored Suo Li specified when these points should not be also generally all hybrid monitorings or in its vicinity.
Specify when determining hybrid monitoring by the measured point of monitored angle, the monitored angledata of Cable Structure is by K in Cable Structure
_{3}individual specified point, cross the L of each specified point
_{3}the H of individual appointment straight line, each appointment straight line
_{3}individual angle coordinate component describes, and the change of Cable Structure angle is exactly the change of all specified points, all appointment straight line, all angle coordinate component of specifying.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}it is an integer being not less than 0.Each specified point can be exactly the fixed endpoint (being such as the stiff end of dragline on bridge floor of cablestayed bridge) of each root support cable or a point near it, this specified point can also be a point near Cable Structure bearing, and the point of monitored angledata generally should all not be chosen as " fixed endpoint of the rope of the monitored Suo Li specified in hybrid monitoring or point in its vicinity " and " point of the monitored strain of specifying in hybrid monitoring or point in its vicinity "; Only can measure one at each specified point and specify an angle coordinate of straight line, such as measuring the Cable Structure surface normal of specified point or the tangent line angle coordinate relative to acceleration of gravity direction, is in fact exactly measurement of dip angle here.
Specify when determining hybrid monitoring by monitored shape data, the monitored shape data of Cable Structure is by K in Cable Structure
_{4}the L of individual specified point and each specified point
_{4}the volume coordinate of individual assigned direction describes, and the change of Cable Structure shape data is exactly K
_{4}the change of all coordinate components of individual specified point.Each total M
_{4}individual coordinates measurements or calculated value characterize Cable Structure shape, M
_{4}for K
_{4}and L
_{4}longpending.M
_{4}it is an integer being not less than 0.Each specified point can be exactly the fixed endpoint (being such as the stiff end of dragline on bridge of cablestayed bridge) of each root support cable, and this specified point can also be a point near Cable Structure bearing, or is exactly directly Cable Structure seat pivot; Here selected point being monitored should all not selected " fixed endpoint of the rope of the monitored Suo Li specified in hybrid monitoring or point in its vicinity ", " point of the monitored strain of specifying in hybrid monitoring or point in its vicinity " and " point of the monitored angledata of specifying in hybrid monitoring or point in its vicinity ".
Comprehensive abovementioned monitored amount, whole Cable Structure is total M monitored amount with regard to hybrid monitoring, and M is Q, M
_{2}, M
_{3}and M
_{4}sum, definition parameter K, K is Q, K
_{2}, K
_{3}and K
_{4}sum, K and M must not be less than the quantity N of evaluation object.Because M monitored amount is dissimilar, so this method is called " during temperature variation the damaged cable of hybrid monitoring and support translation identification method ".For simplicity, in the method by monitored all parameters of Cable Structure " during the hybrid monitoring " listed by this step referred to as " monitored amount ".To M monitored amount serial number, this numbering will be used for generating vector sum matrix in subsequent step.This method represents this numbering, j=1,2,3 with variable j ..., M.
Determine " the temperature survey calculating method of the Cable Structure of this method ", the method concrete steps are as follows:
A walks: inquiry or actual measurement (can be measured by ordinary temperature measuring method, thermal resistance is such as used to measure) obtain the temperature variant thermal conduction study parameter of environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, asconstructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model (such as finite element model) of Cable Structure.Inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day
_{r}, be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation
_{h}, get Δ T for convenience of describing
_{h}unit be DEG C/m, the surface of Cable Structure is got " R Cable Structure surface point ", the Specific Principles getting " R Cable Structure surface point " describes in step b3, the temperature of this R Cable Structure surface point will be obtained below by actual observation record, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ".From the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, if the sea level elevation of such as Cable Structure is between 0m to 200m, so can choose height above sea level 0m, 50m, 100m and height above sea level 200m, crossing with Cable Structure surface with imaginary surface level at the sea level elevation place that each is chosen, obtain intersection, surface level is crossing with Cable Structure obtains cross surface, intersection is the outer edge line of cross surface, 6 points are chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along 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 6 directions of the measurement Cable Structure chosen along 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, Cable Structure surrounding environment etc. determines the sunny slope of Cable Structure and in the shade, the sunny slope of Cable Structure and in the shade face are the parts on the surface of Cable Structure, at the sea level elevation place that each is chosen, aforementioned intersection respectively has one section in sunny slope and in the shade, these two sections of intersection respectively have a mid point, cross the outer normal that these two mid points get Cable Structure, these two outer normals are called the sunny slope outer normal of Cable Structure and in the shade outer normal of Cable Structure by this method, these two outer normal directions are called the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure by this method, the outer normal of obvious sunny slope and the outer normal of in the shade all crossing with aforementioned intersection, also two intersection points are just had, intersection is divided into two line segments by these two intersection points, 2 points are got respectively on two line segments, totally 4 points, each line segment in two of intersection line segments is divided into equal 3 sections of length by taken point, the outer normal on Cable Structure surface is got at these 4 some places, the outer normal on 6 Cable Structure surfaces is just have chosen so altogether at each selected sea level elevation place, the direction of 6 outer normals is exactly " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ".There are two intersection points on the surface of each " measures the direction of Cable Structure along the Temperature Distribution of wall thickness " line and Cable Structure, if Cable Structure is hollow, these, two intersection points are on Cable Structure outside surface, another on an internal 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, straightline segment is chosen three points again, this straightline segment is divided into four sections by these three points, measure two end points of Cable Structure at these three points chosen and straightline segment, amount to the temperature of 5 points, concrete can first hole in Cable Structure, temperature sensor is embedded in this 5 some places, especially, can not hole in support cable, support cable is only measured to the temperature of support cable surface point, in any case, the temperature recorded all 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 the direction of Cable Structure along the Temperature Distribution of wall thickness " to 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 have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, especially, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", the number temperature profile data at sea level elevation place " identical sea level elevation Cable Structure is along the temperature profile data of thickness " will chosen at each in this method ".Measuring temperature in Cable Structure location according to meteorology to require to choose a position, meeting obtaining in this position actual observation record the temperature that meteorology measures the Cable Structure place environment of temperature requirement, in the onsite spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable (as long as there was sunrise the same day, this position just should by sunlight), at the flat board (square plate that the wide 3mm of such as 30cm is thick) of this position of sound production one piece of carbon steel material (such as 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 measures the wooden thermometer screen required, the one side of this reference plate on the sunny side, be called that sunny slope (such as, time 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 with dark color (being conducive to accepting solar irradiation), the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, the nonsunny slope of reference plate is covered with insulation material (the thick calcium carbonate insulation material of such as 5mm), RealTime Monitoring record is obtained the temperature of the sunny slope of reference plate.
B walks, RealTime Monitoring (can be measured by ordinary temperature measuring method, thermal resistance is such as used to measure, such as every 10 minutes survey records temperature data) record R the Cable Structure surface temperature measured data obtaining abovementioned R Cable Structure surface point, RealTime Monitoring (can be measured by ordinary temperature measuring method simultaneously, thermal resistance is such as used to measure, such as every 10 minutes survey records temperature data) obtain the temperature profile data of previously defined Cable Structure along thickness, RealTime Monitoring (can be measured by ordinary temperature measuring method simultaneously, such as in the wooden thermometer screen meeting meteorology temperature measurement requirement, lay thermal resistance and measure temperature, such as every 10 minutes survey records temperature data) record the temperature record obtaining the Cable Structure place environment meeting the requirement of meteorology measurement temperature, (can be measured by ordinary temperature measuring method by RealTime Monitoring, such as in the wooden thermometer screen meeting meteorology temperature measurement requirement, lay thermal resistance and measure temperature, such as every 10 minutes survey records temperature data) record obtains being carved at sunrise the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes 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 of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be designated as Δ T
_{emax}, (such as first the temperature measured data sequence of Cable Structure place environment is carried out curve fitting by Conventional mathematical calculating by the temperature measured data sequence of Cable Structure place environment, then by asking curve to the derivative of time or by asking on curve each point corresponding to survey record data time by numerical method to the rate of change of time) obtain the rate of change of temperature about the time of Cable Structure place environment, this rate of change is also along with time variations, (can be measured by ordinary temperature measuring method by RealTime Monitoring, such as use the temperature of the dull and stereotyped sunny slope of thermal resistance witness mark, such as every 10 minutes survey records temperature data) obtain being carved at sunrise the same day measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be designated as Δ T
_{pmax}, (can be measured by ordinary temperature measuring method by RealTime Monitoring, thermal resistance is such as used to measure Cable Structure surface point, such as every 10 minutes survey records temperature data) record obtains being carved at sunrise the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes the same day, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein designated as Δ T
_{smax}, (such as first each Cable Structure surface temperature measured data sequence is carried out curve fitting by Conventional mathematical calculating by each Cable Structure surface temperature measured data sequence, then by asking curve to the derivative of time or by asking on curve each point corresponding to survey record data time by numerical method to the rate of change of time) obtain the rate of change of temperature about the time of each Cable Structure surface point, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations.Obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by RealTime Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and minimum temperature at the sea level elevation place that each is chosen, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T
_{tmax}.
C walks, and survey calculation obtains Cable Structure steady temperature data; First, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology; The a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, Δ T
_{pmax}with Δ T
_{smax}all be not more than 5 degrees Celsius; The b condition that Section 2 must meet be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the Δ T that survey calculation obtains above
_{emax}be not more than with reference to temperature difference per day Δ T
_{r}, and the Δ T that survey calculation obtains above
_{pmax}deduct 2 degrees Celsius and be not more than Δ T
_{emax}, and the Δ T that survey calculation obtains above
_{smax}be not more than Δ T
_{pmax}; One of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition; Section 3 condition is when obtaining Cable Structure steady temperature data, and 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; Section 4 condition is when obtaining Cable Structure steady temperature data, and 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; Section 5 condition is when obtaining 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 be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise; Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T
_{tmax}be not more than 1 degree Celsius.This method utilizes abovementioned six conditions, any one in following three kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in abovementioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, the second moment is moment of the Section 6 condition only met in abovementioned " condition relevant to determining to obtain moment of Cable Structure steady temperature data ", simultaneously the third moment meets the moment of the Section 1 in abovementioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition, when the mathematics moment obtaining Cable Structure steady temperature data is exactly one in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data, if the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data, the amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method, this method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly moment of the acquisition Cable Structure steady temperature data of this method, then, according to Cable Structure heat transfer characteristic, utilize R Cable Structure surface temperature measured data and " HBE Cable Structure is along the thickness temperature measured data " in the moment obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model (such as finite element model) of Cable Structure, the Temperature Distribution that (such as finite element method) obtains the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steadystate surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steadystate surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steadystate 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 ".When the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%; Cable Structure surface comprises support cable surface; Second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", 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 DEG C 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 obtained, gets Δ T for convenience of describing
_{h}unit be DEG C/m, be m for convenience of describing the unit getting Δ h; " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not 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; 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 geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshineduration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshineduration the most fully those surface points in Cable Structure.
Second step: set up initial mechanical Calculation Basis model A
_{o}.
When Cable Structure is completed, or before setting up health monitoring (damaged cable identification) system, obtain " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " survey calculation (to measure by ordinary temperature measuring method, thermal resistance is such as used to measure), " Cable Structure steady temperature data " now use vector T
_{o}represent, be called initial Cable Structure steady temperature data vector T
_{o}.T is obtained in actual measurement
_{o}while, namely at the synchronization in the moment of the initial Cable Structure steady temperature data vector of acquisition, use the direct survey calculation of conventional method to obtain the initial value of all monitored amount of Cable Structure, form monitored amount initial value vector C
_{o}.
Specifically the synchronization in moment of soandso Cable Structure steady temperature data vector such as (such as initial or current) can be being obtained according to following method in this method, soandso method survey calculation is used to obtain the data of the monitored amount of soandso measured amount (all monitored amount of such as Cable Structure): (to comprise the temperature of Cable Structure place environment in survey record temperature, the temperature of the sunny slope of reference plate and Cable Structure surface temperature) while, such as every 10 minutes survey records temperature, so simultaneously equally also every 10 minutes the monitored amount of soandso measured amount of survey record (all monitored amount of such as Cable Structure) data.Once determine the moment obtaining Cable Structure steady temperature data, so just be called the synchronization in the moment obtaining Cable Structure steady temperature data, the data of the monitored amount of soandso measured amount using soandso method survey calculation method to obtain with the data of the monitored amount of soandso measured amount (all monitored amount of such as Cable Structure) of the moment synchronization of acquisition Cable Structure steady temperature data.
Conventional method (consult reference materials or survey) is used to obtain temperature variant physical parameter (such as thermal expansivity) and the mechanical property parameters (such as elastic modulus, Poisson ratio) of the various materials that Cable Structure uses; Initial Cable Structure steady temperature data vector T is obtained in actual measurement
_{o}synchronization, direct survey calculation obtains the Initial cable force of all support cables, composition Initial cable force vector F
_{o}; Obtain the length of all support cables when free state and Suo Li are 0 according to Cable Structure design data, completion data, in free state time crosssectional area and in free state time the weight of unit length, and the temperature of all support cables when obtaining these three kinds of data, utilize temperature variant physical function parameter and the mechanical property parameters of all support cables on this basis, conveniently physical computing obtains all support cables at initial Cable Structure steady temperature data vector T
_{o}when when Suo Li under condition is 0, the length of all support cables, Suo Li are 0, the crosssectional area of all support cables and Suo Li are 0, the weight of the unit length of all support cables, forms the initial drift vector l of support cable successively
_{o}, initial free crosssectional area vector A
_{o}with the weight vector ω of initial free unit length
_{o}, the initial drift vector l of support cable
_{o}, initial free crosssectional area vector A
_{o}with the weight vector ω of initial free unit length
_{o}the coding rule of element and Initial cable force vector F
_{o}the coding rule of element identical.Initial Cable Structure steady temperature data vector T is obtained at Actual measurement
_{o}while, namely at the synchronization in the moment of acquisition Cable Structure steady temperature data, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Actual measurement data of Cable Structure comprise the measured data such as data, the initial geometric data of Cable Structure, rope force data, drawbar pull data, initial Cable Structure bearing generalized coordinate data, the initial generalized displacement measurement data of Cable Structure bearing, Cable Structure modal data, Cable Structure strain data, Cable Structure angle measurement data, Cable Structure volume coordinate measurement data that the Nondestructive Testing Data of support cable etc. can express the health status of rope.Initial Cable Structure bearing generalized coordinate data refer to the support coordinate data under Cable Structure design point, and Cable Structure bearing initial generalized displacement measurement data refers to setting up initial mechanical Calculation Basis model A
_{o}time, the generalized displacement that Cable Structure bearing occurs relative to the bearing under Cable Structure design point.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 Cable Structure, and object is the geometric properties according to these coordinate data determination Cable Structure.For cablestayed bridge, initial geometric data can be the spatial data that the spatial data of the end points of all ropes adds some points on bridge two ends, socalled bridge type data that Here it is.Utilize the Nondestructive Testing Data etc. of support cable can express the data of the health status of support cable and the initial generalized displacement measurement data of Cable Structure bearing sets up evaluation object initial damage vector d
_{o}(such as formula (1) Suo Shi), uses d
_{o}represent that Cable Structure is (with initial mechanical Calculation Basis model A
_{o}represent) the initial health of evaluation object.If there is no the Nondestructive Testing Data of support cable and other are when can express the data of the health status of support cable, or can think Cable Structure original state be not damaged without relaxed state time, vectorial d
_{o}in each element numerical value relevant to support cable get 0, if when there is no Cable Structure bearing initial generalized displacement measurement data or can think that the initial generalized displacement of Cable Structure bearing is 0, vectorial d
_{o}in each element numerical value relevant to Cable Structure generalized displacement of support get 0.The temperature variant physical and mechanical properties parameter of the various materials utilizing the measured data of the design drawing of Cable Structure, asconstructed drawing and initial Cable Structure, the Nondestructive Testing Data of support cable, Cable Structure to use and initial Cable Structure steady temperature data vector T
_{o}, utilize mechanics method (such as finite element method) to count " Cable Structure steady temperature data " and set up initial mechanical Calculation Basis model A
_{o}.
No matter which kind of method to obtain initial mechanical Calculation Basis model A by
_{o}, count " Cable Structure steady temperature data " (i.e. initial Cable Structure steady temperature data vector T
_{o}), based on A
_{o}the Cable Structure that calculates calculates data must closely its measured data, and error generally must not be greater than 5%.Like this can utility A
_{o}the Suo Li calculated under the analog case of gained calculates data, strain calculation data, Cable Structure shapometer count certificate and displacement meter counts certificate, Cable Structure angledata, Cable Structure spatial data etc., measured data when reliably truly occurring close to institute's analog case.Model A
_{o}the health status evaluation object initial damage vector d of middle evaluation object
_{o}represent, the initial Cable Structure steady temperature data vector T of Cable Structure steady temperature data
_{o}represent.Due to based on A
_{o}the initial value (actual measurement obtains) of the evaluation calculating all monitored amounts closely all monitored amounts, 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 health status with evaluation object initial damage vector d
_{o}represent; Corresponding to A
_{o}the initial value monitored amount initial value vector C of all monitored amount
_{o}represent.T
_{o}and d
_{o}a
_{o}parameter, C
_{o}also can by A
_{o}mechanics Calculation result composition.
3rd step: in the method, alphabetical i is except representing the place of number of steps significantly, and alphabetical i only represents cycle index, i.e. ith circulation; Ith circulation needs the current initial mechanical Calculation Basis model of Cable Structure that is that set up or that set up to be designated as current initial mechanical Calculation Basis model A when starting
^{i} _{o}, A
_{o}and A
^{i} _{o}count temperature parameter, can the Effect on Mechanical Properties of accounting temperature change to Cable Structure; When ith 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 one_to_one corresponding; The current initial damage vector of evaluation object that ith circulation needs when starting is designated as d
^{i} _{o}, d
^{i} _{o}cable Structure A when representing that this circulation starts
^{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 one_to_one corresponding; When ith circulation starts, the initial value of all monitored amounts, with monitored amount current initial value vector 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 one_to_one corresponding, monitored amount current initial value vector C
^{i} _{o}represent and correspond to A
^{i} _{o}the concrete numerical value of all monitored amount; T
^{i} _{o}and d
^{i} _{o}a
^{i} _{o}characterisitic parameter; C
^{i} _{o}by A
^{i} _{o}mechanics Calculation result composition; When first time, circulation started, A
^{i} _{o}be designated as A
^{1} _{o}, set up A
^{1} _{o}method for making A
^{1} _{o}equal A
_{o}; When first time, circulation started, T
^{i} _{o}be designated as T
^{1} _{o}, set up T
^{1} _{o}method for making T
^{1} _{o}equal T
_{o}; When first time, circulation started, d
^{i} _{o}be designated as d
^{1} _{o}, set up d
^{1} _{o}method for making d
^{1} _{o}equal d
_{o}; When first time, circulation started, C
^{i} _{o}be designated as C
^{1} _{o}, set up C
^{1} _{o}method for making C
^{1} _{o}equal C
_{o}.
4th step: the hardware components of pass line structural healthy monitoring system.Hardware components at least comprises: monitored amount monitoring system (such as measuring subsystem, signal conditioner etc. containing measurement of angle subsystem, cable force measurement subsystem, strain measurement subsystem, volume coordinate), Cable Structure temperature monitoring system (containing temperature sensor, signal conditioner etc.) and Cable Structure ambient temperature measurement system (containing temperature sensor, signal conditioner etc.), signal (data) collector, computing machine and communication alert equipment.Each monitored amount, each temperature must arrive by monitored system monitoring, monitoring system by the Signal transmissions that monitors to signal (data) collector; Signal is delivered to computing machine through signal picker; The health monitoring software of the evaluation object running Cable Structure is then responsible for by computing machine, comprises the signal that the transmission of tracer signal collector comes; When monitoring evaluation object health status and changing, computer control communication panalarm is reported to the police to monitor staff, owner and (or) the personnel that specify.
5th step: establishment installation and operation temperature variation angle monitor problem cable and generalized displacement of support progressivetype recognition method system software on computers, the function (i.e. all work that can complete with computing machine in this specific implementation method) such as monitoring, record, control, storage, calculating, notice, warning that this software will complete this method " temperature variation angle monitor problem cable and generalized displacement of support progressivetype recognition method " required by task and wants.
6th step: step starts circulation running thus, in Cable Structure military service process, the current data of Cable Structure steady temperature data is obtained, the current data composition current cable structure steady temperature data vector T of all " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement
^{i}, vector T
^{i}definition mode and vector T
_{o}definition mode identical, T
^{i}element and T
_{o}element one_to_one corresponding; In actual measurement vector T
^{i}while, namely at acquisition current cable structure steady temperature data vector T
^{i}the synchronization in moment, actual measurement obtains the currency of all monitored amounts in Cable Structure, and 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 one_to_one corresponding, represent that identical monitored amount is at not numerical value in the same time.Current cable structure steady temperature data vector T is obtained in actual measurement
^{i}synchronization, actual measurement obtain all M in Cable Structure
_{1}the rope force data of root support cable, all these rope force datas composition current cable force vector F, the element of vectorial F and vectorial F
_{o}the coding rule of element identical; Current cable structure steady temperature data vector T is obtained in actual measurement
^{i}synchronization, Actual measurement obtains all M
_{1}the volume coordinate of two of root support cable supporting end points, the difference of the volume coordinate of two supporting end points component is in the horizontal direction 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 and support end points horizontal range vector l
^{t} _{x}, current support cable two supports end points horizontal range vector l
^{t} _{x}the coding rule of element and Initial cable force vector F
_{o}the coding rule of element identical.
7th step: obtaining current cable structure steady temperature data vector T
^{i}after, compare T
^{i}and T
^{i} _{o}if, T
^{i}equal T
^{i} _{o}, then do not need A
^{i} _{o}and T
^{i} _{o}upgrade, otherwise need current initial mechanical Calculation Basis model A
^{i} _{o}, current initial Cable Structure steady temperature data vector T
^{i} _{o}with monitored amount current initial value vector C
^{i} _{o}upgrade, and evaluation object current initial damage vector d
^{i} _{o}remain unchanged, update method follows these steps to a to step c and carries out:
A. T is calculated
^{i}with T
_{o}difference, T
^{i}with T
_{o}difference be exactly the changes of current cable structure steady temperature data about initial Cable Structure steady temperature data, T
^{i}with T
_{o}difference represent with steady temperature change vector S, S equals T
^{i}deduct T
_{o}, S represents the change of Cable Structure steady temperature data.
B. to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A
_{o}in Cable Structure apply the current initial mechanical Calculation Basis model A that obtains after temperature variation upgrading
^{i} _{o}, upgrade A
^{i} _{o}while, T
^{i} _{o}all elements numerical value also uses T
^{i}all elements numerical value correspondence replace, namely have updated T
^{i} _{o}, so just obtain and correctly correspond to A
^{i} _{o}t
^{i} _{o}; Now d
^{i} _{o}remain unchanged; As renewal A
^{i} _{o}after, A
^{i} _{o}the health status evaluation object of evaluation object current initial damage vector d
^{i} _{o}represent, A
^{i} _{o}cable Structure steady temperature current cable structure steady temperature data vector T
^{i}represent, upgrade C
^{i} _{o}method be: when renewal A
^{i} _{o}after, obtain A by Mechanics Calculation
^{i} _{o}in all monitored amounts, current concrete numerical value, these concrete numerical value composition C
^{i} _{o};
8th step: at current initial mechanical Calculation Basis model A
^{i} _{o}basis on, carry out several times Mechanics Calculation according to step a to steps d, by calculate set up unit damage monitored numerical quantity unit change matrix Δ C
^{i}with unit damage or unit generalized displacement vector D
^{i} _{u}.
A., when ith circulation starts, directly Δ C is obtained by method listed by step b to steps d
^{i}and D
^{i} _{u}; In other moment, when in the 7th step to A
^{i} _{o}after upgrading, Δ C must be regained by method listed by step b to steps d
^{i}and D
^{i} _{u}if, not to A in the 7th step
^{i} _{o}upgrade, then directly proceed to the 9th step herein and carry out followup work.
B. at current initial mechanical Calculation Basis 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, calculation times numerically equals the quantity N of all evaluation objects, has N number of evaluation object just to have N calculating; Calculate hypothesis each time and only have an evaluation object at vectorial d
^{i} _{o}there is unit damage or unit generalized displacement in the basis of the health status of the evaluation object represented, 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}the basis that this support cable represented has a damage there is again unit damage (such as getting 5%, 10%, 20% or 30% equivalent damage is unit damage), if this evaluation object is the generalized displacement component in a direction of a bearing, just suppose this bearing at this sense of displacement at vectorial d
^{i} _{o}there is unit generalized displacement again (if such as this evaluation object is the translational component in the x direction of a bearing in the basis that this bearing represented has a generalized displacement, just suppose that this bearing has unit line displacement in x direction, if this evaluation object is the generalized displacement component around xaxis of a bearing, just suppose that this bearing is around the generalized displacement of xZhou You unit), use D
^{i} _{uk}record this unit damage or unit generalized displacement, wherein k represents the numbering of the evaluation object that unit damage or unit generalized displacement occur, D
^{i} _{uk}unit damage or unit generalized displacement vector D
^{i} _{u}an element, unit damage or unit generalized displacement vector D
^{i} _{u}the coding rule of element and vectorial d
_{o}the coding rule of element identical; Occur in calculating each time that the evaluation object of unit damage or unit generalized displacement is different from during other time calculates the evaluation object occurring unit damage or unit generalized displacement, calculate the current calculated value all utilizing mechanics method to calculate all monitored amount of Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector; When supposing that a kth evaluation object has unit damage or unit generalized displacement, available C
^{i} _{tk}represent corresponding " monitored amount calculation current vector "; When giving the element number of each vector in this step, same coding rule should be used with other vector in this method, to ensure any one element in this step in each vector, with in other vector, number identical element, have 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 one_to_one corresponding.
C. the vectorial C calculated each time
^{i} _{tk}deduct vectorial C
^{i} _{o}obtain a vector, then by each element of this vector divided by the unit damage supposed in this calculating or unit generalized displacement numerical value D
^{i} _{uk}after obtain " the unit change of a monitored amount vector δ C
^{i} _{k}"; N number of evaluation object is had just to have N number of " the unit change vector of monitored amount ".
D. by this is N number of " the unit change vector of monitored amount " according to the coding rule of N number of evaluation object, " the unit damage monitored numerical quantity unit change matrix Δ C having N to arrange is formed successively
^{i}"; Unit damage monitored numerical quantity unit change matrix Δ C
^{i}each row correspond to a monitored amount unit change vector; Unit damage monitored numerical quantity unit change matrix Δ C
^{i}every a line correspond to the different unit change amplitude of same monitored amount when different evaluation object increases unit damage or unit generalized displacement; Unit damage monitored numerical quantity unit change matrix Δ C
^{i}the coding rule of row and vectorial d
_{o}the coding rule of element identical, unit damage monitored numerical quantity unit change matrix Δ C
^{i}the coding rule of coding rule and M monitored amount of row identical.
9th step: set up linear relationship error vector e
^{i}with vectorial g
^{i}.Utilize (" the monitored amount current initial value vector C of data above
^{i} _{o}", " unit damage monitored numerical quantity unit change matrix Δ C
^{i}"); while the 8th step calculates each time; namely calculate each time hypothesis evaluation object in only have an evaluation object increase unit damage or unit generalized displacement while; when hypothesis kth (k=1,2,3; ..., N), when individual evaluation object increases unit damage or unit generalized displacement, calculate composition injury vector each time, use d
^{i} _{tk}represent this injury vector, corresponding monitored amount calculation current vector is C
^{i} _{tk}(see the 8th step), injury vector 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 unit generalized displacement value of the evaluation object of unit damage or unit generalized displacement, d
^{i} _{tk}the numerical value of other element get 0, that be not numbering and the supposition of the element of 0 increase the evaluation object of unit damage or unit generalized displacement corresponding relation, be identical with the element of the same numbering of other vectors with the corresponding relation of this evaluation object; d
^{i} _{tk}with evaluation object initial damage vector d
_{o}element number 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 a linear relationship error vector e each time
^{i} _{k}; e
^{i} _{k}subscript k represent kth (k=1,2,3 ..., N) and individual evaluation object increases unit damage or unit generalized displacement.There is N number of evaluation object just to have N calculating, just have N number of linear relationship error vector e
^{i} _{k}, by this N number of linear relationship error vector e
^{i} _{k}obtaining a vector after addition, is exactly final linear relationship error vector e by each element of this vector divided by the new vector obtained after N
^{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.
Tenth step: define current nominal fatigue vector d
^{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 generalized displacement of support of corresponding evaluation object, d
^{i} _{c}and d
^{i}with evaluation object initial damage vector d
_{o}element number rule identical, d
^{i} _{c}element, d
^{i}element and d
_{o}element be onetoone relationship.
11 step: according to monitored amount current value vector C
^{i}with " monitored amount current initial value vector C
^{i} _{o}", " unit damage monitored numerical quantity unit change matrix Δ C
^{i}" and " current nominal fatigue vector d
^{i} _{c}" between exist linear approximate relationship, this linear approximate relationship can be expressed as formula (11), according to multiobjective optimization algorithm calculate current nominal fatigue vector d
^{i} _{c}noninferior solution, namely can determine the position of damaged cable and the solution of nominal fatigue degree thereof more exactly with reasonable error from all ropes.
The multiobjective optimization algorithm that can adopt has a variety of, such as: 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 algorithm is all conventional algorithm, can realize easily, this implementation step only provides for Objective Programming and solves current nominal fatigue vector d
^{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 Ω limits 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).Formula (35) be meant to the minimum real number γ of searching one, formula (36) is met.G (d in formula (36)
^{i} _{c}) defined by formula (37), the middle G (d of the product representation formula (36) of weighing vector W and γ in formula (36)
^{i} _{c}) and vectorial g
^{i}between allow deviation, g
^{i}definition see formula (17), its value calculates in the 9th step.During actual computation vector W can with vectorial g
^{i}identical.The concrete programming realization of Objective Programming has had universal program directly to adopt.Use Objective Programming just can in the hope of current nominal fatigue vector d
^{i} _{c}.
minimize γ
(35)
γ∈R，
${d}_{c}^{i}\∈\mathrm{\Ω}$
12 step: according to evaluation object current actual damage vector d
^{i}definition (see formula (18)) and the definition (see formula (19)) of its element calculate current actual damage vector d
^{i}each element, thus can by d
^{i}determine the health status of evaluation object.Evaluation object current actual damage vector d
^{i}a kth element d
^{i} _{k}represent the current actual health status of a kth evaluation object in ith circulation, if this evaluation object is support cable, so a d in cable system
^{i} _{k}represent its current actual damage, d
^{i} _{k}represent when being 0 that the support cable of its correspondence is without health problem, d
^{i} _{k}numerical value represents when not being 0 that the support cable of its correspondence is the support cable of unsoundness problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, the degree of the lax or damage of its numerical response; If this evaluation object is generalized displacement component, so a d of a bearing
^{i} _{k}represent its current actual generalized displacement numerical value; So according to evaluation object current actual damage vector d
^{i}which support cable unsoundness problem and order of severity thereof can be defined, define which bearing and there occurs generalized displacement and numerical value thereof; So far the generalized displacement of support identification of Cable Structure is achieved.
By current for evaluation object actual damage vector d
^{i}in the M relevant to support cable
_{1}individual element takes out, composition support cable current actual damage vector d
^{ic}, support cable current actual damage vector d
^{ic}the coding rule of element and Initial cable force vector F
_{o}the coding rule of element identical.Support cable current actual damage vector d
^{ic}h element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., M
_{1}; Current actual damage vector d
^{ic}middle numerical value be not 0 element correspond 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 by removing the protective seam of the support cable of unsoundness problem; visual discriminating is carried out to support cable; or carry out visual discriminating by optical imaging apparatus; also can be differentiated whether support cable is impaired 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 that damage, unsoundness problem is exactly there occurs lax rope, and needing the rope adjusting Suo Li exactly, is exactly slack line, and these ropes that need adjust Suo Li are at support cable current actual damage vector d
^{ic}in corresponding element numerical value (such as one of them element can use d
^{ic} _{h}represent) represent the degree of injury with the relax level mechanic equivalent of these support cables, just determine slack line thus, the computing method of concrete slack are described below.According to support cable current actual damage vector d
^{ic}, identify slack line from the support cable of unsoundness problem after, remaining is exactly damaged cable, and these damaged cables are at support cable current actual damage vector d
^{ic}the numerical value of the element of middle correspondence just represents its degree of injury, the numerical value of corresponding element represents when being 100% that this support cable thoroughly loses loadbearing capacity, represent time between 0 and 100% that this support cable loses the loadbearing capacity of corresponding proportion, so far just have identified damaged cable and degree of injury thereof.
13 step: utilize at current cable structure steady temperature data vector T
^{i}the support cable current actual damage vector d obtained in the 12 step under condition
^{ca}that obtain slack line with current actual equivalent damage degree that is its relax level mechanic equivalent, utilize the 6th step obtain at current cable structure steady temperature data vector T
^{i}current cable force vector F under condition and current support cable two support end points horizontal range vector l
^{t} _{x}, utilize 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}with the weight vector ω of initial free unit length
_{o}, utilize current cable structure steady temperature data vector T
^{i}represent support cable current steady state temperature data, utilize second step obtain at initial Cable Structure steady temperature data vector T
_{o}the support cable initial steady state temperature data represented, the temperature variant physical and mechanical properties parameter of the various materials utilizing the Cable Structure obtained at second step to use, count the impact of temperature variation on support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanic equivalent calculate slack line, with the relax level of current actual equivalent damage degree equivalence, the mechanical condition of equivalence is: one, two equivalences rope without lax identical with the mechanics parameters of initial drift during not damaged, geometrical property parameter, density and material; Two, after lax or damage, the slack line of two equivalences and the Suo Li of damage rope be out of shape after overall length identical.When meeting abovementioned two equivalent conditions, the such mechanics function of two support cables in Cable Structure is exactly identical, if after namely replacing damaged cable with the slack line of equivalence, any change can not occur Cable Structure, and vice versa.Try to achieve according to aforementioned mechanic equivalent condition the relax level that those are judged as slack line, relax level is exactly the knots modification of support cable drift, namely determines the long adjustment amount of those ropes that need adjust the support cable of Suo Li.Particularly can in the hope of the relax level of these ropes (i.e. the long adjustment amount of rope) according to formula (29) or formula (30).So just achieve the lax identification of support cable.So far just damaged cable, slack line and generalized displacement of support is all identified.
14 step: the computing machine in health monitoring systems regularly generates cable system health condition form automatically or by human users's health monitoring systems.
15 step: under specified requirements, the computing machine automatic operation communication alert equipment in health monitoring systems is reported to the police to monitor staff, owner and (or) the personnel that specify.
16 step: set up mark vector B according to formula (31)
^{i}, formula (32) gives mark vector B
^{i}the definition of a kth element; If mark vector B
^{i}element be 0 entirely, then get back to the 6th step and proceed health monitoring to cable system and calculating; If mark vector B
^{i}element be not 0, then after completing subsequent step entirely, enter and circulate next time.
17 step: according to formula (33) calculate next time (namely the ith+1 time, i=1,2,3,4 ...) needed for circulation initial damage vector d
^{i+1} _{o}each element d
^{i+1} _{ok}(k=1,2,3 ..., N); The second, at initial mechanical Calculation Basis model A
_{o}basis on, to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, then makes the health status of evaluation object be d
^{i+1} _{o}after obtain be exactly next time, namely the ith+1 time (i=1,2,3,4 ...) Mechanics Calculation benchmark model A needed for circulation
^{i+1}; Next time (namely the ith+1 time, i=1,2,3,4 ...) current initial Cable Structure steady temperature data vector T needed for circulation
^{i+1} _{o}equal T
^{i} _{o}.Obtain A
^{i+1}, d
^{i+1} _{o}and T
^{i+1} _{o}after, obtain A by Mechanics Calculation
^{i+1}in all monitored amounts, current concrete numerical value, these the monitored amounts of concrete numerical value composition next time, namely needed for the ith+1 time circulation current initial value vector C
^{i+1} _{o}.
18 step: get back to the 6th step, starts by the circulation of the 6th step to the 17.
Claims (1)
1. temperature variation angle monitor problem cable and a generalized displacement of support progressivetype recognition method, is characterized in that described method comprises:
A. for sake of convenience, this method unitedly calls evaluated support cable and generalized displacement of support component to be evaluation object, if the quantity sum of the quantity of evaluated support cable and generalized displacement of support component is N, namely the quantity of evaluation object is N; Determine the coding rule of evaluation object, evaluation objects all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step; This method variable k represents this numbering, k=1,2,3 ..., N; Specify when determining hybrid monitoring by the support cable of monitored Suo Li, if total M in cable system
_{1}root support cable, the monitored rope force data of Cable Structure specifies Q rope force data of support cable to describe by Q in Cable Structure, and the change of Cable Structure Suo Li is exactly the change of the Suo Li of all appointment support cables; Each total Q cable force measurement value or calculated value characterize the rope force information of Cable Structure; Q is one and is not less than 0, is not more than M
_{1}integer; Specify when determining hybrid monitoring by the measured point of monitored strain, the monitored strain data of Cable Structure can by K in Cable Structure
_{2}the L of individual specified point and each specified point
_{2}the strain of individual assigned direction describes, and the change of Cable Structure strain data is exactly K
_{2}the change of all tested strain 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}be be not less than 0 integer; Specify when determining hybrid monitoring by the measured point of monitored angle, the monitored angledata of Cable Structure is by K in Cable Structure
_{3}individual specified point, cross the L of each specified point
_{3}the H of individual appointment straight line, each appointment straight line
_{3}individual angle coordinate component describes, and the change of Cable Structure angle is exactly the change of all specified points, all appointment straight line, all angle coordinate component of specifying; 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}it is an integer being not less than 0; Specify when determining hybrid monitoring by monitored shape data, the monitored shape data of Cable Structure is by K in Cable Structure
_{4}the L of individual specified point and each specified point
_{4}the volume coordinate of individual assigned direction describes, and the change of Cable Structure shape data is exactly K
_{4}the change of all coordinate components of individual specified point; Each total M
_{4}individual coordinates measurements or calculated value characterize Cable Structure shape, M
_{4}for K
_{4}and L
_{4}longpending; M
_{4}it is an integer being not less than 0; The monitored amount of comprehensive abovementioned hybrid monitoring, whole Cable Structure has M monitored amount, and M is Q, M
_{2}, M
_{3}and M
_{4}sum, definition parameter K, K is Q, K
_{2}, K
_{3}and K
_{4}sum, K and M must not be less than the quantity N of evaluation object; Because M monitored amount is dissimilar, so this method is called " during temperature variation the damaged cable of hybrid monitoring and support translation identification method "; For simplicity, in the method by monitored all parameters of Cable Structure " during the hybrid monitoring " listed by this step referred to as " monitored amount "; The quantity sum of all monitored amounts is designated as M, and M must not be less than N; Must not be greater than 30 minutes to the time interval between any twice measurement of same amount RealTime Monitoring in this method, the moment of survey record data is called the physical record data moment;
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 environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, asconstructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model of Cable Structure, inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, in the method daytime can not be seen one of the sun and be called the cloudy day all day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, do not represent that the same day necessarily can see the sun, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day
_{r}, be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation
_{h}, get Δ T for convenience of describing
_{h}unit be DEG C/m, the surface of Cable Structure is got " R Cable Structure surface point ", the Specific Principles getting " R Cable Structure surface point " describes in step b3, the temperature of this R Cable Structure surface point will be obtained below by actual measurement, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ", from the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, at the sea level elevation place that each is chosen, two points are at least chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along 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 in the shade the outer normal direction of the measurement Cable Structure chosen along the sunny slope outer normal direction and Cable Structure that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness, three points are no less than along each measurement Cable Structure along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, along each, Cable Structure is measured for support cable and only gets a point along the direction of the Temperature Distribution of wall thickness, only measure the temperature of the surface point of support cable, measure all temperature be selected a little, the temperature recorded 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 the direction of Cable Structure along the Temperature Distribution of wall thickness " to 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 have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", measure temperature in Cable Structure location according to meteorology to require to choose a position, obtain meeting the temperature that meteorology measures the Cable Structure place environment of temperature requirement by the actual measurement of this position, in the onsite spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable, at the flat board of this position of sound production one piece of carbon steel material, 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 with dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, the nonsunny slope of reference plate is covered with insulation material, RealTime Monitoring is obtained 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 the temperature profile data of previously defined Cable Structure along thickness simultaneously, and RealTime Monitoring obtains meeting the temperature record that meteorology measures the Cable Structure place environment of temperature requirement simultaneously, the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be called environment maximum temperature difference, be designated as Δ T
_{emax}, calculated the rate of change of temperature about the time of Cable Structure place environment by Conventional mathematical by the temperature measured data sequence of Cable Structure place environment, this rate of change is also along with time variations, the measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be called reference plate maximum temperature difference, be designated as Δ T
_{pmax}, the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein called Cable Structure surface maximum temperature difference, be designated as Δ T
_{smax}, calculated the rate of change of temperature about the time of each Cable Structure surface point by Conventional mathematical by each Cable Structure surface temperature measured data sequence, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations, obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by RealTime Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and minimum temperature at the sea level elevation place that each is chosen, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T
_{tmax},
B3: survey calculation obtains Cable Structure steady temperature data, first, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day 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 Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the environment maximum temperature difference Δ T that survey calculation obtains above
_{emax}be not more than with reference to temperature difference per day Δ T
_{r}, and reference plate maximum temperature difference Δ T
_{pmax}Δ T is not more than after deducting 2 degrees Celsius
_{emax}, and Cable Structure surface maximum temperature difference Δ T
_{smax}be not more than Δ T
_{pmax}, one of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition, Section 3 condition is when obtaining Cable Structure steady temperature data, and 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, Section 4 condition is when obtaining Cable Structure steady temperature data, and 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, Section 5 condition is when obtaining 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 be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise, Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T
_{tmax}be not more than 1 degree Celsius, this method utilizes abovementioned six conditions, any one in following three kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in abovementioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, the second moment is moment of the Section 6 condition only met in abovementioned " condition relevant to determining to obtain moment of Cable Structure steady temperature data ", simultaneously the third moment meets the moment of the Section 1 in abovementioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition, when the mathematics moment obtaining Cable Structure steady temperature data is exactly a moment in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data, if the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data, the amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method, this method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, 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 obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution of the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steadystate surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steadystate surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steadystate 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 ", when the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%, Cable Structure surface comprises support cable surface, second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", 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 DEG C 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 obtained, gets Δ T for convenience of describing
_{h}unit be DEG C/m, be m for convenience of describing the unit getting Δ h, " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not 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, 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 geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshineduration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshineduration the most fully those surface points in Cable Structure,
C. the Cable Structure steady temperature data under original state are obtained according to " the temperature survey calculating method of the Cable Structure of this method " direct survey calculation, 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}"; Survey or consult reference materials and obtain the temperature variant physical and mechanical properties parameter of the various materials that Cable Structure uses; Initial Cable Structure steady temperature data vector T is obtained in actual measurement
_{o}synchronization, direct survey calculation obtains the Initial cable force of all support cables, composition Initial cable force vector F
_{o}; Obtain the length of all support cables when free state and Suo Li are 0 according to Cable Structure design data, completion data, in free state time crosssectional area and in free state time the weight of unit length, and the temperature of all support cables when obtaining these three kinds of data, utilize temperature variant physical function parameter and the mechanical property parameters of all support cables on this basis, conveniently physical computing obtains all support cables at initial Cable Structure steady temperature data vector T
_{o}the weight of the unit length of all support cables when the crosssectional area of all support cables and Suo Li are 0 when the length of all support cables, Suo Li are 0 when Suo Li under condition is 0, form the initial drift vector of support cable, the weight vector of the initial free unit length of initial free crosssectional area vector sum successively, the coding rule 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 and Initial cable force vector F
_{o}the coding rule of element identical; T is obtained in actual measurement
_{o}while, namely at the initial Cable Structure steady temperature data vector T of acquisition
_{o}the synchronization in moment, direct survey calculation obtains the measured data of initial Cable Structure, and the measured data of initial Cable Structure comprises the Nondestructive Testing Data of the health status expressing support cable, Cable Structure bearing initial generalized displacement measurement data, the initial value of all monitored amounts, the Initial cable 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; 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; Utilization can express the Nondestructive Testing Data of the health status of support cable and the initial generalized displacement measurement data of Cable Structure bearing sets up evaluation object initial damage vector d
_{o}, vectorial d
_{o}represent with initial mechanical Calculation Basis model A
_{o}the initial health of the evaluation object of the Cable Structure represented; 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 support cable corresponding to this element is intact, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to this element completely loses loadbearing capacity, if its numerical value is between 0 and 100%, then represent that this support cable loses the loadbearing capacity of corresponding proportion; If d
_{o}evaluation object corresponding to some elements be some generalized displacement components of some bearings, so d
_{o}the numerical value of this element represent the initial value of this generalized displacement component of this bearing; If there is no the Nondestructive Testing Data of support cable and other are when can express the data of the health status of support cable, or can think Cable Structure original state be not damaged without relaxed state time, vectorial d
_{o}in each element numerical value relevant to support cable get 0, if when there is no Cable Structure bearing initial generalized displacement measurement data or can think that the initial generalized displacement of Cable Structure bearing is 0, vectorial d
_{o}in each element numerical value relevant to Cable Structure generalized displacement of support get 0; Initial Cable Structure bearing generalized coordinate data refer to the support coordinate data under Cable Structure design point, and Cable Structure bearing initial generalized displacement measurement data refers to setting up initial mechanical Calculation Basis model A
_{o}time, the generalized displacement that Cable Structure bearing occurs relative to the bearing under Cable Structure design point; Bearing generalized coordinate comprises line amount and angular amount two kinds;
Temperature variant physical and mechanical properties parameter, the initial Cable Structure steady temperature data vector T of the various materials d. used according to the measured data of the design drawing of Cable Structure, asconstructed drawing and initial Cable Structure, the Nondestructive Testing Data of support cable, Cable Structure bearing initial generalized displacement measurement data, Cable Structure
_{o}with all Cable Structure data that preceding step obtains, set up the initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data "
_{o}, based on A
_{o}the Cable Structure that calculates calculates data must closely its measured data, 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}evaluation object health status with evaluation object initial damage vector d
_{o}represent; Corresponding to A
_{o}the initial value monitored amount initial value vector C of all monitored amount
_{o}represent; T
_{o}and d
_{o}a
_{o}parameter, by A
_{o}the initial value of all monitored amount that obtains of Mechanics Calculation result and C
_{o}the initial value of all monitored amount represented is identical, therefore alternatively C
_{o}by A
_{o}mechanics Calculation result composition, A in the method
_{o}, C
_{o}, d
_{o}and T
_{o}constant;
E. in the method, alphabetical i is except representing the place of number of steps significantly, and alphabetical i only represents cycle index, i.e. ith circulation; Ith circulation needs the current initial mechanical Calculation Basis model of Cable Structure that is that set up or that set up to be designated as current initial mechanical Calculation Basis model A when starting
^{i} _{o}, A
_{o}and A
^{i} _{o}count temperature parameter, can the Effect on Mechanical Properties of accounting temperature change to Cable Structure; When ith 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 one_to_one corresponding; The current initial damage vector of evaluation object that ith circulation needs when starting is designated as d
^{i} _{o}, d
^{i} _{o}cable Structure A when representing that this circulation starts
^{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 one_to_one corresponding; When ith circulation starts, the initial value of all monitored amounts, with monitored amount current initial value vector 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 one_to_one corresponding, monitored amount current initial value vector C
^{i} _{o}represent and correspond to A
^{i} _{o}the concrete numerical value of all monitored amount; T
^{i} _{o}and d
^{i} _{o}a
^{i} _{o}characterisitic parameter, C
^{i} _{o}by A
^{i} _{o}mechanics Calculation result composition; When first time, circulation started, A
^{i} _{o}be designated as A
^{1} _{o}, set up A
^{1} _{o}method for making A
^{1} _{o}equal A
_{o}; When first time, circulation started, T
^{i} _{o}be designated as T
^{1} _{o}, set up T
^{1} _{o}method for making T
^{1} _{o}equal T
_{o}; When first time, circulation started, d
^{i} _{o}be designated as d
^{1} _{o}, set up d
^{1} _{o}method for making d
^{1} _{o}equal d
_{o}; When first time, circulation started, 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 being walked to s step by f here; In Cable Structure military service process, the current data of Cable Structure steady temperature data is obtained, the current data composition current cable structure steady temperature data vector T of all " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement
^{i}, vector T
^{i}definition mode and vector T
_{o}definition mode identical, T
^{i}element and T
_{o}element one_to_one corresponding; Vector T is obtained in actual measurement
^{i}while, actual measurement obtains at acquisition current cable structure steady temperature data vector T
^{i}moment synchronization Cable Structure 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 one_to_one corresponding, represent that identical monitored amount is at not numerical value in the same time; Current cable structure steady temperature data vector T is obtained in actual measurement
^{i}synchronization, actual measurement obtain all M in Cable Structure
_{1}the rope force data of root support cable, all these rope force datas composition current cable force vector F, the element of vectorial F and vectorial F
_{o}the coding rule of element identical; Current cable structure steady temperature data vector T is obtained in actual measurement
^{i}synchronization, Actual measurement obtains all M
_{1}the volume coordinate of two supporting end points of root support cable, the difference of the volume coordinate component in the horizontal direction of two supporting end points is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables form current support cable two and support end points horizontal range vector, and current support cable two supports coding rule and the Initial cable force vector F of the element of end points horizontal range vector
_{o}the coding rule of element identical;
G. according to current cable structure steady temperature data vector T
^{i}, upgrade current initial mechanical Calculation Basis model A according to step g 1 to g3
^{i} _{o}, monitored amount current initial value vector C
^{i} _{o}with current initial Cable Structure steady temperature data vector T
^{i} _{o}, and evaluation object current initial damage vector d
^{i} _{o}remain unchanged;
G1. T is compared
^{i}with T
^{i} _{o}if, T
^{i}equal T
^{i} _{o}, then A
^{i} _{o}, C
^{i} _{o}and T
^{i} _{o}remain unchanged; Otherwise need to follow these steps to A
^{i} _{o}and T
^{i} _{o}upgrade;
G2. T is calculated
^{i}with T
_{o}difference, T
^{i}with T
_{o}difference be exactly the changes of current cable structure steady temperature data about initial Cable Structure steady temperature data, T
^{i}with T
_{o}difference represent with steady temperature change vector S, S equals T
^{i}deduct T
_{o}, S represents the change of Cable Structure steady temperature data;
G3. to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A
_{o}in Cable Structure apply temperature variation after obtain upgrade current initial mechanical Calculation Basis model A
^{i} _{o}, upgrade A
^{i} _{o}while, T
^{i} _{o}all elements numerical value also uses T
^{i}all elements numerical value correspondence replace, namely have updated T
^{i} _{o}, so just obtain and correctly correspond to A
^{i} _{o}t
^{i} _{o}; Now d
^{i} _{o}remain unchanged; As renewal A
^{i} _{o}after, A
^{i} _{o}the health status evaluation object of rope current initial damage vector d
^{i} _{o}represent, A
^{i} _{o}cable Structure steady temperature current cable structure steady temperature data vector T
^{i}represent, upgrade C
^{i} _{o}method be: when renewal A
^{i} _{o}after, obtain A by Mechanics Calculation
^{i} _{o}in all monitored amounts, current concrete numerical value, these concrete numerical value composition C
^{i} _{o};
H. at current initial mechanical Calculation Basis model A
^{i} _{o}basis on, carry out several times Mechanics Calculation according to step h1 to step h4, by calculate set up unit damage monitored numerical quantity unit change matrix Δ C
^{i}with unit damage or unit generalized displacement vector D
^{i} _{u};
H1., when ith circulation starts, directly Δ C is obtained by method listed by step h2 to step h4
^{i}and D
^{i} _{u}; In other moment, when in step g to A
^{i} _{o}after upgrading, Δ C must be regained by method listed by step h2 to step h4
^{i}and D
^{i} _{u}if, not to A in step g
^{i} _{o}upgrade, then directly proceed to step I herein and carry out followup work;
H2. at the current initial mechanical Calculation Basis model A of Cable Structure
^{i} _{o}basis on carry out several times Mechanics Calculation, calculation times numerically equals the quantity N of all evaluation objects, has N number of evaluation object just to have N calculating; According to the coding rule of evaluation object, calculate successively; Calculating hypothesis each time only has an evaluation object to increase unit damage or unit generalized displacement again on the basis of original damage or generalized displacement, 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 the generalized displacement component in a direction of a bearing, just suppose that this bearing increases unit generalized displacement again at this sense of displacement, use D
^{i} _{uk}record unit damage or the unit generalized displacement of this increase, wherein k represents the numbering of the evaluation object increasing unit damage or unit generalized displacement, D
^{i} _{uk}unit damage or unit generalized displacement vector D
^{i} _{u}an element, unit damage or unit generalized displacement vector D
^{i} _{u}the coding rule of element and vectorial d
_{o}the coding rule of element identical; The evaluation object increasing unit damage or unit generalized displacement in calculating each time is again different from during other time calculates the evaluation object increasing unit damage or unit generalized displacement again, calculate the current calculated value all utilizing mechanics method to calculate all monitored amount of Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector; When supposing that a kth evaluation object increases unit damage or unit generalized displacement again, use C
^{i} _{tk}represent corresponding " monitored amount calculation current vector "; When giving the element number of each vector in this step, same coding rule should be used with other vector in this method, to ensure any one element in this step in each vector, with in other vector, number identical element, have 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 one_to_one corresponding;
H3. the vectorial C calculated each time
^{i} _{tk}deduct vectorial C
^{i} _{o}obtain a vector, then after each element of this vector is calculated divided by this unit damage or unit generalized displacement numerical value supposed, obtain " unit change vector δ a C for monitored amount
^{i} _{k}"; N number of evaluation object is had just to have N number of " the unit change vector of monitored amount ";
H4. by this is N number of " the unit change vector of monitored amount " according to the coding rule of N number of evaluation object, " the unit damage monitored numerical quantity unit change matrix Δ C having N to arrange is formed successively
^{i}"; Unit damage monitored numerical quantity unit change matrix Δ C
^{i}each row correspond to a monitored amount unit change vector; Unit damage monitored numerical quantity unit change matrix Δ C
^{i}every a line correspond to the different unit change amplitude of same monitored amount when different evaluation object increases unit damage or unit generalized displacement; Unit damage monitored numerical quantity unit change matrix Δ C
^{i}the coding rule of row and vectorial d
_{o}the coding rule of element identical, unit damage monitored numerical quantity unit change matrix Δ C
^{i}the coding rule of coding rule and M monitored amount of row identical;
I. current nominal fatigue vector d is defined
^{i} _{c}with current actual damage vector d
^{i}, d
^{i} _{c}and d
^{i}element number equal the quantity of support cable, 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 health status of corresponding evaluation object, d
^{i} _{c}and d
^{i}with evaluation object initial damage vector d
_{o}element number 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 " monitored amount current initial value vector C
^{i} _{o}", " unit damage monitored numerical quantity unit change matrix Δ C
^{i}" and " current nominal fatigue vector d
^{i} _{c}" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, except d in formula 1
^{i} _{c}other outer amount is known, solves formula 1 and just can calculate current nominal fatigue vector d
^{i} _{c};
K. the current actual damage vector d utilizing formula 2 to express
^{i}a kth element d
^{i} _{k}with evaluation object current initial damage vector d
^{i} _{o}a kth element d
^{i} _{ok}with current nominal fatigue vector d
^{i} _{c}a kth element d
^{i} _{ck}between relation, calculate current actual damage vector d
^{i}all elements;
k=1 in formula 2 formula 2,2,3 ..., N; d
^{i} _{k}represent the current actual health status of a kth evaluation object in ith circulation, d
^{i} _{k}represent when being 0 that a kth evaluation object is without health problem, d
^{i} _{k}numerical value represents when not being 0 that a kth evaluation object is the evaluation object of unsoundness problem, if this evaluation object is support cable, so a d in cable system
^{i} _{k}represent the order of severity of its current health problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, d
^{i} _{k}the degree of the lax or damage of this support cable of numerical response, if this evaluation object is generalized displacement component, so a d of a bearing
^{i} _{k}represent its current actual generalized displacement numerical value; So according to evaluation object current actual damage vector d
^{i}which support cable unsoundness problem and order of severity thereof can be defined, define which bearing and there occurs generalized displacement and numerical value thereof;
L. by current for evaluation object actual damage vector d
^{i}in with M
_{1}the M that root support cable is relevant
_{1}individual element takes out, composition support cable current actual damage vector d
^{ca}, support cable current actual damage vector d
^{ca}the coding rule of element and Initial cable force vector F
_{o}the coding rule of element identical; Support cable current actual damage vector d
^{ca}h element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., M
_{1}; Support cable current actual damage vector d
^{ca}middle numerical value be not 0 element correspond to the support cable of unsoundness problem, from the support cable of these unsoundness problems, identify damaged cable, remaining is exactly slack line, support cable current actual damage vector d
^{ca}in correspond to slack line element numerical expression be the current actual equivalent damage degree with slack line relax level mechanic equivalent;
M. utilize at current cable structure steady temperature data vector T
^{i}under condition, walk the slack line that identifies and with support cable current actual damage vector d at l
^{ca}these slack lines of expressing, with the current actual equivalent damage degree of its relax level mechanic equivalent, utilize f step obtain at current cable structure steady temperature data vector T
^{i}current cable force vector F under condition and current support cable two support end points horizontal range vector, utilize 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, the vectorial F of Initial cable force
_{o}, utilize current cable structure steady temperature data vector T
^{t}represent support cable current steady state temperature data, utilize c step obtain at initial Cable Structure steady temperature data vector T
_{o}the support cable initial steady state temperature data represented, the temperature variant physical and mechanical properties parameter of the various materials utilizing the Cable Structure obtained in c step to use, count the impact of temperature variation on support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanic equivalent calculate slack line, with the relax level of current actual equivalent damage degree equivalence, mechanic equivalent condition is: one, two equivalences rope without lax identical with the mechanics parameters of initial drift during not damaged, geometrical property parameter, density and material; Two, after lax or damage, the slack line of two equivalences and the Suo Li of damage rope be out of shape after overall length identical; When meeting abovementioned two equivalent conditions, the such mechanics function of two support cables in Cable Structure is exactly identical, if after namely replacing damaged cable with the slack line of equivalence, any change can not occur Cable Structure, and vice versa; Try to achieve according to aforementioned mechanic equivalent condition the relax level that those are judged as slack line, relax level is exactly the knots modification of support cable drift, namely determines the long adjustment amount of those ropes that need adjust the support cable of Suo Li; So just achieve lax identification and the nondestructive tests of support cable; During calculating, institute's demand power is provided by current cable force vector F corresponding element;
N. current nominal fatigue vector d is tried to achieve
^{i} _{c}after, set up mark vector B according to formula 3
^{i}, formula 4 gives mark vector B
^{i}the definition of a kth element;
formula 4
Element B in formula 4
^{i} _{k}mark vector B
^{i}a kth element, D
^{i} _{uk}unit damage or unit generalized displacement vector D
^{i} _{u}a kth element, d
^{i} _{ck}evaluation object current nominal fatigue vector d
^{i} _{c}a kth element, they all represent the relevant information of a kth evaluation object, k=1 in formula 4,2,3 ..., N;
If o. mark vector B
^{i}element be 0 entirely, then get back to step f continue this circulation; If mark vector B
^{i}element be not 0 entirely, then enter next step, i.e. step p;
P. calculate next time according to formula 5, evaluation object current initial damage vector d namely needed for the ith+1 time circulation
^{i+1} _{o}each element;
d in formula 5 formula 5
^{i+1} _{ok}the evaluation object current initial damage vector d next time, namely needed for the ith+1 time circulation
^{i+1} _{o}a kth element, d
^{i} _{ok}this, i.e. the evaluation object of ith circulation current initial damage vector d
^{i} _{o}a kth element, D
^{i} _{uk}unit damage or the unit generalized displacement vector D of ith circulation
^{i} _{u}a kth element, B
^{i} _{k}the mark vector B of ith circulation
^{i}a kth element, k=1 in formula 5,2,3 ..., N;
Q. take off once, namely the ith+1 time circulation needed for current initial Cable Structure steady temperature data vector T
^{i+1} _{o}equal the current initial Cable Structure steady temperature data vector T of ith circulation
^{i} _{o};
R. at initial mechanical Calculation Basis model A
_{o}basis on, to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, then makes the health status of rope be d
^{i+1} _{o}after obtain be exactly next time, namely the ith+1 time circulation needed for Mechanics Calculation benchmark model A
^{i+1}; Obtain A
^{i+1}after, obtain A by Mechanics Calculation
^{i+1}in all monitored amounts, current concrete numerical value, these the monitored amounts of concrete numerical value composition next time, namely needed for the ith+1 time circulation current initial value vector C
^{i+1} _{o};
S. get back to step f, start to circulate next time.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201210171652.2A CN102735460B (en)  20120530  20120530  Progressive generalized displacement identification method for problem cables and support based on temperature change angle monitoring 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201210171652.2A CN102735460B (en)  20120530  20120530  Progressive generalized displacement identification method for problem cables and support based on temperature change angle monitoring 
Publications (2)
Publication Number  Publication Date 

CN102735460A CN102735460A (en)  20121017 
CN102735460B true CN102735460B (en)  20150610 
Family
ID=46991358
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201210171652.2A CN102735460B (en)  20120530  20120530  Progressive generalized displacement identification method for problem cables and support based on temperature change angle monitoring 
Country Status (1)
Country  Link 

CN (1)  CN102735460B (en) 
Families Citing this family (1)
Publication number  Priority date  Publication date  Assignee  Title 

CN103630417A (en) *  20131209  20140312  东南大学  Method for progressively identifying defective cables, concentrated load and linear displacement by virtue of angle monitoring 
Citations (5)
Publication number  Priority date  Publication date  Assignee  Title 

CN101387600A (en) *  20081017  20090318  东南大学  Cable system health monitoring method in cable structure based on mixed monitor 
CN101699248A (en) *  20091105  20100428  东南大学  Progressive health monitoring method for identifying loose carrying cable based on angle monitoring 
CN101806668A (en) *  20100402  20100818  东南大学  Cable structure health monitoring method based on cable tension monitoring 
CN101832874A (en) *  20100402  20100915  东南大学  Progressive cable structure health monitoring method based on angle monitoring 
CN102221478A (en) *  20110531  20111019  东南大学  Progressive health monitoring method for cable system based on hybrid monitoring during generalized displacement of support 
Family Cites Families (2)
Publication number  Priority date  Publication date  Assignee  Title 

CN101819098B (en) *  20100331  20110907  东南大学  Mixed monitoring based progressive method for identifying damaged cables and support displacement 
CN102323074A (en) *  20110531  20120118  东南大学  The method of going forward one by one based on the generalized displacement of hybrid monitoring identification damaged cable slack line bearing 

2012
 20120530 CN CN201210171652.2A patent/CN102735460B/en not_active IP Right Cessation
Patent Citations (5)
Publication number  Priority date  Publication date  Assignee  Title 

CN101387600A (en) *  20081017  20090318  东南大学  Cable system health monitoring method in cable structure based on mixed monitor 
CN101699248A (en) *  20091105  20100428  东南大学  Progressive health monitoring method for identifying loose carrying cable based on angle monitoring 
CN101806668A (en) *  20100402  20100818  东南大学  Cable structure health monitoring method based on cable tension monitoring 
CN101832874A (en) *  20100402  20100915  东南大学  Progressive cable structure health monitoring method based on angle monitoring 
CN102221478A (en) *  20110531  20111019  东南大学  Progressive health monitoring method for cable system based on hybrid monitoring during generalized displacement of support 
Also Published As
Publication number  Publication date 

CN102735460A (en)  20121017 
Similar Documents
Publication  Publication Date  Title 

CN103913328A (en)  Generalized displacement hybrid monitoring damaged cable load progressive identification method  
CN103913342A (en)  Method for progressively recognizing fault cable, load and generalized displacement based on angle monitoring  
CN103616249A (en)  Strainmonitoring damaged cable centralized load generalized displacement recognition method  
CN102706659B (en)  Defective cable and support angular displacement progressive identification method based on angular monitoring of temperature change  
CN103604549A (en)  Linear displacement progressive type identification method for strain monitoring of damaged cable and intensive load  
CN103852331A (en)  Progressive recognition method for problem cable load generalized displacement based on mixing monitoring  
CN103852309A (en)  Progressive recognition method for problem cable load linear displacement based on mixing monitoring  
CN103852304A (en)  Method for recognizing damaged cables, loads and linear displacement based on cable force monitoring  
CN103616222A (en)  Hybridmonitoring linear displacement progressive identification method for defective cable and concentrated loads  
CN102706627A (en)  Damaged cable and support angular displacement identification method on basis of hybrid monitoring during temperature variation  
CN103604642A (en)  Problematic cable/concentrated load progressive recognition method on basis of generalized displacement cable power monitoring  
CN103604634A (en)  Linear displacement progressive type identification method for space coordinate monitoring of damaged cable and intensive load  
CN103604623A (en)  Progressive type damaged rope and concentrated load identification method based on angular displacement strain monitoring  
CN103852291A (en)  Progressive recognition method for problematic cable loads through generalized displacement and space coordinate monitoring  
CN102735459B (en)  The problem cable generalized displacement of support progressivetype recognition method of temperature variation cable force monitoring  
CN103604640A (en)  Method for progressive recognition of concentrated loads of damaged cables based on generalized displacement angle monitoring  
CN103604628A (en)  Problematic cable/concentrated load progressive recognition method on basis of linear displacement angle monitoring  
CN103604638A (en)  Problematic cable/concentrated load/linear displacement recognition method on basis of strain monitoring  
CN103604649A (en)  Problematic cable/concentrated load/linear displacement progressive recognition method on basis of strain monitoring  
CN103604629A (en)  Damaged cable/concentrated load/bracket angular displacement progressive recognition method on basis of strain monitoring  
CN103616231A (en)  Generalized displacement strain monitoring progressive identification method for damaged cable and concentrated loads  
CN103604644A (en)  Problematic cable/concentrated load progressive recognition method on basis of linear displacement strain monitoring  
CN103604630A (en)  Problematic cable/concentrated load/linear displacement recognition method on basis of hybrid monitoring  
CN103616245A (en)  Progressive concentrated load support angular displacement identification method for damaged cable in process of cable force monitoring  
CN103616215A (en)  Defective cable, concentrated load and angular displacement identification method based on space coordinate monitoring 
Legal Events
Date  Code  Title  Description 

C06  Publication  
PB01  Publication  
C10  Entry into substantive examination  
SE01  Entry into force of request for substantive examination  
C14  Grant of patent or utility model  
GR01  Patent grant  
CF01  Termination of patent right due to nonpayment of annual fee  
CF01  Termination of patent right due to nonpayment of annual fee 
Granted publication date: 20150610 Termination date: 20180530 