Calculation method of shared tower load based on finite element analysis

Abstract

The present application relates to the field of electrical data processing technology, and specifically to a method for calculating a shared tower load based on finite element analysis, the method comprising: establishing a finite element model of a tower structure by finite element analysis, obtaining the displacement vector and stress data of each node on the tower at each moment, and the height of each node based on the simulation process of the finite element model; determining the force consistency coefficient of each node at each moment; obtaining the relative difference coefficient of any node at each moment; calculating the load transfer anomaly of any node at each moment; determining the abnormal deviation coefficient of any node; determining the transfer anomaly factor of any node; determining the nonlinear effect influence index of the tower; obtaining the adjusted wind load adjustment coefficient, and calculating the load of the tower. The present application improves the accuracy of the shared tower load calculation so as to more accurately evaluate the stability and safety of the tower.

Description

Technical Field

[0001] The present application relates to the technical field of electrical data processing, and in particular to a method for calculating a shared tower load based on finite element analysis. Background Art

[0002] Due to the rapid development of communication technology, the demand for communication towers has greatly increased. It is necessary to rationally utilize the existing power tower resources and add communication base stations to the transmission towers to realize the sharing of power transmission and communication. The shared tower model came into being. The wind, ice, and the deadweight load of communication equipment affect the structural strength and stability of the tower. How to accurately calculate its load under various complex working conditions has become a key issue.

[0003] The development of finite element analysis technology has provided a new way to solve these problems. Finite element analysis is a numerical analysis technology that uses mathematical approximation to simulate real physical systems. Due to the strong geometric nonlinearity of the transmission tower line structure and the complexity of wind, ice and the deadweight load of communication equipment, the traditional finite element analysis process fails to fully consider the impact of multiple load conditions on the nonlinear effect of the tower structure, making it difficult to accurately evaluate the bearing performance of the transmission tower under different stress states, resulting in low accuracy in the calculation of shared tower loads. Summary of the invention

[0004] In order to solve the above technical problems, a shared tower load calculation method based on finite element analysis is provided to solve the existing problems.

[0005] The solution to the technical problem of this application is to provide a shared tower load calculation method based on finite element analysis, including the following steps:

[0006] The finite element model of the tower structure is established by the finite element analysis method. Based on the simulation process of the finite element model, the displacement vector and stress data of each node on the tower at each time and the height of each node are obtained;

[0007] Analyze the correlation between the displacement vectors of each node and its adjacent nodes at each moment, and determine the force consistency coefficient of each node at each moment in combination with the difference in stress data;

[0008] Clustering the force consistency coefficients of all nodes at each moment; analyzing the difference between the force consistency coefficients of any node at each moment and the remaining nodes in the cluster to which it belongs, as well as the difference between the force consistency coefficients of any node and the nodes in other clusters, to obtain the relative difference coefficient of any node at each moment; combining the discrete distribution of the corresponding heights of the nodes in the cluster to which any node belongs at each moment, calculating the load transfer abnormality of any node at each moment;

[0009] According to the difference degree of the load transfer anomaly between any node and the remaining nodes at all times, and the concentrated distribution trend of the load transfer anomaly, the abnormal deviation coefficient of any node is determined; based on the change trend of the displacement vector of any node at all times, combined with the abnormal deviation coefficient, the transfer anomaly factor of any node is determined; the nonlinear effect influence index of the tower is determined by fusing the transfer anomaly factors of all nodes; based on the nonlinear effect influence index, the adjusted wind load adjustment coefficient is obtained, and the load of the tower is calculated.

[0010] Preferably, the correlation degree is the cosine similarity between the displacement vector of each node at each moment and the displacement vectors of each of its adjacent nodes at the same moment.

[0011] Preferably, determining the force consistency coefficient of each node at each moment includes:

[0012] Calculate the difference between the stress data of each node at each moment and the stress data of each adjacent node at the same moment, and record it as stress difference;

[0013] The sum of the correlation degree and the value 1 is calculated; the ratio of the sum to the stress difference is calculated, and the average of the ratios of each node to all other adjacent nodes at each moment is taken as the force consistency coefficient of each node at each moment.

[0014] Preferably, obtaining the relative difference coefficient of any node at each moment includes:

[0015] The mean of the difference between the force consistency coefficient of any node in each cluster corresponding to each moment and all other nodes is recorded as the intra-cluster difference corresponding to any node at each moment;

[0016] The mean of the difference in the force consistency coefficient between any node in each cluster corresponding to each moment and all nodes in all other clusters is recorded as the inter-cluster difference corresponding to any node at each moment;

[0017] The ratio of the inter-cluster difference to the intra-cluster difference corresponding to any node at each moment is calculated as the relative difference coefficient of any node at each moment.

[0018] Preferably, the calculation method of the load transfer abnormality of any node at each moment is:

[0019] Calculate the discrete degree of the height of all nodes in each cluster corresponding to each moment;

[0020] The product of the relative difference coefficient of any node at each moment and the discrete degree of the cluster to which the any node belongs is taken as the load transfer abnormality of the any node at each moment.

[0021] Preferably, the difference degree is the DTW distance of the load transfer anomaly between any node and the remaining nodes at all times.

[0022] Preferably, determining the abnormal deviation coefficient of any node includes:

[0023] The average value of the load transfer anomaly degree at all times corresponding to any node and the remaining nodes is recorded as the average anomaly coefficient;

[0024] The product of the difference degree and the average abnormal coefficient is recorded as the load deviation; the sum of the load deviations of any node and the remaining nodes is taken as the abnormal deviation coefficient of any node.

[0025] Preferably, the determining of the transmission anomaly factor of any node includes:

[0026] Calculating the modulus of the displacement vector of any node at each time; calculating the autocorrelation coefficient of the modulus at all times of any node;

[0027] The ratio of the abnormal deviation coefficient of any node to the absolute value of the autocorrelation coefficient is calculated as the transmission abnormality factor of any node.

[0028] Preferably, the nonlinear effect influence index of the tower is the cumulative sum of the transmission anomaly factors of all nodes.

[0029] Preferably, the calculation method of the adjusted wind load adjustment coefficient is: ,in, The height of the tower is The adjusted wind load adjustment factor is The height of the tower is The preset initial wind load adjustment factor is is the nonlinear effect index of the tower, is the normalization function.

[0030] This application has at least the following beneficial effects:

[0031] The present application analyzes the correlation between the displacement vectors of each node and its adjacent nodes at each moment, and determines the force consistency coefficient of each node at each moment in combination with the difference in stress data. The beneficial effect is that the force consistency of adjacent nodes is taken into account to reflect the stability of the tower structure; the force consistency coefficients of all nodes at each moment are clustered; the difference in the force consistency coefficients of any node at each moment and the remaining nodes in the cluster to which it belongs, as well as the difference in the force consistency coefficients between any node and the nodes in other clusters are analyzed to obtain the relative difference coefficient of any node at each moment. The beneficial effect is that the deviation in the force consistency of nodes between different clusters is taken into account; combined with the discrete distribution of the corresponding heights of the nodes in the cluster to which any node belongs at each moment, the load transfer anomaly of any node at each moment is calculated. The beneficial effect is that the distribution of nodes with small force differences on the tower is taken into account, thereby reflecting the force conditions of nodes at different positions on the tower to illustrate the possibility of abnormal load transfer between different nodes; according to the difference in the load transfer anomaly between any node and the remaining nodes at all times The abnormal deviation coefficient of any node is determined based on the change trend of the displacement vector of any node at all times and the concentrated distribution trend of the load transfer abnormality; the transfer abnormal factor of any node is determined based on the change trend of the displacement vector of any node at all times in combination with the abnormal deviation coefficient; the transfer abnormal factors of all nodes are integrated to determine the nonlinear effect influence index of the tower, which has the beneficial effect of considering the change trend of the displacement of different nodes at different times to reflect the regularity of the displacement state change of the node under the influence of wind force. Secondly, the difference in the force of different nodes is analyzed to reflect the load transfer state of different nodes, indicating whether the corresponding nodes are under the same load conditions in the structure, indicating the influence of the nonlinear effect on the shared tower, and reflecting the possible degree of instability of the tower under the action of wind force; based on the nonlinear effect influence index, the adjusted wind load adjustment coefficient is obtained to calculate the load of the tower, which has the beneficial effect of considering the influence of the nonlinear effect on the tower, resulting in the uncertainty and complexity of the wind load on the tower structure, thereby selecting a suitable wind load adjustment coefficient to improve the accuracy of the shared tower load calculation, so as to more accurately evaluate the stability and safety of the tower. BRIEF DESCRIPTION OF THE DRAWINGS

[0032] The shared tower load calculation method based on finite element analysis of the present application is further described in detail below in conjunction with the accompanying drawings.

[0033] Figure 1 A flow chart of the steps of a shared tower load calculation method based on finite element analysis provided in an embodiment of the present application;

[0034] Figure 2A flowchart of the steps of a method for obtaining the load transfer anomaly degree of any node at each moment provided in an embodiment of the present application;

[0035] Figure 3 A flowchart of the steps of a method for obtaining a nonlinear effect influence index of an iron tower provided in an embodiment of the present application. DETAILED DESCRIPTION

[0036] In order to make the purpose, technical solution and advantages of this application more clear, the following is a further detailed description of the shared tower load calculation method based on finite element analysis proposed in this application in conjunction with the accompanying drawings and implementation examples. It should be understood that the specific embodiments described herein are only used to explain this application and are not used to limit this application.

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

[0038] See also Figure 1 , which shows a flow chart of the steps of a shared tower load calculation method based on finite element analysis provided by an embodiment of the present application, the method comprising the following steps:

[0039] Step 1: Establish a finite element model of the tower structure through the finite element analysis method. Based on the simulation process of the finite element model, obtain the displacement vector and stress data of each node on the tower at each time, as well as the height of each node.

[0040] Finite Element Analysis (FEA) is a numerical calculation method that approximates the solution of engineering and physical problems by discretizing a continuous structure or system into a finite number of small units (called finite elements). Finite element analysis discretizes a continuous problem into a finite number of unit combinations that are only connected to each other at a finite number of nodes. This technology can simulate the performance of the tower in a complex environment. Therefore, in the process of analyzing the bearing performance of the shared tower under the influence of different working loads, the finite element model of the shared tower is established by using ANSYS software based on the structural dimensions, material properties and other parameters of the shared tower. Based on the constructed finite element model, each structural unit in the main material and cross diagonal material of the shared tower can be obtained. The structural units are connected through support nodes, so each node on the tower is obtained.

[0041] In order to analyze whether the tower is unstable under different working conditions, the force conditions of the nodes in the tower are simulated by applying wind load, ice load and communication equipment deadweight load. Therefore, in the process of applying load to the finite element model, the displacement vector and stress data of each node of the tower in the finite element model at different times are collected. The collection time interval is t, and the height of each node in the finite element model is collected.

[0042] Preferably, in this embodiment, the collection time interval is 0.1s. As other implementation methods, the implementer can set it according to the actual situation.

[0043] At this point, the displacement vector and stress data of each node of the tower at each moment, as well as the height of each node, are obtained.

[0044] Step 2: Analyze the correlation between the displacement vectors of each node and its adjacent nodes at each moment, and determine the force consistency coefficient of each node at each moment in combination with the difference in stress data.

[0045] Shared towers can be divided into two parts, transmission towers and communication equipment. The load of shared towers can be regarded as the load of transmission towers plus the load of communication equipment. The loads on shared towers in actual projects are relatively complex, and the loads they are subjected to mainly include wind loads, ice loads, and wire loads. In the load calculation of shared towers, ice loads and the deadweight loads of communication equipment are mainly determined by analyzing the pressure they exert on the nodes. In contrast, the calculation of wind loads is more complicated because it needs to consider the correlation characteristics between the spatial structures of the towers, and the ice layer covered in the tower structure and the communication equipment will increase the area of ​​action of the wind force. In addition, wind loads are also the main factor causing the instability of the tower structure. Therefore, the wind loads of shared towers are calculated by analyzing the relevant characteristics of the stress and displacement states of the nodes.

[0046] Normally, the magnitude or direction of wind acting on the tower changes rapidly, causing frequent eddies near the tower. In addition, due to the influence of the geometry, material or boundary conditions of the structural unit, it is easy to cause nonlinear effects in the tower structure, thus affecting the accuracy of load calculation. The force difference characteristics between different nodes can accurately reflect the degree of influence of the tower by nonlinear effects.

[0047] Based on the above analysis, the nodes are connected through structural units, so that the stress data and displacement between adjacent nodes can affect each other. When the tower structure is stable, the displacement, displacement direction and stress magnitude between adjacent nodes should be relatively consistent at the same time. Therefore, by analyzing the force conditions between adjacent nodes, the force consistency coefficient is determined, which is:

[0048] Calculate the cosine similarity between the displacement vector of each node at each time and the displacement vector of each of its adjacent nodes at the same time;

[0049] It should be noted that the calculation of cosine similarity is a well-known technique and will not be described in detail here.

[0050] Calculate the difference between the stress data of each node at each moment and the stress data of each adjacent node at the same moment, and record it as stress difference;

[0051] Preferably, in this embodiment, the absolute value of the difference between the stress data of each node at each moment and the stress data of each adjacent node at the same moment is calculated and recorded as the stress difference.

[0052] Calculate the sum of the cosine similarity and the value 1;

[0053] Calculate the ratio of the sum to the stress difference, and take the average of the ratios of each node to all other adjacent nodes at each moment as the force consistency coefficient of each node at each moment;

[0054] Preferably, in this embodiment, the calculation method of the force consistency coefficient of each node at each moment is: ,in, For the The node in The force consistency coefficient at the moment, For the The node in The displacement vector at time, For the The node adjacent to The node in The displacement vector at time, For the The node in The stress data at the moment, For the The node adjacent to The node in The stress data at the moment, For the The number of other nodes adjacent to a node, To calculate the cosine similarity, To preset a value greater than 0, avoid the denominator being 0. The value range is , in this embodiment, The value is 0.1. As other implementation methods, the implementer can set it according to the actual situation. is the stress difference.

[0055] It should be noted that cosine similarity The closer it is to 1, the more consistent the displacement distance and direction of the two nodes. Adding 1 to the numerator is to avoid the numerator being a negative number, so the sum is The value of , the larger the sum is, the closer the displacement direction and magnitude between the two nodes are; secondly, the stress difference The smaller it is, the smaller the stress difference between the two nodes is, and the larger the force consistency coefficient is, which means that the stress state of the corresponding node is more consistent with that of its adjacent nodes, reflecting that the structure of the tower is more stable.

[0056] At this point, the force consistency coefficient of each node at each moment is obtained.

[0057] Step 3, clustering the force consistency coefficients of all nodes at each moment; analyzing the difference in the force consistency coefficients of any node at each moment and the remaining nodes in the cluster to which it belongs, as well as the difference in the force consistency coefficients between any node and the nodes in other clusters, to obtain the relative difference coefficient of any node at each moment; combining the discrete distribution of the corresponding heights of the nodes in the cluster to which any node belongs at each moment, calculate the load transfer anomaly of any node at each moment.

[0058] Furthermore, the displacement degree and stress magnitude of the nodes of the shared tower are different at different heights, but the difference in the force consistency relationship corresponding to different nodes is small. For example, the displacement of the nodes corresponding to the bottom and top positions of the tower is quite different. However, the load between these two nodes is gradually transmitted through multiple nodes in the middle. Therefore, when the influence of nonlinear effects is relatively small, the difference in the force consistency coefficients corresponding to different nodes in the entire tower should be small, and the load transfer state in the tower is relatively good.

[0059] However, different parts of the tower will use different metal materials. Due to the different molecular structures of the metal materials themselves, especially when affected by nonlinear effects, the force consistency coefficients corresponding to the nodes in different parts of the tower will have different degrees of deviation. At this time, the load transfer state inside the tower is relatively poor, which makes the force consistency coefficients of all nodes show different degrees of category difference characteristics.

[0060] Furthermore, the flowchart of the method for obtaining the load transfer abnormality of any node at each time provided in the embodiment of the present application is as follows: Figure 2 shown.

[0061] First, cluster analysis is performed on the force consistency coefficients between different nodes, specifically:

[0062] Performing cluster analysis on the force consistency coefficients of all nodes at each moment to obtain multiple clusters;

[0063] Preferably, in this embodiment, an AP (Affinity propagation) clustering algorithm is used to perform clustering analysis, wherein the AP clustering algorithm is a well-known technology and will not be described in detail herein.

[0064] It should be noted that the greater the number of clusters, the greater the difference in force consistency coefficients between different nodes.

[0065] Furthermore, the differences in the force consistency coefficients between the nodes in each cluster and the height distribution of the nodes are analyzed to determine the load transfer state to reflect the load transfer state inside the tower of the node, specifically:

[0066] The mean of the difference between the force consistency coefficient of any node in each cluster and all other nodes at each moment is recorded as the intra-cluster difference corresponding to any node at each moment;

[0067] The mean of the difference in the force consistency coefficient between any node in each cluster corresponding to each moment and all nodes in all other clusters is recorded as the inter-cluster difference corresponding to any node at each moment;

[0068] Preferably, in this embodiment, the average of the absolute values ​​of the difference in the force consistency coefficient between any node in each cluster and all other nodes corresponding to each moment is recorded as the intra-cluster difference corresponding to any node at each moment; the average of the absolute values ​​of the difference in the force consistency coefficient between any node in each cluster and all other nodes in all clusters corresponding to each moment is recorded as the inter-cluster difference corresponding to any node at each moment.

[0069] Calculating the ratio of the inter-cluster difference to the intra-cluster difference corresponding to any node at each moment as the relative difference coefficient of any node at each moment;

[0070] Calculate the discrete degree of the height of all nodes in each cluster corresponding to each moment;

[0071] Preferably, in this embodiment, the degree of dispersion is measured by calculating the variance of the heights of all nodes in each cluster corresponding to each moment.

[0072] The product of the relative difference coefficient of any node at each moment and the discrete degree of the cluster to which the any node belongs is used as the load transfer abnormality of any node at each moment;

[0073] It should be noted that the smaller the intra-cluster difference is, the smaller the force consistency difference between the corresponding node and the other nodes in the same cluster is; the larger the inter-cluster difference is, the greater the force consistency difference between the corresponding node and the nodes in other clusters is, and the larger the relative difference coefficient is, the greater the force difference between the corresponding node and other nodes is; the larger the degree of discreteness is, the more discrete the heights of the nodes in the corresponding cluster are, and the more likely the tower is to become unstable; the larger the load transfer anomaly is, the greater the difference in force consistency coefficients between different nodes is, reflecting that the possibility of load transfer anomaly at the corresponding node on the tower is greater, that is, the load transfer state inside the tower is relatively poor.

[0074] At this point, the load transfer anomaly degree of any node at each moment is obtained.

[0075] Step 4, determining the abnormal deviation coefficient of any node according to the difference degree of the load transfer abnormality between any node and the remaining nodes at all times, and the concentrated distribution trend of the load transfer abnormality; determining the transfer abnormality factor of any node based on the change trend of the displacement vector of any node at all times and the abnormal deviation coefficient; fusing the transfer abnormality factors of all nodes to determine the nonlinear effect influence index of the tower; based on the nonlinear effect influence index, obtaining the adjusted wind load adjustment coefficient, and calculating the load of the tower.

[0076] Furthermore, when the iron tower is subjected to wind load, each node will more or less have a certain degree of displacement. For an iron tower in a good stable state, the displacement of each node at different times is not uniform, but shows a regular ups and downs state, and the load transfer anomaly corresponding to different nodes is small, and the load transfer state between nodes is relatively close. Therefore, the displacement changes of nodes at different times and the differences in load transfer anomaly between different nodes are analyzed to determine the nonlinear effect influence index to reflect the possible degree of instability of the iron tower under the action of wind. The step flow chart of the method for obtaining the nonlinear effect influence index of the iron tower provided in the embodiment of the present application is as follows: Figure 3 As shown, including:

[0077] Calculating the DTW distance of the load transfer anomaly degree between any node and the remaining nodes at all times;

[0078] Calculating the modulus of the displacement vector of any node at each time; calculating the autocorrelation coefficient of the modulus at all times of any node;

[0079] It should be noted that the calculation of DTW distance, module length and autocorrelation coefficient are all well-known technologies and will not be described in detail here. The time delay of the autocorrelation coefficient is 3. As other implementation methods, the implementer can set it according to the actual situation.

[0080] The average value of the load transfer anomaly degree at all times corresponding to any node and the remaining nodes is recorded as the average anomaly coefficient;

[0081] The product of the DTW distance and the average anomaly coefficient is recorded as the load deviation; the sum of the load deviations of any node and the rest of the nodes is recorded as the anomaly deviation coefficient of any node;

[0082] Calculating a ratio of the abnormal deviation coefficient of any node to the absolute value of the autocorrelation coefficient as a transmission abnormality factor of any node;

[0083] The sum of the transmission anomaly factors of all nodes is used as the nonlinear effect influence index of the tower;

[0084] Preferably, in this embodiment, the calculation method of the nonlinear effect influence index of the iron tower is: ,in, is the nonlinear effect index of the tower, For the The node and The average abnormal coefficient of nodes, No. The node and The DTW distance of the load transfer anomaly between nodes at all times, For the The autocorrelation coefficient of the nodes, is the number of all nodes on the tower; secondly, is the abnormal deviation coefficient, To transmit abnormal factors.

[0085] It should be noted that, the larger the absolute value of the autocorrelation coefficient is, the more regular the displacement state of the corresponding node is under the influence of wind load, the smaller the DTW distance is, the closer the load transfer state between different nodes is, the smaller the average anomaly coefficient is, the smaller the possibility of load transfer anomaly between different nodes is, the smaller the obtained transfer anomaly factor is, the more consistent the load transfer state between the corresponding node and other nodes is; conversely, the larger the obtained transfer anomaly factor is, the more inconsistent the load transfer state between the corresponding node and other nodes is, reflecting that the corresponding node is subjected to unique load conditions in the structure; secondly, the nonlinear effect influence index reflects the possible degree of instability of the tower under the action of wind force, the smaller the obtained nonlinear effect influence index is, the less the shared tower is affected by the nonlinear effect, and the smaller the possibility of instability of the tower is, conversely, the larger the nonlinear effect influence index is, the more affected the tower is by the nonlinear effect, and the greater the possibility of instability of the tower is.

[0086] Furthermore, when the tower is more affected by the nonlinear effect, it means that the behavior of the tower structure in the force process is no longer a simple linear relationship. The nonlinear structure may produce a greater dynamic response under the action of wind load, which may increase the complexity of the structural response and the possibility of instability of the tower. When the wind load of the tower is calculated by the wind load calculation formula, the wind load adjustment coefficient is a coefficient used to correct the wind load. It reflects the uncertainty and complexity of the wind load. When the nonlinear effect has a greater impact, this uncertainty and complexity increases. By selecting a larger wind load adjustment coefficient, the impact of these nonlinear factors on the structural response is taken into account, thereby more accurately evaluating the performance and safety of the tower under extreme wind loads; conversely, the smaller the nonlinear effect influence index, the smaller the impact of the nonlinear effect on the tower, and a relatively small wind load adjustment coefficient can be selected.

[0087] Based on the above analysis and the nonlinear effect influence index, the wind load adjustment coefficient is adjusted as follows:

[0088] The calculation method of the adjusted wind load adjustment factor is: ,in, The height of the tower is The adjusted wind load adjustment factor is The height of the tower is The preset initial wind load adjustment factor is is the nonlinear effect index of the tower, is a normalization function. In this embodiment, the tanh function is used for normalization processing, wherein the tanh function is a well-known technology and will not be described in detail here.

[0089] Preferably, in this embodiment, the preset initial wind load adjustment coefficient is set to 1.5, wherein the value of the preset initial wind load adjustment coefficient is related to the tower height and the type of tower. The initial wind load adjustment coefficient is selected in accordance with the "Overhead Transmission Line Load Regulations". As other implementation methods, the implementer can set it according to actual conditions.

[0090] Based on the adjusted wind load adjustment coefficient and according to the wind load calculation formula, the wind load of the shared tower is determined.

[0091] It should be noted that the adjusted wind load calculation formula is: ,in, is the standard value of wind load perpendicular to the tower surface, ; is the standard value of the base wind pressure, ; is the wind pressure height variation coefficient, which is selected according to the "Overhead Transmission Line Load Regulations" and is determined according to the height of the tower from the ground or sea level and the roughness of the ground; is the body coefficient of the component, which is selected according to the "Overhead Transmission Line Load Regulations". For the tower composed of steel structure, the value of this coefficient is ,in is the wind load reduction factor on the leeward side of the tower, The value of is selected in accordance with the "Overhead Transmission Line Load Regulations". For towers with mixed structures composed of different types of cross-sections such as steel pipes and angle steels, it should be calculated separately according to the windward areas of different types of rods or determined by averaging the windward areas of the towers. is the ice-covered wind load increase coefficient of the tower components. The value of this coefficient is related to the geographical location. If the ice thickness is 5mm, it is 1.1; if the ice thickness is 10mm, it is 1.2; if the ice thickness is 15mm, it is 1.6; if the ice thickness is 20mm, it is 1.8; if the ice thickness is more than 20mm, it is 2.0-2.5; is the adjusted wind load adjustment factor; is the projected area of ​​the component subjected to wind pressure, .

[0092] It should be noted that the wind load calculation formula is a well-known technology and will not be repeated here. The calculation formula for wind load and the wind pressure height variation coefficient are detailed in the "Regulations on Overhead Transmission Line Loads". , the shape coefficient of the component selection regulations.

[0093] Based on the calculated wind load on the tower, other loads on the shared tower are calculated, including ice load and deadweight load of communication equipment. Based on all the loads on the shared tower, the shared tower can be stabilized and reinforced later.

[0094] It should be understood that although Figure 1 The steps in the flowchart are shown in sequence as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. Moreover, Figure 1 At least part of the steps may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least part of the sub-steps or stages of other steps.

[0095] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0096] The above-mentioned embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the present application. It should be pointed out that, for those of ordinary skill in the art, several modifications and improvements can be made without departing from the concept of the present application. Therefore, any simple modification, equivalent change and modification made to the above embodiments according to the technical essence of the present application without departing from the content of the technical solution of the present application, shall fall within the protection scope of the technical solution of the present application.

Claims

1. The shared tower load calculation method based on finite element analysis is characterized by: The method comprises the following steps: The finite element model of the tower structure is established by the finite element analysis method. Based on the simulation process of the finite element model, the displacement vector and stress data of each node on the tower at each time and the height of each node are obtained; Analyze the correlation between the displacement vectors of each node and its adjacent nodes at each moment, and determine the force consistency coefficient of each node at each moment in combination with the difference in stress data; Clustering the force consistency coefficients of all nodes at each moment; analyzing the difference between the force consistency coefficients of any node at each moment and the remaining nodes in the cluster to which it belongs, as well as the difference between the force consistency coefficients of any node and the nodes in other clusters, to obtain the relative difference coefficient of any node at each moment; combining the discrete distribution of the corresponding heights of the nodes in the cluster to which any node belongs at each moment, calculating the load transfer abnormality of any node at each moment; According to the difference degree of the load transfer anomaly between any node and the remaining nodes at all times, and the concentrated distribution trend of the load transfer anomaly, the abnormal deviation coefficient of any node is determined; based on the change trend of the displacement vector of any node at all times, combined with the abnormal deviation coefficient, the transfer anomaly factor of any node is determined; the nonlinear effect influence index of the tower is determined by fusing the transfer anomaly factors of all nodes; based on the nonlinear effect influence index, the adjusted wind load adjustment coefficient is obtained, and the load of the tower is calculated.

2. The shared tower load calculation method based on finite element analysis according to claim 1, characterized in that: The correlation degree is the cosine similarity between the displacement vector of each node at each time and the displacement vectors of each of its adjacent nodes at the same time.

3. The shared tower load calculation method based on finite element analysis according to claim 1, characterized in that: Determining the force consistency coefficient of each node at each moment includes: Calculate the difference between the stress data of each node at each moment and the stress data of each adjacent node at the same moment, and record it as stress difference; The sum of the correlation degree and the value 1 is calculated; the ratio of the sum to the stress difference is calculated, and the average of the ratios of each node to all other adjacent nodes at each moment is taken as the force consistency coefficient of each node at each moment.

4. The shared tower load calculation method based on finite element analysis according to claim 1, characterized in that: The obtaining of the relative difference coefficient of any node at each moment includes: The mean of the difference between the force consistency coefficient of any node in each cluster corresponding to each moment and all other nodes is recorded as the intra-cluster difference corresponding to any node at each moment; The mean of the difference in the force consistency coefficient between any node in each cluster corresponding to each moment and all nodes in all other clusters is recorded as the inter-cluster difference corresponding to any node at each moment; The ratio of the inter-cluster difference to the intra-cluster difference corresponding to any node at each moment is calculated as the relative difference coefficient of any node at each moment.

5. The shared tower load calculation method based on finite element analysis according to claim 1, characterized in that: The calculation method of the load transfer abnormality of any node at each time is: Calculate the discrete degree of the height of all nodes in each cluster corresponding to each moment; The product of the relative difference coefficient of any node at each moment and the discrete degree of the cluster to which the any node belongs is taken as the load transfer abnormality of the any node at each moment.

6. The shared tower load calculation method based on finite element analysis according to claim 1, characterized in that: The difference degree is the DTW distance of the load transfer anomaly between any node and the remaining nodes at all times.

7. The shared tower load calculation method based on finite element analysis according to claim 1, characterized in that: The determining of the abnormal deviation coefficient of any node includes: The average value of the load transfer anomaly degree at all times corresponding to any node and the remaining nodes is recorded as the average anomaly coefficient; The product of the difference degree and the average abnormal coefficient is recorded as the load deviation; the sum of the load deviations of any node and the remaining nodes is taken as the abnormal deviation coefficient of any node.

8. The shared tower load calculation method based on finite element analysis according to claim 1, characterized in that: The determining of the transmission anomaly factor of any node includes: Calculating the modulus of the displacement vector of any node at each time; calculating the autocorrelation coefficient of the modulus at all times of any node; The ratio of the abnormal deviation coefficient of any node to the absolute value of the autocorrelation coefficient is calculated as the transmission abnormality factor of any node.

9. The shared tower load calculation method based on finite element analysis according to claim 1, characterized in that: The nonlinear effect influence index of the tower is the cumulative sum of the transmission anomaly factors of all nodes.

10. The shared tower load calculation method based on finite element analysis according to claim 1, characterized in that: The calculation method of the adjusted wind load adjustment coefficient is: ,in, The height of the tower is The adjusted wind load adjustment factor is The height of the tower is The preset initial wind load adjustment factor is is the nonlinear effect index of the tower, is the normalization function.

Images