Parametric modeling method and system for hybrid tower cable system based on multi-body dynamics

Abstract

The present invention discloses a parameterized modeling method and system of a mixed tower steel cable system based on multi-body dynamics, and relates to the technical field of wind power generation. The method comprises the steps of: creating parameter points at the center of each end face of a concrete segment according to geometric parameters, and establishing a three-dimensional model of a mixed tower; correcting the mass of each segment of the mixed tower according to the material density of the mixed tower; simulating the pressure-resistant and tension-intolerant characteristics of the mixed tower, and establishing contact and bonding relationships between each concrete segment of the three-dimensional model of the mixed tower; obtaining the design parameters of the steel cables in the mixed tower, constructing a parameterized model of the steel cables, and combining the parameterized model of the steel cables with the three-dimensional model of the mixed tower to obtain a parameterized model of the mixed tower steel cable system; performing dynamic simulation on the operating conditions of the wind turbine, and obtaining the time domain response results of the parameterized model of the mixed tower steel cable system. The present invention realizes the response characteristic analysis of the mixed tower under multiple working conditions and multiple structures.

Description

Technical Field

[0001] The present invention relates to the technical field of wind power generation, and in particular to a parameterized modeling method and system for a hybrid tower cable system based on multi-body dynamics. Background Art

[0002] The statements in this section merely provide background information related to the present invention and do not necessarily constitute prior art.

[0003] In recent years, as the power generation capacity and rotor diameter of wind turbines have increased year by year, the tower height has increased year by year. Concrete-steel hybrid towers (referred to as hybrid towers) are increasingly favored by the market due to their large system stiffness, convenient prefabrication and installation, and certain economic efficiency. Due to the short time of hybrid towers, the low proportion of in-service, the diversity of hybrid tower structures, and the differences in wind resources in wind farms, there is little operating data on existing hybrid tower wind turbines, especially a serious lack of response characteristic data in the time and frequency domains of hybrid tower wind turbines under extreme working conditions, making it difficult to effectively draw on the response characteristic data of existing hybrid towers in the time and frequency domains to guide the design of hybrid tower cable systems for other wind farms. With the development of computer-aided engineering technology, simulation models can effectively predict the characteristics of hybrid tower cable systems. Although there are some simulation calculations of hybrid tower steel cable structures at present, these simulation studies focus on the static strength and modal vibration of the hybrid tower structure, and lack the study of the dynamic characteristics of the hybrid tower steel cable system; focus on analyzing the characteristics of a single given hybrid tower structure, and lack the study of parametric modeling methods for hybrid tower steel cable systems; focus on analyzing the response of the structure at the load moment of a given working condition (maximum working condition), and lack the time domain dynamic characteristics of the hybrid tower steel cable system in continuous time under working conditions such as fan startup and yaw.

[0004] In summary, the existing simulation model of the hybrid tower cable structure cannot fully reflect the time domain response characteristics under various complex working conditions, and cannot show the research on parametric modeling methods under various structural characteristics. Summary of the invention

[0005] In view of the deficiencies in the prior art, the purpose of the present invention is to provide a parametric modeling method and system for a hybrid tower cable system based on multi-body dynamics, which realizes the response characteristic analysis of the hybrid tower under multiple working conditions and multiple structures.

[0006] In order to achieve the above object, the present invention is implemented through the following technical solutions:

[0007] A first aspect of the present invention provides a parameterized modeling method for a hybrid tower cable system based on multi-body dynamics, comprising the following steps:

[0008] Obtain the geometric parameters of the concrete tower, create parameter points at the center of each end surface of the concrete segment according to the geometric parameters, and establish a three-dimensional model of the concrete tower;

[0009] The mass of each section of the mixing tower is corrected according to the material density of the mixing tower until the deviation between the mass of the mixing tower and the actual mass is less than the preset range;

[0010] Simulate the characteristics of the concrete tower that it is resistant to pressure but not to tension, and establish the contact and bonding relationship between the concrete segments of the concrete tower 3D model;

[0011] The design parameters of the steel cables in the hybrid tower are obtained, and a parametric model of the steel cables is constructed. The parametric model of the steel cables is combined with the three-dimensional model of the hybrid tower to obtain a parametric model of the hybrid tower steel cable system.

[0012] The operating conditions of the wind turbine are dynamically simulated to obtain the time domain response results of the parameterized model of the hybrid tower cable system.

[0013] Furthermore, the geometric parameters of the mixing tower include the height of the mixing tower, the number of concrete sections, the radius of the mixing tower and the wall thickness of the mixing tower.

[0014] Furthermore, the specific steps of establishing the contact and bonding relationship between the concrete segments of the three-dimensional model of the concrete tower are as follows:

[0015] Setting contact relationships between concrete sections, between concrete sections and wind turbine foundations, and between concrete sections and transition sections;

[0016] Based on the Hertz contact theory, the contact stiffness is determined according to the contact radius, the elastic modulus of the material and the Poisson's ratio, and the concrete segments are arranged together due to the effect of gravity.

[0017] Furthermore, the specific steps to simulate the pressure-resistant and tensile-resistant characteristics of the mixing tower are as follows:

[0018] According to the Young's modulus, strain and cross-sectional area of ​​the bonding material, the bonding force and bonding moment in the six directions of freedom of each bonding surface of the concrete segment are calculated;

[0019] Use spline curves and function expressions to express the changing relationship between bonding force and displacement, bonding moment and rotation angle, and accurately simulate the compressive but not tensile properties of the mixed tower;

[0020] Set the critical position to simulate the fracture failure characteristics of the mixing tower.

[0021] Furthermore, the design parameters of the steel cable include the cable arrangement position, quantity, diameter, axial stiffness, damping, anchor point position and preload.

[0022] Furthermore, the specific steps of dynamic simulation of the operating conditions of the fan are as follows:

[0023] Construct the rotation constraints of the wind rotor hub and the nacelle, and use the rotation drive to simulate the wind turbine startup condition. In this case, the aerodynamic load of the wind turbine rotor is equivalently calculated according to the power, blade aerodynamic characteristics and speed information of the wind turbine to be analyzed, and the aerodynamic load is applied to the wind turbine rotor using spline curves and function expressions to simulate the wind load of the wind turbine from startup to rated operating speed.

[0024] The rotation constraints of the nacelle and tower are constructed, and the rotation drive is used to simulate the torsional effect of the yaw-to-wind motion on the wind turbine tower.

[0025] Furthermore, the specific steps to obtain the time domain response results of the parameterized model of the hybrid tower cable system are as follows:

[0026] Carry out multi-condition simulation of normal, extreme and fault conditions of hybrid tower cable system;

[0027] The time domain response results of the four-degree-of-freedom swing angle, angular velocity, angular acceleration and cable preload of the hybrid tower are obtained.

[0028] The second aspect of the present invention provides a parameterized modeling system for a hybrid tower cable system based on multi-body dynamics, comprising:

[0029] A mixed tower model building module is configured to obtain geometric parameters of the mixed tower, create parameter points at the center of each end surface of the concrete segment according to the geometric parameters, and establish a three-dimensional model of the mixed tower;

[0030] The mass correction module is configured to correct the mass of each section of the mixing tower according to the material density of the mixing tower until the deviation between the mass of the mixing tower and the actual mass is less than a preset range;

[0031] The characteristic simulation module is configured to simulate the compressive and tensile properties of the concrete tower, and to establish contact and bonding relationships between concrete segments of the concrete tower three-dimensional model;

[0032] A cable model building module is configured to obtain design parameters of the cables in the mixed tower, build a cable parametric model, and combine the cable parametric model with the mixed tower three-dimensional model to obtain a mixed tower cable system parametric model;

[0033] The dynamic simulation module is configured to perform dynamic simulation on the operating conditions of the wind turbine and obtain the time domain response results of the parameterized model of the hybrid tower cable system.

[0034] A third aspect of the present invention provides a medium having a program stored thereon, which, when executed by a processor, implements the steps in the power demand response method based on a data matching algorithm as described in the first aspect of the present invention.

[0035] The fourth aspect of the present invention provides a device, including a memory, a processor, and a program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps in the parametric modeling method of a hybrid tower cable system based on multi-body dynamics as described in the first aspect of the present invention are implemented.

[0036] One or more of the above technical solutions have the following beneficial effects:

