Method for quickly modifying tire complex pattern three-dimensional finite element model

Abstract

The invention relates to a method for rapidly modifying a tire complex pattern three-dimensional finite element model, and relates to the technical field of tire pattern finite element modelling, and the method comprises the following steps: taking a tire circumferential direction as an X axis, a tire width direction as a Y axis, and a tire radial direction as a Z axis, carrying out three-dimensional grid division on a tire pattern three-dimensional model to obtain an initial grid; the center O point of the initial grid on the XY plane is the origin of coordinates of the X axis and the Y axis; and performing angle change or mirroring on the initial grid, wherein the angle change is to establish mathematical mapping of a pattern angle and node displacement through definite integral or interpolation. The mirroring is a mirrored grid generated by mirroring the initial grid relative to the XZ surface or the YZ surface; normal information of the mirrored grid is updated, so that the determinant of the Jacobian matrix of the grid unit is positive; according to the method, the steps of adjusting the pattern three-dimensional model and reconstructing the finite element mesh are omitted, the pattern research and development period is greatly shortened, and the labor cost of pattern research and development is reduced.

Description

Technical Field

[0001] This invention relates to the field of tire tread finite element modeling technology, and in particular to a method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern. Background Technology

[0002] Tires are a crucial component of automobiles, and the arrangement of tire tread patterns affects a vehicle's noise characteristics and grip on wet, slippery surfaces. As tire performance requirements increase, tread pattern designs need frequent iterations to balance multiple performance objectives. Currently, tire simulation based on finite element analysis is the core method for verifying the performance of tread pattern designs.

[0003] Currently, the mainstream method for finite element modeling of tire tread patterns involves first constructing a 3D model of the tread pattern, including complex grooves and angular features, using 3D modeling software. This 3D model is then imported into finite element analysis software to generate a finite element mesh, which is then used for finite element simulation analysis. If the tread pattern needs adjustment based on the simulation results, the 3D model must be adjusted using the 3D modeling software. The adjusted 3D model is then imported back into the finite element analysis software to re-generate a finite element mesh, and the simulation analysis is performed again. Existing meshes cannot be reused. Each adjustment requires starting from scratch, resulting in repetitive work, significant time consumption, a long tread pattern development cycle, and high labor costs.

[0004] Therefore, it is necessary to develop a method that can directly modify the pattern angle on the finite element mesh, eliminating the need to adjust the 3D pattern model and reconstruct the finite element mesh, thus shortening the pattern development cycle and reducing labor costs. Summary of the Invention

[0005] To address the shortcomings of related technologies, this invention provides a method for rapidly modifying a three-dimensional finite element model of a tire's complex tread pattern. This method directly modifies the tread pattern angle on the finite element mesh, eliminating the need for adjusting the three-dimensional tread pattern model and reconstructing the finite element mesh. This solves the technical problems of long tread pattern development cycles and high labor costs.

[0006] This invention provides a method for rapidly modifying a three-dimensional finite element model of a tire with complex tread patterns, comprising: S1. Using the tire circumference as the X-axis, the tire width direction as the Y-axis, and the tire radial direction as the Z-axis, perform three-dimensional meshing on the tire tread model to obtain the initial mesh. The center point O of the initial mesh on the XY plane is the origin of the X-axis and Y-axis coordinates. S2. Shift or mirror the initial mesh to generate a modified mesh: Variable angle: The initial grid has multiple sampling points along the Y-axis, with two sampling points located at the two ends of the initial grid and the other sampling points located at the angle transitions of the initial grid pattern; a variable angle value τ is assigned to each sampling point. Calculate the X-axis change for each sampling bit. x, based on two adjacent sample bits x is calculated using a definite integral algorithm or an interpolation algorithm to determine the change in the X-axis of each grid cell between these two sampling points. The formula for calculating x is:

[0007]

[0008] Where y is the y-axis coordinate. x is the change along the X-axis. The curvature of the XY plane; Based on the calculated changes in the X-axis of the nodes, adjust the X-axis coordinates of each node in the grid; Mirroring: Mirrors the initial mesh about the XZ or YZ plane to generate a mirrored mesh; updates the normal information of the mirrored mesh to make the determinant of the Jacobian matrix of the mesh elements positive; S3, Output the modified mesh.

