Integrated design method based on reamer assembly three-dimensional simulation model

Abstract

The invention relates to the technical field of three-dimensional simulation model design, in particular to an integrated design method based on a reamer assembly three-dimensional simulation model. The method comprises the following steps: performing multi-view three-dimensional scanning on a target reamer to generate three-dimensional model grid data of the reamer; performing model attribute setting according to the reamer three-dimensional model grid data, and generating an optimized three-dimensional reamer attribute model; cutting load simulation is carried out according to the optimized three-dimensional reamer attribute model, structural geometric parameter influence analysis is carried out, and reamer key influence parameters are generated; structural variable disturbance processing is conducted according to the reamer key influence parameters, finite element material distribution optimization is conducted on the optimized three-dimensional reamer attribute model, three-dimensional simulation model reconstruction is conducted, and optimized reamer three-dimensional simulation model data are generated. Through dynamic characteristic analysis of the reamer assembly, design optimization of the three-dimensional model of the reamer assembly is achieved, and the problem that the machining quality and efficiency are affected due to chatter generated in the machining process of the reamer assembly is solved.

Description

Technical Field

[0001] The present invention relates to the technical field of three-dimensional simulation model design, and in particular to an integrated design method for a three-dimensional simulation model based on a reamer assembly. Background Art

[0002] Among numerous cutting processes, reaming, as a precise cutting method, is widely used in the finish machining of holes to obtain precise dimensions, shapes, and good surface quality. As a key tool for reaming, the structure and performance of the reamer directly affect the machining quality and efficiency. To meet the machining requirements of complex parts, adjustable or combined reamer assemblies have emerged. Such reamer assemblies are usually composed of multiple blades or tool heads and a tool shank combined through specific connection structures (such as threads, wedge structures, or hinge structures), and can achieve the machining of holes with different diameters or special shapes. During the cutting process, when the natural frequency of the reamer assembly approaches the cutting frequency or its multiple frequencies, resonance will occur. Resonance will significantly amplify the tiny disturbances during the cutting process, and then trigger chatter. Chatter is a self-excited vibration, manifested as a violent relative vibration between the tool and the workpiece, resulting in defects such as ripples and scratches on the machined surface, seriously affecting the machining accuracy and surface quality. In addition, the violent vibration will also accelerate the tool wear, reduce the tool life, and generate harsh noises, affecting the health and safety of the operator. However, most of the traditional three-dimensional simulation model designs of reamer assemblies focus on static or quasi-static analysis, lacking accurate simulation of the dynamic response of the reamer assembly during the actual cutting process, unable to consider the dynamic characteristic changes of the reamer assembly under complex cutting conditions, unable to effectively suppress chatter, and difficult to achieve the optimal design. Summary of the Invention

[0003] Based on this, the present invention provides an integrated design method for a three-dimensional simulation model based on a reamer assembly to solve at least one of the above technical problems.

[0004] To achieve the above object, an integrated design method for a three-dimensional simulation model based on a reamer assembly includes the following steps:

[0005] Step S1: Perform multi-view three-dimensional scanning on the target reamer to obtain the original point cloud data of the reamer assembly; reconstruct a three-dimensional reamer model based on the original point cloud data of the reamer assembly to generate an initial three-dimensional reamer geometric model; perform mesh division on the geometric model according to the initial three-dimensional reamer geometric model to generate the mesh data of the three-dimensional reamer model;

[0006] Step S2: Identify reamer components based on the mesh data of the reamer 3D model to obtain reamer component identification data; set the model attributes of the initial 3D reamer geometric model through the reamer component identification data to generate an optimized 3D reamer attribute model; perform cutting load simulation on the optimized 3D reamer attribute model to generate reamer simulation modal analysis data;

[0007] Step S3: Analyze the influence of structural geometric parameters on the optimized 3D reamer attribute model through the reamer simulation modal analysis data to generate key influencing parameters of the reamer; perform structural variable perturbation processing based on the key influencing parameters of the reamer and conduct 3D grid node spatial mapping to obtain reamer structure sensitivity mapping data;

[0008] Step S4: Optimize the finite element material distribution of the optimized 3D reamer attribute model using the reamer structure sensitivity mapping data and reconstruct the 3D simulation model to generate optimized 3D reamer simulation model data.

[0009] The present invention realizes the accurate digital acquisition of the reamer geometry through multi-view 3D scanning of the reamer assembly, constructs a high-fidelity initial 3D geometric model, avoids the errors introduced by simplification or approximation in traditional modeling methods, and can more realistically reflect the actual structure of the reamer assembly. Next, by meshing the initial 3D model, the applicability of the model in numerical simulation is ensured. Through reamer component identification, the complex reamer structure is decomposed into independent components, and corresponding physical properties are assigned to each component. This decomposition and attribute setting enable the simulation analysis to more precisely simulate the real behavior of the reamer during the cutting process, avoiding the errors caused by treating the reamer as a homogeneous body as a whole, and especially can accurately consider the influence of different materials and connection structures on the overall dynamic performance. On this basis, cutting load simulation is performed to obtain the dynamic response of the reamer under cutting conditions, providing quantitative data for understanding the modal characteristics of the reamer assembly and potential chatter risks. This simulation-based analysis method replaces the traditional empirical design, can more deeply reveal the dynamic characteristics of the reamer assembly under complex working conditions, and the sensitivity mapping realizes a leap from qualitative analysis to quantitative analysis. Using these sensitivity information, the reamer structure can be adjusted more precisely to suppress chatter.

[0010] Through the reasonable allocation of material distribution, the overall performance of the reamer is improved. The goal of this optimization is to be able to increase the natural frequency of the reamer as much as possible and stay away from the cutting frequency while meeting the cutting performance, thereby reducing the possibility of resonance and chatter. Therefore, an integrated design method based on a three-dimensional simulation model of a reamer assembly of the present invention reconstructs an initial three-dimensional reamer geometric model by performing a multi-view three-dimensional scan of the target reamer, meshes the initial model, and performs cutting load simulation. The modal analysis data is used to perform structural geometric parameter influence analysis, determine key influencing parameters, and perform structural variable perturbation processing to obtain the reamer structure sensitivity mapping data, achieving a leap from qualitative analysis to quantitative analysis. Finite element material distribution optimization is performed based on sensitivity mapping data, and an optimized three-dimensional simulation model of the reamer is reconstructed, which effectively improves the dynamic performance and vibration resistance of the reamer and realizes the precise design and manufacture of the reamer.

[0011] Preferably, step S1 comprises the following steps:

[0012] Step S11: performing surface cleaning treatment on the target reamer to obtain a cleaned target reamer component sample;

[0013] Step S12: performing a multi-view three-dimensional scan on the cleaning target reamer component sample to obtain original point cloud data of the reamer component;

[0014] Step S13: performing point cloud data filtering processing according to the original point cloud data of the reamer assembly, and reconstructing a three-dimensional reamer model to generate an initial three-dimensional reamer geometric model;

[0015] Step S14: performing a blade area geometric feature line analysis based on the initial three-dimensional reamer geometric model to generate component blade area feature data;

[0016] Step S15: Setting blade size control parameters for the initial three-dimensional reamer geometric model based on the component blade area feature data to generate blade size control parameters;

[0017] Step S16: meshing the initial three-dimensional reamer geometric model using blade size control parameters to generate reamer three-dimensional model mesh data.

[0018] By performing surface cleaning on the target reamer, the present invention can remove impurities such as dirt and grease on the surface of the reamer, ensuring the accuracy and reliability of the point cloud data obtained by subsequent three-dimensional scanning. Multi-view three-dimensional scanning can comprehensively obtain the geometric information of the reamer assembly, avoiding model missing or deformation caused by single-view scanning. Based on the initial model, geometric feature line analysis of the cutting edge area is carried out to extract key geometric features of the cutting edge area, such as the cutting edge profile, chamfer, etc., and generate the feature data of the cutting edge area of the assembly. According to the set cutting edge size control parameters, mesh generation is performed on the initial three-dimensional reamer geometric model to generate the mesh data of the reamer three-dimensional model that meets the requirements of simulation analysis. Guided by the cutting edge size control parameters, mesh generation can ensure the generation of high-quality meshes in the cutting edge area, improving the accuracy and reliability of the simulation results. Especially in the analysis of key parameters such as stress concentration and deformation in the cutting edge area, more accurate simulation results can be provided.

[0019] Preferably, step S14 includes the following steps:

[0020] Step S141: Identify the cutting edge area of the initial three-dimensional reamer geometric model to generate the reamer cutting edge area data;

[0021] Step S142: Detect the edge contour of the reamer cutting edge area data and extract geometric feature lines to obtain the cutting edge feature line data;

[0022] Step S143: Perform cutting edge surface curvature analysis based on the cutting edge feature line data to generate the cutting edge line curvature data;

[0023] Step S144: Mark the key points of the cutting edge features on the reamer cutting edge area data based on the cutting edge line curvature data to obtain the candidate cutting edge feature point data;

[0024] Step S145: Connect the cutting edge line segments according to the candidate cutting edge feature point data and perform geometric feature parameterization to generate the feature data of the cutting edge area of the assembly.

[0025] By identifying the cutting edge area of the initial three-dimensional reamer geometric model, the present invention can simplify the complex reamer model and focus on the cutting edge area that plays a key role in cutting performance, improving the efficiency and pertinence of subsequent analysis. By performing surface curvature analysis on the cutting edge feature line, the degree of bending of the cutting edge surface can be quantified, and cutting edge line curvature data can be generated. The curvature data can reflect the geometric shape changes of the cutting edge, such as the sharpness of the cutting edge and the size of the chamfer. This information is crucial for evaluating the cutting performance of the tool. Based on the cutting edge line curvature data, key points of the cutting edge features of the reamer cutting edge area can be marked, such as the tool tip and the inflection point of the cutting edge. These key points are the key features of the cutting edge geometry and can more precisely describe the geometry of the cutting edge. Connecting the candidate cutting edge feature points with cutting edge line segments and performing geometric feature parameterization can convert the discrete feature point data into continuous line segments and express the geometric features of the cutting edge in a parametric way, such as the length and angle of the line segments.