[0037] The present invention discloses a parameterized modeling method and system of a mixed tower cable system based on multi-body dynamics. Based on the theory of multi-body dynamics, a parameterized dynamic model of a mixed tower cable system is constructed, and the bonding force of each bonding surface of a concrete segment is simulated by means of spline curves and function expressions, and the compressive and non-tensile properties and fracture properties of concrete are equivalently simulated. Referring to the historical data information of aerodynamic loads acting on the wind turbine rotor of a wind turbine with the same power level, the aerodynamic load of the wind turbine rotor is equivalently calculated according to the power and speed information of the wind turbine to be analyzed, and applied to the wind turbine rotor using spline curves and function expressions, and the yaw condition of the wind turbine rotor is constructed by means of the driving constraint of the rotating pair, which effectively simulates the dynamic response characteristics of the wind turbine in continuous time under various complex working conditions.

[0038] The present invention proposes a parametric modeling method for a hybrid tower cable system based on multi-body dynamics, which realizes the response characteristic analysis of a hybrid tower under multiple working conditions and multiple structures. It can improve the efficiency and safety of structural design. The rationality of different design parameters can be verified through simulation, thereby guiding the design plan. The simulation model can accurately predict the response of the system under various complex working conditions, help discover potential structural failures and formulate maintenance plans.

[0039] Response data is also obtained by simulating complex working conditions. The hybrid tower cable system usually works in a changing environment. Dynamic simulation can simulate external disturbances under real conditions and provide data support for structural design under complex working conditions. It improves construction and operation and maintenance efficiency. Through simulation analysis, the optimal sequence and tension distribution of cable tensioning can be determined to improve construction efficiency. The dynamic model can be used as a basic tool for operating status monitoring and maintenance decision-making to evaluate operating status and formulate maintenance strategies.

[0040] Advantages of additional aspects of the present invention will be given in part in the following description, and in part will become obvious from the following description, or will be learned through practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] The accompanying drawings in the specification, which constitute a part of the present invention, are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations on the present invention.

[0042] Figure 1It is a scheme diagram of a method for creating a parameterized dynamic model of a hybrid tower cable system according to a first embodiment of the present invention;

[0043] Figure 2 This is a simplified model diagram of a mixing tower according to a first embodiment of the present invention;

[0044] Among them, 1. Wind turbine foundation, 2. Concrete section, 3. Transition section, 4. Steel tower section, 5. Steel strand.

[0045] Figure 3 This is a graph showing the relationship between the adhesive force and displacement of a certain section of a mixing tower in a certain direction in Example 1 of the present invention;

[0046] Figure 4 It is a schematic diagram of the pitch angle, roll angle, torsion angle and yaw angle of a mixed tower in accordance with the first embodiment of the present invention;

[0047] Figure 5 This is a dynamic response diagram of a mixing tower fracture failure mode according to the first embodiment of the present invention;

[0048] Figure 6 This is a dynamic response diagram of a mixed tower in Example 1 of the present invention when the pre-tensioning force of the steel cable is normal;

[0049] Figure 7 This is a graph showing the variation trend of the inclination angles of the concrete tower in the roll, pitch and torsion directions over time when the pre-tensioning force of the steel cables is normal according to the first embodiment of the present invention;

[0050] Figure 8 This is a graph showing the variation trend of the angular velocity of the rolling, pitching and torsion directions of the first hybrid tower in the embodiment of the present invention when the pre-tensioning force of the steel cables is normal over time;

[0051] Fig. 9 This is a graph showing the variation trend of the angular acceleration in the roll, pitch and torsion directions of the first hybrid tower in the embodiment of the present invention when the pre-tensioning force of the steel cables is normal over time;

[0052] Fig.10 is a graph showing a change trend of the preload force of each steel cable over time when the preload force of the steel cable is normal in the first embodiment of the present invention;

[0053] Fig.11 is a graph showing the variation trend of the inclination angles of the concrete tower in the roll, pitch and torsion directions over time when the pre-tightening force of the steel cables is insufficient in the first embodiment of the present invention; DETAILED DESCRIPTION

[0054] It should be noted that the following detailed descriptions are exemplary and are intended to provide further explanation of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meanings as those commonly understood by those skilled in the art to which the present invention belongs.

[0055] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present invention. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "include" and / or "include" are used in this specification, it indicates the presence of features, steps, operations, devices, components and / or their combinations;

[0056] Embodiment 1:

[0057] Due to the short time of hybrid towers on the market, the low proportion of in-service hybrid towers, the diversity of hybrid tower structures and the differences in wind resources in wind farms, there is little operating data on existing hybrid tower wind turbines. In particular, there is a serious lack of time and frequency domain response characteristic data on hybrid tower wind turbines under extreme operating conditions. This makes it difficult to effectively draw on the time and frequency domain response characteristic data of existing hybrid towers to guide the design of hybrid tower cable systems for other wind farms.

[0058] In view of the above problems, the present invention provides a parameterized modeling method of a hybrid tower cable system based on multi-body dynamics theory, which realizes the time domain response characteristic analysis of the hybrid tower under multiple working conditions and multiple structures. Figure 1 As shown, first, create a parametric model of the concrete segment, set the material density of the mixed tower for correction, establish contact between the joint surfaces of the concrete segment, and use the six-degree-of-freedom bushing force to construct the pressure-resistant and non-tensile properties of the bonding material. Then, create a parametric model of the steel cable. Create the rotation constraints of the wind turbine hub-nacelle and nacelle-tower, apply wind loads, perform dynamic simulation on the operating conditions of the wind turbine, and obtain the time domain response results of the parametric model of the mixed tower steel cable system.

[0059] The specific steps include:

[0060] Step 1: Obtain the geometric parameters of the concrete tower, create parameter points at the center of each end surface of the concrete segment according to the geometric parameters, and establish a three-dimensional model of the concrete tower.

[0061] In a specific embodiment, the geometric parameters of the mixing tower include the height of the mixing tower, the number of concrete sections, the radius of the mixing tower and the wall thickness of the mixing tower. Figure 2 As shown. The mixed tower is mainly composed of a wind turbine foundation 1, a concrete segment 2, a transition segment 3, a steel tower segment 4 and a steel strand 5. The mixed tower segment consists of three sections: a concrete segment, a transition segment and a steel tower segment. The concrete segment is assembled from a plurality of prefabricated annular segments. Therefore, it is necessary to consider the influence of the bonding surface of the plurality of prefabricated annular segments of the concrete segment. Specifically, after the modeling is completed, a parameter point is created at the center of each end face of the concrete segment. The mixed tower height, the number of concrete segments, the mixed tower radius and the mixed tower wall thickness are used as design variables to construct a full-scale model of the mixed tower cable system to meet the needs of analyzing mixed towers of different structural sizes.

[0062] Step 2: Correct the mass of each section of the mixing tower according to the material density of the mixing tower until the deviation between the mass of the mixing tower and the actual mass is less than a preset range.

[0063] In a specific implementation, the density property of the material is set based on experience, the mass of the three-dimensional model of the mixing tower is calculated by software simulation, and the calculated result is compared with the actual mass of the mixing tower, so as to check the deviation between the mass of the mixing tower and the actual mass, so that it is within the allowable deviation range, and ensure the accuracy of the model. At the same time, the density property of concrete is used as a design variable to meet the needs of analyzing mixing towers of different materials.

[0064] Step 3: Simulate the compressive but tensile properties of the concrete tower and establish contact and bonding relationships between the concrete segments of the concrete tower 3D model.

[0065] Step 3.1: Simulate the characteristics of the mixing tower that it is resistant to pressure but not to tension.

[0066] Step 3.1.1: Calculate the bonding force and bonding moment in the six directions of each bonding surface of the concrete segment based on the Young's modulus, strain and cross-sectional area of ​​the bonding material.

[0067] In a specific embodiment, the six-degree-of-freedom bushing force is used to equalize the bonding material between the concrete segments, and the bonding force and displacement, and the bonding torque and rotation angle relationships in the six degrees of freedom of each bonding surface of the concrete segments are obtained based on parameters such as Young's modulus, strain and concrete cross-sectional area of ​​the bonding material such as epoxy resin.

[0068] Among them, the calculation formula of the bonding strength of the bonding material is as follows:

[0069] ,

[0070] ,

[0071] .

[0072] in, is the bonding strength of the bonding material, is the stress of the bonding material, is the cross-sectional area of ​​the mixing tower, is the Young's modulus of the bonding material, is the strain of the bonding material, is the deformation length of the bonding material, is the original length of the bonding material.

[0073] In this embodiment, the bonding material is epoxy resin. In other possible embodiments, other bonding materials may also be used for bonding.

[0074] Step 3.1.2: Use spline curves and function expressions to express the relationship between bonding force and displacement, bonding torque and rotation angle. The relationship is as follows: Figure 3 As shown, the X-axis represents the displacement of a certain section of the mixing tower, and the Y-axis represents the bonding force. It can accurately simulate the compressive but not tensile properties of the mixing tower.

[0075] Step 3.1.3: Set the critical position to simulate the fracture failure characteristics of the mixing tower.

[0076] In a specific implementation, a critical position is set, such as when the displacement or rotation angle exceeds this value, the bonding force suddenly changes to 0, simulating the fracture failure characteristics of the mixing tower, such as Figure 4 (a) Figure 4 (b) and Figure 5 shown.

[0077] in, Figure 4 (a) is a schematic diagram of the pitch angle of the mixing tower. Figure 4 (b) is a schematic diagram of the yaw angle of the mixed tower. Figure 5 This is the dynamic response diagram of the mixed tower under the failure form of fracture. Figure 4 (a) Figure 4 As shown in (b), the origin of the following coordinate system (X'Y'Z') coincides with the origin of the cabin coordinate system (X''Y''Z''), the projection line of the following coordinate system X' axis to the XY plane of the geodetic coordinate system (XYZ) coincides with the projection line of the X'' axis of the cabin coordinate system to the XY plane of the geodetic coordinate system, and the X'Y' plane of the following coordinate system is always parallel to the XY plane of the geodetic coordinate system. The rotation angle of the following coordinate system relative to the geodetic coordinate system around the Z axis is the mixed tower yaw angle, as shown in Figure 4 As shown in (b) in the figure. The rotation angle of the cabin coordinate system relative to the follower coordinate system around the X' axis is the mixed tower roll angle. The rotation angle of the cabin coordinate system relative to the follower coordinate system around the Y' axis is the mixed tower pitch angle, as shown in Figure 4 As shown in (a) in the figure, the rotation angle of the cabin coordinate system relative to the follower coordinate system around the Z' axis is the mixed tower torsion angle.

[0078] Step 3.2: Establish contact and bonding relationships between concrete segments in the three-dimensional model of the concrete tower.

[0079] Step 3.2.1: Set contact relationships between concrete segments, between concrete segments and wind turbine foundations, and between concrete segments and transition segments.

[0080] In a specific implementation, contact is established between each concrete segment, between the concrete segment and the wind turbine foundation, and between the concrete segment and the transition segment. A fixing pair is set between the transition segment and the steel tower segment to simulate the pressure-resistant but tension-intolerant properties of the concrete tower, and fix the wind turbine foundation to the ground.

[0081] Step 3.2.2: Based on Hertz contact theory, the contact stiffness is determined according to the contact radius, elastic modulus of the material and Poisson's ratio. In order to express the contact between the end faces of each concrete tower segment, the contact force between the end faces provides support force and forms a balancing force with the gravity of each concrete tower segment. Therefore, the modeling simulates the actual situation so that the concrete segments are arranged together due to the effect of gravity.

[0082] In a specific implementation, the contact stiffness is determined based on the Hertz contact theory, according to the contact radius, the elastic modulus of the material, and the Poisson's ratio.

[0083] The contact stiffness is calculated as follows:

[0084] ,

[0085] ,

[0086] .

[0087] in, is the elastic modulus of concrete; is the Poisson’s ratio of concrete; is the contact radius. Since it is a planar contact, .

[0088] Step 4: Obtain the design parameters of the steel cables in the hybrid tower, construct a parametric model of the steel cables, and combine the parametric model of the steel cables with the three-dimensional model of the hybrid tower to obtain a parametric model of the hybrid tower steel cable system.

[0089] Step 4.1: Obtain the design parameters of the steel cables in the hybrid tower. The design parameters of the steel cables include the cable arrangement position, quantity, diameter, axial stiffness, damping, anchor point position and preload.

[0090] Step 4.2: Construct a parametric model of the steel cable. The parametric model of the steel cable is combined with the three-dimensional model of the mixed tower to obtain a parametric model of the mixed tower steel cable system.

[0091] In a specific implementation, due to the large number of steel cables in the tower, modeling is time-consuming and laborious. Therefore, in the modeling process, this embodiment simplifies the steel cable modeling and equates multiple strands of steel cables to one steel cable. In order to solve the problem of how to calculate the stiffness of the steel cable, its preload is also simplified by equivalently equating the steel cable preload with the axial force, and using the limit function in the dynamic simulation software to express the characteristic that the steel cable preload can only be generated when stretched, and at the same time, ensure that the steel cable preload changes linearly with the deformation.

[0092] Specifically, the moment generated by the wind load is calculated, and the cable preload generates a moment to offset this moment, thereby calculating the preload. The allowable inclination angle of the concrete tower in the design criteria is evenly divided into each concrete section, and the deformation of each cable is calculated, and then the cable stiffness is derived.

[0093] The calculation formula of cable stiffness is as follows:

[0094] ,

[0095] ,

[0096] ,

[0097] .

[0098] in, is the torque generated by the cable preload; is the torque produced by wind load; is the wind load; It is mixed Tata Gao; It is the preload of two cables symmetrical about the axis of the tower; It is the distance between two steel cables that are symmetrical about the axis of the tower; is the number of steel cables; is the cable preload; is the cable stiffness; is the elongation of the cable.

[0099] At the same time, the cable preload is also simplified equivalently, and the axial force is used to express the preload with the help of the limit function in the dynamic simulation software to simulate the characteristic that the cable preload can only be generated when stretched, and ensure that the preload changes linearly with the deformation.

[0100] Step 5: Perform dynamic simulation on the operating conditions of the wind turbine to obtain the time domain response results of the parameterized model of the hybrid tower cable system.

[0101] Step 5.1: Perform dynamic simulation on the operating conditions of the fan.

[0102] In a specific implementation, a wind rotor hub-nacelle and nacelle-tower rotation constraint is constructed, and a rotation drive is used to simulate wind turbine startup, wind yaw and other working conditions.

[0103] Step 5.1.1: Construct the rotation constraints of the wind rotor hub and nacelle, and use the rotation drive to simulate the wind turbine startup condition.

[0104] Referring to the historical data information of aerodynamic load acting on the wind rotor of wind turbines with the same power level, the wind rotor includes three blades. The aerodynamic load of the wind turbine rotor is equivalently calculated according to the aerodynamic characteristics, power and speed information of the blades of the wind turbine to be analyzed. The relationship among the torque coefficient, thrust coefficient and lateral force coefficient of the wind turbine rotor and the angle of attack is expressed by using spline curve. The axial force is used to apply the horizontal thrust and lateral force to the wind turbine rotor, and the calculation formula of the horizontal thrust and lateral force is expressed with the help of function expression.

[0105] The details are as follows:

[0106] From the fan speed information, the fan rotor rotation angular velocity can be calculated according to the following formula:

[0107] .

[0108] in, is the angular velocity of the wind turbine rotor, is the rotation speed, the unit is revolutions per minute.

[0109] Based on the wind turbine power generation and the wind turbine rotor rotation angular velocity information, the wind turbine rotor torque can be calculated according to the following formula:

[0110] .

[0111] The wind wheel torque T and torque coefficient C M , air density , the swept area S of the wind rotor, and the rotation radius R of the wind rotor, the outflow wind speed can be calculated according to the following formula :

[0112] .

[0113] According to the fan thrust coefficient C T Combined with the wind speed information calculated in the previous step, the horizontal thrust of the fan is calculated according to the following formula F T :

[0114] .

[0115] According to the fan lateral force coefficient C L Combined with the wind speed information calculated in the previous step, the fan lateral force is calculated according to the following formula: F L :

[0116]

[0117] in, C M , CT and C L are the torque coefficient, thrust coefficient and side force coefficient of the fan blade, which are related to the aerodynamic characteristics of the fan blade.

[0118] Step 5.1.2: Construct the rotation constraints of the nacelle and tower, and use the rotation drive to simulate the torsion effect of the yaw movement on the wind turbine tower. In order to effectively simulate the dynamic response characteristics of the mixed tower under different wind directions, the present invention sets the drive constraint of the rotation pair between the nacelle and the tower to construct the yaw working condition of the wind rotor.

[0119] Step 5.2: Obtain the time domain response results of the parameterized model of the hybrid tower cable system.

[0120] Step 5.2.1: Carry out multi-operating condition simulations of the hybrid tower cable system, including normal, extreme and fault conditions.

[0121] Step 5.2.2: Obtain the time domain response results of the four-degree-of-freedom swing angle, angular velocity, angular acceleration and cable preload of the hybrid tower.

[0122] In a specific embodiment, Figure 6 It represents the dynamic response of the mixed tower when the cable preload is normal. Figure 7 Indicates the change of mixing tower inclination angle under normal working conditions of preload. Figure 8 Indicates the change of angular velocity under normal working conditions of preload. Fig. 9 Indicates the change of angular acceleration under normal working conditions of preload, Fig.10 Indicates the change of preload force of the steel cable under normal working conditions. Fig.11 Indicates the change of mixing tower inclination angle under the condition of insufficient preload. Figure 6 , 7 , 8, 9, 10 and 11. The comparative analysis shows that the inclination angle of the mixed tower when the preload is insufficient is greater than that when the preload is normal, which leads to an increase in the deformation of the steel cable and an increase in the preload. If the allowable stress of the steel cable is exceeded, the steel cable will break and the tower will collapse, which is in line with the actual situation and the model is reasonable. It provides guidance for the design and later maintenance of the mixed tower.

[0123] Embodiment 2:

[0124] Embodiment 2 of the present invention provides a parameterized modeling system for a hybrid tower cable system based on multi-body dynamics, including:

[0125] A mixed tower model building module is configured to obtain geometric parameters of the mixed tower, create parameter points at the center of each end surface of the concrete segment according to the geometric parameters, and establish a three-dimensional model of the mixed tower;

[0126] The mass correction module is configured to correct the mass of each section of the mixing tower according to the material density of the mixing tower until the deviation between the mass of the mixing tower and the actual mass is less than a preset range;

[0127] The characteristic simulation module is configured to simulate the pressure-resistant and tension-resistant characteristics of the mixing tower, and establish contact and bonding relationships between the mixing tower sections in the three-dimensional model of the mixing tower;

[0128] A cable model building module is configured to obtain design parameters of the cables in the mixed tower, build a cable parametric model, and combine the cable parametric model with the mixed tower three-dimensional model to obtain a mixed tower cable system parametric model;

[0129] The dynamic simulation module is configured to perform dynamic simulation on the operating conditions of the wind turbine and obtain the time domain response results of the parameterized model of the hybrid tower cable system.

[0130] Embodiment three:

[0131] Embodiment 3 of the present invention provides a medium on which a program is stored. When the program is executed by a processor, the steps in the parameterized modeling method of a hybrid tower cable system based on multi-body dynamics as described in Embodiment 1 of the present invention are implemented.

[0132] Embodiment 4:

[0133] Embodiment 4 of the present invention provides a device, including a memory, a processor, and a program stored in the memory and executable on the processor. When the processor executes the program, the steps in the parametric modeling method of a mixed tower cable system based on multi-body dynamics as described in Embodiment 1 of the present invention are implemented.

[0134] The steps involved in the above embodiments 2, 3 and 4 correspond to the method embodiment 1. For the specific implementation methods, please refer to the relevant description part of embodiment 1.

[0135] Those skilled in the art should understand that the modules or steps of the present invention described above can be implemented by a general-purpose computer device, or alternatively, they can be implemented by a program code executable by a computing device, so that they can be stored in a storage device and executed by the computing device, or they can be made into individual integrated circuit modules, or multiple modules or steps therein can be made into a single integrated circuit module for implementation. The present invention is not limited to any specific combination of hardware and software.

[0136] Although the above describes the specific implementation mode of the present invention in conjunction with the accompanying drawings, it is not intended to limit the scope of protection of the present invention. Those skilled in the art should understand that various modifications or variations that can be made by those skilled in the art on the basis of the technical solution of the present invention without creative work are still within the scope of protection of the present invention.

Claims

1. A parametric modeling method for a hybrid tower cable system based on multi-body dynamics, characterized in that: The following steps are involved: Obtain the geometric parameters of the mixing tower, create parameter points at the center of each end surface of the mixing tower section according to the geometric parameters, and establish a three-dimensional model of the mixing tower; The mass of each section of the mixing tower is corrected according to the material density of the mixing tower until the deviation between the mass of the mixing tower and the actual mass is less than the preset range; Simulating the characteristics of the concrete tower that it is resistant to pressure but not to tension, a contact and bonding relationship is established between the concrete segments of the concrete tower three-dimensional model. Specifically, a contact relationship is set between the concrete segments, between the concrete segment and the wind turbine foundation, and between the concrete segment and the transition segment. Based on the Hertz contact theory, the contact stiffness is determined according to the contact radius, the elastic modulus of the material and the Poisson's ratio, and the concrete segments are arranged together due to the effect of gravity; According to the Young's modulus, strain and cross-sectional area of ​​the bonding material, the bonding force and bonding moment in the six directions of freedom of each bonding surface of the concrete segment are calculated; Use spline curves and function expressions to express the changing relationship between bonding force and displacement, bonding moment and rotation angle, and accurately simulate the compressive but not tensile properties of the mixed tower; Set critical positions to simulate the fracture failure characteristics of the mixing tower; The design parameters of the steel cables in the hybrid tower are obtained, and a parametric model of the steel cables is constructed. The parametric model of the steel cables is combined with the three-dimensional model of the hybrid tower to obtain a parametric model of the hybrid tower steel cable system. The operating conditions of the wind turbine are dynamically simulated to obtain the time domain response results of the parameterized model of the hybrid tower cable system.

2. The parameterized modeling method of a hybrid tower cable system based on multi-body dynamics according to claim 1 is characterized in that: The geometric parameters of the mixing tower include the mixing tower height, the number of concrete sections, the mixing tower radius and the mixing tower wall thickness.

3. The parameterized modeling method of a hybrid tower cable system based on multi-body dynamics according to claim 1 is characterized in that: The design parameters of the steel cables include cable arrangement position, number, diameter, axial stiffness, damping, anchor point position and preload.

4. The parameterized modeling method of a hybrid tower cable system based on multi-body dynamics according to claim 1 is characterized in that: The specific steps for dynamic simulation of the fan's operating conditions are as follows: Construct the rotation constraints of the wind rotor hub and the nacelle, and use the rotation drive to simulate the wind turbine startup condition. In this case, the aerodynamic load of the wind turbine rotor is equivalently calculated according to the power, blade aerodynamic characteristics and speed information of the wind turbine to be analyzed, and the aerodynamic load is applied to the wind turbine rotor using spline curves and function expressions to simulate the wind load of the wind turbine from startup to rated operating speed. The rotation constraints of the nacelle and tower are constructed, and the rotation drive is used to simulate the wind yaw conditions of the wind turbine under different wind directions.

5. The parameterized modeling method of a hybrid tower cable system based on multi-body dynamics according to claim 4 is characterized in that: The specific steps to obtain the time domain response results of the parameterized model of the hybrid tower cable system are as follows: Carry out multi-condition simulation of normal, extreme and fault conditions of hybrid tower cable system; The time domain response results of the four-degree-of-freedom swing angle, angular velocity, angular acceleration and cable preload of the hybrid tower are obtained.

6. A parametric modeling system for a hybrid tower cable system based on multi-body dynamics, characterized in that: include: A mixed tower model building module is configured to obtain geometric parameters of the mixed tower, create parameter points at the center of each end surface of the concrete segment according to the geometric parameters, and establish a three-dimensional model of the mixed tower; The mass correction module is configured to correct the mass of each section of the mixing tower according to the material density of the mixing tower until the deviation between the mass of the mixing tower and the actual mass is less than a preset range; The characteristic simulation module is configured to simulate the pressure-resistant and tension-resistant characteristics of the concrete tower, and to establish contact and bonding relationships between concrete segments of the three-dimensional model of the concrete tower. Specifically, contact relationships are set between concrete segments, between concrete segments and wind turbine foundations, and between concrete segments and transition segments. Based on the Hertz contact theory, the contact stiffness is determined according to the contact radius, the elastic modulus of the material and the Poisson's ratio, and the concrete segments are arranged together due to the effect of gravity; According to the Young's modulus, strain and cross-sectional area of ​​the bonding material, the bonding force and bonding moment in the six directions of freedom of each bonding surface of the concrete segment are calculated; Use spline curves and function expressions to express the changing relationship between bonding force and displacement, bonding moment and rotation angle, and accurately simulate the compressive but not tensile properties of the mixed tower; Set critical positions to simulate the fracture failure characteristics of the mixing tower; A cable model building module is configured to obtain design parameters of the cables in the mixed tower, build a cable parametric model, and combine the cable parametric model with the mixed tower three-dimensional model to obtain a mixed tower cable system parametric model; The dynamic simulation module is configured to perform dynamic simulation on the operating conditions of the wind turbine and obtain the time domain response results of the parameterized model of the hybrid tower cable system.

7. A computer-readable storage medium, characterized in that: A plurality of instructions are stored therein, and the instructions are suitable for being loaded by a processor of a terminal device and executing the parameterized modeling method of a hybrid tower cable system based on multi-body dynamics as described in any one of claims 1-5.

8. A terminal device, characterized in that: It includes a processor and a computer-readable storage medium, the processor is used to implement various instructions; the computer-readable storage medium is used to store multiple instructions, and the instructions are suitable for being loaded by the processor and executed by the parametric modeling method of the hybrid tower cable system based on multi-body dynamics described in any one of claims 1-5.

Images