[0009] By adopting the above technical solution, the present invention can realize the local or overall angle adjustment of the finite element model of the tire pattern, and can also realize the overall mirroring of the finite element model of the tire pattern. It skips the steps of adjusting the three-dimensional model of the pattern and reconstructing the finite element mesh, dynamically links the adjustment of design parameters with mesh modification, directly modifies the existing mesh, greatly shortens the pattern development cycle, and reduces the labor cost of pattern development.

[0010] In some embodiments, a pattern unit is defined between two adjacent sampling positions, and multiple adjacent pattern units form a continuous pattern. Before calculating the change in the X-axis of the node in S2, all sampling positions in all pattern units involved in the continuous pattern are sequentially numbered along the Y-axis. Each sampling position... x superimposed on the previous sample bit x obtain Based on two adjacent sampling bits The change in the X-axis of each grid cell between these two sampling points is calculated using a definite integral algorithm or an interpolation algorithm.

[0011] By adopting the above technical solution, the problem of misalignment between adjacent pattern units in a continuous pattern can be avoided after the angle of some pattern units in the continuous pattern is adjusted by changing the angle. This avoids the impact of local angle adjustment on the overall continuity of the pattern.

[0012] In some embodiments, the entire initial mesh is viewed as a single continuous pattern, or as two continuous patterns symmetrical about the XZ plane.

[0013] By adopting the above technical solution, for all types of tire tread patterns, the entire initial grid can be regarded as a continuous pattern to achieve unidirectional integration along the positive or negative Y-axis, thereby improving the overall pattern continuity. For tire tread patterns symmetrical about the XZ plane, in addition to unidirectional integration, the entire initial grid can also be regarded as two continuous patterns symmetrical about the XZ plane. The two continuous patterns are integrated along the Y-axis towards point O, or the two continuous patterns are integrated along the Y-axis away from point O, to achieve bidirectional integration.

[0014] In some embodiments, after obtaining the change in the X-axis of a node in S2, the change in the X-axis of the node at Y=0 is used as... Subtract the X-axis coordinate of each node This yields the corrected X-axis coordinates of the nodes.

[0015] By adopting the above technical solution, the change in the X-axis at Y=0 of the corrected grid is 0, so that the X-axis grid coordinates at Y=0 can be used as the coordinate system correction reference for the modified grid.

[0016] In some embodiments, the definite integral algorithm type in S2 is any one of Simpson integral, trapezoidal rule integral, Gaussian integral, and adaptive step size control algorithm; the interpolation algorithm type is any one of piecewise quadratic interpolation and polynomial interpolation.

[0017] In some embodiments, in S2, mirroring about the XZ plane changes the node coordinates of all mesh cells from (x, y, z) to (x, -y, z); mirroring about the YZ plane changes the node coordinates of all mesh cells from (x, y, z) to (-x, y, z).

[0018] In some embodiments, the normal information in S2 is updated using finite element software through any of the following methods: flipping the normal direction of the mirrored element; displaying and storing the normal direction of the elements in the initial mesh and updating the normal direction of the mirrored element accordingly; swapping the output order of the nodes on the top and bottom surfaces of the mirrored element; and adjusting the node number corresponding to the node coordinates.

[0019] By adopting the above technical solution, the cell normal error caused by mirroring operation was corrected, ensuring the normal consistency of the cells.

[0020] In some embodiments, multiple sampling points are set at the angle transition of a single pattern.

[0021] By adopting the above technical solution, the curvature continuity of the overall mesh after local angle change operation is improved.

[0022] In some embodiments, the Y-axis interval in S1 is applicable to tire tread patterns of all width specifications, the Y-axis interval occupied by the initial grid is user data, and the X-axis offset and τ outside the user data are both zero.

[0023] By adopting the above technical solution, the boundary of the finite element mesh is prevented from abruptly changing or breaking during local angle changes, thus improving the stability and fault tolerance of the angle change operation.

[0024] In some embodiments, the initial mesh is divided using three-dimensional six-node elements or three-dimensional eight-node elements in S1.

[0025] Compared with the prior art, the beneficial effects of the present invention are as follows: 1. This invention can achieve local or overall angle adjustment of the finite element model of tire tread pattern, and can also achieve overall mirroring of the finite element model of tire tread pattern. It skips the steps of adjusting the 3D model of the tread pattern and reconstructing the finite element mesh, dynamically links the adjustment of design parameters with mesh modification, and directly modifies the existing mesh, which greatly shortens the tread pattern development cycle and reduces the labor cost of tread pattern development.

[0026] 2. This invention preserves the seed point distribution of the mirrored mesh and adjusts the node order while keeping the node coordinates unchanged, thus ensuring the normal consistency of the elements and avoiding mesh distortion caused by mirroring.

[0027] 3. The modified mesh based on this invention has a Jacobian value greater than 0.9, and the modified mesh has higher quality. Attached Figure Description

[0028] The accompanying drawings, which are included to provide a further understanding of the invention and form part of this application, illustrate exemplary embodiments of the invention and, together with their description, serve to explain the invention and do not constitute an undue limitation thereof. In the drawings: Figure 1 This is a flowchart illustrating a method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to the present invention. Figure 2 This is the initial grid used in the embodiments of the present invention; Figure 3 for Figure 2 A schematic diagram illustrating the parameter changes of sampling points during angle-changing operations in the medium grid. Figure 4 for Figure 2 Mesh based Figure 3 The modified mesh after the content has undergone an angle-shifting operation; Figure 5 for Figure 2 Modified mesh after mirroring the XZ plane in the medium mesh; Figure 6 for Figure 2Modified mesh after mirroring the YZ plane in the medium mesh; Figure 7 This is a schematic diagram illustrating the node order adjustment of a three-dimensional eight-node mesh unit in an embodiment of the present invention. Detailed Implementation

[0029] The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not all of them. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative effort are within the scope of protection of the present invention.

[0030] In the description of this invention, it should be understood that the terms "center", "lateral", "longitudinal", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", etc., indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings. They are only for the convenience of describing this invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on this invention.

[0031] The terms "first," "second," and "third" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Therefore, a feature defined as "first," "second," or "third" may explicitly or implicitly include one or more of that feature.

[0032] In the description of this invention, it should be noted that, unless otherwise explicitly specified and limited, the terms "installation," "connection," and "joining" should be interpreted broadly. For example, they can refer to a fixed connection, a detachable connection, or an integral connection; they can refer to a direct connection or an indirect connection through an intermediate medium; and they can refer to the internal communication between two components. Those skilled in the art can understand the specific meaning of the above terms in this invention based on the specific circumstances.

[0033] like Figure 1-4 As shown in the illustrative embodiment of a method for rapidly modifying a three-dimensional finite element model of a tire with complex tread patterns according to the present invention, the method includes at least the following: S1. Using the tire circumference as the X-axis, the tire width direction as the Y-axis, and the tire radial direction as the Z-axis, perform three-dimensional meshing on the tire tread model to obtain the initial mesh. The center point O of the initial mesh on the XY plane is the origin of the X-axis and Y-axis coordinates. S2, Variable Angle: The initial grid is set with multiple sampling points along the Y-axis, with two sampling points located at the two ends of the initial grid and the other sampling points located at the turning points of the pattern angle of the initial grid; a variable angle value τ is set for each sampling point. Calculate the X-axis change for each sampling bit. x, based on two adjacent sample bits x is calculated using a definite integral algorithm or an interpolation algorithm to determine the change in the X-axis of each grid cell between these two sampling points. The formula for calculating x is:

[0034]

[0035] Where y is the y-axis coordinate. x is the change along the X-axis. The curvature of the XY plane; Based on the calculated changes in the X-axis of the nodes, adjust the X-axis coordinates of each node in the grid; S3, Output the modified mesh.

[0036] The angle variation method described above can achieve local or overall angle adjustment of the finite element model of the tire tread pattern; the method eliminates the steps of adjusting the 3D model of the tread pattern and reconstructing the finite element mesh, which helps to shorten the tread pattern development cycle and reduce the labor cost of tread pattern development.

[0037] For tread patterns that cycle periodically along the tire circumference, the above operations can be performed on the finite element model of a single pitch. Then, the mesh output by S3 can be rotated and arrayed around the tire axis using the rotation function in the finite element software to generate a complete tire tread mesh, further improving the modification speed of the finite element model. Here, "pitch" refers to the length of a complete tread pattern cycle element along the tire circumference.

[0038] In some embodiments, a pattern unit is defined between two adjacent sampling positions, and multiple adjacent pattern units form a continuous pattern. Before calculating the change in the X-axis of the node in S2, all sampling positions in all pattern units involved in the continuous pattern are sequentially numbered along the Y-axis. Each sampling position... x superimposed on the previous sample bit x obtain Based on two adjacent sampling bits The change in the X-axis of each grid cell between these two sampling points is calculated using a definite integral algorithm or an interpolation algorithm.

[0039] By following the above steps, we can avoid the problem of misalignment between adjacent pattern units in a continuous pattern after the angle of some pattern units is adjusted by changing the angle. This avoids the impact of local angle adjustment on the overall continuity of the pattern.

[0040] In some embodiments, for all types of tire tread patterns, the entire initial mesh can be considered as a single continuous pattern to achieve unidirectional integration along the positive or negative Y-axis. For tire tread patterns symmetrical about the XZ plane, in addition to unidirectional integration, the entire initial mesh can also be considered as two continuous patterns symmetrical about the XZ plane. The two continuous patterns are integrated along the Y-axis towards point O, or the two continuous patterns are integrated along the Y-axis away from point O to achieve bidirectional integration.

[0041] Furthermore, all initial meshes for tire tread patterns are constructed using one-way integration to simplify the calculation process and avoid the complexity introduced by two-way integration.

[0042] In some embodiments, after obtaining the change in the X-axis of a node in S2, the change in the X-axis of the node at Y=0 is used as... Subtract the X-axis coordinate of each node The corrected node X-axis coordinates are obtained to ensure that the change in the X-axis at Y=0 of the corrected mesh is 0, so that the X-axis mesh coordinates at Y=0 can be used as the coordinate system correction reference for the modified mesh.

[0043] In some embodiments, the definite integral algorithm in S2 can be any one of Simpson's integral, trapezoidal rule integral, Gaussian integral, or adaptive step size control algorithm; the interpolation algorithm can be any one of piecewise quadratic interpolation or polynomial interpolation. Piecewise quadratic interpolation and polynomial interpolation can be considered low-precision versions of Simpson's integral. Using Simpson's integral to calculate the X-axis change of nodes in S2 is preferable, as it helps ensure the continuity of the integration path's values, achieves smooth transitions between mesh nodes, and reduces model modification time from hours to minutes.

[0044] In some embodiments, multiple sampling points are set at the turning points of a single pattern angle to avoid curvature discontinuities caused by definite integral algorithms or interpolation algorithms.

[0045] The aforementioned Y-axis range ensures that this embodiment can modify the angle of tire tread patterns for all width specifications. Operations where the X-axis offset and τ are both zero (excluding the user data) prevent abrupt changes or breaks in the finite element mesh boundaries during the angle adjustment operation in S2, thus improving the stability and fault tolerance of this embodiment.

[0046] Furthermore, the Y-axis range is from -10000mm to 10000mm.

[0047] In some embodiments, the Y-axis interval in S1 applies to tire tread patterns of all width specifications; for example, the Y-axis interval is from negative infinity to positive infinity. The initial mesh occupies the Y-axis interval of user data, and the X-axis offset and τ outside the user data are both zero.

[0048] like Figure 1 and Figure 5-7 As shown, in some other embodiments, a method for rapidly modifying a three-dimensional finite element model of a tire with complex tread patterns includes at least the following: S1. Using the tire circumference as the X-axis, the tire width direction as the Y-axis, and the tire radial direction as the Z-axis, perform three-dimensional meshing on the tire tread model to obtain the initial mesh. The center point O of the initial mesh on the XY plane is the origin of the X-axis and Y-axis coordinates. S2, Mirror: Mirror the initial mesh about the XZ or YZ plane to generate a mirrored mesh; update the normal information of the mirrored mesh to make the determinant of the Jacobian matrix of the mesh elements positive; S3, Output the modified mesh.

[0049] The mirroring method described above can achieve a complete mirroring of the finite element model of the tire tread pattern. This method eliminates the need to adjust the 3D tread pattern model and reconstruct the finite element mesh, which helps to shorten the tread pattern development cycle and reduce the labor costs of tread pattern development.

[0050] For tread patterns that cycle periodically along the tire circumference, the above operations can be performed on the finite element model of a single pitch. Then, the mesh output by S3 can be rotated and arrayed around the tire axis using the rotation function in the finite element software to generate a complete tire tread mesh, further improving the modification speed of the finite element model. Here, "pitch" refers to the length of a complete tread pattern cycle element along the tire circumference.

[0051] In some embodiments, in S2, mirroring about the XZ plane changes the node coordinates of all mesh cells from (x, y, z) to (x, -y, z); mirroring about the YZ plane changes the node coordinates of all mesh cells from (x, y, z) to (-x, y, z).

[0052] In some embodiments, updating the normal information in S2 is achieved using finite element software through any of the following methods: flipping the normal direction of the mirrored element; displaying and storing the normal directions of the elements in the initial mesh and updating the normal directions of the mirrored element accordingly; swapping the node output order of the top and bottom surfaces in the mirrored element; and adjusting the node numbers corresponding to the node coordinates. Among these, swapping the node output order of the top and bottom surfaces in the mirrored element is optimal, as it does not change the mesh topology, only adjusts the node order, and avoids negative impacts on mesh quality.

[0053] In some embodiments, the initial mesh is divided using three-dimensional six-node elements or three-dimensional eight-node elements in S1.

[0054] The following is combined Figure 1-4 A preferred embodiment of the angle-changing operation in S2 of the method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to the present invention will be described: S1. Establish a coordinate system with the tire circumferential direction as the X-axis, the tire width direction as the Y-axis, and the tire radial direction as the Z-axis. The Y-axis range is -10000mm to 10000mm. Based on this coordinate system, use three-dimensional eight-node elements to perform three-dimensional mesh generation on the single-pitch tire tread pattern 3D model to obtain the initial mesh, such as... Figure 2 As shown; the grid cell node numbering order follows the right-hand rule, as follows: Figure 3 As shown; the initial mesh is located at point O, the center of the XY plane, which is the origin of the X and Y axes. S2, Angle Variable: The initial grid is set with four sampling points along the Y-axis, namely a, b, c, and d, where a and d are... Located at both ends of the initial grid, b and c are located at the turning points of the tread angle of the initial grid. A variable angle value τ is set for each sampling position. The τ corresponding to a is -Q, and Q is a positive value. The τ corresponding to b and c is 0, and the τ corresponding to d is Q. That is, only the tread angle at the two tire shoulders is adjusted, and the tread angle at other positions remains unchanged.

[0055] The initial grid occupies the Y-axis interval of user data, and the X-axis offset and τ outside the user data are both zero.

[0056] Calculate the X-axis change for each sampling bit. x, The formula for calculating x is:

[0057]

[0058] Where y is the y-axis coordinate. x is the change along the X-axis. The curvature of the XY plane; The area between two adjacent sampling points is a pattern unit. The entire initial grid is considered as a continuous pattern. All sampling points in all pattern units are numbered sequentially along the Y-axis. Each sampling point... x superimposed on the previous sample bit x obtain Based on two adjacent sampling bits The Simpson integral is used to calculate the change in the X-axis of each grid cell between the two sampling points. A mathematical mapping between the pattern angle and the node displacement is established through one-way integration. The Simpson algorithm is used to ensure the continuity of the numerical values ​​of the integration path, achieving a smooth transition of grid nodes. This can reduce the model modification time from hours to minutes.

[0059] The change in the X-axis at the node where Y=0 is taken as... Subtract the X-axis coordinate of each node This yields the corrected X-axis coordinates of the nodes.

[0060] S3, Output the modified mesh.

[0061] This embodiment can achieve local or overall angle adjustment of the finite element model of tire tread pattern, skipping the steps of adjusting the 3D model of the tread pattern and reconstructing the finite element mesh. It dynamically links the adjustment of design parameters with mesh modification, directly modifies the existing mesh, significantly shortens the tread pattern development cycle, and reduces the labor cost of tread pattern development.

[0062] The modified mesh, obtained using the above method, has a Jacobian value greater than 0.9, indicating high mesh quality. The Jacobian is used in mesh quality checks to evaluate the difference between the actual and ideal shape of a cell. Specifically, the Jacobian is the ratio of the minimum to the maximum value of the Jacobian matrix at each integration point within the cell.

[0063] Furthermore, multiple sampling points can be set at the angle transition of a single pattern to improve the curvature continuity of the overall mesh after local angle changes. For example... Figure 2 If multiple sampling points are set at point b, the impact of the angle change at point a on the curvature continuity at the junction of segments ab and bc can be reduced.

[0064] The following is combined Figure 1 and Figure 5-7 A preferred embodiment of the mirroring operation in S2 of the method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to the present invention will be described below: S1. Establish a coordinate system with the tire circumferential direction as the X-axis, the tire width direction as the Y-axis, and the tire radial direction as the Z-axis. The Y-axis range is -10000mm to 10000mm. Based on this coordinate system, use three-dimensional eight-node elements to perform three-dimensional mesh generation on the single-pitch tire tread pattern 3D model to obtain the initial mesh, such as... Figure 2 As shown; the direction of the element normal is determined by following the right-hand rule based on the node numbering order of the mesh elements. Figure 7 N1, N2, N3, N4, N5, N6, N7, and N8 are the eight nodes of the grid cell. Figure 7 In This is a schematic diagram of the nodes of the original mesh. The bottom and top faces are numbered counterclockwise, and the normals of the mesh elements point upwards. The center point O of the initial mesh on the XY plane is the origin of the X and Y axes.

[0065] S2, Mirroring: The initial mesh occupies the Y-axis interval of user data, and the X-axis offset outside the user data is zero. The initial mesh is mirrored about the XZ or YZ plane to generate a mirrored mesh; mirroring about the XZ plane changes the node coordinates of all mesh elements from (x, y, z) to (x, -y, z), such as... Figure 7 In As shown. Mirroring the YZ plane changes the node coordinates of all mesh elements from (x, y, z) to (-x, y, z), as... Figure 7 In As shown.

[0066] Due to mirroring operations, Figure 7 middle and The numbering order of the bottom and top surfaces is changed to clockwise, therefore... Figure 7 middle and The mesh elements are output with node information in the order of N1, N2, N3, N4, N5, N6, N7, N8. This can lead to problems such as downward-facing element normals, disordered stress and strain directions, and impaired material property definitions or accuracy of contact analysis simulation results. To ensure that the Jacobian matrix determinant is positive and meets the requirements of finite element analysis, the element orientation needs to be reversed so that the element normal direction points from the bottom surface to the top surface, thus automatically correcting the element normal direction errors caused by mirroring operations.

[0067] Based on the right-hand rule, when calculating the normal direction of a mesh element, if the Z-axis coordinate of the bottom surface is greater than the Z-axis coordinate of the top surface, it will cause the element normal direction to reverse. Therefore, the finite element software outputs the values ​​in the order of N5, N6, N7, N8, N1, N2, N3, and N4. Figure 7 middle and By analyzing the information of the grid cell nodes, the cell normal direction can be reversed, making the cell normal direction upward.

[0068] S3, Output the modified mesh.

[0069] This embodiment also enables the complete mirroring of the finite element model of the tire tread pattern, skipping the steps of adjusting the 3D tread pattern model and reconstructing the finite element mesh. It dynamically links design parameter adjustments with mesh modifications, directly modifying the existing mesh, significantly shortening the tread pattern development cycle and reducing the labor costs of tread pattern development. This embodiment retains the seed point distribution of the mirrored mesh, adjusting the node order while keeping the node coordinates unchanged, ensuring the consistency of the element normals and avoiding mesh distortion caused by mirroring.

[0070] The modified mesh, obtained using the above method, has a Jacobian value greater than 0.9, indicating high mesh quality. The Jacobian is used in mesh quality checks to evaluate the difference between the actual and ideal shape of a cell. Specifically, the Jacobian is the ratio of the minimum to the maximum value of the Jacobian matrix at each integration point within the cell.

[0071] Through the description of several embodiments of the method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to the present invention, it can be seen that the embodiments of the method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to the present invention have at least one or more of the following advantages: 1. This invention enables local or overall angle adjustment of the finite element model of tire tread patterns, and also allows for the complete mirroring of the finite element model of tire tread patterns. It skips the steps of adjusting the 3D tread pattern model and reconstructing the finite element mesh, dynamically linking design parameter adjustments with mesh modifications, directly modifying the existing mesh, significantly shortening the tread pattern development cycle and reducing the labor costs of tread pattern development. This invention is particularly suitable for scenarios requiring frequent iterations in the early stages of tread pattern development, such as tire tread angle optimization, groove style adjustment, and symmetry generation, significantly reducing finite element mesh reconstruction time.

[0072] 3. This invention preserves the seed point distribution of the mirrored mesh and adjusts the node order while keeping the node coordinates unchanged, thus ensuring the normal consistency of the elements and avoiding mesh distortion caused by mirroring.

[0073] 4. The modified mesh based on this invention has a Jacobian value greater than 0.9, and the modified mesh has higher quality.

[0074] Finally, it should be noted that the various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. The same or similar parts between the various embodiments can be referred to each other.

[0075] The above embodiments are only used to illustrate the technical solutions of the present invention and not to limit them; although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications can still be made to the specific implementation of the present invention or equivalent substitutions can be made to some technical features without departing from the spirit of the technical solutions of the present invention, and all such modifications and substitutions should be covered within the scope of the technical solutions claimed in the present invention.

Claims

1. A method for rapidly modifying a three-dimensional finite element model of a tire with complex tread patterns, characterized in that, include: S1. Using the tire circumference as the X-axis, the tire width direction as the Y-axis, and the tire radial direction as the Z-axis, perform three-dimensional meshing on the tire tread model to obtain the initial mesh. The center point O of the initial mesh on the XY plane is the origin of the X-axis and Y-axis coordinates. S2. Shift or mirror the initial mesh to generate a modified mesh: Variable angle: The initial grid has multiple sampling points along the Y-axis, with two sampling points located at the two ends of the initial grid and the other sampling points located at the angle transitions of the initial grid pattern; a variable angle value τ is assigned to each sampling point. Calculate the X-axis change for each sampling bit. x, based on two adjacent sample bits x is calculated using a definite integral algorithm or an interpolation algorithm to determine the change in the X-axis of each grid cell between these two sampling points. The formula for calculating x is: Where y is the y-axis coordinate. x is the change along the X-axis. The curvature of the XY plane; Based on the calculated changes in the X-axis of the nodes, adjust the X-axis coordinates of each node in the grid; Mirroring: Mirrors the initial mesh about the XZ or YZ plane to generate a mirrored mesh; updates the normal information of the mirrored mesh to make the determinant of the Jacobian matrix of the mesh elements positive; S3, Output the modified mesh.

2. The method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to claim 1, characterized in that, Two adjacent sampling positions constitute a pattern unit. Multiple adjacent pattern units form a continuous pattern. Before calculating the change in the X-axis of the node in S2, all sampling positions in all pattern units involved in the continuous pattern are sequentially numbered along the Y-axis. Each sampling position... x superimposed on the previous sample bit x obtain Based on two adjacent sampling bits The change in the X-axis of each grid cell between these two sampling points is calculated using a definite integral algorithm or an interpolation algorithm.

3. The method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to claim 2, characterized in that, The entire initial mesh can be viewed as a single continuous pattern, or as two continuous patterns symmetrical about the XZ plane.

4. The method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to claim 3, characterized in that, After obtaining the change in the X-axis of the node in S2, the change in the X-axis of the node at Y=0 is taken as... Subtract the X-axis coordinate of each node This yields the corrected X-axis coordinates of the nodes.

5. A method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to any one of claims 1-4, characterized in that, In S2, the definite integral algorithm type can be any one of Simpson's integral, trapezoidal rule integral, Gaussian integral, or adaptive step size control algorithm; the interpolation algorithm type can be any one of piecewise quadratic interpolation or polynomial interpolation.

6. A method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to any one of claims 1-4, characterized in that, In S2, mirroring the XZ plane changes the node coordinates of all mesh elements from (x, y, z) to (x, -y, z); mirroring the YZ plane changes the node coordinates of all mesh elements from (x, y, z) to (-x, y, z).

7. A method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to any one of claims 1-4, characterized in that, In S2, updating the normal information is achieved using finite element software through any of the following methods: flipping the normal direction of the mirrored element; displaying and storing the normal direction of the elements in the initial mesh and updating the normal direction of the mirrored element accordingly; swapping the output order of the nodes on the top and bottom surfaces of the mirrored element; and adjusting the node number corresponding to the node coordinates.

8. A method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to any one of claims 1-4, characterized in that, Multiple sampling points are set at the angle transition of a single pattern.

9. A method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to any one of claims 1-4, characterized in that, The Y-axis interval in S1 is applicable to tire tread patterns of all width specifications. The Y-axis interval occupied by the initial mesh is user data, and the X-axis offset and τ outside the user data are both zero.

10. A method for rapidly modifying a three-dimensional finite element model of a complex tire tread pattern according to any one of claims 1-4, characterized in that, In S1, the initial mesh is generated using three-dimensional six-node elements or three-dimensional eight-node elements.

Images