A vortex vibration response prediction method, device and equipment of a cantilever column structure and a medium

Abstract

This application relates to the field of vortex-induced vibration technology, and particularly to a method, apparatus, device, and medium for predicting the vortex-induced vibration response of a cantilever column structure. The method includes: dividing the target cantilever column structure into structural elements to obtain several series-connected structural elements; constructing a model of a constant cross-section structure for each structural element to obtain a constant cross-section structure, and establishing a finite element model to be tested based on each constant cross-section structure; determining whether all meet the requirements of a preset characteristic size evaluation system, a preset average wind load evaluation system, and a preset vortex-induced force evaluation system, until the finite element model to be tested is determined to be the target finite element model; obtaining the target average wind load and target vortex-induced force parameters; obtaining the overall average wind load and overall vortex-induced force parameters of the target cantilever column structure in a preset coordinate system; constructing a dynamic equilibrium equation, and predicting the vortex-induced vibration response based on the dynamic equilibrium equation. This application can improve the accuracy of vortex-induced vibration response prediction.

Description

Technical Field

[0001] This application relates to the field of vortex-induced vibration technology, and in particular to a method, apparatus, equipment and medium for predicting the vortex-induced vibration response of a cantilever column structure. Background Technology

[0002] Many engineering structures are cantilevered column structures with similar and gradually varying cross-sections, such as wind turbine towers, chimneys, bridge towers, and power poles. These structures have cross-sectional types such as circular, elliptical, and polygonal, with cross-sectional dimensions gradually changing as the structure's height increases—sometimes decreasing, sometimes increasing, sometimes decreasing. The towering shape of these cantilevered column structures presents new challenges to safety. The taller and more flexible the structure, the higher the probability of vortex-induced vibration (VEM). While VEM does not directly cause structural dynamic instability and failure, large-amplitude VEM can reduce assembly precision during construction, increase construction difficulty, endanger the safety of construction personnel, affect structural safety during operation, reduce fatigue life, and ultimately impact the normal use of the structure. Therefore, a predictive method for structural VEM is needed to improve the analytical accuracy and precision of VEM problems. Summary of the Invention

[0003] This application provides a method, apparatus, device, and medium for predicting the vortex-induced vibration response of a cantilever column structure, in order to solve one or more technical problems existing in the prior art, and at least provide a beneficial option or create conditions.

[0004] Other features and advantages of this application will become apparent from the following detailed description, or may be learned in part from practice of this application.

[0005] According to one aspect of the embodiments of this application, a method for predicting the vortex-induced vibration response of a cantilever column structure is provided, the method comprising:

[0006] S1, Divide the target cantilever column structure into structural units to obtain several connected structural units;

[0007] S2, construct a model of a constant cross-section structure for each of the structural units to obtain a constant cross-section structure corresponding to each of the structural units, and establish a finite element model to be tested corresponding to the target cantilever column structure based on each of the constant cross-section structures;

[0008] S3, determine whether the finite element models to be tested all meet the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system and the preset vortex-induced force evaluation system.

[0009] S4. If the finite element model to be tested does not meet the requirements of the characteristic size evaluation system, the average wind load evaluation system, or the vortex-induced force evaluation system, then repeat steps S1-S2 until the finite element model to be tested meets the requirements of the characteristic size evaluation system, the average wind load evaluation system, and the vortex-induced force evaluation system, and determine the finite element model to be tested as the target finite element model.

[0010] S5, obtain the structural characteristic parameters of the target cantilever column structure, and determine the target average wind load and the target vortex-induced force parameters of the target finite element model under the average wind load evaluation system based on the structural characteristic parameters.

[0011] S6. Based on the target average wind load and the target vortex-induced force parameters, the coordinate system transformation of the element transformation matrix of each target structural unit in the target finite element model is performed to obtain the overall average wind load and overall vortex-induced force parameters of the target cantilever column structure under the preset coordinate system.

[0012] S7. Determine the structural vibration acceleration, structural vibration velocity, and structural vibration displacement of the target cantilever column structure based on the structural characteristic parameters, and construct the dynamic equilibrium equation corresponding to the target cantilever column structure based on the structural vibration acceleration, structural vibration velocity, structural vibration displacement, overall average wind load, overall vortex-induced force parameters, and preset weak shear flow effect correction coefficient.

[0013] S8. Based on the dynamic equilibrium equation, predict the vortex-induced vibration response of the target cantilever column structure when subjected to the incoming wind.

[0014] In one embodiment of this application, based on the foregoing scheme, a model of a uniform cross-section structure is constructed for a single structural unit, including the following steps:

[0015] Obtain any target point in the height direction of the structural unit, and construct a cross-section structure with the same height as the structural unit using the cross-section where the target point is located;

[0016] In the constant cross-section structure, the shape of each cross section is consistent with the shape of the cross section where the target point is located.

[0017] In one embodiment of this application, based on the foregoing scheme, whether the finite element model to be tested meets the requirements of the preset feature size evaluation system is determined by the following steps:

[0018] Obtain the characteristic dimension parameters of each of the structural units, the characteristic dimension parameters including the upper end face cross-sectional dimension, the lower end face cross-sectional dimension and the target face cross-sectional dimension of the structural unit, the target face cross-sectional dimension corresponding to the dimension of the cross-section where the target point is located;

[0019] The characteristic dimension discrimination value of each structural unit is determined based on the cross-sectional dimensions of each upper end face, each lower end face, and each target surface.

[0020] If all the discriminant values ​​of the aforementioned feature dimensions are lower than the preset feature dimension threshold, then the finite element model to be tested is determined to meet the requirements of the preset feature dimension evaluation system.

[0021] The characteristic size discrimination value of each of the structural units satisfies the following formula:

[0022]

[0023] In the formula, The feature size discriminant value, The dimensions of the upper end face cross-section are as follows: The dimensions of the lower end face cross-section are as follows: The target surface cross-sectional dimensions are... , , .

[0024] In one embodiment of this application, based on the foregoing scheme, whether the finite element model to be tested meets the requirements of the average wind load evaluation system is determined by the following steps:

[0025] The average wind load parameters on the upper end face of each structural unit, the average wind load parameters on the lower end face of each structural unit, and the average wind load parameters on the target face of each structural unit are obtained.

[0026] The average wind load parameters of each of the upper end faces, the average wind load parameters of each of the lower end faces, and the average wind load parameters of each of the target faces are input into the preset average wind load calculation formula to obtain the first average wind load of each of the upper end faces, the second average wind load of each of the lower end faces, and the third average wind load of each of the target faces.

[0027] The average wind load discrimination value of each structural unit is determined based on each of the first average wind load, each of the second average wind load, and each of the third average wind load.

[0028] If all the average wind load discrimination values ​​are lower than the preset average wind load threshold, then the finite element model to be tested is determined to meet the requirements of the average wind load evaluation system.

[0029] The average wind load parameters include air density, drag coefficient, and wind speed at any height. The cross-sectional dimensions and arbitrary height of the structural unit at the location The average incoming air velocity;

[0030] any height The average wind load can be obtained using the following formula for calculating the average wind load:

[0031]

