A parametric migration method for isogeometric analysis based on discrete geometric mapping

Abstract

The present invention discloses a parameterized migration method for isogeometric analysis based on discrete geometric mapping, which first classifies given feature constraints. Secondly, according to the classification results, if it is impossible to directly move the boundary of the facet where the constraint is located to generate high-quality parameterization, the facet where the constraint is located is divided into multiple sub-facets so that the constraint is located at the boundary of the new sub-facet. Then radial basis function interpolation is performed to move the endpoints of the constraint point or constraint line to the facet boundary. Finally, for the situation where cracks appear on the facet due to feature points or lines, an optimization method is used to optimize the local parameterization quality to ensure that the constraint position remains unchanged, and obtain a high-quality parameterized result that meets the constraint. The present invention optimizes the facets with severe quality degradation, and obtains a parameterized result with higher quality while satisfying the feature constraints, which meets the application of subsequent analysis.

Description

Technical Field

[0001] The invention relates to the technical field of geometric model processing, and in particular to an isogeometric analysis parameterization migration method based on discrete geometric mapping. Background Art

[0002] In 2005, Professor Hughe first proposed isogeometric analysis. Its basic idea is to use the same mathematical language in modeling and analysis, and directly apply the geometric spline model designed in CAD (Computer Aided Design) software to the physical analysis in CAE (Computer Aided Engineering), so as to achieve a seamless combination of CAD and CAE. Whether it is a geometric model designed by CAD software or a model reconstructed by three-dimensional scanning, it usually only has a spline boundary of the region. In order to obtain the physical quantities distributed inside the model, it is necessary to construct a spline representation of the internal area of ​​the model, which is the parameterization problem in isogeometric analysis. Parameterization is an important step that plays a cornerstone role in isogeometric analysis, and its quality greatly affects the accuracy and stability of subsequent analysis.

[0003] At present, there are many computational domain parameterization methods for isogeometric analysis in academia, but these parameterization methods for isogeometric analysis all generate the overall spline representation of the computational domain from the boundary model, and do not consider the situation with internal constraints. In isogeometric analysis, for planar models with feature points and line constraints, when parameterizing, it is hoped that the feature points and lines are located at the corners or boundaries of the patch, so that constraints can be imposed on the features in subsequent isogeometric analysis. How to generate high-quality patches near the feature constraints and ensure that the boundaries of the patches are well attached to the feature constraints so that specific constraints are met in simulation analysis is a key problem that needs to be solved in the field of isogeometric analysis parameterization, which has certain theoretical value and application prospects. Summary of the invention

[0004] In order to solve the above problems, the present invention proposes a parameterized migration method of isogeometric analysis based on discrete geometric mapping. For plane geometric models and internal point features or line features, the present invention provides a high-quality parameterization method that makes the features located at the internal facet boundary, which facilitates the subsequent imposition of constraints at the features.

[0005] In order to solve the above technical problems, the technical solution of the present invention is:

[0006] A parameterized migration method for isogeometric analysis based on discrete geometric mapping includes the following steps:

[0007] Step (1) classifying the given feature constraints;

[0008] Step (2) according to the classification result of step (1), if it is not possible to directly move the boundary of the face where the constraint is located to generate a high-quality parameterization, then split the face where the constraint is located into multiple sub-faces, so that the constraint can be located at the boundary of the new sub-face, which is convenient for applying the constraint;

[0009] Step (3) performing radial basis function (RBF) interpolation to move the constraint point or the endpoint of the constraint line to the boundary of the patch;

[0010] 3-1. Generate a triangular mesh of a plane model;

[0011] 3-2. Perform RBF interpolation movement;

[0012] Step (4) When a feature point or line has a small crack in the patch due to its special position, an optimization method is used to optimize the local parameterization quality to ensure that the constraint position remains unchanged, so as to obtain a high-quality parameterization result that satisfies the constraint;

[0013] 4-1. Construct a quality optimization function that preserves constraints;

[0014] 4-2. Fit the plane model boundary and obtain the final optimization result through interpolation movement.

[0015] Preferably, the step (1) comprises the following sub-steps:

[0016] 1-1. For a given plane model, perform parameterization based on segmentation;

[0017] 1-2. For feature point constraints, the corner points of the patch closest to the feature points on the plane model are used as the reference, and the feature points are divided into three categories according to their positions relative to the reference points;

[0018] 1-3. For feature line constraints, the feature lines are divided into 7 categories based on the position of the feature line endpoints in the patch, taking the patch corner point closest to the feature line on the plane model as the reference.

[0019] Preferably, in step (2), the method for dividing the constrained patch into multiple sub-patches is:

[0020] Uniformly sample the patch where the feature constraint is located, calculate the sampling point closest to the feature point or the first and last endpoints of the feature line, and the node value to be inserted is the node value corresponding to the selected sampling point. Assuming that the order of the patch where the constraint is located in both directions is k, insert the node value until its repetition is k, so that the patch is divided into multiple small patches.

[0021] Preferably, in step 3-1, the method for generating a triangular mesh of a plane model is as follows: for a plane model with known parameterization results, a Delaunay triangulation (a classic triangulation algorithm that can maximize the minimum angle of the mesh triangle after triangulation) is performed to generate a background mesh according to the positions of its internal control points and boundary control polygons.

[0022] Preferably, in step 3-2, the RBF interpolation moving method is:

[0023] According to the constraint classification, the boundary control points of the patch to be moved are determined, and the corresponding m mesh vertices are marked in the background mesh. According to the positions of the m points to be moved and the target point, the offset of the m points is determined. Based on the RBF interpolation method, the positions of other control points inside the plane geometric model are determined to obtain the initial parameterization that satisfies the feature constraints.

[0024] Preferably, the step (4) comprises the following sub-steps:

[0025] 4-1. For geometric models with known parametric results and internal feature point constraints or line constraints, the above process can generate parametric results that accurately meet the constraints. However, for some point and line constraints with special locations, although the above method can avoid the decline of internal parametric quality as much as possible while meeting the constraints, the patches near the features may have tiny cracks due to segmentation and discretization interpolation movement. The quality of parameterization can be improved by optimizing the following objective function:

[0026]

[0027] in,

[0028] Q l (P i,j )=∫ Ω ||s u || 2 +||s v || 2 dΩ

[0029] represents the parametric line length function,

[0030] Q u (P i,j )=∫ Ω ||s uu || 2 +2||s uv || 2 +||s vv || 2 dΩ

[0031] represents the uniformity function of the parameter line,

[0032] Q o (P i,j )=∫ Ω (s u ·s v ) 2 dΩ

[0033] represents the orthogonality function,

[0034]

[0035] represents the skewness function,

[0036] Q a (P i,j )=∫ Ω |s u ,s v | 2 dΩ

[0037] represents the area function,

[0038]

[0039] represents the eccentricity function. l ,w u ,w o ,w s ,w a ,w e is a non-negative weight value that determines the influence of the related function. i,j are the control points in the region C formed by all the thinned and split patches in the initial parameterization, and s is the patch in region C. v represents the first-order partial derivative of the patch with respect to the v direction, s u represents the first-order partial derivative of the patch about the u direction, s uu represents the second-order partial derivative of the patch about the u direction, s uv represents the mixed partial derivatives of the patch about the u and v directions, s vv represents the second-order partial derivative of the patch with respect to the v direction, s v Represents the first-order partial derivative of the patch with respect to the v direction.

[0040] 4-2. Since the boundary of the geometric model is not fixed during quality optimization, the shape of the model changes. To solve this problem, after optimization, the boundary curve basis function of C is used to approximate the curve of this part before, and the final optimization result can be obtained by interpolation movement.

[0041] The present invention has the following characteristics and beneficial effects:

[0042] The present invention proposes a constrained parameterization method for a plane model for isogeometric analysis, which can obtain a parameterization result that satisfies the internal point constraint or line constraint under the premise of the known initial parameterization of the input model. For feature constraints close to the endpoints of the facet boundary, the coordinates of the facet boundary control points can be moved based on RBF interpolation to obtain a parameterization result that satisfies the constraint. For some point constraints or line constraints with special positions, the facets where the point constraints or line constraints are located are subdivided using heavy node interpolation before RBF interpolation, and then subsequent operations are performed to ensure the quality of the facets to a certain extent. The present invention includes a local area quality optimization method, which optimizes the facets with severe quality degradation, and obtains a parameterization result with higher quality while satisfying the feature constraints, which meets the application of subsequent analysis. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] In order to more clearly illustrate the embodiments of the present invention 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 of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0044] Figure 1 is a flowchart of steps of an embodiment of the present invention;

[0045] Figure 2 is p located inside the model, p s Located at the circle centered at p p Parameterized results inside the circle;

[0046] Figure 3 is that p is on the boundary of the model, or p s Not located in the circle centered on p p Parameterized results inside the circle;

[0047] Figure 4 is the parameterization result corresponding to classification 3;

[0048] Figure 5 is the parameterization result corresponding to classification 4;

[0049] Figure 6 is the parameterization result corresponding to classification 5;

[0050] Figure 7 is the parameterization result corresponding to classification 6;

[0051] Figure 8 It is the result of direct RBF interpolation movement of the constraint point located in the center of the patch;

[0052] Fig. 9It is the parameterized result that satisfies the constraints after local refinement;

[0053] Fig.10 is the parameterization result when the constraint line is on the diagonal;

[0054] Fig.11 It is the parameterized result that satisfies the constraints after local quality optimization. DETAILED DESCRIPTION

[0055] It should be noted that, in the absence of conflict, the embodiments of the present invention and the features in the embodiments may be combined with each other.

[0056] In the description of the present invention, it should be understood that the terms "center", "longitudinal", "lateral", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside" and the like indicate positions or positional relationships based on the positions or positional relationships shown in the accompanying drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as limiting the present invention. In addition, the terms "first", "second", and the like are only used for descriptive purposes, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, features defined as "first", "second", and the like may explicitly or implicitly include one or more of the features. In the description of the present invention, unless otherwise specified, "multiple" means two or more.

[0057] In the description of the present invention, it should be noted that, unless otherwise clearly specified and limited, the terms "installed", "connected", and "connected" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be indirectly connected through an intermediate medium, or it can be the internal communication of two components. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood by specific circumstances.

[0058] The present invention provides an isogeometric parameterization method for a plane model with feature constraints, such as Figure 1 As shown, the following steps are included:

[0059] Step (1) classifies the given feature constraints.

[0060] Specific:

[0061] 1-1. For a given plane model, use the method in "Constructing IGA-suitable planarparameterization from complex CAD boundary by domain partition and global / local optimization" to generate partition-based parameterization;

[0062] 1-2. For feature point constraints, define the distance from the feature point p on the plane model. s The nearest patch corner point is p, and the arc lengths of the patch boundary curves C1 and C2 on both sides of p are dis1 and dis2 respectively. The circle with p as the center and min(dis1 / 2, dis2 / 2) as the radius is the circle Cir p According to the distance p s The nearest corner point p is divided into two categories based on whether it is on the outer boundary of the input model. If the corner point p is on the outer boundary of the model, moving the corner point p to meet the feature constraint will change the overall shape of the model, which does not meet the requirements of simulation analysis in actual engineering, and is marked as category 1. s Located on the circle Cir with p as the center p Direct interpolation of moving corner points can obtain high-quality parameterized results while satisfying the constraints. Based on this, when the corner point p is located inside the model, according to p s Is it located in the circle Cir with p as the center? p It is divided into two categories, namely category 2 and category 3. Figure 2 and Figure 3 ;

[0063] 1-3. For the characteristic constraint line, its first and last endpoints are p a ,p b , from p a ,p b The nearest patch corners are p1 and p2, and the boundary curves C on both sides of p1 a1 ,C a2 The arc lengths are dis a1 ,dis a2 , the boundary curves C on both sides of p2 b1 ,C b2 The arc lengths are dis b1 ,dis b2 , with p1 as the center and min(dis a1 / 2,dis a2 / 2) is a circle with a radius Similarly, with p2 as the center and min(dis b1 / 2,dis b2 / 2) is a circle with a radius If the leading and trailing endpoints p of the feature line constraint a ,p b The nearest corner points p1 and p2 are points on the diagonal line of the patch. Directly moving the internal patch boundary of the model to meet the feature constraint will greatly reduce the quality of the patch around the constraint, and it is marked as category 1; if the corner points p1 and p2 closest to the feature line constraint are the two endpoints of the patch boundary, and p1 and p2 are both located on the external boundary of the input model, similar to the feature point constraint, directly moving the patch boundary will change the external shape of the model, and it is marked as category 2; if p1 and p2 are the two endpoints of the patch boundary and one is located at the model boundary and the other is located inside the model, assuming that p1 is located at the model boundary, at this time, determine the constraint point p corresponding to p2 b Is it in the circle with p2 as the center? If it is not in the circle If p1 and p2 are the two endpoints of the patch boundary and both p1 and p2 are inside the model, the model is classified as category 3, otherwise it is classified as category 4. If p1 and p2 are the two endpoints of the patch boundary and both p1 and p2 are inside the model, the model is classified as category 3 according to whether p1 and p2 are respectively inside the corresponding circle. This situation is further divided into three categories, namely category 5, category 6 and category 7. Figure 4-Figure 7 shown.

[0064] Step (2) divides the constrained patch into multiple sub-patches.

[0065] The specific implementation method is as follows:

[0066] Assume that the face where the constraint is located is

[0067]

[0068] The node vectors are U={u0,u1,…,u n+p+1} and V={v0,v1,…,v m+q+1},N i,p (u) represents the i-th p-th spline basis function in the direction of the node vector U, N j,q (v) represents the j-th q-order spline basis function in the direction of the knot vector V, P i,j Represents the control vertices of the patch. The node values ​​to be inserted are Assume that after interpolation in the u direction, the new node vector is in:

[0069]

[0070] Then the new node vector U1 can define n+2 p-order B-spline basis functions i=0,1,…,n+2, the initial p×q B-spline surface can be obtained by the new B-spline basis function and the new control vertices after refinement express, The definition of is as follows:

[0071]

[0072] in

[0073]

[0074] Where u is a parameter, u i 1 is the node value of the node vector U1.

[0075] To split the original face f into multiple small face units, the above process needs to be repeated p times in the u direction and q times in the v direction, so that the face f is broken into 4 small face units f i ’ ,i=1,…,4.

[0076] Step (3) performs RBF interpolation to move the constraint point or the endpoint of the constraint line to the patch boundary, such as Figure 8 shown.

[0077] 3-1. Generate a triangular mesh of the plane model.

[0078] Specifically, the method for generating a triangular mesh for a plane model with known parameterized results is: based on the Triangle library (a classic mesh processing program library), Delaunay triangulation is performed according to the parameterized internal control point positions and boundary control polygons.

[0079] 3-2. Determine the boundary control points of the patch to be moved according to the constraint classification, and mark the corresponding m mesh vertices in the background mesh. Determine the offset of the m points according to the positions of the m points to be moved and the target point. Determine the positions of other control points inside the plane geometric model based on the RBF interpolation method, and obtain the initial parameterization that satisfies the feature constraints.

[0080] Specifically, RBF consists of m radial basis functions q(t i ), i=1,…,m, t i Represents the distance between the point to be moved and the target point. RBF can provide a smooth interpolation function for the entire discrete space:

[0081]

[0082] where a i, i=1,…,m and c0, c1, c2, c3 are the (m+4) unknown quantities that need to be calculated, thus forming a certain affine transformation. RBF interpolation can calculate the unknown quantities in the above equation according to the equation system formed by the following formula, and determine the offset of other grid points except m known points.

[0083]

[0084] Among them A x ,A y ,A z is an (m+4)-dimensional column vector (a1, a2, …, a m ,c0,c1,c2,c3). G is a (m+4)×(m+4) matrix, as shown below:

[0085]

[0086] where g ij =g(||u i ||),i=1,…,m, The present invention uses the following function to determine g ij :

[0087]

[0088] Usually k=1, so (m+4) coefficients to be determined can be determined, and the RBF interpolation function can be determined. Based on the determined RBF interpolation function, the offsets of other grid points inside the model except the m known points are calculated to determine the initial parameterization results that meet the feature constraints.

[0089] Step (4) For the patches with serious quality degradation, an optimization method is used to perform local parameterized quality optimization while keeping the constraint position unchanged, and a high-quality parameterized result that satisfies the constraints is obtained, such as Fig. 9 shown.

[0090] 4-1. For geometric models with known parametric results and internal feature point constraints or line constraints, the above process can generate parametric results that accurately meet the constraints. However, for some point and line constraints with special locations, although the above method can avoid the decline of internal parametric quality as much as possible while meeting the constraints, the patches near the features may have tiny cracks due to segmentation and discretization interpolation movement. The quality of parameterization can be improved by optimizing the following objective function:

[0091]

[0092] in,

[0093] Q l (Pi,j )=∫ Ω ||s u || 2 +||s v || 2 dΩ

[0094] represents the parametric line length function,

[0095] Q u (P i,j )=∫ Ω ||s uu || 2 +2||s uv || 2 +||s vv || 2 dΩ

[0096] represents the uniformity function of the parameter line,

[0097] Q o (P i,j )=∫ Ω (s u ·s v ) 2 dΩ

[0098] represents the orthogonality function,

[0099]

[0100] represents the skewness function,

[0101] Q a (P i,j )=∫ Ω |s u ,s v | 2 dΩ

[0102] represents the area function,

[0103]

[0104] represents the eccentricity function. l ,w u ,w o ,w s ,w a ,w e is a non-negative weight value that determines the influence of the related function. i,j is the control point in the region C formed by all the thinned and split patches in the initial parameterization, s is the patch in region C, and s v represents the first-order partial derivative of the patch with respect to the v direction, s uurepresents the second-order partial derivative of the patch about the u direction, s uv represents the mixed partial derivatives of the patch about the u and v directions, s vv represents the second-order partial derivative of the patch with respect to the v direction, s v Represents the first-order partial derivative of the patch with respect to the v direction.

[0105] 4-2. The boundary of region C is recorded as The boundary of plane model A is Fixed constraint line C s and the boundary of region C The control points on the curve segment are optimized, and the positions of other control points inside are optimized, and the parameterization quality of the model is greatly improved. However, since the boundary of A is not fixed, the shape of the model changes. To solve this problem, the present invention uses The curve basis function of the curve is approximately optimized before the part of the curve is optimized, and the final optimization result is obtained by interpolation movement, such as Fig.10 and Fig.11 shown.

[0106] The embodiments of the present invention are described in detail above with reference to the accompanying drawings, but the present invention is not limited to the described embodiments. For those skilled in the art, various changes, modifications, substitutions and variations of these embodiments including components are made without departing from the principles and spirit of the present invention, and still fall within the scope of protection of the present invention.

Claims

1. A parameterized migration method for isogeometric analysis based on discrete geometric mapping, characterized in that: The steps include: Step (1) classifies the given feature constraints; the specific process is as follows: 1-1. For a given plane model, perform parameterization based on segmentation; 1-2. For feature point constraints, the corner point of the patch closest to the feature point on the plane model is used as the reference, and the feature points are divided into three categories according to their positions relative to the reference points; specifically: For feature point constraints, define the distance from the feature point p on the plane model. s The nearest patch corner point is p. The arc lengths of the patch boundary curves C1 and C2 on both sides of p are dis1 and dis2 respectively. The circle with p as the center and min(dis1 / 2, dis2 / 2) as the radius is the circle Cir. p , according to the distance p s Whether the nearest corner point p is on the outer boundary of the input model is divided into two categories. If the corner point p is on the outer boundary of the model, moving the corner point p satisfies the feature constraint, which changes the overall shape of the model and does not meet the requirements of simulation analysis in actual engineering. It is marked as category 1. s Located on the circle Cir with p as the center p Direct interpolation moves the corner point to obtain a parameterized result while satisfying the constraints. Based on this, the corner point p is located inside the model. s Is it located in the circle Cir with p as the center? p It is subdivided into two categories, namely category 2 and category 3; 1-3. For feature line constraints, the corner point of the plane model closest to the feature line is used as the reference, and the feature lines are divided into 7 categories according to the position of the feature line endpoints in the patch, specifically: For the feature constraint line, its starting and ending points are p a ,p b , from p a ,p b The nearest patch corners are p1 and p2, and the boundary curves C on both sides of p1 a1 ,C a2 The arc lengths are dis a1 ,dis a2 , the boundary curves C on both sides of p2 b1 ,C b2 The arc lengths are dis b1 ,dis b2 , with p1 as the center and min(dis a1 / 2,dis a2 / 2) is a circle with a radius Similarly, with p2 as the center and min(dis b1 / 2,dis b2 / 2) is a circle with a radius If the leading and trailing endpoints p of the feature line constraint a ,p b The nearest corner points p1 and p2 are points on the diagonal line of the patch. Directly moving the internal patch boundary of the model to meet the feature constraint will greatly reduce the quality of the patch around the constraint, and it is marked as category 1; if the corner points p1 and p2 closest to the feature line constraint are the two endpoints of the patch boundary, and p1 and p2 are both located on the external boundary of the input model, similar to the feature point constraint, directly moving the patch boundary will change the external shape of the model, and it is marked as category 2; if p1 and p2 are the two endpoints of the patch boundary and one is located at the model boundary and the other is located inside the model, assuming that p1 is located at the model boundary, at this time, determine the constraint point p corresponding to p2 b Is it in the circle with p2 as the center? If it is not in the circle If p1 and p2 are the two endpoints of the patch boundary and both p1 and p2 are inside the model, the model is classified as category 3, otherwise it is classified as category 4. If p1 and p2 are the two endpoints of the patch boundary and both p1 and p2 are inside the model, the model is classified as category 3 according to whether p1 and p2 are respectively inside the corresponding circle. This situation is further divided into three categories, namely category 5, category 6 and category 7; Step (2) according to the classification result, if it is not possible to directly move the boundary of the face where the constraint is located to generate parameterization, then split the face where the constraint is located into multiple sub-faces, so that the constraint is located at the boundary of the new sub-face; Step (3) performing radial basis function (RBF) interpolation to move the constraint point or the endpoint of the constraint line to the boundary of the patch; Step (4) In the case where cracks appear in the facet due to feature points or lines, an optimization method is used to optimize the local parameterization quality to ensure that the constraint position remains unchanged, and a parameterization result that satisfies the constraint is obtained. The specific process is as follows: 4-1. Construct a quality optimization function that preserves constraints; The quality of parameterization is improved by optimizing the following objective function: in, Q l (P i,j )=∫ Ω ||s u || 2 +||s v || 2 dΩ represents the parametric line length function, Q u (P i,j )=∫ Ω ||s uu || 2 +2||s uv || 2 +||s vv || 2 dΩ represents the uniformity function of the parameter line, Q o (P i,j )=∫ Ω (s u ·s v ) 2 dΩ represents the orthogonality function, represents the skewness function, Q a (P i,j )=∫ Ω |s u ,s v | 2 dΩ represents the area function, represents the eccentricity function; w l ,w u ,w o ,w s ,w a ,w e is a non-negative weight value that determines the influence of the related function; P i,j is the control point in the region C formed by all the thinned and split patches in the initial parameterization, s is the patch in region C, and s v represents the first-order partial derivative of the patch with respect to the v direction, s u represents the first-order partial derivative of the patch about the u direction, s uu represents the second-order partial derivative of the patch about the u direction, s uv represents the mixed partial derivatives of the patch in the u and v directions, s vv Represents the second-order partial derivative of the patch with respect to the v direction; 4-2. Fit the plane model boundary and obtain the optimization result through interpolation movement; After the optimization, the boundary curve basis function of C is used to optimize the curve of the previous part, and the optimization result is obtained by interpolation movement.

2. The isogeometric analysis parameterization migration method based on discrete geometric mapping according to claim 1 is characterized in that: The specific process of dividing the face where the constraint is located into multiple sub-faces in step (2) is as follows: The patch where the feature constraint is located is uniformly sampled, and the sampling point closest to the feature point or the first and last endpoints of the feature line is calculated. The node value to be inserted is the node value corresponding to the selected sampling point. Assuming that the order of the patch where the constraint is located in both directions is k, the patch is divided into multiple small patches by inserting node values ​​until its repetition degree is k.

3. The isogeometric analysis parameterization migration method based on discrete geometric mapping according to claim 2 is characterized in that: The specific process of step (3) is as follows: 3-1. Generate a triangular mesh of a plane model; For a plane model with known parameterized results, Delaunay triangulation is performed based on the positions of its internal control points and boundary control polygons to generate a background mesh. 3-2. Perform RBF interpolation movement; According to the constraint classification, the boundary control points of the patch to be moved are determined, and the corresponding m mesh vertices are marked in the background mesh; according to the positions of the m points to be moved and the target point, the offsets of the m points are determined; based on the RBF interpolation method, the positions of other control points inside the plane geometric model are determined to obtain the initial parameterization that satisfies the feature constraints.

Images