Finite element modeling method for assembled retaining wall of transformer substation

Abstract

The invention relates to the technical field of finite element modeling, in particular to a finite element modeling method for an assembled retaining wall of a transformer substation, which comprises the following steps: building an experimental environment and the assembled retaining wall, and monitoring the assembled retaining wall to obtain displacement data of a prefabricated part; setting a plurality of load conditions, and obtaining a relative displacement characteristic group by utilizing the relative displacement conditions of the prefabricated parts under different load conditions, so as to divide the prefabricated parts and obtain displacement category values; and on the basis of the difference distribution situation and the relative displacement situation between the displacement category values of the adjacent prefabricated parts, the connection areas and boundary coefficients between the adjacent prefabricated parts are determined, and therefore the finite element size of the corresponding connection areas between the prefabricated parts in the finite element modeling process is adjusted, and a finite element model of the assembly type retaining wall is obtained. According to the method, the precision of finite element modeling of the fabricated retaining wall is improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of finite element modeling, in particular to a finite element modeling method of an assembled retaining wall of a substation. BACKGROUND

[0002] High-quality, high-efficiency, low-carbon and environmentally friendly construction of power transmission and transformation projects is a major development requirement in the field of green construction of power transmission and transformation projects. Especially in the complex geological environment of the Yellow River Basin, the traditional construction period is long, the construction process is complex, and the carbon emission is high. Therefore, it is necessary to develop a green and efficient assembly structure system of power facilities structures such as substation frame (support) and equipment assembly concrete foundation, retaining wall, and to reveal its stress mechanism and establish a structure green design and construction technology system, which has important theoretical significance and application value for ecological protection in the Yellow River Basin.

[0003] In the prior art of finite element analysis modeling of assembled retaining walls, the connection interface between the prefabricated components of the assembled retaining wall is usually simplified. This analysis method can achieve relatively ideal modeling results in the case of simple retaining wall structure and environment, but in the case of complex geological environment and limited construction conditions of the substation, the assembled retaining wall will present a complex three-dimensional structure to adapt to the terrain conditions. Therefore, when using a simplified model for finite element analysis, the friction and slip characteristics of the contact surface are not fully considered, which leads to inaccurate prediction of stress distribution and deformation in the joint area, resulting in unsatisfactory finite element modeling results. SUMMARY

[0004] The present application provides a finite element modeling method for an assembled retaining wall of a substation to solve the existing problems.

[0005] The finite element modeling method for an assembled retaining wall of a substation provided by the present application adopts the following technical scheme: One embodiment of the present application provides a finite element modeling method for an assembled retaining wall of a substation, which comprises the following steps: An experimental environment and an assembled retaining wall are set up, and the assembled retaining wall is monitored in real time by sensors to obtain displacement data corresponding to each prefabricated component in the assembled retaining wall; A plurality of load conditions are set, and the relative displacement characteristics of the prefabricated components are obtained by using the relative displacement of the displacement data of different prefabricated components in the assembled retaining wall under all load conditions; the prefabricated components are classified by the relative displacement characteristics of the prefabricated components in the assembled retaining wall, and the displacement classification value of each prefabricated component is obtained; The connecting area between the adjacent prefabricated components and the boundary coefficient of the connecting area are determined based on the difference distribution between the target prefabricated component and the displacement category value of the adjacent prefabricated component, and in combination with the relative displacement between the adjacent prefabricated components. The finite element size of the corresponding connecting area between the prefabricated components in the finite element modeling process is adjusted by using the boundary coefficient, so as to obtain the finite element model of the assembled retaining wall.

[0006] Optionally, the method comprises the following steps of: A plurality of load conditions are set, and displacement data of all prefabricated components of the assembled retaining wall under each load condition is obtained. Under any load condition, the relative displacement factor of the target prefabricated component under the load condition is calculated according to the relative displacement shown by the displacement data between the target prefabricated component and other prefabricated components connected thereto. The relative displacement characteristic group of the target prefabricated component is obtained through the numerical distribution of the relative displacement factor of the target prefabricated component under all load conditions and the correlation degree with the load conditions.

[0007] Optionally, the method comprises the following steps of: In the experimental environment, a hydraulic loading system is used to apply graded loads, and the horizontal load, the vertical load and the dynamic load with increasing degrees are applied to the assembled retaining wall respectively, so as to obtain a plurality of load conditions.

[0008] Optionally, the relative displacement factor of the target prefabricated component under the load condition comprises the following specific method: Under any load condition, the range of displacement data of any prefabricated component is obtained as the displacement amplitude of the prefabricated component, and the displacement amplitude difference between any prefabricated component connected to the target prefabricated component and the target prefabricated component is obtained; the interval formed by the maximum value and the minimum value of the displacement data of any prefabricated component is recorded as the displacement interval of the prefabricated component, and the intersection-union ratio of the displacement intervals corresponding to any prefabricated component connected to the target prefabricated component and the target prefabricated component is recorded as the displacement interval overlap rate; the Pearson correlation coefficient of the displacement data corresponding to any prefabricated component connected to the target prefabricated component and the target prefabricated component is obtained, and the relative displacement factor of the target prefabricated component under any load condition is calculated in combination with the displacement amplitude difference and the displacement interval overlap rate.

[0009] Optionally, the relative displacement characteristic group of the target prefabricated component comprises the following specific method: According to the numerical distribution of the relative displacement factor of the target precast component under all load conditions, the relative displacement factor of the target precast component is calculated; according to the correlation degree of the relative displacement factor and the load condition, the displacement-load correlation degree is calculated; the array formed by the relative displacement factor and the displacement-load correlation degree is taken as the relative displacement feature group of the target precast component.

[0010] Optionally, the relative displacement feature group of the precast component in the assembled retaining wall is used for classifying the precast components to obtain the displacement category value of each precast component, and the specific method comprises the following steps: The Euclidean distance between the relative displacement feature groups of the precast components is taken as the distance measurement method of the DBSCAN clustering algorithm, and the DBSCAN clustering algorithm is used for clustering processing of all the precast components to obtain a plurality of clustering clusters; The displacement category value of the precast component in each clustering cluster is obtained through the relative displacement level of all the precast components in each clustering cluster, and for any clustering cluster, the product of the relative displacement factor and the displacement-load correlation degree of any precast component in the clustering cluster is taken as the displacement level coefficient of the precast component; The average displacement level coefficient of all the precast components in the clustering cluster is taken as the displacement type parameter of the clustering cluster, the displacement type parameters of all the clustering clusters are normalized by maximum and minimum, and the normalized result is taken as the displacement category value of the precast component in the corresponding clustering cluster.

[0011] Optionally, the difference distribution between the displacement category values of the target precast component and the adjacent precast components is combined with the relative displacement between the adjacent precast components to determine the connection area between the adjacent precast components and the boundary coefficient of the connection area, and the specific method comprises the following steps: The neighborhood conflict coefficient of the target precast component is obtained based on the difference between the displacement category values of the target precast component and the adjacent precast components; According to the relative displacement between any adjacent precast components, a bell-shaped curve is established, and the connection area between the adjacent precast components is obtained by using the bell-shaped curve; Any connection area is taken as a target connection area, two precast components constituting the connection area are obtained, the information entropy of the neighborhood conflict coefficient corresponding to all the connection areas formed by the two precast components is obtained, and the boundary coefficient of the target connection area is obtained by combining the neighborhood conflict coefficient of the target connection area with the information entropy.

[0012] Optionally, the specific method for obtaining the neighborhood conflict coefficient of the target precast component comprises the following steps: The absolute value of the difference between the displacement category values ​​of the target precast component and any adjacent precast component is recorded as the neighborhood conflict factor between the target precast component and the adjacent precast component. The average value of the neighborhood conflict factors between the target precast component and all adjacent precast components is taken as the neighborhood conflict coefficient of the target precast component.

[0013] Optionally, the specific method for establishing a bell-shaped curve based on the relative displacement between any adjacent prefabricated components and using the bell-shaped curve to obtain the connection area between the adjacent prefabricated components includes: Select adjacent precast components and prefabricated components Obtain prefabricated components respectively and prefabricated components Relative displacement data under all load conditions; along the precast component and prefabricated components The normal direction of the interface between the two points divides the interface into several equally spaced measurement segments. On each measurement segment, measurements are evenly distributed. The system uses 10 monitoring points to obtain the relative displacement values ​​of each point under different load conditions; for each monitoring point, it calculates the average relative displacement value under all load conditions. The preset number of monitoring points; All monitoring points are sorted according to their position coordinates on the connection interface; with position coordinates as independent variables and average relative displacement as dependent variables, a bell curve in the form of a Gaussian distribution is fitted. Determine the effective range of the bell curve and define the boundary position of the connecting region as follows. ,in, The default boundary factor is denoted as , where The mathematical expectation of the bell curve; It represents the standard deviation of the bell curve; the region between the two boundary positions is taken as the boundary region.

[0014] Optionally, the specific method for adjusting the finite element size of the corresponding connection area between prefabricated components during the finite element modeling process using boundary coefficients to obtain the finite element model of the prefabricated retaining wall includes: The basic finite element mesh size parameters are preset, and the finite element mesh size in the boundary region is obtained by combining the bell curve and boundary coefficient of the boundary region; A 3D model of the prefabricated retaining wall is created using the 3D point cloud data. The custom mesh control function of the finite element software is then used to perform mesh generation in the prefabricated retaining wall model according to the finite element mesh size distribution, thus obtaining the finite element model of the prefabricated retaining wall.

[0015] The beneficial effects of the technical scheme of the present application are: through self-adaptive mesh division, the fine simulation of the key connection area of the assembled retaining wall is realized, and at the same time, the waste of computing resources caused by high-density mesh division of the overall structure is avoided, compared with the traditional uniform mesh division method, the finite element model of the present application can reduce the number of calculation units under the premise of ensuring the calculation accuracy, significantly improve the finite element analysis efficiency, and is especially suitable for the mechanical property analysis of the assembled retaining wall of the transformer substation with complex geological environment and various structural forms. BRIEF DESCRIPTION OF DRAWINGS

[0016] In order to more clearly illustrate the technical scheme in the embodiments of the present application or the prior art, the drawings needed to be used in the embodiments or the prior art description will be briefly introduced below, and obviously, the drawings in the following description can only be some embodiments of the present application, and other drawings can also be obtained by those skilled in the art without creative labor.

[0017] Figure 1 The step flow chart of the finite element modeling method of the assembled retaining wall of the transformer substation according to the present application; Figure 2 The structural block diagram of the finite element modeling system of the assembled retaining wall of the transformer substation according to the present application. DETAILED DESCRIPTION

[0018] In order to further explain the technical means and effects adopted by the present application to achieve the predetermined invention purpose, the specific implementation, structure, features and effects of the finite element modeling method of the assembled retaining wall of the transformer substation according to the present application will be described in detail below in combination with the drawings and preferred embodiments. In the following description, different "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. In addition, the specific features, structures or characteristics in one or more embodiments can be combined in any suitable form.

[0019] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which the present application belongs.

[0020] The specific scheme of the finite element modeling method of the assembled retaining wall of the transformer substation provided by the present application will be specifically described below in combination with the drawings.

[0021] Please refer to Figure 1 which shows the step flow chart of the finite element modeling method of the assembled retaining wall of the transformer substation according to one embodiment of the present application, and the method comprises the following steps: Step S001: build an experimental environment and an assembled retaining wall, and monitor the assembled retaining wall in real time through a sensor to obtain displacement data corresponding to a plurality of prefabricated components in the assembled retaining wall, respectively.

[0022] It should be noted that in the substation engineering construction under complex geological conditions and limited construction environment, the assembled retaining wall needs to adapt to the variable terrain and bear various external load actions. In such an environment, due to different connection methods between the precast components of the retaining wall, different degrees of relative displacement and interface friction slip phenomena often occur in the actual stress process. As shown in the figure, the traditional finite element modeling method is too simplified in handling the connection interface, and it is difficult to accurately reflect the stress distribution and deformation characteristics of the joint area. The present application focuses on such engineering practice, and for the assembled retaining wall system with complex three-dimensional structure form and multi-component splicing form in the substation, the interaction behavior between the components under real service conditions is considered, combined with experimental monitoring and load response characteristics, a refined finite element model that can truly reflect the mechanical properties of the contact surface is constructed, to improve the structural analysis accuracy and design reliability.

[0023] Specifically, in order to realize the finite element modeling method of the assembled retaining wall of the substation proposed in this embodiment, first of all, an experimental environment and an assembled retaining wall need to be built, and displacement data of precast fittings in the assembled retaining wall need to be collected, and the specific process is as follows: Step S101, building an assembled retaining wall to form an experimental environment.

[0024] As an example, the method of building an assembled retaining wall to form an experimental environment includes: assembling an assembled retaining wall physical model according to the design drawing in a controlled experimental site simulating the actual working condition of the substation. The experimental environment needs to reproduce the typical geological characteristics of the Yellow River basin (such as slope terrain with slope ≥15°, sandy clay layer distribution), and strictly control the construction process parameters (such as precast component assembly precision error ≤3mm), to ensure the consistency of the model and the actual engineering.

[0025] Step S102, using a high-precision three-dimensional laser scanner to non-contact scan all precast components of the retaining wall.

[0026] As an example, the method of using a high-precision three-dimensional laser scanner to non-contact scan all precast components of the retaining wall includes: the model of the high-precision three-dimensional laser scanner is Faro Focus S350, the ranging accuracy is ±1mm, and the point cloud density is ≥50 points / mm 2; The original point cloud data is denoised, registered and encapsulated by a matching software (such as Geomagic Control X) to generate a high-fidelity three-dimensional digital model (.stl format) with a resolution controlled within 1mm. The prefabricated component at least includes a wall panel, a bottom plate and a buttress column, and the geometric size data covers key parameters: wall panel thickness (150-300mm), bottom plate width (1.5-3.0 times the wall height), buttress column spacing (3-6m), and micro-topography of the connecting interface (such as keyway depth, surface roughness Ra value).

[0027] As an example, when the retaining wall is scanned, the scanning process is performed in regions, specifically: First, fix the scanner at a preset reference point (the origin of the coordinate system is aligned with the design coordinate system of the substation), and perform full-surface scanning of the wall panel (including vertical or inclined panels), focusing on capturing the panel curvature, reserved connecting hole position and edge chamfer details.

[0028] Secondly, move the scanner to the bottom plate area, and scan the bottom surface roughness of the bottom plate (including anti-skid keyway and drainage hole) and the contact interface with the wall panel.

[0029] Finally, the buttress column (including column foot embedded parts and column body reinforcing ribs) is scanned around to ensure that the geometric characteristics of the column-plate connecting node are completely captured.

[0030] Step S103, collecting displacement data of the prefabricated component.

[0031] As an example, the specific method of collecting monitoring data is: first, install high-sensitivity displacement sensors (model: Micro-Epsilon eddyNCDT 3300, range ±50mm, accuracy 0.01mm) at key positions of the wall panel (50mm away from the connecting interface), the bottom plate (keyway area), and the buttress column (column foot and column top).

[0032] Then, the sensor arrangement density is graded according to the importance of the component: ≥4 / m 2 (high density area) in the connecting interface area, and ≥1 / m 2 (low density area) in the non-critical area, to ensure that the relative displacement sensitive points are covered, and the real-time displacement data of each sensor is synchronously collected under the set load condition (sampling frequency 10Hz).

[0033] At this point, the three-dimensional point cloud data of the assembled retaining wall and the displacement data corresponding to each prefabricated component in the assembled retaining wall are obtained by the above method.

[0034] Step S002: setting several load conditions, obtaining the relative displacement characteristics group of the prefabricated components by using the relative displacement conditions shown by the displacement data of different prefabricated components in the assembled retaining wall under all load conditions; and classifying the prefabricated components by using the relative displacement characteristics group of the prefabricated components in the assembled retaining wall, to obtain the displacement category value of each prefabricated component.

[0035] It should be noted that the monitoring data of the assembled retaining wall under the corresponding load condition is obtained by adjusting the load conditions to test the assembled retaining wall in the experimental environment. Since the relative displacement conditions between the prefabricated components of the assembled retaining wall under different load conditions are different, in order to facilitate subsequent finite element modeling, the prefabricated components are respectively represented by the matched model accuracy, the embodiment of the present application selects to classify the displacement of the prefabricated components under different load conditions.

[0036] Specifically, in step S201, several load conditions are set, and the relative displacement characteristics group of the prefabricated components is obtained by using the relative displacement conditions shown by the displacement data of different prefabricated components in the assembled retaining wall under all load conditions.

[0037] It should be noted that since the assembled retaining wall is combined and spliced by different prefabricated components through corresponding connection modes, under different load conditions, the assembled retaining wall composed of prefabricated components with different connection modes will have corresponding system-level responses, that is, the different prefabricated components will have a certain degree of displacement, and the direction and degree of displacement will be different, thereby presenting the relative displacement condition. Therefore, in order to effectively describe the displacement condition of the prefabricated component, the embodiment of the present application selects to obtain the relative displacement characteristics group of the prefabricated component.

[0038] As a preferred embodiment, the specific method for obtaining the relative displacement characteristics group of the prefabricated component includes: First, several load conditions are set, and the displacement data of all prefabricated components of the assembled retaining wall under each load condition is obtained.

[0039] As an optional embodiment, the specific method for setting several load conditions includes: in the experimental environment, a hydraulic loading system is used to apply graded loads, and the assembled retaining wall is respectively applied with increasing horizontal load, vertical load and dynamic load, to obtain several load conditions.

[0040] As an optional embodiment, the horizontal load, vertical load and dynamic load are respectively: the horizontal load is simulated soil pressure, which is increased by three levels of 0.2H, 0.5H and 0.8H (H is the wall height), and each level is held for 10 minutes; the vertical load is simulated equipment additional load, and the applied range is 50-200kN / m2 Dynamic load: simulate earthquake working condition by vibration table, input El-Centro wave (peak acceleration 0.3g).

[0041] Then, under any load condition, according to the relative displacement shown by the displacement data between the target precast component and other precast components connected thereto, the relative displacement factor of the target precast component under the load condition is calculated.

[0042] As an optional embodiment, the specific method for obtaining the relative displacement factor of the target precast component under any load condition is as follows: under any load condition, the range of displacement data of any precast component is obtained as the displacement amplitude of the precast component, and the difference between the displacement amplitude of any precast component connected to the target precast component and the displacement amplitude of the target precast component is obtained; the interval formed by the maximum value and the minimum value of the displacement data of any precast component is recorded as the displacement interval of the precast component, and the intersection-union ratio of the displacement intervals corresponding to any precast component connected to the target precast component and the target precast component is recorded as the displacement interval overlap rate; the Pearson correlation coefficient of the displacement data corresponding to any precast component connected to the target precast component and the target precast component is obtained, and the relative displacement factor of the target precast component under any load condition is calculated in combination with the displacement amplitude difference and the displacement interval overlap rate.

[0043] As an optional embodiment, the specific calculation method for the relative displacement factor of the target precast component under any load condition is as follows: wherein, represents the relative displacement factor of the target precast component under the i-th load condition; represents the relative displacement factor of the target precast component under the i-th load condition; represents the number of other precast components connected to the target precast component; represents the displacement amplitude difference between the i-th precast component connected to the target precast component and the target precast component under the i-th load condition; represents the displacement amplitude difference between the i-th precast component connected to the target precast component and the target precast component under the i-th load condition; represents the displacement interval overlap rate between the i-th precast component connected to the target precast component and the target precast component under the i-th load condition; represents the displacement interval overlap rate between the i-th precast component connected to the target precast component and the target precast component under the i-th load condition; represents the displacement interval overlap rate between the i-th precast component connected to the target precast component and the target precast component under the i-th load condition; represents the displacement interval overlap rate between the i-th precast component connected to the target precast component and the target precast component under the i-th load condition; represents the Pearson correlation coefficient of the displacement data between the i-th precast component connected to the target precast component and the target precast component under the i-th load condition. represents the Pearson correlation coefficient of the displacement data between the i-th precast component connected to the target precast component and the target precast component under the i-th load condition.

[0044] ​It should be noted that the relative displacement factor is used to describe the degree of relative displacement between the corresponding prefabricated component and the connected prefabricated component after the corresponding load condition is applied. The greater the value of the relative displacement factor, the greater the degree of relative displacement between the corresponding prefabricated component and the connected prefabricated component.

[0045] Finally, the relative displacement characteristics group of the target prefabricated component is obtained through the value distribution of the relative displacement factor of the target prefabricated component under all load conditions and the correlation degree with the load conditions.

[0046] As an optional embodiment, the specific acquisition method of the relative displacement characteristics group of the target prefabricated component is: calculating the relative displacement coefficient of the target prefabricated component according to the value distribution of the relative displacement factor of the target prefabricated component under all load conditions; calculating the displacement-load correlation degree according to the correlation degree of the relative displacement factor and the load condition; and taking the array formed by the relative displacement coefficient and the displacement-load correlation degree as the relative displacement characteristics group of the target prefabricated component.

[0047] As an optional embodiment, the specific calculation method of the relative displacement characteristics group of the target prefabricated component is: wherein, represents the relative displacement coefficient of the target prefabricated component; represents the displacement-load correlation degree of the target prefabricated component; represents the exponential function with the natural constant as the base number; represents the standard deviation of the relative displacement factor of the target prefabricated component under all load conditions; represents the average value of the relative displacement factor of the target prefabricated component under all load conditions; is the sequence of the relative displacement factor and the sequence of the load intensity ; is the standard deviation of the relative displacement factor; is the standard deviation of the load intensity.

[0048] It should be noted that the relative displacement characteristics group can comprehensively reflect the relative displacement behavior characteristics of the prefabricated component under multiple load conditions, and provide a quantitative basis for subsequent component classification and grid division. It should be noted that the relative displacement characteristics group reflects the corresponding prefabricated component.

[0049] Step S202: classifying the prefabricated components through the relative displacement characteristics group of the prefabricated components in the assembled retaining wall, and obtaining the displacement category value of each prefabricated component.

[0050] It should be noted that in the assembled retaining wall, when the prefabricated components are in similar load environment under the same load condition, the relative movement types of the prefabricated components will show similar characteristics, and therefore, the embodiment of the present application classifies the characteristics shown by the prefabricated components to determine the relative displacement categories of the prefabricated components which are relatively stable and are not affected by various load conditions.

[0051] Firstly, the Euclidean distance between the relative displacement feature groups of the prefabricated components is taken as the distance measurement method of the DBSCAN clustering algorithm, and the DBSCAN clustering algorithm is used to cluster all the prefabricated components to obtain a plurality of clustering clusters.

[0052] Then, the displacement category value of the prefabricated components in each clustering cluster is obtained through the relative displacement level of all the prefabricated components in each clustering cluster. For any clustering cluster, the product of the relative displacement coefficient and the displacement-load correlation degree of any prefabricated component in the clustering cluster is taken as the displacement level coefficient of the prefabricated component.

[0053] Finally, the average displacement level coefficient of all the prefabricated components in the clustering cluster is taken as the displacement type parameter of the clustering cluster, the displacement type parameters of all the clustering clusters are normalized by maximum and minimum, and the normalized result is taken as the displacement category value of the prefabricated components in the corresponding clustering cluster.

[0054] It should be noted that the displacement category value reflects the stability and sensitivity of the relative displacement behavior of the prefabricated components under various loads. The larger the value, the higher the possibility of significant, load-sensitive relative displacement between the component and its adjacent components during the stress process, and the more unstable the connection interface. Conversely, it indicates that the displacement behavior is coordinated and the load response is stable. Therefore, the displacement category value provides a quantitative basis for subsequent judgment of whether the connection interface needs to be refined modeling, and is a key input parameter for realizing adaptive mesh division.

[0055] At this point, a plurality of categories of prefabricated components and the displacement category value of each prefabricated component are obtained by the above method.

[0056] Step S003: based on the difference between the displacement category value of the target prefabricated component and the displacement category value of the adjacent prefabricated component, obtaining the neighborhood conflict coefficient of the target prefabricated component; establishing a bell-shaped curve according to the relative displacement between any adjacent prefabricated components, and using the bell-shaped curve to obtain the connection area between the adjacent prefabricated components.

[0057] It should be noted that in the assembled retaining wall, due to the different connection modes between the prefabricated components, the relative displacement between the prefabricated components under the influence of load conditions is usually restricted by the connection mode, so that even if a large degree of relative displacement occurs, it may not cause a large degree of deformation, and therefore when the finite element modeling of the area corresponding to the connection relationship is performed, high-resolution finite element units are not needed to form the corresponding model, so as to avoid wasting computing resources when using the finite element model of the assembled retaining wall for subsequent finite element analysis.

[0058] Specifically, in step S301, the neighborhood conflict coefficient of the target prefabricated component is obtained based on the difference between the displacement category value of the target prefabricated component and the displacement category value of the adjacent prefabricated component.

[0059] As an optional embodiment, the specific calculation method of the neighborhood conflict coefficient of the target prefabricated component is that the absolute value of the difference between the displacement category value of the target prefabricated component and the displacement category value of any adjacent prefabricated component is denoted as the neighborhood conflict factor between the target prefabricated component and the adjacent prefabricated component, and the average value of the neighborhood conflict factors between the target prefabricated component and all adjacent prefabricated components is taken as the neighborhood conflict coefficient of the target prefabricated component.

[0060] It should be noted that in the complex terrain and the substation scene coupled with multiple loads, some connection nodes (such as the intersection between the buttress column and the bottom plate) may be weak in structure or have significant interface slip, resulting in significant difference in displacement mode of the components on both sides; while other areas (such as the whole section of the bottom plate) have strong connection or good stiffness matching, and the displacement is highly synchronized; the neighborhood conflict coefficient is a spatial heterogeneity index that captures the local deformation incoordination. Incoordination of displacement between components will cause stress redistribution in the connection area, especially at the keyway, embedded part or contact surface, so the higher the conflict coefficient, the greater the risk; if the overall structure is composed of a large number of high-conflict areas, it indicates that the integrity of the assembled system is poor, and it is closer to a loose assembly; otherwise, it reflects good assembly coordination.

[0061] In step S302, a bell-shaped curve is established according to the relative displacement between any adjacent prefabricated components, and the connection area between the adjacent prefabricated components is obtained by using the bell-shaped curve.

[0062] As a preferred embodiment, the specific acquisition method of the connection area between the adjacent prefabricated components is: First, select the adjacent prefabricated components and the prefabricated component , respectively, to obtain the relative displacement data of the prefabricated component and the prefabricated component under all load conditions; along the prefabricated component And prefabricated components The normal direction of the connecting interface between the prefabricated components, divides the connecting interface into several equidistant measuring line segments, on each measuring line segment, uniformly arranges monitoring points, obtains the relative displacement values of the monitoring points under different load conditions; for each monitoring point, calculates the average relative displacement value of the monitoring point under all load conditions, wherein is a preset monitoring point quantity parameter.

[0063] It should be noted that the monitoring point quantity parameter is empirically preset to 5, which can be adjusted according to actual conditions, and the embodiments of the application are not specifically limited.

[0064] Then, all the monitoring points are sorted according to the position coordinates on the connecting interface; a bell-shaped curve in Gaussian distribution form is fitted with the position coordinates as the independent variable and the average relative displacement value as the dependent variable.

[0065] Finally, the effective range of the bell-shaped curve is determined, and the boundary position of the connecting region is defined as , wherein is a preset boundary factor, wherein is the mathematical expectation of the bell-shaped curve; represents the standard deviation of the bell-shaped curve; the region between the two boundary positions is regarded as the boundary region.

[0066] It should be noted that the value of the boundary factor is empirically preset to 2.5, which can be adjusted according to actual conditions, and the embodiments of the application are not specifically limited.

[0067] Step S303, taking any connecting region as a target connecting region, obtaining two prefabricated components constituting the connecting region, obtaining the information entropy of the domain conflict coefficient between the two prefabricated components corresponding to all connecting regions formed by the two prefabricated components, obtaining the boundary coefficient of the target connecting region through the domain conflict coefficient of the target connecting region and combining the information entropy.

[0068] As an optional embodiment, the specific calculation method of the boundary coefficient of the target connecting region is: , wherein is the boundary coefficient of the target connecting region; is the domain conflict coefficient of the target connecting region; is the information entropy of the domain conflict coefficient between the two prefabricated components corresponding to all connecting regions formed by the two prefabricated components; represents a sigmoid normalization function.

[0069] It should be noted that the boundary coefficient reflects the difference degree of the mechanical behavior between adjacent prefabricated components, the greater the value, the easier the connection area is to cause relative slip or stress concentration under the action of load, and finer mesh division needs to be used in finite element modeling; on the contrary, the smaller the boundary coefficient, the more coordinated the mechanical behavior of the connection area is, and relatively rough mesh division can be used. Through the scientific determination of the boundary coefficient, the adaptive division of the mesh in the connection area of the assembled retaining wall finite element model can be realized, which significantly improves the calculation efficiency while ensuring the calculation accuracy.

[0070] At this point, the boundary coefficient of each connection area in the assembled retaining wall is obtained by the above method.

[0071] Step S004: adjusting the finite element size of the corresponding connection area between the prefabricated components in the finite element modeling process by using the boundary coefficient, thereby obtaining the finite element model of the assembled retaining wall.

[0072] It should be noted that based on the obtained connection area boundary coefficient, the adaptive mesh division technology is used to implement differentiated mesh density control for different connection areas of the assembled retaining wall, thereby effectively reducing the calculation resource consumption while ensuring the calculation accuracy, improving the finite element modeling efficiency and the accuracy of the analysis results.

[0073] First, the basic finite element mesh size parameter is preset, and the finite element mesh size in the boundary area is obtained in combination with the bell-shaped curve of the boundary area and the boundary coefficient.

[0074] It should be noted that in the embodiment of the application, the value range of the basic finite element mesh size parameter is 10-50mm, and the specific value is determined according to the overall size of the retaining wall and the calculation accuracy requirement.

[0075] It should be noted that in the embodiment of the application, the value range of the mesh adjustment coefficient is 0.5-2.0, and the preferred value is 1.2 according to experience, which can be adjusted according to the actual situation, and the embodiment of the application is not limited specifically.

[0076] Then, the three-dimensional point cloud data of the assembled retaining wall is modeled, and the self-defined mesh control function of the finite element software is used to implement mesh division in the assembled retaining wall model according to the finite element mesh size distribution, thereby obtaining the finite element model of the assembled retaining wall.

[0077] As an optional embodiment, the relationship between the boundary coefficient and the finite element mesh size is defined as: wherein, represents the finite element mesh size at the corresponding position of the data point in the boundary area of the three-dimensional point cloud data; represents the finite element mesh size at the corresponding position of the data point in the boundary area of the three-dimensional point cloud data; represents the preset basic finite element mesh size parameter; This is the preset grid adjustment coefficient; This is the boundary coefficient of the connected region, with a value range of [0,1]. Represents the natural constant; Represents data points within the boundary region of 3D point cloud data. The corresponding value in the bell curve.

[0078] In addition, during the analysis using the finite element model, the finite element analysis results are compared and verified with the collected monitoring data. If the displacement error exceeds 15% or the stress error exceeds 20%, the boundary coefficients of the corresponding areas are fine-tuned according to the error distribution until the model accuracy meets the engineering analysis requirements.

[0079] This concludes the embodiment.

[0080] It should be noted that the embodiments used in this example The model is only used to represent negative correlations and the results of the constraint model output are in Within this range, in specific implementations, other models with the same purpose can be substituted; this embodiment is merely an example. The description will be based on a model, without making specific limitations on it. This refers to the input of the model.

[0081] Please see Figure 2 The diagram shows a structural block diagram of a finite element modeling system for a prefabricated retaining wall in a substation according to an embodiment of the present invention. The system includes a memory 202, a processor 201, and a computer program 2021 stored in the memory 202 and executable on the processor. When the processor 201 executes the computer program 2021, it implements steps S001 to S004 of the finite element modeling method for a prefabricated retaining wall in a substation.

[0082] Furthermore, in an optional embodiment, the memory 202 described above may include read-only memory and random access memory, and provide instructions and data to the processor. The memory 202 may also include non-volatile random access memory. For example, the memory may also store device type information.

[0083] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for finite element modeling of a substation fabricated retaining wall, characterized in that, The method comprises the following steps: An experimental environment is built, and the prefabricated retaining wall is assembled and monitored in real time by sensors to obtain displacement data corresponding to each prefabricated component in the prefabricated retaining wall; A plurality of load conditions are set, and the relative displacement characteristics of the prefabricated components are obtained by using the relative displacement of the displacement data of different prefabricated components in the prefabricated retaining wall under all load conditions; the prefabricated components are classified according to the relative displacement characteristics of the prefabricated components in the prefabricated retaining wall, and a displacement category value of each prefabricated component is obtained; Any prefabricated component is taken as a target prefabricated component, and the connection area between adjacent prefabricated components and the boundary coefficient of the connection area are determined based on the difference distribution between the displacement category values of the target prefabricated component and adjacent prefabricated components and the relative displacement between the adjacent prefabricated components; The finite element size of the connection area between the prefabricated components is adjusted by using the boundary coefficient in the process of finite element modeling, so as to obtain a finite element model of the prefabricated retaining wall.

2. The finite element modeling method of a substation assembled retaining wall according to claim 1, characterized in that, The specific method for setting a plurality of load conditions and obtaining the relative displacement characteristics of the prefabricated components by using the relative displacement of the displacement data of different prefabricated components in the prefabricated retaining wall under all load conditions comprises the following steps: A plurality of load conditions are set, and the displacement data of all prefabricated components of the prefabricated retaining wall under each load condition is obtained; Under any load condition, the relative displacement factor of the target prefabricated component under the load condition is calculated according to the relative displacement of the displacement data between the target prefabricated component and other prefabricated components connected thereto; The relative displacement characteristics of the target prefabricated component are obtained by the numerical distribution of the relative displacement factor of the target prefabricated component under all load conditions and the correlation degree with the load conditions.

3. The finite element modeling method of a substation assembled retaining wall according to claim 2, characterized in that, The specific method for setting a plurality of load conditions comprises the following steps: In the experimental environment, a hydraulic loading system is used to apply graded loads, and horizontal load, vertical load and dynamic load with increasing degrees are applied to the prefabricated retaining wall to obtain a plurality of load conditions.

4. The finite element modeling method of a substation assembled retaining wall according to claim 2, characterized in that, The specific method for obtaining the relative displacement factor of the target prefabricated component under the load condition comprises the following steps: Under any load condition, the range of the displacement data of any prefabricated component is obtained as the displacement amplitude of the prefabricated component, and the displacement amplitude difference between any prefabricated component connected to the target prefabricated component and the target prefabricated component is obtained; The interval formed by the maximum and minimum values of the displacement data of any prefabricated component is recorded as the displacement interval of the prefabricated component, and the intersection-union ratio of the displacement intervals corresponding to any prefabricated component connected to the target prefabricated component and the target prefabricated component is recorded as the displacement interval overlap rate; The Pearson correlation coefficient of the displacement data corresponding to any prefabricated component connected to the target prefabricated component and the target prefabricated component is obtained, and the relative displacement factor of the target prefabricated component under any load condition is calculated in combination with the displacement amplitude difference and the displacement interval overlap rate.

5. The finite element modeling method of a substation assembled retaining wall according to claim 2, wherein, The specific method for obtaining the relative displacement characteristics of the target prefabricated component comprises the following steps: According to the numerical distribution of the relative displacement factor of the target prefabricated component under all load conditions, the relative displacement coefficient of the target prefabricated component is calculated; according to the correlation degree of the relative displacement factor and the load condition, the displacement-load correlation degree is calculated; The array formed by the relative displacement coefficient and the displacement-load correlation degree is taken as the relative displacement feature group of the target prefabricated component.

6. The finite element modeling method of a substation assembled retaining wall according to claim 1, wherein, The specific method for classifying the prefabricated components by the relative displacement feature group of the prefabricated components in the assembled retaining wall and obtaining the displacement category value of each prefabricated component comprises: Taking the Euclidean distance between the relative displacement feature groups of the prefabricated components as the distance measurement method of the DBSCAN clustering algorithm, and using the DBSCAN clustering algorithm to cluster all the prefabricated components to obtain a plurality of clustering clusters; According to the relative displacement level of all prefabricated components in each clustering cluster, the displacement category value of the prefabricated components in each clustering cluster is obtained, and for any clustering cluster, the product of the relative displacement coefficient and the displacement-load correlation degree of any prefabricated component in the clustering cluster is taken as the displacement level coefficient of the prefabricated component; Taking the average displacement level coefficient of all prefabricated components in the clustering cluster as the displacement type parameter of the clustering cluster, and performing maximum and minimum normalization on the displacement type parameters of all clustering clusters, the normalized result is taken as the displacement category value of the prefabricated components in the corresponding clustering cluster.

7. The finite element modeling method of a substation assembled retaining wall according to claim 1, characterized in that, The specific method for determining the connection area between adjacent prefabricated components and the boundary coefficient of the connection area based on the difference distribution between the displacement category values of the target prefabricated component and the adjacent prefabricated components and combining the relative displacement between the adjacent prefabricated components comprises: Based on the difference between the displacement category values of the target prefabricated component and the adjacent prefabricated components, the neighborhood conflict coefficient of the target prefabricated component is obtained; According to the relative displacement between any adjacent prefabricated components, a bell-shaped curve is established, and the connection area between the adjacent prefabricated components is obtained by using the bell-shaped curve; Taking any connection area as a target connection area, two prefabricated components constituting the connection area are obtained, the information entropy of the neighborhood conflict coefficient corresponding to all connection areas formed by the two prefabricated components is obtained, and the boundary coefficient of the target connection area is obtained by combining the neighborhood conflict coefficient of the target connection area with the information entropy.

8. The finite element modeling method of a substation assembled retaining wall according to claim 7, characterized in that, The specific method for obtaining the neighborhood conflict coefficient of the target prefabricated component comprises: The absolute value of the difference between the displacement category values of the target prefabricated component and any adjacent prefabricated component is taken as the neighborhood conflict factor between the target prefabricated component and the adjacent prefabricated component, and the average value of the neighborhood conflict factors between the target prefabricated component and all adjacent prefabricated components is taken as the neighborhood conflict coefficient of the target prefabricated component.

9. The finite element modeling method of a substation assembled retaining wall according to claim 7, characterized in that, The specific method for establishing a bell-shaped curve according to the relative displacement between any adjacent prefabricated components and obtaining the connection area between the adjacent prefabricated components by using the bell-shaped curve comprises: selecting adjacent precast components and precast components , respectively obtaining relative displacement data of the precast components and precast components under all load conditions; dividing the connecting interface between the precast components and precast components into a plurality of equally-spaced measuring line segments along the normal direction of the connecting interface, and arranging monitoring points on each measuring line segment to obtain relative displacement values of the monitoring points under different load conditions; calculating, for each monitoring point, an average relative displacement value of the monitoring point under all load conditions, wherein is a preset monitoring point quantity parameter; All monitoring points are sorted according to the position coordinates on the connection interface; A bell-shaped curve in Gaussian distribution form is fitted with the position coordinates as independent variables and the average relative displacement values as dependent variables. determining the effective range of the bell curve, defining the boundary position of the connection area as wherein, is a preset boundary factor, wherein is the mathematical expectation of the bell curve; represents the standard deviation of the bell curve; the area between the boundary positions on both sides is taken as the boundary area.

10. The finite element modeling method of a substation assembled retaining wall according to claim 9, wherein, The method comprises the following steps of: presetting a finite element grid size parameter of a foundation, and obtaining a finite element grid size in a boundary region by combining a bell-shaped curve of the boundary region and a boundary coefficient; performing three-dimensional modeling on three-dimensional point cloud data of the assembled retaining wall, and implementing grid division in the assembled retaining wall model according to finite element grid size distribution by using a self-defined grid control function of finite element software, so as to obtain a finite element model of the assembled retaining wall.

Images