A method, apparatus and device for designing a ring-shaped tuned liquid column damper structure
By conducting finite element numerical simulation and modal analysis of the wind turbine tower, the structural parameters of the annular tuned liquid column damper were determined, solving the problem of the inability to accurately tune the liquid frequency and tower frequency in the existing technology, achieving multi-directional vibration reduction effect and meeting the needs of engineering applications.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2023-03-10
- Publication Date
- 2026-04-07
AI Technical Summary
Existing technologies cannot effectively determine the optimal structural parameters of the annular tuned liquid column damper, resulting in the inability to accurately tune the liquid frequency and tower frequency, and thus failing to achieve multi-directional vibration reduction effects.
By establishing a finite element numerical simulation model of the wind turbine tower, modal analysis was performed to determine the first-order bending natural frequency and structural modal mass of the tower. The tuning frequency of the annular tuned liquid column damper was calculated, the liquid column length was initialized, and the liquid depth and other structural parameters were determined through numerical relationships. Bidirectional fluid-structure interaction numerical calculations were performed to optimize the damper structure.
The design accuracy and efficiency of the annular tuned liquid column damper have been improved, ensuring that it can effectively reduce vibration in multiple directions in wind power generation systems and meet the needs of engineering practice.
Smart Images

Figure CN116401778B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of vibration control in wind power generation systems, specifically to a method, apparatus, and equipment for designing a ring-shaped tuned liquid column damper. Background Technology
[0002] Under the influence of external environmental loads and operational loads, the tower of a wind power generation system will experience cyclic vibrations. As the capacity of wind turbine generators increases, the height of the tower also gradually increases. The vibration of the entire structure—comprising the turbine, tower, foundation, and ground—has become a crucial indicator for structural design and safe, stable operation. Installing tuned damping devices on the tower can effectively suppress the vibration of the turbine and tower, ensuring their safety. Commonly used tuned damping devices include tuned mass dampers (TMDs), tuned liquid dampers (TLDs), and tuned liquid column dampers (TLCDs). TLCDs typically consist of a U-shaped rectangular water tank filled with liquid. By adjusting the liquid length, the TLCD frequency is made close to the frequency of the controlled object. During vibration, energy dissipation is achieved through the liquid head loss caused by the liquid movement and the viscous action in the boundary layer. TLCD can only control the vibration response of a structure in a single direction. However, the input directions of the most important loads affecting wind turbine generators, such as wind loads, wave loads, and seismic loads, are uncertain, making it impossible to achieve multi-directional vibration reduction and energy dissipation for wind turbine generators.
[0003] A new type of annular tuned liquid column damper has been proposed in recent years, with its specific structure detailed in documents CN114909010A and CN217501237U. The annular tuned liquid column damper is a novel damping and vibration reduction device capable of multi-directional tuning. Its geometric characteristics are consistent with those of a circular tower, offering the advantage of simultaneous multi-directional tuning and adapting to the random, multi-directional loads borne by the entire wind turbine structure. However, due to the greater complexity of the internal structure of the annular tuned liquid column damper compared to a TLCD (Turbocooled Liquid Column Damper), traditional design methods used for TLCDs are insufficient to meet the design requirements of the annular tuned liquid column damper. Existing TLCD design methods cannot determine the optimal structural parameters of the annular tuned liquid column damper, such as the outer shell radius and inner shell radius, nor can they ensure precise tuning between the liquid frequency and the tower frequency. Therefore, a new structural design method is urgently needed to determine the structural parameters of the annular tuned liquid column damper. Summary of the Invention
[0004] In view of this, the present invention provides a method, apparatus and equipment for designing a ring-shaped tuned liquid column damper, which can obtain structural parameters that make the ring-shaped tuned liquid column damper have better vibration reduction effect.
[0005] According to a first aspect, embodiments of the present invention provide a structural design method for an annular tuned liquid column damper, the method comprising: establishing a finite element numerical simulation model of a wind turbine tower; performing modal analysis on the finite element numerical simulation model to determine the first-order bending natural frequency and the corresponding structural modal mass of the wind turbine tower; calculating the tuning frequency of the annular tuned liquid column damper using a preset mass ratio and the first-order bending natural frequency, wherein the preset mass ratio is the ratio of the liquid mass in the annular tuned liquid column damper to the structural modal mass; determining a first strip liquid column and a second strip liquid column from the vertical liquid column of the annular tuned liquid column damper, wherein the first strip liquid column is a liquid column opposite to the first strip liquid column along the excitation direction on the outer shell diameter, and the second strip liquid column is a relative liquid column in the annular tuned liquid column damper adjacent to and on the same side of the first strip liquid column; initializing a first horizontal length and a second horizontal length between the first strip liquid column. The second horizontal length between the liquid columns is determined; a first numerical relationship is established between the liquid depth, the first horizontal length, the second horizontal length, and the tuning frequency in the annular tuned liquid column damper, and the liquid depth is determined based on the first numerical relationship; a second numerical relationship is established between the liquid mass, the liquid depth, the first horizontal length, the second horizontal length, and other structural parameters, and the other structural parameters are determined based on the second numerical relationship, including the outer shell radius, inner shell radius, horizontal channel height, inner shell height, and outer shell height of the annular tuned liquid column damper; the first horizontal length, the second horizontal length, and the other structural parameters are used as damper structural parameters to perform bidirectional fluid-structure interaction numerical calculations of the wind turbine tower-annular tuned liquid column damper to obtain vibration reduction indexes; it is determined whether the vibration reduction indexes meet the preset vibration reduction conditions, and if they do, the damper structural parameters are output.
[0006] Optionally, the first numerical relationship between the liquid depth, the first horizontal length, the second horizontal length, and the tuning frequency in the generated annular tuned liquid column damper is as follows:
[0007]
[0008] In the formula, χ represents the frequency correction coefficient, with a default value of 1, and η represents the cross-sectional area ratio, with a default value of 1. H1 represents the first horizontal length, H2 represents the second horizontal length, and f f The tuning frequency is represented by , and g represents the gravitational acceleration.
[0009] Optionally, before determining the other structural parameters based on the second numerical relationship, the method further includes: calculating the initial length ratio p using the following formula.
[0010]
[0011] In the formula, H1 represents the first horizontal length, and V represents the liquid depth;
[0012] Determine whether the initial length ratio falls within the range of [0.5, 0.6].
[0013] If it does not fall within the range of [0.5, 0.6], then the first horizontal length and the second horizontal length are reinitialized.
[0014] Optionally, the established second numerical relationship includes: m f =8A V V+A H H1+2A H H2, In the formula, m f A represents the mass of the liquid. V Let A represent the cross-sectional area of a certain strip-shaped liquid column. H R1 and R2 represent the cross-sectional area of the corresponding equivalent TLCD horizontal segment liquid column, respectively, B3 represents the outer shell radius and inner shell radius of the annular tuned liquid column damper, and V represents the liquid depth.
[0015] The step of determining the other structural parameters based on the second numerical relationship includes: substituting the liquid mass, the liquid depth, the first horizontal length, and the second horizontal length into the second numerical relationship to calculate the outer shell radius, inner shell radius, and horizontal channel height of the annular tuned liquid column damper; determining the outer shell height of the annular tuned liquid column damper by a preset multiple of the liquid depth; and calculating the inner shell height of the annular tuned liquid column damper based on the difference between the outer shell height and the horizontal channel height.
[0016] Optionally, the method further includes: creating a hydrodynamic model using the damper structural parameters; adding the liquid depth to a preset compensation to obtain a compensated liquid depth; inputting the compensated liquid depth into the hydrodynamic model to perform free decay vibration simulation to find the simulated liquid depth corresponding to the tuning frequency; updating the value of the liquid depth using the simulated liquid depth; and correcting the damper structural parameters using the relationship between the updated liquid depth and the liquid mass.
[0017] Optionally, the method further includes: substituting the updated liquid depth into the formula. In the formula, a new frequency correction coefficient is output; where χ' represents the new frequency correction coefficient, η represents the cross-sectional area ratio, and η is preset to 1. H1 represents the first horizontal length, H2 represents the second horizontal length, and f f The tuning frequency is represented by , and g represents the gravitational acceleration.
[0018] Optionally, the preset compensation is equal to the horizontal channel height.
[0019] According to a second aspect, embodiments of the present invention provide a structural design device for an annular tuned liquid column damper. The device includes: a tower simulation module for establishing a finite element numerical simulation model of a wind turbine tower, and determining the first-order bending natural frequency and corresponding structural modal mass of the wind turbine tower by performing modal analysis on the finite element numerical simulation model; a tuning frequency calculation module for calculating the tuning frequency of the annular tuned liquid column damper using a preset mass ratio and the first-order bending natural frequency, wherein the preset mass ratio is the ratio of the liquid mass in the annular tuned liquid column damper to the structural modal mass; an equivalent liquid column determination module for determining a first strip liquid column and a second strip liquid column from the vertical liquid column of the annular tuned liquid column damper, wherein the first strip liquid column is a liquid column opposite to the first strip liquid column along the excitation direction on the outer shell diameter, and the second strip liquid column is a relative liquid column in the annular tuned liquid column damper that is adjacent to the first strip liquid column and on the same side; and a horizontal length initialization module for initializing the first horizontal length between the first strip liquid columns and the second horizontal length between the second strip liquid columns. The system includes: a second horizontal length between the liquid columns; a liquid depth calculation module for establishing a first numerical relationship between the liquid depth, the first horizontal length, the second horizontal length, and the tuning frequency in the annular tuned liquid column damper, and determining the liquid depth based on the first numerical relationship; a other structural parameter calculation module for establishing a second numerical relationship between the liquid mass, the liquid depth, the first horizontal length, the second horizontal length, and other structural parameters, and determining the other structural parameters based on the second numerical relationship, wherein the other structural parameters include the outer shell radius, inner shell radius, horizontal channel height, inner shell height, and outer shell height of the annular tuned liquid column damper; a vibration reduction experiment module for performing bidirectional fluid-structure interaction numerical calculations of the wind turbine tower-annular tuned liquid column damper using the first horizontal length, the second horizontal length, and the other structural parameters as damper structural parameters to obtain vibration reduction indicators; and a structural parameter output module for determining whether the vibration reduction indicators meet preset vibration reduction conditions, and if so, outputting the damper structural parameters.
[0020] According to a third aspect, embodiments of the present invention provide a device for designing a ring-shaped tuned liquid column damper structure, comprising: a memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the computer instructions to perform the method described in the first aspect, or any optional embodiment of the first aspect.
[0021] According to a fourth aspect, embodiments of the present invention provide a computer-readable storage medium storing computer instructions for causing a computer to perform the method described in the first aspect, or any optional embodiment of the first aspect.
[0022] The technical solution provided in this application has the following advantages:
[0023] The technical solution provided in this application first establishes a finite element numerical simulation model of the wind turbine tower. Modal analysis is performed on the finite element numerical simulation model to determine the first-order bending natural frequency and the corresponding structural modal mass of the wind turbine tower. Then, the tuning frequency of the annular tuned liquid column damper is calculated using a preset mass ratio and the first-order bending natural frequency. Next, a first and second strip-shaped liquid column are determined from the vertical liquid column of the annular tuned liquid column damper, and the first horizontal length between the first and second strip-shaped liquid columns is initialized. This generates the liquid in the annular tuned liquid column damper. The invention establishes a first numerical relationship between depth, first horizontal length, second horizontal length, and tuning frequency, and determines the liquid depth based on this relationship. Then, it establishes a second numerical relationship between liquid mass, liquid depth, first horizontal length, second horizontal length, and other structural parameters, and calculates these other structural parameters based on this relationship. Finally, it performs bidirectional fluid-structure interaction numerical calculations on the wind turbine tower-annular tuned liquid column damper using the damper structural parameters formed by the first horizontal length, second horizontal length, and other structural parameters to obtain vibration reduction indices. When the vibration reduction indices meet preset vibration reduction conditions, the damper structural parameters are output. This embodiment of the invention closely integrates self-developed theoretical calculations with numerical simulation, improving design efficiency and ensuring design accuracy. Using numerical simulation for the overall calculation of the wind turbine tower and annular tuned liquid column damper allows for a more accurate evaluation of the vibration reduction performance of the annular tuned liquid column damper, meeting the needs of direct application in engineering practice. Attached Figure Description
[0024] The features and advantages of the invention will be more clearly understood by referring to the accompanying drawings, which are schematic and should not be construed as limiting the invention in any way. In the drawings:
[0025] Figure 1 This diagram illustrates the steps of a method for designing a ring-shaped tuned liquid column damper structure according to one embodiment of the present invention.
[0026] Figure 2 A schematic diagram of the structure of an annular tuned liquid column damper according to one embodiment of the present invention is shown;
[0027] Figure 3 A flowchart illustrating a method for designing a ring-shaped tuned liquid column damper structure according to one embodiment of the present invention is shown.
[0028] Figure 4 This diagram illustrates a structural design device for an annular tuned liquid column damper according to one embodiment of the present invention.
[0029] Figure 5 A schematic diagram of a structure design device for an annular tuned liquid column damper according to one embodiment of the present invention is shown. Detailed Implementation
[0030] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not all of them. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.
[0031] Please see Figure 1 In one embodiment, a method for designing a ring-shaped tuned liquid column damper structure specifically includes the following steps:
[0032] Step S101: Establish a finite element numerical simulation model of the wind turbine tower. By performing modal analysis on the finite element numerical simulation model, determine the first-order bending natural frequency of the wind turbine tower and the corresponding structural modal mass.
[0033] Specifically, in this embodiment of the invention, based on the wind turbine tower design drawings, a three-dimensional geometric model of the tower is created using CATIA software. The model is then imported into Hypermesh software, where shell elements are used to mesh the tower. The resulting mesh is then imported into the Modal analysis submodule of ANSYS Workbenck, where material properties (elastic modulus, Poisson's ratio, density) are assigned to the mesh elements. Finally, the first-order bending natural frequency f of the tower can be calculated using the Modal analysis submodule. s and the corresponding modal mass m s .
[0034] Step S102: Calculate the tuning frequency of the annular tuned liquid column damper using the preset mass ratio and the first-order bending natural frequency. The preset mass ratio is the ratio of the liquid mass to the structural modal mass in the annular tuned liquid column damper.
[0035] Specifically, the tuning frequency f of a commonly used toroidal tuned liquid column damper can be used. f and the first-order bending natural frequency f s Calculating the relationship between the tuning frequency f f The formula is as follows:
[0036]
[0037] In the formula, f s denoted as the first-order bending natural frequency of the tower, and μ as the preset mass ratio, representing the ratio of the liquid mass to the structural modal mass in the annular tuned liquid column damper.
[0038] It should be noted that in this embodiment of the invention, the ratio of the liquid mass of the damper to the modal mass of the structure is set to be in the range of 0.01 to 0.05. Thus, when calculating the tuning frequency, the influence of the mass of the annular tuned liquid column damper on the natural frequency of the overall structure is taken into account, so that the structural parameters designed later perform better in terms of vibration reduction performance.
[0039] Step S103: Determine the first strip liquid column and the second strip liquid column from the vertical section liquid column of the annular tuned liquid column damper. The first strip liquid column is the liquid column of the annular tuned liquid column damper that is opposite to the outer shell diameter along the excitation direction. The second strip liquid column is the opposite liquid column of the annular tuned liquid column damper that is adjacent to the first strip liquid column and on the same side.
[0040] Specifically, such as Figure 2 The diagram shows a schematic of a ring-shaped tuned liquid column damper. Since current technology indicates that the ring-shaped tuned liquid column damper exhibits optimal multi-directional vibration control performance when it has eight rectangular baffles, this embodiment of the invention uses a ring-shaped tuned liquid column damper with eight rectangular baffles as an example to describe the steps for generating structural parameters. In practical applications, the structural parameter generation method provided in this embodiment is applicable as long as the number of baffles is an even number greater than or equal to six. This embodiment is not limited to a ring-shaped tuned liquid column damper with eight rectangular baffles.
[0041] Based on the relevant parameters of steps S101 to S102 above, this embodiment proposes a structural parameter calculation formula specifically applicable to annular tuned liquid column dampers through extensive research and simulation. Firstly, in the preparation stage, such as... Figure 2 As shown, in this embodiment, the liquid columns that are opposite each other in the diameter of the outer shell along the excitation direction of the annular tuned liquid column damper are defined as the first strip liquid column, and the opposite liquid columns that are adjacent to the first strip liquid column and on the same side in the annular tuned liquid column damper are defined as the second strip liquid column.
[0042] Step S104: Initialize the first horizontal length between the first strip liquid columns and the second horizontal length between the second strip liquid columns.
[0043] Step S105: Generate a first numerical relationship between the liquid depth, first horizontal length, second horizontal length and tuning frequency in the annular tuned liquid column damper, and determine the liquid depth based on the first numerical relationship.
[0044] Specifically, the first horizontal length H1 between the first strip-shaped liquid columns and the second horizontal length H2 between the second strip-shaped liquid columns are two damper structural parameters that need to be pre-initialized. After the first horizontal length H1 and the second horizontal length H2 are initialized, this embodiment of the invention proposes a separate structural parameter design formula applicable to the annular tuned liquid column damper based on the inherent characteristics of the natural frequency of the annular tuned liquid column damper. That is, the first numerical relationship between the liquid depth, the first horizontal length, the second horizontal length, and the tuning frequency in the annular tuned liquid column damper is proposed, thereby realizing the structural design.
[0045] Specifically, this embodiment lists the kinetic energy, potential energy, and non-conservative dissipative forces of the liquid in each part of the annular tuned liquid column damper under external excitation. Substituting these energies and forces into the Lagrange equations of motion yields the equations of motion for the annular tuned liquid column damper. These equations of motion then determine the natural frequency characteristics of the annular tuned liquid column damper. Based on these natural frequency characteristics, this embodiment proposes the following first numerical relationship between the first horizontal length, the second horizontal length, and the tuning frequency:
[0046]
[0047] In the formula, χ represents the frequency correction coefficient, with a default value of 1, and η represents the cross-sectional area ratio, with a default value of 1. H1 represents the first horizontal length, H2 represents the second horizontal length, and f f Indicates the tuning frequency.
[0048] Specifically, assuming the frequency correction coefficient χ = 1 and the cross-sectional area ratio η = 1 (the cross-sectional area ratio is the ratio of the cross-sectional area of any strip-shaped liquid column to the cross-sectional area of its own horizontal segment liquid column), the first horizontal length H1, the second horizontal length H2, and the tuning frequency f calculated in step one are used. f Substituting into the above formula, we can output the liquid height V that should be filled into the annular tuned liquid column damper.
[0049] Step S106: Establish a second numerical relationship between liquid mass, liquid depth, first horizontal length, second horizontal length and other structural parameters, and determine other structural parameters based on the second numerical relationship. Other structural parameters include the outer shell radius, inner shell radius, horizontal channel height, inner shell height and outer shell height of the annular tuned liquid column damper.
[0050] Specifically, in steps S101 to S102, the preset mass ratio value μ, representing the ratio of liquid mass to structural modal mass, has been initialized and set, and the structural modal mass m... s The liquid mass m has been obtained through modeling in step S101. f Structural modal mass ms It is obtained by multiplying by the preset mass ratio μ. Since the liquid mass m f Certainly, the liquid depth required to reach the tuning frequency of the annular tuned liquid column damper is also obtained in step S105. To allow a fixed liquid mass to be injected into the damper to reach a fixed liquid depth, a second numerical relationship between the liquid depth and the liquid mass is established. In addition to the first horizontal length H1 and the second horizontal length H2, other structural parameters that meet the above conditions can also be determined according to certain calculation and analysis methods, thereby forming a complete structural parameter scheme.
[0051] Step S107: Using the first horizontal length, the second horizontal length, and other structural parameters as damper structural parameters, perform bidirectional fluid-structure interaction numerical calculations on the wind turbine tower-annular tuned liquid column damper to obtain vibration reduction indexes.
[0052] Specifically, this embodiment of the invention utilizes the FLUENT and Transientstructural submodules in the ANSYS Workbench application to perform bidirectional fluid-structure interaction numerical calculations on a wind turbine tower-annular tuned liquid column damper, thereby verifying the performance of the annular tuned liquid column damper. After determining the liquid depth in the annular tuned liquid column damper as V0′, the dynamic response is calculated when the liquid depth is V=V0′ and V=0, respectively. Using the standard deviation of the displacement response of the top of the wind turbine tower (i.e., the nacelle) in the downwind and vertical wind directions as indicators, the attenuation rate is statistically analyzed. The vibration reduction performance of the annular tuned liquid column damper is evaluated using the attenuation rate; the larger the attenuation rate, the better the vibration reduction performance. The specific steps include: 1. Using CATIA software, determine the location where the annular tuned liquid column damper needs to be installed in the wind turbine tower established in step S101; 2. Using the damper structural parameters of the annular tuned liquid column damper determined in steps S102 to S106, establish a three-dimensional geometric model of the annular tuned liquid column damper and import it into Hypermesh for mesh generation; 3. Import the mesh models of the wind turbine tower and the annular tuned liquid column damper into the Transient structural submodule of ANSYS Workbench, and set the annular tuned liquid column damper as the interface for fluid-structure interaction calculation; 4. After calculating the wind load (and wave load) of the wind turbine under the conditions of cut-in wind speed, rated operating wind speed, cut-out wind speed, and extreme wind speed, apply it as an external load to the Transient structural submodule of ANSYS Workbench. 5. In the structural submodule, the corresponding mesh element of the wind turbine tower is used, where the wind load is applied to the top of the tower and the wave load is applied below the water level; 6. The wall in the computational fluid dynamics model of the annular tuned liquid column damper is set as the fluid-structure interaction interface; 7. Two-way fluid-structure interaction numerical calculations are carried out for the wind turbine tower-annular tuned liquid column damper. The dynamic response is calculated when the liquid depth is V=V0′ and V=0 respectively. The attenuation rate is calculated using the standard deviation of the displacement response of the top of the wind turbine tower (i.e., the nacelle) in the downwind and vertical directions as indicators to evaluate the vibration reduction performance of the annular tuned liquid column damper.
[0053] Step S108: Determine whether the vibration reduction index meets the preset vibration reduction conditions. If it does, output the damper structural parameters; otherwise, return to step S103 to redetermine the damper structural parameters.
[0054] Specifically, the vibration reduction conditions of the preset vibration reduction index in this embodiment of the invention characterize whether the structural parameters of the annular tuned liquid column damper can meet the vibration reduction requirements. For example, it determines whether the attenuation rate output in step S107 exceeds a preset value. If the attenuation rate exceeds the preset value, it indicates that the damper structural parameters designed in the above steps are effective, and thus the damper structural parameters are output. If the attenuation rate does not exceed the preset value, step S103 is re-executed to determine new damper structural parameters, and the vibration reduction index is recalculated through step S107 until the damper structural parameters can meet the vibration reduction requirements.
[0055] The solution provided by this invention closely integrates theoretical calculations and numerical simulations, improving design efficiency while ensuring design accuracy. Using numerical simulation for the overall calculation of the wind turbine tower and the annular tuned liquid column damper allows for a more accurate evaluation of the vibration reduction performance of the annular tuned liquid column damper, meeting the needs for direct application in engineering practice.
[0056] Specifically, in one embodiment, before step S106 described above, the following step is also included:
[0057] Step 1: Calculate the initial length ratio p using the following formula.
[0058]
[0059] In the formula, H1 represents the first horizontal length, and V represents the liquid depth.
[0060] Step 2: Determine whether the initial length ratio falls within the range of [0.5, 0.6].
[0061] Step 3: If it does not fall within the range of [0.5, 0.6], then reinitialize the first and second horizontal lengths.
[0062] Specifically, based on the research results of this application in the application scenario, through formula
[0063]
[0064] The calculated initial length ratio p characterizes the control performance of the annular tuned liquid column damper. Typically, the initial length ratio p ranges from 0.2 to 0.8; the larger p is, the better the control performance and the larger the liquid response amplitude. In this embodiment of the invention, since the tower imposes certain limitations on the outer shell radius of the annular tuned liquid column damper, a value of p between 0.5 and 0.6 is optimal. Therefore, before step S106, this embodiment of the invention requires calculating the initial length ratio using the first horizontal length H1 and determining whether the initial length ratio falls within the range of [0.5, 0.6]. If the initial length ratio does not fall within the range of [0.5, 0.6], the first and second horizontal lengths are re-initialized. When the initial length ratio meets the condition of [0.5, 0.6], step S106 is executed, thereby enabling the generated structural parameters to further improve the control performance of the annular tuned liquid column damper. If the outer shell radius exceeds the limitation of the inner radius of the wind turbine tower, the value of p should be reduced, and the structural parameters of the annular tuned liquid column damper should be recalculated.
[0065] Specifically, in one embodiment, step S106 above includes the following steps:
[0066] Step 4: Establish the following second numerical relationship, and calculate the outer shell radius, inner shell radius, and horizontal channel height of the annular tuned liquid column damper using the following formula.
[0067] m f =8A V V+A H H1+2A H H2,
[0068] In the formula, m f A represents the mass of the liquid. V Let A represent the cross-sectional area of a certain strip-shaped liquid column. H R1 and R2 represent the cross-sectional area of the corresponding equivalent TLCD horizontal segment liquid column, respectively, and B3 represents the outer shell radius and inner shell radius of the annular tuned liquid column damper.
[0069] Specifically, for the annular tuned liquid column damper, based on extensive theoretical analysis and simulation studies, this embodiment of the invention proposes that the second numerical relationship satisfied by the liquid mass and liquid depth is: m f =8A V V+A H H1+2A H H2. In this embodiment, the cross-sectional area ratio η = 1 is set. When the cross-sectional area ratio η = 1, the control performance of the annular tuned liquid column damper is determined to be optimal. Therefore, the cross-sectional area A of a certain strip-shaped liquid column is calculated using the above formula. VThe cross-sectional area A of a certain equivalent TLCD horizontal segment liquid column H Then, substituting the liquid mass, liquid depth, first horizontal length, and second horizontal length into the above formula, the outer shell radius R1, inner shell radius R2, and horizontal channel height B3 of the annular tuned liquid column damper are calculated.
[0070] Step 5: Determine the housing height of the annular tuned liquid column damper by using a preset multiple of the liquid depth.
[0071] Specifically, the present invention determines the shell height B1 of the annular tuned liquid column damper by a preset multiple of the liquid depth V. The research of this application shows that the optimal value range of the preset multiple is 2 to 3 times.
[0072] Step 6: Calculate the inner shell height of the annular tuned liquid column damper based on the difference between the outer shell height and the horizontal channel height.
[0073] Specifically, the inner shell height of the annular tuned liquid column damper is further calculated as B2 = B1 - B3. The thickness of the annular tuned liquid column damper is R3 = R1 - R2. The outer shell radius, inner shell radius, horizontal channel height, inner shell height, and outer shell height are then output as damper structural parameters. After vibration reduction performance verification, these parameters can be used as the optimal design scheme for the damper, improving the vibration reduction effect of the annular tuned liquid column damper.
[0074] Specifically, in one embodiment, such as Figure 3 As shown in the embodiment of the present invention, the method for determining the structural parameters of an annular tuned liquid column damper further includes the following steps:
[0075] Step 7: Create a fluid dynamics model using the damper structural parameters.
[0076] Step 8: Add the liquid depth to the preset compensation to obtain the compensated liquid depth.
[0077] Step 9: Input the compensated liquid depth into the fluid dynamics model to simulate free decay vibration and find the simulated liquid depth corresponding to the tuning frequency.
[0078] Step 10: Update the liquid depth value using simulated liquid depth.
[0079] Step 11: Correct the damper structural parameters by updating the relationship between the fluid depth and fluid mass.
[0080] Specifically, in step S105, the assumed frequency correction coefficient χ = 1 means that the vertical section of the annular tuned liquid column damper is assumed to participate entirely in vertical oscillation without any horizontal swaying. This is an ideal state, which is usually unattainable in real-world environments. This results in the calculated liquid depth in step S105 being too small, thus requiring further improvement in the accuracy of the damper's structural parameters. In this embodiment of the invention, a free decay vibration simulation using a fluid dynamics model is performed to correct the liquid depth of the damper, thereby further improving the accuracy of subsequent structural parameter schemes.
[0081] First, a three-dimensional computational fluid dynamics model is established using the damper structural parameters output in step S106. The liquid depth in the model is set to be slightly larger than the liquid depth calculated in step S105. Specifically, the liquid depth is added to a preset compensation to obtain a slightly larger compensation liquid depth. In this embodiment, the preset compensation is set to be the optimal value of the horizontal channel height B3. Then, the compensation liquid depth is input into the fluid dynamics model for free decay vibration simulation. The dynamic response of the liquid is extracted and subjected to Fourier spectrum analysis. Through multiple trials, the design value f of the liquid frequency is found. f The corresponding simulated liquid depth V0′. The specific steps for attenuated vibration simulation include: obtaining the damper structural parameters output in step S105; using CATIA software to create a three-dimensional geometric model of the annular tuned liquid column damper; importing the model into Hypermesh; using a three-dimensional hexahedral mesh to divide the annular tuned liquid column damper into mesh elements; subsequently importing the mesh elements into the FLUENT analysis submodule of ANSYS Workbenck for parameter setting; and setting different initial liquid depths V0′. i Self-degrading vibration simulation was conducted, and the time history of the self-degrading vibration was Fourier transformed to determine the initial liquid depth V. i The corresponding frequency f i Finally, the design value f was determined. f The corresponding simulated liquid depth V0′ is used to correct the liquid depth.
[0082] It is important to note that when using a three-dimensional hexahedral mesh to generate the mesh elements for the annular tuned liquid column damper in the damping vibration simulation, it is necessary to simultaneously select hexahedral meshes of different sizes for the same mesh generation method. For different mesh sizes, the same initial liquid depth should be set in the FLUENT submodule, and free damping vibration calculations should be performed separately. Only by comparing the calculation results can it be ensured that the mesh generation does not affect the calculation results.
[0083] Furthermore, in this embodiment, when setting the parameters of the FLUENT submodule, a pressure-based solver is selected for transient solution; a multiphase flow model of Volume of Fluid is selected, a Realizable k-ε turbulence model is adopted, and a standard wall function is used; material parameters for the liquid and ambient air are set; the container surface of the annular tuned liquid column damper is set as a wall, and the surfaces above the cylindrical outer shell and inner shell of the annular tuned liquid column damper that are connected to the atmosphere are set as pressure outlet surfaces; by writing a user-defined function (UDF), the motion conditions of the fluid domain inside the annular tuned liquid column damper are given, preferably a sinusoidal motion signal with a frequency equal to the design natural frequency f of the annular tuned liquid column damper. f With an amplitude of 0.1 times R1 and a motion period of 10 cycles; dividing the gas-liquid two phases: taking a certain V i The values are used to divide the gas phase and liquid phase, with the upper part being the gas phase and the lower part being the liquid phase; a monitoring point for the liquid level height of the annular tuned liquid column damper is created, and the output results are saved during the calculation process; the fixed time step for the calculation is preferably 0.01s, and the calculation time length can be selected as 20 cycles of a sinusoidal motion signal.
[0084] Finally, this embodiment re-executes steps four to six based on the corrected liquid depth to update the values of the damper structural parameters, thereby outputting more accurate damper structural parameters.
[0085] Specifically, in one embodiment, the method for determining the structural parameters of an annular tuned liquid column damper provided by the present invention further includes the following steps:
[0086] Step 12: Substitute the updated liquid depth into the formula In the formula, a new frequency correction coefficient is output; where χ' represents the new frequency correction coefficient, η represents the cross-sectional area ratio, and η is preset to 1. H1 represents the first horizontal length, H2 represents the second horizontal length, and f f Indicates the tuning frequency.
[0087] Specifically, in this embodiment of the invention, based on the corrected liquid depth output in step thirteen above, and combined with the aforementioned cross-sectional area ratio, first horizontal length, second horizontal length, and tuning frequency series parameters, the values are substituted back into the following formula.
[0088]
[0089] This yields a more accurate new frequency correction coefficient χ'. Based on this new frequency correction coefficient, it can be directly used in the design of the next annular tuned liquid column damper structural parameters. Compared to the calculation process assuming a frequency correction coefficient of 1, this further improves the generation efficiency and accuracy of the next annular tuned liquid column damper structural parameters.
[0090] Specifically, in one embodiment, the present invention performs the above steps S101 to step 12 for various extreme operating conditions, establishes the design space between the tuning frequency, structural parameters and vibration reduction effect of the annular tuned liquid column damper, and saves the new frequency correction coefficient χ' for each extreme operating condition. Thus, in the subsequent optimization process, the frequency correction coefficient for various calculated operating conditions can be directly read and used, thereby improving the parameter design and optimization efficiency of the annular tuned liquid column damper for wind turbine towers of different unit models and tower sizes in the same wind farm.
[0091] Through the above steps, the technical solution provided in this application first establishes a finite element numerical simulation model of the wind turbine tower. Modal analysis is then performed on the finite element numerical simulation model to determine the first-order bending natural frequency and the corresponding structural modal mass of the wind turbine tower. Next, the tuning frequency of the annular tuned liquid column damper is calculated using a preset mass ratio and the first-order bending natural frequency. Then, the first and second strip liquid columns are determined from the vertical section of the annular tuned liquid column damper, and the first horizontal length between the first and second strip liquid columns is initialized. This generates the tuning frequency of the annular tuned liquid column damper. The invention establishes a first numerical relationship between the liquid depth, first horizontal length, second horizontal length, and tuning frequency, and determines the liquid depth based on this relationship. Then, it establishes a second numerical relationship between the liquid mass, liquid depth, first horizontal length, second horizontal length, and other structural parameters, and calculates these other structural parameters based on this relationship. Finally, it performs bidirectional fluid-structure interaction numerical calculations on the wind turbine tower-annular tuned liquid column damper using the damper structural parameters formed by the first horizontal length, second horizontal length, and other structural parameters to obtain vibration reduction indicators. When the vibration reduction indicators meet preset vibration reduction conditions, the damper structural parameters are output. The solution provided in this embodiment of the invention closely integrates self-developed theoretical calculations with numerical simulation, improving design efficiency and ensuring design accuracy. Using numerical simulation for the overall calculation of the wind turbine tower and annular tuned liquid column damper allows for a more accurate evaluation of the vibration reduction performance of the annular tuned liquid column damper, meeting the needs of direct application in engineering practice.
[0092] like Figure 4 As shown, this embodiment also provides a design device for a ring-shaped tuned liquid column damper structure, the device comprising:
[0093] The tower simulation module 101 is used to establish a finite element numerical simulation model of the wind turbine tower. By performing modal analysis on the finite element numerical simulation model, the first-order bending natural frequency and the corresponding structural modal mass of the wind turbine tower are determined. For details, please refer to the relevant description of step S101 in the above method embodiment, which will not be repeated here.
[0094] The tuning frequency calculation module 102 is used to calculate the tuning frequency of the annular tuned liquid column damper using a preset mass ratio and the first-order bending natural frequency. The preset mass ratio is the ratio of the liquid mass to the structural modal mass in the annular tuned liquid column damper. For details, please refer to the relevant description of step S102 in the above method embodiments, which will not be repeated here.
[0095] The equivalent liquid column determination module 103 is used to determine a first strip liquid column and a second strip liquid column from the vertical section liquid column of the annular tuned liquid column damper. The first strip liquid column is a liquid column opposite to the first strip liquid column in the diameter direction of the outer shell of the annular tuned liquid column damper, and the second strip liquid column is a liquid column in the annular tuned liquid column damper that is adjacent to the first strip liquid column and on the same side. For details, please refer to the relevant description of step S103 in the above method embodiment, which will not be repeated here.
[0096] The horizontal length initialization module 104 is used to initialize the first horizontal length between the first strip liquid columns and the second horizontal length between the second strip liquid columns. For details, please refer to the relevant description of step S104 in the above method embodiment, which will not be repeated here.
[0097] The liquid depth calculation module 105 is used to establish a first numerical relationship between the liquid depth, the first horizontal length, the second horizontal length, and the tuning frequency in the annular tuned liquid column damper, and to determine the liquid depth based on the first numerical relationship. For details, please refer to the relevant description of step S105 in the above method embodiment, which will not be repeated here.
[0098] The other structural parameter calculation module 106 is used to establish a second numerical relationship between liquid mass, liquid depth, first horizontal length, second horizontal length, and other structural parameters, and to determine other structural parameters based on the second numerical relationship. These other structural parameters include the outer shell radius, inner shell radius, horizontal channel height, inner shell height, and outer shell height of the annular tuned liquid column damper. For details, please refer to the relevant description of step S106 in the above method embodiment; it will not be repeated here.
[0099] The vibration reduction experimental module 107 is used to perform bidirectional fluid-structure interaction numerical calculations on the wind turbine tower-annular tuned liquid column damper using the first horizontal length, the second horizontal length, and other structural parameters as damper structural parameters, to obtain vibration reduction indicators. For details, please refer to the relevant description of step S107 in the above method embodiment, which will not be repeated here.
[0100] The structural parameter output module 108 determines whether the vibration reduction index meets the preset vibration reduction conditions. If it does, it outputs the damper structural parameters. For details, please refer to the relevant description of step S108 in the above method embodiment, which will not be repeated here.
[0101] The present invention provides a device for designing a ring-shaped tuned liquid column damper structure, which is used to execute a method for designing a ring-shaped tuned liquid column damper structure provided in the above embodiments. Its implementation method and principle are the same. For details, please refer to the relevant description of the above method embodiments, which will not be repeated here.
[0102] Through the collaborative efforts of the aforementioned components.
[0103] The technical solution provided in this application first establishes a finite element numerical simulation model of the wind turbine tower. Modal analysis is performed on the finite element numerical simulation model to determine the first-order bending natural frequency and the corresponding structural modal mass of the wind turbine tower. Then, the tuning frequency of the annular tuned liquid column damper is calculated using a preset mass ratio and the first-order bending natural frequency. Next, a first and second strip-shaped liquid column are determined from the vertical liquid column of the annular tuned liquid column damper, and the first horizontal length between the first and second strip-shaped liquid columns is initialized. This establishes the liquid structure within the annular tuned liquid column damper. The invention establishes a first numerical relationship between depth, first horizontal length, second horizontal length, and tuning frequency, and determines the liquid depth based on this relationship. Then, it establishes a second numerical relationship between liquid mass, liquid depth, first horizontal length, second horizontal length, and other structural parameters, and calculates these other structural parameters based on this relationship. Finally, it performs bidirectional fluid-structure interaction numerical calculations on the wind turbine tower-annular tuned liquid column damper using the damper structural parameters formed by the first horizontal length, second horizontal length, and other structural parameters to obtain vibration reduction indices. When the vibration reduction indices meet preset vibration reduction conditions, the damper structural parameters are output. This embodiment of the invention closely integrates self-developed theoretical calculations with numerical simulation, improving design efficiency and ensuring design accuracy. Using numerical simulation for the overall calculation of the wind turbine tower and annular tuned liquid column damper allows for a more accurate evaluation of the vibration reduction performance of the annular tuned liquid column damper, meeting the needs of direct application in engineering practice.
[0104] Figure 5 This invention illustrates a design device for a ring-shaped tuned liquid column damper structure according to an embodiment of the present invention. The device includes a processor 901 and a memory 902, which can be connected via a bus or other means. Figure 5 Taking the example of a connection between China and Israel via a bus.
[0105] Processor 901 can be a Central Processing Unit (CPU). Processor 901 can also be other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or combinations of the above types of chips.
[0106] The memory 902, as a non-transitory computer-readable storage medium, can be used to store non-transitory software programs, non-transitory computer-executable programs, and modules, such as the program instructions / modules corresponding to the methods in the above method embodiments. The processor 901 executes various functional applications and data processing of the processor by running the non-transitory software programs, instructions, and modules stored in the memory 902, thereby implementing the methods in the above method embodiments.
[0107] The memory 902 may include a program storage area and a data storage area. The program storage area may store the operating system and applications required for at least one function; the data storage area may store data created by the processor 901, etc. Furthermore, the memory 902 may include high-speed random access memory and may also include non-transitory memory, such as at least one disk storage device, flash memory device, or other non-transitory solid-state storage device. In some embodiments, the memory 902 may optionally include memory remotely located relative to the processor 901, and these remote memories may be connected to the processor 901 via a network. Examples of such networks include, but are not limited to, the Internet, corporate intranets, local area networks, mobile communication networks, and combinations thereof.
[0108] One or more modules are stored in memory 902, and when executed by processor 901, they perform the methods described in the above method embodiments.
[0109] The specific details of the equipment for determining the structural parameters of the annular tuned liquid column damper can be understood by referring to the relevant descriptions and effects in the above method embodiments, and will not be repeated here.
[0110] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The implemented program can be stored in a computer-readable storage medium. When the program is executed, it can include the processes of the embodiments of the above methods. The storage medium can be a magnetic disk, optical disk, read-only memory (ROM), random access memory (RAM), flash memory, hard disk drive (HDD), or solid-state drive (SSD), etc.; the storage medium can also include combinations of the above types of memory.
[0111] Although embodiments of the invention have been described in conjunction with the accompanying drawings, those skilled in the art can make various modifications and variations without departing from the spirit and scope of the invention, and such modifications and variations all fall within the scope defined by the appended claims.
Claims
1. A method for designing a ring-shaped tuned liquid column damper structure, characterized in that, The method includes: A finite element numerical simulation model of the wind turbine tower is established. Modal analysis is performed on the finite element numerical simulation model to determine the first-order bending natural frequency and the corresponding structural modal mass of the wind turbine tower. The tuning frequency of the annular tuned liquid column damper is calculated by using a preset mass ratio and the first-order bending natural frequency. The preset mass ratio is the ratio of the liquid mass in the annular tuned liquid column damper to the structural modal mass. A first strip liquid column and a second strip liquid column are determined from the vertical section liquid column of the annular tuned liquid column damper. The first strip liquid column is the liquid column of the annular tuned liquid column damper that is opposite to the outer shell diameter along the excitation direction. The second strip liquid column is the opposite liquid column of the annular tuned liquid column damper that is adjacent to the first strip liquid column and on the same side. Initialize the first horizontal length between the first and second liquid columns; A first numerical relationship is established between the liquid depth, the first horizontal length, the second horizontal length, and the tuning frequency in the annular tuned liquid column damper, and the liquid depth is determined based on the first numerical relationship. A second numerical relationship is established between the liquid mass, the liquid depth, the first horizontal length, the second horizontal length, and other structural parameters, and the other structural parameters are determined based on the second numerical relationship. The other structural parameters include the outer shell radius, inner shell radius, horizontal channel height, inner shell height, and outer shell height of the annular tuned liquid column damper. Using the first horizontal length, the second horizontal length, and the other structural parameters as damper structural parameters, a two-way fluid-structure interaction numerical calculation of the wind turbine tower-annular tuned liquid column damper is carried out to obtain the vibration reduction index. Determine whether the vibration reduction index meets the preset vibration reduction conditions. If it does, output the damper structural parameters.
2. The method according to claim 1, characterized in that, The first numerical relationship between the liquid depth, the first horizontal length, the second horizontal length, and the tuning frequency in the generated annular tuned liquid column damper is as follows: In the formula, χ represents the frequency correction coefficient, with a default value of 1, and η represents the cross-sectional area ratio, with a default value of 1. H1 represents the first horizontal length, H2 represents the second horizontal length, and f f The tuning frequency is represented by , and g represents the gravitational acceleration.
3. The method according to claim 2, characterized in that, Before determining the other structural parameters based on the second numerical relationship, the method further includes: The initial length ratio p is calculated using the following formula. In the formula, H1 represents the first horizontal length, and V represents the liquid depth; Determine whether the initial length ratio falls within the range of [0.5, 0.6]. If it does not fall within the range of [0.5, 0.6], then the first horizontal length and the second horizontal length are reinitialized.
4. The method according to claim 2, characterized in that, The second numerical relationship established includes: m f =8A V V+A H H1+2A H H2, In the formula, m f A represents the mass of the liquid. V Let A represent the cross-sectional area of a certain strip-shaped liquid column. H R1 and R2 represent the cross-sectional area of the corresponding equivalent TLCD horizontal segment liquid column, respectively, B3 represents the outer shell radius and inner shell radius of the annular tuned liquid column damper, and V represents the liquid depth. Determining the other structural parameters based on the second numerical relationship includes: Substitute the liquid mass, the liquid depth, the first horizontal length, and the second horizontal length into the second numerical relationship to calculate the outer shell radius, inner shell radius, and horizontal channel height of the annular tuned liquid column damper; The housing height of the annular tuned liquid column damper is determined by a preset multiple of the liquid depth; The inner shell height of the annular tuned liquid column damper is calculated based on the difference between the outer shell height and the horizontal channel height.
5. The method according to claim 4, characterized in that, The method further includes: A fluid dynamics model is created using the damper structural parameters; The liquid depth is added to the preset compensation to obtain the compensated liquid depth; The compensated liquid depth is input into the fluid dynamics model to simulate free decay vibration, and the simulated liquid depth corresponding to the tuning frequency is found. The simulated liquid depth is used to update the numerical value of the liquid depth; The damper structural parameters are corrected by updating the relationship between the liquid depth and the liquid mass.
6. The method according to claim 5, characterized in that, The method further includes: Substitute the updated liquid depth into the formula In the formula, a new frequency correction coefficient is output; where χ' represents the new frequency correction coefficient, η represents the cross-sectional area ratio, and η is preset to 1. H1 represents the first horizontal length, H2 represents the second horizontal length, and f f The tuning frequency is represented by , and g represents the gravitational acceleration.
7. The method according to claim 5, characterized in that, The preset compensation is equal to the horizontal channel height.
8. A structural design device for a ring-shaped tuned liquid column damper, characterized in that, The device includes: The tower simulation module is used to establish a finite element numerical simulation model of the wind turbine tower. By performing modal analysis on the finite element numerical simulation model, the first-order bending natural frequency and the corresponding structural modal mass of the wind turbine tower are determined. The tuning frequency calculation module is used to calculate the tuning frequency of the annular tuned liquid column damper by using a preset mass ratio and the first-order bending natural frequency. The preset mass ratio is the ratio of the liquid mass in the annular tuned liquid column damper to the structural modal mass. The equivalent liquid column determination module is used to determine a first strip liquid column and a second strip liquid column from the vertical section liquid column of the annular tuned liquid column damper. The first strip liquid column is the liquid column of the annular tuned liquid column damper that is opposite to the outer shell diameter along the excitation direction. The second strip liquid column is the opposite liquid column of the annular tuned liquid column damper that is adjacent to the first strip liquid column and on the same side. A horizontal length initialization module is used to initialize the first horizontal length between the first strip liquid columns and the second horizontal length between the second strip liquid columns; The liquid depth calculation module is used to establish a first numerical relationship between the liquid depth, the first horizontal length, the second horizontal length and the tuning frequency in the annular tuned liquid column damper, and to determine the liquid depth based on the first numerical relationship. The other structural parameter calculation module is used to establish a second numerical relationship between the liquid mass, the liquid depth, the first horizontal length, the second horizontal length and other structural parameters, and to determine the other structural parameters based on the second numerical relationship. The other structural parameters include the outer shell radius, inner shell radius, horizontal channel height, inner shell height and outer shell height of the annular tuned liquid column damper. The vibration reduction test module is used to perform bidirectional fluid-structure interaction numerical calculations on the wind turbine tower-annular tuned liquid column damper using the first horizontal length, the second horizontal length, and the other structural parameters as damper structural parameters, and to obtain vibration reduction indicators. The structural parameter output module determines whether the vibration reduction index meets the preset vibration reduction conditions. If it does, it outputs the damper structural parameters.
9. A device for designing a ring-shaped tuned liquid column damper structure, characterized in that, include: A memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, the processor executing the computer instructions to perform the method as described in any one of claims 1-7.
10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions for causing the computer to perform the method as described in any one of claims 1-7.
Citation Information
Patent Citations
Tuned air pressure liquid column damper and tower drum
CN114909010A
Tuning liquid column damper and tower drum
CN217501237U
Multilayer circular ring-shaped tuned liquid column damper
CN112761401A
Tuning liquid column damper with adjustable liquid level and tower drum
CN114809765A