[0032] In the formula, The air density is mentioned. The drag coefficient is mentioned above; To be at any height The cross-sectional dimensions of the structural unit at that location. For the arbitrary height The average incoming air velocity; For any height The average wind load experienced at the location.

[0033] In one embodiment of this application, based on the foregoing scheme, the average wind load discrimination value can be obtained by the following formula:

[0034]

[0035] In the formula, The average wind load discrimination value is... The first average wind load, This is the second average wind load. The third average wind load; The average incoming wind speed at the height of the upper end face. The average incoming wind speed at the height of the lower end face. The average incoming wind speed at the height of the target surface.

[0036] In one embodiment of this application, based on the foregoing scheme, whether the finite element model to be tested meets the requirements of the vortex-induced force evaluation system is determined by the following steps:

[0037] Obtain the first vortex-induced force parameter of the upper end face of each structural unit, the second vortex-induced force parameter of the lower end face of each structural unit, and the third vortex-induced force parameter of the target surface of each structural unit.

[0038] The vortex-induced force discrimination value of each structural unit is determined based on each of the first vortex-induced force parameters, each of the second vortex-induced force parameters, and each of the third vortex-induced force parameters.

[0039] If all the vortex-induced force discrimination values ​​are lower than the preset vortex-induced force threshold, then the finite element model to be tested is determined to meet the requirements of the vortex-induced force evaluation system.

[0040] In one embodiment of this application, based on the foregoing scheme, the structural characteristic parameters include the average mass per meter, average damping per meter, and average stiffness per meter of the target cantilever column structure; the dynamic equilibrium equation is specifically the following formula:

[0041]

[0042] In the formula, The average mass per meter of the target cantilever column structure. The average damping per meter of the target cantilever column structure is given. The average stiffness per meter of the target cantilever column structure; The vibration acceleration of the structure, The vibration velocity of the structure is... The vibration displacement of the structure; The overall average wind load is... The overall vortex-induced force parameters are as follows: This is the correction coefficient for the eddy-induced weak shear flow effect.

[0043] According to one aspect of the embodiments of this application, a device for predicting the vortex-induced vibration response of a cantilever column structure is provided, the device comprising:

[0044] The unit division is used to divide the target cantilever column structure into structural units, resulting in several connected structural units;

[0045] The construction unit is used to construct a model of a constant cross-section structure for each of the structural units, thereby obtaining a constant cross-section structure that corresponds one-to-one with each of the structural units, and to establish a finite element model to be tested corresponding to the target cantilever column structure based on each of the constant cross-section structures.

[0046] The judgment unit is used to determine whether the finite element model to be tested meets the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system, and the preset vortex-induced force evaluation system.

[0047] An iterative unit is used to repeat the above steps if the finite element model to be tested does not meet the requirements of the characteristic size evaluation system, the average wind load evaluation system, or the vortex-induced force evaluation system, until the finite element model to be tested meets the requirements of the characteristic size evaluation system, the average wind load evaluation system, and the vortex-induced force evaluation system, and then determines the finite element model to be tested as the target finite element model.

[0048] The acquisition unit is used to acquire the structural characteristic parameters of the target cantilever column structure, and determine the target average wind load and the target vortex-induced force parameters of the target finite element model under the average wind load evaluation system based on the structural characteristic parameters.

[0049] The transformation unit is used to perform coordinate system transformation of the element transformation matrix of each target structural element in the target finite element model according to the target average wind load and the target vortex-induced force parameters, so as to obtain the overall average wind load and overall vortex-induced force parameters of the target cantilever column structure under the preset coordinate system.

[0050] The dynamic equilibrium unit is used to determine the structural vibration acceleration, structural vibration velocity, and structural vibration displacement of the target cantilever column structure based on the structural characteristic parameters, and to construct the dynamic equilibrium equation corresponding to the target cantilever column structure based on the structural vibration acceleration, structural vibration velocity, structural vibration displacement, overall average wind load, overall vortex-induced force parameters, and a preset weak shear flow effect correction coefficient.

[0051] The prediction unit is used to predict the vortex-induced vibration response of the target cantilever column structure when subjected to incoming wind, based on the dynamic equilibrium equation.

[0052] According to one aspect of the embodiments of this application, a computer-readable storage medium is provided that stores a computer program thereon, the computer program including executable instructions that, when executed by a processor, implement the method described in the above embodiments.

[0053] According to one aspect of the embodiments of this application, an electronic device is provided, including: one or more processors; and a memory for storing executable instructions of the processors, which, when executed by the one or more processors, cause the one or more processors to perform the method as described in the above embodiments.

[0054] The beneficial effects of this application are as follows: First, by dividing the target cantilever column structure into structural elements, several connected structural elements are obtained. The model of the divided structural elements with equal cross-sections can be constructed, and then a finite element model to be tested corresponding to the target cantilever column structure can be established.

[0055] By determining whether the finite element models to be tested all meet the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system, and the preset vortex-induced force evaluation system, the final target finite element models all meet the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system, and the preset vortex-induced force evaluation system, thereby equivalently simulating the vortex-induced vibration that may occur in the target cantilever column structure in the actual scenario.

[0056] Furthermore, by obtaining the structural characteristic parameters of the target cantilever column structure, that is, by obtaining the target average wind load and the target vortex-induced force parameters under the average wind load evaluation system and the target vortex-induced force evaluation system, the coordinate system transformation of the element transformation matrix of each target structural unit is performed to equivalently simulate the overall average wind load and overall vortex-induced force parameters of the target cantilever column structure under the preset coordinate system.

[0057] Furthermore, by determining the structural vibration acceleration, vibration velocity, and vibration displacement of the target cantilever column structure through structural characteristic parameters, and correcting them according to a preset weak shear flow effect correction coefficient, the influence of weak shear flow in the vertical direction of the cantilever column structure's cross-section due to its gradual change characteristics can be eliminated, thereby improving the accuracy and precision of the analysis of vortex-induced vibration problems. Simultaneously, the constructed dynamic equilibrium equations can be used to predict the vortex-induced vibration response of the target cantilever column structure under the influence of incoming wind. By monitoring the parameters of the incoming wind, it is possible to predict in real time whether the cantilever column structure will experience vortex-induced vibration response.

[0058] Therefore, the vortex-induced vibration response prediction method for cantilever column structures provided in this application can eliminate the influence of weak shear flow and establish dynamic equilibrium equations to accurately predict the vortex-induced vibration response. Attached Figure Description

[0059] To more clearly illustrate the technical solutions in the embodiments of this application, the accompanying drawings used in the description of the embodiments will be briefly explained below. Obviously, the described drawings are only a part of the embodiments of this application, and not all of them. Those skilled in the art can obtain other design schemes and drawings based on these drawings without creative effort.

[0060] Figure 1 This is a schematic diagram illustrating the target cantilever column structure under the influence of incoming wind, according to an embodiment of this application.

[0061] Figure 2 This is a flowchart illustrating a method for predicting the vortex-induced vibration response of a cantilever column structure according to an embodiment of this application;

[0062] Figure 3This is a schematic diagram illustrating a single structural unit (left side) and a single uniform cross-section structure (right side) according to embodiments of this application;

[0063] Figure 4 This is a logic flowchart illustrating the entire vortex-induced vibration response prediction according to an embodiment of this application;

