CPC A01K 61/70 (2017.01) [G01P 5/001 (2013.01); G01P 5/241 (2013.01)] | 2 Claims |
1. A method for optimization design of an artificial reef structure, comprising:
step 1, arranging an artificial reef model to be tested on a test platform;
step 2, testing a flow field to obtain flow field data and a pull force of an artificial reef;
step 3, analyzing the flow field data to obtain a flow velocity reference point; and
step 4, carrying out optimization analysis in combination with the flow field data and the pull force, wherein
the testing a flow field to obtain flow field data and a pull force of an artificial reef in the step 2 comprises:
with the test platform comprising a water tank, wherein one end of the water tank is connected to a flow generator, the flow generator generates water flow having a fixed velocity, the generated water flow passes through the artificial reef, and the artificial reef is suspended on the water tank by means of a six-component force meter, such that a bottom of the artificial reef does not make contact with the water tank,
setting a measurement region, wherein the measurement region comprises a periphery of the artificial reef, and gridding the measurement region to obtain a plurality of grids;
measuring a water flow velocity around the artificial reef through an acoustic Doppler flow meter; and
measuring a gravity of the artificial reef by means of the six-component force meter, and recording the gravity as the pull force of the artificial reef under different flow velocities;
the testing a flow field to obtain flow field data and a pull force of an artificial reef further comprises:
starting, by the flow generator, to work, and setting the water flow velocity V0, wherein the water flow enters the water tank and flows through the artificial reef;
placing the acoustic Doppler flow meter in each grid of the measurement region to measure the water flow velocity, and using an average value of flow velocity data continuously recorded for 10 s at each point as the water flow velocity; and
with a size of the grid being 0.1 L×0.1 L, wherein L is a side length of the artificial reef, using a length of a side of the artificial reef perpendicular to a water flow direction as L under the condition that a length and width of the artificial reef are unequal, and recording a side of each artificial reef on a horizontal plane parallel to the water flow direction as a side A, and a side perpendicular to the water flow as a side B;
the analyzing the flow field data to obtain a flow velocity reference point in step 3 comprises:
obtaining the water flow velocity in each grid;
constructing a velocity contour of a cross section according to the water flow velocity; and
changing spacing between the artificial reefs to obtain velocity contours of the cross section under different spacing, and recording the velocity contours of the cross section as a velocity contour set VP, wherein
the set spacing at least comprises a value less than L, and the water flow velocity refers to a component of a velocity of the water flow in a direction parallel to a water surface;
the analyzing the flow field data to obtain a flow velocity reference point in step 3 specifically comprises:
obtaining a center position of a grid having a maximum water flow velocity and a corresponding water flow velocity value from the velocity contour set VP under different transverse spacing of the artificial reef, and recording the center position of the grid having the maximum water flow velocity and the corresponding water flow velocity value as maxPi and maxVi, wherein i is a variable, maxPi and maxVi represent a position of a point having the maximum water flow velocity obtained in the velocity contour VPi and the corresponding water flow velocity under i-th spacing respectively, i∈[1, len(VP)], len(VP) is a size of VP, i.e. the number of the set spacing; and VPi is an i-th velocity contour, such that a value of i is 1;
step 3.1, proceeding to step 3.2 under the condition that a value of max Vi/V0 is greater than a ratio of L(VPi) to L, and otherwise, proceeding to step 3.3 to obtain the flow velocity reference point; and recording L(VPi) as spacing between the artificial reefs under the i-th velocity contour;
step 3.2, increasing a value of i by 1 by using the position of maxPi as a reference point, proceeding to step 3.4 under the condition that i is greater than len(VP), and otherwise restarting step 3.1;
step 3.3, recording two opposite sides A in two artificial reefs as sides A1 of each artificial reef respectively, using end points of one ends of the two sides A1 close to a water flow incoming direction to make a line segment, increasing the value of i by 1 by using a midpoint of the line segment and a midpoint of the point maxPi as reference points, proceeding to step 3.4 under the condition of i>len(VP), and otherwise restarting step 3.1; and
step 3.4, obtaining flow velocity reference points under all artificial reef spacing; and
the carrying out optimization analysis in combination with the flow field data and the pull force in step 4 comprises:
in one velocity contour, constructing a closed region as a turbulence region according to a straight line in which two sides A1 are located and a straight line in which two sides B of one artificial reef are located separately;
using a closed region formed by a straight line of a side B of one artificial reef away from the water flow and two sides A and an edge of the water tank as a wake region of the flow artificial reef;
selecting a first velocity contour obtained from a flow field test in the velocity contour set VP when the spacing between the artificial reefs is less than L, obtaining a value of a maximum flow velocity point in the turbulence region as V1max, using maxV1 as a turbulence reference value Vref under the condition of V1max<maxV1, and otherwise, using a value of V1max as a turbulence reference value Vref, wherein max V1 is a maximum flow velocity value obtained in the first velocity contour;
recording the number of the spacing of the artificial reefs in the velocity contour set VP less than L as n1, and setting the value of i as n1+1; and initializing an empty set Vmean, and a point set Pmax;
proceeding to step 4.1 under the condition that a flow velocity reference point is provided in the turbulence region in the velocity contour VPi, and otherwise, proceeding to step 4.2;
step 4.1, recording a maximum point of a flow velocity in the turbulence region in the velocity contour VPi as Pmaxi, and recording a closed region formed by a minimum point of a flow velocity in two wake regions and Pmaxi as a first wake region;
step 4.1.1, removing an overlapping region from the first wake region as a new first wake region under the condition that the constructed first wake region overlaps any artificial reef, recording the first wake region as Ai and an average flow velocity of all sampling points in the first wake region as Vmeani, and putting Vmeani into Vmean and Pmaxi into a point set Pmax; and proceeding to step 4.3;
step 4.2, selecting the maximum point of the flow velocity from the velocity contour VPi, recording the maximum point as Pmaxi, constructing the first wake region according to an end point closest to Pmaxi from two artificial reefs and Pmaxi separately, recording the first wake region as Ai and the average flow velocity of all sampling points in the region as Vmeani, putting Vmeani into Vmean, and proceeding to step 4.3;
step 4.3, sequentially traversing points of a reference point set, and proceeding to step 4.1 or 4.2 to obtain the first wake region and the average flow velocity in the region;
step 4.4, setting an upwelling resistance limiting condition:
and
Vmeank<Vhk, wherein
Vmeank is an average flow velocity value of all measurement regions at k-th spacing, Vref is a turbulence reference value, UPF(Ak) is an upwelling volume of a region Ak, Wk is a pull force of the artificial reef at the k-th spacing, W0 is a pull force of the artificial reef at the water flow velocity of 0, exp( ) is an exponential function having a natural constant e as a base, Lk is a spacing value of the artificial reef at the k-th spacing, H is a height of the artificial reef, and Vhk is a vertical maximum velocity of the water flow in all the measurement areas at the k-th spacing; and
step 4.5, determining the upwelling resistance limiting condition when the spacing between the artificial reefs is greater than L, and obtaining a structure optimization direction of the artificial reef under the condition that the number of the artificial reefs satisfying the upwelling resistance limiting condition is greater than half len(VP): optimizing the flow velocity of the artificial reef structure; and otherwise, obtaining a structure optimization direction of the artificial reef: optimizing a resistance of the artificial reef structure.
|