Large-span bridge pneumatic admittance identification method and device based on operation measured data

Abstract

The invention provides a long-span bridge pneumatic admittance identification method and device based on operation measured data, and relates to the technical field of bridge safety, and the method comprises the steps: building a generalized motion equation of a long-span bridge under the action of buffeting force; establishing a bridge finite element model, and updating the bridge finite element model through operation measured data to obtain an optimal bridge finite element model; calculating a fluctuating wind coherence function of the large-span bridge girder; obtaining a three-component coefficient, a three-component force and a generalized buffeting force of the long-span bridge at different wind attack angles based on the optimal bridge finite element model; performing frequency domain analysis on the generalized motion equation to obtain simulated power spectral density of each point of the large-span bridge girder; and based on the simulated power spectrum density, identifying pneumatic admittance parameters in a preset pneumatic admittance function, and completing pneumatic admittance identification of the long-span bridge. According to the method, the problem that the pneumatic admittance of the long-span bridge in the operation period cannot be identified by an existing method is solved.

Description

Technical Field

[0001] The present invention relates to the field of bridge safety technology, and in particular to a method and device for identifying the aerodynamic admittance of a long-span bridge based on operational measured data. Background Art

[0002] Buffeting refers to a random vibration caused by fluctuating wind on engineering structures. Buffeting is inevitable for bridges exposed to natural fluctuating wind fields. Although buffeting usually does not cause catastrophic damage, large-scale buffeting responses can have a negative impact on driving safety, comfort, and fatigue performance of structural components. In the buffeting analysis of long-span bridges, the aerodynamic admittance function plays an important role as a transfer function that links the aerodynamic force acting on the bridge section with the fluctuating wind field of the incoming flow. Therefore, the identification of the aerodynamic admittance function has become an important part of the buffeting analysis of long-span bridges.

[0003] At present, the identification of aerodynamic admittance function is mostly determined by the pulsating wind speed and pulsating force obtained by wind tunnel test measurement. However, the wind field simulated by the wind tunnel often does not have the large-scale eddy current of the actual wind field, which often leads to low aerodynamic admittance values ​​in the low-frequency range. At present, more and more long-span bridges are equipped with long-term health monitoring systems that can measure wind speed and the flutter response of the bridge on site. Therefore, some scholars have identified the aerodynamic admittance by measuring the wind speed and the pulsating force obtained by the on-site pressure measurement method. However, if the bridge is in the operation stage, the on-site pressure measurement method is difficult to implement, which will make it impossible to identify the aerodynamic admittance of the large-span bridge during operation. Summary of the invention

[0004] The purpose of the present invention is to provide a method and device for identifying the aerodynamic admittance of a long-span bridge based on operational measured data, so as to improve the above-mentioned problem. In order to achieve the above-mentioned purpose, the technical solution adopted by the present invention is as follows:

[0005] In a first aspect, the present application provides a method for identifying aerodynamic admittance of a long-span bridge based on operational measured data, comprising:

[0006] Based on the parameters of the long-span bridge, the generalized motion equation of the long-span bridge under the action of buffeting force is established;

[0007] Establishing a bridge finite element model of the long-span bridge, updating the bridge finite element model through operational measured data of the long-span bridge, and obtaining an optimal bridge finite element model;

[0008] Based on the operational measured data, the pulsating wind coherence function of the main beam of the long-span bridge is calculated;

[0009] Obtaining the three-force coefficients of the long-span bridge at different wind attack angles based on the optimal bridge finite element model, and calculating the three-force and generalized buffeting forces of the main beam of the long-span bridge based on the three-force coefficients;

[0010] Based on the pulsating wind coherence function, the three-part force and the generalized buffeting force, the generalized motion equation is analyzed in the frequency domain to obtain the simulated power spectrum density of each point of the main beam of the long-span bridge;

[0011] Based on the simulated power spectrum density, the aerodynamic admittance parameters in the preset aerodynamic admittance function are identified, the aerodynamic admittance function is obtained, and the aerodynamic admittance identification of the long-span bridge is completed.

[0012] In the second aspect, the present application also provides a long-span bridge aerodynamic admittance identification device based on operational measured data, comprising:

[0013] The first building block is used to establish the generalized motion equation of the long-span bridge under the buffeting force based on the parameters of the long-span bridge;

[0014] The second construction module is used to establish a bridge finite element model of the long-span bridge, update the bridge finite element model through the operational measured data of the long-span bridge, and obtain an optimal bridge finite element model;

[0015] A first calculation module is used to calculate the pulsating wind coherence function of the main beam of the long-span bridge based on the operational measured data;

[0016] A second calculation module is used to obtain the three-force coefficients of the long-span bridge at different wind attack angles based on the optimal bridge finite element model, and calculate the three-force and generalized buffeting force of the main beam of the long-span bridge based on the three-force coefficients;

[0017] A third calculation module is used to perform frequency domain analysis on the generalized motion equation based on the pulsating wind coherence function, the three-part force and the generalized buffeting force to obtain a simulated power spectrum density of each point of the main beam of the long-span bridge;

[0018] The identification module is used to identify the aerodynamic admittance parameters in the preset aerodynamic admittance function based on the simulated power spectrum density, obtain the aerodynamic admittance function, and complete the aerodynamic admittance identification of the long-span bridge.

