141 results about "Discrete points" patented technology

The invention provides a crack identification method and system based on multi-modal data fusion, electronic equipment and a storage medium, and the method comprises the steps: firstly obtaining multi-modal data of a to-be-detected structure surface, specifically including an image sequence obtained by an acquisition module, full-field stress time sequence data and preset key part discrete point stress time sequence data; performing time synchronization calibration on the multi-modal data, mapping a stress field corresponding to full-field stress time sequence data to a pixel coordinate system of an image sequence through camera calibration parameters, and associating discrete point stress time sequence data to the pixel coordinate system to form a comprehensive data set of unified space coordinates; and finally, inputting the comprehensive data set into a trained multi-modal crack recognition model, and performing model recognition to obtain crack parameters and stress analysis results, thereby realizing crack recognition driven by multi-modal data fusion. Multi-modal data collaborative input overcomes the problems of misinformation, missing report and information loss of a single data source, and improves the accuracy and reliability of crack recognition.

The invention relates to a self-adaptive abrasive jet strengthening method and system for a complex curved surface. The method comprises the steps that a three-dimensional point cloud data set of a workpiece is acquired; analyzing the three-dimensional point cloud data set to obtain a first principal component vector; a section plane is set, an intersection line of the section plane and the three-dimensional point cloud data set is calculated based on the first principal component vector, and the intersection line is an initial path; a normal vector of each discrete point on the initial path is calculated, calculation is carried out based on the normal vectors and the custom viewpoints, newest unit outer normal vectors of the discrete points are obtained, and the newest unit outer normal vectors are all perpendicular to the complex curved surface on the workpiece; and calculating corresponding actual position coordinates of the nozzle relative to each discrete point based on the unit outer normal vectors of the discrete points, and forming a continuous strengthening path based on the actual position coordinates. According to the method and the device, the deviation between the CAD model and the real object can be eliminated, the strengthening path planning of the nozzle has the curved surface self-adaptive capability, and the jet flow energy distribution is uniform.

The invention discloses a model discretization method and device, electronic equipment and a computer readable storage medium. Comprising the steps that a to-be-discretized CAD model is obtained, the to-be-discretized CAD model is analyzed, geometric features of the to-be-discretized CAD model are obtained, and the geometric features at least comprise faces and edges; for any adjacent edge A and edge B, discretizing the edge A and the edge B according to the relationship between the edge A and the edge B, and discretizing the edge A and the edge B respectively to obtain a discrete point column of the edge A and a discrete point column of the edge B; and for any adjacent surface C and surface D, discretizing the surface C and the surface D according to the relationship between the surface C and the surface D, and discretizing the surface C and the surface D respectively to obtain a discrete grid of the surface C and a discrete grid of the surface D. According to the method, when the model is converted into the net-shaped structure, it is guaranteed that selfing of the grids does not occur.

The invention discloses a plane target contour reconstruction method and device suitable for geographic information data processing, which are used for recovering the contour shape of a target geographic element from a disordered discrete point set, and comprises the following steps: obtaining a plane discrete point set representing the target geographic element, the plane discrete point set being disordered coordinate data; the plane discrete point set is processed, the plane discrete point set is structured into a group of line segment sets representing local contour features of the plane discrete point set through iterative splitting and merging operation based on an elliptic model, and the splitting and merging operation is controlled by a linearity threshold value and a disconnected connectivity threshold value which are generated in a self-adaptive mode; processing the line segment set, and integrating the line segment set into closed polygon contour data representing the target geographic element by establishing a spatial index for end points of the line segment set and executing nearest neighbor connection; and outputting closed polygon contour data. According to the method and the device, the accuracy and the automation degree of contour reconstruction and the robustness of an algorithm are improved.

The invention discloses a variable grid method grid automatic division method and device based on vector field guidance. The method comprises the following steps: extracting boundary feature points according to geometric features of an input model boundary so as to construct a feature polygon of a model, and dispersing the feature polygon into a triangular mesh model; constructing a vector field based on a Laplacian operator by using the triangular mesh model; starting from the geometric center of the model, setting an initial quadrilateral grid as an initial paving layer, paving the quadrilateral grid outwards layer by layer, and if the paving layer exceeds the boundary of the model, stopping paving; when each layer of quadrilateral grid deforms, uniform transition is achieved through a size transition template; and finally, generating a boundary quadrilateral mesh of the model by adopting a mapping method, mapping discrete points on the boundary of the feature polygon back to the boundary of the input model, and outputting a generated quadrilateral mesh model. According to the method, full-automatic generation of the quadrilateral mesh can be realized, and the method is particularly suitable for mesh generation in finite element simulation pretreatment.

The embodiment of the invention discloses a metal product modeling method and system based on three-dimensional point cloud. The method comprises the steps that original discrete point cloud of a target metal product is acquired, and point cloud regularization processing is carried out; identifying a density mutation region and a normal vector jump region through global scanning, and carrying out region correlation analysis to generate a key structure feature boundary; dividing local point clouds and mining a structure dependency relationship to determine a structure association rule; constructing a structural section connection relation network and creating a topological structure; and aligning the topological structure with the corresponding discrete point cloud segments, inputting the aligned topological structure and the corresponding discrete point cloud segments into a graph neural network to execute reasoning analysis, and realizing structural segment association reconstruction based on a result to obtain a three-dimensional model. Thus, through cooperation of multi-link technical means, the geometrical morphology of the metal product is accurately reproduced, meanwhile, the internal correlation and topological logic of the structure are reserved, and the precision of the three-dimensional model and the structural expression integrity are improved.

The invention discloses an optimization method and system for compressed point cloud data, equipment and a medium. The method comprises the following steps: acquiring first data and second data; extracting local geometric features and global structural features from the first data; adapting a compression distortion type of the first data according to the coding parameter; based on the local geometric features, determining a deviation range of each discrete point in the first data and a point missing region in the first data; restoring the first data by compressing a restoration network corresponding to the distortion type based on the global structure feature, the deviation range and a point missing region in the first data to obtain second data; the first data and the second data are merged to obtain optimized data, distortion types under different compression strengths can be adapted through the coding parameters, and targeted restoration and optimization of the compressed point cloud data are realized in combination with local and global features, so that the geometric fidelity and visual quality of the compressed point cloud data are improved.

The invention relates to the technical field of bolt connection fatigue life prediction, in particular to a high-strength bolt connection fatigue life prediction method and system, and the method comprises the steps: obtaining a bolt entity morphology discrete point cloud, and generating a continuous geometric topology function through non-uniform rational B-spline interpolation reconstruction; discretizing the standard thread parameters to construct a standard finite element model, calculating a normal spatial position difference value based on a continuous geometric topology function, and packaging to generate a geometric defect mapping tensor; according to the converted node coordinates, constructing a real morphology finite element model; solving the nonlinear equilibrium equation by utilizing a penalty function method, and generating a contact stiffness characteristic matrix; calculating an axial load distribution sequence and executing gradient convolution operation, and extracting a local field variable peak value; and calculating an accumulated damage index in combination with a damage evolution equation and outputting a predicted life. By mapping the entity discrete point cloud into a finite element geometric defect tensor and reconstructing a real morphology model, the accurate prediction of the fatigue life of the bolt considering the micromachining deviation influence is realized.

The invention discloses a workpiece contour projection extraction method and system based on ordered point cloud, and is applied to the technical field of three-dimensional precision measurement. The method comprises the following steps of obtaining high-precision ordered point cloud data of a target workpiece and performing preprocessing; geometric analysis is carried out through the three-dimensional visual point cloud model, and section parameters are defined; performing rapid gridding processing on the ordered point cloud data through an efficient triangulation algorithm; calculating to obtain an intersection discrete point set of the cross section and the triangular mesh; and establishing a section plane local coordinate system and implementing two-dimensional projection transformation to obtain the contour projection of the target workpiece. According to the method, high-precision, high-efficiency and automatic extraction of the section contour of the workpiece can be realized, the actual requirements of industrial field online detection and quality control are met, and the method can be applied to scenes such as part size detection, reverse engineering, quality control and contour analysis in industrial manufacturing.

The invention discloses an aircraft grid generation method and device based on conformal mapping and topological segmentation, equipment and a storage medium, and relates to the technical field of computational fluid mechanics, and the method comprises the steps: generating a surface unstructured grid according to an aircraft model, carrying out the conformal mapping to a parameter domain plane grid through Ricci flow, and calculating the compression ratio of the surface unstructured grid; if the compression ratio is smaller than a first threshold value, extracting a contour line corresponding to a second threshold value; positioning a triangular grid corresponding to the isoline discrete points in the parameter domain plane grid by using a spatial partitioning algorithm, inserting the discrete points into the grid according to the relationship between the vertex distance and the edge distance, forming a closed internal boundary through a one-dimensional boundary recovery operation, and determining a corresponding region; the area is inversely mapped back to the three-dimensional surface, the surface grid is updated, the mapping and judging steps are repeated until the compression ratio meets the requirement, and the parameter domain plane grid is set as a target aircraft grid, so that the efficiency of generating the aircraft grid is improved.

The invention belongs to the field of robot trajectory optimization, and particularly relates to a trajectory optimization method and system considering redundancy characteristics of a robot macro-micro force control grinding and polishing system. The method comprises the following steps: S1, generating an initial grinding and polishing track based on a CAD model, and providing a discrete point curvature estimation method to estimate the curvature of each path point; s2, on the basis of the curvature of the discrete points, a mixed optimization target of a curve energy item and a curvature change rate about curvature square is constructed, arc length and redundancy direction optimization constraints are considered, and an adjusted grinding and polishing track is generated through a nonlinear optimization method; s3, on the basis of the optimized initial points, an initial dominant point selection method based on an angle-chord height-chord length joint criterion is adopted, and the trajectory precision is kept while the number of dominant points is reduced; and S4, on the basis of the optimized curved surface dominant points, an improved least square iterative approximation method is provided to generate a parameterized B-spline grinding and polishing path generation method, and final grinding and polishing track generation is achieved.

The invention provides a shot blasting path generation method for an irregular profile double-curvature integral wallboard, which comprises the following steps of: firstly, extracting geometric information of the integral wallboard, and discretizing an appearance surface of the wallboard by utilizing a discretization algorithm to obtain a uniformly distributed and fully covered discrete point set; extracting the geometric information of the wallboard at the discrete points, solving the matrix thickness of the discrete points at the ribs, and storing the matrix thickness in an XML (Extensible Markup Language) data structure, so that the extracted geometric information can be conveniently used for shot blasting path planning and process parameter calculation; and finally, based on the discrete point two-dimensional array, performing shot blasting path planning through the extreme curvature value and the direction of the discrete point, and combining an iterative adjustment strategy to ensure that the generated path meets the requirements of equipment parameters and coverage rate, so that the adaptability to complex geometry and accurate range control of spacing are realized, and the path spacing meets the range requirement. Manual adjustment is reduced, and the path generation process is stable and reliable.

The invention discloses a motion trail smooth optimization method, device and equipment and a storage medium, and the method comprises the steps: collecting the distance data from a Z axis to a glass surface, and obtaining m groups of two-dimensional discrete point data through combining the lineation distance data of an X axis and the sampling density rho, inputting a plurality of feature vectors established by combining the two-dimensional discrete point data (Xpos, Zpos), the sampling density rho, the lineation speed V, the track included angle theta and the normalized parameter k into a regression model to obtain a recommended parameter beta, and calculating a tension control parameter alpha according to the recommended parameter beta; and finally, calculating a tangential vector of each two-dimensional discrete point and establishing a plurality of sections of Hermite spline curves to obtain a trajectory optimization curve so as to finish smooth optimization of the slave axis motion curve, and optimizing each section of slave axis motion trajectory according to a trajectory included angle theta. The following oscillation of the slave shaft caused by a high-frequency instruction during high-speed operation of the main shaft can be inhibited, and the coincidence degree of the actual track and the theoretical track of the slave shaft is effectively improved.

This invention belongs to the technical field of multi-scale particle transport simulation, and discloses a kinetic solution method and system for multi-scale particle transport simulation. The method includes: a particle kinetic model based on discrete velocity form; selecting discrete points in velocity space; calculating the distribution function along the direction of the discrete points in velocity space at the center of the grid according to a compact scheme; accumulating the contribution values ​​of the distribution function to macroscopic physical quantities; scanning the downstream grid sequentially based on the direction of the discrete points in velocity space; changing the discrete points in velocity space and repeating the operation until the contribution values ​​of the grid at all discrete points in velocity space are accumulated, and obtaining the final macroscopic physical quantities; outputting the calculation results of macroscopic physical quantities when the convergence condition is met. This invention discloses a kinetic solution method and system for multi-scale particle transport simulation that balances computational memory, accuracy, and efficiency, and is particularly suitable for multi-scale particle transport simulation with a large number of discrete points in velocity space and strong heterogeneity.

The invention belongs to the technical field of underwater sound channels and physical fields, and relates to a numerical mixing calculation method for a forward scattering sound field of an underwater regular-shape object. The method comprises the following steps: acquiring an initial projection contour of a three-dimensional geometric model of a regular shape object in a specified projection plane and an incident direction; discrete points of the initial projection contour are divided into upper boundary points and lower boundary points; performing interpolation completion on the upper boundary points and the lower boundary points, so that the upper boundary points and the lower boundary points are the same in number and share abscissas; discretizing the supplemented projection contour by using a trapezoidal plane element; calculating a forward scattering far field morphological function based on a discretization result; and calculating the forward scattering target intensity of the object. The method provided by the invention can cover the shape range of the scattering object suitable for the analysis method, does not need to calculate special functions such as Bessel and spherical harmonic and a complex series expression formed by combining the special functions, and obviously improves the calculation efficiency.

The application discloses an aircraft grid generation method and device based on conformal mapping and topological segmentation, equipment and a storage medium, relates to the technical field of computational fluid dynamics, and comprises the following steps: generating a surface unstructured grid according to an aircraft model, conformally mapping the surface unstructured grid to a parameter domain plane grid through Ricci flow, and calculating the compression rate of the parameter domain plane grid; if the compression rate is less than a first threshold value, extracting an isoline corresponding to a second threshold value; locating a triangular grid corresponding to an isoline discrete point in the parameter domain plane grid by using a spatial partition algorithm, inserting the discrete point into the grid according to the vertex distance and edge distance relationship, forming a closed internal boundary through a one-dimensional boundary recovery operation, and determining a corresponding region; inversely mapping the region back to a three-dimensional surface, updating the surface grid, and repeating the mapping and judging steps until the compression rate meets the requirement, and setting the parameter domain plane grid as a target aircraft grid, so that the efficiency of generating the aircraft grid is improved.

The invention belongs to the technical field of optical lenses, and particularly relates to a discrete point source bus set envelope fitting and optimization design method for expanding a light source to realize uniform illumination. Dispersing the expansion source into point light sources, and obtaining a cluster bus set based on an accurate mathematical model of the point light sources; according to an envelope principle, selecting an outmost contour including all features as much as possible, and fitting the outmost contour by using a polynomial under a polar coordinate system to serve as an initial generatrix of the light distribution lens; a root mean square error of an irradiance distribution numerical value obtained by running a two-dimensional ray tracing program is used as an evaluation index to optimize and fit various coefficients of a polynomial, finally, the lens suitable for achieving uniform illumination through extended source light distribution is obtained, and the problem that mathematical modeling of an extended source light distribution light path is difficult is solved; through verification, the performance is good, the irradiation uniformity is gradually improved to 92.7% from 76.8% of an initial structure, the edge irradiance gradient is improved to 0.45 from the initial 0.27, and the edge dispersion problem is effectively weakened.

This invention discloses an energy-saving color matching method and apparatus for electronic maps based on semantic relationships. The method includes: spatializing the colors of the original map; modeling the map colors as discrete point clouds in the CIELab color space according to perceptual uniformity rules; analyzing the discrete point clouds using a clustering algorithm based on color difference rules to extract semantic relationships between map colors, including differences, associations, and orders; establishing a color semantic relationship evaluation model based on the extracted semantic relationships; reorganizing the colors in the RGB color space and establishing an energy consumption estimation model based on the red, green, and blue component values ​​of each color; establishing the map energy-saving problem as a multi-objective, multi-constraint optimization problem; using a heuristic search algorithm to find energy-saving colors; and obtaining the optimal energy-saving color matching scheme by analyzing the Pareto front solution set. This invention is applicable to vector and raster electronic maps of different scales and themes, automatically extracting semantic relationships to reduce map energy consumption.

The present application relates to the technical field of data analysis, and specifically relates to a user portrait label value analysis method and system, comprising obtaining a first subject user portrait label value mapping to form a first user discrete point group, generating virtual user data to be respectively input into a first subject and a second subject clustering model to form a first simulated discrete point group and a second simulated discrete point group, transforming a vector through space fitting records and identifying an unmatched discrete point group, performing cutting and scaling processing on the unmatched discrete point group and recording a cutting and scaling vector, and obtaining a second user discrete point group through coordinate transformation of the first user discrete point group through the transformation vector and the cutting and scaling vector, and obtaining a user portrait label value of the second subject. The present application can realize user portrait conversion across label systems without sharing original user data, and solve the label matching problem caused by loss of normalized information.

The application discloses a mixed layer thickness identification method based on a carrier coordinate system, and the identification method comprises the following steps: importing flow field data; uniformly processing the flow field data, and marking fluid domains and non-fluid domains; identifying active flow and passive flow, and respectively tracking characteristic flow lines of the active flow and the passive flow to determine the boundary of a mixed area; establishing a carrier coordinate system, and determining an integral path set of the mixed area; integrating all the integral paths to obtain the mixed layer thickness of each integral path, i.e., to obtain the mixed layer thickness along the flow direction of the mixed area; screening out integral path discrete points with abnormal growth or reduction of the mixed layer thickness before and after the mixed layer thickness, and removing the integral path discrete points; drawing a curve graph of the mixed layer thickness along the flow direction, and completing the identification of the mixed layer thickness. The application is applied to the field of flow field mixing, can meet the calculation requirements of the mixed layer thickness under more general flow field conditions, and provides more accurate data for flow field analysis, mixed supercharging evaluation and other applications.

The application discloses a jump line generation method based on profile laser cutting contour geometric information, extracts the start and end points of corresponding geometric information based on given geometric information, then obtains the coordinates on the circumscribed circle in the direction of the corresponding normal vector, then generates corresponding three-dimensional space jump lines, simultaneously calculates the discrete point coordinates of the jump lines and the direction vectors of the discrete points generated by circular arc discrete interpolation, and can introduce intermediate points for special scenes, so that the machining path meets the actual scene requirements, for example, the path is always located on the upper part, the intermediate points can be arranged on the upper part, and collision and interference can be avoided.

This application discloses a method, apparatus, electronic device, and storage medium for character stroke extraction. The method includes: acquiring the glyph curve of the character to be processed; performing discrete sampling on the glyph curve to obtain a discrete point set; selecting at least two stroke key points from the discrete point set to obtain a closed curve containing at least two stroke key points, thereby obtaining the character strokes. Through the above method, this application can improve the speed and accuracy of character stroke extraction.

The present disclosure provides a phase configuration method and electronic equipment, which relates to the technical field of communication, obtains a discrete point set corresponding to a phase interval of a first optical MZI array system, calibrates a first phase set obtained according to the discrete point set by using a test optical signal corresponding to the first optical MZI array system, adjusts each phase in the first optical MZI array system by using a second phase set, obtains a second optical MZI array system, and the second phase set includes at least one phase to be adjusted in the calibrated first phase set; obtains an optimal phase configuration corresponding to the second optical MZI array system according to a first output light intensity of each port in the second optical MZI array system under a continuous phase condition and a heuristic algorithm, and configures the phase of the second optical MZI array system by using the optimal phase configuration. The pre-calibration operation can be reduced, the power consumption can be reduced, the expressiveness and robustness of the system can be improved, the accuracy of the phase configuration can be improved, and resources can be saved.

The application relates to the technical field of image processing, and provides an image comparison method and device, electronic equipment and a storage medium. The method comprises the following steps: obtaining a texture feature matrix of an image in a horizontal axis direction and a texture feature matrix of the image in a vertical axis direction; performing feature fusion on the texture feature matrix in the horizontal axis direction and the texture feature matrix in the vertical axis direction to obtain a fused texture feature matrix corresponding to the image, so that the fused texture feature matrix can accurately express the features of the image in the horizontal axis and the vertical axis; removing discrete points in the fused texture feature matrix to obtain a target texture matrix corresponding to the image, so that the target texture matrix can more accurately express the texture features of the image, thereby avoiding the influence of the discrete points in the image on the image comparison result; and finally, performing difference comparison on the target texture matrix of the image and a target texture matrix of a to-be-compared image, so that the obtained difference comparison is accurate.

The invention discloses a PVT synthesis-based ultrasonic inspection control method for a welding seam of a connecting pipe of a reactor pressure vessel. The method comprises the following steps: S1, setting parameters; s2, obtaining a spatial trajectory discrete point sequence; s3, the shaft speed is constant; s4, any two continuous spatial trajectory discrete points in the spatial trajectory discrete point sequence are selected; according to the movement time t between the two points, the axis and the component velocity of the axis are obtained through calculation; s5, based on the three-axis component, PVT interpolation operation between two points is carried out, and a smooth track is generated; s6, repeatedly executing the steps S4-S5; and S7, performing incremental motion according to the stepping interval, and repeating the steps S2-S6. Aiming at the special geometrical characteristics of the saddle surface shape of the connecting pipe, a PVT synthesis algorithm is designed by virtue of a linear rotating device, a scanning motion track is obtained, the fitting state of the ultrasonic probe and the detected surface is accurately controlled, and automatic nondestructive detection of the welding seam in the inner fillet connecting pipe area is realized.

The application belongs to the technical field of power systems, and discloses an EIGD time-delay power system stability analysis method and device, wherein the method comprises the following steps: a linear model of the time-delay power system is established, and the differential equation in the model is converted into an abstract Cauchy problem by using an infinitesimal generator; a group of discrete points is selected for each time-delay interval of the time-delay power system, a discrete function space is established according to the discrete points, and the time-delay variable is discretized to generate a low-order partial infinitesimal generator discretization matrix; displacement and inverse transformation are performed on the discretization matrix to obtain an inverse matrix, and the characteristic value of the system is obtained according to the inverse matrix; the characteristic value is checked by using the Newton method to obtain the accurate characteristic value and the characteristic vector, and the stability of the time-delay power system is analyzed and judged. The method solves the problem of large matrix LU decomposition calculation amount in the existing EIGD time-delay power system characteristic value sparse calculation method based on the DDAE model.

The application discloses a large model access control system and method based on calling intention analysis, relates to the technical field of natural language data processing, and standardizes, cleans and unifies the format of user input text, and obtains an intention vector based on an intention encoder constructed by a large model. Then, the intention vector continuously generated is arranged in sequence to form a discrete point sequence with a dialogue round as a discrete time axis, and the end points are connected to form a dialogue moving track through sliding window maintenance. The curvature of each point of the track and the maximum curvature in the sliding window are calculated by processing the obtained digital data. A history normal curvature library is established for each user, a personalized curvature threshold is obtained, and an abnormal curvature segment is marked, a first-level early warning is triggered through a continuous abnormality counter, a second-level early warning is triggered and a graded response is executed through accumulated abnormal curvature integration and variance determination. The application achieves the purposes of reducing the false positive rate and improving the accuracy and reliability of large language model security detection.

The invention relates to a multi-beam sounding data quality real-time evaluation method based on a Delaunay triangulation network model, and the method comprises the steps: carrying out the thinning of an original discrete data set, generating a thinned data set, and constructing a Delaunay triangulation network; performing fragmentation processing on the original data set according to the Delaunay triangulation network, and extracting an accurate data range of the fragmented data set; and respectively carrying out gridding processing on each fragmented data set to generate fragmented grids, and carrying out multi-beam sounding data quality real-time evaluation. According to the method, the gridding processing of mass and large-area discrete point data of the multi-beam sounding data can be quickly and accurately realized, the problems of real-time quality evaluation of the mass multi-beam sounding data and the like are solved, and the method can be widely applied to the field of multi-beam sounding data processing.

The application provides a driving track control method of a vehicle in an S-shaped curve, a vehicle and a storage medium, and relates to the technical field of intelligent driving. A discrete point set is fitted into a target curve of an Nth order polynomial, and since N is greater than or equal to 5, the shape of the real lane boundary can be more accurately and completely represented, and the reliability is high. Since the shape of the real lane boundary can be more accurately and completely represented, the target curve can be accurately identified as an S-shaped curve. In the case of identifying that the target curve is an S-shaped curve, the discrete point set of the S-shaped curve and the current motion state data of the vehicle are input into a preconfigured model predictive control algorithm to output a target driving track that minimizes the sum of the steering angles when the vehicle drives in the lane of the S-shaped curve, so that the lateral swing of the vehicle is small, the driving stability is high, and the driving safety and the driving experience of the user are further improved.

The application discloses a grid division method for a closed cylindrical surface with complex boundary constraints, and comprises the following steps: obtaining a plurality of boundary rings of the closed cylindrical surface and a plurality of parameters from a solid model, wherein the parameters comprise a coordinate system of the closed cylindrical surface, a radius of the closed cylindrical surface and an axial direction of the closed cylindrical surface; sequentially performing discrete calculation on each edge in each boundary ring of the closed cylindrical surface to obtain a plurality of discrete points; all the discrete points corresponding to the boundary ring form a discrete point linked list corresponding to the boundary ring; the arrangement order of the discrete points in the discrete point linked list is consistent with the arrangement order of the discrete points in the edge; each discrete point in the discrete point linked list corresponding to each boundary ring is subjected to mapping transformation to obtain a corresponding plane mapping point, and a mapping relationship between the plane mapping point and the corresponding discrete point is established. The application can solve the technical problem that there is no grid division method for a closed cylindrical surface with complex boundary constraints.