Method, device and system for determining flutter critical wind speed of flexible support

Abstract

The invention discloses a method, device and system for determining the flutter critical wind speed of a flexible support. And considering the interaction of the supporting structures among the rows of the flexible bracket array, and constructing a fluid-solid coupling model of the flexible bracket array covering multiple rows of flexible brackets and photovoltaic panels on the flexible brackets. The calculation limitation of a single-row flexible support is broken, and the practical application value is high. A fluid-solid coupling model is divided into a solid domain model and a fluid domain model, a bidirectional fluid-solid coupling numerical simulation method is introduced, solving is performed based on an FEA solver and a CFD solver, real-time change of cable force of a cable system is considered, deformation of a photovoltaic panel is dynamically monitored, and accuracy is improved. And obtaining a torsion angle of the photovoltaic panel according to a solving result, and adjusting an inlet wind speed of the fluid domain model until the torsion angle is equal to a preset value. And obtaining the inlet wind speed at the moment as the flutter critical wind speed. Actual measurement and gas bomb wind tunnel tests are not needed, operation is easy and convenient, cost is low, the calculation period is short, and application in actual engineering design is facilitated.

Description

Technical Field

[0001] The present application relates to the technical field of photovoltaic brackets, and in particular to a method, device and system for determining a flutter critical wind speed of a flexible bracket. Background Art

[0002] Flexible brackets in the field of photovoltaic brackets have the characteristics of large span and strong terrain adaptability. Under the dual influence of static and dynamic wind loads, the flexible bracket will show wind-induced vibration effects on the structure, and will flutter in extreme weather and be easily damaged. In order to ensure the safety of the flexible bracket, the site wind speed needs to be lower than the flutter critical wind speed of the flexible bracket. Therefore, how to determine the flutter critical wind speed of the flexible bracket is crucial.

[0003] In the prior art, the critical flutter wind speed of the flexible bracket can be determined by field measurement and aeroelastic wind tunnel test methods, but field measurement is greatly affected by environmental factors, the results are not accurate, and the measurement cycle is long and the process is complicated. Aeroelastic wind tunnel test is costly, and the difficulty of making the test model and the requirements for test conditions are high. Therefore, how to accurately and simply determine the critical flutter wind speed of the flexible bracket has become an urgent problem to be solved. Summary of the invention

[0004] Based on the above problems, the present application provides a method, device and system for determining the critical wind speed of flutter of a flexible support, so as to accurately and simply determine the critical wind speed of flutter of the flexible support.

[0005] The present application discloses a method for determining a flutter critical wind speed of a flexible support, the method comprising:

[0006] Establishing a fluid-solid coupling model of a flexible support array; the flexible support array includes multiple rows of flexible supports and photovoltaic panels arranged on the flexible supports; the fluid-solid coupling model includes a solid domain model and a fluid domain model;

[0007] Solving the fluid-solid coupling model based on a finite element analysis (FEA) solver and a computational fluid dynamics (CFD) solver to obtain a torsion angle of the photovoltaic panel in the fluid-solid coupling model;

[0008] When the torsion angle is not equal to the preset value, adjusting the inlet wind speed of the fluid domain model until the torsion angle is equal to the preset value;

[0009] The inlet wind speed when the torsion angle is equal to the preset value is obtained as the flutter critical wind speed.

[0010] Optionally, the establishing of a fluid-structure coupling model of the flexible support array includes:

[0011] The solid domain model is set to be a plurality of rows of flexible brackets, each row of flexible brackets being provided with a plurality of photovoltaic panels;

[0012] According to the working scenario of the flexible support array, the inlet wind speed of the fluid domain model is set;

[0013] The solid domain model and the fluid domain model are meshed respectively to obtain a plurality of finite element units.

[0014] Optionally, solving the fluid-solid coupling model based on a finite element analysis (FEA) solver and a computational fluid dynamics (CFD) solver includes:

[0015] The FEA solver and the CFD solver form a bidirectional data transmission system, wherein the solution result of the solid domain model by the FEA solver and the solution result of the fluid domain model by the CFD solver are bidirectionally transmitted to realize bidirectional fluid-solid coupling.

[0016] Optionally, obtaining the torsion angle of the photovoltaic panel in the fluid-structure coupling model includes:

[0017] The bidirectional data transmission system outputs deformation data of the photovoltaic panel in each time step;

[0018] The displacement of the target feature point in the deformation data is obtained, and the torsion angle of the photovoltaic panel is obtained according to the displacement of the target feature point; the target feature point is located on the photovoltaic panel.

[0019] Optionally, the solution result of the solid domain model obtained by the FEA solver and the solution result of the fluid domain model obtained by the CFD solver are bidirectionally transmitted, including:

