Material ablation interface grid reconstruction method, device, equipment and medium

Abstract

The invention provides a material ablation interface grid reconstruction method, device and equipment and a medium, and the method comprises the steps: defining the normal direction of a material matrix phase plane, carrying out horizontal expansion pretreatment on a curved surface space quadrilateral grid of the initial shape of a single cylindrical fiber ablation interface to obtain a plane grid, and according to the normal direction of the material matrix phase plane, carrying out horizontal expansion pretreatment on the curved surface space quadrilateral grid; the ablation interface grid points are updated through interpolation, and ablation interface horizontal grids are obtained; and based on a single cylindrical fiber ablation interface horizontal grid, performing grid distortion judgment and reconstruction updating on each layer of grid points. And distortion judgment and reconstruction updating of all grid layers and grid reconstruction of each single cylindrical fiber are completed according to the steps, so that material ablation interface grid reconstruction is realized. The method has the advantages of being high in grid distortion judgment efficiency, high in distortion reconstruction processing speed and small in numerical calculation amount, and is suitable for optimizing and reconstructing the interface grid in the ablation interface evolution tracking process based on the Lagrange 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 environments such as aerodynamic heating, different 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 to the simulation of the thermal response of the material and the evaluation of the thermal performance. The simulation of the evolution of the ablation shape belongs to the tracking problem of complex moving interfaces. The simulation methods that can be used for such problems include the Euler method and the Lagrangian method. The Euler method generally requires that the grid division of the calculation area is fine enough in the process of tracking 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, and regenerates the ablation surface with all the newly generated grid points at the new moment. The grid size requirement is relatively small, but it is necessary to solve the grid distortion and non-physical distortion of the grid caused by the evolution of the interface grid.

[0003] For the Lagrangian method of ablation interface evolution, if the ablation interface mesh is distorted or even distorted, the ablation calculation process will not converge or even fail in simulation calculation, and adaptive meshing or mesh reconstruction is required. Among them, the adaptive meshing method is generally used to deal with relatively small ablation deformation displacement, but once 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 a good simulation effect for typical problems, it involves the processing of volume grids, the amount of numerical simulation calculations is relatively large, and the implementation process of the method is also relatively complicated.

[0005] In view of this, there is an urgent need to provide a volume mesh analysis and reconstruction method with high material mesh distortion judgment efficiency 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 material ablation interface mesh reconstruction method, comprising the following steps:

[0008] Step S1, defining the normal direction of the material matrix phase plane, performing 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, and updating 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 matrix plane grid lines, thereby obtaining a horizontal grid 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 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 less numbers, and distribute them at equal distances, and finally move the grid point closest to the intersection point to cover the intersection point, so as 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 realizing the mesh reconstruction of the material ablation interface.

[0011] One or more embodiments of the present specification provide a material ablation interface mesh 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 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;

[0013] A grid distortion reconstruction module is used to judge and reconstruct each layer of grid points based on the horizontal grid of the ablation interface of a single cylindrical fiber obtained by the horizontal 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 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 less numbers, and distribute them at equal distances, and finally move the grid points closest to the intersection points to cover the intersection points, so as 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 the present specification provide a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the material ablation interface mesh reconstruction method as described above when executing the computer program.

[0016] One or more embodiments of the present 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 invention provides a material ablation interface grid reconstruction method, device, equipment and medium, which have the advantages of being a non-physical distortion judgment and reconstruction method of ablation morphology grid under the synergistic effect of strategies such as "horizontal preprocessing of ablation interface grid" and "intersection judgment and inductive update of in-plane distortion grid". 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 a grid line segment intersection judgment problem in the plane. The grid coordinate update is realized by combining the inductive method to quickly judge and correct the intersection type in the plane. 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 optimizing and reconstructing the interface grid in the process of tracking the evolution of the ablation interface 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 drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. 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 paying 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 an initial material matrix provided for one or more embodiments of this specification;

[0021] Figure 3A schematic diagram of an ablation interface grid obtained after horizontal processing provided by one or more embodiments of this specification;

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

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

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

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

[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] Fig. 9 A schematic diagram of the structure of a computer device provided for one or more embodiments of this specification. DETAILED DESCRIPTION

[0028] In order to enable those skilled in the art to 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 creative work should fall within the scope of protection of the present invention.

[0029] The present invention is described in detail below in conjunction with specific implementation methods and the accompanying drawings.

[0030] Method Embodiment

[0031] According to an embodiment of the present invention, a material ablation interface mesh reconstruction method is provided, such as Figure 1 As shown, it 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 space quadrilateral grid of the initial shape of the ablation interface of a single cylindrical fiber to obtain a plane grid, and updating 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 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.

[0033] Step S2, based on the horizontal grid of the ablation interface of a single cylindrical fiber, the grid distortion judgment and reconstruction update are performed on each layer of grid points; 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 turn according to the traversal direction, and the number of grid points between the first and last traversed two 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 numbers are moved to between the two grid points with less numbers, and are distributed at equal distances, and finally the grid points closest to the intersection are moved and cover the intersection to complete the grid reconstruction of the corresponding grid layer.

[0034] Step S3, according to step S2, complete the distortion judgment and reconstruction update of all grid layers, and according to step S1 to step S2, complete the grid reconstruction of each single cylindrical fiber, so as to achieve the reconstruction of the material ablation interface grid. Thus, the coordinates of the spatial grid points of the material ablation interface after reconstruction are obtained, and finally the ablation shape grid information without physical distortion and distortion at the current calculation time is obtained.

[0035] The material ablation interface grid reconstruction method provided in this embodiment is 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 "intersection judgment and inductive update of the distorted grid in the plane". By performing a specific preprocessing of "horizontal expansion" on the spatial grid of the ablation interface of the composite material, the judgment of the distortion of the ablation interface grid is converted into a problem of grid line segment intersection judgment in the plane. The grid coordinate update is realized by combining the inductive method to quickly judge and correct the intersection type in the plane. 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 optimizing and reconstructing the interface grid in the process of tracking the evolution of the ablation interface based on the Lagrangian method, which is beneficial to improving the comprehensive analysis efficiency of the simulation problem of the evolution of the composite material ablation interface grid.

[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, it is a schematic diagram of the initial material matrix provided in this embodiment. It can be seen that before material ablation occurs, the fiber 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 direction of the defined material matrix phase plane, the initial shape of the ablation interface of the single cylindrical fiber and the curved surface space quadrilateral grid are horizontalized to obtain n r ×n s A 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 grids, 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, with 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 interpolated and updated based on the coordinate values ​​of the grid points with the corresponding numbers in the current column and the 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, determine the mesh layer that needs to be reconstructed, refer to Figure 4 As shown, it is a schematic diagram of the horizontal grid after the initial fiber is expanded provided in this embodiment, so after step S1 updates the ablation interface grid points by interpolation or before step S13, the following steps are also 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 , then it 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 crossing. The grid after the "horizontalization" processing in step S1 will also result in grid line crossing. 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 are crossed.

[0044] In a specific embodiment, reference Figure 5-6 As shown, Figure 5 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, specifically, based on the horizontal mesh of the ablation interface of a single cylindrical fiber, the mesh distortion judgment and reconstruction and updating 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, determine the grid lines P in turn i P i+1 and the 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 there is a mesh segment intersection, the four mesh points corresponding to the crossed mesh 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 the calculation points ② and ③ are horizontally moved and evenly distributed between the two grid points corresponding to the calculation points ① and ④, and the grid point closest to Pc is moved to cover the Pc point.

[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 point 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 FIG. 1 is a schematic diagram of a single cylindrical fiber ablation interface reconstruction provided by an embodiment, wherein Figure 7 (a) is the initial schematic diagram of the ablation interface. Figure 7 (b) Schematic diagram of the reconstruction of the ablation interface.

[0053] The method provided in this embodiment is applicable to, but not limited to, the interface mesh reconstruction of 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 "horizontalization", so that the judgment of whether the ablation interface grid is distorted only needs to be performed within the plane grid points after the "horizontalization" of each layer, which solves the problem of difficulty in judging the distortion of the composite material ablation interface grid under large ablation deformation conditions. Practice has shown that it has high accuracy and analysis efficiency.

[0056] (2) Compared with other mesh reconstruction methods, the method of this embodiment attributes the mesh 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 situations and reconstructs them separately. It has the characteristics of high mesh distortion judgment efficiency and fast distortion reconstruction processing speed, which is beneficial to improving the comprehensive simulation efficiency of the mesh 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 As shown, it is a block diagram of the material ablation interface mesh reconstruction device provided in this embodiment. The complex material ablation interface mesh reconstruction device according to the 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 and reconstruct each layer of grid points based on the horizontal grid of the ablation interface of a single cylindrical fiber 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 less numbers, and distribute them at equal distances, and 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.

[0061] The material ablation interface mesh reconstruction module 30 is used to complete the mesh reconstruction of each single cylindrical fiber through the horizontal processing module 10 and the mesh distortion reconstruction module 20, so as to achieve the material ablation interface mesh reconstruction. In this way, the coordinates of the spatial mesh points of the material ablation interface after the reconstruction processing are obtained, and finally the ablation shape mesh information without physical distortion and distortion at the current calculation time is obtained.

[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 judgment of the ablation interface grid distortion is converted into a grid line segment intersection judgment problem in the plane. The grid coordinate update is realized by combining the inductive method to quickly judge and correct the intersection type in the plane. 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 evolution of the ablation interface 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 horizontal 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 mesh of the curved surface according to the normal direction of the defined material matrix phase plane to obtain n r ×n s A 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, with Z k The value is the interpolation target height. For each column (1, 2, 3, ..., n r) The calculation points in the grid points with the same sequence number as the current layer number are updated by interpolation 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 , then it 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 judgment and reconstruction update, which are 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, determine the grid lines P in turn i P i+1 and the 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 there is a mesh segment intersection, the four mesh points corresponding to the crossed mesh 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 the calculation points ② and ③ are horizontally moved and evenly distributed between the two grid points corresponding to the calculation points ① and ④, and the grid point closest to Pc is moved to cover the Pc point.

[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 point 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, which will not be repeated here.

[0077] like Fig. 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-mentioned 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 can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and 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. As an illustration and not limitation, RAM is available in many 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, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the device or system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment. 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 may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[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 by equivalents. These modifications or replacements do not make the essence of the corresponding technical solutions deviate 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 belong to the common knowledge of those skilled in the art.

Claims

1. A material ablation interface mesh 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 space quadrilateral grid of the initial shape of the ablation interface of a single cylindrical fiber to obtain a plane grid, and updating 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 matrix plane grid lines, thereby obtaining a horizontal grid 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 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 less numbers, and distribute them at equal distances, and finally move the grid point closest to the intersection point to cover the intersection point, so as 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 realizing 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 surface space quadrilateral mesh 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, with 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 the 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 grid reconstruction 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, determine the grid lines P in turn i P i+1 Respectively and grid lines Whether there is 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 grid point; Step S23, if there is a mesh segment intersection, the four mesh points corresponding to the crossed mesh lines are numbered in ascending order as calculation point ①, calculation point ②, calculation point ③ and calculation point ④, and the intersection point is numbered as 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 the Pc point; 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 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 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; The grid distortion reconstruction module is used to judge and reconstruct each layer of grid points based on the horizontal grid of the ablation interface of a single cylindrical fiber obtained by the horizontal 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 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 less numbers, and distribute them at equal distances, and finally move the grid points closest to the intersection points to cover the intersection points 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 comprises: Material matrix phase plane definition submodule, used to define the normal direction of the material matrix phase plane; A horizontal processing submodule is used to horizontalize the initial shape of the ablation interface of a single cylindrical fiber and the surface space quadrilateral mesh based on the normal 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, with 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 the 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 horizontal 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 that needs to be reconstructed 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, determine the grid lines P in turn i P i+1 Respectively and grid lines Whether there is 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 grid point; Step S23, if there is a mesh segment intersection, the four mesh points corresponding to the crossed mesh lines are numbered in ascending order as calculation point ①, calculation point ②, calculation point ③ and calculation point ④, and the intersection point is numbered as 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 grid points between the two grid points corresponding to the calculation points ② and ③ are horizontally moved and evenly distributed between the two grid points corresponding to the calculation points ① and ④, and the grid point closest to Pc is moved to cover the Pc point; if N 23 >N 14 , then 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 point 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, characterized in that: 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