[0026] Preferably, step S2 includes the following steps:

[0027] Step S21: Identify the reamer components based on the mesh data of the reamer three-dimensional model to obtain reamer component identification data;

[0028] Step S22: Set the model attributes of the initial three-dimensional reamer geometric model through the reamer component identification data to generate an optimized three-dimensional reamer attribute model;

[0029] Step S23: Perform finite element processing on the optimized three-dimensional reamer attribute model to obtain a finite element three-dimensional reamer simulation model;

[0030] Step S24: Perform cutting load simulation on the finite element three-dimensional reamer simulation model and perform multi-order natural frequency / vibration mode analysis to generate reamer simulation modal analysis data.

[0031] By identifying the components of the reamer three-dimensional model mesh data, the present invention can decompose the reamer assembly into different functional components, such as the tool body, tool shank, threads, etc., and assign corresponding attribute information to each component, such as material properties, boundary conditions, etc., providing a more refined model for subsequent simulation analysis. Using the identified component information, the model attributes of the initial three-dimensional reamer geometric model can be set, such as assigning different material properties to different components (for example, the tool body is made of cemented carbide and the tool shank is made of steel), and defining the connection relationship between the components. Based on the finite element model, cutting load simulation can be performed to simulate the stress state, deformation conditions, and temperature distribution of the tool during actual cutting. Performing multi-order natural frequency / vibration mode analysis on the finite element model can obtain the natural frequencies and vibration modes of the reamer assembly under different vibration modes. This information is the key to analyzing the vibration characteristics of the tool and predicting chatter phenomena.

[0032] Preferably, step S22 includes the following steps:

[0033] Step S221: Query the component material properties of the reamer component recognition data through a preset reamer component material property library to generate component material property data;

[0034] Step S222: Map the component material property data to the initial three-dimensional reamer geometric model to obtain a three-dimensional reamer property model;

[0035] Step S223: Analyze the component contact relationship of the reamer component recognition data through the three-dimensional reamer property model to generate component contact relationship data;

[0036] Step S224: Set the contact constraint parameters for the three-dimensional reamer property model based on the component contact relationship data to generate component contact constraint condition data;

[0037] Step S225: Set the model constraints for the three-dimensional reamer property model through the component contact constraint condition data to generate an optimized three-dimensional reamer property model.

[0038] In the present invention, by querying the component material properties of the reamer component recognition data through a preset reamer component material property library, the corresponding material properties, such as density, elastic modulus, Poisson's ratio, etc., can be quickly and accurately assigned to each component, avoiding the cumbersome operation of manually inputting material properties, improving the modeling efficiency, and reducing human errors. Based on the three-dimensional reamer property model, the component contact relationship of the reamer component recognition data is analyzed, and the contact relationship between components, such as the connection between the tool body and the tool holder, the fit between the blade and the tool body, etc., can be automatically recognized, and the component contact relationship data is generated. According to the generated component contact relationship data, the contact constraint parameters for the three-dimensional reamer property model can be set, such as setting the contact type (such as bonded contact, frictional contact), friction coefficient, etc., and the component contact constraint condition data is generated. By setting the model constraints for the three-dimensional reamer property model through the component contact constraint condition data, the interaction between components, such as the fixed connection between the tool body and the tool holder, the sliding friction between the blade and the tool body, etc., can be simulated.

[0039] Preferably, step S223 includes the following steps:

[0040] Mark the component meshes of the reamer component recognition data through the three-dimensional reamer property model to obtain the reamer component mesh marking data;

[0041] Perform geometric proximity search based on the reamer component mesh marking data and perform component surface mesh node processing to generate component surface mesh node data;

[0042] Calculate the spatial node distance on the component surface based on the mesh node data of the component surface to generate the spatial node distance data of the component surface;

[0043] Screen the contact surface node pairs of the component surface mesh node data through the spatial node distance data of the component surface to obtain the potential component contact surface data;

[0044] Analyze the actual reamer assembly relationship of the three-dimensional reamer attribute model using the potential component contact surface data to generate the actual reamer assembly data;

[0045] Process the relative motion form of the contact surface for the potential component contact surface data through the actual reamer assembly data, and determine the type of contact relationship to generate the component contact relationship data.

[0046] The present invention performs geometric proximity search based on the marked mesh data, can quickly find the component surface mesh nodes that are relatively close in space, and perform node processing to generate the component surface mesh node data. Using the component surface mesh node data to calculate the spatial node distance on the component surface can quantify the distance between different component surface nodes and generate the spatial node distance data of the component surface. By setting an appropriate distance threshold, screening the contact surface node pairs of the component surface mesh node data can preliminarily determine the potential component contact surface and obtain the potential component contact surface data. Using the potential component contact surface data to analyze the actual reamer assembly relationship of the three-dimensional reamer attribute model, such as analyzing the assembly relationship between the tool body and the tool handle, the mating relationship between the blade and the tool body, etc., and generating the actual reamer assembly data. Through the actual reamer assembly data, process the relative motion form of the contact surface for the potential component contact surface data, such as determining whether the contact surfaces are fixedly connected, slidingly connected or rollingly connected, etc., and determine the type of contact relationship to generate the final component contact relationship data.

[0047] Preferably, step S24 includes the following steps:

[0048] Step S241: Obtain the actual working frequency range of the target reamer;

[0049] Step S242: Process the resonance natural frequency value according to the actual working frequency range of the target reamer, and perform mode vector analysis to generate the natural frequency / mode vector data;

[0050] Step S243: Perform modal mode / resonance risk modal processing according to the natural frequency / mode vector data, and set multi-order modal parameters to generate a mode / resonance multi-order modal parameter list;

[0051] Step S244: Perform energy contribution analysis according to the mode / resonance multi-order modal parameter list, and superimpose the target modal simulation parameters to obtain the modal superposition simulation parameters;

[0052] Step S245: Perform cutting load simulation using the finite element three-dimensional reamer simulation model based on the modal superposition simulation parameters to obtain dynamic cutting load response data;

[0053] Step S246: Perform multi-order modal response superposition based on the dynamic cutting load response data, and perform frequency-domain response analysis on the reamer part to obtain reamer simulation modal analysis data.

[0054] By obtaining the actual working frequency range of the target reamer, the present invention can clarify the key frequency range of the simulation analysis, avoid unnecessary calculations, and improve the simulation efficiency. Processing the resonance natural frequency values according to the actual working frequency range and performing mode shape vector analysis can extract the natural frequencies and mode shape vectors related to the working frequency range. Based on the natural frequency / mode shape vector data, modal mode shape / resonance risk modal processing can be performed to identify potential resonance risk modes. Using the modal superposition simulation parameters and the finite element three-dimensional reamer simulation model to perform cutting load simulation can obtain dynamic cutting load response data, which contains the response information of different modes to the cutting load. Performing multi-order modal response superposition on the dynamic cutting load response data and performing frequency-domain response analysis on the reamer part can obtain reamer simulation modal analysis data, which can more accurately simulate the dynamic response of the reamer under complex cutting conditions, effectively improve the accuracy and efficiency of the simulation analysis, and provide more reliable guidance for the design optimization of the reamer. Through frequency-domain response analysis, the vibration response of the reamer at different frequencies can be clearly identified, providing key information for suppressing chatter.

[0055] Preferably, step S3 includes the following steps:

[0056] Step S31: Analyze the influence of structural geometric parameters on the optimized three-dimensional reamer attribute model through the reamer simulation modal analysis data to generate key influence parameters of the reamer;

[0057] Step S32: Perform reamer structure feature mapping according to the key influence parameters of the reamer, and perform variable geometric structure quantization processing to generate key structure variable data for reamer adjustment;

[0058] Step S33: Perform structural variable perturbation processing on the key structure variable data for reamer adjustment using the preset structural variable perturbation rules, and perform finite difference calculation to generate structural variable sensitivity data;

[0059] Step S34: Perform three-dimensional grid node space mapping on the optimized three-dimensional reamer attribute model through the structural variable sensitivity data to obtain reamer structure sensitivity mapping data.

[0060] Through the analysis of the influence of structural geometric parameters on the optimized three-dimensional reamer attribute model using the reamer simulation modal analysis data, it is possible to determine which structural geometric parameters have a greater impact on the dynamic characteristics of the reamer (such as natural frequency, vibration mode), and generate the key influence parameters of the reamer. Using the preset structural variable perturbation rule to perform structural variable perturbation processing on the key structural variable data of the reamer adjustment, and performing finite difference calculation, the influence degree of each key structural variable on the dynamic characteristics of the reamer can be quantified, and the structural variable sensitivity data is generated. This data reflects the sensitivity of each variable to the objective function (such as natural frequency). According to the sensitivity mapping result, the structural parameters of the reamer are modified targeted, such as increasing the material thickness or changing the shape of the high-sensitivity area, so as to improve the dynamic performance and stability of the reamer, and effectively suppress the occurrence of chatter.

[0061] Preferably, step S4 includes the following steps:

[0062] Step S41: Identify the high-sensitivity area of the reamer by presetting the sensitivity threshold for the reamer structure sensitivity mapping data to obtain the high-sensitivity area of the reamer;

[0063] Step S42: Calculate the minimum response amplitude frequency according to the high-sensitivity area of the reamer to generate the minimum target response frequency value;

[0064] Step S43: Set the structural response target for the minimum target response frequency value based on the high-sensitivity area of the reamer to obtain the reamer structure response type data;

[0065] Step S44: Set the response optimization expression for the high-sensitivity area of the reamer through the minimum target response frequency value and the reamer structure response type data to generate the reamer structure response optimization expression;

[0066] Step S45: Use the high-sensitivity area of the reamer to perform volume constraint processing on the optimized three-dimensional reamer attribute model to generate the reamer volume constraint condition data; obtain the reamer manufacturing constraint condition data;

[0067] Step S46: Based on the reamer manufacturing constraint condition data and the reamer volume constraint condition data, use the reamer structure response optimization expression to optimize the finite element material distribution of the optimized three-dimensional reamer attribute model, and perform three-dimensional simulation model reconstruction to generate the optimized reamer three-dimensional simulation model data.

[0068] Based on the high-sensitivity region of the reamer and the minimum target response frequency value, the present invention sets the structural response target, which can clarify that the optimization objective is to reduce the vibration response at this frequency and obtain the data of the reamer structural response type. By setting the response optimization expression for the high-sensitivity region of the reamer through the minimum target response frequency value and the reamer structural response type data, an objective function can be established to guide the subsequent optimization process and generate the reamer structural response optimization expression. Before optimization, the volume constraint processing of the optimized three-dimensional reamer attribute model is carried out using the high-sensitivity region of the reamer, which can ensure that the optimized reamer meets the volume requirements and generate the reamer volume constraint condition data. By optimizing the material distribution, such as changing the material density or thickness of the high-sensitivity region, the vibration response of the reamer at the target frequency is reduced, thereby improving the dynamic performance and stability of the reamer. This optimization method based on sensitivity analysis and multiple constraint conditions can more effectively improve the vibration resistance performance of the reamer and achieve the optimal design of the reamer structure on the premise of meeting the actual manufacturing requirements.

[0069] Preferably, step S45 includes the following steps:

[0070] Step S451: Evaluate the influence degree of material change on the optimized three-dimensional reamer attribute model through the high-sensitivity region of the reamer to obtain the structural region influence degree data;

[0071] Step S452: Extract the initial volume of the structure according to the optimized three-dimensional reamer attribute model to obtain the optimized structural initial volume value;

[0072] Step S453: Analyze the volume reduction ratio of the optimized structural initial volume value through the structural region influence degree data to generate the regional volume reduction ratio data;

[0073] Step S454: Calculate the maximum material topological removal volume for the high-sensitivity region of the reamer using the regional volume reduction ratio data to generate the regional maximum removal volume data;

[0074] Step S455: Carry out structural volume constraint processing according to the regional maximum removal volume data to obtain the reamer volume constraint condition data.

[0075] The present invention extracts the initial volume of the structure according to the optimized three-dimensional reamer attribute model, can obtain the volume of the reamer before optimization, and get the initial volume value of the optimized structure as the reference value for volume constraint. Through the analysis of the volume reduction ratio of the initial volume value of the optimized structure by the structure area influence degree data, the volume ratio of the material allowed to be removed in different areas can be determined, and the area volume reduction ratio data can be generated. The material ratio allowed to be removed in the high influence degree area is small, while the material ratio allowed to be removed in the low influence degree area is large. By using the area volume reduction ratio data to calculate the maximum material topology removal volume in the high-sensitivity area of the reamer, the maximum material volume allowed to be removed in each area can be determined, and the area maximum removal volume data can be generated. According to the area maximum removal volume data, the structure volume constraint processing can be carried out, the volume constraint condition can be converted into specific numerical limits, and the reamer volume constraint condition data can be obtained. This data will be used to constrain the subsequent material distribution optimization process to ensure that the optimized reamer meets the volume requirements. Brief Description of the Drawings

[0076] Figure 1 It is a schematic flow chart of the steps of the integrated design method based on the three-dimensional simulation model of the reamer assembly of the present invention;

[0077] Figure 2 is Figure 1 a schematic detailed implementation step flow chart of step S3 in

[0078] Figure 3 is Figure 1 a schematic detailed implementation step flow chart of step S4 in

[0079] The realization, functional features and advantages of the object of the present invention will be further described with reference to the embodiments and the accompanying drawings. Detailed Embodiments

[0080] The technical method of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all of the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0081] In addition, the accompanying drawings are only schematic diagrams of the present invention and are not necessarily drawn to scale. The same reference numerals in the drawings represent the same or similar parts, and thus the repeated description thereof will be omitted. Some of the block diagrams shown in the drawings are functional entities and do not necessarily correspond to physically or logically independent entities. The functional entities can be implemented in software form, or in one or more hardware modules or integrated circuits, or in different networks and / or processor methods and / or microcontroller methods.

[0082] It should be understood that although terms such as "first", "second", etc. may be used herein to describe various units, these units should not be limited by these terms. These terms are only used to distinguish one unit from another. For example, without departing from the scope of the exemplary embodiments, the first unit may be referred to as the second unit, and similarly, the second unit may be referred to as the first unit. The term "and / or" used herein includes any and all combinations of one or more of the listed associated items.

[0083] To achieve the above object, please refer to Figures 1 to 3 , the present invention provides an integrated design method based on a three-dimensional simulation model of a reamer assembly, comprising the following steps:

[0084] Step S1: Perform multi-view three-dimensional scanning on the target reamer to obtain the original point cloud data of the reamer assembly; reconstruct a three-dimensional reamer model based on the original point cloud data of the reamer assembly to generate an initial three-dimensional reamer geometric model; perform geometric model mesh division on the initial three-dimensional reamer geometric model to generate reamer three-dimensional model mesh data;

[0085] Step S2: Identify reamer components based on the reamer three-dimensional model mesh data to obtain reamer component identification data; set model attributes for the initial three-dimensional reamer geometric model through the reamer component identification data to generate an optimized three-dimensional reamer attribute model; perform cutting load simulation on the optimized three-dimensional reamer attribute model to generate reamer simulation modal analysis data;

[0086] Step S3: Analyze the influence of structural geometric parameters on the optimized three-dimensional reamer attribute model through the reamer simulation modal analysis data to generate key influence parameters of the reamer; perform structural variable perturbation processing based on the key influence parameters of the reamer and perform three-dimensional grid node spatial mapping to obtain reamer structure sensitivity mapping data;

[0087] Step S4: Use the reamer structure sensitivity mapping data to optimize the finite element material distribution of the optimized three-dimensional reamer attribute model and perform three-dimensional simulation model reconstruction to generate optimized reamer three-dimensional simulation model data.

[0088] In an embodiment of the present invention, the integrated design method based on the three-dimensional simulation model of the reamer assembly comprises the following steps:

[0089] Step S1: Perform multi-view three-dimensional scanning on the target reamer to obtain the original point cloud data of the reamer assembly; reconstruct a three-dimensional reamer model based on the original point cloud data of the reamer assembly to generate an initial three-dimensional reamer geometric model; perform geometric model mesh division on the initial three-dimensional reamer geometric model to generate reamer three-dimensional model mesh data;

[0090] In the embodiments of the present invention, the target reamer is placed on the scanning platform of a three-dimensional scanner, and the scanner parameters are adjusted, such as laser power, scanning speed, scanning accuracy, etc. The reamer is scanned from different angles to obtain multi-viewpoint cloud data. Each scan ensures that there is at least a 30% overlapping area for subsequent point cloud registration. During the registration process, the iterative closest point (ICP) algorithm is used to unify the point cloud data from different viewpoints into the same coordinate system. The denoising operation removes the spurious points and outliers generated during the scanning process, simplifies the operation, reduces the amount of point cloud data, and improves the subsequent processing efficiency. After the preprocessing is completed, the point cloud data is converted into a three-dimensional mesh model by using the surface reconstruction function in the software. An appropriate reconstruction algorithm can be selected, such as triangulation, Poisson reconstruction, etc., according to the geometric complexity and accuracy requirements of the reamer. The reconstructed three-dimensional mesh model is the initial three-dimensional reamer geometric model. The initial three-dimensional reamer geometric model is imported into a finite element analysis software (such as ANSYS, ABAQUS) for mesh generation of the geometric model. An appropriate mesh type is selected, such as tetrahedral mesh, hexahedral mesh, and the mesh size parameters are set. The mesh density should be adjusted according to the geometric complexity of the reamer and the simulation accuracy requirements. For example, for regions with a larger curvature, a denser mesh should be used. After the mesh generation is completed, the mesh data of the three-dimensional reamer model is generated.

[0091] Step S2: Identify the reamer components based on the mesh data of the three-dimensional reamer model to obtain the reamer component identification data; set the model attributes of the initial three-dimensional reamer geometric model through the reamer component identification data to generate an optimized three-dimensional reamer attribute model; perform a cutting load simulation based on the optimized three-dimensional reamer attribute model to generate the reamer simulation modal analysis data;

[0092] In an embodiment of the present invention, the three-dimensional model grid data of the reamer generated in step S1 is imported into CAD software. Using the feature recognition function of the software, different components of the reamer, such as the tool shank, tool body, cutting edge, etc., are automatically or manually recognized. During the recognition process, judgments can be made based on information such as the geometric features and topological relationships of the components. The recognition results are marked with different colors or layers to generate reamer component recognition data. According to the reamer component recognition data, model attribute settings are performed on the initial three-dimensional reamer geometric model. For example, the material property of the tool shank is set to 40Cr steel, the elastic modulus is 200 GPa, and the Poisson's ratio is 0.3; the material property of the tool body is set to cemented carbide, the elastic modulus is 600 GPa, and the Poisson's ratio is 0.22. At the same time, according to the actual processing conditions, cutting parameters are set, such as the cutting speed is 100 m / min, the feed rate is 0.1 mm / r, and the cutting depth is 2 mm. The optimized three-dimensional reamer attribute model with set attributes is imported into the finite element analysis software. According to the cutting parameters and material properties, a cutting load simulation model is established. Cutting forces, cutting heat and other loads are applied in the model, and boundary conditions are set. Cutting load simulation is carried out to obtain reamer simulation modal analysis data, including stress distribution, deformation amount, and natural frequency.

[0093] Step S3: Analyze the influence of the structural geometric parameters on the optimized three-dimensional reamer attribute model through the reamer simulation modal analysis data to generate the key influence parameters of the reamer; perform structural variable perturbation processing according to the key influence parameters of the reamer, and perform three-dimensional grid node space mapping to obtain the reamer structural sensitivity mapping data;

[0094] In an embodiment of the present invention, the influence of different structural geometric parameters (such as tool shank diameter, tool body length, cutting edge angle) on the simulation results, such as the changes in stress, deformation, and natural frequency, is analyzed for the reamer simulation modal analysis data generated in step S2. Through comparative analysis, the key influence parameters that have a greater impact on the performance of the reamer are determined, such as the tool body length and the cutting edge angle. Three values of the tool body length are selected as 50 mm, 60 mm, and 70 mm, and three values of the cutting edge angle are selected as 10°, 15°, and 20° as the structural variable perturbation values. Structural variable perturbation processing is performed on the optimized three-dimensional reamer attribute model to generate 9 different reamer models. The grid node coordinates of each group of models are compared with the grid node coordinates of the initial model, and the node displacement change amount is calculated. The node displacement change amount is associated with the corresponding structural variable perturbation value to establish a reamer structural sensitivity mapping relationship and obtain the reamer structural sensitivity mapping data. For example, when the tool body length changes by 1 mm, the displacement change amount of each node can be calculated to obtain the influence of the tool body length on the reamer structural sensitivity.

[0095] Step S4: Use the reamer structure sensitivity mapping data to optimize the finite element material distribution of the optimized 3D reamer property model, and perform 3D simulation model reconstruction to generate optimized reamer 3D simulation model data.

[0096] In the embodiment of the present invention, with minimizing the reamer deformation or maximizing the natural frequency of the reamer as the objective function and the material volume as the constraint condition, the finite element material distribution is optimized using the reamer structure sensitivity mapping data. During the optimization process, the software will automatically adjust the material density of each part of the reamer according to the sensitivity information, and finally obtain the optimized material distribution result. Import the optimized material distribution result into CAD software, and reconstruct the 3D model of the reamer according to the optimization result. For example, the areas with lower material density can be removed, and the areas with higher material density can be retained or their shapes can be adjusted. The reconstructed 3D model is the optimized reamer 3D simulation model. Import the optimized reamer 3D simulation model data into the finite element analysis software for performance verification. Through simulation analysis, verify whether the optimized reamer model meets the performance requirements, such as strength, stiffness, vibration characteristics, etc. If the requirements are met, the optimization design is completed; if the requirements are not met, it is necessary to return to the previous steps for adjustment until the requirements are met.

[0097] Preferably, step S1 includes the following steps:

[0098] Step S11: Perform surface cleaning on the target reamer to obtain a cleaned target reamer component sample;

[0099] Step S12: Perform multi-view 3D scanning on the cleaned target reamer component sample to obtain the original point cloud data of the reamer component;

[0100] Step S13: Perform point cloud data filtering on the original point cloud data of the reamer component and perform 3D reamer model reconstruction to generate an initial 3D reamer geometric model;

[0101] Step S14: Analyze the geometric feature lines of the cutting edge area based on the initial 3D reamer geometric model to generate component cutting edge area feature data;

[0102] Step S15: Set the cutting edge size control parameters for the initial 3D reamer geometric model based on the component cutting edge area feature data to generate cutting edge size control parameters;

[0103] Step S16: Perform geometric model mesh division on the initial 3D reamer geometric model through the cutting edge size control parameters to generate reamer 3D model mesh data.

[0104] In an embodiment of the present invention, a target reamer assembly sample is placed in an ultrasonic cleaner, and a suitable cleaning agent is selected, such as industrial alcohol or acetone. The ultrasonic cleaning time is set to 5 minutes, and the power is set to 100W. After cleaning, a lint-free cloth is used to wipe the surface of the reamer to remove the residual cleaning agent and dirt. Ensure that the surface of the reamer is clean, dry, and free of impurities such as oil stains and dust to improve the accuracy and quality of 3D scanning. If there is rust or an oxide layer on the surface of the reamer, it is necessary to use sandpaper or a polishing machine for polishing until the surface is smooth. The cleaned target reamer assembly sample is fixed on the rotating table of the 3D scanner. Set the scanning parameters, such as the laser line width is 0.1mm, the scanning speed is 200mm / s, and the point cloud density is 1000 points / mm 2 ³. Start the scanner and perform a 360° omnidirectional scan of the reamer to obtain multiple sets of point cloud data at different angles. There should be a certain overlapping area between each set of scan data, for example, the overlapping rate is 30%. Use the filtering function in the software to remove the noise points and outliers in the point cloud data. Gaussian filtering, median filtering and other algorithms can be selected for filtering. The filtering parameters are adjusted according to the quality and noise situation of the point cloud data. For example, set the Gaussian filtering radius to 0.5mm. After filtering, perform point cloud registration to unify multiple sets of point cloud data into the same coordinate system. Then, use the surface reconstruction function in the software to convert the point cloud data into a 3D mesh model. Triangular meshes, quadrilateral meshes or other types of meshes can be selected for reconstruction. Import the initial 3D reamer geometric model after reconstruction into CAD software (such as SolidWorks). Use the curve extraction function in the software to extract the geometric feature lines in the cutting edge area. For example, extract the contour line and chamfer line of the cutting edge. The extracted feature lines should accurately reflect the geometric shape and dimensional information of the cutting edge. The angle, width, and length parameters of the cutting edge can be calculated based on the extracted feature lines. According to the extracted feature data of the cutting edge area of the component, such as the cutting edge width and angle, combined with the design requirements and machining accuracy of the reamer, set the cutting edge size control parameters. For example, set the cutting edge width tolerance to ±0.01mm and the angle tolerance to ±0.1°. Import the initial 3D reamer geometric model into finite element analysis software (such as ANSYS). According to the set cutting edge size control parameters, perform mesh division on the geometric model. In the cutting edge area, a finer mesh division is adopted, such as a mesh size of 0.1mm. In other areas, a coarser mesh division can be adopted, such as a mesh size of 0.5mm. The mesh type can be tetrahedral mesh or hexahedral mesh. After mesh division, generate the mesh data of the 3D reamer model.

[0105] Preferably, step S14 includes the following steps:

[0106] Step S141: Identify the cutting edge area of the initial 3D reamer geometric model to generate reamer cutting edge area data;

[0107] Step S142: Perform edge contour detection on the data of the reamer blade region, and extract geometric feature lines to obtain blade feature line data;

[0108] Step S143: Analyze the curvature of the blade surface based on the blade feature line data to generate blade line curvature data;

[0109] Step S144: Based on the blade line curvature data, label the key points of the blade features for the reamer blade region data to obtain candidate blade feature point data;

[0110] Step S145: Connect the blade line segments according to the candidate blade feature point data, and perform geometric feature parameterization to generate component blade region feature data.

[0111] In the embodiment of the present invention, the reamer model is segmented into different regions by using the segmentation function of the software. According to the geometric features of the reamer, such as curvature, normal vector, etc., the blade region is identified. A curvature-based segmentation algorithm can be adopted, with the curvature threshold set to 0.5, and the region with a curvature greater than the threshold is identified as the blade region. The recognition result is marked with different colors or labels to generate the reamer blade region data. The three-dimensional model is converted into a two-dimensional image, for example, projected onto the XOY plane. The Canny edge detection algorithm is used to perform edge contour detection on the blade region. The Gaussian filter kernel size is set to 5x5, the high threshold is set to 100, and the low threshold is set to 50. The detected edge contour points are connected into line segments, and curve fitting is performed, such as using cubic spline curve fitting, to obtain the blade feature line data. Using the curvature analysis function of the software, the curvature of the blade surface is calculated. Gaussian curvature, mean curvature, principal curvature, etc. can be calculated. The sampling point spacing is set to 0.1 mm, and the curvature value of each sampling point is calculated. The calculation result is displayed in the form of a color cloud map to intuitively show the distribution of the curvature of the blade surface, generating the blade line curvature data. Based on the blade line curvature data generated in step S143, key points of the blade features are labeled for the reamer blade region data. For example, the points with curvature values reaching local maximum or minimum values can be marked as candidate blade feature points. Feature points can also be identified according to the curvature change rate, such as the curvature gradient. For example, the curvature gradient threshold is set to 0.2, and the points with a gradient greater than the threshold are marked as candidate blade feature points to obtain candidate blade feature point data. According to the spatial position relationship of the feature points, adjacent feature points are connected into line segments to form the contour line of the blade. Then, geometric feature parameterization is performed on the connected line segments. For example, a B-spline curve can be used to parameterize the blade contour line. Geometric parameters such as the length, width, angle, and curvature radius of the blade can be extracted to generate component blade region feature data.

[0112] Preferably, step S2 includes the following steps:

[0113] Step S21: Identify the reamer components based on the mesh data of the reamer 3D model to obtain the reamer component identification data;

[0114] Step S22: Set the model attributes of the initial 3D reamer geometric model through the reamer component identification data to generate an optimized 3D reamer attribute model;

[0115] Step S23: Perform finite element processing on the optimized 3D reamer attribute model to obtain a finite element 3D reamer simulation model;

[0116] Step S24: Conduct cutting load simulation on the finite element 3D reamer simulation model and perform multi-order natural frequency / vibration mode modal analysis to generate reamer simulation modal analysis data.

[0117] In the embodiments of the present invention, for example, the tool holder can be identified based on the cylindrical surface feature, and the cutting edge can be identified based on the helical surface feature. The recognition results are marked with different colors or layers, and a unique ID is assigned to each component to generate the recognition data of the reamer components. For example, the ID of the tool holder is 1, the ID of the tool body is 2, and the ID of the cutting edge is 3. Analyze geometric features such as the curvature, normal direction, and regional connectivity of the mesh model. For example, the cutting edge region usually has large curvature changes and specific normal directions, the chip flutes have concave geometric features, while the tool body and the shank have relatively smooth surfaces. By setting different feature thresholds, the reamer mesh model is segmented into different parts. For example, the cutting edge is identified using a curvature threshold above 0.7, the chip flutes are identified using a curvature threshold below -0.2, and the tool body and the shank are identified through regional connectivity. The recognition process can also be combined with user interaction, allowing the user to manually adjust the recognition results, such as adding or deleting recognition regions, so as to ensure the accuracy of recognition. The recognition results are stored in the form of reamer component recognition data, which includes the mesh index, component type, and corresponding attribute information of each component. Set the model attributes for the initial three-dimensional reamer geometric model. Different components need to be set with different material attributes. For example, the tool holder is set to 45 steel, with an elastic modulus of 210 GPa and a Poisson's ratio of 0.3; the tool body is set to cemented carbide YG8, with an elastic modulus of 550 GPa and a Poisson's ratio of 0.22. In addition, contact attributes, boundary conditions, etc. also need to be set. For example, the contact between the tool holder and the tool body is set to bonded contact, and the bottom of the reamer is set to fixed constraint. After completing the attribute setting, an optimized three-dimensional reamer attribute model is generated. According to the geometric complexity and accuracy requirements of the model, a suitable element type is selected for mesh division. For example, hexahedral elements or tetrahedral elements can be selected. Set the mesh size parameters, such as the global mesh size is 2 mm and the mesh size in the cutting edge region is 0.5 mm. Conduct a mesh quality check to ensure that the mesh quality meets the requirements. For example, the aspect ratio of the element should be less than 5, and the distortion should be less than 0.5. After completing the mesh division, a finite element three-dimensional reamer simulation model is obtained. According to the actual cutting conditions, cutting loads are applied to the finite element three-dimensional reamer simulation model. For example, tangential force and radial force can be applied, the cutting speed is 100 m / min, the feed rate is 0.1 mm / rev, and the cutting depth is 2 mm. Set the simulation parameters, such as the solution type is transient analysis and the time step is 0.01 s. Conduct a cutting load simulation to obtain data such as stress, strain, and displacement of the reamer during the cutting process. At the same time, conduct multi-order natural frequency / vibration mode analysis. For example, calculate the first 6 natural frequencies and vibration modes. Modal analysis can help understand the vibration characteristics of the reamer and avoid the occurrence of resonance phenomena. Integrate the cutting simulation and modal analysis results to generate the reamer simulation modal analysis data.

[0118] Preferably, step S22 includes the following steps:

[0119] Step S221: Query the component material attributes of the reamer component identification data through a preset reamer component material attribute library to generate component material attribute data;

[0120] Step S222: Map the component material attribute data to the initial three-dimensional reamer geometric model for component attribute mapping to obtain a three-dimensional reamer attribute model;

[0121] Step S223: Analyze the component contact relationship of the reamer component identification data through the three-dimensional reamer attribute model to generate component contact relationship data;

[0122] Step S224: Set the contact constraint parameters for the three-dimensional reamer attribute model based on the component contact relationship data to generate component contact constraint condition data;

[0123] Step S225: Set the model constraints for the three-dimensional reamer attribute model through the component contact constraint condition data to generate an optimized three-dimensional reamer attribute model.

[0124] In the embodiment of the present invention, a reamer component material attribute library is established in advance, which contains the attribute information of various common reamer materials, such as elastic modulus, Poisson's ratio, density, etc. The attribute library can be stored in the form of a database or a table. For example, use an Excel table to store information such as material name, elastic modulus, Poisson's ratio, density, etc., and assign a unique ID to each material. Match the component ID in the reamer component identification data obtained in step S21 with the material ID in the material attribute library to query the corresponding material attributes. For example, if the ID in the tool holder identification data is 1, and the ID of 45 steel in the corresponding material library is also 1, then the elastic modulus of 45 steel is queried as 210 GPa, Poisson's ratio is 0.3, and density is 7850 kg / m 3 . Map the component material attribute data generated in step S221 to the initial three-dimensional reamer geometric model. According to the ID of each component in the component identification data, assign the corresponding material attributes to the corresponding components in the three-dimensional model. For example, set the material attributes of the tool holder to 45 steel, with an elastic modulus of 210 GPa, Poisson's ratio of 0.3, and density of 7850 kg / m 3 ; set the material attributes of the tool body to cemented carbide YG8, with an elastic modulus of 550 GPa, Poisson's ratio of 0.22, and density of 14300 kg / m 3After completing the attribute mapping, a 3D reamer attribute model is obtained. Based on the 3D reamer attribute model, component contact relationship analysis is performed on the reamer component recognition data. Analyze the contact relationships between different components. For example, the contact between the tool shank and the tool body is a bonded contact, and the contact between the cutting edge and the workpiece is a frictional contact. The judgment of the contact relationship can be carried out according to the geometric position relationship and assembly relationship between the components. For example, if there is a surface-to-surface contact between two components, it is considered that there is a contact relationship between them. Represent the analysis results in the form of a table or a graph. For example, use a table to record information such as the ID of the contacting components, the contact type, and the contact surface, and generate component contact relationship data. According to the component contact relationship data generated in step S223, set the contact constraint parameters for the 3D reamer attribute model. Different contact types require different contact parameters to be set. For example, a bonded contact requires setting a bonding constraint, and a frictional contact requires setting a friction coefficient. The values of the contact parameters can be determined based on experience or by referring to relevant materials. For example, the friction coefficient between the tool shank and the tool body is set to 0.2. Store the set contact parameters as component contact constraint condition data. Apply the component contact constraint condition data generated in step S224 to the 3D reamer attribute model. According to the component contact relationship data, add corresponding contact constraints in the model. Add a frictional contact between the cutting edge and the workpiece, and set the friction coefficient to 0.2. In addition, other boundary conditions need to be set according to the actual situation, such as fixed constraints, displacement constraints, etc. For example, set the end face of the tool shank of the reamer as a fixed constraint. After completing all the constraint settings, generate an optimized 3D reamer attribute model.

[0125] Preferably, step S223 includes the following steps:

[0126] Perform component mesh marking on the reamer component recognition data through the 3D reamer attribute model to obtain reamer component mesh marking data;

[0127] Perform geometric proximity search according to the reamer component mesh marking data, and perform component surface mesh node processing to generate component surface mesh node data;

[0128] Perform component surface space node distance calculation according to the component surface mesh node data to generate component surface space node distance data;

[0129] Perform contact surface node pair screening on the component surface mesh node data through the component surface space node distance data to obtain potential component contact surface data;

[0130] Perform actual reamer assembly relationship analysis on the 3D reamer attribute model using the potential component contact surface data to generate actual reamer assembly data;

[0131] Process the relative motion form of the contact surface for the potential component contact surface data through the actual reamer assembly data, and determine the type of contact relationship to generate the component contact relationship data.

[0132] In the embodiments of the present invention, according to the reamer component identification data obtained in step S21, each component in the 3D model is marked with a mesh. For example, the mesh of the tool shank can be marked as 1, the mesh of the tool body can be marked as 2, and so on. The marking process can be realized by using the mesh grouping function of the software. The meshes belonging to the same component are divided into the same group, and the corresponding group ID is assigned to generate the reamer component mesh marking data. According to the reamer component mesh marking data generated in step S2231, geometric proximity search is performed on the meshes between different components. Spatial search algorithms such as K-d tree or Octree can be used to quickly find the relatively close mesh nodes. Set the search radius to 1 mm, and search for the neighboring nodes of each node within the specified radius range. For the searched neighboring nodes, if they belong to different components, their surface mesh nodes are extracted to generate the component surface mesh node data, including information such as node coordinates and the ID of the belonging component. According to the component surface mesh node data generated in step S2232, the spatial distance between the surface mesh nodes of different components is calculated. The Euclidean distance formula can be used to calculate the distance between two points. Calculate the minimum distance from each node to the surface nodes of other components to generate the component surface spatial node distance data, including information such as node ID, minimum distance value, and ID of the nearest neighbor node. According to the component surface mesh node data generated in step S2232, the spatial distance between the surface mesh nodes of different components is calculated. The Euclidean distance formula can be used to calculate the distance between two points. Calculate the minimum distance from each node to the surface nodes of other components to generate the component surface spatial node distance data, including information such as node ID, minimum distance value, and ID of the nearest neighbor node. Set the distance threshold to 0.1 mm. In the component surface spatial node distance data generated in step S2233, the node pairs with a distance less than the threshold are screened out as potential contact surface node pairs. These nodes form the actual contact surface. The screening result is stored as the potential component contact surface data, including information such as the ID of the contact node pair and the ID of the belonging component. Set the distance threshold to 0.1 mm. In the component surface spatial node distance data generated in step S2233, the node pairs with a distance less than the threshold are screened out as potential contact surface node pairs. These nodes form the actual contact surface. The screening result is stored as the potential component contact surface data, including information such as the ID of the contact node pair and the ID of the belonging component. Using the potential component contact surface data obtained in step S2234, combined with the 3D reamer attribute model, the actual reamer assembly relationship analysis is carried out. Analyze the assembly relationship between different components. For example, the interference fit exists between the tool shank and the tool body, and the welding connection exists between the tool body and the blade. The determination of the assembly relationship can be based on information such as design drawings and assembly processes. The analysis result is represented in the form of a tree structure or a table. For example, the tree structure is used to represent the assembly hierarchy of the reamer to generate the actual reamer assembly data.Based on the actual reamer assembly data generated in step S2235, process the relative motion forms of the contact surfaces of the potential component contact surface data. For example, if two components are fixedly connected, the relative motion form of the contact surface between them is no relative motion; if two components are hinged, the relative motion form of the contact surface between them is rotational motion. Determine the type of contact relationship according to the relative motion form of the contact surface. For example, the contact relationship type of no relative motion is bonded contact, and the contact relationship type of rotational motion is hinged contact.

[0133] Preferably, step S24 includes the following steps:

[0134] Step S241: Obtain the actual operating frequency range of the target reamer;

[0135] Step S242: Process the resonance natural frequency values according to the actual operating frequency range of the target reamer, and perform mode shape vector analysis to generate natural frequency / mode shape vector data;

[0136] Step S243: Perform modal mode / resonance risk modal processing according to the natural frequency / mode shape vector data, and set multi-order modal parameters to generate a mode shape / resonance multi-order modal parameter list;

[0137] Step S244: Perform energy contribution analysis according to the mode shape / resonance multi-order modal parameter list, and superimpose the target modal simulation parameters to obtain the modal superposition simulation parameters;

[0138] Step S245: Use the finite element three-dimensional reamer simulation model based on the modal superposition simulation parameters to perform cutting load simulation to obtain dynamic cutting load response data;

[0139] Step S246: Perform multi-order modal response superposition according to the dynamic cutting load response data, and perform frequency-domain response analysis of the reamer part to obtain the reamer simulation modal analysis data.

[0140] In the embodiments of the present invention, according to information such as the design specifications of the reamer, machining parameters (such as rotational speed, feed rate), and actual working conditions, the actual working frequency range of the target reamer is determined. For example, if the rotational speed of the reamer is 1000 rpm and the number of cutting edges is 4, the passing frequency of the teeth can be calculated as 66.67 Hz. Considering factors such as vibration during the machining process, the actual working frequency range is set to 0 Hz to 200 Hz. Perform modal analysis on the finite element three-dimensional reamer simulation model obtained in step S23. Calculate the natural frequencies and vibration modes of the reamer. Compare the calculated natural frequencies with the actual working frequency range of the target reamer obtained in step S241 to identify the natural frequency values at which resonance occurs. For example, if the calculated first-order natural frequency is 150 Hz and falls within the actual working frequency range, it is considered that there is a resonance risk at this frequency. For each natural frequency, extract the corresponding vibration mode vector, which represents the vibration pattern of the reamer at this frequency. Store the natural frequency and vibration mode vector data as natural frequency / vibration mode vector data. Based on the natural frequency / vibration mode vector data generated in step S242, analyze the modal vibration modes to identify the modes with higher resonance risks. For example, the modes with larger amplitudes and located at key parts such as the cutting edges have higher resonance risks. Select multiple modes that need to be analyzed in detail, such as the first six modes. For each mode, set corresponding modal parameters, such as damping ratio, excitation frequency, etc. Based on the vibration mode / resonance multiple-mode parameter list generated in step S243, perform an energy contribution analysis on each mode. Indicators such as modal participation factors can be used to evaluate the contribution degree of each mode to the overall vibration response. Select the modes with higher contribution degrees as the target modes, such as the modes with a contribution degree greater than 10%. Superimpose the simulation parameters of the target modes, such as excitation frequency, amplitude, etc., to obtain the modal superposition simulation parameters. The superposition method can be selected according to the actual situation, such as linear superposition or sum-of-squares superposition. Based on the modal superposition simulation parameters obtained in step S244, use the finite element three-dimensional reamer simulation model to perform a cutting load simulation. Apply dynamic cutting loads in the simulation model, such as considering the fluctuations and impacts of the cutting forces. Perform transient dynamics analysis to calculate the dynamic response of the reamer under the action of the cutting load, such as displacement, velocity, and acceleration. Based on the dynamic cutting load response data obtained in step S245, perform a multi-mode response superposition. Superimpose the responses of each mode according to the superposition method set in step S244 to obtain the total response of the reamer. Perform a frequency-domain analysis on the total response, such as using the fast Fourier transform (FFT) to convert the time-domain signal into a frequency-domain signal. Analyze the vibration responses of different parts of the reamer at different frequencies, such as the amplitudes and phases of the handle, tool body, cutting edges, etc.

[0141] As an example of the present invention, refer to Figure 2 shown in Figure 1Schematic diagram of the detailed implementation steps of step S3. In this example, step S3 includes:

[0142] Step S31: Analyze the influence of structural geometric parameters on the optimized three-dimensional reamer attribute model through the reamer simulation modal analysis data, and generate the key influence parameters of the reamer.

[0143] In the embodiment of the present invention, the reamer simulation modal analysis data generated in step S2, such as natural frequency, vibration mode, etc., is imported into data analysis software (such as MATLAB). Key indicators related to the performance of the reamer are extracted, such as the first-order natural frequency. Analyze the influence of structural geometric parameters on the optimized three-dimensional reamer attribute model, such as the influence of the shank diameter, tool body length, cutting edge angle, etc. on the key performance indicators. Methods such as the response surface method and orthogonal experimental design can be used to analyze the contribution degree of different geometric parameters to the performance indicators. For example, by changing the shank diameter, observe the change of the first-order natural frequency and calculate its sensitivity. Through analysis, determine the key influence parameters that have the greatest influence on the reamer performance, such as the tool body length and the rake angle of the cutting edge, and generate the key influence parameter data of the reamer.

[0144] Step S32: Map the reamer structural characteristics according to the key influence parameters of the reamer, and perform variable geometric structure quantization processing to generate the key structural variable data for reamer adjustment.

[0145] In the embodiment of the present invention, according to the key influence parameters of the reamer generated in step S31, such as the tool body length and the rake angle of the cutting edge, map these parameters to the specific geometric features in the optimized three-dimensional reamer attribute model. For example, the tool body length corresponds to the axial length of the tool body part in the three-dimensional model, and the rake angle of the cutting edge corresponds to the angle between the cutting edge and the tool body axis. Using the parametric modeling function of CAD software (such as SolidWorks), parameterize these key geometric features to make them variable geometric structure variables. For example, the tool body length can be defined as a variable parameter L, and its value range is from 50 mm to 70 mm; the rake angle of the cutting edge can be defined as a variable parameter α, and its value range is from 10° to 20°. Store these variable geometric structure variables and their value ranges as the key structural variable data for reamer adjustment.

[0146] Step S33: Perform structural variable perturbation processing on the key structural variable data for reamer adjustment using the preset structural variable perturbation rules, and perform finite difference calculation to generate the structural variable sensitivity data.

[0147] In the embodiment of the present invention, using a preset structural variable perturbation rule, such as the central difference method, the key structural variable data of the reamer adjustment generated in step S32 is subjected to structural variable perturbation processing. For example, for the cutter body length L, the perturbation step size is set to 1 mm, and the corresponding reamer performance index values, such as the first natural frequency, are calculated when L + 1 mm and L - 1 mm respectively. Then, the structural variable sensitivity is calculated using the finite difference method. For example, the central difference formula is used to calculate the sensitivity of the cutter body length L to the first natural frequency: sensitivity = (frequency(L + 1 mm) - frequency(L - 1 mm)) / (2×1 mm). Perturbation and sensitivity calculation are performed for each key structural variable to generate structural variable sensitivity data.

[0148] Step S34: Perform three-dimensional grid node space mapping on the optimized three-dimensional reamer attribute model through the structural variable sensitivity data to obtain reamer structure sensitivity mapping data.

[0149] In the embodiment of the present invention, the perturbation amount and sensitivity value of each structural variable are obtained. Then, according to the reamer structure feature mapping relationship established in step S32, the grid area affected by each structural variable is determined. For example, the cutting edge width variable mainly affects the grids near the cutting edge, and the chip flute depth variable mainly affects the grids near the chip flute. For each affected grid node, according to the perturbation amount of the structural variable, its corresponding spatial displacement amount is calculated. For example, for the case where the cutting edge width increases by 0.005 mm, the grid nodes near the cutting edge are translated 0.005 mm along the direction of the cutting edge width. The specific calculation of the displacement vector needs to be calculated according to the spatial position of the node and the perturbation direction of the structural variable. If the structural variable is an angle, the grid nodes are rotated according to the angle change to ensure that the mapped grid does not distort. After completing the spatial mapping of all grid nodes, reamer structure sensitivity mapping data is generated, which includes the original coordinates, mapped coordinates, and corresponding structural variable information of each grid node. For example, for grid node 10, the original coordinates are (x1, y1, z1), the perturbed coordinates are (x1', y1', z1'), and the corresponding structural variable is the cutting edge width; for grid node 20, the original coordinates are (x2, y2, z2), the perturbed coordinates are (x2', y2', z2'), and the corresponding structural variable is the chip flute depth. The structural variable sensitivity information is associated with the geometric structure of the three-dimensional model to obtain reamer structure sensitivity mapping data.

[0150] As an example of the present invention, refer to Figure 3 shown, which is Figure 1 a schematic diagram of the detailed implementation step flow of step S4 in

[0151] Step S41: Identify the high-sensitivity region of the reamer structure by using a preset sensitivity threshold for the reamer structure sensitivity mapping data, and obtain the high-sensitivity region of the reamer;

[0152] In the embodiment of the present invention, the reamer structure sensitivity mapping data is read, which contains the original coordinates, the mapped coordinates, and the corresponding structural variable information of each grid node. Then, the sensitivity value of each grid node is calculated. The sensitivity of a node can be calculated according to the change in node displacement before and after mapping. The change in node displacement is related to the perturbation of the corresponding structural variable. The calculation formula is, for example: sensitivity = |mapped node coordinate - original node coordinate| / perturbation of structural variable. A preset sensitivity threshold is set, for example, 0.8. Traverse the sensitivity values of all grid nodes, and mark the grid nodes with sensitivity values greater than this threshold as high-sensitivity nodes. These high-sensitivity nodes constitute the high-sensitivity region of the reamer, and their grid indices and coordinate information are stored as the reamer high-sensitivity region data.

[0153] Step S42: Calculate the minimum response amplitude frequency based on the high-sensitivity region of the reamer, and generate the minimum target response frequency value;

[0154] In the embodiment of the present invention, modal analysis is performed on the optimized three-dimensional reamer property model to calculate the first several natural frequencies and vibration modes of the reamer. The modal analysis results include the vibration mode patterns corresponding to each natural frequency. For each natural frequency, calculate the average amplitude of the high-sensitivity region. The calculation method is, for example: for all grid nodes in the high-sensitivity region, calculate their amplitude values at this frequency, and then calculate the average amplitude. By calculating the average amplitude corresponding to all natural frequencies, find the natural frequency corresponding to the minimum average amplitude. This natural frequency is defined as the minimum response amplitude frequency. This frequency will be used as the target frequency for subsequent structural responses.

[0155] Step S43: Set the structural response target for the minimum target response frequency value based on the high-sensitivity region of the reamer, and obtain the reamer structural response type data;

[0156] In the embodiments of the present invention, according to the design objective, the type of structural response is determined. Since the minimum response amplitude frequency calculated in step S42 is obtained, and it is desired to optimize the reamer so that its response at this frequency is smaller, the structural response type is defined as "minimizing the response amplitude". At the same time, the minimum target response frequency value is used as the target frequency of the structural response, clarifying that the goal of subsequent optimization is to minimize the response amplitude at this frequency. In addition to the response target, it is also necessary to set the constraint conditions in the optimization process. For example, volume constraints, mass constraints, etc. of the reamer can be set, and these constraint conditions will limit the distribution of materials in the optimization process. For example, it can be set that the volume of the reamer cannot exceed 105% of the initial volume. The structural response target type, the minimum target response frequency value, and the relevant constraint conditions are stored as reamer structural response type data. For example, the data records that the structural response type is "minimizing the response amplitude", the minimum target response frequency value is 3200 Hz, and the volume constraint is 105% of the initial volume.

[0157] Step S44: Set the response optimization expression for the high-sensitivity region of the reamer through the minimum target response frequency value and the reamer structural response type data to generate the reamer structural response optimization expression;

[0158] In the embodiments of the present invention, according to the structural response type, a suitable optimization objective function is selected. Since the response type is "minimizing the response amplitude", the objective function is selected as minimizing the average response amplitude of the high-sensitivity region at the target frequency. This average response amplitude can be calculated through finite element analysis. To meet the constraint conditions, a penalty term needs to be introduced into the optimization objective function. For example, if the volume of the reamer exceeds the preset volume constraint, a penalty term is added to the objective function to reduce the tendency of the volume to increase. The response optimization expression needs to include optimization variables. For example, the material density of the mesh elements. By changing the material density, the distribution of materials in the reamer structure can be controlled, thereby changing the response characteristics of the structure. The form of the optimization expression is as follows: Minimize J = average response amplitude + penalty term, where the average response amplitude is the average response value of the mesh nodes in the high-sensitivity region at the target frequency, and the penalty term is a term defined according to the constraint conditions. If the volume constraint is exceeded, the penalty term will increase, guiding the optimization process to reduce the excess volume constraint. The reamer structural response optimization expression is stored in the form of a mathematical formula, including the objective function, the constraint conditions, and the optimization variables.

[0159] Step S45: Process the volume constraint of the optimized three-dimensional reamer attribute model using the high-sensitivity region of the reamer to generate the reamer volume constraint condition data; obtain the reamer manufacturing constraint condition data;

[0160] In an embodiment of the present invention, a volume constraint condition is set. For example, the maximum volume increase allowed for the reamer is 5%. The initial volume and the maximum volume increase are stored as reamer volume constraint condition data, which includes the initial volume value, the maximum volume limit value, and the corresponding unit in mm 3 . For example, the initial volume is 2500 mm 3 , and the maximum volume constraint is 2625 mm 3 . While obtaining the reamer volume constraint condition data, it is also necessary to obtain the reamer manufacturing constraint condition data. The manufacturing constraint conditions are determined according to the actual manufacturing process of the reamer. For example, if machining is used, the minimum cutting size limit needs to be satisfied; if additive manufacturing is used, the minimum wall thickness limit, the maximum overhang angle limit, etc. need to be satisfied. These manufacturing constraint conditions will limit the optimized structural design to ensure its manufacturability. For example, the minimum cutting size limit is 0.2 mm, the minimum wall thickness limit is 0.5 mm, and the maximum overhang angle is 45 degrees. The reamer manufacturing constraint condition data is stored in the form of a list, which includes the type and value of each manufacturing constraint.

[0161] Step S46: Based on the reamer manufacturing constraint condition data and the reamer volume constraint condition data, use the reamer structure response optimization expression to perform finite element material distribution optimization on the optimized three-dimensional reamer property model, and perform three-dimensional simulation model reconstruction to generate optimized reamer three-dimensional simulation model data.

[0162] In an embodiment of the present invention, the reamer structure response optimization expression, the reamer manufacturing constraint condition data, and the reamer volume constraint condition data are imported into the finite element analysis optimization software. Using the topology optimization method, perform material distribution optimization on the optimized three-dimensional reamer property model. According to the optimization expression, iteratively adjust the material density of each grid cell. During the optimization process, the software continuously adjusts the material density of each grid cell according to the objective function. For regions where the volume exceeds the constraint conditions, the material density will be reduced; for regions with high sensitivity, the material density will be adjusted to minimize the response of the reamer at the target frequency. At the same time, the optimization process will satisfy the manufacturing constraint conditions to avoid generating unmanufacturable structures. After multiple iterations, the material distribution gradually converges to obtain the optimal material distribution scheme. According to the optimized material distribution scheme, reconstruct the three-dimensional simulation model. This process remaps the optimized material density information back to the three-dimensional grid model to generate optimized reamer three-dimensional simulation model data, which includes the geometric information, material distribution information, and manufacturing constraint information of the reamer.

[0163] Preferably, step S45 includes the following steps:

[0164] Step S451: Evaluate the influence degree of material change on the optimized three-dimensional reamer property model through the high-sensitivity region of the reamer to obtain structural region influence degree data;

[0165] Step S452: Extract the initial volume of the structure according to the optimized three-dimensional reamer attribute model to obtain the optimized initial structure volume value;

[0166] Step S453: Analyze the volume reduction ratio of the optimized initial structure volume value through the structure region influence degree data to generate region volume reduction ratio data;

[0167] Step S454: Use the region volume reduction ratio data to calculate the maximum material topology removal volume for the high-sensitivity region of the reamer to generate region maximum removal volume data;

[0168] Step S455: Perform structural volume constraint processing according to the region maximum removal volume data to obtain the reamer volume constraint condition data.

[0169] In the embodiment of the present invention, within the high-sensitivity region, a small change (such as a 5% reduction) in the simulated material property (such as elastic modulus) is made, and finite element analysis is performed to observe its influence on the target performance (such as the minimum target response frequency value). The influence degree is quantified, for example, expressed as a percentage of performance change, to generate the structure region influence degree data. The influence of material changes in different regions on the performance is different, and the influence degree data reflects the sensitivity of different regions to material changes. The total volume of the directly measured model is measured, or the model is exported in STL format, and then its volume value is calculated using volume calculation software. The extracted volume value is recorded as the optimized initial structure volume value, for example, the initial volume is 10000 mm 3 . The high-sensitivity region is divided into several sub-regions, and according to the influence degree data of each sub-region, a reasonable volume reduction ratio is determined. For regions with a high influence degree, a larger volume reduction ratio is allowed; for regions with a low influence degree, the volume reduction ratio needs to be restricted to avoid excessive weakening of the structural strength. For example, the high-sensitivity region is divided into three sub-regions with influence degrees of 0.9, 0.7, and 0.5 respectively, and the corresponding volume reduction ratios are set to 20%, 15%, and 10% respectively to generate region volume reduction ratio data. Multiply the initial volume of each sub-region by the corresponding volume reduction ratio to obtain the maximum volume that can be removed from the sub-region. For example, the initial volume of a certain sub-region is 2000 mm 3 and the volume reduction ratio is 20%, then the maximum volume that can be removed from this sub-region is 400 mm 3 . Summarize the maximum removal volumes of all sub-regions to generate region maximum removal volume data. Use the region maximum removal volume value as the upper limit of the volume constraint, for example: total removal volume <= sum of the maximum removal volumes of all sub-regions. This means that during the topology optimization process, the volume of the removed material cannot exceed this upper limit value.

[0170] This application aims to accurately digitize the geometric shape of a reamer by performing multi-view 3D scanning on the target reamer to obtain the original point cloud data of the reamer assembly and reconstructing the initial 3D reamer geometric model based on this data. Compared with traditional manual modeling or simplified modeling methods, 3D scanning technology can more realistically and comprehensively reflect the actual structure of the reamer assembly, avoiding errors introduced by simplification or approximation and ensuring the accuracy and reliability of the design results. Based on the initial 3D model, its applicability in numerical simulation is ensured. Mesh generation is a key step in numerical simulation methods such as finite element analysis. Reasonable mesh generation can improve the accuracy and efficiency of simulation calculations. Through mesh generation, the complex geometric model is discretized into a finite number of elements. By identifying components from the mesh data of the reamer 3D model, the complex reamer structure is decomposed into independent components, and corresponding physical properties are assigned to each component. This decomposition and property setting enable the simulation analysis to more precisely simulate the real behavior of the reamer during the cutting process, avoiding errors caused by treating the entire reamer as a homogeneous body, and particularly accurately considering the influence of different materials and connection structures on the overall dynamic performance. For example, the cutter body and the tool shank may be made of different materials with different elastic moduli and densities. Through component identification and property setting, simulation analysis can be performed on these components made of different materials respectively, so as to more accurately predict the stress distribution, deformation conditions, and vibration characteristics of the reamer under cutting loads. Through simulation of cutting loads, the dynamic response of the reamer under cutting conditions can be obtained, and the dynamic response of the reamer under different cutting parameters, different materials, and different structures can be simulated, and the natural frequency, vibration mode, and the relationship with the cutting frequency of the reamer can be analyzed, providing a scientific basis for suppressing chatter and optimizing the reamer design. For example, through simulation analysis, it can be found that when the natural frequency of the reamer is close to the cutting frequency, resonance is likely to occur, resulting in defects such as ripples and scratches on the machined surface, seriously affecting the machining accuracy and surface quality. According to the simulation results, the structural parameters of the reamer, such as the cutter body length and the tool shank diameter, can be adjusted specifically to change the natural frequency of the reamer and avoid the resonance frequency range, thus effectively suppressing the occurrence of chatter. Through reasonable allocation of material distribution, the overall performance of the reamer is improved. For example, increasing the material density or changing the material properties in high-sensitivity areas can improve the stiffness and vibration resistance of the reamer in these areas, thereby improving the cutting performance of the reamer.

[0171] Therefore, in any aspect, the embodiments should be regarded as exemplary and non-limiting. The scope of the present invention is defined by the appended claims rather than the above description. Therefore, all changes falling within the meaning and scope of the equivalent elements of the application documents are intended to be encompassed within the present invention.

[0172] The above are only specific embodiments of the present invention, enabling those skilled in the art to understand or implement the present invention. Various modifications to these embodiments will be obvious to those skilled in the art, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to these embodiments shown herein, but rather will be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. An integrated design method based on a three-dimensional simulation model of a reamer assembly, characterized in that: The following steps are involved: Step S1: Perform multi-view 3D scanning on the target reamer to obtain original point cloud data of the reamer component; reconstruct a 3D reamer model according to the original point cloud data of the reamer component to generate an initial 3D reamer geometric model; perform geometric model meshing according to the initial 3D reamer geometric model to generate 3D model mesh data of the reamer; Step S2: Reamer component identification is performed according to the reamer three-dimensional model mesh data to obtain reamer component identification data; model attributes of the initial three-dimensional reamer geometric model are set according to the reamer component identification data to generate an optimized three-dimensional reamer attribute model; cutting load simulation is performed according to the optimized three-dimensional reamer attribute model to generate reamer simulation modal analysis data; Step S3: Perform structural geometric parameter influence analysis on the optimized three-dimensional reamer attribute model through reamer simulation modal analysis data to generate key influencing parameters of the reamer; perform structural variable disturbance processing according to the key influencing parameters of the reamer, and perform three-dimensional grid node space mapping to obtain reamer structural sensitivity mapping data; Step S4: Utilize the reamer structure sensitivity mapping data to perform finite element material distribution optimization on the optimized three-dimensional reamer property model, and perform three-dimensional simulation model reconstruction to generate optimized reamer three-dimensional simulation model data.

2. The integrated design method based on the three-dimensional simulation model of the reamer assembly according to claim 1 is characterized in that: Step S1 includes the following steps: Step S11: performing surface cleaning treatment on the target reamer to obtain a cleaned target reamer component sample; Step S12: performing a multi-view three-dimensional scan on the cleaning target reamer component sample to obtain original point cloud data of the reamer component; Step S13: performing point cloud data filtering processing according to the original point cloud data of the reamer assembly, and reconstructing a three-dimensional reamer model to generate an initial three-dimensional reamer geometric model; Step S14: performing a blade area geometric feature line analysis based on the initial three-dimensional reamer geometric model to generate component blade area feature data; Step S15: Setting blade size control parameters for the initial three-dimensional reamer geometric model based on the component blade area feature data to generate blade size control parameters; Step S16: meshing the initial three-dimensional reamer geometric model using blade size control parameters to generate reamer three-dimensional model mesh data.

3. The integrated design method based on the three-dimensional simulation model of the reamer assembly according to claim 2 is characterized in that: Step S14 includes the following steps: Step S141: performing blade area recognition on the initial three-dimensional reamer geometric model to generate reamer blade area data; Step S142: performing edge contour detection on the reamer blade area data and extracting geometric feature lines to obtain blade feature line data; Step S143: performing blade surface curvature analysis according to blade characteristic line data to generate blade line curvature data; Step S144: marking blade feature key points of the reamer blade area data based on the blade line curvature data to obtain candidate blade feature point data; Step S145: Connect blade line segments according to the candidate blade feature point data, and perform geometric feature parameterization processing to generate component blade area feature data.

4. The integrated design method based on the three-dimensional simulation model of the reamer assembly according to claim 1 is characterized in that: Step S2 includes the following steps: Step S21: Reamer component identification is performed according to the reamer three-dimensional model mesh data to obtain reamer component identification data; Step S22: setting model attributes of the initial three-dimensional reamer geometric model through the reamer component identification data to generate an optimized three-dimensional reamer attribute model; Step S23: performing finite element processing on the optimized three-dimensional reamer property model to obtain a finite element three-dimensional reamer simulation model; Step S24: performing cutting load simulation according to the finite element three-dimensional reamer simulation model, and performing multi-order natural frequency / vibration mode analysis to generate reamer simulation modal analysis data.

5. The integrated design method based on the three-dimensional simulation model of the reamer assembly according to claim 4 is characterized in that: Step S22 includes the following steps: Step S221: querying the component material attributes of the reamer component identification data through a preset reamer component material attribute library to generate component material attribute data; Step S222: mapping the component material attribute data to the initial three-dimensional reamer geometric model to obtain a three-dimensional reamer attribute model; Step S223: performing component contact relationship analysis on the reamer component identification data through the three-dimensional reamer attribute model to generate component contact relationship data; Step S224: setting contact constraint parameters for the three-dimensional reamer attribute model based on the component contact relationship data to generate component contact constraint condition data; Step S225: setting model constraints on the three-dimensional reamer property model through the component contact constraint condition data to generate an optimized three-dimensional reamer property model.

6. The integrated design method based on the three-dimensional simulation model of the reamer assembly according to claim 5 is characterized in that: Step S223 includes the following steps: Perform component mesh marking on reamer component identification data through a three-dimensional reamer attribute model to obtain reamer component mesh marking data; Perform geometric proximity search based on the mesh tag data of the reamer component, process the component surface mesh nodes, and generate the component surface mesh node data; Calculate the spatial node distance of the component surface according to the component surface mesh node data to generate the spatial node distance data of the component surface; The contact surface node pairs are screened for the component surface mesh node data by using the component surface spatial node distance data to obtain potential component contact surface data; The potential component contact surface data is used to analyze the actual reamer assembly relationship of the three-dimensional reamer attribute model, and the actual reamer assembly data is generated; The contact surface data of potential components are processed in the form of relative motion of the contact surface through the actual reamer assembly data, and the contact relationship type is determined to generate component contact relationship data.

7. The integrated design method based on the three-dimensional simulation model of the reamer assembly according to claim 4 is characterized in that: Step S24 includes the following steps: Step S241: obtaining the actual operating frequency range of the target reamer; Step S242: Processing the resonance natural frequency value according to the actual working frequency range of the target reamer, and performing mode vector analysis to generate natural frequency / mode vector data; Step S243: performing modal mode / resonance risk mode processing according to the natural frequency / mode vector data, and performing multi-order modal parameter setting to generate a mode mode / resonance multi-order modal parameter list; Step S244: performing energy contribution analysis according to the vibration mode / resonance multi-order modal parameter list, and superimposing target modal simulation parameters to obtain modal superposition simulation parameters; Step S245: performing cutting load simulation using a finite element three-dimensional reamer simulation model based on modal superposition simulation parameters to obtain dynamic cutting load response data; Step S246: performing multi-order modal response superposition according to the dynamic cutting load response data, and performing frequency domain response analysis of the reamer part, thereby obtaining reamer simulation modal analysis data.

8. The integrated design method based on the three-dimensional simulation model of the reamer assembly according to claim 1 is characterized in that: Step S3 includes the following steps: Step S31: performing structural geometric parameter influence analysis on the optimized three-dimensional reamer attribute model through reamer simulation modal analysis data to generate key influencing parameters of the reamer; Step S32: mapping the structural characteristics of the reamer according to the key influencing parameters of the reamer, and performing variable processing on the variable geometric structure to generate key structural variable data for adjusting the reamer; Step S33: using preset structural variable disturbance rules to perform structural variable disturbance processing on the key structural variable data of the reamer adjustment, and performing finite difference calculation to generate structural variable sensitivity data; Step S34: Perform three-dimensional grid node space mapping on the optimized three-dimensional reamer attribute model through the structural variable sensitivity data to obtain reamer structure sensitivity mapping data.

9. The integrated design method based on the three-dimensional simulation model of the reamer assembly according to claim 1 is characterized in that: Step S4 includes the following steps: Step S41: identifying the reamer high sensitivity area on the reamer structure sensitivity mapping data by using a preset sensitivity threshold to obtain the reamer high sensitivity area; Step S42: Calculate the minimum response amplitude frequency according to the high sensitivity area of ​​the reamer to generate a minimum target response frequency value; Step S43: setting a structural response target for the minimum target response frequency value based on the reamer high sensitivity area to obtain reamer structural response type data; Step S44: setting a response optimization expression for the reamer high sensitivity region according to the minimum target response frequency value and the reamer structure response type data, and generating a reamer structure response optimization expression; Step S45: using the reamer high sensitivity area to perform volume constraint processing on the optimized three-dimensional reamer attribute model to generate reamer volume constraint condition data; and obtaining reamer manufacturing constraint condition data; Step S46: Based on the reamer manufacturing constraint data and the reamer volume constraint data, the optimized three-dimensional reamer property model is optimized by finite element material distribution optimization using the reamer structure response optimization expression, and the three-dimensional simulation model is reconstructed to generate optimized reamer three-dimensional simulation model data.

10. The integrated design method based on the three-dimensional simulation model of the reamer assembly according to claim 9, characterized in that: Step S45 includes the following steps: Step S451: evaluating the influence of material changes on the optimized three-dimensional reamer property model through the reamer high sensitivity area to obtain structural area influence data; Step S452: extracting the initial volume of the structure according to the optimized three-dimensional reamer attribute model to obtain the initial volume value of the optimized structure; Step S453: performing volume reduction ratio analysis on the initial volume value of the optimized structure through the structural region influence data to generate regional volume reduction ratio data; Step S454: Calculate the maximum material topological removal volume of the reamer high sensitivity area using the regional volume reduction ratio data to generate regional maximum removal volume data; Step S455: Perform structural volume constraint processing according to the regional maximum removal volume data to obtain reamer volume constraint condition data.

Images