[0019] The beneficial effects of the present invention are:

[0020] The present invention collects data of the bridge during operation to obtain measured operation data, and identifies the pulsating wind coherence function of the long-span bridge, the three-force coefficient of the main beam and the generalized buffeting force through the measured operation data, which can indicate the effect of the buffeting force during the operation of the long-span bridge. At the same time, the aerodynamic admittance function of the long-span bridge in different directions during operation can be accurately identified through the optimization algorithm, which is convenient for the later buffeting analysis of the long-span bridge based on the aerodynamic admittance function identified by the present invention.

[0021] Other features and advantages of the present invention will be set forth in the following description, and in part will become apparent from the description, or may be understood by practicing embodiments of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for use in the embodiments are briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present invention and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without creative work.

[0023] Figure 1 It is a schematic flow chart of a method for identifying aerodynamic admittance of a long-span bridge based on operational measured data described in an embodiment of the present invention;

[0024] Figure 2 It is a schematic diagram of the structure of the aerodynamic admittance identification equipment for a large-span bridge based on operational measured data described in an embodiment of the present invention.

[0025] Markings in the figure: 800, aerodynamic admittance identification equipment for large-span bridges based on operational measured data; 801, processor; 802, memory; 803, multimedia component; 804, I / O interface; 805, communication component. DETAILED DESCRIPTION

[0026] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, rather than all of the embodiments. The components of the embodiments of the present invention generally described and shown in the drawings here can be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present invention provided in the drawings is not intended to limit the scope of the claimed invention, but merely represents selected embodiments of the present invention. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0027] It should be noted that similar reference numerals and letters represent similar items in the following drawings, so once an item is defined in one drawing, it does not need to be further defined and explained in the subsequent drawings. At the same time, in the description of the present invention, the terms "first", "second", etc. are only used to distinguish the description and cannot be understood as indicating or implying relative importance.

[0028] Embodiment 1:

[0029] This embodiment provides a method for identifying the aerodynamic admittance of a long-span bridge based on actual operational measured data.

[0030] See also Figure 1 , the figure shows that the method includes step S100, step S200, step S300, step S400, step S500, and step S600.

[0031] Step S100: establishing a generalized motion equation of the long-span bridge under buffeting force based on the long-span bridge parameters;

[0032] In this embodiment, the expression of the generalized motion equation is:

[0033]

[0034]

[0035] In the formula, I i represents the generalized mass of the i-th mode, K i represents the reduced frequency of the i-th mode, ω i represents the angular frequency of the i-th mode, D represents the main beam height, U represents the average wind speed, ζ i represents the damping ratio of the i-th mode, ξ i (·), and They represent the generalized displacement, generalized velocity and generalized acceleration of the i-th mode, B represents the width of the main beam, Q i (·) represents the generalized buffeting force of the i-th mode, and s represents the dimensionless time.

[0036] In this embodiment, the damping ratio can be obtained by processing the measured buffeting acceleration time history data using a fast Bayesian FFT method, wherein the measured buffeting acceleration time history data is operational measured data of a large-span bridge.

[0037] Step S200: establishing a bridge finite element model of the long-span bridge, updating the bridge finite element model through operational measured data of the long-span bridge, and obtaining an optimal bridge finite element model;

[0038] In this embodiment, measuring points are arranged at one-quarter span, one-half span and three-quarter span of the long-span bridge, and multiple sensors are arranged at each measuring point to collect actual operational data of the long-span bridge during operation.

[0039] The step S200 specifically includes:

[0040] Step S201: Based on the parameters of the long-span bridge, a three-dimensional beam unit is selected to simulate the main beam and the bridge tower of the long-span bridge, a spatial rod unit is selected to simulate the cable and the hanger of the long-span bridge, and a mass unit is used to simulate the auxiliary structure and the second-phase dead load of the long-span bridge to obtain a finite element model of the bridge;

[0041] Step S202: collecting operational measured data of the long-span bridge, wherein the operational measured data includes the measured vertical natural frequency, measured lateral natural frequency, measured torsional natural frequency, measured displacement data, measured environmental parameters and measured power spectrum density of the long-span bridge;

[0042] Step S203: defining a first objective function based on the operational measured data, and selecting the main beam material density, main beam elastic modulus, main cable material density and main cable initial stress in the bridge finite element model as update parameters;

[0043] In this embodiment, the first objective function is constructed using measured natural frequency data, wherein the measured natural frequency data includes the measured vertical natural frequency, the measured transverse natural frequency and the measured torsional natural frequency of the bridge. The expression of the first objective function is:

[0044]

[0045] Where OBJ Bridge represents the first objective function, f i sim represents the natural frequency of the i-th mode of the updated bridge finite element model, f i mea represents the natural frequency of the i-th mode of the measured frequency data, and n represents the total number of natural frequencies.

[0046] Step S204: updating the update parameters in the bridge finite element model by using a particle swarm algorithm until the first objective function converges to obtain an optimal bridge finite element model.

[0047] In this embodiment, a particle swarm algorithm is used to minimize the first objective function, wherein the performance of the particles is associated with the first objective function, and the positions of the particles correspond to the update parameters.

[0048] Step S300: Calculating the fluctuating wind coherence function of the main beam of the long-span bridge based on the operational measured data;

[0049] In this embodiment, based on the conditional simulation method, an iterative approach is adopted to obtain the pulsating wind coherence function of the entire long-span bridge main beam.

[0050] The Davenport coherence function is selected, and the attenuation constant in the Davenport coherence function is fitted based on the wind speed time history data of the measuring point to obtain a preliminary attenuation constant, and a preliminary coherence function is obtained through the preliminary attenuation constant.

[0051] Based on the preliminary coherence function and the wind spectrum information of the measuring point, the main beam wind field of the long-span bridge is conditionally simulated, and then the attenuation constant in the Davenport coherence function is fitted based on the conditionally simulated main beam wind field to obtain the updated attenuation constant. Repeat the above steps to iteratively determine the optimal attenuation constant and obtain the pulsating wind coherence function of the entire main beam, where the wind speed time history data and wind spectrum information are both operational measured data during the operation of the long-span bridge.

[0052] Step S400: obtaining the three-force coefficients of the long-span bridge at different wind attack angles based on the optimal bridge finite element model, and calculating the three-force and generalized buffeting forces of the main beam of the long-span bridge based on the three-force coefficients;

[0053] In this embodiment, the three-force coefficient of the main beam is determined by finite element analysis and on-site measured displacement data. The difference between the measured average displacement response of the main beam and the average displacement response of the bridge deck is used as the second objective function. The three-force coefficient is used as the identification parameter, and an optimization algorithm is used to identify the three-force coefficient under different wind attack angles.

[0054] And when the average wind speed is low, the flutter derivative of the bridge main beam will be very small, and the corresponding self-excited force will have a negligible effect on the bridge buffeting. Therefore, in order to eliminate the effect of the self-excited force on the identification of the aerodynamic admittance function, this embodiment will identify the aerodynamic admittance function at low wind speed.

[0055] The step S400 specifically includes:

[0056] Step S401: Calculating the measured average displacement response of the main beam based on the measured displacement data;

[0057] Step S402: calculating the average displacement response of the bridge deck based on the optimal bridge finite element model;

[0058] Step S403: defining a second objective function based on the measured average displacement response of the main beam and the average displacement response of the bridge deck;

[0059] In this embodiment, the difference between the measured average displacement response of the main beam and the average displacement response of the bridge deck is used as the second objective function.

[0060] Step S404: Using the three-force coefficient as the identification parameter of the optimal bridge finite element model, optimizing the identification parameters in the optimal bridge finite element model through the optimization algorithm until the second objective function converges, and obtaining the three-force coefficient under different wind attack angles, wherein the three-force coefficient includes a lift coefficient, a drag coefficient and a moment coefficient.

[0061] The step S400 further includes:

[0062] Step S405: Calculating the slope of the three-force coefficient based on the three-force coefficient;

[0063] Step S406: Calculating the three-force of the main beam of the long-span bridge based on the measured environmental parameters, the three-force coefficient and the slope of the three-force coefficient, wherein the measured environmental parameters include the average wind speed, the fluctuating wind speed in the downwind direction and the fluctuating wind speed in the vertical direction;

[0064] In this embodiment, the three forces of the main beam of the long-span bridge include lift, resistance and moment. The calculation formulas of the lift, resistance and moment are:

[0065]

[0066]

[0067]

[0068] Where, L b , D b and M b They represent lift, drag and torque respectively, C L , C D and C M They represent lift coefficient, drag coefficient and moment coefficient respectively, ρ represents air density, U represents average wind speed, C′ L , C′ D and C′ M They represent the slope of lift coefficient, drag coefficient and moment coefficient respectively, B represents the main beam width, u′ represents the fluctuating wind speed along the wind, w′ represents the fluctuating wind speed in the vertical direction, χ L , χ D and χ M They represent the aerodynamic admittance functions of the lift, drag, and moment to be identified respectively.

[0069] Step S407: Obtain the length of the long-span bridge and the vertical mode, lateral mode and torsional mode of the main beam of the long-span bridge;

[0070] Step S408: Calculate the generalized buffeting force of the main beam of the long-span bridge based on the length of the long-span bridge, the vertical mode, the lateral mode, the torsional direction mode and the three-part force.

[0071] In this embodiment, the calculation formula of the generalized buffeting force of the main beam of the long-span bridge is:

[0072]

[0073] In the formula, Q i (·) represents the generalized buffeting force of the i-th mode, s represents dimensionless time, l represents the length of the long-span bridge, and h i (·), p i (·) and α i (·) represents the i-th vertical mode, transverse mode and torsional mode of the main beam, x represents the longitudinal coordinate along the length of the main beam, L b (·),D b (·) and M b (·) represent the lift, drag and moment on the main beam, respectively, and B represents the width of the main beam.

[0074] Step S500: Based on the pulsating wind coherence function, the three-component force and the generalized buffeting force, the generalized motion equation is analyzed in the frequency domain to obtain the simulated power spectrum density of each point of the main beam of the long-span bridge;

[0075] The step S500 specifically includes:

[0076] Step S501: Substituting the generalized buffeting force into the generalized motion equation and performing Fourier transform to obtain a frequency domain equation;

[0077] In this embodiment, the generalized motion equation after substituting the generalized buffeting force is subjected to Fourier transform, and the expression of the frequency domain equation is obtained as follows:

[0078]

[0079] In the formula, K i represents the reduced frequency of the i-th order mode, K represents the reduced frequency, ζ i represents the damping ratio of the i-th mode, ρ represents the air density, B represents the width of the main beam, I i represents the generalized mass of the i-th mode, h i (·), p i (·) and α i (·) represents the i-th vertical mode, transverse mode and torsional mode of the main beam, x represents the longitudinal coordinate along the length of the main beam, and They represent the Fourier transform of lift, drag and moment on the main beam, Represents the Fourier transform of the generalized displacement of the ith mode.

[0080] Step S502: Calculate the first power spectrum density of the generalized displacement according to the random vibration theory and the frequency domain equation;

[0081] In this embodiment, the calculation formula of the first power spectrum density is:

[0082]

[0083] In the formula, represents the first power spectral density of the generalized displacement of the i-th mode, K represents the reduction frequency, K i represents the reduced frequency of the i-th order mode, ρ represents the air density, B represents the width of the main beam, I i represents the generalized mass of the i-th mode, ζ i represents the damping ratio of the i-th mode, l represents the length of the long-span bridge, S F (x A ,x B ,K) represents the generalized buffeting force at point x under the reduced frequency K A and point x B The cross-spectral density of .

[0084] Step S503: calculating three-dimensional generalized displacement based on the pulsating wind coherence function, the three-dimensional force and the first power spectrum density, wherein the three-dimensional generalized displacement includes generalized displacement in vertical, lateral and torsional directions;

[0085] In this embodiment, the expression of the three-way generalized displacement is:

[0086]

[0087]

[0088]

[0089]

[0090] r=h,p,α

[0091] a=u,w

[0092] In the formula, and They represent the generalized displacements in the vertical, lateral and torsional directions of the i-th mode, K represents the reduced frequency, K i represents the reduced frequency of the i-th order mode, ρ represents the air density, B represents the width of the main beam, I i represents the generalized mass of the i-th mode, ζ i represents the damping ratio of the i-th mode, χ L , χ D and χ MThey represent the aerodynamic admittance functions of the lift, drag and moment to be identified, respectively, and C L , C D and C M They represent lift coefficient, drag coefficient and moment coefficient respectively, C′ L , C′ D and C′ M They represent the slope of lift coefficient, drag coefficient and moment coefficient respectively, S uu (K) and S ww (K) represents the power spectrum density of the fluctuating wind in the downwind direction and vertical direction w under the reduced frequency K, represents the intermediate vector in direction a and direction r at the reduced frequency K, a and r both represent direction parameters, h represents vertical, p represents horizontal, α represents torsional direction, u represents the downwind direction of pulsating wind, w represents the vertical direction of pulsating wind, coh a (·) represents the fluctuating wind coherence function, l represents the length of the long-span bridge, h i (·), p i (·) and α i (·) represents the i-th vertical mode, lateral mode and torsional mode of the main beam, respectively, and |·| represents the absolute value.

[0093] Step S504: Based on the three-dimensional generalized displacement and the generalized motion equation, the simulated power spectrum density of each point on the main beam of the long-span bridge is calculated.

[0094] In this embodiment, the expression of the simulated power spectrum density of each point of the main beam of the long-span bridge is:

[0095]

[0096]

[0097]

[0098] In the formula, S h (·), S p (·) and S α (·) respectively represent the vertical simulated power spectrum density, lateral simulated power spectrum density and torsional simulated power spectrum density of each point on the main beam of the long-span bridge, x represents the longitudinal coordinate along the length of the main beam, K represents the reduced frequency, h i (·), p i (·) and α i (·) represents the i-th vertical mode, transverse mode and torsional mode of the main beam, B represents the width of the main beam, and They represent the generalized displacements in the vertical, lateral and torsional directions of the i-th mode, respectively.

[0099] Step S600: Based on the simulated power spectrum density, the aerodynamic admittance parameters in the preset aerodynamic admittance function are identified to obtain the aerodynamic admittance function, thereby completing the aerodynamic admittance identification of the long-span bridge.

[0100] The step S600 specifically includes:

[0101] Step S601: defining an aerodynamic admittance objective function based on the measured power spectral density and the simulated power spectral density;

[0102] In this embodiment, the aerodynamic admittance objective function focuses on the difference between the measured power spectrum density and the simulated power spectrum density at the peak point corresponding to the main vibration mode at the measuring point, wherein the expression of the aerodynamic admittance objective function is:

[0103]

[0104] Where OBJ r represents the aerodynamic admittance objective function in the r direction, x j represents the jth measured point, K k represents the reduced frequency of the kth concern, represents the simulated power spectral density corresponding to the kth reduced frequency of interest at the jth measured point in the r direction, It represents the measured power spectral density corresponding to the kth reduced frequency of interest at the jth measured point in the r direction, where r represents the direction parameter, r=h,p,α, h represents the vertical direction, p represents the lateral direction, and α represents the torsional direction.

[0105] Step S602: optimizing the parameters to be fitted in the preset aerodynamic admittance function by using a genetic algorithm until the aerodynamic admittance objective function converges to obtain aerodynamic admittance parameters;

[0106] In this embodiment, the preset aerodynamic admittance function is represented by the Larose model, wherein the expression of the preset aerodynamic admittance function is:

[0107]

[0108] In the formula, χ Lu (·) represents the aerodynamic admittance function in the downwind direction of the lift pulsating wind, χ Lw (·) represents the aerodynamic admittance function in the vertical direction of the lift pulsating wind, χ Du (·) represents the aerodynamic admittance function in the downwind direction of the resistance pulsating wind, χ Dw (·) represents the aerodynamic admittance function in the vertical direction of the resistance pulsating wind, χ Mu (·) represents the aerodynamic admittance function in the downwind direction of the moment pulsating wind, χ Mw(·) represents the aerodynamic admittance function in the vertical direction of the moment pulsating wind, a represents the direction parameter, f represents the force parameter, K represents the reduced frequency, γ and β represent the parameters to be fitted, and |·| represents the absolute value.

[0109] Step S603: Substitute the aerodynamic admittance parameter into the preset aerodynamic admittance function to obtain the aerodynamic admittance function of the long-span bridge.

[0110] In summary, the present invention can accurately identify the aerodynamic admittance function of the main beam of a long-span bridge in the vertical, lateral and torsional directions based on the measured data. At the same time, based on the aerodynamic admittance function identified by the present invention, a buffeting analysis of a long-span bridge is performed, and the power spectrum density of the measured point response in the numerical simulation matches well with the power spectrum density of the measured response.

[0111] Embodiment 2:

[0112] This embodiment provides a device for identifying aerodynamic admittance of a long-span bridge based on operational measured data, the device comprising:

[0113] The first building block is used to establish the generalized motion equation of the long-span bridge under the buffeting force based on the parameters of the long-span bridge;

[0114] The second construction module is used to establish a bridge finite element model of the long-span bridge, update the bridge finite element model through the operational measured data of the long-span bridge, and obtain an optimal bridge finite element model;

[0115] A first calculation module is used to calculate the pulsating wind coherence function of the main beam of the long-span bridge based on the operational measured data;

[0116] A second calculation module is used to obtain the three-force coefficients of the long-span bridge at different wind attack angles based on the optimal bridge finite element model, and calculate the three-force and generalized buffeting force of the main beam of the long-span bridge based on the three-force coefficients;

[0117] A third calculation module is used to perform frequency domain analysis on the generalized motion equation based on the pulsating wind coherence function, the three-part force and the generalized buffeting force to obtain a simulated power spectrum density of each point of the main beam of the long-span bridge;

[0118] The identification module is used to identify the aerodynamic admittance parameters in the preset aerodynamic admittance function based on the simulated power spectrum density, obtain the aerodynamic admittance function, and complete the aerodynamic admittance identification of the long-span bridge.

[0119] The second building block comprises:

[0120] A construction unit is used to select a three-dimensional beam unit to simulate the main beam and the bridge tower of the long-span bridge based on the parameters of the long-span bridge, select a spatial rod unit to simulate the cable and the hanger of the long-span bridge, and simulate the auxiliary structure and the second-phase dead load of the long-span bridge through the mass unit to obtain a finite element model of the bridge;

[0121] A data acquisition unit, used to collect operational measured data of the long-span bridge, wherein the operational measured data includes the measured vertical natural frequency, measured lateral natural frequency, measured torsional natural frequency, measured displacement data, measured environmental parameters and measured power spectrum density of the long-span bridge;

[0122] A first definition unit is used to define a first objective function based on the operational measured data, and select the main beam material density, the main beam elastic modulus, the main cable material density and the main cable initial stress in the bridge finite element model as update parameters;

[0123] The first updating unit is used to update the update parameters in the bridge finite element model by using a particle swarm algorithm until the first objective function converges to obtain an optimal bridge finite element model.

[0124] The second calculation module includes:

[0125] A first calculation unit is used to calculate the measured average displacement response of the main beam based on the measured displacement data;

[0126] A second calculation unit is used to calculate the average displacement response of the bridge deck based on the optimal bridge finite element model;

[0127] A second definition unit is used to define a second objective function based on the measured average displacement response of the main beam and the average displacement response of the bridge deck;

[0128] A second updating unit is used to use the three-force coefficient as an identification parameter of the optimal bridge finite element model, optimize the identification parameter in the optimal bridge finite element model by an optimization algorithm, until the second objective function converges, and obtain the three-force coefficient under different wind attack angles, wherein the three-force coefficient includes a lift coefficient, a drag coefficient, and a moment coefficient;

[0129] A third calculation unit, configured to calculate a slope of a three-force coefficient based on the three-force coefficient;

[0130] A fourth calculation unit is used to calculate the three-force of the main beam of the long-span bridge based on the measured environmental parameters, the three-force coefficient and the slope of the three-force coefficient, wherein the measured environmental parameters include the average wind speed, the fluctuating wind speed in the downwind direction and the fluctuating wind speed in the vertical direction;

[0131] A data acquisition unit, used to acquire the length of the long-span bridge and the vertical mode, lateral mode and torsional mode of the main beam of the long-span bridge;

[0132] The fifth calculation unit is used to calculate the generalized flutter force of the main beam of the long-span bridge based on the length of the long-span bridge, the vertical mode, the lateral mode, the torsional direction mode and the three-part force.

[0133] The third computing module comprises:

[0134] a sixth calculation unit, configured to substitute the generalized buffeting force into the generalized motion equation and perform Fourier transform to obtain a frequency domain equation;

[0135] a seventh calculation unit, used for calculating a first power spectrum density of generalized displacement according to random vibration theory and the frequency domain equation;

[0136] an eighth calculation unit, configured to calculate a three-dimensional generalized displacement based on the pulsating wind coherence function, the three-component force and the first power spectrum density, wherein the three-dimensional generalized displacement includes generalized displacements in vertical, lateral and torsional directions;

[0137] The ninth calculation unit is used to calculate the simulated power spectrum density of each point on the main beam of the long-span bridge based on the three-dimensional generalized displacement and the generalized motion equation.

[0138] In this embodiment, the identification module includes:

[0139] A third definition unit, configured to define an aerodynamic admittance objective function based on the measured power spectral density and the simulated power spectral density;

[0140] A third updating unit is used to optimize the parameters to be fitted in the preset aerodynamic admittance function by using a genetic algorithm until the aerodynamic admittance target function converges to obtain the aerodynamic admittance parameters;

[0141] The identification unit is used to substitute the aerodynamic admittance parameter into the preset aerodynamic admittance function to obtain the aerodynamic admittance function of the long-span bridge.

[0142] It should be noted that, regarding the device in the above embodiment, the specific manner in which each module performs the operation has been described in detail in the embodiment of the method, and will not be elaborated here.

[0143] Embodiment 3:

[0144] Corresponding to the above method embodiment, this embodiment also provides a large-span bridge aerodynamic admittance identification device based on operational measured data. The large-span bridge aerodynamic admittance identification device based on operational measured data described below and the large-span bridge aerodynamic admittance identification method based on operational measured data described above can be referenced to each other.

[0145] Figure 2FIG. 8 is a block diagram of a long-span bridge aerodynamic admittance identification device 800 based on operational measured data according to an exemplary embodiment. Figure 2 As shown, the long-span bridge aerodynamic admittance identification device 800 based on operational measured data may include: a processor 801, a memory 802. The long-span bridge aerodynamic admittance identification device 800 based on operational measured data may also include one or more of a multimedia component 803, an I / O interface 804, and a communication component 805.

[0146] The processor 801 is used to control the overall operation of the large-span bridge pneumatic admittance identification device 800 based on operational measured data, so as to complete all or part of the steps in the large-span bridge pneumatic admittance identification method based on operational measured data. The memory 802 is used to store various types of data to support the operation of the large-span bridge pneumatic admittance identification device 800 based on operational measured data, and these data may include, for example, instructions for any application or method operated on the large-span bridge pneumatic admittance identification device 800 based on operational measured data, as well as application-related data, such as contact data, sent and received messages, pictures, audio, video, etc. The memory 802 can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, disk or optical disk. The multimedia component 803 may include a screen and an audio component. The screen may be, for example, a touch screen, and the audio component is used to output and / or input audio signals. For example, the audio component may include a microphone, which is used to receive external audio signals. The received audio signal may be further stored in the memory 802 or sent via the communication component 805. The audio component also includes at least one speaker for outputting audio signals. The I / O interface 804 provides an interface between the processor 801 and other interface modules, and the above-mentioned other interface modules can be keyboards, mice, buttons, etc. These buttons can be virtual buttons or physical buttons. The communication component 805 is used for wired or wireless communication between the large-span bridge aerodynamic admittance identification device 800 based on operational measured data and other devices. Wireless communication, such as Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G or 4G, or a combination of one or more of them, so the corresponding communication component 805 can include: Wi-Fi module, Bluetooth module, NFC module.

[0147] In an exemplary embodiment, the large-span bridge aerodynamic admittance identification device 800 based on operational measured data can be implemented by one or more application-specific integrated circuits (ASIC), digital signal processors (DSP), digital signal processing devices (DSPD), programmable logic devices (PLD), field programmable gate arrays (FPGA), controllers, microcontrollers, microprocessors or other electronic components to execute the above-mentioned large-span bridge aerodynamic admittance identification method based on operational measured data.

[0148] In another exemplary embodiment, a computer-readable storage medium including program instructions is also provided, and when the program instructions are executed by a processor, the steps of the above-mentioned method for identifying the aerodynamic admittance of a large-span bridge based on operational measured data are implemented. For example, the computer-readable storage medium may be the above-mentioned memory 802 including program instructions, and the above-mentioned program instructions may be executed by the processor 801 of the device 800 for identifying the aerodynamic admittance of a large-span bridge based on operational measured data to complete the above-mentioned method for identifying the aerodynamic admittance of a large-span bridge based on operational measured data.

[0149] Embodiment 4:

[0150] Corresponding to the above method embodiment, a readable storage medium is also provided in this embodiment. The readable storage medium described below and the aerodynamic admittance identification method of a large-span bridge based on operational measured data described above can refer to each other.

[0151] A readable storage medium stores a computer program, which, when executed by a processor, implements the steps of the aerodynamic admittance identification method for a large-span bridge based on operational measured data of the above method embodiment.

[0152] The readable storage medium may specifically be a U disk, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, or other readable storage medium that can store program codes.

[0153] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

[0154] The above description is only a specific implementation mode of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can easily think of changes or substitutions within the technical scope disclosed by the present invention, which should be covered by the protection scope of the present invention.

Claims

1. A method for identifying aerodynamic admittance of a long-span bridge based on operational measured data, characterized in that: include: Based on the parameters of the long-span bridge, the generalized motion equation of the long-span bridge under the action of buffeting force is established; Establishing a bridge finite element model of the long-span bridge, updating the bridge finite element model through operational measured data of the long-span bridge, and obtaining an optimal bridge finite element model; Based on the operational measured data, the pulsating wind coherence function of the main beam of the long-span bridge is calculated; Obtaining the three-force coefficients of the long-span bridge at different wind attack angles based on the optimal bridge finite element model, and calculating the three-force and generalized buffeting forces of the main beam of the long-span bridge based on the three-force coefficients; Based on the pulsating wind coherence function, the three-part force and the generalized buffeting force, the generalized motion equation is analyzed in the frequency domain to obtain the simulated power spectrum density of each point of the main beam of the long-span bridge; Based on the simulated power spectrum density, the aerodynamic admittance parameters in the preset aerodynamic admittance function are identified, the aerodynamic admittance function is obtained, and the aerodynamic admittance identification of the long-span bridge is completed.

2. The aerodynamic admittance identification method of a long-span bridge based on operational measured data according to claim 1 is characterized in that The step of establishing the bridge finite element model of the long-span bridge, updating the bridge finite element model through the operational measured data of the long-span bridge, and obtaining the optimal bridge finite element model comprises: Based on the parameters of the long-span bridge, a three-dimensional beam unit is selected to simulate the main beam and the bridge tower of the long-span bridge, a spatial rod unit is selected to simulate the cable and the hanger of the long-span bridge, and a mass unit is used to simulate the auxiliary structure and the second-phase dead load of the long-span bridge to obtain a finite element model of the bridge; Collecting operational measured data of the long-span bridge, wherein the operational measured data includes the measured vertical natural frequency, measured lateral natural frequency, measured torsional natural frequency, measured displacement data, measured environmental parameters and measured power spectrum density of the long-span bridge; A first objective function is defined based on the operational measured data, and the main beam material density, the main beam elastic modulus, the main cable material density and the main cable initial stress in the bridge finite element model are selected as update parameters; The update parameters in the bridge finite element model are updated by using a particle swarm algorithm until the first objective function converges, thereby obtaining an optimal bridge finite element model.

3. The aerodynamic admittance identification method of a long-span bridge based on operational measured data according to claim 2 is characterized in that The method of obtaining the three-force coefficient of the long-span bridge under different wind attack angles based on the optimal bridge finite element model includes: Based on the measured displacement data, calculating the measured average displacement response of the main beam; Based on the optimal bridge finite element model, calculating the average displacement response of the bridge deck; Defining a second objective function based on the measured average displacement response of the main beam and the average displacement response of the bridge deck; The three-force coefficient is used as the identification parameter of the optimal bridge finite element model. The identification parameters in the optimal bridge finite element model are optimized by an optimization algorithm until the second objective function converges, and the three-force coefficient under different wind attack angles is obtained. The three-force coefficient includes a lift coefficient, a drag coefficient and a moment coefficient.

4. The method for identifying aerodynamic admittance of a long-span bridge based on operational measured data according to claim 2 is characterized in that , the calculation of the three-force and generalized buffeting force of the main beam of the long-span bridge based on the three-force coefficient includes: Based on the three-force coefficient, calculating the three-force coefficient slope; Calculating the three-force of the main beam of the long-span bridge based on the measured environmental parameters, the three-force coefficient and the slope of the three-force coefficient, wherein the measured environmental parameters include the average wind speed, the fluctuating wind speed in the downwind direction and the fluctuating wind speed in the vertical direction; Obtain the length of the long-span bridge and the vertical mode, lateral mode and torsional mode of the main beam of the long-span bridge; The generalized buffeting force of the main beam of the long-span bridge is calculated based on the length of the long-span bridge, the vertical mode, the lateral mode, the torsional direction mode and the three-part force.

5. The method for identifying aerodynamic admittance of a long-span bridge based on operational measured data according to claim 1 is characterized in that Based on the pulsating wind coherence function, the three-part force and the generalized buffeting force, the generalized motion equation is analyzed in the frequency domain to obtain the simulated power spectrum density of each point of the main beam of the long-span bridge, including: Substituting the generalized buffeting force into the generalized motion equation and performing Fourier transformation to obtain a frequency domain equation; According to the random vibration theory and the frequency domain equation, the first power spectrum density of the generalized displacement is calculated; Based on the pulsating wind coherence function, the three-part force and the first power spectrum density, calculating the three-dimensional generalized displacement, the three-dimensional generalized displacement including the generalized displacement in the vertical, lateral and torsional directions; Based on the three-dimensional generalized displacement and the generalized motion equation, the simulated power spectrum density of each point on the main beam of the long-span bridge is calculated.

6. The method for identifying aerodynamic admittance of a long-span bridge based on operational measured data according to claim 2 is characterized in that , based on the simulated power spectrum density, identifying the aerodynamic admittance parameters in the preset aerodynamic admittance function to obtain the aerodynamic admittance function, including: Based on the measured power spectral density and the simulated power spectral density, defining an aerodynamic admittance objective function; Optimizing the parameters to be fitted in the preset aerodynamic admittance function by using a genetic algorithm until the aerodynamic admittance objective function converges, thereby obtaining the aerodynamic admittance parameters; The aerodynamic admittance parameter is substituted into the preset aerodynamic admittance function to obtain the aerodynamic admittance function of the long-span bridge.

7. A device for identifying aerodynamic admittance of a long-span bridge based on operational measured data, characterized in that: include: The first building block is used to establish the generalized motion equation of the long-span bridge under the buffeting force based on the parameters of the long-span bridge; The second construction module is used to establish a bridge finite element model of the long-span bridge, update the bridge finite element model through the operational measured data of the long-span bridge, and obtain an optimal bridge finite element model; A first calculation module is used to calculate the pulsating wind coherence function of the main beam of the long-span bridge based on the operational measured data; A second calculation module is used to obtain the three-force coefficients of the long-span bridge at different wind attack angles based on the optimal bridge finite element model, and calculate the three-force and generalized buffeting force of the main beam of the long-span bridge based on the three-force coefficients; A third calculation module is used to perform frequency domain analysis on the generalized motion equation based on the pulsating wind coherence function, the three-part force and the generalized buffeting force to obtain a simulated power spectrum density of each point of the main beam of the long-span bridge; The identification module is used to identify the aerodynamic admittance parameters in the preset aerodynamic admittance function based on the simulated power spectrum density, obtain the aerodynamic admittance function, and complete the aerodynamic admittance identification of the long-span bridge.

8. The aerodynamic admittance identification device for a long-span bridge based on operational measured data according to claim 7 is characterized in that: The second building block comprises: A construction unit is used to select a three-dimensional beam unit to simulate the main beam and the bridge tower of the long-span bridge based on the parameters of the long-span bridge, select a spatial rod unit to simulate the cable and the hanger of the long-span bridge, and simulate the auxiliary structure and the second-phase dead load of the long-span bridge through the mass unit to obtain a finite element model of the bridge; A data acquisition unit, used to collect operational measured data of the long-span bridge, wherein the operational measured data includes the measured vertical natural frequency, measured lateral natural frequency, measured torsional natural frequency, measured displacement data, measured environmental parameters and measured power spectrum density of the long-span bridge; A first definition unit is used to define a first objective function based on the operational measured data, and select the main beam material density, the main beam elastic modulus, the main cable material density and the main cable initial stress in the bridge finite element model as update parameters; The first updating unit is used to update the update parameters in the bridge finite element model by using a particle swarm algorithm until the first objective function converges to obtain an optimal bridge finite element model.

9. The aerodynamic admittance identification device for a long-span bridge based on operational measured data according to claim 8 is characterized in that: The second calculation module includes: A first calculation unit is used to calculate the measured average displacement response of the main beam based on the measured displacement data; A second calculation unit is used to calculate the average displacement response of the bridge deck based on the optimal bridge finite element model; A second definition unit is used to define a second objective function based on the measured average displacement response of the main beam and the average displacement response of the bridge deck; A second updating unit is used to use the three-force coefficient as an identification parameter of the optimal bridge finite element model, optimize the identification parameter in the optimal bridge finite element model by an optimization algorithm, until the second objective function converges, and obtain the three-force coefficient under different wind attack angles, wherein the three-force coefficient includes a lift coefficient, a drag coefficient, and a moment coefficient; A third calculation unit, configured to calculate a slope of a three-force coefficient based on the three-force coefficient; A fourth calculation unit is used to calculate the three-force of the main beam of the long-span bridge based on the measured environmental parameters, the three-force coefficient and the slope of the three-force coefficient, wherein the measured environmental parameters include the average wind speed, the fluctuating wind speed in the downwind direction and the fluctuating wind speed in the vertical direction; A data acquisition unit, used to acquire the length of the long-span bridge and the vertical mode, lateral mode and torsional mode of the main beam of the long-span bridge; The fifth calculation unit is used to calculate the generalized flutter force of the main beam of the long-span bridge based on the length of the long-span bridge, the vertical mode, the lateral mode, the torsional direction mode and the three-part force.

10. The aerodynamic admittance identification device for a long-span bridge based on operational measured data according to claim 7, characterized in that: The third computing module comprises: a sixth calculation unit, configured to substitute the generalized buffeting force into the generalized motion equation and perform Fourier transform to obtain a frequency domain equation; A seventh calculation unit, used for calculating a first power spectrum density of generalized displacement according to random vibration theory and the frequency domain equation; an eighth calculation unit, configured to calculate a three-dimensional generalized displacement based on the pulsating wind coherence function, the three-component force and the first power spectrum density, wherein the three-dimensional generalized displacement includes generalized displacements in vertical, lateral and torsional directions; The ninth calculation unit is used to calculate the simulated power spectrum density of each point on the main beam of the long-span bridge based on the three-dimensional generalized displacement and the generalized motion equation.

Images