[0020] In a first time step, the solid domain model is solved by a FEA solver to obtain the displacement of the flexible support array as deformation data of the first time step;

[0021] Inputting the displacement into a CFD solver, solving the fluid domain, obtaining the force of the flexible support array, and inputting the force into a FEA solver;

[0022] In the second time step, the solid domain model is solved by the FEA solver based on the forces obtained in the first time step to obtain updated displacements of the flexible support array as deformation data of the second time step.

[0023] Optionally, after obtaining the torsion angle of the photovoltaic panel in the fluid-structure coupling model, the method further includes:

[0024] Continuously monitoring the displacement of the target feature point to obtain a displacement time history;

[0025] A torsion angle time history curve is output according to the displacement time history.

[0026] Optionally, before solving the fluid-solid coupling model based on a finite element analysis (FEA) solver and a computational fluid dynamics (CFD) solver, the method further includes:

[0027] The parameters of the FEA solver, the CFD solver, and the coupling of the bidirectional data transmission system are set.

[0028] Optionally, when the torsion angle is not equal to a preset value, adjusting the inlet wind speed of the fluid domain model includes:

[0029] When the torsion angle is greater than a preset value, reducing the inlet wind speed of the fluid domain model;

[0030] When the torsion angle is less than a preset value, the inlet wind speed of the fluid domain model is increased.

[0031] Based on the above method for determining the critical wind speed of flutter of a flexible support, the present application also discloses a device for determining the critical wind speed of flutter of a flexible support, including: a modeling unit, a solving unit and an adjusting unit;

[0032] The modeling unit is used to establish a fluid-solid coupling model of a flexible support array; the flexible support array includes multiple rows of flexible supports and photovoltaic panels arranged on the flexible supports; the fluid-solid coupling model includes a solid domain model and a fluid domain model;

[0033] The solving unit is used to solve the fluid-solid coupling model based on a finite element analysis (FEA) solver and a computational fluid dynamics (CFD) solver to obtain a torsion angle of the photovoltaic panel in the fluid-solid coupling model;

[0034] The adjustment unit is used to adjust the inlet wind speed of the fluid domain model until the torsion angle is equal to the preset value when the torsion angle is not equal to the preset value;

[0035] It is also used to obtain the inlet wind speed when the torsion angle is equal to the preset value as the flutter critical wind speed.

[0036] Optionally, the modeling unit includes:

[0037] A solid setting subunit, used to set the solid domain model to be a plurality of rows of flexible brackets, each row of flexible brackets being provided with a plurality of photovoltaic panels;

[0038] A wind speed setting subunit, used to set the inlet wind speed of the fluid domain model according to the working scenario of the flexible support array;

[0039] The division subunits are used to perform mesh division on the solid domain model and the fluid domain model respectively to obtain a plurality of finite element units.

[0040] Optionally, the solving unit includes:

[0041] The system forms a subunit, which is used for the FEA solver and the CFD solver to form a bidirectional data transmission system, wherein the solution results of the FEA solver on the solid domain model and the solution results of the CFD solver on the fluid domain model are bidirectionally transmitted to realize bidirectional fluid-solid coupling.

[0042] Optionally, the solving unit includes:

[0043] A deformation output subunit, used for the bidirectional data transmission system to output deformation data of the photovoltaic panel in each time step;

[0044] The angle acquisition subunit is used to acquire the displacement of the target feature point in the deformation data, and obtain the torsion angle of the photovoltaic panel according to the displacement of the target feature point; the target feature point is located on the photovoltaic panel.

[0045] Optionally, the system forms a subunit, including:

[0046] A first displacement calculation subunit is used to solve the solid domain model by using a FEA solver in a first time step to obtain the displacement of the flexible support array as deformation data of the first time step;

[0047] A force calculation subunit, used for inputting the displacement into a CFD solver, solving the fluid domain, obtaining the force of the flexible support array, and inputting the force into a FEA solver;

[0048] The second displacement calculation subunit is used to solve the solid domain model in a second time step through an FEA solver based on the force obtained in the first time step, and obtain the updated displacement of the flexible support array as the deformation data of the second time step.

[0049] Optionally, the device further comprises:

[0050] A time history calculation unit, used for continuously monitoring the displacement of the target feature point to obtain a displacement time history;

[0051] The curve acquisition unit is used to output a torsion angle time history curve according to the displacement time history.

[0052] Optionally, the device further comprises:

[0053] The setting unit is used to set the parameters of the FEA solver, the CFD solver, and the coupling of the bidirectional data transmission system.

[0054] Optionally, the adjustment unit includes:

[0055] A reduction subunit, used for reducing the inlet wind speed of the fluid domain model when the torsion angle is greater than a preset value;

[0056] The increasing subunit is used to increase the inlet wind speed of the fluid domain model when the torsion angle is less than a preset value.

[0057] Based on the above method for determining the flutter critical wind speed of a flexible support, the present application also discloses a system for determining the flutter critical wind speed of a flexible support, which is applied to the above method, and the system includes: a modeling module, a meshing module, a solver setting module, a solving module, a parameter acquisition module, a flutter instability determination module and a flutter critical wind speed output module;

[0058] The modeling module is used to establish a fluid-solid coupling model of the flexible support array; the fluid-solid coupling model includes a solid domain model and a fluid domain model;

[0059] The meshing module is used to mesh the solid domain model and the fluid domain model;

[0060] The solver setting module is used to set the coupling between the FEA solver, the CFD solver and the bidirectional data transmission system; the FEA solver and the CFD solver constitute the bidirectional data transmission system;

[0061] The solution module is used for solving the bidirectional data transmission system;

[0062] The parameter acquisition module is used to monitor the solution results, obtain the displacement of the characteristic points on the photovoltaic panel, and output the torsion angle according to the displacement of the characteristic points;

[0063] The flutter instability determination module is used to determine whether the torsion angle is equal to a preset value;

[0064] The flutter critical wind speed output module is used to output the inlet wind speed of the fluid domain model when the torsion angle is equal to a preset value as the flutter critical wind speed.

[0065] The present application discloses a method, device and system for determining the critical wind speed of flutter of a flexible support. Considering the integrity of the flexible support array, that is, the interaction of the support structure between rows, the present application constructs a fluid-solid coupling model of a flexible support array covering multiple rows of flexible supports and photovoltaic panels installed thereon. It surpasses the limitation of calculating the critical wind speed of flutter for a single row of flexible supports only, and has high practical application value. In addition, the fluid-solid coupling model is divided into a solid domain model and a fluid domain model, and a bidirectional fluid-solid coupling numerical simulation method is introduced. The fluid-solid coupling model is solved based on the FEA solver and the CFD solver, taking into account the real-time change of the cable force of the cable system under the action of wind load on the flexible support, and dynamically monitoring the deformation of the photovoltaic panel to improve accuracy. The torsion angle of the photovoltaic panel is obtained according to the solution result. When the torsion angle is not equal to the preset value, the inlet wind speed of the fluid domain model is adjusted until the torsion angle is equal to the preset value. The inlet wind speed at this time is obtained as the critical wind speed of flutter. The present application does not require actual measurement and aeroelastic wind tunnel test, is easy to operate, low in cost, and has a short calculation cycle, which is more conducive to application in actual engineering design. BRIEF DESCRIPTION OF THE DRAWINGS

[0066] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying any creative work.

[0067] Figure 1a A schematic flow chart of a method for determining a flutter critical wind speed of a flexible support disclosed in an embodiment of the present application;

[0068] Figure 1b It is a schematic diagram of the structure of the solid domain model disclosed in the embodiment of the present application;

[0069] Figure 1c A schematic diagram of the structure of the fluid domain model disclosed in the embodiment of the present application;

[0070] Figure 1d A schematic diagram of a torsion angle time history curve disclosed in an embodiment of the present application;

[0071] Figure 2 A schematic flow chart of another method for determining a flutter critical wind speed of a flexible support disclosed in an embodiment of the present application;

[0072] Figure 3 A schematic diagram of the structure of a device for determining a critical flutter wind speed of a flexible support disclosed in an embodiment of the present application;

[0073] Figure 4A schematic diagram of the structure of a system for determining a flutter critical wind speed of a flexible support disclosed in an embodiment of the present application. DETAILED DESCRIPTION

[0074] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0075] Embodiment 1: The present application discloses a method for determining the critical wind speed of flutter of a flexible support.

[0076] For details, please refer to Figure 1a A method for determining a flutter critical wind speed of a flexible support disclosed in this embodiment includes the following steps:

[0077] Step 101: Establish a fluid-structure interaction model of the flexible support array.

[0078] In the method described in this embodiment, the flexible bracket array includes multiple rows of flexible brackets and photovoltaic panels arranged on the flexible brackets. In the field of photovoltaic brackets, flexible brackets usually arrange inter-row support structures between rows to increase the flutter critical wind speed of the entire flexible bracket array. Based on the flexible bracket array composed of multiple rows of flexible brackets, the interaction between multiple rows, that is, the role of the inter-row support structure, can be considered. Compared with the modeling based only on a single row of flexible brackets, the method described in this embodiment is more comprehensive and has higher practical value.

[0079] In the method described in this embodiment, the constructed fluid-solid coupling model includes a solid domain model and a fluid domain model. Among them, the solid domain model includes the material properties, geometric cross-sections, boundary conditions, spatial positions and contact relationships between components of each component. These parameters are consistent with the flexible bracket array that the model needs to test the critical flutter wind speed in actual engineering. As an optional method, the solid domain model is set to multiple rows of flexible brackets, and multiple photovoltaic panels are arranged on each row of flexible brackets. For example, Figure 1b This is a schematic diagram of the structure of the solid domain model disclosed in the embodiment of this application, such as Figure 1bAs shown in the figure, the solid domain model is established according to the flexible bracket array, including end bracket 1, photovoltaic panel 2, panel frame 3, inter-row support structure 4, upper cable 5 and lower cable 6. The solid domain model is a five-row single-span structure, with 30 photovoltaic panels arranged in a single row and a total of 150 photovoltaic panels arranged in five rows. The overall east-west length is 50m, the north-south length is 20m, the height of the component panel from the ground is 5m, and the inclination angle of the component panel is 15°. Among them, the single row adopts a two-layer three-cable structure. The three flexible cables and the rigid internal support between the cables form a three-cable-truss structure system, which enhances the vertical bending stiffness of the structure. At the same time, the inter-row support structure between the rows connects the front and rear rows in series as a whole, greatly improving the wind resistance stability.

[0080] In the method described in this embodiment, the fluid domain model adopts the velocity inlet boundary condition. As an optional method, the inlet wind speed of the fluid domain model is set according to the working scenario of the flexible support array, that is, a certain wind speed value is input. For example, Figure 1c This is a schematic diagram of the structure of the fluid domain model disclosed in the embodiment of this application, such as Figure 1c As shown, the fluid domain model includes an inlet 7, an outlet 8, a bottom surface 9, a top surface 10, a front side surface 11, a rear side surface 12 and a fluid-solid coupling surface (the interface between the solid domain and the fluid domain) 13. The geometric dimensions are 200m×150m×50m, the inlet is 40m away from the flexible support array, the outlet is 90m away from the flexible support array, the front and rear side surfaces are both 75m away from the flexible support array, and the top surface is 45m away from the flexible support array. The blockage rate of the fluid domain model is 2.5%, which meets the blockage rate requirements of fluid simulation.

[0081] In the method described in this embodiment, meshing software can be used to mesh the solid domain model and the fluid domain model respectively, and the continuous geometric model can be discretized into a series of discrete units to obtain multiple finite element units, which is convenient for numerical solution. Among them, the appropriate mesh type and mesh density can be selected according to the requirements of bidirectional fluid-solid coupling.

[0082] Step 102: Solve the fluid-solid coupling model based on a finite element analysis (FEA) solver and a computational fluid dynamics (CFD) solver to obtain a torsion angle of the photovoltaic panel in the fluid-solid coupling model.

[0083] In the method described in this embodiment, the parameters of the finite element analysis FEA (Fi nite Element Analysis) solver and the computational fluid dynamics CFD (Computational Fluid Dynamics) solver are set, and the coupling of the two-way data transmission system formed by the FEA solver and the CFD solver is performed. Among them, the solid domain model is solved by the FEA solver, that is, the finite element method, and the fluid domain model is solved by the CFD solver, that is, the finite volume method.

[0084] In the method described in this embodiment, the solution results of the solid domain model by the FEA solver and the solution results of the fluid domain model by the CFD solver are bidirectionally transmitted to realize bidirectional fluid-solid coupling. Specifically, in the first time step, the solid domain model is solved by the FEA solver to obtain the displacement of the flexible support array as the deformation data of the first time step. The displacement is input to the CFD solver, and the fluid domain is solved to obtain the force of the flexible support array. The force is then input to the FEA solver, and the operation in the first time step is completed. Next, in the second time step, the solid domain model is solved again by the FEA solver based on the force obtained in the first time step to obtain the updated displacement of the flexible support array as the deformation data of the second time step.

[0085] As an optional method, each component of the flexible support array selects appropriate units for discretization in the FEA solver. For example, when computing resources are sufficient, each component can be solved using Solid units with large computational complexity and long computational time. When computing resources need to be saved, corresponding units can be selected according to the characteristics of each component. Figure 1b In the figure, the photovoltaic panel 2 can be solved by using Shell units, the end bracket 1, the component frame 3 and the inter-row support structure 4 can be solved by using Beam units, and the upper cable 5 and the lower cable 6 can be solved by using Link units. The selection of the above units is for example only, and appropriate solving units can be selected according to actual needs. The selection and type of solving units are not specifically limited here, and it is sufficient to complete the solution.

[0086] In addition, to ensure the wind resistance of the flexible support array, a certain amount of prestress is usually applied to the upper cable 5 and the lower cable 6 during installation. In the FEA solver, the same prestress as that of the actual flexible support array needs to be applied to the cable system.

[0087] As an optional method, CFD solver settings include solver selection, turbulence model settings, fluid material settings, boundary condition settings, dynamic grid settings, pressure-velocity coupling mode settings, discrete format settings, relaxation factor settings, residual settings, initialization settings, time step number and time step size settings, etc. Figure 1c In the figure, the inlet 7 adopts the velocity inlet boundary condition, with a preset velocity (inlet wind speed), the outlet 8 adopts the pressure outlet boundary condition, the bottom surface 9 and the fluid-solid coupling surface 13 adopt the no-slip wall boundary condition, and the top surface 10, the front side surface 11 and the rear side surface 12 all adopt the symmetric boundary condition. The selection of the above boundary conditions is for example only, and appropriate boundary conditions can be selected according to actual needs. The selection and type of boundary conditions are not specifically limited here, and the solution can be completed.

[0088] In the method described in this embodiment, the FEA solver and the CFD solver perform a bidirectional data transmission in each time step, and the bidirectional data transmission system outputs deformation data of a photovoltaic panel. As an optional method, a target feature point is selected on the photovoltaic panel, the displacement of the target feature point in the deformation data is obtained, and the torsion angle of the photovoltaic panel is obtained according to the displacement of the target feature point.

[0089] In addition, as an optional method, the displacement of the target feature point can be continuously monitored to obtain the displacement time history. The displacement time history is processed by Matlab to obtain the torsion angle time history curve of the photovoltaic panel.

[0090] Step 103: When the torsion angle is not equal to the preset value, adjusting the inlet wind speed of the fluid domain model until the torsion angle is equal to the preset value.

[0091] In the field of aeroelasticity, flutter instability refers to the phenomenon that when a structure is subjected to the action of a fluid, the vibration amplitude of the structure increases continuously due to the interaction between the fluid and the structure, and eventually the structure loses stability. Flutter instability is a self-excited vibration, that is, the vibration of the structure is not caused by external forced vibration, but is the result of the combined action of the structure's own dynamic response and fluid excitation.

[0092] In the method described in this embodiment, as an optional method, when the torsion angle is greater than a preset value, the flexible support array is deemed to have experienced flutter instability, and the inlet wind speed of the fluid domain model is reduced. When the torsion angle is less than the preset value, the flexible support array is deemed not to have experienced flutter instability, and the inlet wind speed of the fluid domain model is increased. This adjustment is continued until the torsion angle is equal to the preset value, and the flexible support array is deemed to have just experienced flutter instability.

[0093] As another optional method, it is also possible to determine whether the flexible support array has vibrated and become unstable by using the torsion angle time history curve obtained in step 102 . Figure 1dSchematic diagram of the torsion angle time history curve disclosed in the embodiment of the present application, such as Figure 1d As shown, the torsion angle shows a decreasing trend under the action of wind load, and it can be considered that the flexible support array has not experienced flutter instability at the current inlet wind speed. When the torsion angle shows an increasing trend under the action of wind load and reaches a preset angle (preset value), it can be considered that the flexible support array has experienced flutter instability at the current inlet wind speed. Figure 1d The torsion angle in the vibration changes periodically, and the torsion angle amplitude of each vibration cycle remains basically unchanged, about 12.5°. It can be considered that the flexible support array has reached the critical state of flutter instability.

[0094] The preset value can be set according to actual needs, and the specific value of the preset value is not limited here, as long as the method described in this embodiment can be implemented. For the adjustment (increase / decrease) of the inlet wind speed, the plus / minus value of each adjustment can be set according to actual needs, and the specific value of the plus / minus value is not limited here, as long as the method described in this embodiment can be implemented.

[0095] Step 104: obtaining the inlet wind speed when the torsion angle is equal to the preset value as the flutter critical wind speed.

[0096] In the method described in this embodiment, when the torsion angle is equal to the preset value, the flexible support array reaches the critical state of flutter instability, and the inlet wind speed at this time is the system flutter critical wind speed. When the inlet wind speed is less than the flutter critical wind speed, the flexible support system has no risk of flutter damage. When the inlet wind speed is greater than the flutter critical wind speed, the flexible support system will be damaged by flutter. Figure 1d As shown, the critical wind speed for flutter of the flexible bracket array is 40m / s. When the wind speed in the working scene is less than 40m / s, the structural form of the flexible bracket array can be adopted. When the wind speed in the working scene is greater than 40m / s, the flexible bracket array needs to be further improved to enhance its structural form.

[0097] The method described in this embodiment does not require cumbersome aeroelastic wind tunnel tests, which saves time and reduces the consumption of manpower and material resources, and is less costly. For each project, this method can provide customized calculations of the critical flutter wind speed, which not only has a short calculation cycle, but also has high results accuracy. The inter-row support structure between multiple rows of flexible brackets is taken into account, and is not limited to a single row of flexible brackets, and has high practical application value. It also incorporates considerations of real-time cable force changes in the flexible bracket cable system, performs dynamic detection, and has high calculation accuracy. At the same time, the meshing step can effectively improve the calculation efficiency while maintaining the calculation accuracy.

[0098] Embodiment 2: This application discloses another method for determining the critical wind speed of the flutter of the flexible support, please refer to Figure 2The method described in this embodiment introduces the whole process of determining the critical flutter wind speed of the flexible support.

[0099] Step 201: Based on the target flexible support array, a fluid-solid coupling model including a solid domain model and a fluid domain model is constructed.

[0100] Step 202: Setting the inlet wind speed of the fluid domain model.

[0101] Step 203: setting parameters of a bidirectional data transmission system formed by the FEA solver and the CFD solver, and solving the fluid-structure coupling model therethrough to obtain deformation data of the photovoltaic panel.

[0102] Step 204: Obtain the torsion angle of the photovoltaic panel according to the deformation data of the photovoltaic panel.

[0103] Step 205: Determine whether the torsion angle is not equal to a preset value. If yes, proceed to step 206. If no, proceed to step 207.

[0104] Step 206: Determine whether the torsion angle is greater than a preset value. If yes, proceed to step 208. If no, proceed to step 209.

[0105] Step 207: Output the current inlet wind speed as the flutter critical wind speed.

[0106] Step 208: Reduce the inlet wind speed of the fluid domain model, and return to step 203.

[0107] Step 209: Increase the inlet wind speed of the fluid domain model, and return to step 203.

[0108] Based on the method for determining the critical wind speed of the flutter of the flexible support disclosed in the above embodiment, this embodiment correspondingly discloses a device for determining the critical wind speed of the flutter of the flexible support. Figure 3 , the device for determining the flutter critical wind speed of the flexible support comprises: a modeling unit 301, a solving unit 302 and an adjusting unit 303;

[0109] The modeling unit 301 is used to establish a fluid-solid coupling model of a flexible support array; the flexible support array includes multiple rows of flexible supports and photovoltaic panels arranged on the flexible supports; the fluid-solid coupling model includes a solid domain model and a fluid domain model;

[0110] The solving unit 302 is used to solve the fluid-solid coupling model based on a finite element analysis (FEA) solver and a computational fluid dynamics (CFD) solver to obtain a torsion angle of the photovoltaic panel in the fluid-solid coupling model;

[0111] The adjusting unit 303 is used to adjust the inlet wind speed of the fluid domain model when the torsion angle is not equal to the preset value until the torsion angle is equal to the preset value;

[0112] It is also used to obtain the inlet wind speed when the torsion angle is equal to the preset value as the flutter critical wind speed.

[0113] Optionally, the modeling unit 301 includes:

[0114] A solid setting subunit, used to set the solid domain model to be a plurality of rows of flexible brackets, each row of flexible brackets being provided with a plurality of photovoltaic panels;

[0115] A wind speed setting subunit, used to set the inlet wind speed of the fluid domain model according to the working scenario of the flexible support array;

[0116] The division subunits are used to perform mesh division on the solid domain model and the fluid domain model respectively to obtain a plurality of finite element units.

[0117] Optionally, the solving unit 302 includes:

[0118] The system forms a subunit, which is used for the FEA solver and the CFD solver to form a bidirectional data transmission system, wherein the solution results of the FEA solver on the solid domain model and the solution results of the CFD solver on the fluid domain model are bidirectionally transmitted to realize bidirectional fluid-solid coupling.

[0119] Optionally, the solving unit 302 includes:

[0120] A deformation output subunit, used for the bidirectional data transmission system to output deformation data of the photovoltaic panel in each time step;

[0121] The angle acquisition subunit is used to acquire the displacement of the target feature point in the deformation data, and obtain the torsion angle of the photovoltaic panel according to the displacement of the target feature point; the target feature point is located on the photovoltaic panel.

[0122] Optionally, the system forms a subunit, including:

[0123] A first displacement calculation subunit is used to solve the solid domain model by using a FEA solver in a first time step to obtain the displacement of the flexible support array as deformation data of the first time step;

[0124] A force calculation subunit, used for inputting the displacement into a CFD solver, solving the fluid domain, obtaining the force of the flexible support array, and inputting the force into a FEA solver;

[0125] The second displacement calculation subunit is used to solve the solid domain model in a second time step through an FEA solver based on the force obtained in the first time step, and obtain the updated displacement of the flexible support array as the deformation data of the second time step.

[0126] Optionally, the device further comprises:

[0127] A time history calculation unit, used for continuously monitoring the displacement of the target feature point to obtain a displacement time history;

[0128] The curve acquisition unit is used to output a torsion angle time history curve according to the displacement time history.

[0129] Optionally, the device further comprises:

[0130] The setting unit is used to set the parameters of the FEA solver, the CFD solver, and the coupling of the bidirectional data transmission system.

[0131] Optionally, the adjusting unit 303 includes:

[0132] A reduction subunit, used for reducing the inlet wind speed of the fluid domain model when the torsion angle is greater than a preset value;

[0133] The increasing subunit is used to increase the inlet wind speed of the fluid domain model when the torsion angle is less than a preset value.

[0134] Based on the above method for determining the critical wind speed of flutter of a flexible support, the present application also discloses a system for determining the critical wind speed of flutter of a flexible support, which is applied to the above method. Figure 3 , the system for determining the flutter critical wind speed of the flexible support comprises: a modeling module 401, a meshing module 402, a solver setting module 403, a solving module 404, a parameter acquisition module 405, a flutter instability determination module 406 and a flutter critical wind speed output module 407;

[0135] The modeling module 401 is used to establish a fluid-solid coupling model of the flexible support array; the fluid-solid coupling model includes a solid domain model and a fluid domain model;

[0136] The meshing module 402 is used to mesh the solid domain model and the fluid domain model;

[0137] The solver setting module 403 is used to set the coupling between the FEA solver, the CFD solver and the bidirectional data transmission system; the FEA solver and the CFD solver constitute the bidirectional data transmission system;

[0138] The solution module 404 is used for solving the bidirectional data transmission system;

[0139] The parameter acquisition module 405 is used to monitor the solution results, obtain the displacement of the characteristic point on the photovoltaic panel, and output the torsion angle according to the displacement of the characteristic point; it is also used for the displacement time history of the characteristic point and output the torsion angle time history curve;

[0140] The flutter instability determination module 406 is used to determine whether the torsion angle is equal to a preset value; and is also used to analyze the torsion angle time history curve to determine whether the flexible support array has flutter instability at the current inlet wind speed;

[0141] The flutter critical wind speed output module 407 is used to output the inlet wind speed of the fluid domain model when the torsion angle is equal to a preset value as the flutter critical wind speed.

[0142] In the critical wind speed prediction method of the flexible support array based on bidirectional fluid-solid coupling applied by the system described in this embodiment, all operation steps can be integrated and embedded into the corresponding system module. This approach not only simplifies the operation process and reduces the difficulty of operation, but also significantly improves the intelligence level of the system, making it easier to apply and promote in actual engineering design.

[0143] The embodiments in this specification are described in a progressive manner. As for the devices disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the description is relatively simple, and the relevant parts can be referred to the description of the methods.

[0144] It should also be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the presence of other identical elements in the process, method, article or device including the elements.

[0145] The steps of the method or algorithm described in conjunction with the embodiments disclosed herein may be implemented directly using hardware, a software module executed by a processor, or a combination of the two. The software module may be placed in a random access memory (RAM), a memory, a read-only memory (ROM), an electrically programmable ROM, an electrically erasable programmable ROM, a register, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

[0146] The features described in the embodiments of this specification can be replaced or combined with each other, so that professional and technical personnel in this field can implement or use the present application.

[0147] The above description of the disclosed embodiments enables those skilled in the art to implement or use the present application. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present application. Therefore, the present application will not be limited to the embodiments shown herein, but will conform to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for determining the critical wind speed of flutter of a flexible support, characterized in that: include: Establishing a fluid-solid coupling model of a flexible support array; the flexible support array includes a plurality of rows of flexible supports and photovoltaic panels arranged on the flexible supports; The fluid-solid coupling model includes a solid domain model and a fluid domain model; Solving the fluid-solid coupling model based on a finite element analysis (FEA) solver and a computational fluid dynamics (CFD) solver to obtain a torsion angle of the photovoltaic panel in the fluid-solid coupling model; When the torsion angle is not equal to the preset value, adjusting the inlet wind speed of the fluid domain model until the torsion angle is equal to the preset value; The inlet wind speed when the torsion angle is equal to the preset value is obtained as the flutter critical wind speed.

2. The method according to claim 1, characterized in that The method of establishing a fluid-solid coupling model of a flexible support array includes: The solid domain model is set to be a plurality of rows of flexible brackets, each row of flexible brackets being provided with a plurality of photovoltaic panels; According to the working scenario of the flexible support array, the inlet wind speed of the fluid domain model is set; The solid domain model and the fluid domain model are meshed respectively to obtain a plurality of finite element units.

3. The method according to claim 1, characterized in that The method of solving the fluid-solid coupling model based on a finite element analysis (FEA) solver and a computational fluid dynamics (CFD) solver includes: The FEA solver and the CFD solver form a bidirectional data transmission system, wherein the solution result of the solid domain model by the FEA solver and the solution result of the fluid domain model by the CFD solver are bidirectionally transmitted to realize bidirectional fluid-solid coupling.

4. The method according to claim 3, characterized in that The obtaining of the torsion angle of the photovoltaic panel in the fluid-solid coupling model comprises: The bidirectional data transmission system outputs deformation data of the photovoltaic panel in each time step; The displacement of the target feature point in the deformation data is obtained, and the torsion angle of the photovoltaic panel is obtained according to the displacement of the target feature point; the target feature point is located on the photovoltaic panel.

5. The method according to claim 3, characterized in that: The solution result of the FEA solver on the solid domain model and the solution result of the CFD solver on the fluid domain model are bidirectionally transmitted, including: In a first time step, the solid domain model is solved by a FEA solver to obtain the displacement of the flexible support array as deformation data of the first time step; Inputting the displacement into a CFD solver, solving the fluid domain, obtaining the force of the flexible support array, and inputting the force into a FEA solver; In the second time step, the solid domain model is solved by the FEA solver based on the forces obtained in the first time step to obtain updated displacements of the flexible support array as deformation data of the second time step.

6. The method according to claim 4, characterized in that After obtaining the torsion angle of the photovoltaic panel in the fluid-solid coupling model, the method further includes: Continuously monitoring the displacement of the target feature point to obtain a displacement time history; A torsion angle time history curve is output according to the displacement time history.

7. The method according to claim 3, characterized in that Before solving the fluid-solid coupling model based on a finite element analysis (FEA) solver and a computational fluid dynamics (CFD) solver, the method further includes: The parameters of the FEA solver, the CFD solver, and the coupling of the bidirectional data transmission system are set.

8. The method according to claim 1, characterized in that When the torsion angle is not equal to a preset value, adjusting the inlet wind speed of the fluid domain model includes: When the torsion angle is greater than a preset value, reducing the inlet wind speed of the fluid domain model; When the torsion angle is less than a preset value, the inlet wind speed of the fluid domain model is increased.

9. A device for determining the critical wind speed of flutter of a flexible support, characterized in that: include: Modeling unit, solving unit and adjustment unit; The modeling unit is used to establish a fluid-solid coupling model of a flexible support array; the flexible support array includes multiple rows of flexible supports and photovoltaic panels arranged on the flexible supports; the fluid-solid coupling model includes a solid domain model and a fluid domain model; The solving unit is used to solve the fluid-solid coupling model based on a finite element analysis (FEA) solver and a computational fluid dynamics (CFD) solver to obtain a torsion angle of the photovoltaic panel in the fluid-solid coupling model; The adjustment unit is used to adjust the inlet wind speed of the fluid domain model until the torsion angle is equal to the preset value when the torsion angle is not equal to the preset value; It is also used to obtain the inlet wind speed when the torsion angle is equal to the preset value as the flutter critical wind speed.

10. A system for determining the critical wind speed of flutter of a flexible support, characterized in that: For implementing the method described in any one of claims 1 to 8, the system comprises: a modeling module, a meshing module, a solver setting module, a solving module, a parameter acquisition module, a flutter instability determination module and a flutter critical wind speed output module; The modeling module is used to establish a fluid-solid coupling model of the flexible support array; the fluid-solid coupling model includes a solid domain model and a fluid domain model; The meshing module is used to mesh the solid domain model and the fluid domain model; The solver setting module is used to set the coupling between the FEA solver, the CFD solver and the bidirectional data transmission system; the FEA solver and the CFD solver constitute the bidirectional data transmission system; The solution module is used for solving the bidirectional data transmission system; The parameter acquisition module is used to monitor the solution results, obtain the displacement of the characteristic points on the photovoltaic panel, and output the torsion angle according to the displacement of the characteristic points; The flutter instability determination module is used to determine whether the torsion angle is equal to a preset value; The flutter critical wind speed output module is used to output the inlet wind speed of the fluid domain model when the torsion angle is equal to a preset value as the flutter critical wind speed.

Images