A calculation method, device and equipment for simultaneously solving cable force and length of cable rod
By obtaining the cable bending stiffness design value and power spectrum diagram of the cable rod, and using iterative rules to calculate the cable rod length and cable force, the problem of difficult calculating the cable rod length in the existing technology is solved, and accurate cable rod length determination is achieved.
Patent Information
- Application Number
- CN202510045433.7
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2025-01-13
- Publication Date
- 2025-05-13
- Estimated Expiration
- 2045-01-13
AI Technical Summary
Existing cable rods are difficult to determine cable rod length based on non-rigid constraints or support requirements.
By obtaining the cable bending stiffness design value of the cable rod and the power spectrum diagram of the cable rod vibration, the length and cable force of the cable rod are calculated using iterative rules until the convergence conditions and threshold range are met.
The accurate determination of the cable rod length based on non-rigid constraints or support needs is achieved, and the problem of difficult calculating the cable rod length in the prior art is solved.
Smart Images

Figure CN119441708B_ABST
Abstract
Description
Technical Field
[0001] The present application relates to the technical field of cable rods, and in particular to a calculation method, device and equipment for synchronously solving the cable force and length of a cable rod. Background Art
[0002] Cable-stayed structure is a common form of building structure, which uses cables (or cable rods) suspended from fulcrums to support the gravity of the building, forming a unique appearance and function.
[0003] In order to avoid cable vibration, the anchorage sections of prestressed cables in buildings and other structural projects are anchored more reliably, which is closer to consolidation. When the length of prestressed cables and cables is not too long, the cable rods do not sag significantly, and the consolidation formula can be used for analysis. However, on bridge cables and prestressed cables, in order to reduce the damage of vibration to the anchorage and reduce the overall vibration of the structure, different forms of non-rigid constraints or supports are placed on the cables or at both ends of the cables, which brings adjustments to the determination of the length of the cables. The length of existing cables is difficult to calculate. Summary of the invention
[0004] The present application provides a calculation method, device and equipment for simultaneously solving the cable force and length of a cable rod, which is used to solve the technical problem that it is difficult to determine the length of the cable rod according to non-rigid constraints or support requirements in existing cable rods.
[0005] In order to achieve the above objectives, this application provides the following technical solutions:
[0006] On the one hand, a method for calculating the cable force and length of a cable rod simultaneously is provided, comprising the following steps:
[0007] Obtaining a design value of the bending stiffness of the cable of the cable rod and a power spectrum diagram of the vibration of the cable rod, and obtaining n frequency data in a series and a frequency order corresponding to each of the frequency data from the power spectrum diagram;
[0008] Acquire k+1 initial assumed length data of the cable rod, and perform iterative calculation using an iterative rule according to each of the initial assumed length data, the n frequency data and the n frequency orders, to obtain initial cable bending stiffness data and initial tension data corresponding to each of the initial assumed length data and satisfying a convergence condition;
[0009] According to the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data, the k+2th length update data is obtained; according to the k+2th length update data, the nth frequency data and the nth frequency orders, an iterative calculation is performed using an iterative rule to obtain the k+2th cable bending stiffness data and the k+2th tension data that meet the convergence condition;
[0010] If the ratio of the k+2th cable bending stiffness data to the cable bending stiffness design value is within the threshold data range, the k+2th length update data is used as the length of the cable rod and the k+2th tension data is used as the cable force of the cable rod; if the ratio of the k+2th cable bending stiffness data to the cable bending stiffness design value is not within the threshold data range, the k+2th length update data and the k+2th cable bending stiffness data are used as initial assumed length data to recalculate the length update data until the ratio of the iteratively calculated cable bending stiffness data to the cable bending stiffness design value is within the threshold data range and the corresponding length update data is used as the length of the cable rod;
[0011] Wherein, k is a natural number not equal to 0, n is a natural number greater than 3, and the convergence condition is that the change rate of two adjacent tension data is less than a set threshold.
[0012] Preferably, the content of the iteration rule includes:
[0013] Obtaining the cable line density of the cable rod, and performing fitting inference using a first vibration least squares fitting formula according to the initial assumed length data or the length update data, the cable line density, the n frequency data, and the n frequency orders, to obtain a first tension and a first cable bending stiffness;
[0014] Calculating using a frequency angle formula according to the cable line density, the first tension, the first cable bending stiffness, and the frequency data corresponding to the frequency order to obtain a frequency angle corresponding to the frequency order;
[0015] According to the initial assumed length data or the length update data, the cable line density, the n frequency data, the n frequency orders and the n frequency angles, a second vibration least square method fitting formula is used to perform fitting inference to obtain a second tension and a second cable bending stiffness;
[0016] Determine whether the convergence condition is met according to the first tension and the second tension. If so, use the second tension and the second cable bending stiffness as the corresponding initial cable bending stiffness data and initial tension data, or use the second tension and the second cable bending stiffness as the corresponding cable bending stiffness data and tension data; if not, use the second tension and the second cable bending stiffness as the first tension and the first cable bending stiffness to update the frequency angle, tension and cable bending stiffness until the updated tension meets the convergence condition and uses the updated tension and cable bending stiffness as the corresponding initial cable bending stiffness data and initial tension data, or uses the updated tension and cable bending stiffness as the corresponding cable bending stiffness data and tension data.
[0017] Preferably, the first vibration least squares fitting formula is:
[0018]
[0019] The second vibration least squares fitting formula is:
[0020]
[0021] Where EI0 is the bending stiffness of the first cable, L is the initial assumed length data or length update data, ρ is the cable line density, i and f i are the frequency order and the frequency data of the i-th order, T0 is the first tension, EI k is the second cable bending stiffness or the cable bending stiffness calculated for the kth update, T k is the second tension or the tension calculated for the kth update, is the frequency angle of the i-th order.
[0022] Preferably, the frequency angle formula is:
[0023]
[0024] In the formula, i and f i are the frequency order and the frequency data of the i-th order, T is the tension, EI is the bending stiffness of the cable, ρ is the linear density of the cable, is the frequency angle of the i-th order.
[0025] Preferably, the calculation method for simultaneously solving the cable tension and length of the cable rod comprises: using the length update formula to calculate according to the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data, to obtain the k+2th length update data; the length update formula is: L k+2 =((L k+1 -L k ) / (EI k+1 -EI k ))EI+(L k -EI k (L k+1 -L k ) / (EI k+1 -EI k )), where L k+2 Update data for the k+2th length, L k+1 is the k+1th initial assumption length data, L k is the kth initial hypothesis length data, EI k+1is the k+1th initial cable bending stiffness data, EI k is the kth initial cable bending stiffness data, and EI is the design value of the cable bending stiffness.
[0026] In another aspect, a computing device for synchronously solving the cable force and length of a cable rod is provided, comprising a data acquisition module, a first iterative computing module, a second iterative computing module and a length module;
[0027] The data acquisition module is used to obtain the design value of the bending stiffness of the cable of the cable rod and the power spectrum diagram of the cable rod vibration, and obtain n frequency data in a series and the frequency order corresponding to each frequency data from the power spectrum diagram;
[0028] The first iterative calculation module is used to obtain k+1 initial assumed length data of the cable rod, and perform iterative calculation using an iterative rule according to each of the initial assumed length data, the n frequency data and the n frequency orders to obtain initial cable bending stiffness data and initial tension data corresponding to each of the initial assumed length data and satisfying a convergence condition;
[0029] The second iterative calculation module is used to calculate according to the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data to obtain the k+2th length update data; iteratively calculate according to the k+2th length update data, the nth frequency data and the nth frequency orders using an iterative rule to obtain the k+2th cable bending stiffness data and the k+2th tension data that meet the convergence condition;
[0030] The length module is used to use the k+2th length update data as the length of the cable rod and the k+2th tension data as the cable force of the cable rod according to the ratio of the k+2th cable bending stiffness data to the cable bending stiffness design value being within the threshold data range; or to use the k+2th length update data and the k+2th cable bending stiffness data as the initial assumed length data to recalculate the length update data according to the ratio of the k+2th cable bending stiffness data to the cable bending stiffness design value not being within the threshold data range, until the ratio of the iteratively calculated cable bending stiffness data to the cable bending stiffness design value is within the threshold data range and the corresponding length update data is used as the length of the cable rod;
[0031] Wherein, k is a natural number not equal to 0, n is a natural number greater than 3, and the convergence condition is that the change rate of two adjacent tension data is less than a set threshold.
[0032] Preferably, the content of the iteration rule includes:
[0033] Obtaining the cable line density of the cable rod, and performing fitting inference using a first vibration least squares fitting formula according to the initial assumed length data or the length update data, the cable line density, the n frequency data, and the n frequency orders, to obtain a first tension and a first cable bending stiffness;
[0034] Calculating using a frequency angle formula according to the cable line density, the first tension, the first cable bending stiffness, and the frequency data corresponding to the frequency order to obtain a frequency angle corresponding to the frequency order;
[0035] According to the initial assumed length data or the length update data, the cable line density, the n frequency data, the n frequency orders and the n frequency angles, a second vibration least square method fitting formula is used to perform fitting inference to obtain a second tension and a second cable bending stiffness;
[0036] Determine whether the convergence condition is met based on the first tension and the second tension. If so, use the second tension and the second cable bending stiffness as the corresponding initial cable bending stiffness data and initial tension data, or use the second tension and the second cable bending stiffness as the corresponding cable bending stiffness data and tension data; if not, use the second tension and the second cable bending stiffness as the first tension and the first cable bending stiffness to update the frequency angle, tension and cable bending stiffness until the updated tension meets the convergence condition and uses the updated tension and cable bending stiffness as the corresponding initial cable bending stiffness data and initial tension data, or uses the updated tension and cable bending stiffness as the corresponding cable bending stiffness data and tension data.
[0037] Preferably, the first vibration least squares fitting formula is:
[0038]
[0039] The second vibration least squares fitting formula is:
[0040]
[0041] The frequency angle formula is:
[0042]
[0043] Where T is the tension, EI is the cable bending stiffness, EI0 is the first cable bending stiffness, L is the initial assumed length data or length update data, ρ is the cable linear density, i and f i are the frequency order and the frequency data of the i-th order, T0 is the first tension, EI kis the second cable bending stiffness or the cable bending stiffness calculated for the kth update, T k is the second tension or the tension calculated for the kth update, is the frequency angle of the i-th order.
[0044] Preferably, the second iterative calculation module is further used to calculate the k+2th length update data by using the length update formula according to the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data; the length update formula is: L k+2 =((L k+1 -L k ) / (EI k+1 -EI k ))EI+(L k -EI k (L k+1 -L k ) / (EI k+1 -EI k )), where L k+2 Update data for the k+2th length, L k+1 is the k+1th initial assumption length data, L k is the kth initial hypothesis length data, EI k+1 is the k+1th initial cable bending stiffness data, EI k is the kth initial cable bending stiffness data, and EI is the design value of the cable bending stiffness.
[0045] In another aspect, a terminal device is provided, comprising a processor and a memory;
[0046] The memory is used to store program codes and transmit the program codes to the processor;
[0047] The processor is used to execute the above-mentioned calculation method for synchronously solving the cable rod tension and length according to the instructions in the program code.
[0048] The method, device and equipment for calculating the cable tension and length of a cable rod synchronously include obtaining the cable bending stiffness design value of the cable rod and the power spectrum diagram of the cable rod vibration, obtaining n frequency data in a series and the frequency order corresponding to each frequency data from the power spectrum diagram; obtaining k+1 initial assumed length data of the cable rod, performing iterative calculation according to each initial assumed length data, n frequency data and n frequency orders using an iterative rule, and obtaining initial cable bending stiffness data and initial tension data corresponding to each initial assumed length data and satisfying a convergence condition; calculating according to the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data, and obtaining the k+2th length update data; and calculating according to the k+2th length update data. The k+2th length update data, n frequency data and n frequency orders are iteratively calculated using an iteration rule to obtain the k+2th cable bending stiffness data and the k+2th tension data that meet the convergence condition; if the ratio of the k+2th cable bending stiffness data to the design value of the cable bending stiffness is within the threshold data range, the k+2th length update data is used as the length of the cable rod and the k+2th tension data is used as the cable force of the cable rod; if the ratio of the k+2th cable bending stiffness data to the design value of the cable bending stiffness is not within the threshold data range, the k+2th length update data and the k+2th cable bending stiffness data are used as the initial assumed length data to recalculate the length update data until the ratio of the iteratively calculated cable bending stiffness data to the design value of the cable bending stiffness is within the threshold data range and the corresponding length update data is used as the length of the cable rod.
[0049] It can be seen from the above technical scheme that the present application has the following advantages: the calculation method for synchronously solving the cable tension and length of the cable rod first calculates the corresponding initial cable bending stiffness data, initial tension data and length update data through n frequency data in an approximate geometric progression and the set initial assumed length data, and then calculates the k+2th cable bending stiffness data according to the length update data, and finally determines the corresponding k+2th length update data as the length of the cable rod according to whether the ratio of the calculated k+2th cable bending stiffness data to the cable bending stiffness design value is within the threshold data range. The length of the cable rod meets the non-rigid constraints or support requirements of the cable rod engineering requirements, which solves the technical problem that it is difficult to determine the length of the cable rod according to the non-rigid constraints or support requirements of the existing cable rods. BRIEF DESCRIPTION OF THE DRAWINGS
[0050] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.
[0051] Figure 1 A flowchart of the steps of the method for calculating the cable force and length of the cable rod simultaneously according to an embodiment of the present application;
[0052] Figure 2 A framework flow chart of a method for calculating the cable force and length of a cable rod simultaneously as described in an embodiment of the present application;
[0053] Figure 3 A framework flow chart of the iteration rules in the method for calculating the cable force and length of the cable rod simultaneously described in the embodiment of the present application;
[0054] Figure 4 It is a schematic diagram of the framework of the calculation device for simultaneously solving the cable force and length of the cable rod according to the embodiment of the present application;
[0055] Figure 5 A schematic diagram of a terminal device described in an embodiment of the present application. DETAILED DESCRIPTION
[0056] In order to make the purpose, features, and advantages of the invention of this application more obvious and easy to understand, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the drawings in the embodiments of this application. Obviously, the embodiments described below are only part of the embodiments of this application, not all of them. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.
[0057] In the description of the embodiments of the present application, the terms "first" and "second" are used for descriptive purposes only and cannot be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of the features. In the description of the embodiments of the present application, the meaning of "multiple" is two or more, unless otherwise clearly and specifically defined.
[0058] In the embodiments of the present application, unless otherwise clearly specified and limited, the terms "installed", "connected", "connected", "fixed" and the like should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral one; it can be a mechanical connection or an electrical connection; it can be a direct connection or an indirect connection through an intermediate medium, it can be the internal connection of two elements or the interaction relationship between two elements. For those of ordinary skill in the art, the specific meanings of the above terms in the embodiments of the present application can be understood according to the specific circumstances.
[0059] The embodiments of the present application provide a calculation method, device and equipment for simultaneously solving the cable force and length of a cable rod, which solves the technical problem that it is difficult to determine the length of the cable rod according to non-rigid constraints or support requirements in existing cable rods.
[0060] Embodiment 1:
[0061] Figure 1 This is a flowchart of the steps of the method for calculating the cable force and length of the cable rod simultaneously according to an embodiment of the present application. Figure 2 The present invention is a framework flow chart of a method for calculating the cable tension and length of a cable rod simultaneously as described in an embodiment of the present application.
[0062] like Figure 1 and Figure 2 As shown, the embodiment of the present application provides a method for calculating the cable force and length of a cable rod simultaneously, comprising the following steps:
[0063] S1. Obtain the design value of the bending stiffness of the cable of the cable rod and the power spectrum diagram of the cable rod vibration, and obtain n frequency data in a series and the frequency order corresponding to each frequency data from the power spectrum diagram.
[0064] It should be noted that in step S1, first, the cable bending stiffness design value EI of the cable rod and the power spectrum of the cable rod vibration are obtained; second, n frequency data presenting an approximate geometric progression and the corresponding n frequency orders are obtained from the power spectrum. In this embodiment, the cable bending stiffness design value EI can be set according to demand, and the specific parameters of the cable bending stiffness design value are not limited here.
[0065] In an embodiment of the present application, obtaining a power spectrum diagram of cable rod vibration includes: collecting vibration signals of various orders by a signal collection element disposed at a middle position of the cable rod; and analyzing and processing the vibration signals of various orders using an autocorrelation power spectrum to obtain a power spectrum diagram of cable rod vibration.
[0066] It should be noted that the signal acquisition element can be selected as an acceleration sensor. In this embodiment, the acceleration sensor is fixed on the cable rod, and the fixing position needs to avoid the stationary point of the first 10th order vibration of the cable rod, and is preferably located in the middle of two adjacent stationary points. For a cable rod with a cable length of L, the stationary point position within the 10th order is: x k0,n0 =k0L / n0, k0=0, 1, ..., n0, n0=1, 2, ....., 10. The vibration signal of the forced free vibration section of the cable rod after being struck is collected by an acceleration sensor, and the collected vibration signal is subjected to autocorrelation power spectrum analysis to obtain a power spectrum diagram of the vibration. In other embodiments, the signal collection element can also be a collection element (such as a vibration sensor) that can collect vibration signals. The vibration sensor includes an inductive vibration sensor, an eddy current vibration sensor, a capacitive vibration sensor, a resistive strain vibration sensor, a piezoelectric vibration sensor, etc. The autocorrelation power spectrum analysis technology is a relatively mature technology in this field and will not be elaborated here.
[0067] In an embodiment of the present application, in step S1 of the method for simultaneously solving the cable tension and length of the cable rod, if n frequency data in a series cannot be obtained from the power spectrum diagram, the installation position of the signal acquisition element on the cable rod is adjusted, and it is checked whether the signal acquisition element is firmly installed on the cable rod, and the vibration signal is recollected to obtain the power spectrum diagram until n frequency data in a series are obtained from the power spectrum diagram.
[0068] S2. Obtain k+1 initial assumed length data of the cable rod, and perform iterative calculations using an iterative rule based on each initial assumed length data, n frequency data, and n frequency orders to obtain initial cable bending stiffness data and initial tension data corresponding to each initial assumed length data and satisfying the convergence condition. Wherein, k is a natural number not equal to 0, n is a natural number greater than 3, and the convergence condition is that the rate of change of two adjacent tension data is less than a set threshold.
[0069] It should be noted that in step S2, first, k+1 initial assumed length data of the cable rod are obtained, and second, initial cable bending stiffness data and initial tension data corresponding to each initial assumed length data and satisfying the convergence condition are obtained by iterative calculation based on the k+1 initial assumed length data and the frequency data and frequency order obtained in step S1. In this embodiment, the threshold is set to 0.1%, and in other embodiments, the threshold can be set as required.
[0070] In the embodiment of the present application, the initial assumed length data of the method for calculating the cable rod tension and length simultaneously includes the first initial assumed length data L1=L when k=1 and the second initial assumed length data L2=L / 2, where L is the length assumption value. Figure 2 As shown, L k , Lk+1 As the initial assumed length data, EI k , EI k+1 As the corresponding initial cable bending stiffness data, T k , T k+1 As the corresponding initial tension data.
[0071] S3. Obtain the k+2th length update data based on the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data; perform iterative calculations based on the k+2th length update data, n frequency data and n frequency orders using an iterative rule to obtain the k+2th cable bending stiffness data and the k+2th tension data that meet the convergence conditions.
[0072] It should be noted that in step S3, the k+2th length update data is calculated first, and then the k+2th cable bending stiffness data and tension data are calculated according to the length update data using an iteration rule to provide data for subsequent steps.
[0073] S4. If the ratio of the k+2th cable bending stiffness data to the design value of the cable bending stiffness is within the threshold data range, the k+2th length update data is used as the length of the cable rod and the k+2th tension data is used as the cable force of the cable rod; if the ratio of the k+2th cable bending stiffness data to the design value of the cable bending stiffness is not within the threshold data range, the k+2th length update data and the k+2th cable bending stiffness data are used as the initial assumed length data to recalculate the length update data until the ratio of the iteratively calculated cable bending stiffness data to the design value of the cable bending stiffness is within the threshold data range and the corresponding length update data is used as the length of the cable rod.
[0074] It should be noted that in step S4, it is first determined whether the ratio of the cable bending stiffness data calculated in step S3 to the cable bending stiffness design value is within the threshold data range, and whether to use the k+2th length update data as the length of the cable rod and the k+2th tension data as the cable force of the cable rod. In this embodiment, if the ratio of the k+2th cable bending stiffness data calculated in step S3 to the cable bending stiffness design value is not within the threshold data range, the k+2th length update data and the k+2th cable bending stiffness data are used as the initial assumed length data to recalculate the length update data and repeat steps S3 and S4 until the ratio of the iteratively calculated cable bending stiffness data to the cable bending stiffness design value is within the threshold data range and the corresponding length update data is used as the length of the cable rod. The threshold data is [0.99, 1.01].
[0075] In an embodiment of the present application, the calculation method for simultaneously solving the cable tension and length of the cable rod first obtains n frequency data in an approximate geometric progression from the power spectrum diagram of the cable rod, and then calculates the corresponding initial cable bending stiffness data, initial tension data and length update data according to the set initial assumed length data, and then calculates the k+2th cable bending stiffness data and the k+2th tension data according to the length update data, and determines the length of the cable rod based on whether the ratio of the calculated k+2th cable bending stiffness data to the design value of the cable bending stiffness is not within the threshold data range.
[0076] The present application provides a calculation method for synchronously solving the cable force and length of a cable rod, comprising obtaining a design value of the cable bending stiffness of the cable rod and a power spectrum diagram of the cable rod vibration, obtaining n frequency data in a series and a frequency order corresponding to each frequency data from the power spectrum diagram; obtaining k+1 initial assumed length data of the cable rod, performing iterative calculations using an iteration rule according to each initial assumed length data, n frequency data and n frequency orders, and obtaining initial cable bending stiffness data and initial tension data corresponding to each initial assumed length data and satisfying a convergence condition; calculating according to the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data, and obtaining the k+2th length update data; calculating according to the k+2th length update data, the nth frequency data and the nth frequency orders, and obtaining the initial cable bending stiffness design value and the corresponding initial cable bending stiffness data; and obtaining the initial cable bending stiffness design value and the initial cable bending stiffness design value according to the k+2th length update data, the nth frequency data and the nth frequency orders. The frequency data and n frequency orders are iteratively calculated using an iteration rule to obtain the k+2th cable bending stiffness data and the k+2th tension data that meet the convergence conditions; if the ratio of the k+2th cable bending stiffness data to the design value of the cable bending stiffness is within the threshold data range, the k+2th length update data is used as the length of the cable rod and the k+2th tension data is used as the cable force of the cable rod; if the ratio of the k+2th cable bending stiffness data to the design value of the cable bending stiffness is not within the threshold data range, the k+2th length update data and the k+2th cable bending stiffness data are used as the initial assumed length data to recalculate the length update data until the ratio of the iteratively calculated cable bending stiffness data to the design value of the cable bending stiffness is within the threshold data range and the corresponding length update data is used as the length of the cable rod. The calculation method for simultaneously solving the cable force and length of the cable rod first calculates the corresponding initial cable bending stiffness data, initial tension data and length update data through n frequency data in an approximate geometric progression and the set initial assumed length data, then calculates the k+2th cable bending stiffness data according to the length update data, and finally determines the corresponding k+2th length update data as the length of the cable rod according to whether the ratio of the calculated k+2th cable bending stiffness data to the cable bending stiffness design value is within the threshold data range. The length of the cable rod meets the non-rigid constraints or support requirements of the cable rod engineering requirements, which solves the technical problem that it is difficult to determine the length of the cable rod according to the non-rigid constraints or support requirements of the existing cable rods.
[0077] Figure 3 This is a framework flow chart of the iteration rules in the calculation method for simultaneously solving the cable force and length of the cable rod described in an embodiment of the present application.
[0078] like Figure 3 As shown, in one embodiment of the present application, the content of the iteration rule includes:
[0079] Obtain the cable line density of the cable rod, and use the first vibration least square method fitting formula to perform fitting inference according to the initial assumed length data or the length update data, the cable line density, the n frequency data, and the n frequency orders to obtain the first tension and the first cable bending stiffness;
[0080] The frequency angle corresponding to the frequency order is calculated by using a frequency angle formula according to the cable line density, the first tension, the first cable bending stiffness and the frequency data corresponding to the frequency order;
[0081] According to the initial assumed length data or the length update data, the cable line density, the n frequency data, the n frequency orders and the n frequency angles, the second vibration least square method fitting formula is used for fitting and inference to obtain the second tension and the second cable bending stiffness;
[0082] Determine whether the convergence condition is met according to the first tension and the second tension. If so, use the second tension and the second cable bending stiffness as the corresponding initial cable bending stiffness data and initial tension data, or use the second tension and the second cable bending stiffness as the corresponding cable bending stiffness data and tension data; if not, use the second tension and the second cable bending stiffness as the first tension and the first cable bending stiffness to update the frequency angle, tension and cable bending stiffness until the updated tension meets the convergence condition and uses the updated tension and cable bending stiffness as the corresponding initial cable bending stiffness data and initial tension data, or uses the updated tension and cable bending stiffness as the corresponding cable bending stiffness data and tension data.
[0083] It should be noted that in the iterative rules of the calculation method for simultaneously solving the cable rod tension and length, in the process of calculating the first tension and the first cable bending stiffness, the frequency order is first used as the horizontal coordinate and the frequency data as the vertical coordinate, and the n frequency data and the corresponding frequency order are used to construct a first scatter plot. The first vibration least squares fitting formula is used to fit the first scatter plot, and the fitting coefficients are the inferred first tension T1 and the first cable bending stiffness EI1. Similarly, in the process of calculating the second tension and the second cable bending stiffness, the ratio of the frequency angle to the pi is first used as a coefficient, the difference between the frequency order and the coefficient is used as the horizontal coordinate, and the frequency data is used as the vertical coordinate. The n frequency data and the corresponding difference are used to construct a second scatter plot. The second vibration least squares fitting formula is used to fit in the second scatter plot. The fitting coefficient is the inferred second tension T2 and the second cable bending stiffness EI2. Under the same initial assumed length data or length update data, according to whether the change rate between the first tension and the second tension is less than the set threshold, it is judged whether the calculated second tension and the second cable bending stiffness can be used as the initial cable bending stiffness data and initial tension data or the cable bending stiffness data and tension data corresponding to the initial assumed length data or the length update data. If the calculated second tension and the second cable bending stiffness do not meet the convergence condition, the frequency angle formula is re-calculated based on the second tension and the second cable bending stiffness to obtain an updated frequency angle, thereby updating the second scatter plot, and iteratively updating until the second tension and the second cable bending stiffness that meet the convergence condition are obtained.
[0084] In an embodiment of the present application, the first vibration least squares fitting formula is:
[0085]
[0086] The second vibration least squares fitting formula is:
[0087]
[0088] Where EI0 is the bending stiffness of the first cable, L is the initial assumed length data or length update data, ρ is the cable line density, i and f i are the frequency order and the frequency data of the i-th order, T0 is the first tension, EI k is the second cable bending stiffness or the cable bending stiffness calculated for the kth update, T k is the second tension or the tension calculated for the kth update, is the frequency angle of the i-th order.
[0089] In the embodiment of the present application, the frequency angle formula is:
[0090]
[0091] In the formula, i and f i are the frequency order and the frequency data of the i-th order, T is the tension, EI is the bending stiffness of the cable, ρ is the linear density of the cable, is the frequency angle of the i-th order.
[0092] In one embodiment of the present application, the calculation method for simultaneously solving the cable tension and length of the cable rod includes: using the length update formula to calculate according to the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data, to obtain the k+2th length update data; the length update formula is: L k+2 =((L k+1 -L k ) / (EI k+1 -EI k ))EI+(L k -EI k (L k+1 -L k ) / (EI k+1 -EI k )), where L k+2 Update data for the k+2th length, L k+1 is the k+1th initial assumption length data, L k is the kth initial hypothesis length data, EI k+1 is the k+1th initial cable bending stiffness data, EI k is the kth initial cable bending stiffness data, and EI is the design value of the cable bending stiffness.
[0093] In the embodiment of the present application, the calculation method for solving the cable force and length of the cable rod simultaneously assumes the length value of the cable rod first, adopts the tension and cable bending stiffness corresponding to the iteration rule, and when the ratio of the k+2nd cable bending stiffness data to the cable bending stiffness design value is calculated to be within the threshold data range, the k+2nd length update data is used as the length of the cable rod that meets the engineering requirements. When the ratio of the k+2nd cable bending stiffness data to the cable bending stiffness design value is calculated to be not within the threshold data range, the k+2nd length update data and the k+2nd cable bending stiffness data are used to infer the next length update data using the length update formula, and the contents of step S3 and step S4 are repeatedly executed until the k+2nd cable bending stiffness data is iteratively calculated to be within the threshold data range.
[0094] Embodiment 2:
[0095] Figure 4 A schematic diagram of the framework of a calculation device for simultaneously solving the cable tension and length of a cable rod as described in an embodiment of the present application.
[0096] like Figure 4 As shown, an embodiment of the present application provides a calculation device for synchronously solving the cable force and length of a cable rod, including a data acquisition module 10, a first iterative calculation module 20, a second iterative calculation module 30 and a length module 40;
[0097] The data acquisition module 10 is used to obtain the design value of the bending stiffness of the cable of the cable rod and the power spectrum diagram of the cable rod vibration, and obtain n frequency data in a series and the frequency order corresponding to each frequency data from the power spectrum diagram;
[0098] The first iterative calculation module 20 is used to obtain k+1 initial assumed length data of the cable rod, and perform iterative calculation using an iterative rule according to each initial assumed length data, n frequency data and n frequency orders to obtain initial cable bending stiffness data and initial tension data corresponding to each initial assumed length data and satisfying a convergence condition;
[0099] The second iterative calculation module 30 is used to calculate according to the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data to obtain the k+2th length update data; iteratively calculate according to the k+2th length update data, the n frequency data and the n frequency orders using an iterative rule to obtain the k+2th cable bending stiffness data and the k+2th tension data that meet the convergence condition;
[0100] The length module 40 is used to use the k+2th length update data as the length of the cable rod and the k+2th tension data as the cable force of the cable rod according to the ratio of the k+2th cable bending stiffness data to the cable bending stiffness design value being within the threshold data range; or to use the k+2th length update data and the k+2th cable bending stiffness data as the initial assumed length data to recalculate the length update data according to the ratio of the k+2th cable bending stiffness data to the cable bending stiffness design value not being within the threshold data range, until the ratio of the iteratively calculated cable bending stiffness data to the cable bending stiffness design value is within the threshold data range and the corresponding length update data is used as the length of the cable rod;
[0101] Wherein, k is a natural number not equal to 0, n is a natural number greater than 3, and the convergence condition is that the change rate of two adjacent tension data is less than the set threshold.
[0102] It should be noted that the contents of the modules in the calculation device for synchronously solving the cable rod tension and length in the second embodiment correspond to the steps of the calculation method for synchronously solving the cable rod tension and length in the first embodiment. The steps of the calculation method for synchronously solving the cable rod tension and length have been described in the first embodiment, and the module contents of the calculation device for synchronously solving the cable rod tension and length are not described in detail in this embodiment. In this embodiment, the calculation device for synchronously solving the cable rod tension and length realizes the calculation of determining the cable rod length according to the design value of the cable bending stiffness through the data acquisition module 10, the first iterative calculation module 20, the second iterative calculation module 30 and the length module 40.
[0103] In the embodiment of the present application, the content of the iteration rule includes:
[0104] Obtain the cable line density of the cable rod, and use the first vibration least square method fitting formula to perform fitting inference according to the initial assumed length data or the length update data, the cable line density, the n frequency data, and the n frequency orders to obtain the first tension and the first cable bending stiffness;
[0105] The frequency angle corresponding to the frequency order is calculated by using a frequency angle formula according to the cable line density, the first tension, the first cable bending stiffness and the frequency data corresponding to the frequency order;
[0106] According to the initial assumed length data or the length update data, the cable line density, the n frequency data, the n frequency orders and the n frequency angles, the second vibration least square method fitting formula is used for fitting and inference to obtain the second tension and the second cable bending stiffness;
[0107] Determine whether the convergence condition is met according to the first tension and the second tension. If so, use the second tension and the second cable bending stiffness as the corresponding initial cable bending stiffness data and initial tension data, or use the second tension and the second cable bending stiffness as the corresponding cable bending stiffness data and tension data; if not, use the second tension and the second cable bending stiffness as the first tension and the first cable bending stiffness to update the frequency angle, tension and cable bending stiffness until the updated tension meets the convergence condition and uses the updated tension and cable bending stiffness as the corresponding initial cable bending stiffness data and initial tension data, or uses the updated tension and cable bending stiffness as the corresponding cable bending stiffness data and tension data.
[0108] In the embodiment of the present application, the first vibration least squares fitting formula is:
[0109]
[0110] The second vibration least squares fitting formula is:
[0111]
[0112] The frequency angle formula is:
[0113]
[0114] Where T is the tension, EI is the cable bending stiffness, EI0 is the first cable bending stiffness, L is the initial assumed length data or length update data, ρ is the cable linear density, i and f i are the frequency order and the frequency data of the i-th order, T0 is the first tension, EI k is the second cable bending stiffness or the cable bending stiffness calculated for the kth update, T k is the second tension or the tension calculated for the kth update, is the frequency angle of the i-th order.
[0115] In the embodiment of the present application, the second iterative calculation module is further used to calculate the k+2th length update data using the length update formula according to the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data; the length update formula is: L k+2 =((L k+1 -L k ) / (EI k+1 -EI k ))EI+(L k -EI k (L k+1 -L k ) / (EI k+1 -EI k )), where L k+2 Update data for the k+2th length, L k+1 is the k+1th initial assumption length data, L k is the kth initial hypothesis length data, EI k+1 is the k+1th initial cable bending stiffness data, EI k is the kth initial cable bending stiffness data, and EI is the design value of the cable bending stiffness.
[0116] Embodiment three:
[0117] Figure 5 A schematic diagram of a terminal device described in an embodiment of the present application.
[0118] like Figure 5 As shown, an embodiment of the present application provides a terminal device, including a processor and a memory;
[0119] A memory, used for storing program codes and transmitting the program codes to a processor;
[0120] The processor is used to execute the above-mentioned calculation method for synchronously solving the cable rod tension and length according to the instructions in the program code.
[0121] It should be noted that the processor is used to execute the steps of the above-mentioned embodiment of a method for calculating the cable force and length of the cable rod synchronously according to the instructions in the program code. Alternatively, the processor implements the functions of each module / unit in the above-mentioned system / device embodiments when executing the computer program.
[0122] Exemplarily, the computer program may be divided into one or more modules / units, one or more modules / units are stored in a memory and executed by a processor to complete the present application. One or more modules / units may be a series of computer program instruction segments capable of completing a specific function, and the instruction segments are used to describe the execution process of the computer program in a terminal device.
[0123] The terminal device may be a computing device such as a desktop computer, a notebook, a PDA, or a cloud server. The terminal device may include, but is not limited to, a processor and a memory. Those skilled in the art will appreciate that this does not constitute a limitation on the terminal device, and may include more or fewer components than shown in the figure, or a combination of certain components, or different components. For example, the terminal device may also include input and output devices, network access devices, buses, etc.
[0124] The processor may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc.
[0125] The memory may be an internal storage unit of the terminal device, such as a hard disk or memory of the terminal device. The memory may also be an external storage device of the terminal device, such as a plug-in hard disk, a smart memory card (SMC), a secure digital (SD) card, a flash card, etc. equipped on the terminal device. Furthermore, the memory may include both an internal storage unit of the terminal device and an external storage device. The memory is used to store computer programs and other programs and data required by the terminal device. The memory may also be used to temporarily store data that has been output or is to be output.
[0126] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.
[0127] In the several embodiments provided in the present application, it should be understood that the disclosed systems, devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or units, which can be electrical, mechanical or other forms.
[0128] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
[0129] In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.
[0130] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application is essentially or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product. The computer software product is stored in a storage medium, including several instructions for a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present application. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM), random access memory (RAM), disk or optical disk, etc., various media that can store program codes.
[0131] As described above, the above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application.
Claims
1. A method for calculating the cable force and length of a cable rod simultaneously, characterized in that: The following steps are involved: Obtaining a design value of the bending stiffness of the cable of the cable rod and a power spectrum diagram of the vibration of the cable rod, and obtaining n frequency data in a series and a frequency order corresponding to each of the frequency data from the power spectrum diagram; Acquire k+1 initial assumed length data of the cable rod, and perform iterative calculation using an iterative rule according to each of the initial assumed length data, the n frequency data and the n frequency orders, to obtain initial cable bending stiffness data and initial tension data corresponding to each of the initial assumed length data and satisfying a convergence condition; According to the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data, the k+2th length update data is obtained; according to the k+2th length update data, the nth frequency data and the nth frequency orders, an iterative calculation is performed using an iterative rule to obtain the k+2th cable bending stiffness data and the k+2th tension data that meet the convergence condition; If the ratio of the k+2th cable bending stiffness data to the cable bending stiffness design value is within the threshold data range, the k+2th length update data is used as the length of the cable rod and the k+2th tension data is used as the cable force of the cable rod; if the ratio of the k+2th cable bending stiffness data to the cable bending stiffness design value is not within the threshold data range, the k+2th length update data and the k+2th cable bending stiffness data are used as initial assumed length data to recalculate the length update data until the ratio of the iteratively calculated cable bending stiffness data to the cable bending stiffness design value is within the threshold data range and the corresponding length update data is used as the length of the cable rod; Wherein, k is a natural number not equal to 0, n is a natural number greater than 3, and the convergence condition is that the change rate of two adjacent tension data is less than a set threshold; The content of the iteration rules includes: Obtaining the cable line density of the cable rod, and performing fitting inference using a first vibration least squares fitting formula according to the initial assumed length data or the length update data, the cable line density, the n frequency data, and the n frequency orders, to obtain a first tension and a first cable bending stiffness; Calculating using a frequency angle formula according to the cable line density, the first tension, the first cable bending stiffness, and the frequency data corresponding to the frequency order to obtain a frequency angle corresponding to the frequency order; According to the initial assumed length data or the length update data, the cable line density, the n frequency data, the n frequency orders and the n frequency angles, a second vibration least square method fitting formula is used to perform fitting inference to obtain a second tension and a second cable bending stiffness; Determine whether the convergence condition is met based on the first tension and the second tension. If so, use the second tension and the second cable bending stiffness as the corresponding initial cable bending stiffness data and initial tension data, or use the second tension and the second cable bending stiffness as the corresponding cable bending stiffness data and tension data; if not, use the second tension and the second cable bending stiffness as the first tension and the first cable bending stiffness to update the frequency angle, tension and cable bending stiffness until the updated tension meets the convergence condition and uses the updated tension and cable bending stiffness as the corresponding initial cable bending stiffness data and initial tension data, or uses the updated tension and cable bending stiffness as the corresponding cable bending stiffness data and tension data.
2. The method for calculating the cable force and length of a cable rod synchronously according to claim 1, characterized in that: The first vibration least squares fitting formula is: ; The second vibration least squares fitting formula is: ; Where EI0 is the bending stiffness of the first cable, L is the initial assumed length data or length update data, ρ is the cable line density, i and f i are the frequency order and the frequency data of the i-th order, T0 is the first tension, EI k is the second cable bending stiffness or the cable bending stiffness calculated for the kth update, T k is the second tension or the tension calculated for the kth update, is the frequency angle of the i-th order.
3. The method for calculating the cable force and length of a cable rod synchronously according to claim 1, characterized in that: The frequency angle formula is: ; In the formula, i and f i are the frequency order and the frequency data of the i-th order, T is the tension, EI is the bending stiffness of the cable, ρ is the linear density of the cable, is the frequency angle of the i-th order.
4. The method for calculating the cable force and length of a cable rod synchronously according to claim 1, characterized in that: include: Calculate the k+2th length update data using a length update formula according to the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data; The length update formula is: L k+2 =((L k+1 -L k ) / (EI k+1 -EI k ))EI+(L k -EI k (L k+1 -L k ) / (EI k+1 -EI k )), where L k+2 Update data for the k+2th length, L k+1 is the k+1th initial assumption length data, L k is the kth initial hypothesis length data, EI k+1 is the k+1th initial cable bending stiffness data, EI k is the kth initial cable bending stiffness data, and EI is the design value of the cable bending stiffness.
5. A calculation device for simultaneously solving the cable force and length of a cable rod, characterized in that: It includes a data acquisition module, a first iterative calculation module, a second iterative calculation module and a length module; The data acquisition module is used to obtain the design value of the bending stiffness of the cable of the cable rod and the power spectrum diagram of the cable rod vibration, and obtain n frequency data in a series and the frequency order corresponding to each frequency data from the power spectrum diagram; The first iterative calculation module is used to obtain k+1 initial assumed length data of the cable rod, and perform iterative calculation using an iterative rule according to each of the initial assumed length data, the n frequency data and the n frequency orders to obtain initial cable bending stiffness data and initial tension data corresponding to each of the initial assumed length data and satisfying a convergence condition; The second iterative calculation module is used to calculate according to the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data to obtain the k+2th length update data; iteratively calculate according to the k+2th length update data, the nth frequency data and the nth frequency orders using an iterative rule to obtain the k+2th cable bending stiffness data and the k+2th tension data that meet the convergence condition; The length module is used for using the k+2th length update data as the length of the cable rod and the k+2th tension data as the cable force of the cable rod according to the ratio of the k+2th cable bending stiffness data to the cable bending stiffness design value being within the threshold data range; or according to the ratio of the k+2th cable bending stiffness data to the cable bending stiffness design value not being within the threshold data range, the k+2th length update data and the k+2th cable bending stiffness data are used as the initial assumed length data to recalculate the length update data, until the ratio of the iteratively calculated cable bending stiffness data to the cable bending stiffness design value is within the threshold data range and the corresponding length update data is used as the length of the cable rod; Wherein, k is a natural number not equal to 0, n is a natural number greater than 3, and the convergence condition is that the change rate of two adjacent tension data is less than a set threshold; The content of the iteration rules includes: Obtaining the cable line density of the cable rod, and performing fitting inference using a first vibration least squares fitting formula according to the initial assumed length data or the length update data, the cable line density, the n frequency data, and the n frequency orders, to obtain a first tension and a first cable bending stiffness; Calculating using a frequency angle formula according to the cable line density, the first tension, the first cable bending stiffness, and the frequency data corresponding to the frequency order to obtain a frequency angle corresponding to the frequency order; According to the initial assumed length data or the length update data, the cable line density, the n frequency data, the n frequency orders and the n frequency angles, a second vibration least square method fitting formula is used to perform fitting inference to obtain a second tension and a second cable bending stiffness; Determine whether the convergence condition is met based on the first tension and the second tension. If so, use the second tension and the second cable bending stiffness as the corresponding initial cable bending stiffness data and initial tension data, or use the second tension and the second cable bending stiffness as the corresponding cable bending stiffness data and tension data; if not, use the second tension and the second cable bending stiffness as the first tension and the first cable bending stiffness to update the frequency angle, tension and cable bending stiffness until the updated tension meets the convergence condition and uses the updated tension and cable bending stiffness as the corresponding initial cable bending stiffness data and initial tension data, or uses the updated tension and cable bending stiffness as the corresponding cable bending stiffness data and tension data.
6. The device for calculating the cable force and length of a cable rod synchronously according to claim 5, characterized in that: The first vibration least squares fitting formula is: ; The second vibration least squares fitting formula is: ; The frequency angle formula is: ; Where T is the tension, EI is the cable bending stiffness, EI0 is the first cable bending stiffness, L is the initial assumed length data or length update data, ρ is the cable linear density, i and f i are the frequency order and the frequency data of the i-th order, T0 is the first tension, EI k is the second cable bending stiffness or the cable bending stiffness calculated for the kth update, T k is the second tension or the tension calculated for the kth update, is the frequency angle of the i-th order.
7. The device for calculating the cable force and length of a cable rod synchronously according to claim 5, characterized in that: The second iterative calculation module is further used to calculate using a length update formula according to the kth initial assumed length data, the k+1th initial assumed length data, the cable bending stiffness design value and the corresponding initial cable bending stiffness data to obtain the k+2th length update data; The length update formula is: L k+2 =((L k+1 -L k ) / (EI k+1 -EI k ))EI+(L k -EI k (L k+1 -L k ) / (EI k+1 -EI k )), where L k+2 Update data for the k+2th length, L k+1 is the k+1th initial assumption length data, L k is the kth initial hypothesis length data, EI k+1 is the k+1th initial cable bending stiffness data, EI k is the kth initial cable bending stiffness data, and EI is the design value of the cable bending stiffness.
8. A terminal device, characterized in that: including a processor and a memory; The memory is used to store program codes and transmit the program codes to the processor; The processor is used to execute the calculation method for synchronously solving the cable rod tension and length as described in any one of claims 1 to 4 according to the instructions in the program code.
Citation Information
Patent Citations
Bridge inhaul cable bending rigidity recognition method
CN106096178A
Method, system and device for measuring cable force of existing inhaul cable of building structure
CN117928807A