A material ablation interface grid reconstruction method, device, equipment and medium

Abstract

The present application provides a material ablation interface mesh reconstruction method, device, equipment and medium, the method includes defining the normal of the material matrix phase plane, horizontally unfolding the curved space quadrilateral mesh of the initial shape of the ablation interface of a single cylindrical fiber to obtain a plane mesh, updating the ablation interface mesh points by interpolation according to the normal of the material matrix phase plane, and obtaining the ablation interface horizontal mesh; based on the horizontal mesh of the ablation interface of a single cylindrical fiber, performing mesh distortion judgment and reconstruction update on each layer of mesh points. According to the above steps, the distortion judgment and reconstruction update of all mesh layers and the mesh reconstruction of each single cylindrical fiber are completed, thereby realizing the material ablation interface mesh reconstruction. The present application has the characteristics of high mesh distortion judgment efficiency, fast distortion reconstruction processing speed and small amount of numerical calculation, and is suitable for processing the optimization and reconstruction of the interface mesh in the process of ablation interface evolution tracking based on the Lagrangian method.

Description

Technical Field

[0001] The present invention relates to the field of simulation technology, and in particular to a material ablation interface grid reconstruction method, device, equipment and medium. Background Art

[0002] Under thermal boundary conditions such as aerodynamic heating, varying degrees of ablation will occur on the surface of the material. The accuracy of the simulation of the evolution of the ablation shape is crucial for the simulation of the material's thermal response and the evaluation of its thermal performance. The simulation of the evolution of the ablation shape belongs to the problem of tracking a complex moving interface. The simulation methods that can be used for such problems include the Euler method and the Lagrangian method. The Euler method generally requires a sufficiently fine mesh division in the calculation area during the tracking of the ablation surface, so the amount of calculation required is relatively large. The Lagrangian method calculates the normal velocity of each grid point on the ablation surface and generates new grid points. At the new moment, all newly generated grid points are used to regenerate the ablation surface. The mesh size requirement is relatively small, but it is necessary to solve the mesh distortion and non-physical mesh distortion problems that may arise due to the evolution of the interface mesh.

[0003] For Lagrangian methods of ablation interface evolution, if the ablation interface mesh is distorted or even distorted, the ablation calculation process will not converge or even fail, requiring the use of adaptive meshing or mesh reconstruction. Adaptive meshing is generally used for relatively small ablation deformation displacements, but if large ablation deformation occurs, mesh reconstruction is usually required.

[0004] Most of the current grid reconstruction methods are highly targeted. Among them, the grid reconstruction method for ablation problems mainly analyzes and reconstructs the discrete grids of hexahedral and tetrahedral units of thermal protection structures. Although it has good simulation effects for typical problems, due to the processing of volume grids, the numerical simulation calculation amount is relatively large and the implementation process of the method is relatively complicated.

[0005] In view of this, there is an urgent need to provide a volume mesh analysis and reconstruction method with high efficiency in judging material mesh distortion and fast distortion reconstruction processing speed. Summary of the Invention

[0006] In order to overcome the problems existing in the related art, the present disclosure provides a material ablation interface mesh reconstruction method, device, equipment and medium to solve the technical problems of surface mesh distortion and distortion occurring during the numerical simulation of composite material ablation morphology evolution in the related art.

[0007] One or more embodiments of this specification provide a method for reconstructing a material ablation interface mesh, comprising the following steps:

[0008] Step S1, defining the normal direction of the material matrix phase plane, performing horizontal expansion preprocessing on the curved spatial quadrilateral mesh of the initial shape of the ablation interface of a single cylindrical fiber to obtain a plane mesh, and updating the ablation interface mesh points by interpolation according to the normal direction of the material matrix phase plane, so that all radial mesh lines of the ablation interface mesh are parallel to the mesh lines of the matrix plane, thereby obtaining a horizontal mesh of the ablation interface after horizontal processing, wherein the normal direction of the material matrix phase plane is the z-coordinate direction;

[0009] Step S2: Based on the horizontal grid of the ablation interface of a single cylindrical fiber, perform grid distortion judgment and reconstruction update on each layer of grid points; based on the grid layer, traverse and determine the intersecting line segments in the grid layer, determine the four grid points corresponding to the intersecting line segments in sequence according to the traversal direction, and determine the number of grid points between the first and last traversed grid points, as well as the number of grid points between the other two grid points. Move the grid points on the side with a larger number to between the two grid points with a smaller number and distribute them at equal distances. Finally, move the grid point closest to the intersection point to cover the intersection point to complete the grid reconstruction of the corresponding grid layer; and

[0010] Step S3, completing the distortion judgment and reconstruction update of all mesh layers according to step S2, and completing the mesh reconstruction of each single cylindrical fiber according to steps S1 to S2, thereby achieving the mesh reconstruction of the material ablation interface.

[0011] One or more embodiments of this specification provide a material ablation interface grid reconstruction device, comprising:

[0012] A horizontal processing module is used to define the normal direction of the material matrix phase plane, perform horizontal expansion preprocessing on the initial shape of the ablation interface of a single cylindrical fiber and the curved surface space quadrilateral grid to obtain a plane grid, and update the ablation interface grid points by interpolation according to the normal direction of the material matrix phase plane, so that all radial grid lines of the ablation interface grid are parallel to the grid lines of the matrix plane, thereby obtaining the horizontal grid of the ablation interface after horizontal processing, wherein the normal direction of the material matrix phase plane is the z coordinate direction;

[0013] A grid distortion reconstruction module is used to determine grid distortion and reconstruct and update each layer of grid points based on the horizontal grid of the single cylindrical fiber ablation interface obtained by the horizontalization processing module; based on the grid layer, traverse and determine the intersecting line segments in the grid layer, determine the four grid points corresponding to the intersecting line segments in sequence according to the traversal direction, and determine the number of grid points between the first and last traversed grid points, as well as the number of grid points between two other grid points; move the grid points on the side with a larger number to between the two grid points with a smaller number, and distribute them at equal distances; finally, move the grid point closest to the intersection point to cover the intersection point to complete the grid reconstruction of the corresponding grid layer; and

[0014] The material ablation interface mesh reconstruction module is used to complete the mesh reconstruction of each single cylindrical fiber through the horizontal processing module and the mesh distortion reconstruction module, thereby realizing the material ablation interface mesh reconstruction.

[0015] One or more embodiments of this specification provide a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the material ablation interface mesh reconstruction method described above is implemented.

[0016] One or more embodiments of this specification provide a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, the material ablation interface mesh reconstruction method as described above is implemented.

[0017] The present disclosure provides a material ablation interface grid reconstruction method, device, equipment and medium, which have the advantages of being a method for judging and reconstructing the non-physical distortion of the ablation morphology grid under the synergistic effect of strategies such as "horizontal preprocessing of the ablation interface grid" and "in-plane distortion grid intersection judgment and inductive update". By performing a specific preprocessing of "horizontal expansion" on the spatial grid of the composite material ablation interface, the judgment of the ablation interface grid distortion is converted into an in-plane grid line segment intersection judgment problem. The inductive method is combined to quickly judge and correct the in-plane intersection type to realize grid coordinate update. The method has the characteristics of high grid distortion judgment efficiency, fast distortion reconstruction processing speed and small numerical calculation amount. It is suitable for processing the optimization and reconstruction of the interface grid in the process of tracking the ablation interface evolution based on the Lagrangian method, which is beneficial to improving the comprehensive analysis efficiency of the composite material ablation interface grid evolution simulation problem. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] In order to more clearly illustrate one or more embodiments of this specification or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments recorded in this specification. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0019] Figure 1 A flow chart of a material ablation interface mesh reconstruction method provided for one or more embodiments of this specification;

[0020] Figure 2 A schematic diagram of the initial material matrix provided for one or more embodiments of this specification;

[0021] Figure 3A schematic diagram of the ablation interface grid obtained after horizontal processing according to one or more embodiments of this specification;

[0022] Figure 4 A schematic diagram of a horizontal grid after the initial fibers are expanded according to one or more embodiments of this specification;

[0023] Figure 5 A schematic diagram of line segment intersections between grid points in a grid layer provided for one or more embodiments of this specification;

[0024] Figure 6 A flow chart for reconstructing the material ablation interface mesh provided for one or more embodiments of this specification;

[0025] Figure 7 A schematic diagram of a single cylindrical fiber ablation interface reconstruction provided in one or more embodiments of this specification, wherein Figure 7 (a) is the initial schematic diagram of the ablation interface. Figure 7 (b) Schematic diagram of the ablated interface after reconstruction;

[0026] Figure 8 A block diagram of a material ablation interface mesh reconstruction device provided for one or more embodiments of this specification; and

[0027] Figure 9 A schematic diagram of the structure of a computer device provided in one or more embodiments of this specification. DETAILED DESCRIPTION

[0028] In order to help those skilled in the art better understand the technical solutions in one or more embodiments of this specification, the technical solutions in one or more embodiments of this specification will be clearly and completely described below in conjunction with the drawings in one or more embodiments of this specification. Obviously, the described embodiments are only part of the embodiments of this specification, not all of the embodiments. Based on one or more embodiments of this specification, all other embodiments obtained by ordinary technicians in this field without making creative efforts should fall within the scope of protection of the present invention.

[0029] The present invention will be described in detail below with reference to specific implementation methods and the accompanying drawings.

[0030] Method Example

[0031] According to an embodiment of the present invention, a material ablation interface mesh reconstruction method is provided, such as Figure 1 FIG. 1 is a flow chart of the material ablation interface mesh reconstruction method provided in this embodiment. This embodiment describes the mesh reconstruction process based on the ablation morphology of a single cylindrical fiber. The material ablation interface mesh reconstruction method according to the embodiment of the present invention includes the following steps:

[0032] Step S1, defining the normal direction of the material matrix phase plane, performing horizontal expansion preprocessing on the curved spatial quadrilateral mesh of the initial shape of the ablation interface of a single cylindrical fiber to obtain a plane mesh, and updating the ablation interface mesh points by interpolation according to the normal direction of the material matrix phase plane, so that all radial mesh lines of the ablation interface mesh are parallel to the matrix plane mesh lines, thereby obtaining the horizontal mesh of the ablation interface after horizontal processing, wherein the normal direction of the material matrix phase plane is the z-coordinate direction.

[0033] In step S2, based on the horizontal grid of the ablation interface of a single cylindrical fiber, the grid distortion of each layer of grid points is judged and reconstructed and updated; based on the grid layer, the intersecting line segments in the grid layer are traversed and determined, and the four grid points corresponding to the intersecting line segments are determined in sequence according to the traversal direction. The number of grid points between the first and last traversed grid points, as well as the number of grid points between the other two grid points, are determined, and the grid points on the side with more grid points are moved to between the two grid points with fewer grid points, and are distributed at equal distances. Finally, the grid point closest to the intersection point is moved to cover the intersection point to complete the grid reconstruction of the corresponding grid layer.

[0034] In step S3, distortion determination and reconstruction updates are completed for all mesh layers based on step S2. Furthermore, mesh reconstruction of each individual cylindrical fiber is completed based on steps S1 and S2, thereby reconstructing the mesh of the material ablation interface. This results in the coordinates of the spatial mesh points of the reconstructed material ablation interface, ultimately yielding ablation shape mesh information free of physical distortion and distortion at the current calculation time.

[0035] The material ablation interface mesh reconstruction method provided in this embodiment is a method for judging and reconstructing the non-physical distortion of the ablation morphology mesh under the synergistic effect of strategies such as "horizontal preprocessing of the ablation interface mesh" and "in-plane distortion mesh intersection judgment and inductive update". By performing a specific preprocessing of "horizontal expansion" on the spatial mesh of the composite material ablation interface, the judgment of the ablation interface mesh distortion is converted into an in-plane mesh line segment intersection judgment problem. The inductive method is combined to quickly judge and correct the in-plane intersection type to realize the update of the mesh coordinates. This method has the characteristics of high mesh distortion judgment efficiency, fast distortion reconstruction processing speed and small numerical calculation amount. It is suitable for processing the optimization and reconstruction of the interface mesh in the process of tracking the ablation interface evolution based on the Lagrangian method, which is beneficial to improving the comprehensive analysis efficiency of the composite material ablation interface mesh evolution simulation problem.

[0036] In this embodiment, step S1 specifically includes the following steps:

[0037] Step S11, defining the normal direction of the material matrix phase plane, wherein the normal direction of the material matrix phase plane is the z coordinate direction. Figure 2As shown, this is a schematic diagram of the initial material matrix provided in this embodiment. It can be seen that before material ablation occurs, the fibers and the surrounding matrix will form an initial plane, and the ablation direction of the matrix will always remain parallel to this matrix plane, that is, the material matrix phase plane is fixed and not arbitrary.

[0038] Step S12: Based on the normal of the defined material matrix phase plane, the initial shape of the ablation interface of the single cylindrical fiber and the quadrilateral mesh of the curved surface are horizontalized to obtain n r ×n s Plane grid, where n r is the number of meshes around the fiber, n s is the mesh number along the fiber axis.

[0039] Step S13: Based on the plane grid, reconstruct the z-direction spacing ΔZ according to the set grid j1 (ΔZ j1 >ΔZ1), where ΔZ1 is the Z-direction spacing of the first layer of grid, from the first layer to the last layer (n s layer), take Z k =Z k-1 +ΔZ j1 , k=1,2,3,…,n s , respectively Z k The value is the interpolation target height, for each column (1,2,3,...,n r ) The calculation points with the same number as the current layer in the grid points are updated based on the coordinate values ​​of the grid points with the corresponding numbers in the current column and current layer and the coordinate values ​​of the original grid points, so that all radial grid lines of the ablation interface grid are parallel to the substrate plane, and the ablation interface grid after "horizontalization" is obtained. Figure 3 As shown, it is a schematic diagram of the ablation interface grid obtained after the horizontal processing provided in this embodiment.

[0040] In this embodiment, n obtained in step S1 or step S12 r ×n s The plane mesh of the plane mesh, in which the part that is not ablated in the z direction does not need to be interpolated and updated in the z coordinate. Therefore, in order to reduce unnecessary "interpolation update" actions and save computer resources, it is necessary to first determine the mesh layer that needs to be "interpolated and updated", that is, to determine the mesh layer that needs to be reconstructed. Figure 4 FIG. 1 is a schematic diagram of a horizontal grid after the initial fiber expansion provided in this embodiment. Therefore, after updating the ablation interface grid points by interpolation in step S1 or before step S13, the following steps are further included:

[0041] Step S121: Based on the plane grid, starting from the bottom grid layer of the fiber plane grid, determine the maximum z coordinate z of the grid in the layer. j,maxIs it higher than the z coordinate value z corresponding to the current position of the matrix phase plane? matrix , if z j,max >z matrix , is the grid layer that needs to be reconstructed, and the grid layer that needs to be reconstructed is determined layer by layer from bottom to top.

[0042] In a specific embodiment, when determining the grid layers to be reconstructed layer by layer from bottom to top, the grid layers that need to be reconstructed are marked in turn. t (t=1,2,3,……,n z ).

[0043] In one embodiment, the material ablation grid evolution process is likely to result in grid line intersections. The grid after the "horizontalization" processing in step S1 will also result in grid line intersections. Therefore, after the z coordinates of all grid points are processed through the above steps S11-S13, the x, y coordinates of each layer of grid need to be judged and reconstructed and updated. The distortion judgment is the question of whether the line segments between the grid points in the plane intersect.

[0044] In one embodiment, reference Figure 5-6 As shown, Figure 5 This is a schematic diagram of the intersection of line segments between grid points in the grid layer provided in this embodiment. Figure 6 The material ablation interface mesh reconstruction flow chart provided in this embodiment is as follows: Based on the horizontal mesh of the ablation interface of a single cylindrical fiber, the mesh distortion judgment and reconstruction update steps for each layer of mesh points are as follows:

[0045] Step S21: Based on any grid layer, all grid point coordinates in the plane are numbered in counterclockwise order.

[0046] Step S22: Count the number of grid points from 1 to n. r -3, judge the grid lines P in turn i P i+1 and grid lines P i+ 1P i+2 , P i+2 P i+3 ,…, Is there an intersection? If there is an intersection, go to step S23. If there is no intersection, perform grid line intersection judgment based on the next (i+1) grid point.

[0047] Step S23: If a grid segment crosses, the four grid points corresponding to the cross lines are numbered in ascending order as calculation point ①, calculation point ②, calculation point ③ and calculation point ④, and the intersection point is numbered Pc.

[0048] Step S24: Count the number of grid points between calculation points ②, ③ and calculation points ①, ④ as N. 23 and N 14 ; For example, refer to Figure 5 , grid lines P4P5 and P8P9 intersect, then grid points P4, P5, P8, P9 are numbered as calculation point ①, calculation point ②, calculation point ③ and calculation point ④ in sequence, and the number of grid points between calculation points ② and ③ is determined to be 2 (N 23 =2), the number of grid points between calculation points ① and ④ is 0 (N 14 =0).

[0049] Step S25, if N 23 ≤N 14 , all the grid points between the two grid points corresponding to calculation points ② and ③ are moved horizontally and distributed equally between the two grid points corresponding to calculation points ① and ④, and the grid point closest to Pc is moved to cover Pc.

[0050] Step S26, on the contrary, if N 23 >N 14 , all grid points between the two grid points corresponding to calculation points ① and ④ are horizontally moved and evenly distributed between the two grid points corresponding to calculation points ② and ③, and the grid point closest to Pc is moved to Pc.

[0051] Step S27, repeat steps S23 to S26 until the grid points are completed Distortion judgment and reconstruction processing of grid points.

[0052] In another embodiment, if step S21 numbers all the grid point coordinates in the plane in a clockwise order, it can be determined without a doubt that all the judgments in steps S25 to S26 are opposite to those in the above embodiment. Figure 7 As shown, it is a schematic diagram of the reconstruction of the ablation interface of a single cylindrical fiber provided in the embodiment, wherein Figure 7 (a) is the initial schematic diagram of the ablation interface. Figure 7 (b) Schematic diagram of the ablated interface after reconstruction.

[0053] The method provided in this embodiment is applicable to, but not limited to, the reconstruction of interface meshes during the ablation morphology evolution process of various fiber composite materials, particle-reinforced composite materials at the particle scale, fiber scale, fiber bundle scale, and macro-model scale.

[0054] Compared with the prior art, the method of this embodiment has the following beneficial effects:

[0055] (1) Compared with other grid reconstruction methods, the method of this embodiment pre-processes the spatial quadrilateral grid of the composite material ablation interface by "horizontalizing" it, so that the judgment of whether the ablation interface grid is distorted only needs to be performed within the plane grid points after "horizontalization" of each layer. This solves the problem of difficulty in judging the distortion of the composite material ablation interface grid under large ablation deformation conditions, and practice has shown that it has high accuracy and analysis efficiency.

[0056] (2) Compared with other grid reconstruction methods, the method of this embodiment attributes the grid distortion problem that may occur in the evolution process of the composite material ablation interface to the intersection problem of line segments in the plane, and further summarizes the intersection problem into two cases for reconstruction processing respectively. It has the characteristics of high efficiency in grid distortion judgment and fast distortion reconstruction processing speed, which is beneficial to improving the comprehensive simulation efficiency of the grid evolution process of the composite material ablation interface.

[0057] Device embodiment

[0058] According to an embodiment of the present invention, a material ablation interface grid reconstruction device is provided, such as Figure 8 FIG. 1 is a block diagram of a material ablation interface mesh reconstruction device provided in this embodiment. The complex material ablation interface mesh reconstruction device according to an embodiment of the present invention includes:

[0059] The horizontal processing module 10 is used to define the normal direction of the material matrix phase plane, and to perform horizontal expansion preprocessing on the curved space quadrilateral grid of the initial shape of the ablation interface of a single cylindrical fiber to obtain a plane grid. According to the normal direction of the material matrix phase plane, the ablation interface grid points are updated by interpolation so that all radial grid lines of the ablation interface grid are parallel to the matrix plane grid lines, thereby obtaining the horizontal grid of the ablation interface after horizontal processing, wherein the normal direction of the material matrix phase plane is the z coordinate direction.

[0060] The grid distortion reconstruction module 20 is used to judge the grid distortion and reconstruct and update the grid points of each layer based on the horizontal grid of the single cylindrical fiber ablation interface obtained by the horizontal processing module 10; based on the grid layer, traverse and determine the intersecting line segments in the grid layer, determine the four grid points corresponding to the intersecting line segments in turn according to the traversal direction, and determine the number of grid points between the first and last traversed grid points, as well as the number of grid points between the other two grid points, move the grid points on the side with more numbers to between the two grid points with fewer numbers, and distribute them at equal distances, and finally move the grid points closest to the intersection point to cover the intersection point to complete the grid reconstruction of the corresponding grid layer, and.

[0061] The material ablation interface mesh reconstruction module 30 is used to complete the mesh reconstruction of each individual cylindrical fiber through the leveling processing module 10 and the mesh distortion reconstruction module 20, thereby achieving material ablation interface mesh reconstruction. This obtains the coordinates of the spatial mesh points of the reconstructed material ablation interface, ultimately obtaining ablation shape mesh information without physical distortion or distortion at the current calculation time.

[0062] The material ablation interface grid reconstruction device provided in this embodiment is based on the non-physical distortion judgment and reconstruction method of the ablation morphology grid under the synergistic effect of strategies such as "horizontal preprocessing of the ablation interface grid" and "in-plane distortion grid intersection judgment and inductive update". By performing a specific preprocessing of "horizontal expansion" on the spatial grid of the composite material ablation interface, the ablation interface grid distortion judgment is converted into an in-plane grid line segment intersection judgment problem. The inductive method is combined to quickly judge and correct the in-plane intersection type to realize grid coordinate update. This method has the characteristics of high grid distortion judgment efficiency, fast distortion reconstruction processing speed and small numerical calculation amount. It is suitable for processing the optimization and reconstruction of the interface grid in the process of ablation interface evolution tracking based on the Lagrangian method, which is beneficial to improving the comprehensive analysis efficiency of the composite material ablation interface grid evolution simulation problem.

[0063] In this embodiment, the horizontalization processing module 10 includes:

[0064] The material matrix phase plane definition submodule is used to define the normal direction of the material matrix phase plane, wherein the normal direction of the material matrix phase plane is the z coordinate direction.

[0065] The horizontal processing submodule is used to horizontalize the initial shape of the ablation interface of a single cylindrical fiber and the quadrilateral grid of the curved surface according to the normal direction of the defined material matrix phase plane to obtain n r ×n s Plane grid, where n r is the number of meshes around the fiber, n s is the mesh number along the fiber axis.

[0066] The first grid reconstruction submodule is used to reconstruct the z-direction spacing ΔZ based on the plane grid according to the set grid j1 (ΔZ j1 >ΔZ1), from the first layer to the last layer (n s layer), take Z k =Z k-1 +ΔZ j1 , respectively Z k The value is the interpolation target height, for each column (1,2,3,...,n r) The calculation points with the same sequence number as the current layer number in the grid points are interpolated and updated based on the coordinate values ​​of the grid points with corresponding numbers in the current column and the current layer, so that all radial grid lines of the ablation interface grid are parallel to the substrate plane, and the ablation interface grid after "horizontalization" is obtained.

[0067] In this embodiment, the horizontal processing module 10 further includes a grid layer screening submodule for determining the maximum z coordinate z of the grid in the layer starting from the bottom grid layer of the fiber plane grid. j,max Is it higher than the z coordinate value z corresponding to the current position of the matrix phase plane? matrix , if z j,max >z matrix , is the grid layer that needs to be reconstructed, and the grid layer that needs to be reconstructed is determined layer by layer from bottom to top.

[0068] In this embodiment, the mesh distortion reconstruction module 20 is configured to perform the following steps to implement mesh distortion determination and reconstruction update, as follows:

[0069] Step S21: Based on any grid layer, all grid point coordinates in the plane are numbered in counterclockwise order.

[0070] Step S22: Count the number of grid points from 1 to n. r -3, judge the grid lines P in turn i P i+1 and grid lines P i+ 1P i+2 , P i+2 P i+3 ,…, Is there an intersection? If there is an intersection, go to step S23. If there is no intersection, perform grid line intersection judgment based on the next (i+1) grid point.

[0071] Step S23: If a grid segment crosses, the four grid points corresponding to the cross lines are numbered in ascending order as calculation point ①, calculation point ②, calculation point ③ and calculation point ④, and the intersection point is numbered Pc.

[0072] Step S24: Count the number of grid points between calculation points ②, ③ and calculation points ①, ④ as N. 23 and N 14 ; For example, refer to Figure 5 , grid lines P4P5 and P8P9 intersect, then grid points P4, P5, P8, P9 are numbered as calculation point ①, calculation point ②, calculation point ③ and calculation point ④ in sequence, and the number of grid points between calculation points ② and ③ is determined to be 2 (N 23 =2), the number of grid points between calculation points ① and ④ is 0 (N 14 =0).

[0073] Step S25, if N 23 ≤N 14 , all the grid points between the two grid points corresponding to calculation points ② and ③ are moved horizontally and distributed equally between the two grid points corresponding to calculation points ① and ④, and the grid point closest to Pc is moved to cover Pc.

[0074] Step S26, on the contrary, if N 23 >N 14 , all grid points between the two grid points corresponding to calculation points ① and ④ are horizontally moved and evenly distributed between the two grid points corresponding to calculation points ② and ③, and the grid point closest to Pc is moved to Pc.

[0075] Step S27, repeat steps S23 to S26 until the grid points are completed Distortion judgment and reconstruction processing of grid points.

[0076] The embodiment of the present invention is an apparatus embodiment corresponding to the above-mentioned method embodiment. The specific operations of the processing steps of each module can be understood by referring to the description of the method embodiment, and will not be repeated here.

[0077] like Figure 9 As shown, the present invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the material ablation interface mesh reconstruction method in the above embodiment is implemented.

[0078] The present invention also provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the material ablation interface mesh reconstruction method in the above embodiment is implemented, or when the computer program is executed by a processor, the material ablation interface mesh reconstruction method in the above embodiment is implemented.

[0079] Those skilled in the art will appreciate that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0080] Each embodiment in this specification is described in a progressive manner. The same or similar parts between the embodiments can be referred to each other. Each embodiment focuses on the differences from other embodiments. In particular, for the device or system embodiments, since they are basically similar to the method embodiments, the description is relatively simple. For the relevant parts, refer to the partial description of the method embodiments. The device and system embodiments described above are merely schematic, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the scheme of this embodiment. A person of ordinary skill in the art can understand and implement it without making any creative efforts.

[0081] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the above embodiments, or replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the scope of the technical solutions of the embodiments of the present invention, and the contents not described in detail in the specification of the present invention are common knowledge to those skilled in the art.

Claims

1. A material ablation interface grid reconstruction method, characterized in that The following steps are involved: Step S1, defining the normal direction of the material matrix phase plane, performing horizontal expansion preprocessing on the curved spatial quadrilateral mesh of the initial shape of the ablation interface of a single cylindrical fiber to obtain a plane mesh, and updating the ablation interface mesh points by interpolation according to the normal direction of the material matrix phase plane, so that all radial mesh lines of the ablation interface mesh are parallel to the mesh lines of the matrix plane, thereby obtaining a horizontal mesh of the ablation interface after horizontal processing, wherein the normal direction of the material matrix phase plane is the z-coordinate direction; Step S2: Based on the horizontal grid of the ablation interface of a single cylindrical fiber, perform grid distortion judgment and reconstruction update on each layer of grid points; based on the grid layer, traverse and determine the intersecting line segments in the grid layer, determine the four grid points corresponding to the intersecting line segments in sequence according to the traversal direction, and determine the number of grid points between the first and last traversed grid points, as well as the number of grid points between the other two grid points. Move the grid points on the side with a larger number to between the two grid points with a smaller number and distribute them at equal distances. Finally, move the grid point closest to the intersection point to cover the intersection point to complete the grid reconstruction of the corresponding grid layer; and Step S3, completing the distortion judgment and reconstruction update of all mesh layers according to step S2, and completing the mesh reconstruction of each single cylindrical fiber according to steps S1 to S2, thereby achieving the mesh reconstruction of the material ablation interface.

2. The material ablation interface mesh reconstruction method according to claim 1, characterized in that: Furthermore, step S1 specifically includes the following steps: The normal direction of the defined material matrix phase plane; The initial shape of the ablation interface of a single cylindrical fiber and the quadrilateral mesh of the curved surface are horizontalized based on the normal direction of the defined material matrix phase plane; Based on the plane grid, reconstruct the z-direction spacing ΔZ according to the set grid j1 , starting from the first layer to the last layer, take Z k =Z k-1 +ΔZ j1 , respectively Z k The value is the interpolation target height. For the calculation points in each column of grid points with the same sequence number as the current layer number, the coordinate values ​​of the grid points with corresponding numbers in the current column and the current layer are interpolated and updated, so that all radial grid lines of the ablation interface grid are parallel to the substrate plane, and the ablation interface grid after horizontal processing is obtained.

3. The material ablation interface mesh reconstruction method according to claim 1, characterized in that: The step of updating the ablation interface grid points by interpolation also includes determining the grid layer to be reconstructed, and the specific steps are as follows: Based on the plane grid, starting from the bottom grid layer of the fiber plane grid, determine whether the maximum z coordinate of the grid in the layer is higher than the z coordinate value corresponding to the matrix phase plane position at the current moment. If it is greater, it is the grid layer that needs to be reconstructed. The grid layer for reconstructing the grid is determined layer by layer from bottom to top.

4. The material ablation interface mesh reconstruction method according to any one of claims 1 to 3, characterized in that: Furthermore, step S2 specifically includes the following steps: Step S21: Based on any grid layer, all grid point coordinates in the plane are numbered in counterclockwise order. Step S22: Count the number of grid points from 1 to n. r -3, judge the grid lines P in turn i P i+1 Respectively and grid lines Is there an intersection? If so, go to step S23. If not, perform a grid line intersection determination based on the next grid point. Step S23: If a grid segment crosses, the four grid points corresponding to the cross are numbered in ascending order as calculation point ①, calculation point ②, calculation point ③, and calculation point ④, and the intersection point is numbered Pc. Step S24: Count the number of grid points between calculation points ②, ③ and calculation points ①, ④ as N. 23 and N 14 ; Step S25, if N 23 ≤N 14 , then all the grid points between the two grid points corresponding to the calculation points ② and ③ are moved and distributed equally between the two grid points corresponding to the calculation points ① and ④, and the grid point closest to Pc is moved to cover Pc point; if N 23 >N 14 , then all the grid points between the two grid points corresponding to calculation points ① and ④ are moved and distributed equally between the two grid points corresponding to calculation points ② and ③, and the grid point closest to Pc is moved to point Pc; Step S26, repeat steps S23 to S25 until the grid points are completed Distortion judgment and reconstruction processing of grid points.

5. A material ablation interface grid reconstruction device, characterized in that: include: A horizontal processing module is used to define the normal direction of the material matrix phase plane, perform horizontal expansion preprocessing on the initial shape of the ablation interface of a single cylindrical fiber and the curved surface space quadrilateral grid to obtain a plane grid, and update the ablation interface grid points by interpolation according to the normal direction of the material matrix phase plane, so that all radial grid lines of the ablation interface grid are parallel to the grid lines of the matrix plane, thereby obtaining the horizontal grid of the ablation interface after horizontal processing, wherein the normal direction of the material matrix phase plane is the z coordinate direction; The grid distortion reconstruction module is used to judge the grid distortion and reconstruct and update the grid points of each layer based on the horizontal grid of the single cylindrical fiber ablation interface obtained by the horizontal processing module. Based on the grid layer, the module traverses and determines the intersecting line segments in the grid layer, determines the four grid points corresponding to the intersecting line segments in sequence according to the traversal direction, and determines the number of grid points between the first and last traversed grid points, as well as the number of grid points between the other two grid points. The module moves the grid points on the side with a larger number to between the two grid points with a smaller number and distributes them at equal distances. Finally, the grid point closest to the intersection is moved to cover the intersection to complete the grid reconstruction of the corresponding grid layer. as well as The material ablation interface mesh reconstruction module is used to complete the mesh reconstruction of each single cylindrical fiber through the horizontal processing module and the mesh distortion reconstruction module, thereby realizing the material ablation interface mesh reconstruction.

6. The material ablation interface grid reconstruction device according to claim 5, characterized in that: The horizontal processing module includes: Material matrix phase plane definition submodule, used to define the normal direction of the material matrix phase plane; The horizontal processing submodule is used to horizontalize the initial shape of the ablation interface of a single cylindrical fiber and the quadrilateral mesh of the curved surface space based on the normal direction of the defined material matrix phase plane; The first grid reconstruction submodule is used to reconstruct the z-direction spacing ΔZ based on the plane grid according to the set grid j1 , starting from the first layer to the last layer, take Z k =Z k-1 +ΔZ j1 , respectively Z k The value is the interpolation target height. For the calculation points in each column of grid points with the same sequence number as the current layer number, the coordinate values ​​of the grid points with corresponding numbers in the current column and the current layer are interpolated and updated, so that all radial grid lines of the ablation interface grid are parallel to the substrate plane, and the ablation interface grid after horizontal processing is obtained.

7. The material ablation interface grid reconstruction device according to claim 5, characterized in that: The horizontalization processing module also includes a grid layer screening submodule, which is used to determine whether the maximum z coordinate of the grid in the layer is higher than the z coordinate value corresponding to the matrix phase plane position at the current moment based on the plane grid, starting from the bottom grid layer of the fiber plane grid. If it is greater, it is a grid layer that needs to be reconstructed, and the grid layer for grid reconstruction is determined layer by layer from bottom to top.

8. The material ablation interface grid reconstruction device according to any one of claims 5 to 7, characterized in that: The mesh distortion reconstruction module is configured to perform the following steps: Step S21: Based on any grid layer, all grid point coordinates in the plane are numbered in counterclockwise order. Step S22: Count the number of grid points from 1 to n. r -3, judge the grid lines P in turn i P i+1 Respectively and grid lines Is there an intersection? If so, go to step S23. If not, perform a grid line intersection determination based on the next grid point. Step S23: If a grid segment crosses, the four grid points corresponding to the cross are numbered in ascending order as calculation point ①, calculation point ②, calculation point ③, and calculation point ④, and the intersection point is numbered Pc. Step S24: Count the number of grid points between calculation points ②, ③ and calculation points ①, ④ as N. 23 and N 14 ; Step S25, if N 23 ≤N 14 , all grid points between the two grid points corresponding to calculation points ② and ③ are horizontally moved and distributed equally between the two grid points corresponding to calculation points ① and ④, and the grid point closest to Pc is moved to cover Pc point; if N 23 >N 14 , all grid points between the two grid points corresponding to calculation points ① and ④ are horizontally moved and distributed equally between the two grid points corresponding to calculation points ② and ③, and the grid point closest to Pc is moved to Pc; Step S26, repeat steps S23 to S25 until the grid points are completed Distortion judgment and reconstruction processing of grid points.

9. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the material ablation interface mesh reconstruction method according to any one of claims 1 to 4 is implemented.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the material ablation interface mesh reconstruction method according to any one of claims 1 to 4 is implemented.

Images