CN116026511B - Cable force automatic identification method, system, computer and readable storage medium - Google Patents
Cable force automatic identification method, system, computer and readable storage medium Download PDFInfo
- Publication number
- CN116026511B CN116026511B CN202310294935.4A CN202310294935A CN116026511B CN 116026511 B CN116026511 B CN 116026511B CN 202310294935 A CN202310294935 A CN 202310294935A CN 116026511 B CN116026511 B CN 116026511B
- Authority
- CN
- China
- Prior art keywords
- peak
- frequency
- peaks
- frequency difference
- calculating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention provides a cable force automatic identification method, a system, a computer and a readable storage medium, wherein the method comprises the steps of carrying out Fourier transformation on original vibration data to obtain a frequency domain diagram; carrying out peak searching in the frequency domain graph to obtain four peaks before the frequency domain, and judging whether the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio or not; if the frequency difference is smaller than the preset value, calculating a four-peak frequency difference value, and calculating a vibration fundamental frequency according to the frequency four-peak frequency difference value; if not, calculating a trimodal frequency difference value, and calculating a vibration fundamental frequency according to the trimodal frequency difference value; and calculating the cable force of the inhaul cable according to the vibration fundamental frequency and a preset cable force calculation formula. According to the invention, through accurately identifying the vibration fundamental frequency of the inhaul cable and combining with a preset cable force calculation formula, the influence of interference signals on cable force identification can be avoided, and automatic identification and monitoring of cable force can be realized.
Description
Technical Field
The invention belongs to the technical field of bridge cable force identification, and particularly relates to an automatic cable force identification method, an automatic cable force identification system, a computer and a readable storage medium.
Background
In the prior art, the cable force of the bridge cable mainly adopts manual monitoring and manual calculation, the vibration fundamental frequency of the cable cannot be automatically identified, and further the cable force cannot be automatically calculated, so that the bridge is greatly puzzled in the construction, use and later maintenance processes.
Disclosure of Invention
In order to solve the technical problems, the invention provides a cable force automatic identification method, a system, a computer and a readable storage medium, which are used for solving the technical problems in the prior art.
In a first aspect, the present invention provides the following technical solutions, and a method for automatically identifying a cable force, where the method includes:
collecting original vibration data of a inhaul cable, and carrying out Fourier transform on the original vibration data to obtain a frequency domain diagram;
carrying out peak searching in the frequency domain graph to obtain four peaks before the frequency domain, wherein the four peaks before the frequency domain comprise a first peak, a second peak, a third peak and a fourth peak, and judging whether the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio or not;
if the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio, calculating a four-peak frequency difference value between two adjacent peaks in the front four peaks of the frequency domain, and calculating a vibration fundamental frequency according to the four-peak frequency difference value;
if the ratio between the height value of the fourth peak and the height value of the third peak is not smaller than the preset ratio, calculating a three-peak frequency difference value between two adjacent peaks in the first peak, the second peak and the third peak, and calculating a vibration fundamental frequency according to the three-peak frequency difference value;
and calculating the cable force of the inhaul cable according to the vibration fundamental frequency and a preset cable force calculation formula.
Compared with the prior art, the beneficial effects of this application are: according to the method, original vibration data of the inhaul cable are collected, and Fourier transformation is conducted on the original vibration data to obtain a frequency domain diagram; and then searching peaks in the frequency domain graph to obtain four peaks before the frequency domain, wherein the four peaks before the frequency domain comprises a first peak, a second peak, a third peak and a fourth peak, judging whether the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio, sequentially converting the original vibration data of the bridge cable into the frequency domain graph so as to identify the vibration fundamental frequency based on the frequency domain graph, calculating the four peak frequency difference between two adjacent peaks in the four peaks before the frequency domain if the ratio between the height value of the fourth peak and the height value of the third peak is smaller than the preset ratio, and calculating the vibration fundamental frequency according to the four peak frequency difference, if the ratio between the height value of the fourth peak and the height value of the third peak is not smaller than the preset ratio, calculating the frequency difference value of three peaks between the first peak, the second peak and two adjacent peaks in the third peak, calculating the vibration fundamental frequency according to the frequency difference value of three peaks, and calculating the cable force of the cable according to the vibration fundamental frequency and a preset cable force calculation formula.
Preferably, the step of performing peak search in the frequency domain map to obtain the first four peaks of the frequency domain includes:
carrying out peak searching in the frequency domain graph to obtain initial first four peaks;
judging whether the frequency interval between two adjacent peaks in the initial front four peaks is larger than a frequency interval threshold value or not;
if the frequency interval between two adjacent peaks in the initial front four peaks is larger than the frequency interval threshold, the frequency domain front four peaks are the initial front four peaks;
if the frequency interval between two adjacent peaks in the initial front four peaks is not greater than the frequency interval threshold, eliminating the peak with the smallest height value in the initial front four peaks, searching for the next peak, re-determining the initial front four peaks according to the next peak, and returning to execute the step of judging whether the frequency interval between two adjacent peaks in the initial front four peaks is greater than the frequency interval threshold so as to obtain the frequency domain front four peaks.
Preferably, the step of calculating a four-peak frequency difference between two adjacent peaks in the four peaks before the frequency domain, and calculating the vibration fundamental frequency according to the four-peak frequency difference includes:
calculating four peak frequency difference values between two adjacent peaks in the front four peaks of the frequency domain, wherein the four peak frequency difference values comprise a first four peak frequency difference value delta f1, a second four peak frequency difference value delta f2 and a third four peak frequency difference value delta f3;
calculating a first fourth peak frequency difference ratio N1 and a second fourth peak frequency difference ratio N2 according to the first fourth peak frequency difference value Δf1, the second fourth peak frequency difference value Δf2, and the third fourth peak frequency difference value Δf3:
and determining a first interval range of the first four-peak frequency difference ratio N1 and the second four-peak frequency difference ratio N2, and calculating a vibration fundamental frequency based on the first interval range.
Preferably, the step of determining a first interval range of the first four-peak frequency difference ratio N1 and the second four-peak frequency difference ratio N2, and calculating the vibration fundamental frequency based on the first interval range includes:
if N1 < M, then:
if M is less than or equal to N1 and less than M+N and N2 is less than M, then:
if M is less than or equal to N1 and less than or equal to M+N and M is less than or equal to N2 and less than or equal to M+N, then:
if M is less than or equal to N1 and less than M+2N and N2 is less than M, then:
if M is less than or equal to N1 and less than or equal to M+2N and M is less than or equal to N2 and less than or equal to M+N, then:
if M is less than or equal to N1 and less than M+2N, and M is less than or equal to N2 and less than M+2N:
wherein f j For the vibration fundamental frequency, M is a first preset value, and N is a second preset value.
Preferably, the step of calculating a trimodal frequency difference between two adjacent peaks among the first peak, the second peak and the third peak, and calculating the vibration fundamental frequency according to the trimodal frequency difference includes:
calculating a three-peak frequency difference value between two adjacent peaks in the first peak, the second peak and the third peak, wherein the three-peak frequency difference value comprises a first three-peak frequency difference value delta f1', and a second three-peak frequency difference value delta f2';
calculating a first trimodal frequency difference ratio N ' from the first trimodal frequency difference value Δf1' and the second trimodal frequency difference value Δf2 ':
and determining a second interval range of the first third peak frequency difference ratio N', and calculating the vibration fundamental frequency based on the second interval range.
Preferably, the step of determining a second interval range of the first peak-to-frequency difference value N' and calculating the vibration fundamental frequency based on the second interval range includes:
if N' < M, then:
if M is less than or equal to N' < M+N, then:
if M is less than or equal to N' < M+2N, then:
wherein f j For the vibration fundamental frequency, M is a first preset value, and N is a second preset value.
Preferably, in the step of collecting the original vibration data of the inhaul cable, performing fourier transform on the original vibration data to obtain a frequency domain map, performing fourier transform on the original vibration data, and converting a time domain of the original vibration data into a frequency domain to generate the frequency domain map.
In a second aspect, the present invention provides a system for automatically identifying a cable force, the system comprising:
the conversion module is used for collecting original vibration data of the inhaul cable and carrying out Fourier conversion on the original vibration data to obtain a frequency domain diagram;
the searching module is used for carrying out peak searching in the frequency domain graph to obtain four peaks before the frequency domain, wherein the four peaks before the frequency domain comprise a first peak, a second peak, a third peak and a fourth peak, and judging whether the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio or not;
the first calculation module is used for calculating a four-peak frequency difference value between two adjacent peaks in the four peaks before the frequency domain if the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio, and calculating a vibration fundamental frequency according to the four-peak frequency difference value;
the second calculation module is used for calculating a three-peak frequency difference value between two adjacent peaks in the first peak, the second peak and the third peak if the ratio between the height value of the fourth peak and the height value of the third peak is not smaller than a preset ratio, and calculating a vibration fundamental frequency according to the three-peak frequency difference value;
and the cable force calculation module is used for calculating the cable force of the inhaul cable according to the vibration fundamental frequency and a preset cable force calculation formula.
In a third aspect, the present invention provides a computer, including a memory, a processor, and a computer program stored in the memory and capable of running on the processor, where the processor implements the method for automatically identifying cable force as described above when executing the computer program.
In a fourth aspect, the present invention provides a readable storage medium, where a computer program is stored, where the computer program is executed by a processor to implement the above-mentioned method for automatically identifying a cable force.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are needed in the embodiments or the description of the prior art will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flowchart of a method for automatically identifying cable force according to a first embodiment of the present invention;
fig. 2 is a detailed flowchart of step S2 in the automatic cable force identification method according to the first embodiment of the present invention;
FIG. 3 is a detailed flowchart of step S3 in the automatic cable force identification method according to the first embodiment of the present invention;
fig. 4 is a detailed flowchart of step S4 in the automatic cable force identification method according to the first embodiment of the present invention;
FIG. 5 is a block diagram of a cable force automatic identification system according to a second embodiment of the present invention;
fig. 6 is a block diagram of a hardware structure of a computer according to another embodiment of the present invention.
Embodiments of the present invention will be further described below with reference to the accompanying drawings.
Detailed Description
Embodiments of the present invention are described in detail below, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to like or similar elements or elements having like or similar functions throughout. The embodiments described below by referring to the drawings are exemplary and intended to illustrate embodiments of the invention and should not be construed as limiting the invention.
In the description of the embodiments of the present invention, it should be understood that the terms "length," "width," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," "outer," and the like indicate orientations or positional relationships based on the orientation or positional relationships shown in the drawings, merely to facilitate description of the embodiments of the present invention and simplify description, and do not indicate or imply that the devices or elements referred to must have a specific orientation, be configured and operated in a specific orientation, and thus should not be construed as limiting the present invention.
Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include one or more such feature. In the description of the embodiments of the present invention, the meaning of "plurality" is two or more, unless explicitly defined otherwise.
In the embodiments of the present invention, unless explicitly specified and limited otherwise, the terms "mounted," "connected," "secured" and the like are to be construed broadly and include, for example, either permanently connected, removably connected, or integrally formed; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communicated with the inside of two elements or the interaction relationship of the two elements. The specific meaning of the above terms in the embodiments of the present invention will be understood by those of ordinary skill in the art according to specific circumstances.
Example 1
As shown in fig. 1, in a first embodiment of the present invention, the present invention provides a method for automatically identifying a cable force, the method comprising:
s1, acquiring original vibration data of a inhaul cable, and performing Fourier transform on the original vibration data to obtain a frequency domain diagram;
specifically, in the step, vibration data of the inhaul cable can be monitored through monitors arranged on the bridge and the inhaul cable, the vibration data is subjected to data preprocessing, vibration data exceeding a vibration threshold range in the vibration data is removed, and the vibration data after the data preprocessing is integrated to obtain original vibration data of the inhaul cable;
specifically, in step S1, fourier transform is performed on the original vibration data, the time domain of the original vibration data is converted into the frequency domain, so as to generate a frequency domain map, and after the time domain of the original vibration data is converted into the frequency domain, the frequency domain map is automatically generated.
S2, carrying out peak search in the frequency domain graph to obtain four peaks before the frequency domain, wherein the four peaks before the frequency domain comprise a first peak, a second peak, a third peak and a fourth peak, and judging whether the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio;
specifically, after the frequency domain map is generated, since the original vibration data has volatility in the frequency domain map, a plurality of peaks are generated in the frequency domain map, in this step, the first four peaks in the frequency domain map are obtained by searching the peaks in the frequency domain map, and are sequentially marked as a first peak, a second peak, a third peak and a fourth peak, and then, whether the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio is determined, so as to determine the vibration fundamental frequency.
As shown in fig. 2, the step S2 includes:
s21, carrying out peak searching in the frequency domain graph to obtain initial front four peaks;
specifically, the first four peaks are the first four peaks in the frequency domain diagram, but there may be peaks with a frequency spacing smaller than the frequency spacing threshold in the first four peaks, so in order to ensure the accuracy of data, it is necessary to determine the frequency spacing between two adjacent peaks in the first four peaks.
S22, judging whether the frequency interval between two adjacent peaks in the initial front four peaks is larger than a frequency interval threshold value.
S23, if the frequency interval between two adjacent peaks in the initial front four peaks is larger than a frequency interval threshold, the frequency domain front four peaks are the initial front four peaks;
specifically, if the frequency interval between two adjacent peaks in the four peaks before the initial peak is greater than the frequency interval threshold, the selected four peaks before the initial peak meets the requirement and has a higher true value, so that the four peaks before the initial peak can be used as the four peaks before the frequency domain.
S24, if the frequency interval between two adjacent peaks in the initial front four peaks is not larger than the frequency interval threshold, eliminating the peak with the smallest height value in the initial front four peaks, searching for the next peak, re-determining the initial front four peaks according to the next peak, and returning to the step of judging whether the frequency interval between two adjacent peaks in the initial front four peaks is larger than the frequency interval threshold or not so as to obtain the front four peaks of the frequency domain;
specifically, if the frequency spacing between two adjacent peaks in the initial previous four peaks is not greater than the frequency spacing threshold, the authenticity of the data cannot be accurately represented, so that the peak with the smallest height value needs to be removed and the step of searching for the next peak is performed, and after searching for the next peak, the step S22 is continuously performed on the new initial previous four peaks to obtain the final previous four peaks in the frequency domain;
it is worth noting that the range of the frequency spacing threshold is [0.1,1].
S3, if the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio, calculating a four-peak frequency difference value between two adjacent peaks in the front four peaks of the frequency domain, and calculating a vibration fundamental frequency according to the four-peak frequency difference value;
specifically, in this embodiment, the preset ratio is 80%, if the ratio between the height value of the fourth peak and the height value of the third peak is smaller than the preset ratio, four peak frequency differences between two adjacent peaks in the four peaks before the frequency domain need to be calculated, and before calculating the four peak frequency differences, four peaks in the four peaks before the frequency domain need to be sequentially arranged from large to small and four peak frequency differences need to be calculated, and because four peaks are included in the four peaks before the frequency domain, the corresponding four peak frequency differences include a first four peak frequency difference Δf1, a second four peak frequency difference Δf2, and a third four peak frequency difference Δf3;
as shown in fig. 3, the step S3 includes:
s31, calculating four peak frequency difference values between two adjacent peaks in the front four peaks of the frequency domain, wherein the four peak frequency difference values comprise a first four peak frequency difference value delta f1, a second four peak frequency difference value delta f2 and a third four peak frequency difference value delta f3.
S32, calculating a first four-peak frequency difference ratio N1 and a second four-peak frequency difference ratio N2 according to the first four-peak frequency difference value Deltaf 1, the second four-peak frequency difference value Deltaf 2 and the third four-peak frequency difference value Deltaf 3:
s33, determining a first interval range of the first four-peak frequency difference ratio N1 and the second four-peak frequency difference ratio N2, and calculating a vibration fundamental frequency based on the first interval range;
specifically, in step S33, if N1 < M, then:
if M is less than or equal to N1 and less than M+N and N2 is less than M, then:
if M is less than or equal to N1 and less than or equal to M+N and M is less than or equal to N2 and less than or equal to M+N, then:
if M is less than or equal to N1 and less than M+2N and N2 is less than M, then:
if M is less than or equal to N1 and less than or equal to M+2N and M is less than or equal to N2 and less than or equal to M+N, then:
if M is less than or equal to N1 and less than M+2N, and M is less than or equal to N2 and less than M+2N:
wherein f j The vibration fundamental frequency is M, wherein M is a first preset value, and N is a second preset value;
wherein M is 1.5, N is 1, and the first interval range of the first peak-to-peak frequency difference ratio N1 and the second peak-to-peak frequency difference ratio N2 is not limited to the above-mentioned determined range, and the corresponding vibration fundamental frequency can be calculated respectively according to the above-mentioned provided first interval range and so on.
S4, if the ratio between the height value of the fourth peak and the height value of the third peak is not smaller than a preset ratio, calculating a three-peak frequency difference value between two adjacent peaks in the first peak, the second peak and the third peak, and calculating a vibration fundamental frequency according to the three-peak frequency difference value;
specifically, in this embodiment, the preset ratio is 80%, if the ratio between the height value of the fourth peak and the height value of the third peak is not smaller than the preset ratio, then the three peak frequency differences between the adjacent two peaks in the first peak, the second peak and the third peak of the four peaks before the frequency domain need to be calculated, and before the three peak frequency differences are calculated, the first peak, the second peak and the third peak need to be arranged in sequence from large to small and three peak frequency differences need to be calculated, and the first peak, the second peak and the third peak have three peaks in total, so the corresponding three peak frequency differences include a first three peak frequency difference Δf1', a second three peak frequency difference Δf2';
as shown in fig. 4, the step S4 includes:
s41, calculating three peak frequency differences among the first peak, the second peak and two adjacent peaks in the third peak, wherein the three peak frequency differences comprise a first three peak frequency difference delta f1', and a second three peak frequency difference delta f2';
s42, calculating a first trimodal frequency difference value N ' according to the first trimodal frequency difference value delta f1' and the second trimodal frequency difference value delta f2 ':
s43, determining a second interval range of the first trimodal frequency difference ratio N', and calculating a vibration fundamental frequency based on the second interval range;
specifically, in the step S43, if N' < M:
if M is less than or equal to N' < M+N, then:
if M is less than or equal to N' < M+2N, then:
wherein f j The vibration fundamental frequency is M, wherein M is a first preset value, and N is a second preset value;
wherein M is 1.5, N is 1, and the second interval range of the first peak-to-peak frequency difference ratio N' is not limited to the above-determined range, and the corresponding vibration fundamental frequency can be calculated according to the above-provided second interval range and the like.
S5, calculating the cable force of the inhaul cable according to the vibration fundamental frequency and a preset cable force calculation formula;
specifically, in steps S3 and S4, the fundamental vibration frequency has been determined, and thus, the fundamental vibration frequency is brought into a preset cable force calculation formula, so that the corresponding cable force of the cable can be calculated, and the preset cable force calculation formula is as follows:
wherein m is the cable linear density, F is the cable force, and l is the cable length;
through the steps, the vibration fundamental frequency can be automatically identified, and the corresponding inhaul cable force is calculated according to the vibration fundamental frequency, so that the automatic monitoring and identification of the inhaul cable force are realized, and meanwhile, the efficiency and the accuracy of inhaul cable force identification are greatly improved.
The first advantage of this embodiment is: according to the method, original vibration data of the inhaul cable are collected, and Fourier transformation is conducted on the original vibration data to obtain a frequency domain diagram; and then searching peaks in the frequency domain graph to obtain four peaks before the frequency domain, wherein the four peaks before the frequency domain comprises a first peak, a second peak, a third peak and a fourth peak, judging whether the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio, sequentially converting the original vibration data of the bridge cable into the frequency domain graph so as to identify the vibration fundamental frequency based on the frequency domain graph, calculating the four peak frequency difference between two adjacent peaks in the four peaks before the frequency domain if the ratio between the height value of the fourth peak and the height value of the third peak is smaller than the preset ratio, and calculating the vibration fundamental frequency according to the four peak frequency difference, if the ratio between the height value of the fourth peak and the height value of the third peak is not smaller than the preset ratio, calculating the frequency difference value of three peaks between the first peak, the second peak and two adjacent peaks in the third peak, calculating the vibration fundamental frequency according to the frequency difference value of three peaks, and calculating the cable force of the cable according to the vibration fundamental frequency and a preset cable force calculation formula.
Example two
As shown in fig. 5, in a second embodiment of the present invention, there is provided an automatic cable force recognition system including:
the transformation module 1 is used for collecting original vibration data of the inhaul cable, and carrying out Fourier transformation on the original vibration data to obtain a frequency domain diagram;
the searching module 2 is configured to perform peak searching in the frequency domain graph to obtain four peaks before the frequency domain, where the four peaks before the frequency domain include a first peak, a second peak, a third peak and a fourth peak, and determine whether a ratio between a height value of the fourth peak and a height value of the third peak is smaller than a preset ratio;
the first calculating module 3 is configured to calculate a four-peak frequency difference value between two adjacent peaks in the front four peaks of the frequency domain if a ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio, and calculate a vibration fundamental frequency according to the four-peak frequency difference value;
the second calculating module 4 is configured to calculate a three-peak frequency difference value between two adjacent peaks among the first peak, the second peak and the third peak if the ratio between the height value of the fourth peak and the height value of the third peak is not less than a preset ratio, and calculate a vibration fundamental frequency according to the three-peak frequency difference value;
and the cable force calculation module 5 is used for calculating the cable force of the inhaul cable according to the vibration fundamental frequency and a preset cable force calculation formula.
Wherein, the transformation module 1 is specifically configured to:
and carrying out Fourier transform on the original vibration data, and converting the time domain of the original vibration data into a frequency domain to generate a frequency domain map.
The search module 2 includes:
the searching sub-module is used for carrying out peak searching in the frequency domain diagram so as to obtain initial first four peaks;
the judging submodule is used for judging whether the frequency interval between two adjacent peaks in the initial front four peaks is larger than a frequency interval threshold value or not;
the first frequency domain peak determining sub-module is used for determining the frequency domain front four peaks as the initial front four peaks if the frequency interval between two adjacent peaks in the initial front four peaks is larger than a frequency interval threshold value;
and the second frequency domain peak determining sub-module is used for eliminating the peak with the smallest height value in the initial front four peaks if the frequency distance between two adjacent peaks in the initial front four peaks is not larger than the frequency distance threshold value, searching the next peak, re-determining the initial front four peaks according to the next peak and returning to the step of executing the judgment on whether the frequency distance between two adjacent peaks in the initial front four peaks is larger than the frequency distance threshold value so as to obtain the frequency domain front four peaks.
The first computing module 3 includes:
the first calculation sub-module is used for calculating four peak frequency difference values between two adjacent peaks in the four peaks before the frequency domain, wherein the four peak frequency difference values comprise a first four peak frequency difference value delta f1, a second four peak frequency difference value delta f2 and a third four peak frequency difference value delta f3;
the second calculating sub-module is configured to calculate a first fourth peak frequency difference ratio N1 and a second fourth peak frequency difference ratio N2 according to the first fourth peak frequency difference value Δf1, the second fourth peak frequency difference value Δf2, and the third fourth peak frequency difference value Δf3:
the first determining submodule is used for determining a first interval range of the first four-peak frequency difference ratio N1 and the second four-peak frequency difference ratio N2 and calculating vibration fundamental frequency based on the first interval range.
The second calculation module 4 includes:
a third calculation sub-module, configured to calculate a trimodal frequency difference value between two adjacent peaks among the first peak, the second peak, and the third peak, where the trimodal frequency difference value includes a first trimodal frequency difference value Δf1', and a second trimodal frequency difference value Δf2';
a fourth calculation sub-module, configured to calculate a first trimodal frequency difference ratio N ' according to the first trimodal frequency difference value Δf1' and the second trimodal frequency difference value Δf2 ':
and the second determining submodule is used for determining a second interval range of the first peak-to-frequency difference value N' and calculating the vibration fundamental frequency based on the second interval range.
In other embodiments of the present invention, a computer is provided in the embodiments of the present invention, and the computer includes a memory 102, a processor 101, and a computer program stored in the memory 102 and capable of running on the processor 101, where the processor 101 implements the cable force automatic identification method described above when executing the computer program.
In particular, the processor 101 may include a Central Processing Unit (CPU), or an application specific integrated circuit (Application Specific Integrated Circuit, abbreviated as ASIC), or may be configured to implement one or more integrated circuits of embodiments of the present application.
The processor 101 reads and executes the computer program instructions stored in the memory 102 to implement the above-described cable force automatic identification method.
In some of these embodiments, the computer may also include a communication interface 103 and a bus 100. As shown in fig. 6, the processor 101, the memory 102, and the communication interface 103 are connected to each other by the bus 100 and perform communication with each other.
The communication interface 103 is used to implement communication between modules, devices, units, and/or units in the embodiments of the present application. The communication interface 103 may also enable communication with other components such as: and the external equipment, the image/data acquisition equipment, the database, the external storage, the image/data processing workstation and the like are used for data communication.
The computer can execute the automatic cable force recognition method based on the obtained cable force automatic recognition system, so that the automatic recognition of the cable force is realized.
In still other embodiments of the present invention, in combination with the above-described cable force automatic identification method, embodiments of the present invention provide a technical solution, a readable storage medium storing a computer program thereon, where the computer program implements the above-described cable force automatic identification method when executed by a processor.
Those of skill in the art will appreciate that the logic and/or steps represented in the flow diagrams or otherwise described herein, e.g., a ordered listing of executable instructions for implementing logical functions, can be embodied in any computer-readable medium for use by or in connection with an instruction execution system, apparatus, or device, such as a computer-based system, processor-containing system, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions. For the purposes of this description, a "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device.
More specific examples (a non-exhaustive list) of the readable medium would include the following: an electrical connection (electronic device) having one or more wires, a portable computer diskette (magnetic device), a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). In addition, the computer readable medium may even be paper or other suitable medium on which the program is printed, as the program may be electronically captured, via, for instance, optical scanning of the paper or other medium, then compiled, interpreted or otherwise processed in a suitable manner, if necessary, and then stored in a computer memory.
It is to be understood that portions of the present invention may be implemented in hardware, software, firmware, or a combination thereof. In the above-described embodiments, the various steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, may be implemented using any one or combination of the following techniques, as is well known in the art: discrete logic circuits having logic gates for implementing logic functions on data signals, application specific integrated circuits having suitable combinational logic gates, programmable Gate Arrays (PGAs), field Programmable Gate Arrays (FPGAs), and the like.
The technical features of the above-described embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above-described embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The above examples merely represent a few embodiments of the present application, which are described in more detail and are not to be construed as limiting the scope of the invention. It should be noted that it would be apparent to those skilled in the art that various modifications and improvements could be made without departing from the spirit of the present application, which would be within the scope of the present application. Accordingly, the scope of protection of the present application is to be determined by the claims appended hereto.
Claims (5)
1. An automatic cable force identification method, which is characterized by comprising the following steps:
collecting original vibration data of a inhaul cable, and carrying out Fourier transform on the original vibration data to obtain a frequency domain diagram;
carrying out peak searching in the frequency domain graph to obtain four peaks before the frequency domain, wherein the four peaks before the frequency domain comprise a first peak, a second peak, a third peak and a fourth peak, and judging whether the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio or not;
if the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio, calculating a four-peak frequency difference value between two adjacent peaks in the front four peaks of the frequency domain, and calculating a vibration fundamental frequency according to the four-peak frequency difference value;
if the ratio between the height value of the fourth peak and the height value of the third peak is not smaller than the preset ratio, calculating a three-peak frequency difference value between two adjacent peaks in the first peak, the second peak and the third peak, and calculating a vibration fundamental frequency according to the three-peak frequency difference value;
calculating the cable force of the inhaul cable according to the vibration fundamental frequency and a preset cable force calculation formula;
the step of performing peak search in the frequency domain map to obtain the first four peaks of the frequency domain comprises the following steps:
carrying out peak searching in the frequency domain graph to obtain initial first four peaks;
judging whether the frequency interval between two adjacent peaks in the initial front four peaks is larger than a frequency interval threshold value or not;
if the frequency interval between two adjacent peaks in the initial front four peaks is larger than the frequency interval threshold, the frequency domain front four peaks are the initial front four peaks;
if the frequency interval between two adjacent peaks in the initial front four peaks is not greater than the frequency interval threshold, eliminating the peak with the smallest height value in the initial front four peaks, searching for the next peak, re-determining the initial front four peaks according to the next peak, and returning to execute the step of judging whether the frequency interval between two adjacent peaks in the initial front four peaks is greater than the frequency interval threshold so as to obtain the frequency domain front four peaks;
the step of calculating the four peak frequency difference between two adjacent peaks in the four peaks before the frequency domain and calculating the vibration fundamental frequency according to the four peak frequency difference comprises the following steps:
calculating four peak frequency difference values between two adjacent peaks in the front four peaks of the frequency domain, wherein the four peak frequency difference values comprise a first four peak frequency difference value delta f1, a second four peak frequency difference value delta f2 and a third four peak frequency difference value delta f3;
calculating a first fourth peak frequency difference ratio N1 and a second fourth peak frequency difference ratio N2 according to the first fourth peak frequency difference value Δf1, the second fourth peak frequency difference value Δf2, and the third fourth peak frequency difference value Δf3:
determining a first interval range of the first four-peak frequency difference ratio N1 and the second four-peak frequency difference ratio N2, and calculating a vibration fundamental frequency based on the first interval range;
the step of determining a first interval range of the first and second peak-to-peak frequency difference ratios N1 and N2, and calculating the vibration fundamental frequency based on the first interval range includes:
if N1 < M, then:
if M is less than or equal to N1 and less than M+N and N2 is less than M, then:
if M is less than or equal to N1 and less than or equal to M+N and M is less than or equal to N2 and less than or equal to M+N, then:
if M+N is less than or equal to N1 and less than M+2N, and N2 is less than M:
if M+N is less than or equal to N1 and less than M+2N, and M is less than or equal to N2 and less than M+N, then:
if M+N is less than or equal to N1 and less than M+2N, and M+N is less than or equal to N2 and less than M+2N:
wherein f j The vibration fundamental frequency is M, wherein M is a first preset value, and N is a second preset value;
the step of calculating the frequency difference of the three peaks between two adjacent peaks in the first peak, the second peak and the third peak, and calculating the vibration fundamental frequency according to the frequency difference of the three peaks comprises the following steps:
calculating a three-peak frequency difference value between two adjacent peaks in the first peak, the second peak and the third peak, wherein the three-peak frequency difference value comprises a first three-peak frequency difference value delta f1', and a second three-peak frequency difference value delta f2';
calculating a first trimodal frequency difference ratio N ' from the first trimodal frequency difference value Δf1' and the second trimodal frequency difference value Δf2 ':
determining a second interval range of the first peak-to-frequency difference ratio N', and calculating a vibration fundamental frequency based on the second interval range;
the step of determining a second interval range of the first trimodal frequency difference ratio value N' and calculating a vibration fundamental frequency based on the second interval range includes:
if N' < M, then:
if M is less than or equal to N' < M+N, then:
if M+N is less than or equal to N' < M+2N, then:
wherein f j For the vibration fundamental frequency, M is a first preset value, N is a second preset valueThe value M is 1.5 and N is 1.
2. The automatic cable force identification method according to claim 1, wherein in the step of acquiring original vibration data of the cable, fourier transforming the original vibration data to obtain a frequency domain map, fourier transforming the original vibration data, and converting a time domain of the original vibration data into a frequency domain to generate the frequency domain map.
3. A cable force automatic identification system applied to the cable force automatic identification method according to claim 1, characterized in that the system comprises:
the conversion module is used for collecting original vibration data of the inhaul cable and carrying out Fourier conversion on the original vibration data to obtain a frequency domain diagram;
the searching module is used for carrying out peak searching in the frequency domain graph to obtain four peaks before the frequency domain, wherein the four peaks before the frequency domain comprise a first peak, a second peak, a third peak and a fourth peak, and judging whether the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio or not;
the first calculation module is used for calculating a four-peak frequency difference value between two adjacent peaks in the four peaks before the frequency domain if the ratio between the height value of the fourth peak and the height value of the third peak is smaller than a preset ratio, and calculating a vibration fundamental frequency according to the four-peak frequency difference value;
the second calculation module is used for calculating a three-peak frequency difference value between two adjacent peaks in the first peak, the second peak and the third peak if the ratio between the height value of the fourth peak and the height value of the third peak is not smaller than a preset ratio, and calculating a vibration fundamental frequency according to the three-peak frequency difference value;
the cable force calculation module is used for calculating the cable force of the inhaul cable according to the vibration fundamental frequency and a preset cable force calculation formula;
the search module comprises:
the searching sub-module is used for carrying out peak searching in the frequency domain diagram so as to obtain initial first four peaks;
the judging submodule is used for judging whether the frequency interval between two adjacent peaks in the initial front four peaks is larger than a frequency interval threshold value or not;
the first frequency domain peak determining sub-module is used for determining the frequency domain front four peaks as the initial front four peaks if the frequency interval between two adjacent peaks in the initial front four peaks is larger than a frequency interval threshold value;
the second frequency domain peak determining sub-module is configured to reject a peak with a smallest height value in the initial preceding four peaks if the frequency spacing between two adjacent peaks in the initial preceding four peaks is not greater than a frequency spacing threshold, search for a next peak, re-determine the initial preceding four peaks according to the next peak, and return to execute the step of determining whether the frequency spacing between two adjacent peaks in the initial preceding four peaks is greater than the frequency spacing threshold, so as to obtain a frequency domain preceding four peaks;
the first computing module includes:
the first calculation sub-module is used for calculating four peak frequency difference values between two adjacent peaks in the four peaks before the frequency domain, wherein the four peak frequency difference values comprise a first four peak frequency difference value delta f1, a second four peak frequency difference value delta f2 and a third four peak frequency difference value delta f3;
the second calculating sub-module is configured to calculate a first fourth peak frequency difference ratio N1 and a second fourth peak frequency difference ratio N2 according to the first fourth peak frequency difference value Δf1, the second fourth peak frequency difference value Δf2, and the third fourth peak frequency difference value Δf3:
the first determining submodule is used for determining a first interval range of the first four-peak frequency difference ratio N1 and the second four-peak frequency difference ratio N2 and calculating vibration fundamental frequency based on the first interval range;
the second computing module includes:
a third calculation sub-module, configured to calculate a trimodal frequency difference value between two adjacent peaks among the first peak, the second peak, and the third peak, where the trimodal frequency difference value includes a first trimodal frequency difference value Δf1', and a second trimodal frequency difference value Δf2';
a fourth calculation sub-module, configured to calculate a first trimodal frequency difference ratio N ' according to the first trimodal frequency difference value Δf1' and the second trimodal frequency difference value Δf2 ':
and the second determining submodule is used for determining a second interval range of the first peak-to-frequency difference value N' and calculating the vibration fundamental frequency based on the second interval range.
4. A computer comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the cable force automatic identification method as claimed in any one of claims 1 to 2 when executing the computer program.
5. A readable storage medium, characterized in that the readable storage medium has stored thereon a computer program which, when executed by a processor, implements the cable force automatic identification method according to any one of claims 1 to 2.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310294935.4A CN116026511B (en) | 2023-03-24 | 2023-03-24 | Cable force automatic identification method, system, computer and readable storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310294935.4A CN116026511B (en) | 2023-03-24 | 2023-03-24 | Cable force automatic identification method, system, computer and readable storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116026511A CN116026511A (en) | 2023-04-28 |
CN116026511B true CN116026511B (en) | 2023-07-04 |
Family
ID=86072734
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310294935.4A Active CN116026511B (en) | 2023-03-24 | 2023-03-24 | Cable force automatic identification method, system, computer and readable storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116026511B (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000028453A (en) * | 1998-07-07 | 2000-01-28 | Hoya Corp | Method and apparatus for measurement of tension |
CN102519651A (en) * | 2011-12-13 | 2012-06-27 | 清华大学 | Method for determining basic frequency of stay cable when testing cable tension of cable stayed bridge by using vibration method |
RU2015134616A (en) * | 2015-08-17 | 2017-02-21 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Сибирский государственный университет путей сообщения" (СГУПС) | The method of determining the tension force of the cable-stayed bridge element |
CN107421672A (en) * | 2017-05-17 | 2017-12-01 | 苏交科集团股份有限公司 | A kind of weighting Suo Li computational methods that peak is searched based on vibration frequency universe |
JP2018013361A (en) * | 2016-07-19 | 2018-01-25 | 株式会社スエロテック | Tensile force measurement device, tensile force measurement method, and tensile force measurement program |
CN109883596A (en) * | 2019-01-16 | 2019-06-14 | 武汉楚云端信息科技有限责任公司 | A kind of Suo Li automation real-time resolving method based on Spectrum Method |
CN111860397A (en) * | 2020-07-28 | 2020-10-30 | 浙江博远电子科技有限公司 | Cable force calculation method and calculation device for inhaul cable |
CN112100713A (en) * | 2020-08-18 | 2020-12-18 | 东南大学 | Cable force automatic identification method based on variable-pitch grid |
-
2023
- 2023-03-24 CN CN202310294935.4A patent/CN116026511B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000028453A (en) * | 1998-07-07 | 2000-01-28 | Hoya Corp | Method and apparatus for measurement of tension |
CN102519651A (en) * | 2011-12-13 | 2012-06-27 | 清华大学 | Method for determining basic frequency of stay cable when testing cable tension of cable stayed bridge by using vibration method |
RU2015134616A (en) * | 2015-08-17 | 2017-02-21 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Сибирский государственный университет путей сообщения" (СГУПС) | The method of determining the tension force of the cable-stayed bridge element |
JP2018013361A (en) * | 2016-07-19 | 2018-01-25 | 株式会社スエロテック | Tensile force measurement device, tensile force measurement method, and tensile force measurement program |
CN107421672A (en) * | 2017-05-17 | 2017-12-01 | 苏交科集团股份有限公司 | A kind of weighting Suo Li computational methods that peak is searched based on vibration frequency universe |
CN109883596A (en) * | 2019-01-16 | 2019-06-14 | 武汉楚云端信息科技有限责任公司 | A kind of Suo Li automation real-time resolving method based on Spectrum Method |
CN111860397A (en) * | 2020-07-28 | 2020-10-30 | 浙江博远电子科技有限公司 | Cable force calculation method and calculation device for inhaul cable |
CN112100713A (en) * | 2020-08-18 | 2020-12-18 | 东南大学 | Cable force automatic identification method based on variable-pitch grid |
Non-Patent Citations (2)
Title |
---|
Study on the Parameter Vibration of the Large-Span Cable-Stayed Bridges;Li Feng-chen et.al;2009 Sixth International Conference on Information Technology: New Generations;全文 * |
斜拉索实测索力计算方法研究;何长江等;城市道桥与防洪(第第12期期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116026511A (en) | 2023-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111681256B (en) | Image edge detection method, image edge detection device, computer equipment and readable storage medium | |
CN111724599B (en) | Method, device, equipment and medium for acquiring safe driving behavior evaluation data | |
CN110706210B (en) | Deep learning-based rebar counting method and device | |
CN116128312B (en) | Dam safety early warning method and system based on monitoring data analysis | |
CN103871420B (en) | The signal processing method of microphone array and device | |
KR101649243B1 (en) | Method and apparatus for detecting correctness of pitch period | |
CN113421586B (en) | Dream recognition method and device and electronic equipment | |
US9570060B2 (en) | Techniques of audio feature extraction and related processing apparatus, method, and program | |
CN104251934A (en) | Harmonic analysis method and apparatus, and method and apparatus for determining clutter in harmonic wave | |
CN114900256B (en) | Communication scene recognition method and device | |
CN110853677B (en) | Drumbeat beat recognition method and device for songs, terminal and non-transitory computer readable storage medium | |
CN110689885B (en) | Machine synthesized voice recognition method, device, storage medium and electronic equipment | |
CN116026511B (en) | Cable force automatic identification method, system, computer and readable storage medium | |
WO2019201029A1 (en) | Candidate box update method and apparatus | |
CN115712029A (en) | Electromagnetic compatibility test evaluation method for train and trackside axle counter and related equipment | |
AU2012244118B2 (en) | Method and System for Identifying Events of Digital Signal | |
CN107527611A (en) | MFCC audio recognition methods, storage medium, electronic equipment and system | |
CN110659620B (en) | Filtering noise reduction method based on fuzzy control | |
CN111897851A (en) | Abnormal data determination method and device, electronic equipment and readable storage medium | |
CN114882912B (en) | Method and device for testing transient defects of time domain of acoustic signal | |
CN113514147B (en) | Vehicle noise identification method, system, equipment and computer readable storage medium | |
CN113569703A (en) | Method and system for judging true segmentation point, storage medium and electronic equipment | |
CN112560970A (en) | Abnormal picture detection method, system, equipment and storage medium based on self-coding | |
CN116796130A (en) | Bridge vibration low-frequency reconstruction denoising method, system, computer and storage medium | |
CN116415195A (en) | Abnormal data identification method, system, computer and readable storage medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |