CN115310172A - Reconstruction method of transmission tower mechanics finite element model - Google Patents

Reconstruction method of transmission tower mechanics finite element model Download PDF

Info

Publication number
CN115310172A
CN115310172A CN202210827822.1A CN202210827822A CN115310172A CN 115310172 A CN115310172 A CN 115310172A CN 202210827822 A CN202210827822 A CN 202210827822A CN 115310172 A CN115310172 A CN 115310172A
Authority
CN
China
Prior art keywords
angle steel
coordinates
transmission tower
node
point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210827822.1A
Other languages
Chinese (zh)
Inventor
赵鹏飞
杜志叶
孟圣淳
龙振华
郝兆扬
岳国华
李�根
修连成
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN202210827822.1A priority Critical patent/CN115310172A/en
Publication of CN115310172A publication Critical patent/CN115310172A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method for reconstructing a transmission tower mechanics finite element model, which comprises the following steps: importing the three-dimensional solid model of the power transmission tower after the laser point cloud scanning processing into finite element software for identifying the angle steel coordinate information; screening out the outermost corner node coordinates in the key point coordinates of the end section of the angle steel as the end point coordinates of the angle steel unit; carrying out adjacent node fusion processing on the end point coordinates of all the angle steel units; identifying the condition that an auxiliary material and a main material are crossed in a transmission tower model, and adding segmented node coordinates on the main material at the crossed point; identifying the condition of crossing of auxiliary materials in the transmission tower model, and calculating node coordinates at the crossing points; and extracting the updated coordinates of the nodes at the two ends and the cross sectional area data of the end parts of the angle steels in the transmission tower, thereby completing the reconstruction of the whole transmission tower mechanical finite element model. The invention can perform equivalent node rigidity conversion on the connecting point of the angle steel rod piece, thereby greatly improving the modeling precision and the modeling efficiency.

Description

Reconstruction method of transmission tower mechanics finite element model
Technical Field
The invention belongs to the technical field of transmission tower digital twinning, and particularly relates to a reconstruction method of a transmission tower mechanics finite element model.
Background
The transmission tower is used as an infrastructure of electric power energy, and the safety and stability of the structure of the transmission tower are important guarantee for the safe operation of the power network. The coverage range of a transmission line in China is wide, and a tower line system is frequently influenced by natural disasters such as icing, strong wind and the like in a complex field environment for a long time, and accidents such as structural damage, tower falling and the like of a transmission tower can be caused under severe conditions, so that research works such as a damage mechanism, disaster prevention, disaster reduction and the like of a transmission tower line system mechanical model are actively developed at home and abroad to ensure the safe and stable operation of a power system, wherein the establishment of a real and accurate transmission tower mechanical finite element model is a basis and key for the simulation research of transmission line mechanical finite element, and the transmission tower mechanical model is a typical space truss structure model, so that at present, a three-dimension 1 is established according to actual structure drawing parameters of transmission tower construction facility man-hour: the method comprises the following steps of 1, establishing a refined transmission tower mechanical finite element model, wherein a large amount of data such as space coordinates of tower angle steel rods and corresponding size parameters need to be manually interpreted in the traditional method for establishing the three-dimensional refined transmission tower mechanical finite element model according to drawing structure parameters, the whole modeling process consumes a large amount of manual modeling time, the modeling efficiency is greatly reduced, and the transmission tower is easily influenced by geological disasters after construction and construction as time goes on, so that the problems of tower inclination, tower foundation slippage and the like can occur, and therefore certain model deviation exists between the transmission tower mechanical finite element model established according to the tower initial drawing parameters and a real transmission tower model under actual operating conditions, and the simulation calculation result of the mechanical finite element model can be inaccurate.
With the continuous development of laser point cloud scanning technology in recent years, the high-precision three-dimensional space information of a tower model can be obtained by applying the laser point cloud scanning technology to the aspect of power transmission tower scanning at present, and the three-dimensional entity model 1 of a real power transmission tower can be realized according to the extracted high-precision three-dimensional point cloud model: 1, reduction, but a mechanical model of a transmission tower is a typical high-flexibility space steel structure, has the characteristics of large deformation, remarkable geometric nonlinearity and the like, structural stress characteristics of various main materials, inclined materials and auxiliary materials in the tower need to be calculated when the transmission tower model is subjected to mechanical finite element simulation, and a transmission tower three-dimensional solid model generated by laser point cloud scanning does not perform equivalent node rigidity conversion on a connecting point of an angle steel rod piece, so that the transmission tower three-dimensional solid model cannot be directly used for mechanical finite element simulation analysis, so that the method for performing digital twin rapid reconstruction modeling on the mechanical finite element model based on the transmission tower solid model subjected to laser point cloud scanning is a key for solving the practical application problem of the transmission tower digital generative technology in the intelligent power grid.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides the reconstruction method of the power transmission tower mechanics finite element model, which can perform equivalent node rigidity conversion on the angle steel rod connecting points, and greatly improves the modeling precision and the modeling efficiency.
In order to solve the technical problems, the invention adopts the following technical scheme:
a method for reconstructing a power transmission tower mechanics finite element model comprises the following steps:
step 1, obtaining a three-dimensional solid model of the transmission tower through laser point cloud scanning, processing the three-dimensional solid model to obtain point cloud coordinate data of each angle steel in the model, identifying the total number of the angle steels in the solid model of the transmission tower and key point coordinate information of two ends of each angle steel, and calculating the cross sectional area of the end part of each angle steel;
step 2, identifying and screening out the key point coordinates of the outermost corner nodes in the end sections of the angle steel according to the key point coordinates of the cross sections of the two ends of each angle steel, and taking the key point coordinates as the node coordinates of the angle steel unit;
step 3, carrying out adjacent node fusion processing on the node coordinates of all the angle steel units;
step 4, identifying the crossing condition of the auxiliary material and the main material in the transmission tower solid model, adding node coordinates on the main material unit at the crossing point, segmenting the main material unit, modifying the node coordinates of the auxiliary material angle steel near the crossing point into the added node coordinates on the main material, and then performing adjacent node fusion processing again;
step 5, identifying the condition that the auxiliary material and the auxiliary material in the transmission tower solid model are crossed in an X mode, calculating node coordinates at the crossed points as newly-added node coordinates, and meanwhile segmenting crossed auxiliary material units;
and step 6, extracting the two-end node coordinates and the cross-sectional area data of the end parts of the angle steels of all the angle steel units in the transmission tower obtained in the step 1-5, and reconstructing a mechanical finite element model of the whole transmission tower according to the extracted node coordinates and the cross-sectional area data of the end parts of the angle steels.
Further, in step 2, a specific method for screening out the coordinates of the key points of the outermost corner nodes in the cross section of the end of the angle steel is as follows:
the method comprises the steps of firstly calculating the distance between any two key points according to key point coordinate data of the end part of the angle steel, then calculating whether the unit vector inner product of two lines connected with a certain key point is close to 0 to screen two edges which are perpendicular to each other, then comparing whether the lengths of the two edges are close to screen the outermost side and the innermost side of the end part of the angle steel, and finally selecting the lengths of the two edges connected with a node and the maximum node coordinate as the outermost side angle node coordinate.
Further, in step 3, when the adjacent node fusion processing is performed, the node coordinate with the large cross-sectional area of the angle steel is selected to be unchanged, and the node coordinate with the small cross-sectional area of the adjacent angle steel is modified, so that the coordinates of the node coordinates are finally fused into the node coordinate of the angle steel with the largest cross-sectional area.
Further, step 4 specifically includes:
step 4.1, extracting the cross-sectional area S of the first angle steel 1 And sequentially connecting the cross section area S of the angle steel with the cross section area S of the rear angle steel i Comparing two by two, when satisfying | S 1 -S i |≥ΔS 1 If so, one of the two angle steels is a main material and the other is an auxiliary material, wherein delta S is 1 Setting according to actual needs, determining an angle steel unit with a large cross-sectional area as a main material, calculating the distance between two end points of an auxiliary material and a main material line segment, and if the distance between one end point and the main material is delta L 1 If the content is within the range, the auxiliary material and the main material are crossed, wherein the content is Delta L 1 Setting according to actual needs;
4.2, on the premise that the auxiliary material and the main material are crossed, calculating the distance between the end points of the angle steel, and setting a judgment condition that the distance between the end point of the auxiliary material and the two end points of the main material is larger than a certain value so as to eliminate the condition that the main material and the auxiliary material are only crossed at the end points;
and 4.3, finally, calculating the coordinates of the newly added nodes on the main material at the cross points, and segmenting the main material according to the newly added nodes.
Further, the method for calculating the coordinates of the newly added nodes on the main material at the intersection in the step 4.3 comprises the following steps:
let the coordinate of the end point of the auxiliary material crossed with the main material be P (x) 0 ,y 0 ,z 0 ) The main line segment is AB, and the coordinate of the A endpoint is (x) 1 ,y 1 ,z 1 ) And the coordinates of the B endpoint are (x) 2 ,y 2 ,z 2 ) Drawing a perpendicular line of the line segment AB from the point P, and calculating a parameter t in a parameter equation of the straight line AB, wherein the perpendicular is Q (x, y, z) 1 Will t 1 The coordinates (X, Y, Z) of the foot drop point can be obtained as the coordinates of the newly added point by substituting the coordinates into a linear equation, wherein,
Figure BDA0003744673970000041
Figure BDA0003744673970000042
further, step 5 specifically includes:
step 5.1, extracting the cross-sectional area S of the first auxiliary material angle steel 1 And the cross sections are sequentially circulated with the cross section S of the auxiliary angle steel at the back i Comparing two by two, when satisfying | S 1 -S i |≤ΔS 2 Then, the cross sectional area of the two angle steels is close to each other, delta S 2 Setting according to actual needs, calculating the distance between every two end points of the angle steel, and if the distance between any two end points is smaller than the length of the angle steel, judging that the auxiliary material and the auxiliary material are crossed;
step 5.2, calculating node coordinates of the auxiliary material intersection points;
and 5.3, modifying the endpoint coordinates of the auxiliary material angle steel crossed at the moment, and segmenting the auxiliary material angle steel from the crossed nodes.
Further, the method for calculating the node coordinates of the auxiliary material intersection point comprises the following steps:
let a point P (X, Y, Z) exist on the first angle iron line segment AB, and the coordinate of the A endpoint is (X) 1 ,y 1 ,z 1 ) And the coordinates of the B endpoint are (x) 2 ,y 2 ,z 2 ) A point Q (U, V, W) exists on the second angle steel line segment CD, and the coordinate of the C endpoint is (x) 3 ,y 3 ,z 3 ) D end point coordinate is (x) 4 ,y 4 ,z 4 ) If the shortest distance between two angle steel line segments under the condition of auxiliary material intersection is PQ, calculating a parameter s in a linear equation 2 And t 2 Partial derivative f(s) of the square of the PQ shortest distance of 2 ,t 2 ) And let the partial derivative f(s) 2 ,t 2 ) Finding s for 0 2 And t 2 A parameter t 2 Substituting into linear equation to calculate node coordinate K (X) of cross point 1 ,Y 1 ,Z 1 ):
Figure BDA0003744673970000043
The distance between two points PQ is:
Figure BDA0003744673970000044
partial derivative f(s) of PQ shortest distance squared 2 ,t 2 ) Comprises the following steps:
Figure BDA0003744673970000051
node coordinate K (X) of the intersection 1 ,Y 1 ,Z 1 ):
Figure BDA0003744673970000052
Further, step 6 specifically includes:
extracting coordinates of end points and cross sectional area data of angle steels in the power transmission tower model obtained in the step 1-5, writing a parametric modeling command stream file required for establishing a mechanical finite element model by utilizing MATLAB, and rapidly completing modeling reconstruction of the power transmission tower mechanical finite element model in finite element software by loading the generated parametric modeling command stream file.
Compared with the prior art, the invention has the beneficial effects that:
1. according to the reconstruction method of the transmission tower mechanical finite element model, the transmission tower solid model is obtained through laser point cloud scanning, equivalent node rigidity conversion is carried out on the angle steel rod connecting points, and the transmission tower mechanical finite element model is rapidly reconstructed, so that the modeling precision and the modeling efficiency are greatly improved, and the obtained mechanical model can be directly used for mechanical finite element simulation analysis and calculation;
2. the method for reconstructing the transmission tower mechanical finite element model is applied to the technical field of power transmission and transformation digital twinning, can realize digital rapid reconstruction modeling from the transmission tower to the mechanical finite element model under a real operating condition, and provides an effective reference value for the practical application of the transmission tower digital twinning technology in an intelligent power grid;
3. the method for reconstructing the transmission tower mechanical finite element model is suitable for digital rapid reconstruction of transmission tower mechanical finite element models of various voltage classes.
Drawings
FIG. 1 is a flow chart of a reconstruction method of a transmission tower mechanics finite element model according to an embodiment of the invention;
FIG. 2 is a physical model of a power transmission tower after laser point cloud scanning processing in an embodiment of the invention;
FIG. 3 is an L-angle unit in the solid model according to the embodiment of the present invention;
FIG. 4 is a diagram illustrating a neighboring node fusion process according to an embodiment of the present invention;
fig. 5 is a power transmission tower mechanical finite element model with fast reconstruction in the embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be described clearly and completely with reference to the following embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without inventive efforts based on the embodiments of the present invention, shall fall within the scope of protection of the present invention.
It should be noted that the embodiments and features of the embodiments may be combined with each other without conflict.
The present invention is further illustrated by the following examples, which are not intended to limit the scope of the invention.
As shown in fig. 1, an embodiment of the present invention provides a method for reconstructing a transmission tower mechanical finite element model, including the following steps:
step 1, obtaining a three-dimensional solid model of the transmission tower through laser point cloud scanning, processing the three-dimensional solid model to obtain point cloud coordinate data of each angle steel in the model, identifying the total number of the angle steels in the solid model of the transmission tower and key point coordinate information of two ends of each angle steel, and calculating the cross sectional area of the end part of each angle steel;
in this embodiment, a three-dimensional solid model of a transmission tower is obtained through laser point cloud scanning, the solid model of the transmission tower subjected to the laser point cloud scanning is imported into finite element software as shown in fig. 2, then, a command stream code is used for circularly recognizing that the total number of angle steels in the solid model of the transmission tower is 74, key point coordinate information and angle steel sectional area of each angle steel endpoint, and data information of each angle steel is gradually recognized and processed and stored in an array.
Then exporting the arrays of the coordinate information and the cross-sectional area of the key points of the stored angle steel to a specified data file, outputting the total number '74' of the angle steel as data of a first row to be used as subsequent cyclic variables, and starting from a second row, the entity number n of each angle steel, the number '14' of the key points of the cross sections at two ends and the cross-sectional area S n The third row is the x, y and z coordinates of the key point at the end of each angle steel, the first 7 rows are the coordinates of the key point at one end of the angle steel, and the last 7 rows are the coordinates of the key point at the other end of the angle steel.
Step 2, identifying and screening out the key point coordinates of the outermost corner nodes in the end sections of the angle steel according to the key point coordinate data of the cross sections of the two ends of each angle steel, and taking the key point coordinates as the node coordinates of the angle steel unit;
the step can be processed in MATALB software, key point coordinate data of two ends of each angle steel are read in a circulating mode, the distance between two points is calculated, whether the unit vector orthogonal inner product of two lines connected with a certain key point of each end of the angle steel is close to 0 is judged in a circulating mode, then whether the two edge lengths are close to equal is judged, so that the outermost side and the innermost side of each end of the angle steel are screened out, and finally the two edge lengths connected with the angle nodes and the maximum coordinate are selected as the outermost side angle node coordinate, namely the No. 15 node coordinate in one angle steel shown in the figure 3.
Step 3, carrying out adjacent node fusion processing on node coordinates of all the angle steel units, circularly calculating the distance between every two angle steel units, and fusing the node coordinates of the angle steel units which meet the requirement of being close to each other within a certain range into a common node coordinate;
in this embodiment, when the adjacent node fusion processing is performed, the end point coordinate with the large cross-sectional area of the angle steel is selected to be unchanged, and the end point coordinate with the small cross-sectional area of the adjacent angle steel is modified, so that the coordinates of the end points are finally fused into the end point coordinate of the angle steel with the largest cross-sectional area. Specifically, coordinate data and sectional area information of two ends of a first angle steel are taken first, then coordinate data and sectional area information of two ends of a subsequent angle steel are taken in a circulating mode, the distance between every two end points of the two angle steels is calculated sequentially, as shown in fig. 4, end point coordinates 1 and 2 of a first steel structure are taken, end point coordinates 3 and 4 (1 (2) (3) (4) in fig. 4) of a second steel structure are taken, the distances between the end points 1, 3,1, 4,2, 3,2 and 4 are calculated respectively, if one of the conditions meets the requirement that the distance is very short, the two end point coordinates are fused into a common end point coordinate, the sectional areas of the two angle steels need to be compared and judged during node fusion, the end points with large sectional areas are unchanged, and the end points of the other angle steel with small sectional areas are modified to enable the coordinate information to be consistent. And then, comparing and fusing the coordinates of the third steel structure and the first steel structure, and repeating the steps, wherein the final intersection points close to the coordinate positions are modified into the end point coordinates of the angle steel with the largest sectional area.
Step 4, identifying the condition that the auxiliary material and the main material are crossed in the transmission tower solid model, adding node coordinates on the main material unit at the crossed point, segmenting the main material unit, modifying the node coordinates of the auxiliary material angle steel near the crossed point into the added node coordinates on the main material, and then performing adjacent node fusion processing again;
in this example, the cross-sectional area S of the first angle steel is taken 1 And the cross sections are sequentially circulated with the cross section S of the rear angle steel i Comparing every two of them, when the absolute value of S is satisfied 1 -S i |≥0.1cm 2 If the distance between one end point and the main material is within 0.01m, the situation that the auxiliary material is crossed with the main material is judged, then the distance between the end points of the angle steels is calculated, and the judgment condition is set as that the end point of the auxiliary material and the two end points of the main material are crossedThe distance between the main material and the auxiliary material is more than 0.1m, so that the condition of crossing the main material and the auxiliary material without newly adding a split point and the condition of crossing the main material and the auxiliary material only at the end points can be eliminated;
then, the coordinates of the newly added nodes on the main material at the intersection points are calculated, and the coordinates of the end points of the auxiliary materials which are intersected with the main material are set as P (x) 0 ,y 0 ,z 0 ) The main line segment is AB, and the coordinate of the A endpoint is (x) 1 ,y 1 ,z 1 ) And the coordinates of the B endpoint are (x) 2 ,y 2 ,z 2 ) Drawing a perpendicular line of the line segment AB from the point P, and calculating the parameter t in the parameter equation of the straight line AB, wherein the vertical line is Q (x, y, z) 1 Will t 1 Substituting into a linear equation to obtain the coordinate (X, Y, Z) of the foot drop point as the coordinate of the newly added point 1 Will t 1 The coordinates (X, Y, Z) of the foot drop point can be obtained as the coordinates of the newly added point by substituting the coordinates into a linear equation, wherein,
Figure BDA0003744673970000081
Figure BDA0003744673970000082
and finally, modifying the node coordinates of the angle steel of the auxiliary material near the intersection point into the calculated newly added node coordinates on the main material, and then performing adjacent node fusion processing again.
Step 5, identifying the crossing condition of the auxiliary materials and the auxiliary materials in the transmission tower solid model, calculating node coordinates at the crossing points as new node coordinates, and segmenting crossed auxiliary material units;
firstly, extracting the cross section area S of a first auxiliary material angle steel 1 And the first circulation is performed with the cross section area S of the rear angle steel i Performing pairwise comparison, when | S is satisfied 1 -S i |≤0.05cm 2 If the distance between any two end points is small, the distance between any two end points is calculatedJudging that the auxiliary material and the auxiliary material are crossed in an X shape according to the length of the angle steel, calculating the node coordinates of the X-shaped crossed point of the auxiliary material, and setting that a point P (X, Y, Z) exists on a first angle steel line segment AB, and the coordinate of an A endpoint is (X) 1 ,y 1 ,z 1 ) And the coordinates of the B endpoint are (x) 2 ,y 2 ,z 2 ) A point Q (U, V, W) exists on the second angle steel line segment CD, and the coordinate of the C endpoint is (x) 3 ,y 3 ,z 3 ) D end point coordinate is (x) 4 ,y 4 ,z 4 ) If the auxiliary material X is crossed, the shortest distance between two angle steel line segments is PQ, and a parameter s in a linear equation is calculated 2 And t 2 Partial derivative f(s) of the square of the PQ shortest distance of 2 ,t 2 ) And let the partial derivative f(s) 2 ,t 2 ) Finding s for 0 2 And t 2 A parameter t 2 Substituting into linear equation to calculate node coordinate K (X) of cross point 1 ,Y 1 ,Z 1 ):
Figure BDA0003744673970000091
The distance between two points of PQ is:
Figure BDA0003744673970000092
partial derivative f(s) of PQ shortest distance squared 2 ,t 2 ) Comprises the following steps:
Figure BDA0003744673970000093
node coordinate K (X) of the intersection 1 ,Y 1 ,Z 1 ):
Figure BDA0003744673970000094
And then modifying the end point coordinates of the auxiliary material angle steel crossed at the moment, and dividing the angle steel into four parts after adding the end point coordinates of the crossed point, so that two entities need to be added to the original matrix, and the entity numbers are sequentially increased to 2 in the cycle.
And 6, extracting the coordinates of the nodes at the two ends of all the angle steel units in the transmission tower and the data of the cross sectional area of the end part of the angle steel obtained in the steps 1-5, and reconstructing a mechanical finite element model of the whole transmission tower according to the extracted coordinates of the nodes at the two ends and the data of the cross sectional area of the end part of the angle steel.
And finally, re-identifying and extracting the coordinates and the cross-sectional area data of the angle iron end points in the transmission tower model in the example of completing the process of the steps, automatically writing a parameterized command stream file required for modeling by utilizing MATLAB codes, and quickly completing modeling reconstruction of the transmission tower mechanics finite element model shown in figure 5 in finite element software by loading the generated parameterized command stream file, wherein the reconstructed transmission tower mechanics finite element model can be directly used for mechanics simulation calculation and analysis, so that the time of manual modeling in the traditional method is greatly reduced, and the modeling efficiency is improved.
While the invention has been described with reference to a preferred embodiment, it will be understood by those skilled in the art that various changes in form and detail may be made therein without departing from the spirit and scope of the invention.

Claims (8)

1. A method for reconstructing a power transmission tower mechanics finite element model is characterized by comprising the following steps:
step 1, obtaining a three-dimensional solid model of a transmission tower through laser point cloud scanning, processing the three-dimensional solid model to obtain point cloud coordinate data of each angle steel in the model, identifying the total number of the angle steels and key point coordinate information of two ends of each angle steel in the solid model of the transmission tower, and calculating the cross sectional area of the end part of each angle steel;
step 2, identifying and screening out the key point coordinates of the outermost corner nodes in the end sections of the angle steel according to the key point coordinates of the cross sections of the two ends of each angle steel, and taking the key point coordinates as the node coordinates of the angle steel unit;
step 3, carrying out adjacent node fusion processing on the node coordinates of all the angle steel units;
step 4, identifying the crossing condition of the auxiliary material and the main material in the transmission tower solid model, adding node coordinates on the main material unit at the crossing point, segmenting the main material unit, modifying the node coordinates of the auxiliary material angle steel near the crossing point into the added node coordinates on the main material, and then performing adjacent node fusion processing again;
step 5, identifying the condition that the auxiliary material and the auxiliary material in the transmission tower solid model are crossed in an X mode, calculating node coordinates at the crossed points as newly-added node coordinates, and meanwhile segmenting crossed auxiliary material units;
and 6, extracting the node coordinates at the two ends of all the angle steel units in the transmission tower and the cross sectional area data of the end parts of the angle steel obtained in the steps 1-5, and reconstructing a mechanical finite element model of the whole transmission tower according to the extracted node coordinates and the cross sectional area data of the end parts of the angle steel.
2. The method for reconstructing a finite element model of transmission tower mechanics according to claim 1, wherein in step 2, the specific method for screening out the coordinates of the key points of the outermost corner nodes in the section of the end of the angle steel comprises the following steps:
the method comprises the steps of firstly calculating the distance between any two key points according to key point coordinate data of the end part of the angle steel, then calculating whether the unit vector inner product of two lines connected with a certain key point is close to 0 to screen two edges which are perpendicular to each other, then comparing whether the lengths of the two edges are close to screen the angle node coordinates of the outermost side and the innermost side of the end part of the angle steel, and finally selecting the length of the two edges connected with the node and the largest node coordinate as the angle node coordinate of the outermost side.
3. The method for reconstructing a finite element model of a power transmission tower mechanism according to claim 1, wherein in the step 3, when the adjacent node fusion processing is performed, the node coordinate with the larger cross-sectional area of the angle steel is selected to be unchanged, and the node coordinate with the smaller cross-sectional area of the adjacent angle steel is modified, so that the coordinates of the node coordinate are finally fused into the node coordinate of the angle steel with the largest cross-sectional area.
4. The method for reconstructing a finite element model of a transmission tower mechanics according to claim 1, wherein step 4 specifically comprises:
step 4.1, extracting the cross-sectional area S of the first angle steel 1 And the cross section area S of the angle steel is sequentially matched with that of the rear angle steel i Comparing two by two, when satisfying | S 1 -S i |≥ΔS 1 If so, one of the two angle steels is a main material and the other is an auxiliary material, wherein delta S is 1 Setting according to actual needs, determining an angle steel unit with a large cross-sectional area as a main material, calculating the distance between two end points of an auxiliary material and a main material line segment, and if the distance between one end point and the main material is delta L 1 If the content is within the range, the auxiliary material and the main material are crossed, wherein the content is Delta L 1 Setting according to actual needs;
4.2, on the premise that the auxiliary material and the main material are crossed, calculating the distance between the end points of the angle steel, and setting a judgment condition that the distance between the end point of the auxiliary material and the two end points of the main material is larger than a certain value so as to eliminate the condition that the main material and the auxiliary material are only crossed at the end points;
and 4.3, finally calculating the coordinates of the newly added nodes on the main material at the intersection points, and segmenting the main material according to the newly added nodes.
5. The method for reconstructing a finite element model of transmission tower mechanics according to claim 1, wherein the method for calculating the coordinates of the newly added nodes on the main material at the intersection in step 4.3 comprises:
let the coordinate of the end point of the auxiliary material crossed with the main material be P (x) 0 ,y 0 ,z 0 ) The main line segment is AB, and the coordinate of the A endpoint is (x) 1 ,y 1 ,z 1 ) And the coordinates of the B endpoint are (x) 2 ,y 2 ,z 2 ) Drawing a perpendicular line of the line segment AB from the point P, and calculating a parameter t in a parameter equation of the straight line AB, wherein the perpendicular is Q (x, y, z) 1 Let t be 1 Substituting into linear equation to obtain foot drop point coordinates (X, Y, Z) as new point coordinates, wherein,
Figure FDA0003744673960000021
Figure FDA0003744673960000022
6. The method for reconstructing a finite element model of a transmission tower mechanics according to claim 1, wherein the step 5 specifically comprises:
step 5.1, extracting the cross-sectional area S of the first auxiliary material angle steel 1 And the cross sections are sequentially circulated with the cross section S of the auxiliary angle steel at the back i Comparing every two of them, when the absolute value of S is satisfied 1 -S i |≤ΔS 2 Then, the cross sectional area of the two angle steels is close to each other, delta S 2 Setting according to actual needs, calculating the distance between every two end points of the angle steel, and if the distance between any two end points is smaller than the length of the angle steel, judging that the auxiliary material and the auxiliary material are crossed;
step 5.2, calculating node coordinates of the auxiliary material intersection points;
and 5.3, modifying the endpoint coordinates of the auxiliary material angle steel crossed at the moment, and segmenting the auxiliary material angle steel from the crossed nodes.
7. The method for reconstructing the finite element model of the power transmission tower mechanics according to claim 1, wherein the method for calculating the node coordinates of the auxiliary material intersection point comprises the following steps:
let a point P (X, Y, Z) exist on the first angle steel line segment AB, and the coordinate of the A endpoint is (X) 1 ,y 1 ,z 1 ) And the coordinates of the B endpoint are (x) 2 ,y 2 ,z 2 ) A point Q (U, V, W) exists on the second angle steel line segment CD, and the coordinate of the C endpoint is (x) 3 ,y 3 ,z 3 ) D end point coordinate is (x) 4 ,y 4 ,z 4 ) The shortest distance exists between two angle steel line sections under the condition of auxiliary material intersectionCalculating the parameter s in the equation of a straight line when the distance is PQ 2 And t 2 Partial derivative f(s) of the square of the PQ shortest distance of 2 ,t 2 ) And let the partial derivative f(s) 2 ,t 2 ) Finding s for 0 2 And t 2 A parameter t 2 Substituting into linear equation to calculate node coordinate K (X) of cross point 1 ,Y 1 ,Z 1 ):
Figure FDA0003744673960000031
The distance between two points PQ is:
Figure FDA0003744673960000032
partial derivative f(s) of PQ shortest distance squared 2 ,t 2 ) Comprises the following steps:
Figure FDA0003744673960000033
node coordinate K (X) of the intersection 1 ,Y 1 ,Z 1 ):
Figure FDA0003744673960000041
8. The method for reconstructing a finite element model of a transmission tower mechanics according to claim 1, wherein step 6 specifically comprises:
extracting the end point coordinates and the cross-sectional area data of the angle steel in the power transmission tower model obtained in the step 1-5, writing a parametric modeling command stream file required for establishing a mechanical finite element model by using MATLAB, and rapidly completing modeling reconstruction of the power transmission tower mechanical finite element model in finite element software by loading the generated parametric modeling command stream file.
CN202210827822.1A 2022-07-13 2022-07-13 Reconstruction method of transmission tower mechanics finite element model Pending CN115310172A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210827822.1A CN115310172A (en) 2022-07-13 2022-07-13 Reconstruction method of transmission tower mechanics finite element model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210827822.1A CN115310172A (en) 2022-07-13 2022-07-13 Reconstruction method of transmission tower mechanics finite element model

Publications (1)

Publication Number Publication Date
CN115310172A true CN115310172A (en) 2022-11-08

Family

ID=83857405

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210827822.1A Pending CN115310172A (en) 2022-07-13 2022-07-13 Reconstruction method of transmission tower mechanics finite element model

Country Status (1)

Country Link
CN (1) CN115310172A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117893711A (en) * 2024-03-14 2024-04-16 大连理工大学 Reconstruction method of entity model of special-shaped curved surface reinforced structure

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117893711A (en) * 2024-03-14 2024-04-16 大连理工大学 Reconstruction method of entity model of special-shaped curved surface reinforced structure

Similar Documents

Publication Publication Date Title
CN112257597B (en) Semantic segmentation method for point cloud data
CN105261068A (en) Micro-CT technology-based reservoir core three-dimensional entity model reconstruction method
CN113379699A (en) Transmission line insulator defect detection method based on deep learning
Qin et al. Automated reconstruction of parametric bim for bridge based on terrestrial laser scanning data
CN113628291B (en) Multi-shape target grid data vectorization method based on boundary extraction and combination
CN106600700B (en) Three-dimensional model data processing system
CN115310172A (en) Reconstruction method of transmission tower mechanics finite element model
CN114708414B (en) Engineering three-dimensional terrain data cutting fusion method and system
CN111652171B (en) Construction method of facial expression recognition model based on double branch network
CN104484527B (en) Uniform load automatic dynamic amending method in a kind of discrete topology process of topology optimization
CN104462414A (en) Topological structure based flow chart similarity method
CN107679332A (en) Electric power pylon BIM model building methods based on Tekla platforms
CN116049941B (en) Method for extracting and analyzing multidimensional state of assembled ring truss structural member before assembly
CN106227919B (en) Manifold learning-based dynamic simulation visualization method for power system
Choi et al. Tetrahedral mesh generation based on advancing front technique and optimization scheme
CN114612751B (en) Whole machine point cloud data down-sampling method based on semantic learning
CN116305418A (en) Method for converting grid structure BIM model and point cloud model into finite element model
CN116664758A (en) Tetrahedral grid self-adaptive hierarchical refinement method based on fault plane
Li et al. Reasoning mechanism: An effective data reduction algorithm for on-line point cloud selective sampling of sculptured surfaces
EP3608841A1 (en) Method for comparing fittingness of two-dimensional drawing in engineering
CN1605961A (en) Direct numerical control processing method based on measuring data points
CN114547761A (en) High-precision model-based prefabricated part steel bar optimization method and system
Sarfraz et al. Visualization of data using genetic algorithm
CN109117496B (en) Three-dimensional simulation design method and system for temporary construction arrangement of transformer substation engineering
CN111444614A (en) Flow field reconstruction method based on graph convolution

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination