Vortex vibration response prediction method and device of cantilever column structure, equipment and medium
By dividing the structural unit of the cantilever column structure and constructing iso-sectional model, combining characteristic size, average wind load and vortex excitation evaluation system, the vortex vibration prediction problem of cantilever column structure is solved, and the analysis accuracy and safety are improved.
Patent Information
- Application Number
- CN202411769622.0
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-12-04
- Publication Date
- 2025-05-06
- Estimated Expiration
- 2044-12-04
AI Technical Summary
Cantilever cylindrical structures are prone to vortex vibration at high places, resulting in increased construction difficulty, increased safety risks and reduced structural fatigue life. It is difficult for the existing technology to effectively predict and analyze vortex vibration problems.
By dividing the structural element of the cantilever column structure, constructing other cross-sectional structural models, establishing a finite element model, and judging and correcting it through the characteristic size, average wind load and vortex excitation evaluation system until the preset requirements are met, a dynamic equilibrium equation is constructed to predict the eddy vibration response.
The analysis accuracy and accuracy of the vortex vibration problem of cantilever column structure is improved, and it can predict in real time whether the structure will undergo vortex vibration response, reducing safety risks and structural damage in construction and operation.
Smart Images

Figure CN119939979A_ABST
Abstract
Description
Background Art
[0002] A large number of engineering structures are cantilever column structures with similar and slowly changing cross-sections, such as wind turbine towers, chimneys, bridge towers, power steel poles, etc. The cross-section types of these structures are circular, elliptical, polygonal, etc. As the height of the structure increases, the cross-section size changes slowly, the size gradually decreases or increases, or increases first and then decreases, decreases first and then increases, etc. This type of cantilever column structure has a towering appearance, which brings new challenges to safety. The taller and softer the structure, the higher the probability of vortex-induced vibration. Although vortex vibration (vortex-induced vibration) will not directly cause structural dynamic instability and damage, vortex vibration with a large amplitude will reduce the assembly accuracy during construction, increase the difficulty of construction, endanger the safety of construction personnel, affect the safety of the structure during operation, reduce fatigue life, and thus affect the normal use of the structure. Therefore, it is necessary to propose a prediction method for structural vortex-induced vibration to improve the analysis accuracy and precision of vortex-induced vibration problems. Summary of the invention
[0003] The present application provides a method, device, equipment and medium for predicting the vortex-vibration response of a cantilever column structure to solve one or more technical problems existing in the prior art and at least provide a beneficial choice or create conditions.
[0004] Other features and advantages of the present application will become apparent from the following detailed description, or may be learned in part by the practice of the present application.
[0005] According to one aspect of an embodiment of the present application, a method for predicting vortex-induced vibration response of a cantilever column structure is provided, the method comprising: S1, dividing the target cantilever column structure into structural units to obtain a number of structural units connected in series; S2, constructing a model of a uniform cross-section structure for each of the structural units, obtaining a uniform cross-section structure corresponding to each of the structural units, and establishing a finite element model to be tested corresponding to the target cantilever column structure according to each of the uniform cross-section structures; S3, judging whether the finite element models to be tested all meet the requirements of a preset characteristic size evaluation system, a preset average wind load evaluation system and a preset vortex excitation force evaluation system; S4, if the finite element model to be tested does not meet the requirements of the characteristic size evaluation system or the average wind load evaluation system or the vortex excitation force evaluation system, 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 excitation force evaluation system, and the finite element model to be tested is determined to be the target finite element model; S5, acquiring structural characteristic parameters of the target cantilever column structure, and determining the target average wind load received by the target finite element model under the average wind load evaluation system and the target vortex excitation force parameters under the vortex excitation force evaluation system according to the structural characteristic parameters; S6, performing a coordinate system transformation of a unit transformation matrix on each target structural unit in the target finite element model according to the target average wind load and the target vortex excitation force parameter, to obtain an overall average wind load and an overall vortex excitation force parameter on the target cantilever column structure in a preset coordinate system; S7, determining the structural vibration acceleration, structural vibration velocity and structural vibration displacement of the target cantilever column structure according to the structural characteristic parameters, and constructing a dynamic equilibrium equation corresponding to the target cantilever column structure according to the structural vibration acceleration, the structural vibration velocity, the structural vibration displacement, the overall average wind load, the overall vortex excitation parameter and a preset weak shear flow effect correction coefficient; S8, predicting the vortex-oscillation response of the target cantilever column structure when it is affected by incoming wind according to the dynamic equilibrium equation.
[0006] In one embodiment of the present application, based on the above scheme, a model of a uniform cross-section structure is constructed for a single structural unit, comprising the following steps: Acquire any target point of the structural unit in the height direction, and construct a uniform cross-sectional structure having the same height as the structural unit based on the cross section where the target point is located; The shape of each cross section in the equal cross-section structure is consistent with the shape of the cross section where the target point is located.
[0007] In one embodiment of the present application, based on the above scheme, whether the finite element model to be tested meets the requirements of the preset characteristic dimension evaluation system is determined by the following steps: Acquire characteristic size parameters of each of the structural units, wherein the characteristic size parameters include the cross-sectional size of the upper end surface, the cross-sectional size of the lower end surface, and the cross-sectional size of the target surface of the structural unit, wherein the cross-sectional size of the target surface corresponds to the size of the cross section where the target point is located; Determine the characteristic size discrimination value of each of the structural units according to each of the upper end surface cross-sectional dimensions, each of the lower end surface cross-sectional dimensions and each of the target surface cross-sectional dimensions; If each of the characteristic dimension discrimination values is lower than a preset characteristic dimension threshold, it is determined that the finite element model to be inspected meets the requirements of the preset characteristic dimension evaluation system; Wherein, the characteristic size discrimination value of each of the structural units satisfies the following formula:
[0008] In the formula, is the characteristic size discrimination value, is the cross-sectional dimension of the upper end face, is the cross-sectional dimension of the lower end face, is the cross-sectional size of the target surface, , , .
[0009] In one embodiment of the present application, based on the above 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: Obtaining an average wind load parameter received by the upper end surface of each of the structural units, an average wind load parameter received by the lower end surface of each of the structural units, and an average wind load parameter received by the target surface of each of the structural units; Input the average wind load parameters received by each of the upper end faces, the average wind load parameters received by each of the lower end faces, and the average wind load parameters received by each of the target faces into a preset average wind load calculation formula to obtain a first average wind load received by each of the upper end faces, a second average wind load received by each of the lower end faces, and a third average wind load received by each of the target faces; Determining an average wind load discrimination value of each of the structural units according to each of the first average wind loads, each of the second average wind loads, and each of the third average wind loads; If each of the average wind load discrimination values is lower than a preset average wind load threshold, it is determined that the finite element model to be tested meets the requirements of the average wind load evaluation system; The average wind load parameters include air density, drag coefficient, The cross-sectional dimensions of the cross section of the structural unit and any height The average incoming wind speed; Any height The average wind load can be calculated using the following average wind load calculation formula:
[0010] In the formula, is the air density; is the drag coefficient; For any height The cross-sectional dimensions of the cross section of the structural unit at For any height The average incoming wind speed; For any height The average wind load on the space.
[0011] In one embodiment of the present application, based on the above solution, the average wind load discrimination value can be obtained by the following formula:
[0012] In the formula, is the average wind load discrimination value, is the first average wind load, is the second average wind load, is the third average wind load; is the average incoming wind speed at the height of the upper end surface, is the average incoming wind speed at the height of the lower end surface, It is the average incoming wind speed at the height of the target surface.
[0013] In one embodiment of the present application, based on the above 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: Acquire a first vortex-induced force parameter of an upper end surface of each of the structural units, a second vortex-induced force parameter of a lower end surface of each of the structural units, and a third vortex-induced force parameter of a target surface of each of the structural units; Determine the vortex excitation force discrimination value of each of the structural units according to each of the first vortex excitation force parameters, each of the second vortex excitation force parameters, and each of the third vortex excitation force parameters; If all the vortex excitation force discrimination values are lower than the preset vortex excitation force threshold, it is determined that the finite element model to be tested meets the requirements of the vortex excitation force evaluation system.
[0014] In one embodiment of the present application, based on the above scheme, the structural characteristic parameters include the average mass per linear meter, the average damping per linear meter and the average stiffness per linear meter of the target cantilever column structure; the dynamic equilibrium equation is specifically the following formula:
[0015] In the formula, is the average mass per linear meter of the target cantilever column structure, is the average damping per linear meter of the target cantilever column structure, is the average stiffness per linear meter of the target cantilever column structure; is the structural vibration acceleration, is the structural vibration velocity, is the structural vibration displacement; is the overall average wind load, is the overall vortex force parameter, is the correction coefficient of the vortex-induced weak shear flow effect.
[0016] According to one aspect of an embodiment of the present application, a device for predicting vortex-induced vibration response of a cantilever column structure is provided, the device comprising: A division unit is used to divide the target cantilever column structure into structural units to obtain a plurality of structural units connected in series; A construction unit, used to construct a model of a uniform cross-section structure for each of the structural units, obtain a uniform 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 according to each of the uniform cross-section structures; A judgment unit, used to judge whether the finite element models to be tested all meet the requirements of a preset characteristic size evaluation system, a preset average wind load evaluation system and a preset vortex excitation force evaluation system; an iterative unit, for repeating the above steps if the finite element model to be tested does not meet the requirements of the characteristic size evaluation system or the average wind load evaluation system or the vortex excitation 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 excitation force evaluation system, and determining that the finite element model to be tested is a target finite element model; an acquisition unit, used for acquiring structural characteristic parameters of the target cantilever column structure, and determining the target average wind load received by the target finite element model under the average wind load evaluation system and the target vortex excitation force parameters under the vortex excitation force evaluation system according to the structural characteristic parameters; A conversion unit, used for performing a coordinate system conversion of a unit conversion matrix on each target structural unit in the target finite element model according to the target average wind load and the target vortex excitation force parameter, so as to obtain the overall average wind load and the overall vortex excitation force parameter on the target cantilever column structure in a preset coordinate system; a dynamic balance unit, for determining the structural vibration acceleration, structural vibration velocity and structural vibration displacement of the target cantilever column structure according to the structural characteristic parameters, and constructing a dynamic balance equation corresponding to the target cantilever column structure according to the structural vibration acceleration, the structural vibration velocity, the structural vibration displacement, the overall average wind load, the overall vortex excitation force parameter and a preset weak shear flow effect correction coefficient; The prediction unit is used to predict the vortex-oscillation response of the target cantilever column structure when it is affected by the incoming wind according to the dynamic equilibrium equation.
[0017] According to one aspect of an embodiment of the present application, a computer-readable storage medium is provided, on which a computer program is stored. The computer program includes executable instructions. When the executable instructions are executed by a processor, the method described in the above embodiment is implemented.
[0018] According to one aspect of an embodiment of the present application, an electronic device is provided, comprising: one or more processors; and a memory for storing executable instructions of the processors, wherein when the executable instructions are executed by the one or more processors, the one or more processors implement the methods described in the above embodiments.
[0019] The beneficial effects of the present application are as follows: first, by dividing the target cantilever column structure into structural units, a number of structural units connected in series are obtained, and a model of an equal-section structure can be constructed for each of the divided structural units, thereby establishing a finite element model to be tested corresponding to the target cantilever column structure.
[0020] By judging whether the finite element models to be tested meet the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system and the preset vortex-excitation force evaluation system, the target finite element models finally obtained meet the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system and the preset vortex-excitation force evaluation system, so as to equivalently simulate the vortex-excitation vibration that may occur in the target cantilever column structure in the actual scenario.
[0021] Furthermore, by obtaining the structural characteristic parameters of the target cantilever column structure, that is, obtaining the target average wind load that the target finite element model is subjected to under the average wind load evaluation system, and the target vortex excitation force parameters under the vortex excitation force evaluation system, and then performing a coordinate system transformation of the unit transformation matrix on each target structural unit, the overall average wind load and overall vortex excitation force parameters that the target cantilever column structure is subjected to in the preset coordinate system are equivalently simulated.
[0022] Furthermore, by determining the structural vibration acceleration, structural vibration velocity and structural vibration displacement of the target cantilever column structure through structural characteristic parameters, and correcting them according to the preset weak shear flow effect correction coefficient, the influence of weak shear flow in the vertical height direction of the cross section of the cantilever column structure due to the gradual characteristics can be eliminated, thereby improving the analysis precision and accuracy of the vortex vibration problem. At the same time, the constructed dynamic equilibrium equation can predict the vortex vibration response of the target cantilever column structure when it is affected by the incoming wind, and the parameters of the incoming wind can be monitored to predict in real time whether the cantilever column structure will have a vortex vibration response.
[0023] Therefore, through the vortex-oscillation response prediction method of the cantilever column structure provided in the present application, the influence of weak shear flow can be eliminated, and a dynamic equilibrium equation can be established to accurately predict the vortex-oscillation response. BRIEF DESCRIPTION OF THE DRAWINGS
[0024] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following is a brief description of the drawings required for use in the description of the embodiments. Obviously, the drawings described are only part of the embodiments of the present application, not all of the embodiments, and those skilled in the art can also obtain other design solutions and drawings based on these drawings without creative work.
[0025] Figure 1 is a schematic diagram of a target cantilever column structure under the influence of incoming wind according to an embodiment of the present application; Figure 2 It is a flow chart of a method for predicting vortex-vibration response of a cantilever column structure according to an embodiment of the present application; Figure 3 It is a schematic diagram of a single structural unit (left side) and a single equal-section structure (right side) according to an embodiment of the present application; Figure 4 A logic flow chart of the entire vortex vibration response prediction according to an embodiment of the present application; Figure 5 A schematic diagram of a concentrated mass additional node for establishing a finite element model of a cantilever column structure according to an embodiment of the present application; Figure 6 is a block diagram of a vortex-vibration response prediction device for a cantilever column structure according to an embodiment of the present application; Figure 7 It is a schematic diagram of the system structure of an electronic device according to an embodiment of the present application. DETAILED DESCRIPTION
[0026] Example embodiments are now described more fully in conjunction with reference to the accompanying drawings. However, example embodiments can be implemented in a variety of forms and should not be construed as limited to the examples set forth herein; rather, these embodiments are provided so that this application will be more comprehensive and complete and fully convey the concept of the example embodiments to those skilled in the art.
[0027] In addition, described feature, structure or characteristic can be combined in one or more embodiments in any suitable manner. In the following description, many specific details are provided so as to provide a full understanding of the embodiments of the present application. However, it will be appreciated by those skilled in the art that the technical scheme of the present application can be put into practice without one or more of the specific details, or other methods, components, devices, steps, etc. can be adopted. In other cases, known methods, devices, realizations or operations are not shown or described in detail to avoid blurring the various aspects of the present application.
[0028] The block diagrams shown in the accompanying drawings are only functional entities and do not necessarily correspond to physically independent entities. These functional entities can be implemented in software form, or in one or more hardware modules or integrated circuits, or in different networks and / or processor devices and / or micro-control node devices.
[0029] The flowcharts shown in the accompanying drawings are only exemplary and do not necessarily include all the contents and operations / steps, nor must they be executed in the order described. For example, some operations / steps can be decomposed, and some operations / steps can be combined or partially combined, so the actual execution order may change according to actual conditions.
[0030] It should be noted that the "multiple" mentioned in this article refers to two or more. "And / or" describes the association relationship of the associated objects, indicating that there can be three relationships. For example, A and / or B can mean: A exists alone, A and B exist at the same time, and B exists alone. The character " / " generally indicates that the associated objects before and after are in an "or" relationship.
[0031] The following is a detailed introduction to the technical background of the embodiments of the present application: A large number of engineering structures are cantilever column structures with similar and slowly changing cross-sections. Figure 1 As shown, Figure 1 It is a cantilever column structure with a circular cross-section, such as wind turbine towers, chimneys, bridge towers, power steel poles, etc. The cross-section types of these structures are circular, elliptical, polygonal, etc. As the height of the structure increases, the cross-section changes slowly, the size gradually decreases, or gradually increases, or increases and then decreases, or first decreases and then increases. This type of cantilever column structure has a towering appearance, which brings new challenges to safety. The taller and more flexible the structure, the higher the probability of vortex vibration. Although vortex vibration will not directly cause dynamic instability and damage to the structure, vortex vibration with a larger amplitude will reduce the assembly accuracy during construction, increase the difficulty of construction, endanger the safety of construction personnel, affect the safety of the structure during operation, reduce fatigue life, and thus affect the normal use of the structure.
[0032] Since the incoming wind has a significant wind profile, the cross-section of the cantilever column structure has a gradual characteristic in the vertical height direction, and there is a weak shear flow in the vertical direction. These factors have a great influence on the analysis accuracy of the structural vortex-induced vibration problem.
[0033] Therefore, in order to solve the above problems, the embodiment of the present application proposes a vortex-induced vibration response prediction method for a cantilever column structure that is applicable to weak shear flow effects, so as to reduce calculation errors and improve the analysis accuracy of vortex-induced vibration problems.
[0034] The implementation details of the technical solution of the embodiment of the present application are described in detail below: According to one aspect of an embodiment of the present application, a method for predicting vortex-induced vibration response of a cantilever column structure is provided. Figure 2 This is a flow chart of a method for predicting vortex-vibration response of a cantilever column structure according to an embodiment of the present application. The method at least includes steps S1 to S8, which are described in detail as follows: In step S1, the target cantilever column structure is divided into structural units to obtain a plurality of structural units connected in series.
[0035] Specifically, by dividing the target cantilever column structure into structural units, several structural units connected in series can be obtained. Figure 1 Taking the cantilever column structure with a circular cross section as an example, the schematic diagram of the structural unit obtained by division can be shown as follows Figure 3 As shown, it is a cross-sectional gradual change structure with a small circular upper end surface and a large circular lower end surface. i Indicates the upper end surface, Indicates the diameter of the upper end surface, j Indicates the lower end face, Indicates the diameter of the lower end surface, e represents the target surface (any cross section in the structural unit can be selected arbitrarily), It should be noted that the target surface selected arbitrarily may not necessarily meet the requirements of the preset characteristic dimension evaluation system described in the embodiment of the present application. Therefore, it is necessary to set a characteristic dimension 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-sectional structure is reconstructed.
[0036] In step S2, a model of a uniform cross-sectional structure is constructed for each of the structural units to obtain a uniform cross-sectional 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 of the uniform cross-sectional structures.
[0037] Specifically, a model of a uniform cross-sectional structure is constructed for each structural unit. First, it should be noted that since the cross-sectional size of the target surface of each structural unit may be different, the sizes of the constructed multiple uniform cross-sectional structures may also be different. The following is a detailed explanation of the model construction of a uniform cross-sectional structure for a single structural unit: In one embodiment of the application, constructing a model of a uniform cross-section structure for a single structural unit includes the following steps: Acquire any target point of the structural unit in the height direction, and construct a uniform cross-sectional structure having the same height as the structural unit based on the cross section where the target point is located; The shape of each cross section in the equal cross-section structure is consistent with the shape of the cross section where the target point is located.
[0038] When dividing the target cantilever column structure, it is divided into equal heights in the height direction, so the height of each structural unit is consistent, and the cross-sectional size has a slow (gradual) relationship. Therefore, the equal cross-sectional construction of a single structural unit is to construct a cross-sectional structure that is consistent with the height of the structural unit and has the cross-sectional size of the selected cross-sectional size of the target point, that is, Figure 3 Of course, the uniform cross-section structure proposed in the embodiment of the present application is not limited to the cylindrical structure. When the cross-section of the target cantilever column structure is rectangular, polygonal, etc., the uniform cross-section structure constructed at this time is also adjusted accordingly according to the cross-section of the target cantilever column structure. Figure 3 This is merely an exemplary illustration of an equal cross-sectional structure.
[0039] The above method can be used to construct equal-section structures corresponding to each structural unit, and these equal-section structures can be connected in series according to the series connection order of each structural unit to construct a finite element model to be tested. The subsequent steps can be to set the characteristic size discrimination value of the finite element model to be tested, that is, to determine whether the requirements of the preset characteristic size evaluation system are met for each equal-section structure in the finite element model to be tested.
[0040] In step S3, it is determined whether the finite element models to be tested all meet the requirements of a preset characteristic dimension evaluation system, a preset average wind load evaluation system and a preset vortex-induced force evaluation system.
[0041] In one embodiment of the present application, whether the finite element model to be tested meets the requirements of the preset characteristic dimension evaluation system is determined by the following steps: Acquire characteristic size parameters of each of the structural units, wherein the characteristic size parameters include the cross-sectional size of the upper end surface, the cross-sectional size of the lower end surface, and the cross-sectional size of the target surface of the structural unit; Determine the characteristic size discrimination value of each of the structural units according to each of the upper end surface cross-sectional dimensions, each of the lower end surface cross-sectional dimensions and each of the target surface cross-sectional dimensions; If each of the characteristic dimension discrimination values is lower than a preset characteristic dimension threshold, it is determined that the finite element model to be inspected meets the requirements of the preset characteristic dimension evaluation system; Wherein, the characteristic size discrimination value of each of the structural units satisfies the following formula: (1) In the formula, is the characteristic size discrimination value, is the cross-sectional dimension of the upper end face, is the cross-sectional dimension of the lower end face, is the cross-sectional size of the target surface, , , .
[0042] It should be noted that any target point can be any point of the structural unit in the height direction, so the target surface can also be the cross section of any target point. However, the cross-sectional size of the target surface needs to meet the requirements of the preset characteristic size evaluation system, which can be further explained by the following method: The cross-sectional dimensions of the upper end face, the lower end face and the target face of the structural unit correspond to Figure 3 In i End face, j End face and e The cross-sectional dimensions of the upper end face, the lower end face and the target face are respectively , and The preset characteristic size threshold can be set arbitrarily in the range of 0.01~0.05. The cross-sectional size of the structural unit changes slowly (gradually) with the height. When finite element modeling is performed, a gradual uniform cross-section is used to simulate the actual slowly changing structure. The end cross-section of any structural unit is i End face (upper end face), j End face (lower end face), the characteristic dimensions are , ; In the finite element model, the characteristic size (Target surface e The cross-sectional dimensions corresponding to the end face) are constructed to form a uniform cross-sectional structure, so , , , you can refer to the above formula (1) to set the characteristic size discrimination value , the preset feature size threshold is used express, The recommended value is 1%~5%. Only when the requirement of the preset characteristic size evaluation system is met, the characteristic size discrimination value is compared with the preset characteristic size threshold to check whether the constructed uniform cross-section structure can equivalently simulate the structural effect of the slowly varying structure of the structural unit.
[0043] If the requirements of the preset characteristic size evaluation system are not met, it is necessary to reselect the target point, that is, reselect the target surface ( eend face), and reconstruct the equal cross-sectional structure with the cross section of the newly selected target surface to perform the next round of feature size judgment value until the reset result is The value is smaller than the preset feature size threshold.
[0044] Among them, it should be noted that , and It represents the general term for the characteristic size parameters of the target cantilever column structure, which can be applied to cantilever column structures of any cross-sectional size type. The following are examples of the cross-sectional size types of the target cantilever column structure: In one embodiment of the present application, if the cross-sectional dimension type of the target cantilever column structure is a circle, the actual structure corresponding to any structural unit is a truncated cone, and the diameter of the circle is used as the characteristic dimension, and the characteristic dimension parameter is , , Substitute into formula (1) to set the characteristic size discrimination value .
[0045] In one embodiment of the present application, if the cross-sectional dimension type of the target cantilever column structure is a rectangle, the actual structure corresponding to any structural unit is a prism, and the long side of the rectangle is used as the characteristic dimension parameter. L i , L j , Substitute into formula (1) to set the characteristic size discrimination value .
[0046] In one embodiment of the present application, if the cross-sectional dimension type of the target cantilever column structure is a square, the actual structure corresponding to any structural unit is a regular prism, and the side length of the square is used as the characteristic dimension parameter. SL i , SL j , Substitute into formula (1) to set the characteristic size discrimination value .
[0047] In one embodiment of the present application, if the cross-sectional dimension type of the target cantilever column structure is a polygon, the actual structure corresponding to any structural unit is a prism, and the equivalent diameter of the polygon is used as the characteristic dimension parameter. ED i , ED j , Substitute into formula (1) to set the characteristic size discrimination value .
[0048] In one embodiment of the present 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: Obtaining an average wind load parameter received by the upper end surface of each of the structural units, an average wind load parameter received by the lower end surface of each of the structural units, and an average wind load parameter received by the target surface of each of the structural units; Input the average wind load parameters received by each of the upper end faces, the average wind load parameters received by each of the lower end faces, and the average wind load parameters received by each of the target faces into a preset average wind load calculation formula to obtain a first average wind load received by each of the upper end faces, a second average wind load received by each of the lower end faces, and a third average wind load received by each of the target faces; Determining an average wind load discrimination value of each of the structural units according to each of the first average wind loads, each of the second average wind loads, and each of the third average wind loads; If each of the average wind load discrimination values is lower than a preset average wind load threshold, it is determined that the finite element model to be tested meets the requirements of the average wind load evaluation system; The average wind load parameters include air density, drag coefficient, The cross-sectional dimensions of the cross section of the structural unit and any height The average incoming wind speed; Any height The average wind load can be calculated using the following average wind load calculation formula: (2) In the formula, is the air density; is the drag coefficient; For any height The cross-sectional dimensions of the cross section of the structural unit at For any height The average incoming wind speed; For any height The average wind load on the space.
[0049] Specifically, when the incoming wind acts on the target cantilever column structure, the target cantilever column structure is at any height along the vertical direction. The average wind load The expression of can be referred to the above formula (2). According to the above formula (2), in any structural unit, i The height, average wind speed and average wind load corresponding to the end face are: , , , j The height, average wind speed and average wind load corresponding to the end face are: , , , target surface e The height, average wind speed and average wind load corresponding to the end face are: , , , and thus set the average wind load judgment value : (3) The preset average wind load threshold is used express, The recommended value is 1%~5%. By comparing the average wind load discrimination value with the preset average wind load threshold, it can be tested whether the equal-section structure can equivalently simulate the wind load effect of the cantilever column structure in the actual scene.
[0050] The distribution of the average incoming wind speed along the height direction obeys the power law or logarithmic law and is expressed by an exponential model or logarithmic model.
[0051] First, when the exponential model is used to calculate the target average wind load, any height Average wind speed The expression is (4) In the formula, is the reference height; for Reference wind speed at is the surface roughness coefficient.
[0052] When the logarithmic model is used to calculate the target average wind load, any height Average wind speed The expression is: (5) In the formula, is the von Karman constant, which is generally taken as 0.4; is the zero plane displacement, which is a function of the properties, height and distribution of the surface roughness elements. The zero plane displacement can be understood as the difference between the zero wind speed height and the roughness height. For urban landforms, , is the ground roughness length, is the average height of the buildings; is the friction wind speed.
[0053] make , , , when the exponential model is used to calculate the average wind load, , , , at this time the average wind load judgment value The expression is: (6) Setting the average wind load threshold , The recommended value is 1%~5%. The requirements are met only when the equal-section structure is used, and it can be verified whether the wind load effect of the cantilever column structure in the actual scene can be equivalently simulated.
[0054] When the logarithmic model is used to calculate the average wind load, , , , at this time the average wind load judgment value The expression is: (7) Setting the average wind load threshold , The recommended value is 1%~5%. The requirements are met only when the equal-section structure is used, and it can be verified whether the wind load effect of the cantilever column structure in the actual scene can be equivalently simulated.
[0055] In one embodiment of the present 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: Acquire a first vortex-induced force parameter of an upper end surface of each of the structural units, a second vortex-induced force parameter of a lower end surface of each of the structural units, and a third vortex-induced force parameter of a target surface of each of the structural units; Determine the vortex excitation force discrimination value of each of the structural units according to each of the first vortex excitation force parameters, each of the second vortex excitation force parameters, and each of the third vortex excitation force parameters; If each of the vortex excitation force discrimination values is lower than a preset vortex excitation force threshold, it is determined that the finite element model to be tested meets the requirements of the vortex excitation force evaluation system.
[0056] Specifically, the first vortex excitation force parameter, the second vortex excitation force parameter and the third vortex excitation force parameter can all be obtained by the following formula (8): (8) In the formula, is the vortex excitation parameter; is the structural vibration displacement; is the structural vibration velocity; represents aerodynamic parameters; M , N is the order of the Taylor expanded polynomial; is the vortex shedding frequency; is the initial phase; is the parameter to be identified, is the converted circular frequency, is the natural circular frequency of the structure.
[0057] In any structural unit, the end section of the structural unit ( i End face, j The vortex force parameters corresponding to the end surface are , The vortex force parameters at the same position of the characteristic size parameters of the actual structural unit and the medium cross-section structure in the finite element model are taken as ,make , , , set the vortex force discrimination value : (9) Setting the vortex force threshold , The recommended value is 1%~5%. The requirements are met only when the equal-section structure is tested to see whether it can simulate the vortex-induced force effect of the cantilever column structure in actual scenarios.
[0058] Commonly used vortex excitation models include Scanlan empirical linear model, Scanlan empirical nonlinear model, and fifth-order polynomial refined simplified model.
[0059] When the vortex excitation force model adopts the Scanlan empirical linear model, the expression of the target vortex excitation force parameter is: (10) From formula (10), we can get , , , substituting into formula (9), we can get the vortex force discrimination value : (11) In the formula, , , , .
[0060] Setting the vortex force threshold , The recommended value is 1%~5%. The requirements are met only when the equal-section structure is tested to see whether it can simulate the vortex-induced force effect of the cantilever column structure in actual scenarios.
[0061] When the vortex excitation force model adopts the Scanlan empirical nonlinear model, the expression of the target vortex excitation force parameter is: (12) From formula (12), we can get , , , substituting into formula (9), we can get the vortex force discrimination value : (13) Setting the vortex force threshold , The recommended value is 1%~5%. The requirements are met only when the equal-section structure is tested to see whether it can simulate the vortex-induced force effect of the cantilever column structure in actual scenarios.
[0062] When the vortex excitation force model is refined and simplified using a fifth-order polynomial, the expression of the target vortex excitation force parameter is: (14) Among them, from formula (14) we can get , , Substituting it into equation (9), we get the vortex force discrimination value: : (15) Setting the vortex force threshold , The recommended value is 1%~5%. The requirements are met only when the equal-section structure is tested to see whether it can simulate the vortex-induced force effect of the cantilever column structure in actual scenarios.
[0063] In step S4, if the finite element model to be tested does not meet the requirements of the characteristic size evaluation system or the average wind load evaluation system or the vortex excitation force evaluation system, 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 excitation force evaluation system, and the finite element model to be tested is determined to be the target finite element model.
[0064] Specifically, please refer to Figure 4 As shown, Figure 4It is a logical flow chart for the entire vortex-induced vibration response prediction. It is necessary for the finite element models to be tested to meet the requirements of the preset characteristic dimension evaluation system, the preset average wind load evaluation system and the preset vortex-induced force evaluation system. That is, the three conditions that each of the characteristic dimension judgment values is lower than the preset characteristic dimension threshold, each of the average wind load judgment values is lower than the preset average wind load threshold and each of the vortex-induced force judgment values is lower than the preset vortex-induced force threshold are met. Only in this way can the finite element model be able to fully simulate the target cantilever column structure.
[0065] In step S5, the structural characteristic parameters of the target cantilever column structure are obtained, and the target average wind load received by the target finite element model under the average wind load evaluation system and the target vortex excitation force parameters under the vortex excitation force evaluation system are determined based on the structural characteristic parameters.
[0066] Specifically, since the average wind load parameters and vortex excitation parameters of each structural unit can be obtained in the above steps, the target average wind load received by the target finite element model under the average wind load evaluation system can be determined. , and the target vortex force parameters under the vortex force evaluation system .
[0067] In step S6, the coordinate system of the unit transformation matrix is transformed for each target structural unit in the target finite element model 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 on the target cantilever column structure in the preset coordinate system.
[0068] Specifically, the unit force in the local coordinate system is converted into the force in the global coordinate system through the unit conversion matrix, that is, the overall average wind load on the target cantilever column structure. and the target cantilever column structure , the dynamic equilibrium equation is established, and the vortex vibration response of the structure is solved by numerical analysis method.
[0069] The kinetic equilibrium equation is: (16) In the formula, m , c , k They are the average mass, damping and stiffness per linear meter of the cantilever column structure respectively.
[0070] 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, and a dynamic equilibrium equation corresponding to the target cantilever column structure is constructed according to the structural vibration acceleration, the structural vibration velocity, the structural vibration displacement, the overall average wind load, the overall vortex excitation parameter and a preset weak shear flow effect correction coefficient.
[0071] In one embodiment of the present application, the structural characteristic parameters include the average mass per linear meter, the average damping per linear meter and the average stiffness per linear meter of the cantilever column structure; the dynamic equilibrium equation is specifically the above formula (16), where, is the average mass per linear meter of the target cantilever column structure, is the average damping per linear meter of the target cantilever column structure, is the average stiffness per linear meter of the target cantilever column structure; is the structural vibration acceleration, is the structural vibration velocity, is the structural vibration displacement; is the overall average wind load, is the overall vortex force parameter, is the correction coefficient of the vortex-induced weak shear flow effect.
[0072] Specifically, the cantilever column is affected by weak shear flow in the vertical direction, which will greatly reduce the accuracy of the constructed finite element model in equivalent simulation of the cantilever column. Therefore, the vortex-induced weak shear flow effect correction coefficient can be used. To correct, is a function of time, frequency, and axial coordinates, It can be obtained through experiments or numerical simulation methods. Among them, the experiment can use a series of narrow-band aeroelastic models to synchronously test the vortex-induced forces of each narrow band and obtain it through correlation analysis; it can also be obtained by testing the vortex-induced force parameters of different sections and through experimental analysis.
[0073] In step S8, a vortex-oscillation response prediction is performed on the target cantilever column structure when it is affected by incoming wind according to the dynamic equilibrium equation.
[0074] Specifically, the constructed dynamic equilibrium equation can accurately predict the vortex-oscillation response of the cantilever column structure according to the specific parameters of the incoming wind.
[0075] In one embodiment of the present application, Figure 5As shown, the finite element model uses finite element software to establish a concentrated mass model. Taking the tower structure as an example, several units are divided from the top to the bottom elevation according to a fixed length. The cantilever column is established as a finite element model of several units in series using the "sugar-coated haws string" model idea. In terms of unit division, a gradual equal cross-section is used to simulate the actual slowly varying structure. The model unit is a beam unit, and the cross-sectional size of the tower changes with height by defining different cross-sectional areas (cross-sectional dimensions) and moments of inertia (integrals related to area and size) for each unit, simplifying the tower into an additional concentrated mass added to the node at the top. superior.
[0076] In summary, the embodiments of the present application divide the target cantilever column structure into structural units to obtain a number of structural units connected in series, and can construct a model of an equal-section structure for each of the divided structural units, thereby establishing a finite element model to be tested corresponding to the target cantilever column structure.
[0077] By judging whether the finite element models to be tested meet the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system and the preset vortex-excitation force evaluation system, the target finite element model finally obtained meets the requirements of the preset characteristic size evaluation system, the preset average wind load evaluation system and the preset vortex-excitation force evaluation system, thereby equivalently simulating the vortex-excitation vibration that may occur in the target cantilever column structure in the actual scenario.
[0078] 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 to which the target finite element model is subjected under the average wind load evaluation system, and the target vortex excitation force parameters under the vortex excitation force evaluation system, and then performing a coordinate system transformation of the unit transformation matrix on each target structural unit, the overall average wind load and overall vortex excitation force parameters to which the target cantilever column structure is subjected in the preset coordinate system are equivalently simulated.
[0079] Furthermore, by determining the structural vibration acceleration, structural vibration velocity and structural vibration displacement of the target cantilever column structure through structural characteristic parameters, and correcting them according to the preset weak shear flow effect correction coefficient, the influence of weak shear flow in the vertical height direction of the cross section of the cantilever column structure due to the gradual characteristics can be eliminated, thereby improving the analysis precision and accuracy of the vortex vibration problem. At the same time, the constructed dynamic equilibrium equation can predict the vortex vibration response of the target cantilever column structure when it is affected by the incoming wind, and the parameters of the incoming wind can be monitored to predict in real time whether the cantilever column structure will have a vortex vibration response.
[0080] Therefore, through the vortex-oscillation response prediction method of the cantilever column structure provided in the present application, the influence of weak shear flow can be eliminated, and a dynamic equilibrium equation can be established to accurately predict the vortex-oscillation response.
[0081] Figure 6 This is a block diagram of a vortex-induced vibration response prediction device 300 for a cantilever column structure according to an embodiment of the present application. According to a vortex-induced vibration response prediction device 300 for a cantilever column structure according to an embodiment of the present application, the device 300 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. A division unit 301 is used to divide the target cantilever column structure into structural units to obtain a plurality of structural units connected in series; The construction unit 302 is used to construct a model of a uniform cross-section structure for each of the structural units, obtain a uniform 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 according to each of the uniform cross-section structures; A judgment unit 303 is used to judge whether the finite element model to be tested meets the requirements of a preset characteristic size evaluation system, a preset average wind load evaluation system and a preset vortex excitation force evaluation system; The iteration unit 304 is used for repeating the above steps if the finite element model to be tested does not meet the requirements of the characteristic size evaluation system or the average wind load evaluation system or the vortex excitation 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 excitation force evaluation system, and determining that the finite element model to be tested is a target finite element model; An acquisition unit 305 is used to acquire structural characteristic parameters of the target cantilever column structure, and determine the target average wind load received by the target finite element model under the average wind load evaluation system and the target vortex excitation force parameters under the vortex excitation force evaluation system according to the structural characteristic parameters; The conversion unit 306 is used to perform a coordinate system conversion of a unit conversion matrix on each target structural unit in the target finite element model according to the target average wind load and the target vortex excitation force parameter, so as to obtain the overall average wind load and the overall vortex excitation force parameter on the target cantilever column structure in a preset coordinate system; A 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 a dynamic balance equation corresponding to the target cantilever column structure according to the structural vibration acceleration, the structural vibration velocity, the structural vibration displacement, the overall average wind load, the overall vortex excitation parameter and a preset weak shear flow effect correction coefficient; The prediction unit 308 is used to predict the vortex-oscillation response of the target cantilever column structure when it is affected by the incoming wind according to the dynamic equilibrium equation.
[0082] As another aspect, the present application further provides a computer-readable storage medium on which a program product capable of implementing the method provided above in this specification is stored. In some possible implementations, various aspects of the present application may also be implemented in the form of a program product, which includes a program code, and when the program product is run on a terminal device, the program code is used to cause the terminal device to execute the steps according to various exemplary implementations of the present application described in the above "Embodiment Method" section of this specification.
[0083] According to the program product for implementing the above method in the embodiment of the present application, it can adopt a portable compact disk read-only memory (CD-ROM) and include program code, and can be run on a terminal device, such as a personal computer. However, the program product of the present application is not limited thereto. In this document, a readable storage medium can be any tangible medium containing or storing a program, which can be used by or in combination with an instruction execution system, an apparatus or a device.
[0084] The program product may use any combination of one or more readable media. The readable medium may be a readable signal medium or a readable storage medium. The readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared or semiconductor system, device or device, or any combination of the above. More specific examples of readable storage media (a non-exhaustive list) include: an electrical connection with one or more wires, a portable disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above.
[0085] Computer readable signal media may include a data signal propagated in baseband as part of a carrier wave, wherein readable program code is carried. Such propagated data signals may take a variety of forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. Readable signal media may also be any readable medium other than a readable storage medium, which may send, propagate, or transmit a program for use by or in conjunction with an instruction execution system, apparatus, or device.
[0086] The program code embodied on the readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wired, optical cable, RF, etc., or any suitable combination of the foregoing.
[0087] Program code for performing the operations of the present application may be written in any combination of one or more programming languages, including object-oriented programming languages, such as Java, C++, etc., and conventional procedural programming languages, such as "C" or similar programming languages. The program code may be executed entirely on the user computing device, partially on the user device, as a separate software package, partially on the user computing device, partially on the remote computing device, or entirely on the remote computing device or server. In the case of a remote computing device, the remote computing device may be connected to the user computing device through any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computing device (e.g., using an Internet service provider to connect through the Internet).
[0088] As another aspect, the present application also provides an electronic device capable of implementing the above method.
[0089] Those skilled in the art will appreciate that various aspects of the present application may be implemented as a system, method or program product. Therefore, various aspects of the present application may be specifically implemented in the following forms, that is, a complete hardware implementation, a complete software implementation (including firmware, microcode, etc.), or a combination of hardware and software, which may be collectively referred to as "circuit", "module" or "system" herein.
[0090] Refer to the following Figure 7 The electronic device 400 according to this embodiment of the present application is described. Figure 7 The electronic device 400 shown is merely an example and should not bring any limitation to the functions and scope of use of the embodiments of the present application.
[0091] like Figure 7As shown, the electronic device 400 is in the form of a general 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 the storage unit 420 and the processing unit 410).
[0092] The storage unit stores program codes, which can be executed by the processing unit 410, so that the processing unit 410 executes the steps described in the above “Example Method” section of this specification according to various exemplary implementations of the present application.
[0093] The storage unit 420 may include a readable medium in the form of a volatile storage unit, such as a random access memory unit (RAM) 421 and / or a cache memory unit 422 , and may further include a read-only memory unit (ROM) 423 .
[0094] The 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 which or some combination may include an implementation of a network environment.
[0095] Bus 430 may represent one or more of several types of bus structures, including a memory unit bus or memory unit controller node, a peripheral bus, an accelerated graphics port, a processing unit, or a local bus using any of a variety of bus architectures.
[0096] The electronic device 400 may also communicate with one or more external devices 1200 (e.g., keyboards, pointing devices, Bluetooth devices, etc.), may also communicate with one or more devices that enable a user to interact with the electronic device 400, and / or communicate with any device that enables the electronic device 400 to communicate with one or more other computing devices (e.g., routers, modems, etc.). Such communication may be performed via an input / output (I / O) interface 450. In addition, the electronic device 400 may also communicate with one or more networks (e.g., local area networks (LANs), wide area networks (WANs), and / or public networks, such as the Internet) via a network adapter 460. As shown, the network adapter 460 communicates with other modules of the electronic device 400 via a bus 430. It should be understood that, although not shown in the figure, other hardware and / or software modules may be used in conjunction with the 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, etc.
[0097] Through the description of the above implementation methods, it is easy for those skilled in the art to understand that the example implementation methods described here can be implemented by software, or by combining software with necessary hardware. Therefore, the technical solution according to the implementation method of the present application can be embodied in the form of a software product, which can be stored in a non-volatile storage medium (which can be a CD-ROM, a USB flash drive, a mobile hard disk, etc.) or on a network, including several instructions to enable a computing device (which can be a personal computer, a server, a terminal device or a network device, etc.) to execute the method according to the implementation method of the present application.
[0098] In addition, the above-mentioned figures are only schematic illustrations of the processes included in the method according to the exemplary embodiment of the present application, and are not intended to be limiting. It is easy to understand that the processes shown in the above-mentioned figures do not indicate or limit the time sequence of these processes. In addition, it is also easy to understand that these processes can be performed synchronously or asynchronously in multiple modules.
[0099] It should be understood that the present application is not limited to the precise structures that have been described above and shown in the drawings, and that various modifications and changes may be performed without departing from the scope thereof. The scope of the present application is limited only by the appended claims.
Claims
1. A method for predicting vortex-induced vibration response of a cantilever column structure, characterized in that: The method comprises: S1, dividing the target cantilever column structure into structural units to obtain a number of structural units connected in series; S2, constructing a model of a uniform cross-section structure for each of the structural units, obtaining a uniform cross-section structure corresponding to each of the structural units, and establishing a finite element model to be tested corresponding to the target cantilever column structure according to each of the uniform cross-section structures; S3, judging whether the finite element models to be tested all meet the requirements of a preset characteristic size evaluation system, a preset average wind load evaluation system and a preset vortex excitation force evaluation system; S4, if the finite element model to be tested does not meet the requirements of the characteristic size evaluation system or the average wind load evaluation system or the vortex excitation force evaluation system, 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 excitation force evaluation system, and the finite element model to be tested is determined to be the target finite element model; S5, acquiring structural characteristic parameters of the target cantilever column structure, and determining the target average wind load received by the target finite element model under the average wind load evaluation system and the target vortex excitation force parameters under the vortex excitation force evaluation system according to the structural characteristic parameters; S6, performing a coordinate system transformation of a unit transformation matrix on each target structural unit in the target finite element model according to the target average wind load and the target vortex excitation force parameter, to obtain an overall average wind load and an overall vortex excitation force parameter on the target cantilever column structure in a preset coordinate system; S7, determining the structural vibration acceleration, structural vibration velocity and structural vibration displacement of the target cantilever column structure according to the structural characteristic parameters, and constructing a dynamic equilibrium equation corresponding to the target cantilever column structure according to the structural vibration acceleration, the structural vibration velocity, the structural vibration displacement, the overall average wind load, the overall vortex excitation parameter and a preset weak shear flow effect correction coefficient; S8, predicting the vortex-oscillation response of the target cantilever column structure when it is affected by incoming wind according to the dynamic equilibrium equation.
2. The vortex-induced vibration response prediction method of a cantilever column structure according to claim 1, characterized in that: The model construction of a single structural unit with a uniform cross-section structure comprises the following steps: Acquire any target point of the structural unit in the height direction, and construct a uniform cross-sectional structure having the same height as the structural unit based on the cross section where the target point is located; The shape of each cross section in the equal cross-section structure is consistent with the shape of the cross section where the target point is located.
3. The method for predicting vortex-induced vibration response of a cantilever column structure according to claim 2, characterized in that: Whether the finite element model to be tested meets the requirements of the preset characteristic dimension evaluation system is determined by the following steps: Acquire characteristic size parameters of each of the structural units, wherein the characteristic size parameters include the cross-sectional size of the upper end surface, the cross-sectional size of the lower end surface, and the cross-sectional size of the target surface of the structural unit, wherein the cross-sectional size of the target surface corresponds to the size of the cross section where the target point is located; Determine the characteristic size discrimination value of each of the structural units according to each of the upper end surface cross-sectional dimensions, each of the lower end surface cross-sectional dimensions and each of the target surface cross-sectional dimensions; If each of the characteristic dimension discrimination values is lower than a preset characteristic dimension threshold, it is determined that the finite element model to be inspected meets the requirements of the preset characteristic dimension evaluation system; Wherein, the characteristic size discrimination value of each of the structural units satisfies the following formula: In the formula, is the characteristic size discrimination value, is the cross-sectional dimension of the upper end surface, is the cross-sectional dimension of the lower end face, is the cross-sectional size of the target surface, , , .
4. The vortex-induced vibration response prediction method of a cantilever column structure according to claim 3 is 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: Obtaining an average wind load parameter received by the upper end surface of each of the structural units, an average wind load parameter received by the lower end surface of each of the structural units, and an average wind load parameter received by the target surface of each of the structural units; Input the average wind load parameters received by each of the upper end faces, the average wind load parameters received by each of the lower end faces, and the average wind load parameters received by each of the target faces into a preset average wind load calculation formula to obtain a first average wind load received by each of the upper end faces, a second average wind load received by each of the lower end faces, and a third average wind load received by each of the target faces; Determining an average wind load discrimination value of each of the structural units according to each of the first average wind loads, each of the second average wind loads, and each of the third average wind loads; If each of the average wind load discrimination values is lower than a preset average wind load threshold, it is determined that the finite element model to be tested meets the requirements of the average wind load evaluation system; The average wind load parameters include air density, drag coefficient, The cross-sectional dimensions of the cross section of the structural unit and any height The average incoming wind speed; Any height The average wind load can be calculated using the following average wind load calculation formula: In the formula, is the air density; is the drag coefficient; For any height The cross-sectional dimensions of the cross section of the structural unit at For any height The average incoming wind speed; For any height The average wind load on the space.
5. The method for predicting vortex-induced vibration response of a cantilever column structure according to claim 4, characterized in that: The average wind load discrimination value can be obtained by the following formula: In the formula, is the average wind load discrimination value, is the first average wind load, is the second average wind load, is the third average wind load, is the average incoming wind speed at the height of the upper end surface, is the average incoming wind speed at the height of the lower end surface, It is the average incoming wind speed at the height of the target surface.
6. The method for predicting vortex-induced vibration response of a cantilever column structure according to claim 5, characterized in that: Whether the finite element model to be tested meets the requirements of the vortex excitation force evaluation system is determined by the following steps: Acquire a first vortex-induced force parameter of an upper end surface of each of the structural units, a second vortex-induced force parameter of a lower end surface of each of the structural units, and a third vortex-induced force parameter of a target surface of each of the structural units; Determine the vortex excitation force discrimination value of each of the structural units according to each of the first vortex excitation force parameters, each of the second vortex excitation force parameters, and each of the third vortex excitation force parameters; If each of the vortex excitation force discrimination values is lower than a preset vortex excitation force threshold, it is determined that the finite element model to be tested meets the requirements of the vortex excitation force evaluation system.
7. The method for predicting vortex-induced vibration response of a cantilever column structure according to claim 6, characterized in that: The structural characteristic parameters include the average mass per linear meter, the average damping per linear meter and the average stiffness per linear meter of the target cantilever column structure; the dynamic equilibrium equation is specifically the following formula: In the formula, is the average mass per linear meter of the target cantilever column structure, is the average damping per linear meter of the target cantilever column structure, is the average stiffness per linear meter of the target cantilever column structure; is the structural vibration acceleration, is the structural vibration velocity, is the structural vibration displacement; is the overall average wind load, is the overall vortex force parameter, is the correction coefficient of the vortex-induced weak shear flow effect.
8. A vortex-vibration response prediction device for a cantilever column structure, characterized in that: The device comprises: A division unit is used to divide the target cantilever column structure into structural units to obtain a plurality of structural units connected in series; A construction unit, used to construct a model of a uniform cross-section structure for each of the structural units, obtain a uniform 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 according to each of the uniform cross-section structures; A judgment unit, used to judge whether the finite element models to be tested all meet the requirements of a preset characteristic size evaluation system, a preset average wind load evaluation system and a preset vortex excitation force evaluation system; an iterative unit, for repeating the above steps if the finite element model to be tested does not meet the requirements of the characteristic size evaluation system or the average wind load evaluation system or the vortex excitation 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 excitation force evaluation system, and determining that the finite element model to be tested is a target finite element model; an acquisition unit, used for acquiring structural characteristic parameters of the target cantilever column structure, and determining the target average wind load received by the target finite element model under the average wind load evaluation system and the target vortex excitation force parameters under the vortex excitation force evaluation system according to the structural characteristic parameters; A conversion unit, used for performing a coordinate system conversion of a unit conversion matrix on each target structural unit in the target finite element model according to the target average wind load and the target vortex excitation force parameter, so as to obtain the overall average wind load and the overall vortex excitation force parameter on the target cantilever column structure in a preset coordinate system; a dynamic balance unit, for determining the structural vibration acceleration, structural vibration velocity and structural vibration displacement of the target cantilever column structure according to the structural characteristic parameters, and constructing a dynamic balance equation corresponding to the target cantilever column structure according to the structural vibration acceleration, the structural vibration velocity, the structural vibration displacement, the overall average wind load, the overall vortex excitation force parameter and a preset weak shear flow effect correction coefficient; The prediction unit is used to predict the vortex-oscillation response of the target cantilever column structure when it is affected by the incoming wind according to the dynamic equilibrium equation.
9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores at least one program code, and the at least one program code is loaded and executed by a processor to implement the operations performed by the method according to any one of claims 1 to 7.
10. An electronic device, characterized in that: The electronic device includes one or more processors and one or more memories, wherein the one or more memories store at least one program code, and the at least one program code is loaded and executed by the one or more processors to implement the operations performed by the method according to any one of claims 1 to 7.
Citation Information
Patent Citations
Bridge vortex-induced vibration amplitude prediction method based on decision tree and recurrent neural network
CN113642068A
Vortex excitation fluid force data calculation program for linear structure, design program for linear structure considering vortex excitation and vortex excitation fluid force test method for linear structure
JP2014106552A
Enhancement of vortex induced forces and motion through surface roughness control
US20090114001A1