[0064] Figure 5 A schematic diagram of the lumped mass additional nodes for establishing a finite element model of the cantilever column structure shown in the embodiments of this application;

[0065] Figure 6 This is a block diagram of a vortex-induced vibration response prediction device for a cantilever column structure according to an embodiment of this application;

[0066] Figure 7 This is a schematic diagram of the system structure of an electronic device according to an embodiment of this application. Detailed Implementation

[0067] Exemplary embodiments will now be described more fully with reference to the accompanying drawings. However, these exemplary embodiments can be implemented in many forms and should not be construed as limited to the examples set forth herein; rather, they are provided so that this application will be more complete and comprehensive, and will fully convey the concept of the exemplary embodiments to those skilled in the art.

[0068] Furthermore, the described features, structures, or characteristics can be combined in any suitable manner in one or more embodiments. In the following description, numerous specific details are provided to give a thorough understanding of embodiments of this application. However, those skilled in the art will recognize that the technical solutions of this application can be practiced without one or more of the specific details, or other methods, components, apparatuses, steps, etc., can be employed. In other instances, well-known methods, apparatuses, implementations, or operations are not shown or described in detail to avoid obscuring various aspects of this application.

[0069] The block diagrams shown in the accompanying drawings are merely functional entities and do not necessarily correspond to physically independent entities. These functional entities can be implemented in software, in one or more hardware modules or integrated circuits, or in different network and / or processor devices and / or microcontroller node devices.

[0070] The flowcharts shown in the accompanying drawings are merely illustrative and do not necessarily include all content and operations / steps, nor do they necessarily have to be performed in the described order. For example, some operations / steps can be broken down, while others can be combined or partially combined; therefore, the actual execution order may change depending on the specific circumstances.

[0071] It should be noted that "multiple" in this article refers to two or more. "And / or" describes the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A alone, A and B simultaneously, or B alone. The character " / " generally indicates that the preceding and following related objects have an "or" relationship.

[0072] The technical background of the embodiments of this application is described in detail below:

[0073] Many engineering structures are cantilever column structures with similar and gradually varying cross-sections. Cantilever column structures can be like... Figure 1 As shown, Figure 1 Cantilever column structures with circular cross-sections, such as wind turbine towers, chimneys, bridge towers, and power poles, are examples. These structures have cross-sections of circular, elliptical, or polygonal shapes, which gradually change in size as the structure's height increases—the dimensions may gradually decrease, gradually increase, increase and then decrease, or decrease first and then increase. The towering shape of these cantilever column structures presents new challenges to safety. The taller and more flexible the structure, the higher the probability of vortex-induced vibration. While vortex-induced vibration does not directly cause structural dynamic instability and failure, larger amplitude vortex-induced vibrations can reduce assembly precision during construction, increase construction difficulty, endanger the safety of construction personnel, affect the structural safety during operation, reduce fatigue life, and ultimately affect the normal use of the structure.

[0074] Because the incoming wind has a significant wind profile, the cross-section of the cantilever column structure has a gradual change in the vertical height direction, and there is a weak shear flow effect in the vertical direction. These factors have a significant impact on the analysis accuracy of the vortex-induced vibration problem of the structure.

[0075] Therefore, in order to solve the above problems, this application proposes a method for predicting the vortex-induced vibration response of cantilever column structures that is applicable to weak shear flow effects, so as to reduce calculation errors and improve the analysis accuracy of vortex-induced vibration problems.

[0076] The implementation details of the technical solutions in the embodiments of this application are described in detail below:

[0077] According to one aspect of the embodiments of this application, a method for predicting the vortex-induced vibration response of a cantilever column structure is provided. Figure 2 The flowchart below shows a method for predicting the vortex-induced vibration response of a cantilever column structure according to an embodiment of this application. The method includes at least steps S1 to S8, which are described in detail below:

[0078] In step S1, the target cantilever column structure is divided into structural units to obtain several connected structural units.

[0079] Specifically, by dividing the target cantilever column structure into structural units, several connected structural units can be obtained. Among them, Figure 1 Taking a cantilever column structure with a circular cross-section as an example, the schematic diagram of the resulting structural units can be shown as follows: Figure 3 As shown, this is a cross-sectional structure with a gradually changing cross-section, characterized by a small circle at the top and a large circle at the bottom. i Indicates the top surface. Indicates the diameter of the upper end face. j Indicates the lower end face. Indicates the diameter of the lower end face. e The target surface can be any cross-section in the structural unit. This represents the diameter of the target surface. It should be noted that an arbitrarily selected target surface may not meet the requirements of the preset feature size evaluation system described in the embodiments of this application. Therefore, it is necessary to set a feature size discrimination value to determine whether the requirements are met. If the requirements are not met, the target surface is reselected, that is, the equal cross-section structure is reconstructed.

[0080] In step S2, a model of a constant cross-section structure is constructed for each of the structural units to obtain a constant cross-section structure corresponding to each of the structural units, and a finite element model to be tested corresponding to the target cantilever column structure is established based on each constant cross-section structure.

[0081] Specifically, a model of a constant cross-section structure is constructed for each structural unit. First, it should be noted that since the cross-sectional dimensions of the target surface of each structural unit may differ, the dimensions of the constructed multiple constant cross-section structures may also differ. The following section will explain in detail the model construction of a constant cross-section structure for a single structural unit:

[0082] In one embodiment of the application, the model construction of a uniform cross-section structure for a single structural unit includes the following steps:

[0083] Obtain any target point in the height direction of the structural unit, and construct a cross-section structure with the same height as the structural unit using the cross-section where the target point is located;

[0084] In the constant cross-section structure, the shape of each cross section is consistent with the shape of the cross section where the target point is located.

[0085] When dividing the target cantilever column structure, it is divided into equal-height sections along the height direction. Therefore, the height of each structural unit is consistent, and the cross-sectional dimensions exhibit a gradual change. Thus, the equal-section construction of a single structural unit involves using the cross-section of any target point within the structural unit along the height direction as a reference, constructing a cross-section with the same height as the structural unit and dimensions matching those of the selected target point. This results in a uniform-section structure, as shown below. Figure 3 The cylindrical structure shown on the right. Of course, the constant cross-section structure proposed in this application is not limited to a cylindrical structure. When the cross-section of the target cantilever column structure is rectangular, polygonal, or other shapes, the constructed constant cross-section structure is also adjusted accordingly based on the cross-section of the target cantilever column structure. Figure 3 This is merely an exemplary illustration of a structure with equal cross-sections.

[0086] The above method can be used to construct a uniform cross-section structure that corresponds one-to-one with each structural unit. These uniform cross-section structures are then connected in series according to the series sequence of each structural unit to construct a finite element model to be tested. Subsequent steps can be taken to set the characteristic dimension discrimination values ​​of each uniform cross-section structure in the finite element model to be tested and to determine whether they meet the requirements of the preset characteristic dimension evaluation system.

[0087] In step S3, it is determined whether the finite element model to be tested meets the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system, and the preset vortex-induced force evaluation system.

[0088] In one embodiment of this application, whether the finite element model to be tested meets the requirements of the preset feature size evaluation system is determined by the following steps:

[0089] Obtain the characteristic dimension parameters of each of the structural units, including the upper end face cross-sectional dimensions, lower end face cross-sectional dimensions, and target face cross-sectional dimensions of the structural unit;

[0090] The characteristic dimension discrimination value of each structural unit is determined based on the cross-sectional dimensions of each upper end face, each lower end face, and each target surface.

[0091] If all the discriminant values ​​of the aforementioned feature dimensions are lower than the preset feature dimension threshold, then the finite element model to be tested is determined to meet the requirements of the preset feature dimension evaluation system.

[0092] The characteristic size discrimination value of each of the structural units satisfies the following formula:

[0093] (1)

[0094] In the formula, The feature size discriminant value, The dimensions of the upper end face cross-section are as follows: The dimensions of the lower end face cross-section are as follows: The target surface cross-sectional dimensions are... , , .

[0095] It should be noted that any target point can be any point in the height direction of a structural unit, and therefore the target surface can be any cross-section where the target point is located. However, the dimensions of the target surface cross-section must meet the requirements of the preset characteristic dimension evaluation system, which can be further explained by the following method:

[0096] The upper end face cross-sectional dimensions, lower end face cross-sectional dimensions, and target surface cross-sectional dimensions of the structural unit respectively correspond to Figure 3 In i end face, j end face and e The dimensions of the end face, the upper end face cross-section, the lower end face cross-section, and the target surface cross-section are respectively represented by... , and The preset feature size threshold can be arbitrarily selected between 0.01 and 0.05. The cross-sectional dimensions of the structural element gradually change with height. In finite element modeling, a gradually changing uniform cross-section is used to simulate the actual gradually changing structure. The end cross-section of any structural element is... i End face (top end face) j The end face (lower end face) has the following feature dimensions: , In the finite element model, the characteristic dimension is used as the basis for the design. (Target surface) e Construct a uniform cross-section structure (based on the cross-sectional dimensions corresponding to the end face), let... , , The feature size discrimination value can be set by referring to the above formula (1). The preset feature size threshold is used express, A value of 1% to 5% is recommended when modeling. Only when the requirements of the preset characteristic size evaluation system are met can the characteristic size discrimination value be compared with the preset characteristic size threshold to verify whether the constructed uniform cross-section structure can equivalently simulate the structural effect of the gradually changing structure of the structural unit.

[0097] If the requirements of the preset feature size evaluation system are not met, the target point needs to be reselected, that is, the target surface needs to be reselected. e (End face), and reconstruct a uniform cross-section structure based on the cross-section where the newly selected target face is located, and perform the next round of feature size discrimination value. Reset until the reset is achieved. The value is less than the preset feature size threshold.

[0098] It should be noted that, , and This refers to the collective characteristic dimensional parameters of the target cantilever column structure, and can be applied to cantilever column structures of any cross-sectional size. The following examples illustrate the different cross-sectional size types of the target cantilever column structures:

[0099] In one embodiment of this application, if the cross-sectional dimension of the target cantilever column structure is circular, and the actual structure corresponding to any structural unit is a frustum, the diameter of the circle is used as the characteristic dimension, and the characteristic dimension parameter is... , , Substitute into formula (1) to set the feature size discrimination value. .

[0100] In one embodiment of this application, if the cross-sectional dimension type of the target cantilever column structure is rectangular, and the actual structure corresponding to any structural unit is a frustum, the long side of the rectangle is used as the characteristic dimension parameter, and the characteristic dimension parameter is... L i , L j , Substitute into formula (1) to set the feature size discrimination value .

[0101] In one embodiment of this application, if the cross-sectional dimension type of the target cantilever column structure is square, and the actual structure corresponding to any structural unit is a regular frustum, the side length of the square is used as a characteristic dimension parameter. SL i , SL j , Substitute into formula (1) to set the feature size discrimination value .

[0102] In one embodiment of this application, if the cross-sectional dimension type of the target cantilever column structure is polygonal, and the actual structure corresponding to any structural unit is a frustum, the equivalent diameter of the polygon is used as a characteristic dimension parameter. ED i , ED j , Substitute into formula (1) to set the feature size discrimination value .

[0103] In one embodiment of this application, whether the finite element model to be tested meets the requirements of the average wind load evaluation system is determined by the following steps:

[0104] The average wind load parameters on the upper end face of each structural unit, the average wind load parameters on the lower end face of each structural unit, and the average wind load parameters on the target face of each structural unit are obtained.

[0105] The average wind load parameters of each of the upper end faces, the average wind load parameters of each of the lower end faces, and the average wind load parameters of each of the target faces are input into the preset average wind load calculation formula to obtain the first average wind load of each of the upper end faces, the second average wind load of each of the lower end faces, and the third average wind load of each of the target faces.

[0106] The average wind load discrimination value of each structural unit is determined based on each of the first average wind load, each of the second average wind load, and each of the third average wind load.

[0107] If all the average wind load discrimination values ​​are lower than the preset average wind load threshold, then the finite element model to be tested is determined to meet the requirements of the average wind load evaluation system.

[0108] The average wind load parameters include air density, drag coefficient, and wind speed at any height. The cross-sectional dimensions and arbitrary height of the structural unit at the location The average incoming air velocity;

[0109] any height The average wind load can be obtained using the following formula for calculating the average wind load:

[0110] (2)

[0111] In the formula, The air density is mentioned. The drag coefficient is mentioned above; To be at any height The cross-sectional dimensions of the structural unit at that location. For the arbitrary height The average incoming air velocity; For any height The average wind load on the location.

[0112] Specifically, when the incoming airflow acts on the target cantilever structure, the target cantilever structure at any height in the vertical direction Average wind load at the location The expression can be found in formula (2) above. According to formula (2), in any structural unit, i The height, average wind speed, and average wind load corresponding to the end face are: , , , j The corresponding height, average wind speed, and average wind load of the end face are: , , Target surface e The corresponding height, average wind speed, and average wind load of the end face are: , , This is used to set the average wind load discrimination value. :

[0113] (3)

[0114] The preset average wind load threshold is used express, A value of 1% to 5% is recommended when modeling the target finite element model. Only when the preset average wind load evaluation system is met can the requirements be verified. By comparing the average wind load discrimination value with the preset average wind load threshold, it can be verified whether the uniform cross-section structure can equivalently simulate the wind load effect of the cantilever column structure in the actual scene.

[0115] The distribution of average incoming wind speed along the height direction follows a power law or a logarithmic law, and can be expressed by an exponential or logarithmic model.

[0116] First, when using an exponential model to calculate the target mean wind load, at any height Average incoming wind speed The expression is

[0117] (4)

[0118] In the formula, For reference height; for Reference wind speed at the location; This is the surface roughness coefficient.

[0119] When using a logarithmic model to calculate the target mean wind load, at any height Average incoming wind speed The expression is:

[0120] (5)

[0121] In the formula, This is the von Karman constant, typically taken as 0.4; Zero plane displacement is a function of the properties, height, and distribution of the surface roughness element. It can be understood as the difference between the zero wind speed height and the roughness height. For urban topography, , The length of the ground roughness. The average height of the building; This refers to the frictional wind speed.

[0122] make , , When using an exponential model to calculate the average wind load, , , At this time, the average wind load discriminant value The expression is:

[0123] (6)

[0124] Set average wind load threshold , A value of 1% to 5% is recommended when modeling. Only when the requirements are met can we verify whether a uniform cross-section structure can effectively simulate the wind load effect of a cantilever column structure in a real-world scenario.

[0125] When using a logarithmic model to calculate the average wind load , , At this time, the average wind load discriminant value The expression is:

[0126] (7)

[0127] Set average wind load threshold , A value of 1% to 5% is recommended when modeling. Only when the requirements are met can we verify whether a uniform cross-section structure can effectively simulate the wind load effect of a cantilever column structure in a real-world scenario.

[0128] In one embodiment of this application, whether the finite element model to be tested meets the requirements of the vortex-induced force evaluation system is determined by the following steps:

[0129] Obtain the first vortex-induced force parameter of the upper end face of each structural unit, the second vortex-induced force parameter of the lower end face of each structural unit, and the third vortex-induced force parameter of the target surface of each structural unit.

[0130] The vortex-induced force discrimination value of each structural unit is determined based on each of the first vortex-induced force parameters, each of the second vortex-induced force parameters, and each of the third vortex-induced force parameters.

[0131] If all the vortex-induced force discrimination values ​​are lower than the preset vortex-induced force threshold, then the finite element model to be tested is determined to meet the requirements of the vortex-induced force evaluation system.

[0132] Specifically, the first vortex excitation parameter, the second vortex excitation parameter, and the third vortex excitation parameter can all be obtained by the following formula (8):

[0133] (8)

[0134] In the formula, These are the parameters of the vortex-induced force. This refers to the structural vibration displacement; The structural vibration velocity; Represents aerodynamic parameters; M , N Let the order of the polynomial in the Taylor expansion be . The frequency of vortex shedding; This is the initial phase; For the parameters to be identified, To calculate the angular frequency, The natural angular frequency of the structure.

[0135] In any structural element, the end section of the structural element ( i end face, j The vortex-induced force parameters corresponding to the end face are respectively , The vortex-induced force parameters are taken at the same positions as the characteristic dimension parameters of the isosection structure in the actual structural element and the finite element model. ,make , , Set the vortex-induced force discrimination value :

[0136] (9)

[0137] Set vortex-induced force threshold , A value of 1% to 5% is recommended when modeling. Only when the requirements are met can we verify whether a uniform cross-section structure can simulate the vortex-induced force effect of a cantilever column structure in a real-world scenario.

[0138] Commonly used eddy-induced force models include the Scanlan empirical linear model, the Scanlan empirical nonlinear model, and the fifth-order polynomial refined simplification model.

[0139] When the vortex-induced force model adopts the Scanlan empirical linear model, the expression for the target vortex-induced force parameters is:

[0140] (10)

[0141] From equation (10) we get , , Substituting into equation (9), we obtain the vortex-induced force discrimination value. :

[0142] (11)

[0143] In the formula, , ,

[0144] ,

[0145] .

[0146] Set vortex-induced force threshold , A value of 1% to 5% is recommended when modeling. Only when the requirements are met can we verify whether a uniform cross-section structure can simulate the vortex-induced force effect of a cantilever column structure in a real-world scenario.

[0147] When the vortex-induced force model adopts the Scanlan empirical nonlinear model, the expression for the target vortex-induced force parameters is:

[0148] (12)

[0149] From equation (12) we get , , Substituting into equation (9), we obtain the vortex-induced force discrimination value. :

[0150] (13)

[0151] Set vortex-induced force threshold , A value of 1% to 5% is recommended when modeling. Only when the requirements are met can we verify whether a uniform cross-section structure can simulate the vortex-induced force effect of a cantilever column structure in a real-world scenario.

[0152] When the vortex-induced force model adopts a fifth-order polynomial refined simplification model, the expression for the target vortex-induced force parameters is:

[0153] (14)

[0154] Wherein, from equation (14) we get , , Substituting into equation (9), we obtain the vortex-induced force discrimination value. :

[0155] (15)

[0156] Set vortex-induced force threshold , A value of 1% to 5% is recommended when modeling. Only when the requirements are met can we verify whether a uniform cross-section structure can simulate the vortex-induced force effect of a cantilever column structure in a real-world scenario.

[0157] In step S4, if the finite element model to be tested does not meet the requirements of the characteristic size evaluation system, the average wind load evaluation system, or the vortex-induced force evaluation system, then steps S1-S2 are repeated until the finite element model to be tested meets the requirements of the characteristic size evaluation system, the average wind load evaluation system, and the vortex-induced force evaluation system, and the finite element model to be tested is determined to be the target finite element model.

[0158] Specifically, please refer to Figure 4 As shown, Figure 4 The logical flowchart for the entire vortex-induced vibration response prediction requires that the finite element model to be tested fully meet the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system, and the preset vortex-induced force evaluation system. In other words, the finite element model can completely and equivalently simulate the target cantilever column structure only if the three conditions are met: the discrimination value of each of the aforementioned characteristic sizes is lower than the preset characteristic size threshold, the discrimination value of each of the aforementioned average wind loads is lower than the preset average wind load threshold, and the discrimination value of each of the aforementioned vortex-induced force is lower than the preset vortex-induced force threshold.

[0159] In step S5, the structural characteristic parameters of the target cantilever column structure are obtained, and the target average wind load and the target vortex-induced force parameters under the average wind load evaluation system are determined based on the structural characteristic parameters.

[0160] Specifically, since the average wind load parameters and vortex-induced force parameters of each structural unit can be obtained in the above steps, the target average wind load on the target finite element model under the average wind load evaluation system can be determined. And the target vortex-induced force parameters under the vortex-induced force evaluation system. .

[0161] In step S6, the coordinate system transformation of each target structural element in the target finite element model is performed according to the target average wind load and the target vortex-induced force parameters to obtain the overall average wind load and overall vortex-induced force parameters of the target cantilever column structure under the preset coordinate system.

[0162] Specifically, the element forces in the local coordinate system are converted into forces in the global coordinate system through the element transformation matrix, which is the overall average wind load on the target cantilever column structure. And the target cantilever column structure subjected to Based on this, the dynamic equilibrium equations are established, and the structural vortex-induced vibration response is solved using numerical analysis methods.

[0163] The dynamic equilibrium equations are:

[0164] (16)

[0165] In the formula, m , c , k These represent the average mass per meter of the cantilever column structure, its damping, and its stiffness, respectively.

[0166] In step S7, the structural vibration acceleration, structural vibration velocity, and structural vibration displacement of the target cantilever column structure are determined according to the structural characteristic parameters. The dynamic equilibrium equation corresponding to the target cantilever column structure is constructed based on the structural vibration acceleration, the structural vibration velocity, the structural vibration displacement, the overall average wind load, the overall vortex-induced force parameters, and the preset weak shear flow effect correction coefficient.

[0167] In one embodiment of this application, the structural characteristic parameters include the average mass per meter, average damping per meter, and average stiffness per meter of the cantilever column structure; the dynamic equilibrium equation is specifically the above formula (16), where, The average mass per meter of the target cantilever column structure. The average damping per meter of the target cantilever column structure is given. The average stiffness per meter of the target cantilever column structure; The vibration acceleration of the structure, The vibration velocity of the structure is... The vibration displacement of the structure; The overall average wind load is... The overall vortex-induced force parameters are as follows: This is the correction coefficient for the eddy-induced weak shear flow effect.

[0168] Specifically, the weak shear flow effect in the vertical direction of the cantilever column significantly reduces the accuracy of the constructed finite element model when simulating the cantilever column. Therefore, a correction factor for the weak shear flow effect of vortex-induced force can be used. To correct, It is a function of time, frequency, and axial coordinates. It can be obtained through experiments or numerical simulations. Experiments can be conducted by using a series narrow-band aeroelastic model to simultaneously test the vortex-induced force of each narrow band and obtaining the result through correlation analysis; alternatively, it can be obtained by testing the vortex-induced force parameters of different cross-sections and then analyzing the results through experiments.

[0169] In step S8, the vortex-induced vibration response of the target cantilever column structure under the action of incoming wind is predicted according to the dynamic equilibrium equation.

[0170] Specifically, the constructed dynamic equilibrium equations can accurately predict the vortex-induced vibration response of the cantilever column structure based on the specific parameters of the incoming wind.

[0171] In one embodiment of this application, such as Figure 5 As shown, the finite element model is constructed using finite element software to create a lumped mass model. Taking a tower structure as an example, the model is divided into several elements of a fixed length from top to bottom elevation. The "skewer of candied hawthorns" model concept is used to construct a finite element model of several elements connected in series. In terms of element division, a gradually changing uniform cross-section is used to simulate the actual slowly varying structure. The model elements are beam elements. By defining different cross-sectional areas (section dimensions) and moments of inertia (integrals related to area and dimension) for each element, the change in the cross-sectional dimensions of the tower with height is reflected, simplifying the tower into an additional lumped mass added to the top node. superior.

[0172] In summary, the embodiments of this application divide the target cantilever column structure into structural units to obtain several connected structural units. The model of each divided structural unit with an equal cross-section can be constructed, thereby establishing a finite element model to be tested corresponding to the target cantilever column structure.

[0173] By determining whether the finite element model to be tested meets the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system, and the preset vortex-induced force evaluation system, the final target finite element model meets the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system, and the preset vortex-induced force evaluation system, thereby equivalently simulating the vortex-induced vibration that may occur in the target cantilever column structure in the actual scenario.

[0174] Furthermore, by obtaining the structural characteristic parameters of the target cantilever column structure and determining that the finite element model to be tested is the target average wind load under the average wind load evaluation system and the target vortex-induced force parameters under the vortex-induced force evaluation system, the coordinate system transformation of the element transformation matrix of each target structural unit is performed to equivalently simulate the overall average wind load and overall vortex-induced force parameters of the target cantilever column structure under the preset coordinate system.

[0175] Furthermore, by determining the structural vibration acceleration, vibration velocity, and vibration displacement of the target cantilever column structure through structural characteristic parameters, and correcting them according to a preset weak shear flow effect correction coefficient, the influence of weak shear flow in the vertical direction of the cantilever column structure's cross-section due to its gradual change characteristics can be eliminated, thereby improving the accuracy and precision of the analysis of vortex-induced vibration problems. Simultaneously, the constructed dynamic equilibrium equations can be used to predict the vortex-induced vibration response of the target cantilever column structure under the influence of incoming wind. By monitoring the parameters of the incoming wind, it is possible to predict in real time whether the cantilever column structure will experience vortex-induced vibration response.

[0176] Therefore, the vortex-induced vibration response prediction method for cantilever column structures provided in this application can eliminate the influence of weak shear flow and establish dynamic equilibrium equations to accurately predict the vortex-induced vibration response.

[0177] Figure 6 The diagram shows a block diagram of a vortex-induced vibration response prediction device 300 for a cantilever column structure according to an embodiment of this application. The device 300 according to an embodiment of this application includes: a division unit 301, a construction unit 302, a judgment unit 303, an iteration unit 304, an acquisition unit 305, a conversion unit 306, a dynamic balance unit 307, and a prediction unit 308.

[0178] Dividing unit 301 is used to divide the target cantilever column structure into structural units, resulting in several connected structural units;

[0179] Construction unit 302 is used to construct a model of a constant cross-section structure for each of the structural units, to obtain a constant cross-section structure corresponding to each of the structural units, and to establish a finite element model to be tested corresponding to the target cantilever column structure based on each of the constant cross-section structures.

[0180] The judgment unit 303 is used to judge whether the finite element model to be tested meets the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system and the preset vortex-induced force evaluation system.

[0181] The iteration unit 304 is used to repeat the above steps if the finite element model to be tested does not meet the requirements of the characteristic size evaluation system, the average wind load evaluation system, or the vortex-induced force evaluation system, until the finite element model to be tested meets the requirements of the characteristic size evaluation system, the average wind load evaluation system, and the vortex-induced force evaluation system, and then determines the finite element model to be tested as the target finite element model.

[0182] The acquisition unit 305 is used to acquire the structural characteristic parameters of the target cantilever column structure, and determine the target average wind load and the target vortex-induced force parameters of the target finite element model under the average wind load evaluation system based on the structural characteristic parameters.

[0183] The transformation unit 306 is used to perform coordinate system transformation of the element transformation matrix of each target structural element in the target finite element model according to the target average wind load and the target vortex-induced force parameters, so as to obtain the overall average wind load and overall vortex-induced force parameters of the target cantilever column structure under the preset coordinate system.

[0184] The dynamic balance unit 307 is used to determine the structural vibration acceleration, structural vibration velocity, and structural vibration displacement of the target cantilever column structure according to the structural characteristic parameters, and to construct the dynamic balance equation corresponding to the target cantilever column structure according to the structural vibration acceleration, structural vibration velocity, structural vibration displacement, overall average wind load, overall vortex-induced force parameters, and preset weak shear flow effect correction coefficient.

[0185] The prediction unit 308 is used to predict the vortex-induced vibration response of the target cantilever column structure when subjected to the incoming wind, based on the dynamic equilibrium equation.

[0186] In another aspect, this application also provides a computer-readable storage medium storing a program product capable of implementing the methods provided above in this specification. In some possible implementations, various aspects of this application may also be implemented as a program product comprising program code that, when run on a terminal device, causes the terminal device to perform the steps described in the "Embodiment Methods" section of this specification according to various exemplary embodiments of this application.

[0187] The program product for implementing the above-described method according to the embodiments of this application may employ a portable compact disc read-only memory (CD-ROM) and include program code, and may run on a terminal device, such as a personal computer. However, the program product of this application is not limited thereto. In this document, the readable storage medium may be any tangible medium containing or storing a program that may be used by or in conjunction with an instruction execution system, apparatus, or device.

[0188] The program product may employ any combination of one or more readable media. A readable medium may be a readable signal medium or a readable storage medium. A readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of readable storage media (a non-exhaustive list) include: an electrical connection having one or more wires, a portable disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.

[0189] Computer-readable signal media may include data signals propagated as part of a carrier wave in baseband, carrying readable program code. Such propagated data signals may take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A readable signal medium may also be any readable medium other than a readable storage medium, capable of sending, propagating, or transmitting programs for use by or in connection with an instruction execution system, apparatus, or device.

[0190] The program code contained on the readable medium may be transmitted using any suitable medium, including but not limited to wireless, wired, optical fiber, RF, etc., or any suitable combination thereof.

[0191] Program code for performing the operations of this application can be written in any combination of one or more programming languages, including object-oriented programming languages ​​such as Java and C++, and conventional procedural programming languages ​​such as C or similar languages. The program code can execute entirely on the user's computing device, partially on the user's device, as a standalone software package, partially on the user's computing device and partially on a remote computing device, or entirely on a remote computing device or server. In cases involving remote computing devices, the remote computing device can be connected to the user's computing device via any type of network, including a local area network (LAN) or a wide area network (WAN), or it can be connected to an external computing device (e.g., via the Internet using an Internet service provider).

[0192] In another respect, this application also provides an electronic device capable of implementing the above-described method.

[0193] Those skilled in the art will understand that various aspects of this application can be implemented as a system, method, or program product. Therefore, various aspects of this application can be specifically implemented in the following forms: a completely hardware implementation, a completely software implementation (including firmware, microcode, etc.), or a combination of hardware and software implementations, collectively referred to herein as a "circuit," "module," or "system."

[0194] The following reference Figure 7 To describe an electronic device 400 according to this embodiment of the present application. Figure 7 The electronic device 400 shown is merely an example and should not impose any limitations on the functionality and scope of use of the embodiments of this application.

[0195] like Figure 7 As shown, the electronic device 400 is manifested in the form of a general-purpose computing device. The components of the electronic device 400 may include, but are not limited to: at least one processing unit 410, at least one storage unit 420, and a bus 430 connecting different system components (including storage unit 420 and processing unit 410).

[0196] The storage unit stores program code that can be executed by the processing unit 410, causing the processing unit 410 to perform the steps described in the "Embodiment Methods" section above according to various exemplary embodiments of this application.

[0197] Storage unit 420 may include a readable medium in the form of a volatile storage unit, such as random access memory (RAM) 421 and / or cache memory 422, and may further include a read-only memory (ROM) 423.

[0198] Storage unit 420 may also include a program / utility 424 having a set (at least one) of program modules 425, such program modules 425 including but not limited to: an operating system, one or more application programs, other program modules and program data, each of these examples or some combination thereof may include an implementation of a network environment.

[0199] Bus 430 can represent one or more of several types of bus structures, including a memory cell bus or memory cell control node, a peripheral bus, a graphics acceleration port, a processing unit, or a local bus using any of the various bus structures.

[0200] Electronic device 400 can also communicate with one or more external devices 1200 (e.g., keyboard, pointing device, Bluetooth device, etc.), and with one or more devices that enable a user to interact with electronic device 400, and / or with any device that enables electronic device 400 to communicate with one or more other computing devices (e.g., router, modem, etc.). This communication can be performed via input / output (I / O) interface 450. Furthermore, electronic device 400 can also communicate with one or more networks (e.g., local area network (LAN), wide area network (WAN), and / or public networks, such as the Internet) via network adapter 460. As shown, network adapter 460 communicates with other modules of electronic device 400 via bus 430. It should be understood that, although not shown in the figures, other hardware and / or software modules can be used in conjunction with electronic device 400, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage systems.

[0201] Through the above description of the embodiments, those skilled in the art will readily understand that the exemplary embodiments described herein can be implemented by software or by combining software with necessary hardware. Therefore, the technical solutions according to the embodiments of this application can be embodied in the form of a software product, which can be stored in a non-volatile storage medium (such as a CD-ROM, USB flash drive, external hard drive, etc.) or on a network, including instructions to cause a computing device (such as a personal computer, server, terminal device, or network device, etc.) to execute the methods according to the embodiments of this application.

[0202] Furthermore, the above figures are merely illustrative of the processes included in the method according to exemplary embodiments of this application, and are not intended to be limiting. It is readily understood that the processes shown in the above figures do not indicate or limit the temporal order of these processes. Additionally, it is readily understood that these processes may be executed synchronously or asynchronously in multiple modules.

[0203] It should be understood that this application is not limited to the precise structures described above and shown in the accompanying drawings, and various modifications and changes can be made without departing from its scope. The scope of this application is limited only by the appended claims.

Claims

1. A method for predicting the vortex-induced vibration response of a cantilever column structure, characterized in that, The method includes: S1, Divide the target cantilever column structure into structural units to obtain several connected structural units; S2, construct a model of a constant cross-section structure for each of the structural units to obtain a constant cross-section structure corresponding to each of the structural units, and establish a finite element model to be tested corresponding to the target cantilever column structure based on each of the constant cross-section structures; S3, determine whether the finite element models to be tested all meet the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system and the preset vortex-induced force evaluation system. S4. If the finite element model to be tested does not meet the requirements of the characteristic size evaluation system, the average wind load evaluation system, or the vortex-induced force evaluation system, then repeat steps S1-S3 until the finite element model to be tested meets the requirements of the characteristic size evaluation system, the average wind load evaluation system, and the vortex-induced force evaluation system, and determine the finite element model to be tested as the target finite element model. S5, obtain the structural characteristic parameters of the target cantilever column structure, and determine the target average wind load and the target vortex-induced force parameters of the target finite element model under the average wind load evaluation system based on the structural characteristic parameters. S6. Based on the target average wind load and the target vortex-induced force parameters, the coordinate system transformation of the element transformation matrix of each target structural unit in the target finite element model is performed to obtain the overall average wind load and overall vortex-induced force parameters of the target cantilever column structure under the preset coordinate system. S7. Determine the structural vibration acceleration, structural vibration velocity, and structural vibration displacement of the target cantilever column structure based on the structural characteristic parameters, and construct the dynamic equilibrium equation corresponding to the target cantilever column structure based on the structural vibration acceleration, structural vibration velocity, structural vibration displacement, overall average wind load, overall vortex-induced force parameters, and preset vortex-induced force weak shear flow effect correction coefficient. S8. Based on the dynamic equilibrium equation, predict the vortex-induced vibration response of the target cantilever column structure when subjected to the incoming wind. The process of constructing a model of a uniform cross-section structure for a single structural unit includes the following steps: Obtain any target point in the height direction of the structural unit, and construct a cross-section structure with the same height as the structural unit using the cross-section where the target point is located; In the constant cross-section structure, the shape of each cross section is consistent with the shape of the cross section where the target point is located; Whether the finite element model to be tested meets the requirements of the preset feature size evaluation system is determined by the following steps: Obtain the characteristic dimension parameters of each of the structural units, the characteristic dimension parameters including the upper end face cross-sectional dimension, the lower end face cross-sectional dimension and the target face cross-sectional dimension of the structural unit, the target face cross-sectional dimension corresponding to the dimension of the cross-section where the target point is located; The characteristic dimension discrimination value of each structural unit is determined based on the cross-sectional dimensions of each upper end face, each lower end face, and each target surface. If all the discriminant values ​​of the aforementioned feature dimensions are lower than the preset feature dimension threshold, then the finite element model to be tested is determined to meet the requirements of the preset feature dimension evaluation system. The characteristic size discrimination value of each of the structural units satisfies the following formula: In the formula, The feature size discriminant value, The dimensions of the upper end face cross-section are as follows: The dimensions of the lower end face cross-section are as follows: The target surface cross-sectional dimensions are... , , .

2. The method for predicting the vortex-induced vibration response of a cantilever column structure according to claim 1, characterized in that, Whether the finite element model to be tested meets the requirements of the average wind load evaluation system is determined by the following steps: The average wind load parameters on the upper end face of each structural unit, the average wind load parameters on the lower end face of each structural unit, and the average wind load parameters on the target face of each structural unit are obtained. The average wind load parameters of each of the upper end faces, the average wind load parameters of each of the lower end faces, and the average wind load parameters of each of the target faces are input into the preset average wind load calculation formula to obtain the first average wind load of each of the upper end faces, the second average wind load of each of the lower end faces, and the third average wind load of each of the target faces. The average wind load discrimination value of each structural unit is determined based on each of the first average wind load, each of the second average wind load, and each of the third average wind load. If all the average wind load discrimination values ​​are lower than the preset average wind load threshold, then the finite element model to be tested is determined to meet the requirements of the average wind load evaluation system. The average wind load parameters include air density, drag coefficient, and wind speed at any height. The cross-sectional dimensions and arbitrary height of the structural unit at the location The average incoming air velocity; any height The average wind load is obtained using the following formula for calculating the average wind load: In the formula, The air density is mentioned. The drag coefficient is mentioned above; To be at any height The cross-sectional dimensions of the structural unit at that location. For the arbitrary height The average incoming air velocity; For any height The average wind load experienced at the location.

3. The method for predicting the vortex-induced vibration response of a cantilever column structure according to claim 2, characterized in that, The average wind load discrimination value is obtained by the following formula: In the formula, The average wind load discrimination value is... The first average wind load, This is the second average wind load. The third average wind load, The average incoming wind speed at the height of the upper end face. The average incoming wind speed at the height of the lower end face. The average incoming wind speed at the height of the target surface.

4. The method for predicting the vortex-induced vibration response of a cantilever column structure according to claim 3, characterized in that, Whether the finite element model to be tested meets the requirements of the vortex-induced force evaluation system is determined by the following steps: Obtain the first vortex-induced force parameter of the upper end face of each structural unit, the second vortex-induced force parameter of the lower end face of each structural unit, and the third vortex-induced force parameter of the target surface of each structural unit. The vortex-induced force discrimination value of each structural unit is determined based on each of the first vortex-induced force parameters, each of the second vortex-induced force parameters, and each of the third vortex-induced force parameters. If all the vortex-induced force discrimination values ​​are lower than the preset vortex-induced force threshold, then the finite element model to be tested is determined to meet the requirements of the vortex-induced force evaluation system.

5. The method for predicting the vortex-induced vibration response of a cantilever column structure according to claim 4, characterized in that, The structural characteristic parameters include the average mass per meter, average damping per meter, and average stiffness per meter of the target cantilever column structure; the dynamic equilibrium equation is specifically the following formula: In the formula, The average mass per meter of the target cantilever column structure. The average damping per meter of the target cantilever column structure is given. The average stiffness per meter of the target cantilever column structure; The vibration acceleration of the structure, The vibration velocity of the structure is... The vibration displacement of the structure; The overall average wind load is... The overall vortex-induced force parameters are as follows: This is the preset correction coefficient for the weak shear flow effect of vortex-induced force.

6. A device for predicting the vortex-induced vibration response of a cantilever column structure, characterized in that, The device includes: The unit division is used to divide the target cantilever column structure into structural units, resulting in several connected structural units; The construction unit is used to construct a model of a constant cross-section structure for each of the structural units, thereby obtaining a constant cross-section structure that corresponds one-to-one with each of the structural units, and to establish a finite element model to be tested corresponding to the target cantilever column structure based on each of the constant cross-section structures. The judgment unit is used to determine whether the finite element model to be tested meets the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system, and the preset vortex-induced force evaluation system. An iterative unit is used to repeatedly execute the execution steps of the partitioning unit, the construction unit, and the judgment unit if the finite element model to be tested does not meet the requirements of the characteristic size evaluation system, the average wind load evaluation system, or the vortex-induced force evaluation system, until the finite element model to be tested meets the requirements of the characteristic size evaluation system, the average wind load evaluation system, and the vortex-induced force evaluation system, and then determines the finite element model to be tested as the target finite element model. The acquisition unit is used to acquire the structural characteristic parameters of the target cantilever column structure, and determine the target average wind load and the target vortex-induced force parameters of the target finite element model under the average wind load evaluation system based on the structural characteristic parameters. The transformation unit is used to perform coordinate system transformation of the element transformation matrix of each target structural element in the target finite element model according to the target average wind load and the target vortex-induced force parameters, so as to obtain the overall average wind load and overall vortex-induced force parameters of the target cantilever column structure under the preset coordinate system. The dynamic equilibrium unit is used to determine the structural vibration acceleration, structural vibration velocity, and structural vibration displacement of the target cantilever column structure based on the structural characteristic parameters, and to construct the dynamic equilibrium equation corresponding to the target cantilever column structure based on the structural vibration acceleration, structural vibration velocity, structural vibration displacement, overall average wind load, overall vortex-induced force parameters, and a preset vortex-induced force weak shear flow effect correction coefficient. The prediction unit is used to predict the vortex-induced vibration response of the target cantilever column structure when subjected to the incoming wind, based on the dynamic equilibrium equation. The process of constructing a model of a uniform cross-section structure for a single structural unit includes the following steps: Obtain any target point in the height direction of the structural unit, and construct a cross-section structure with the same height as the structural unit using the cross-section where the target point is located; In the constant cross-section structure, the shape of each cross section is consistent with the shape of the cross section where the target point is located; Whether the finite element model to be tested meets the requirements of the preset feature size evaluation system is determined by the following steps: Obtain the characteristic dimension parameters of each of the structural units, the characteristic dimension parameters including the upper end face cross-sectional dimension, the lower end face cross-sectional dimension and the target face cross-sectional dimension of the structural unit, the target face cross-sectional dimension corresponding to the dimension of the cross-section where the target point is located; The characteristic dimension discrimination value of each structural unit is determined based on the cross-sectional dimensions of each upper end face, each lower end face, and each target surface. If all the discriminant values ​​of the aforementioned feature dimensions are lower than the preset feature dimension threshold, then the finite element model to be tested is determined to meet the requirements of the preset feature dimension evaluation system. The characteristic size discrimination value of each of the structural units satisfies the following formula: In the formula, The feature size discriminant value, The dimensions of the upper end face cross-section are as follows: The dimensions of the lower end face cross-section are as follows: The target surface cross-sectional dimensions are... , , .

7. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores at least one piece of program code, which is loaded and executed by a processor to perform the operations performed by the method as described in any one of claims 1 to 5.

8. An electronic device, characterized in that, The electronic device includes one or more processors and one or more memories, wherein at least one piece of program code is stored in the one or more memories, and the at least one piece of program code is loaded and executed by the one or more processors to perform the operation performed by the method as described in any one of claims 1 to 5.

Images