Method, system and computer software product for rapid partitioning of tire complex pattern grid
Technical Field
The invention relates to the technical field of simulation design of tires, in particular to a method, a system and a computer software product for rapidly dividing complex pattern grids of tires.
Background
With the increasing demands on tire performance, particularly in terms of running safety, durability, and comfort, finite Element Analysis (FEA) has become a key tool for evaluating tire designs. The design of the complex pattern of the tire has important influence on the performances such as grip, wear resistance, noise control and the like, so that the characteristics of the complex pattern of the tire need to be fully considered when various performances of the tire are accurately analyzed. However, conventional meshing methods face challenges of modeling complexity, time consumption, and reliance on third party business software when dealing with tire complex patterns.
In the prior art, a three-dimensional model of a tire tread is generally generated by three-dimensional geometric modeling software, and then is subjected to meshing by means of third-party meshing software (such as HYPERMESH and the like). For example, the Chinese patent application (publication No. CN116796604A, publication No. 2023-09-22) filed by the applicant can be used for ensuring that the coordinates of nodes on the adjacent surfaces behind the tread array with certain thickness are completely consistent after grid division by taking a tread with certain thickness (the thickness can be defined independently according to the requirement) from the intercept tread pattern in the radial direction. Not only does this approach require significant computational resources and specialized software support, it can take 8 to 20 hours in practical applications, resulting in significant increases in cost and time overhead. Therefore, this partitioning approach, which relies on three-dimensional modeling, presents an inefficient bottleneck in industrial applications, which makes it difficult to meet the requirements of rapid iterative design and performance evaluation.
To solve these problems, current technical research is mainly focused on the aspects of improving the grid division efficiency through an optimization algorithm and developing alternative simplified methods to reduce the dependence on three-dimensional modeling. However, most of these methods cannot achieve both high precision and low cost, and it is still difficult to meet the application requirements of the tire pattern in high-precision simulation analysis. Therefore, there is an urgent need for a tire complex pattern grid division method that does not rely on three-dimensional modeling, is easy and efficient to operate, to reduce the computational cost, improve the division efficiency, and accelerate the design feedback.
Disclosure of Invention
In order to solve the technical problems, the invention provides a method for rapidly dividing a tire complex pattern grid based on plane unfolding projection, according to the method, the two-dimensional plane graph of the tire pattern is directly subjected to grid division, and then plane grid nodes and units are projected to the three-dimensional tread curve of the tire, so that a three-dimensional grid structure is constructed. Compared with the traditional method, the method does not depend on three-dimensional geometric modeling software, the whole dividing process is completed in about 1 hour, the grid dividing efficiency is greatly improved, and the cost of design and analysis is reduced. In addition, the plane grid division and projection scheme adopted by the method can accurately capture the geometric characteristics of the tire pattern, so that the geometric characteristics can be kept high in the finite element simulation.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
a method for rapid meshing of a complex tread grid of a tire, said method comprising the steps of:
1) Grid division is carried out on the tire intercept or the plane graph of all patterns, the patterns are divided into quadrilateral or triangle units, and a plane coordinate system is constructed by an x axis and a y axis;
2) Drawing a contour curve, generating a tread contour curve, a groove bottom curve and a pattern bottom curve according to design parameters of the tread contour of the tire, and adding smooth transition curves between different curves;
3) Generating node and unit information of a pattern plan grid, defining node coordinates as (xi, yi), wherein a unit comprises a triangle or quadrilateral grid;
4) Expanding nodes and units of the pattern plane grid to a tread curve, and defining node projection coordinates on the tread curve;
5) Projecting the nodes into the patterns, and generating internal projection points according to the intersection points of the normal lines and the profile curves;
6) And forming a three-dimensional unit structure and generating a final three-dimensional grid division.
Preferably, the profile curve drawing in the step 2) comprises drawing smooth curves according to different pattern depths to form tread pattern bottom profile lines with gradual change depths so as to increase grid density and optimize simulation precision, and sequentially naming each curve as C i and i= 0~n according to the height removing positions, wherein n is the number of the curves minus 1, C 0 is a crown curve, C n is a pattern bottom curve, the curves are required to be positioned according to actual positions, and a tire center point is taken as an origin.
Preferably, the node and cell information generation in step 3) includes defining node numbers and cell numbers to ensure that the nodes of the triangle and quadrilateral mesh have unique identifications and that node locations are trackable.
Preferably, the coordinate system of the node in step 4) uses the center of the tire as the origin, the radial direction as the z-coordinate, the wheel axis direction as the y-coordinate, and the direction perpendicular to the two directions as the x-axis. The tire circumferential angle occupied by the pattern is θ (as shown in fig. 3), the x coordinate of the node T c0i corresponding to the node N i on the tread curve is from the left end of the tread curve, the x coordinate and the y coordinate of the point corresponding to the point with the curve length equal to x i are x ri,yri, and the x, y and z coordinate values of the projection point T c0i are respectively:
The constituent node numbers of the projected unit M i are the same as E i.
Preferably, in step 5), a node T c0i on the tread curve is projected into the tread, where i is the node number. For the nodes inside the pattern block, calculating the normal line of the pattern contour corresponding to the position of the node, marking as n i, calculating the intersection point of the normal line n i and other contours C v, marking as T cvi, wherein v=1-n is the internal curve number, i is the node number, and if the normal line and the contour C v have no intersection point, the intersection point coordinates are the same as C v-1.
Preferably, the initial projection curve C b of the cell M i is specified in step 6), where b= 0~n, which means that a three-dimensional cell is formed from below the curve C b, and the node corresponding to the three-dimensional cell S vi (where v represents the layer number, corresponds to the initial contour line, and i represents the cell number) is (T cvm,Tcvk,Tcvp,Tc(v+1)m,Tc(v+1)k,Tc(v+1)p) or (T cvm,Tcvk,Tcvp,Tcvd,Tc(v+1)m,Tc(v+1)k,Tc(v+1)p,Tc(v+1)d), where v= 0~n-1. Thereby forming a three-dimensional unit.
Further, the invention also provides a system for rapidly dividing the complex tread grid of the tire, which realizes the method and comprises the following modules:
the data input module is used for receiving the tire pattern plan and the profile parameters thereof;
the processing module is used for executing grid division and three-dimensional node generation of the tire pattern plan, and comprises contour curve drawing, node and unit information generation, node unfolding projection and internal projection generation;
the storage module is used for storing the generated node and unit data, including node plane coordinates and projection coordinates;
and the output module is used for outputting the three-dimensional grid data and providing a structured file for finite element simulation analysis.
Further, the invention also provides a computer device comprising a memory, a processor and a computer program stored on the memory, the processor executing the computer program to implement the method.
Further, the present invention also provides a computer-readable storage medium having stored thereon a computer program or instructions which, when executed by a processor, implement the method.
Further, the invention also provides a computer program product comprising a computer program or instructions which, when executed by a processor, implements the method.
By adopting the technical scheme, the two-dimensional plane graph of the tire pattern is directly subjected to grid division, and then plane grid nodes and units are projected to the three-dimensional tread curve of the tire, so that the three-dimensional grid structure is constructed. Compared with the traditional method, the method does not depend on three-dimensional geometric modeling software, the whole dividing process is completed in about 1 hour, the grid dividing efficiency is greatly improved, and the cost of design and analysis is reduced. In addition, the plane grid division and projection scheme adopted by the method can accurately capture the geometric characteristics of the tire pattern, so that the geometric characteristics can be kept high in the finite element simulation. The following beneficial effects are embodied:
1) The invention directly carries out grid division on the two-dimensional plane graph of the complex pattern, and then projects the plane grid to the three-dimensional tread curve, thereby omitting the step of modeling the three-dimensional geometry of the pattern in the traditional method. The method obviously reduces the time of grid division, the whole process can be completed within 1 hour, the time is shortened by about 4 times compared with the traditional method, and the working efficiency is greatly improved.
2) The invention does not need to rely on third-party three-dimensional modeling software, realizes grid division by own method and system, greatly reduces the dependence and use cost on expensive software, and reduces the cost of enterprises in terms of computing resources and business software authorization.
3) The method is suitable for complex pattern structures, based on plane graph grid division and unfolding projection, can accurately capture and reproduce complex tire pattern structures, and supports layered projection processing of patterns with different depths through multi-level contour curve design, so that a high-precision three-dimensional grid structure conforming to the geometric characteristics of actual tire patterns is formed.
4) Optimizing simulation precision, namely, increasing node density control between a tread curve and a groove bottom curve, so that grid division has enough resolution on different depths, geometric details of tire patterns can be better represented, the precision of simulation results is improved, and the method is suitable for high-precision finite element analysis of tire performance.
5) The method is simple and convenient to operate and easy to integrate, and the method provides a simplified dividing flow, does not need complex three-dimensional modeling and adjustment, and is simpler and more convenient to operate. In addition, through systematic module design, the method is easy to integrate into the existing finite element analysis workflow, and is convenient for engineering personnel to carry out rapid pattern grid division and performance evaluation.
6) The method and the system can realize automatic processing, and the whole grid dividing process can be automatically completed through standardized input, processing and output module design, so that manual intervention is reduced, the production efficiency and consistency are improved, and the method and the system are suitable for industrial scenes needing rapid iterative design and performance evaluation.
In summary, the invention provides a high-efficiency and low-cost method and a system for dividing complex pattern grids of tires, which can remarkably improve the efficiency and the precision of pattern grid division and are suitable for being applied to rapid iteration of tire design and performance analysis.
Drawings
FIG. 1 is a schematic diagram of a pattern intercept meshing and coordinate system;
FIG. 2 is a schematic diagram of a pattern profile;
FIG. 3 is a schematic diagram of a three-dimensional pattern coordinate system;
FIG. 4 is an embodiment of a grid division;
FIG. 5 is an embodiment pattern profile;
FIG. 6 is a node projection view;
fig. 7 is a three-dimensional patterned grid pattern generated.
Detailed Description
The technical solutions in the embodiments are clearly and completely described below in connection with the embodiments of the present invention, and it is obvious that the described embodiments are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The invention provides a method for rapidly dividing a complex pattern grid of a tire based on plane unfolding projection, which directly divides the grid by a two-dimensional plane graph of the tire pattern, and then the plane grid nodes and the units are projected to a three-dimensional tread curve of the tire, so that a three-dimensional grid structure is constructed. The method specifically comprises the following steps of:
The first step is to divide the plane view of the tire intercept or all patterns into quadrangles or triangles, wherein the width direction of the patterns is x-axis, the length direction of the patterns is y-axis, the leftmost end of the patterns is x-axis origin, the right direction is positive direction, the bottommost position of the patterns is y-axis origin, and the upward direction is positive direction, as shown in figure 1.
And step two, drawing a contour curve. And drawing a tread profile curve according to the designed tread profile of the tire. Since tire tread patterns have patterns of different depths, it is necessary to draw smooth curves through the bottoms of the grooves according to the different pattern depths. And drawing a smooth pattern bottom curve according to the position of the deepest pattern groove bottom. To increase the mesh density, a smooth curve may be arbitrarily added between the tread curve, the groove bottom curve, and the pattern bottom curve. Each curve is named as C i and i= 0~n in sequence according to the height removing position, wherein n is the number of curves minus 1, C 0 is a crown curve, C n is a pattern bottom curve, the curves need to be positioned according to the actual position, and the center point of the tire is taken as the origin, as shown in fig. 2.
And thirdly, generating node and unit information of the pattern plan grid. The plane coordinate of the node N i is (x i,yi), and the unit E j is formed by the node (N m,Nk,Np) or (N m,Nk,Np,Nd) and corresponds to the triangular grid and the quadrilateral grid respectively, wherein m, k, p and d are node numbers.
And fourthly, expanding the nodes and the units of the pattern plane grid to the tread curve. The coordinate system of the node takes the center of the tire as an origin, the radial direction as a z coordinate, the wheel axis direction as a y coordinate and the direction perpendicular to the two directions as an x axis. The tire circumferential angle occupied by the pattern is θ (as shown in fig. 3), the x coordinate of the node T c0i corresponding to the node N i on the tread curve is from the left end of the tread curve, the x coordinate and the y coordinate of the point corresponding to the point with the curve length equal to x i are x ri,yri, and the x, y and z coordinate values of the projection point T c0i are respectively:
The constituent node numbers of the projected unit M i are the same as E i.
And fifthly, projecting a node T c0i on the tread curve into the pattern, wherein i is the node number. For the nodes inside the pattern block, calculating the normal line of the pattern contour corresponding to the position of the node, marking as n i, calculating the intersection point of the normal line n i and other contours C v, marking as T cvi, wherein v=1-n is the internal curve number, i is the node number, and if the normal line and the contour C v have no intersection point, the intersection point coordinates are the same as C v-1.
And sixth, forming a three-dimensional unit. The initial projection curve C b of cell M i is specified, where b= 0~n, indicating that a three-dimensional cell is formed from below the C b curve. The node corresponding to the three-dimensional cell S vi (where v represents a layer number, corresponds to a starting contour, and i represents a cell number) is (T cvm,Tcvk,Tcvp,Tc(v+1)m,Tc(v+1)k,Tc(v+1)p) or (T cvm,Tcvk,Tcvp,Tcvd,Tc(v+1)m,Tc(v+1)k,Tc(v+1)p,Tc(v+1)d), where v= 0~n-1. Thereby forming a three-dimensional unit.
The three-dimensional grid division of the tire complex pattern is completed through the steps, the node set T cvi and the unit set S i are obtained, a three-dimensional model is not required to be established in the whole calculation process of generating the grid, and the grid division time is greatly shortened only by about 1 hour in the whole process.
The following is a specific example of a common tire pattern, which is an intercept:
the first step is to divide the plane view of the tire intercept or all patterns into quadrangles or triangles, wherein the width direction of the patterns is x-axis, the length direction of the patterns is y-axis, the leftmost end of the patterns is x-axis origin, the right direction is positive direction, the bottommost position of the patterns is y-axis origin, and the upward direction is positive direction, as shown in figure 4.
And step two, drawing a contour curve. And drawing a tread profile curve according to the designed tread profile of the tire. Since tire tread patterns have patterns of different depths, it is necessary to draw smooth curves through the bottoms of the grooves according to the different pattern depths. And drawing a smooth pattern bottom curve according to the position of the deepest pattern groove bottom. To increase the mesh density, a smooth curve may be arbitrarily added between the tread curve, the groove bottom curve, and the pattern bottom curve. The 5 curves are sequentially named as C i, i= 0~n according to the height removing position, wherein n=4, C 0 is a crown curve, C 4 is a pattern bottom curve, the curves are positioned according to the actual position, the position of the highest contour point is 400mm, and the center point of the tire is taken as the origin, as shown in fig. 5.
Third, generating nodes (such as table 1) and unit information (such as table 2) of the pattern plan grid.
Table 1 grid node coordinates of planar patterns
Node ID |
Coordinate X |
Coordinate Y |
1 |
10.457295 |
-84.078747 |
2 |
13.658201 |
-84.078747 |
3 |
13.89152 |
-87.500364 |
4 |
18.054475 |
-87.500364 |
5 |
17.682803 |
-84.078747 |
6 |
13.38896 |
-80.130324 |
7 |
17.253906 |
-80.130324 |
8 |
10.126329 |
-80.130324 |
9 |
9.571526 |
-75.76002 |
10 |
12.578976 |
-75.76002 |
11 |
16.482978 |
-75.76002 |
12 |
8.97809 |
-72.238201 |
Table 2 planar pattern element information
Element ID |
Node list |
1 |
1, 2, 3 |
2 |
2, 4, 3 |
3 |
5, 4, 2 |
4 |
6, 7, 5, 2 |
5 |
1, 8, 6, 2 |
6 |
8, 9, 10, 6 |
7 |
10, 11, 7, 6 |
8 |
9, 12, 13, 10 |
9 |
13, 14, 11, 10 |
10 |
12, 15, 16, 13 |
11 |
16, 17, 14, 13 |
12 |
15, 18, 19, 16 |
13 |
19, 20, 17, 16 |
Node N i plane coordinates are (x i,yi), unit E j is defined by node (N m,Nk,Np) or (N m,Nk,Np,Nd), corresponding to triangle mesh and quadrilateral mesh, respectively, where m, k, p and d are node numbers, e.g., node N 1 plane coordinates are (x 1=10.457295,y1 = -84.078747), node of unit E 1 is defined as (N 1,N2,N3), and node of unit E 4 is defined as (N 6,N5,N7,N2).
And fourthly, expanding the nodes and the units of the pattern plane grid to the tread curve. The coordinate system of the node takes the center of the tire as an origin, the radial direction as a z coordinate, the wheel axis direction as a y coordinate and the direction perpendicular to the two directions as an x axis. The pattern occupies a tire circumferential angle θ=10° (as shown in fig. 3). Taking the node N 15 as an example, the x coordinate of the node T c015 corresponding to the plane coordinate of the node N 15 being (x 15=30.5,y15 =0) on the tread curve is from the left end of the tread curve, the x coordinate and the y coordinate of the point corresponding to the point with the curve length equal to x 15 =30.5 are denoted as x r15=-70.6,yr15 = 398.5, and the x, y and z coordinate values of the projection point T c015 are respectively:
The projected element M i has the same constituent node number as E i, i.e., element M 1 has the same constituent node number (N TC01,NTC02,NTC03) as E 1 (N 1,N2,N3).
Taking the node T c015 as an example, projecting the node T c015 on the tread curve into the pattern, wherein i is the node number. For the nodes inside the pattern block, calculating the normal line of the pattern contour corresponding to the position of the node, namely n 15, calculating the intersection point of the normal line n 15 and other contours C v, namely T cvi, wherein v=1-n is the internal curve number, i is the node number, and obtaining the coordinates of T c115,Tc215,Tc315,Tc415 as shown in fig. 6.
And sixth, forming a three-dimensional unit. Taking M 1 as an example, the initial projection curve C 1 of the designated cell M 1 indicates that a three-dimensional cell is formed from below the curve C 1. The node corresponding to the three-dimensional unit S 11 is (T c11,Tc12,Tc13,Tc21,Tc22,Tc23), the node corresponding to the three-dimensional unit S 21 is (T c21,Tc22,Tc23,Tc31,Tc32,Tc33), and the node corresponding to the three-dimensional unit S 31 is (T c31,Tc32,Tc33,Tc41,Tc42,Tc43). Thereby forming a three-dimensional unit as shown in fig. 7.
The three-dimensional grid division of the tire complex pattern is completed through the steps, the node set T cvi and the unit set S vi are obtained, a three-dimensional model is not required to be established in the whole calculation process of generating the grid, the whole process only needs about 45 minutes, and the grid division time is greatly shortened.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art. The generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.