Special-shaped plate stamping die reset process simulation method and system combined with data modeling

Abstract

The invention provides a special-shaped plate stamping die reset process simulation method and system combined with data modeling, and relates to the technical field of computer modeling. Firstly, a special-shaped plate reset dynamic topology model is constructed based on special-shaped plate three-dimensional shape data and die reset component motion feature data; associating the contact area with the acting force transmission path through the topological node; then, shape deformation data and mold reset force real-time data in the reset process are collected to generate a shape force coupling data chain; inputting the morphological force coupling data chain into a special-shaped plate reset dynamic topology model to drive topology nodes to dynamically recombine and generate topology recombination data; generating dynamic adaptive reset path data based on the topology recombination data; and finally, reversely inputting the dynamic adaptive reset path data into the special-shaped plate reset dynamic topology model for iterative optimization, and generating simulation complete data. According to the invention, the reset process of the special-shaped plate stamping die can be accurately simulated, and the simulation accuracy and reliability are improved.

Description

Technical Field

[0001] This invention relates to the field of computer modeling technology, and more specifically, to a method and system for simulating the resetting process of irregular plate stamping dies by combining data modeling. Background Technology

[0002] In the production and application of irregular-shaped sheet stamping dies, the reset process after stamping directly affects the quality and production efficiency of the stamped products. Traditional research on the reset process of irregular-shaped sheet stamping dies has focused on the mechanical motion principles of die components and simple mechanical analysis. For example, theoretical calculations and empirical formulas are used to estimate the force required for die reset and the motion trajectory of components. However, these methods often neglect the influence of the complex morphological characteristics of the irregular-shaped sheet itself on the reset process.

[0003] Furthermore, while existing technologies employ finite element analysis (FEM) to simulate the mold reset process, FEM requires the construction of a fine mesh model. For complex and irregularly shaped components like the aforementioned irregularly shaped plates, mesh generation is challenging, computationally intensive, and it's difficult to accurately capture the dynamic coupling relationship between the shape and force of the irregularly shaped plate during reset. Simultaneously, existing methods cannot reflect the synchronous changes between the deformation of the irregularly shaped plate and the mold reset force in real time, making it difficult to generate a precise dynamically adapted reset path. This results in low accuracy and reliability in simulating the mold reset process, failing to meet the demands of modern high-precision stamping production. Summary of the Invention

[0004] In view of the aforementioned problems, and in conjunction with the first aspect of the present invention, embodiments of the present invention provide a method for simulating the resetting process of a stamping die for irregularly shaped plates by incorporating data modeling, the method comprising: Based on the three-dimensional morphological data of the irregular plate and the motion characteristic data of the mold reset component, a dynamic topology model for irregular plate reset is constructed. The dynamic topology model for irregular plate reset is associated with the contact area of ​​the irregular plate and the force transmission path of the mold reset through topology nodes. Collect morphological deformation data and real-time mold reset force data during the resetting process of the irregular plate, and generate a morphological force coupling data chain. The morphological force coupling data chain records the synchronous change information of deformation and corresponding reset force in each contact area of ​​the irregular plate. The morphological force coupling data chain is input into the dynamic topology model of the irregular plate reset, which drives the dynamic topology model of the irregular plate reset to perform dynamic reorganization of topology nodes and generate topology reorganization data. The topology reorganization data reflects the node connection relationship after the deformation of the irregular plate and the reset force are adapted. Based on the topology reorganization data, dynamic adaptation reset path data is generated. The dynamic adaptation reset path data includes the motion trajectory and action sequence of each component of the mold after adjustment according to the topology node relationship. The dynamic adaptation reset path data is input in reverse into the dynamic topology model of the irregular plate reset, and the association strength of the topology nodes and the reset force transmission coefficient are iteratively optimized to generate the complete simulation data after iterative optimization.

[0005] Furthermore, embodiments of the present invention also provide a simulation system for the reset process of a stamping die for irregularly shaped plates, incorporating data modeling, characterized in that it includes: A processor; a machine-readable storage medium for storing machine-executable instructions of the processor; wherein the processor is configured to execute the above-described simulation method for resetting a non-standard stamping die incorporating data modeling by executing the machine-executable instructions.

[0006] In another aspect, embodiments of the present invention also provide a computer program product, the computer program product including machine-executable instructions stored in a computer-readable storage medium, wherein a processor of a simulation system for the resetting process of a shaped plate stamping die in conjunction with data modeling reads the machine-executable instructions from the computer-readable storage medium, and the processor executes the machine-executable instructions, causing the simulation system for the resetting process of a shaped plate stamping die in conjunction with data modeling to execute the above-described simulation method for the resetting process of a shaped plate stamping die in conjunction with data modeling.

[0007] Based on the above, a dynamic topology model for irregular plate resetting is constructed using 3D morphological data of the irregular plate and motion characteristic data of the mold resetting components. Topological nodes cleverly link the contact areas of the irregular plate with the force transmission path of the mold resetting action. Real-time data on the morphological deformation and mold resetting force during the resetting process are collected and a morphological-force coupling data chain is generated. This chain completely records the synchronous changes in deformation and corresponding resetting force in each contact area of ​​the irregular plate, reflecting the dynamic interaction between form and force during resetting. The morphological-force coupling data chain is input into the dynamic topology model to drive dynamic reorganization of topology nodes, generating topology reorganization data reflecting the node connection relationship after the irregular plate deformation and resetting force are adapted. This, in turn, generates dynamically adapted resetting path data containing the motion trajectory and action sequence of each mold component adjusted according to the topology node relationship, achieving precise planning of the mold resetting process. Finally, the dynamically adapted resetting path data is input back into the model for iterative optimization, improving the accuracy of the topology node association strength and the resetting force transmission coefficient. Ultimately, iteratively optimized complete simulation data is generated, significantly improving the accuracy, reliability, and comprehensiveness of the simulation of the irregular plate stamping mold resetting process. Attached Figure Description

[0008] Figure 1 This is a schematic diagram of the execution flow of the simulation method for resetting irregular plate stamping die by combining data modeling provided in the embodiments of the present invention.

[0009] Figure 2This is a schematic diagram of exemplary hardware and software components of a simulation system for the resetting process of a non-circular plate stamping die, which combines data modeling, provided in an embodiment of the present invention. Detailed Implementation

[0010] Figure 1 This is a flowchart illustrating a method for simulating the resetting process of a non-circular plate stamping die based on data modeling, as provided in one embodiment of the present invention. A detailed description follows.

[0011] Step S110: Based on the three-dimensional morphological data of the irregular plate and the motion characteristic data of the mold reset component, construct a dynamic topology model for the irregular plate reset. The dynamic topology model for the irregular plate reset is associated with the contact area of ​​the irregular plate and the force transmission path of the mold reset through topology nodes.

[0012] This embodiment uses the stamping die reset process of an irregularly shaped side panel of an automobile body as an application scenario. This irregularly shaped panel has complex curved surface features such as wheel arch protrusions, waistline recesses, and window frame edges. The die reset system includes multiple sets of wedge slider assemblies, return push rods, and nitrogen cylinders. During the reset stage after stamping, the irregularly shaped panel forms a dynamic contact relationship with the die components due to elastic deformation. A topological model is needed to map the contact area to the force transmission path. Specifically, the key areas on the surface of the irregularly shaped panel that may contact the die are first identified, and the areas within the die components where reset forces may be applied are determined. Then, by establishing the relationships between these areas, a topological structure reflecting the interaction between the two is constructed. This topological model includes not only the spatial location information of the nodes but also the force transmission rules between the nodes.

[0013] Step S111: Collect the three-dimensional contour data, thickness distribution data and material elasticity data of the irregular plate, and generate the three-dimensional morphological data of the irregular plate through three-dimensional modeling technology. The three-dimensional morphological data of the irregular plate is used to present the spatial structure and physical parameters of the irregular plate.

[0014] For irregularly shaped automotive side panels, a multi-source data acquisition scheme was adopted. A 3D optical scanning device was used to scan the entire surface of the panel, acquiring the spatial coordinates of various points on the surface to form 3D contour data. An ultrasonic thickness gauge was used to measure the thickness at multiple points on the panel surface, obtaining thickness distribution data. A materials testing machine was used to perform elasticity tests on the panel material, determining parameters such as the elastic modulus and Poisson's ratio. The above data was imported into professional 3D modeling software, and data fusion technology was used to integrate the 3D contour, thickness distribution, and material elastic parameters into a single model. During the modeling process, it was necessary to ensure that data from different sources remained consistent in spatial coordinates; coordinate transformation and registration techniques were used to eliminate systematic errors generated during data acquisition. The generated 3D morphological data of the irregularly shaped panel fully presented its spatial geometric structure and physical property parameters in digital form.

[0015] Step S1111: Use a 3D laser scanning device to perform a full-surface scan of the irregular plate, record the spatial coordinate data of each surface point of the irregular plate, and form the original 3D contour data.

[0016] The irregularly shaped side panel of the car was fixed to a measuring platform equipped with pneumatic clamps to ensure that the panel would not shift during scanning. A 3D laser scanning device was used, which emitted a laser beam onto the surface of the irregularly shaped panel and calculated the spatial coordinates of each surface point by receiving the reflected beam. The device was calibrated before scanning, including angle and distance calibration of the laser emitter, to ensure scanning accuracy. The scanning path followed a pre-planned trajectory, covering the entire surface of the irregularly shaped panel, including complex feature areas such as wheel arches and window frames. During scanning, the device automatically recorded the X, Y, and Z axis coordinates of each sampling point, and these coordinates were stored sequentially to form point cloud data. To improve data integrity, areas with significant curvature changes were scanned multiple times to ensure that no key geometric features were missed. The original 3D contour data contained a large amount of point coordinate information, which directly reflected the surface shape of the irregularly shaped panel.

[0017] Step S1112: Measure the thickness of different areas of the irregularly shaped plate using an ultrasonic thickness gauge, record the thickness values ​​according to the measurement location, and form the original thickness distribution data.

[0018] Based on the structural characteristics of the irregularly shaped side panel of the automobile, multiple measurement areas were divided on its surface. These areas were divided considering the geometric features of the panel and the stress conditions during subsequent stamping. Within each measurement area, multiple measurement points were selected, and the thickness was measured using an ultrasonic thickness gauge. During measurement, the probe was placed in close contact with the panel surface, and a coupling agent was applied to ensure good sound wave transmission. Multiple measurements were performed at each measurement point, and the average value of the results was taken as the thickness value at that point. The thickness values ​​were recorded sequentially according to the coordinates of the measurement points to form the original thickness distribution data. For areas that were difficult to measure directly, the probe angle was adjusted or a specially shaped probe was used to ensure the comprehensiveness of the thickness data. The original thickness distribution data reflects the thickness variation of the irregularly shaped panel at different locations.

[0019] Step S1113: Use a materials mechanics testing device to perform elasticity testing on the material of the irregular plate, record the elastic modulus data and Poisson's ratio data of the material, and form the original material elasticity data.

[0020] Standard tensile specimens were cut from the scrap area of ​​the irregularly shaped sheet, with the specimen shape and dimensions conforming to the material mechanics testing specifications. The specimens were mounted on the clamps of the material mechanics testing equipment, ensuring uniform stress during the tensile process. The loading rate of the testing equipment was set to apply a constant tensile force to the specimens. During the test, the elongation (longitudinal strain) and contraction (transverse strain) of the specimens perpendicular to the force direction were measured using an extensometer, while the magnitude of the applied tensile force was recorded. Based on the measured stress (tensile force divided by the specimen's cross-sectional area) and strain data, a stress-strain curve was plotted. The elastic modulus, i.e., the ratio of stress to strain, was determined from the linear phase of the curve. Poisson's ratio was calculated using the ratio of transverse strain to longitudinal strain. Multiple specimens were tested, and the elastic modulus and Poisson's ratio data obtained from each test were recorded to form the original material elasticity data. This data characterizes the elastic properties of the irregularly shaped sheet material and is a key parameter for mechanical analysis and simulation.

[0021] Step S1114: Denoise the original 3D contour data, remove abnormal point data generated during the scanning process, and form purified 3D contour data.

[0022] The original 3D contour data may contain outliers caused by factors such as scanning equipment errors, environmental interference, or surface reflections on the board. These outliers can affect the accuracy of subsequent modeling. A statistical filtering method is used to denoise the original data. First, for each data point, a certain number of neighboring points are identified, and the average coordinates of these neighboring points are calculated. Then, the distance between the data point and the average coordinates is calculated. If the distance exceeds a set threshold, the point is identified as an outlier and removed. The threshold is set based on the overall distribution of the data points and determined by analyzing the standard deviation of the data. After removing outliers, the data is smoothed using a weighted average method to adjust the coordinates of each point, making the point cloud data smoother while preserving the geometric features of the board. The processed and purified 3D contour data eliminates noise interference and more accurately reflects the actual surface shape of the irregularly shaped board.

[0023] Step S1115: Interpolate the original thickness distribution data to supplement the thickness values ​​of the measurement gap area and form continuous thickness distribution data.

[0024] The original thickness distribution data was acquired at discrete measurement points, with gaps existing between these points where the thickness values ​​were unknown. To obtain a continuous thickness distribution across the entire irregularly shaped plate surface, interpolation is required. A distance-weighted interpolation method is used to estimate the thickness of the gaps based on the known thickness values ​​of the measurement points. For each interpolation point, known measurement points within a certain radius are considered, and different weights are assigned based on their distance from the interpolation point, with closer points receiving greater weights. The thickness value of the interpolation point is then calculated using a weighted average. During interpolation, the number of neighboring measurement points is appropriately increased for the edge regions and areas with drastic feature changes in the irregularly shaped plate to improve interpolation accuracy. After interpolation, continuous thickness distribution data covering the entire irregularly shaped plate surface is generated. This continuous thickness distribution data is stored in a grid format, with each grid node corresponding to a thickness value, visually reflecting the spatial variation of the plate thickness.

[0025] Step S1116: Average the original material elasticity data to eliminate numerical fluctuations caused by experimental errors and form standard material elasticity data.

[0026] The original material elasticity data were obtained through testing multiple samples. Due to factors such as sample preparation and testing environment, the test results may fluctuate to some extent. To obtain stable and reliable material elasticity parameters, the original data needs to be mean-processed. The elastic modulus and Poisson's ratio data obtained from multiple tests are statistically analyzed, and the arithmetic mean of each parameter is calculated. Simultaneously, the standard deviation of the data is calculated to analyze the degree of data dispersion. If the deviation of a test data point from the mean exceeds a set range (e.g., twice the standard deviation), the data is considered an outlier, discarded, and the mean is recalculated. After mean-processing, standard material elasticity data is obtained, including the average elastic modulus and average Poisson's ratio. This data eliminates numerical fluctuations caused by experimental errors and can more accurately represent the elastic properties of the irregularly shaped plate material.

[0027] Step S1117: Import the purified 3D contour data, continuous thickness distribution data, and standard material elasticity data into the 3D modeling software.

[0028] Save the processed cleaned 3D contour data, continuous thickness distribution data, and standard material elasticity data in a format recognizable by the 3D modeling software. The cleaned 3D contour data is typically stored as a point cloud, continuous thickness distribution data can be stored as texture maps or mesh data, and standard material elasticity data is stored as a parameter file. Launch the 3D modeling software and import the above data into the software environment sequentially using the software's import function. During the import process, ensure that the coordinate systems of all data are consistent; use coordinate transformation tools to align data from different sources to the same reference coordinate system. After importing, check the data in the software to confirm its integrity and accuracy, ensuring that no data is lost or corrupted.

[0029] Step S1118: Construct a three-dimensional solid model of the irregular plate using three-dimensional modeling software, and assign thickness parameters and material elasticity parameters to the corresponding areas in the three-dimensional solid model.

[0030] In 3D modeling software, based on imported cleaned 3D contour data, a surface model of the irregularly shaped plate is constructed using surface fitting or mesh reconstruction techniques. Then, based on continuous thickness distribution data, the surface model is thickened to generate a 3D solid model with a certain thickness. During model construction, complex feature areas of the irregularly shaped plate, such as wheel arches and window frame edges, are meticulously modeled to ensure the model accurately reflects the actual geometry. Next, material elastic parameters are assigned to the 3D solid model, allocating parameters such as the elastic modulus and Poisson's ratio from standard material elasticity data to the model. A material property library is created in the software, associating material parameters with the model to give it actual physical properties. The software's rendering function allows for a visual inspection of the model's geometry and material distribution.

[0031] Step S1119: Perform mesh generation on the 3D solid model to form mesh model data containing spatial coordinates and physical parameters.

[0032] To facilitate subsequent numerical simulation analysis, the 3D solid model needs to be meshed. In the 3D modeling software, a suitable mesh type (such as tetrahedral or hexahedral mesh) and mesh size are selected. The choice of mesh size needs to comprehensively consider the complexity of the model and the required computational accuracy. For areas with complex geometric features, a smaller mesh size is used to improve computational accuracy; for relatively simple areas, the mesh size can be appropriately increased to reduce computational load. During meshing, mesh quality control parameters, such as mesh distortion and aspect ratio, are set to ensure that the generated mesh quality meets the requirements of numerical computation. After meshing, each mesh element contains spatial coordinate information and corresponding physical parameters (such as thickness and elastic modulus). The mesh model data is stored in a specific format, containing node information, element information, and the association information between elements and physical parameters.

[0033] Step S11110: Extract spatial structure information and physical parameter information from the mesh model data, integrate them to generate three-dimensional morphological data of the irregular plate, and the three-dimensional morphological data of the irregular plate fully presents the spatial structure and physical parameters of the irregular plate.

[0034] Spatial coordinate information of all nodes is extracted from the mesh model data. This coordinate information describes the spatial structure of the irregularly shaped plate. Simultaneously, physical parameter information, such as thickness, elastic modulus, and Poisson's ratio, is extracted from each mesh cell. The node coordinate information and physical parameter information are organized and correlated according to certain rules to ensure that each physical parameter accurately corresponds to its corresponding position in the model. Through data integration, the above information is merged into a unified data structure, forming the three-dimensional morphological data of the irregularly shaped plate. This three-dimensional morphological data of the irregularly shaped plate not only contains its geometric shape information but also its material physical properties information, comprehensively presenting the spatial structure and physical parameters of the irregularly shaped plate.

[0035] Step S112: Collect motion stroke data, rotation angle data, and component surface contact parameter data of each reset component of the mold to form motion characteristic data of the mold reset component. The motion characteristic data of the mold reset component is used to reflect the motion capability and contact parameter information of the reset component.

[0036] For each reset component of the mold, such as the wedge slider and return ejector, corresponding sensors are installed for data acquisition. For linearly moving components, displacement sensors are used to measure their stroke and record the positional changes during the reset process. For rotary components, angle sensors are used to measure their rotation angle and obtain their rotational attitude information. Simultaneously, contact parameters of the component surfaces, including surface roughness and coefficient of friction, are collected. Surface roughness is measured using a surface roughness meter to obtain the microscopic morphology parameters of the component surface. The coefficient of friction is determined using a dedicated friction testing device to simulate the contact friction of the components under actual working conditions. The collected stroke data, rotation angle data, and surface contact parameter data are categorized and organized by component to form the mold reset component motion characteristic data. This mold reset component motion characteristic data reflects the motion range, motion accuracy, and contact characteristics of each reset component.

[0037] Step S113: Extract the data of the protruding area, the concave area, and the edge contour from the three-dimensional shape data of the irregular plate. Determine the surface area of ​​the irregular plate corresponding to the data of the protruding area, the concave area, and the edge contour from the three-dimensional shape data of the irregular plate as the candidate area of ​​the topology node. Each candidate area of ​​the topology node corresponds to the key position where the irregular plate may contact the mold.

[0038] Feature analysis was performed on the 3D morphological data of the irregularly shaped plate to identify convex, concave, and edge contours on the surface. By calculating the curvature values ​​at various points on the model surface, areas with curvature values ​​greater than a set threshold were identified as convex areas, and areas with curvature values ​​less than a negative threshold were identified as concave areas. Edge contours were obtained by extracting the boundary lines of the model surface; these boundary lines are typically the edges of the plate or the intersections of different curved surfaces. The surface areas corresponding to the identified convex, concave, and edge contours were determined as candidate regions for topology nodes. When selecting candidate regions, the area size and geometric features of the regions were considered, prioritizing regions with larger areas and obvious geometric features, as these regions are more likely to come into contact with the mold components during the resetting process. Each candidate region for topology nodes was represented by the spatial coordinates of its geometric center for subsequent association and matching.

[0039] Step S114: Extract the component contact end face data, motion direction data and force output data from the motion characteristic data of the mold reset component. Determine the component action area corresponding to the component contact end face data, motion direction data and force output data from the motion characteristic data of the mold reset component as the topology node matching area. Each topology node matching area corresponds to the key position where the mold may apply a reset force.

[0040] Information related to contact is extracted from the motion characteristic data of the mold reset components. The component contact end face data includes information such as the geometry, dimensions, and spatial position of the contact end face, obtained through analysis of the mold's 3D model. Motion direction data represents the direction vector of the reset component during movement, determined based on the component's motion mechanism and design parameters. Force output data reflects the magnitude and direction of the force that the component can apply, calculated based on the performance parameters of the component's drive device (e.g., a nitrogen cylinder). By comprehensively considering the above data, the key locations where each reset component may apply reset force are determined, i.e., the topology node matching region. The matching region is typically located on the contact end face of the component, and its position is determined based on the geometry and motion direction of the contact end face. Each topology node matching region is represented by the spatial coordinates of its center and associated with the corresponding motion direction and force parameters.

[0041] Step S115: Establish the initial association between the candidate region of the topology node and the matching region of the topology node, and allocate the initial association weight according to the overlap between the spatial position of the irregular plate and the motion trajectory of the mold to form the initial topology association matrix.

[0042] The candidate regions of the topological nodes of the irregularly shaped plate and the matching regions of the topological nodes of the mold are compared in spatial position on the same coordinate system. The spatial distance between each candidate region and each matching region is calculated; the smaller the distance, the closer they are in space, and the greater the probability of contact during the resetting process. The association weight is determined based on the spatial distance, with a larger weight for smaller distances, and the weight value ranges from 0 to 1. Simultaneously, considering the motion trajectory of the mold components, the association weight is appropriately increased for matching regions whose motion trajectories pass near the candidate regions. These association weights are organized in matrix form, with rows representing candidate regions of the topological nodes, columns representing matching regions of the topological nodes, and the matrix elements representing the corresponding association weights, forming an initial topological association matrix. This initial topological association matrix describes the initial degree of association between the candidate regions and the matching regions.

[0043] Step S116: Based on the initial topology association matrix, define the candidate region of each topology node and the corresponding matching region of the topology node as the basic topology node of the irregular plate reset dynamic topology model. Each basic topology node includes spatial coordinate parameters and initial association weight parameters.

[0044] Based on the initial topological correlation matrix, candidate regions and matching regions with correlation weights greater than a set threshold are selected, and each pair of regions is defined as a basic topological node. Each basic topological node contains the spatial coordinates of the candidate region, the spatial coordinates of the matching region, and the initial correlation weight between them. The spatial coordinate parameter is used to determine the node's position in three-dimensional space, and the initial correlation weight parameter reflects the initial correlation strength between the nodes. All basic topological nodes are numbered to facilitate subsequent identification and processing. These basic topological nodes are the basic units constituting the dynamic topological model for irregular plate repositioning, and their attributes and interrelationships will be dynamically updated in the model.

[0045] Step S117: Construct a force transmission path model between basic topology nodes, and set the path transmission coefficient according to the elastic data of the irregular plate material and the contact parameter data of the mold component. The path transmission coefficient is directly related to the force transmission efficiency of adjacent basic topology nodes.

[0046] Analyze the spatial relationships between basic topological nodes to determine potential force transmission paths. For spatially adjacent basic topological nodes, a potential force transmission path is considered to exist between them. Based on the material elasticity data (such as elastic modulus) of the irregular plate and the contact parameter data (such as friction coefficient) of the mold components, a transmission coefficient is set for each path. A higher elastic modulus indicates less material deformation, higher force transmission efficiency, and a larger path transmission coefficient; a higher friction coefficient indicates greater force loss during transmission, and a smaller path transmission coefficient. The path transmission coefficient ranges from 0 to 1, directly reflecting the force transmission efficiency between adjacent basic topological nodes. Integrate the above paths and their transmission coefficients into the model to form a force transmission path model, which describes the transmission method and efficiency of forces between different nodes.

[0047] Step S118: Integrate the basic topology nodes, the initial topology association matrix, and the force transmission path model to form the core structure of the irregular plate reset dynamic topology model. The core structure of the irregular plate reset dynamic topology model is used to define the node composition and association rules.

[0048] The basic topology nodes, initial topology association matrix, and force transmission path model are integrated according to a certain logical relationship. The basic topology nodes serve as the fundamental elements of the model, the initial topology association matrix describes the initial association strength between nodes, and the force transmission path model specifies the rules for force transmission between nodes. During the integration process, the attributes of the nodes (spatial coordinates, association weights), the association relationships between nodes (represented by the association matrix), and the force transmission paths and coefficients are clearly defined. The resulting core structure defines the basic framework of the dynamic topology model, including the way nodes are constructed and the association rules between nodes.

[0049] Step S119: Input the material elasticity parameters in the three-dimensional morphology data of the irregular plate and the force output parameters in the motion characteristic data of the mold reset component, fill in the node physical properties and path transmission parameters in the core structure of the irregular plate reset dynamic topology model, and form a preliminary irregular plate reset dynamic topology model.

[0050] Material elastic parameters, such as elastic modulus and Poisson's ratio, are extracted from the 3D morphological data of the irregularly shaped plate and assigned to the basic topological nodes as their physical properties. Simultaneously, force output parameters, such as maximum force and force direction, are extracted from the motion characteristic data of the mold reset components and associated with the corresponding topological node matching regions. Based on the physical properties of the nodes and the force parameters, the path transmission coefficients in the force transmission path model are adjusted to reflect the actual force transmission. By filling these parameters into the core structure, the model acquires practical physical meaning and mechanical properties, forming a preliminary dynamic topological model for irregularly shaped plate reset. This preliminary model can simulate the basic interactions between nodes during the reset process.

[0051] Step S1110: Using the static contact test data between the irregular plate and the mold, adjust the association weights of the basic topology nodes and the force transmission path coefficients in the preliminary dynamic topology model for irregular plate resetting, and generate the final dynamic topology model for irregular plate resetting.

[0052] Static contact tests were conducted between the irregularly shaped plate and the mold. The plate was placed in the mold, and pressure was applied. The force and deformation data at different contact points were measured. The experimental data were compared with the simulation results of the preliminary dynamic topology model to analyze the differences. Based on the differences, the association weights of the basic topology nodes and the force transmission path coefficients in the model were adjusted. For contact areas with large forces in the experiment, the association weights of the corresponding nodes were increased; for paths whose force transmission efficiency did not match the experimental results, their transmission coefficients were adjusted. Through multiple iterative adjustments, the simulation results of the model were made as consistent as possible with the experimental data. The final dynamic topology model for resetting the irregularly shaped plate more accurately reflects the actual contact relationship and force transmission.

[0053] Step S120: Collect morphological deformation data and real-time mold reset force data during the resetting process of the irregular plate, and generate a morphological force coupling data chain. The morphological force coupling data chain records the synchronous change information of deformation and corresponding reset force in each contact area of ​​the irregular plate.

[0054] In the resetting test of the stamping die for an automotive side panel, it is necessary to collect real-time data on the morphological deformation of the panel and the resetting force of the die. Miniature deformation sensors are installed in the candidate topological node areas of the panel, and pressure sensors are installed in the matching topological node areas of the die. The sensor signals are transmitted to a computer in real time via a data acquisition system. The acquisition system is set with a uniform sampling frequency to ensure that the deformation data and resetting force data are synchronized in time. During the resetting process, the sensors continuously collect data, recording the deformation of each contact area of ​​the panel and the corresponding die resetting force at each moment. This data is organized chronologically to form a morphological-force coupled data chain. Each data unit in this morphological-force coupled data chain contains the deformation data and resetting force data at the corresponding moment, as well as a timestamp and node identifier, enabling a complete record of the synchronous changes of both during the resetting process.

[0055] Step S121: Install miniature deformation sensors in the candidate regions of topological nodes on the surface of the irregular plate, and fix the sensor acquisition points according to the spatial distribution of each candidate region of topological nodes, so that each acquisition point corresponds precisely to a candidate region of topological nodes.

[0056] Based on the distribution of candidate regions for topological nodes on the irregularly shaped plate surface, a miniature deformation sensor is installed at the center of each candidate region. A specialized adhesive is used to install the sensors, ensuring a tight fit between the sensor and the plate surface and preventing relative displacement during deformation. Before installation, the candidate region surface is cleaned to remove oil and impurities to ensure a firm bond. Each sensor's acquisition point corresponds one-to-one with its respective candidate region for topological nodes, identified by a number. The sensor leads are connected to the data acquisition device via flexible wires, and the wire arrangement avoids interfering with the normal resetting movement of the irregularly shaped plate. After installation, the sensors are calibrated to ensure their measurement accuracy meets the experimental requirements.

[0057] Step S122: Install pressure sensors in the topology node matching area of ​​the mold, and arrange the sensor detection surfaces according to the force output direction of each topology node matching area so that the detection surfaces are completely aligned with the contact direction of the irregular plate.

[0058] For each topological node matching area of ​​the mold, the installation direction of the pressure sensor is determined according to its force output direction, ensuring that the sensor's detection surface is perpendicular to the contact direction of the irregular plate. The sensor is fixed to the mold using mechanical clamps, the design of which ensures that the sensor will not shift its position under pressure. During installation, a level and angle gauge are used to adjust the sensor's orientation, ensuring that the detection surface is accurately aligned with the center of the matching area. The sensor's signal line is connected to the data acquisition system through a wiring channel inside the mold, preventing damage to the signal line during mold movement. After installation, a load test is performed on the sensor to check if its output signal is normal, ensuring accurate measurement of the reset force.

[0059] Step S123: Synchronously start the acquisition program of the deformation sensor and the pressure sensor, set a uniform acquisition time interval, the acquisition time interval is set according to the mold reset movement speed, and capture synchronous data at each moment of movement.

[0060] Before conducting the reset test, the data acquisition programs for the deformation and pressure sensors were configured. A uniform data acquisition time interval was set, determined based on the speed of the mold's reset movement; the faster the movement, the shorter the acquisition time interval, ensuring data is captured at every instant of movement. The data acquisition system's synchronous trigger function enabled all sensors to begin acquiring data simultaneously. During acquisition, the system monitored the sensor's operating status in real time to ensure the continuity and accuracy of data acquisition. The acquired data was stored in real time on the computer's hard drive, with each data sample containing a corresponding timestamp for subsequent data synchronization processing.

[0061] Step S124: Collect real-time deformation data, deformation rate data and deformation recovery data of each candidate region of the topology node through deformation sensors to form a set of deformation data of the irregular plate shape. The deformation information of each candidate region of the topology node is arranged in the order of collection time in the set of deformation data of the irregular plate shape.

[0062] In the candidate topological nodes of the irregularly shaped side panel of the automobile, miniature deformation sensors installed in each area continuously collect data at a preset frequency. Before data acquisition, the sensors are zero-point calibrated, and the deformation in the initial state is set as the reference value. Real-time deformation data is acquired through strain gauges built into the sensors. When the irregularly shaped panel deforms, the change in the resistance value of the strain gauge is converted into an electrical signal output, which is then converted into digital deformation data after analog-to-digital conversion. Deformation rate data is calculated by subtracting the real-time deformation data from two consecutive acquisition cycles and then dividing by the acquisition time interval, reflecting the amount of change in deformation per unit time. Deformation recovery data is acquired after the mold reset action has completely stopped and is continuously recorded for a preset duration until the deformation tends to stabilize. During this process, the sensors continue to collect data at the original frequency to capture the dynamic changes in the elastic recovery process of the irregularly shaped panel.

[0063] Step S1241: Set the acquisition frequency of the deformation sensor to be consistent with that of the pressure sensor to keep the data acquisition synchronized.

[0064] Based on the reset speed of the irregularly shaped side panel of the automobile, the sampling frequency of the deformation sensor is set. This sampling frequency must be the same as the sampling frequency of the pressure sensor installed in the matching area of ​​the mold topology node to ensure a one-to-one correspondence between deformation data and reset force data in the time dimension. The clock synchronization function of the data acquisition system ensures that the sampling start time and sampling interval of all sensors are strictly consistent, avoiding data misalignment caused by time asynchrony. The sampling frequency setting must satisfy the Nyquist sampling theorem to ensure complete capture of the dynamic changes in deformation and force.

[0065] Step S1242: The initial morphological data of the candidate region of the topological node is captured by the detection unit of the deformation sensor as the deformation comparison benchmark data.

[0066] Before the mold begins its reset process, the deformation sensor's detection unit acquires the initial state data of the candidate topological node region. At this time, the irregularly shaped plate is in a free state, without any external force, and the deformation recorded by the sensor is zero. This initial shape data is stored as the baseline value, and all subsequent real-time deformation data are calculated using this as a reference: real-time deformation = current acquired value - initial baseline value. The acquisition of baseline data must be carried out under conditions of stable ambient temperature and no vibration interference to avoid external factors affecting the accuracy of the baseline value.

[0067] Step S1243: During the resetting process of the irregular plate, continuously capture the real-time morphological data of each candidate region of the topology node, calculate the difference with the deformation comparison benchmark data, and obtain the real-time deformation data.

[0068] After the mold initiates its reset action, the deformation sensor enters continuous acquisition mode. Within each sampling cycle, the detection unit acquires the current morphological data of the candidate region of the topological node, subtracts it from the stored deformation comparison benchmark data, and obtains the real-time deformation data for that moment. The sign of the real-time deformation data represents the deformation direction (tension is positive, compression is negative), and the absolute value represents the deformation magnitude. The data acquisition system independently records the real-time deformation of each sensor channel and adds timestamp information.

[0069] Step S1244: Based on the real-time deformation data collected in two consecutive transactions, calculate the deformation change value per unit time to obtain the deformation rate data.

[0070] For each candidate region of the topology node, real-time deformation data ε1 and ε2 are acquired at two consecutive adjacent acquisition times t1 and t2 (t2>t1). The deformation rate v is calculated using the formula v=(ε2-ε1) / (t2-t1), with units of deformation units / time units. When ε2>ε1, the deformation rate is positive, indicating that the deformation is increasing; when ε2<ε1, the deformation rate is negative, indicating that the deformation is decreasing. The deformation rate data reflects how quickly the deformation changes over time and is used to analyze the dynamic response characteristics of the irregularly shaped plate during the reset process.

[0071] Step S1245: After the reset action stops, continuously collect the morphological recovery data of the candidate region of the topology node, record the process data from the end of the reset to the morphological stabilization, and obtain the deformation recovery data.

[0072] After the mold reset action completely stops, the data acquisition system continues to collect deformation data for a preset duration. During this period, the irregularly shaped plate begins to recover due to elastic deformation, and the deformation recorded by the sensors gradually decreases and tends to stabilize. The deformation recovery data includes the real-time deformation during the recovery process, the deformation rate, and the residual deformation after stabilization. The acquisition duration is determined based on the elastic recovery characteristics of the irregularly shaped plate material, and typically needs to continue until the change in deformation is less than a set threshold over several consecutive sampling periods.

[0073] Step S1246: Assign a unique data identifier to each candidate region of the topology node, and associate the data identifier with the corresponding real-time deformation data, deformation rate data, and deformation recovery data.

[0074] Each candidate region of the topology node is assigned a unique numerical code as a data identifier based on its spatial location on the irregular plate. Each data identifier corresponds one-to-one with the sensor number of that region, ensuring accurate data attribution. During data storage, real-time deformation data, deformation rate data, and deformation recovery data are all bound to their corresponding data identifiers, forming a dataset at the node level. For example, the dataset for node N01 contains all deformation-related data for that node throughout the entire reset process.

[0075] Step S1247: Sort the real-time deformation data, deformation rate data and deformation recovery data corresponding to each data identifier according to the acquisition time order to form a deformation data sequence of a single node.

[0076] For each node corresponding to a data identifier, the real-time deformation data, deformation rate data, and deformation recovery data collected during the reset process are arranged in chronological order of the collection timestamps to form a one-dimensional time series. Each element in the data series contains three data items: timestamp, real-time deformation, and deformation rate (the data in the deformation recovery phase also contains these three items). The length of the time series is determined by the total collection duration and collection frequency, ensuring complete coverage of the entire reset process from the start to complete deformation stabilization.

[0077] Step S1248: Collect the individual node deformation data sequence of all candidate topological nodes to form a preliminary set of irregular plate morphological deformation data.

[0078] The deformation data sequences of individual nodes in all candidate topological node regions are aggregated and organized according to the order of node data identifiers to form a preliminary irregular plate morphological deformation dataset. This dataset adopts a two-dimensional array structure, where rows represent different node data identifiers, columns represent time series, and array elements are the deformation data items (timestamp, real-time deformation, deformation rate) of the corresponding node at the corresponding time. The preliminary dataset contains deformation information of all key contact areas of the irregular plate, but does not yet include spatial location and other related information.

[0079] Step S1249: Supplement the collection timestamp and node location information of each data sequence, improve the association information of the irregular plate morphological deformation data set, and generate a complete irregular plate morphological deformation data set. The irregular plate morphological deformation data set arranges the deformation information of each topological node candidate region in the order of collection time.

[0080] Based on the initial dataset, detailed acquisition timestamps (accurate to milliseconds) and node spatial location information (3D coordinates) are added to each data sequence. The timestamps are used for synchronization with the mold reset force data, and the node location information is used for subsequent spatial correlation analysis of topological nodes. The refined irregular plate morphology deformation dataset is stored in a structured data format, with each data record containing five fields: node identifier, spatial coordinates, timestamp, real-time deformation amount, and deformation rate. The dataset is arranged in ascending order of acquisition time, and the deformation process data for a specific area can be quickly indexed using the node identifier.

[0081] Step S125: Collect real-time force data, force change rate data, and force transmission duration data of each topology node matching area through pressure sensors to form a real-time mold reset force data set. The real-time mold reset force data set is arranged in the order of collection time, with the force information of each topology node matching area arranged in order of collection time.

[0082] During the reset process, the pressure sensor collects force data from the matching area of ​​the topology node. Real-time force data is a direct measurement from the sensor, representing the magnitude of the force applied by the mold to the irregular plate. Force change rate data is obtained by calculating the time derivative of the real-time force data, reflecting the rate of change of force over time. Force transmission duration data is obtained by recording the time interval from the start to the end of the force application. The above data are arranged in chronological order of acquisition, with each topology node matching area corresponding to a set of time-series data. The data from all matching areas are integrated to form a real-time mold reset force dataset. This real-time mold reset force dataset records the changes in force at key locations of the mold during the reset process.

[0083] Step S126: Extract deformation data of all candidate regions of topological nodes at the same acquisition time point from the deformation data set of the irregular plate shape, and generate instantaneous deformation data frames. Each instantaneous deformation data frame contains the complete deformation parameters of all candidate regions of topological nodes at the acquisition time point.

[0084] From the deformation data set of the irregularly shaped plate, deformation data of all candidate regions of topological nodes at the same acquisition time point are extracted according to the timestamp. This data includes parameters such as the real-time deformation and deformation rate of each candidate region. These parameters are organized into a data structure according to the candidate region's numbering order, generating instantaneous deformation data frames. Each instantaneous deformation data frame corresponds to a specific time point and contains the deformation state information of all key regions of the irregularly shaped plate at that moment. This method converts the time-series deformation data into a data structure based on time frames, facilitating subsequent synchronization with the reset force data.

[0085] Step S127: Extract the force data of all topology node matching areas at the same acquisition time point from the real-time data set of mold reset force, and generate instantaneous force data frames. Each instantaneous force data frame contains the complete force parameters of all topology node matching areas at that time point.

[0086] Similarly, force data for all matching regions of topological nodes at the same acquisition time point are extracted from the real-time data set of mold reset force. This data includes parameters such as real-time force, force change rate, and force transmission duration. These parameters are organized into a data structure according to the matching region numbering order to generate instantaneous force data frames. Each instantaneous force data frame corresponds to a time point and contains the force status information of all key locations of the mold at that moment. The instantaneous force data frames and instantaneous deformation data frames have the same timestamp.

[0087] Step S128: Associate and bind the instantaneous deformation data frame with the corresponding instantaneous force data frame, and connect all the associated and bound frame data in sequence according to the acquisition time to form a preliminary morphological force coupling data chain.

[0088] Instantaneous deformation data frames and instantaneous force data frames are associated and bound one-to-one according to the timestamp. Each bound data unit contains deformation data and force data at the same time point. The bound data units are arranged sequentially according to the acquisition time to form a preliminary shape-force coupled data chain. Each element in the data chain is a comprehensive data unit containing timestamp, deformation data, and force data, which can reflect the deformation state of the irregular plate and the corresponding mold force state at a certain moment during the resetting process.

[0089] Step S129: Supplement the timestamp information and node corresponding identifier of each frame of data in the preliminary morphological force coupling data chain, mark the node to which each deformation data and force data belongs and the acquisition time, and generate a complete morphological force coupling data chain. The morphological force coupling data chain continuously records the synchronous change information of deformation and corresponding reset force of each contact area of ​​the irregular plate in chronological order. Each data node includes spatial position, deformation parameters, force parameters and time parameters.

[0090] Based on the initial morphoforce coupling data chain, detailed timestamp information, accurate to the millisecond level, is added to each frame of data to ensure time accuracy. Simultaneously, node-specific identifiers are added to each deformation and force data point, clearly identifying the candidate or matching region of the topological node to which the data belongs. By adding this information, each data point in the data chain can be traced back to a specific node and acquisition time. After supplementation and annotation, a complete morphoforce coupling data chain is generated. This morphoforce coupling data chain continuously records the synchronous changes in deformation and corresponding reset force in each contact area of ​​the irregularly shaped plate in chronological order. Each data node includes spatial location (node ​​coordinates), deformation parameters (deformation amount, deformation rate, etc.), force parameters (force magnitude, rate of change, etc.), and time parameters (timestamp).

[0091] Step S130: Input the morphological force coupling data chain into the irregular plate reset dynamic topology model, drive the irregular plate reset dynamic topology model to perform dynamic reorganization of topology nodes, and generate topology reorganization data. The topology reorganization data reflects the node connection relationship after the deformation of the irregular plate and the reset force are adapted.

[0092] The complete morpho-force coupling data chain is imported into the dynamic topology model for resetting the irregularly shaped plate. The model updates the state of the topology nodes time-by-time based on the time-series data in the data chain. For each time point, the model analyzes the deformation and force data to determine whether the relationships between the topology nodes need adjustment. When the deformation or force exceeds a set threshold, the model dynamically reassembles the corresponding topology nodes, including adjusting the node association weights and reconstructing the force transmission path. During the reassembly process, the model comprehensively considers the material properties of the irregularly shaped plate, the motion characteristics of the mold, and the interaction between the two to ensure that the reassembled node connections accurately reflect the actual contact and force transmission. Through processing the entire morpho-force coupling data chain, the model completes the dynamic reassembly of the topology nodes, generating topology reassembly data. This topology reassembly data records the changes in the connection relationships between the topology nodes at different times during the resetting process, reflecting the dynamic node associations after the irregularly shaped plate deformation and resetting force are adapted.

[0093] Step S131: Analyze each instantaneous deformation data frame and instantaneous force data frame in the morphological force coupling data chain in chronological order, and extract the actual deformation of each candidate region of the topology node and the actual force of the corresponding matching region of the topology node.

[0094] The morphological-force coupled data chain is analyzed in time sequence, processing each instantaneous deformation data frame and its corresponding instantaneous force data frame. From the instantaneous deformation data frames, the actual deformation variables of each candidate region of the topological node are extracted. These deformation variables are actual data measured by sensors, reflecting the deformation state of the irregularly shaped plate at that moment. From the instantaneous force data frames, the actual forces of the matching region of each topological node are extracted. These forces are the actual force values ​​applied to the irregularly shaped plate by the mold. The extracted actual deformation variables and actual forces are associated by node number to ensure that the deformation variable of each candidate region corresponds to the force of its corresponding matching region. Through the above analysis process, the comprehensive data in the data chain is decomposed into the specific deformation and force data of each node.

[0095] Step S132: Calculate the deformation force adaptation coefficient of each basic topology node. The deformation force adaptation coefficient is the ratio of the actual deformation of the candidate region of the corresponding topology node to the actual force of the matching region of the topology node, reflecting the degree of adaptation between deformation and force at the node level.

[0096] For each basic topology node, the deformation-force adaptation coefficient is obtained by dividing the actual deformation of its corresponding candidate region by the actual force applied to its matching region. This coefficient reflects the deformation per unit force at the node level and is an indicator of the degree of adaptation between deformation and force. If the coefficient is within a reasonable range, it indicates that the relationship between deformation and force at that node meets expectations; if the coefficient is too large or too small, it suggests that the node may be abnormal and requires further analysis and adjustment. By calculating the deformation-force adaptation coefficients of all basic topology nodes, the adaptation status of each node can be comprehensively understood.

[0097] Step S133: Input the deformation force adaptation coefficients of all basic topology nodes into the irregular plate reset dynamic topology model, compare them with the adaptation coefficient threshold initially set in the irregular plate reset dynamic topology model, and filter out abnormal nodes whose adaptation coefficients exceed the threshold range.

[0098] The calculated deformation force adaptation coefficients of all basic topological nodes are input into the dynamic topology model for irregular plate resetting. A threshold range for the adaptation coefficients is pre-defined in the model, determined based on the material properties of the irregular plate, the mold design parameters, and historical experimental data. The adaptation coefficient of each node is compared with the threshold range; if the coefficient exceeds the range, the node is identified as an abnormal node. Abnormal nodes may occur due to poor contact, sensor malfunction, or unreasonable model parameter settings. By filtering abnormal nodes, the model can be adjusted and optimized in a targeted manner.

[0099] Step S134: For abnormal nodes, extract continuous multi-frame deformation and force change data corresponding to the abnormal nodes in the morphological force coupling data chain, analyze the deformation trend and force change trend of the abnormal nodes, and locate the core cause of the node abnormality.

[0100] For the identified anomalous nodes, deformation and force variation data for multiple consecutive time frames corresponding to the node are extracted from the morpho-force coupling data chain. Trend analysis is performed on these data to observe the changes in deformation over time (deformation trend) and the changes in force over time (force variation trend). The consistency between the deformation trend and the force variation trend is compared. If the trends do not match (e.g., force increases while deformation decreases), there may be an abnormal contact state; if the trends are consistent but the values ​​exceed the normal range, it may be a problem with the sensor or model parameters. Through in-depth analysis of these data, the core cause of the node anomaly is identified.

[0101] Step S135: Based on the deformation trend and force change trend of abnormal nodes, adjust the association weight between abnormal nodes and adjacent nodes, and adjust the node association strength according to the degree of fit between deformation and force. When the degree of fit between deformation and force meets the preset standard, increase the association strength; when it does not meet the preset standard, decrease the association strength.

[0102] Based on the analysis results of the deformation and force change trends of abnormal nodes, the association weights between these nodes and their adjacent nodes are adjusted. If the compatibility between deformation and force meets the preset standards (e.g., both trends are consistent and values ​​are within a reasonable range), the association weights between these nodes and their adjacent nodes are appropriately increased to strengthen the association. If the compatibility does not meet the preset standards, the association weights are decreased to weaken the association. During the adjustment process, the magnitude of the change in association weights is determined by the degree of deviation from the compatibility; the greater the deviation, the greater the adjustment. This method ensures that the association relationships between nodes better reflect the actual deformation and force transmission.

[0103] Step S136: Reconstruct the force transmission path corresponding to the abnormal node, adjust the path transmission direction according to the deformation trend, so that the force transmission path corresponds to the actual deformation direction of the irregular plate, and optimize the path transmission coefficient to match the actual force transmission efficiency.

[0104] For base topology nodes identified as anomalous, their deformation trends during the reset process are first analyzed to determine the main direction and characteristics of the deformation. Based on the main deformation direction, the force transmission path between this node and its adjacent nodes is adjusted to ensure that the path direction aligns with the deformation direction, avoiding force transmission efficiency loss due to conflict between the two. Simultaneously, based on the relationship between the actual force and deformation recorded in the form-force coupling data chain, the actual force transmission efficiency of this path is calculated. The path transmission coefficient is adjusted to ensure that the model-calculated transmission efficiency matches the actual value.

[0105] Step S1361: Extract deformation trend data of abnormal nodes and determine the main direction and magnitude of deformation.

[0106] Real-time deformation data sequences of abnormal nodes throughout the entire reset process are extracted from the morpho-force coupled data chain. The deformation direction vectors at different times are calculated using the sliding window method. Principal component analysis of the direction vectors determines the main deformation directions (e.g., positive X-axis, negative Y-axis, etc.). The magnitude of change data is obtained by calculating the difference between the maximum and minimum values ​​of the deformation data sequence, reflecting the total range of deformation. Deformation trend data also includes the time point when the deformation reaches its peak and the deformation value during the stable phase, used to analyze the dynamic process of deformation.

[0107] Step S1362: Based on the main deformation direction data, plan a new force transmission path direction so that the new path direction forms a complementary angle with the main deformation direction.

[0108] Based on the dominant deformation direction of the abnormal node, a new force transmission path direction is designed. When the dominant deformation direction is tension, the path direction should form a preset complementary angle (e.g., 45 degrees) with the tension direction to effectively transmit the force resisting tension; when the dominant deformation direction is compression, the path direction forms another preset complementary angle with the compression direction. The planning of the path direction needs to refer to the material elastic parameters of the irregular plate to ensure that the path direction conforms to the force transmission law in material mechanics and avoids stress concentration.

[0109] Step S1363: Based on the deformation change amplitude data, set the number of branches in the force transmission path. The deformation change amplitude data and the number of path branches are allocated according to a preset correspondence, and the number of branches is determined according to the actual deformation requirements.

[0110] A pre-defined table maps deformation amplitudes to the number of path branches. For example, a deformation amplitude within the range of A1 corresponds to one branch, within the range of A2 corresponds to two branches, and so on. Based on the deformation amplitude data of abnormal nodes, the required number of path branches is determined by looking up the table. The larger the deformation amplitude, the more branches are needed to distribute the transmitted force and avoid overloading a single path. The starting point of each branch path is an abnormal node, and the ending point is a different adjacent normal node, forming a radial transmission network.

[0111] Step S1364: Extract the position data of adjacent normal nodes of the abnormal node, and use the adjacent normal nodes as the connection nodes of the new force transmission path.

[0112] By using the topological node association matrix, neighboring nodes with initial associations to the anomalous node are identified, and normal nodes with fit coefficients within the normal range are selected. The spatial coordinate data of these normal nodes are extracted and used as the endpoint connection nodes for the new force transmission path. The selection of connection nodes must meet spatial distance requirements, prioritizing normal nodes closest to the anomalous node to shorten the force transmission path length and improve transmission efficiency.

[0113] Step S1365: Construct path segments between abnormal nodes and adjacent normal nodes to form a preliminary reconstruction force transmission path.

[0114] In a three-dimensional coordinate system, starting from an abnormal node and ending at each adjacent normal node, straight line segments are drawn as the initial force transmission path. The attributes of the path segments include the starting coordinates, ending coordinates, path length (calculated using the distance formula between two points), and path direction vector (obtained by subtracting the starting coordinates from the ending coordinates). The initially reconstructed force transmission path does not consider the intersection interference between paths and is constructed solely based on spatial positional relationships.

[0115] Step S1366: Calculate the length data and path resistance data of the initial reconfiguration force transmission path, and set the initial path transmission coefficient based on the length data and path resistance data of the initial reconfiguration force transmission path.

[0116] Path length data is calculated using the distance formula between two points in space. Path resistance data is calculated based on the thickness distribution data and material elasticity data of the irregularly shaped plates in the area traversed by the path; the greater the thickness and the higher the elastic modulus, the greater the path resistance. The initial path transfer coefficient is set as a function inversely proportional to both path length and path resistance; that is, the shorter the path and the lower the resistance, the larger the initial transfer coefficient, and vice versa. The initial transfer coefficient is limited to a value between 0 and 1.

[0117] Step S1367: Extract the actual force transmission efficiency data corresponding to the abnormal nodes in the morphological force coupling data chain.

[0118] Force and deformation data of anomalous nodes and their corresponding matching regions are extracted from the force-morphocoupled data chain. The actual force transmission efficiency is calculated using the formula: "Actual transmission efficiency = (Anomalous node deformation × Adjacent node deformation) / (Anomalous node force × Adjacent node force)," reflecting the actual force transmission effect between the anomalous node and its adjacent nodes. During the calculation, it is crucial to ensure that the timestamps of the force and deformation are synchronized.

[0119] Step S1368: Adjust the initial path transmission coefficient so that the transmission efficiency corresponding to the adjusted path transmission coefficient corresponds to the actual force transmission efficiency data.

[0120] Substitute the initial path transfer coefficient into the dynamic topology model to calculate the theoretical force transmission efficiency. Compare the theoretical value with the actual force transmission efficiency data. If the theoretical value is lower than the actual value, increase the transfer coefficient; if the theoretical value is higher than the actual value, decrease the transfer coefficient. The adjustment range is determined based on the deviation between the two values; the larger the deviation, the larger the adjustment range, until the deviation between the theoretical transmission efficiency and the actual data is within the set error range.

[0121] Step S1369: Record the direction data, branch data, connection node data and optimized path transmission coefficient of the reconstructed force transmission path, and generate complete abnormal node force transmission path reconstruction data. The abnormal node force transmission path reconstruction data includes the adjusted path transmission direction and the path transmission coefficient that matches the actual force transmission efficiency.

[0122] The reconstructed force transmission path's direction vector, number of branches, coordinates of each branch's connecting nodes, and optimized path transmission coefficients are organized into structured data. Each path branch corresponds to one record, containing six fields: path ID, starting node ID, ending node ID, direction vector (X, Y, Z), path length, and transmission coefficient. All branch data are aggregated to form the abnormal node force transmission path reconstruction data, which is stored as a component of the topology reorganization data.

[0123] Step S137: For normal nodes that do not exceed the adaptation coefficient threshold, fine-tune the associated weights and path transmission coefficients of the nodes according to the actual data in the morphological force coupling data chain to make the node parameters fit the actual working conditions.

[0124] For normal nodes with adaptation coefficients within the threshold range, parameter fine-tuning is also necessary based on the actual data in the morpho-force coupling data chain. The actual deformation and force data of the nodes are compared with the simulation results of the model, and the deviation between the two is calculated. Based on the magnitude and direction of the deviation, the association weights and path transmission coefficients of the nodes are slightly adjusted to better reflect actual working conditions. This fine-tuning improves the overall accuracy of the model, ensuring that all nodes accurately reflect the actual contact and force transmission.

[0125] Step S138: Record the association weight adjustment data, force transmission path reconstruction data and parameter fine-tuning data of all nodes to form a node adjustment detail. Integrate the node adjustment details to generate the node association matrix and force transmission path model after dynamic reorganization of topology nodes, and mark the connection relationship and transmission rules of the reorganized basic topology nodes.

[0126] During the adjustment and optimization of abnormal and normal nodes, detailed records are kept of the association weight adjustment data, force transmission path reconstruction data, and parameter fine-tuning data for each node, forming a node adjustment detail. This detail includes information such as node number, parameter values ​​before adjustment, parameter values ​​after adjustment, and the reason for adjustment. The node adjustment details are then integrated to update the initial topology association matrix and force transmission path model, generating a new matrix and model after dynamic reorganization of the topology nodes. The new model annotates the connection relationships and transmission rules between the reorganized basic topology nodes, clarifying the associated objects, force transmission paths, and coefficients for each node.

[0127] Step S139: Integrate the recombined node association matrix, force transmission path model and all node adjustment parameters to generate topology recombined data. The topology recombined data is used to reflect the node connection relationship after the deformation of the irregular plate and the reset force are adapted.

[0128] The recombined node association matrix, force transmission path model, and all node adjustment parameters (such as association weight adjustment values ​​and path transmission coefficient optimization values) are integrated and organized according to a certain format to form topology recombined data. The topology recombined data comprehensively reflects the connection relationships, force transmission paths, and related parameters between each basic topology node after the deformation and reset force of the irregular plate are adapted.

[0129] Step S140: Based on the topology reorganization data, generate dynamic adaptation reset path data, which includes the motion trajectory and action sequence of each component of the mold after adjustment according to the topology node relationship.

[0130] By utilizing the node connectivity and force transmission path information in the topology reconfiguration data, motion trajectories and action sequences are planned for each reset component of the mold. The target position of the component is determined based on the spatial coordinate changes of the nodes, and the motion priority and force magnitude of the component are determined based on the force transmission path and associated weights. Considering the time requirements of the reset process, the action sequence and time intervals of each component are rationally arranged to ensure that the mold's reset motion adapts to the deformation of the irregular plate. By comprehensively considering these factors, dynamically adapted reset path data is generated, which details the motion trajectory, velocity, acceleration, and sequence of actions of each component of the mold during the reset process.

[0131] Step S141: Analyze the node association matrix in the topology reorganization data, extract the spatial coordinate change data of the candidate area of ​​the irregular plate topology node corresponding to the matching area of ​​each mold topology node, and obtain the target position adjustment information of the action area of ​​each component.

[0132] The node association matrix in the topology reorganization data is analyzed to determine the correspondence between the matching area of ​​each mold topology node and the candidate area of ​​the irregular plate topology node. These correspondences are tracked over time, and the spatial coordinate change data of the candidate areas corresponding to the matching areas is extracted. Based on the coordinate change data, the target position to which the functional area (i.e., the matching area) of each mold component needs to be adjusted is calculated. The target position adjustment information includes parameters such as the change in position coordinates and the adjustment direction.

[0133] Step S142: Extract the force transmission path model from the topology reorganization data, determine the force transmission priority of each mold component, and sort the components according to the path transmission coefficient values ​​to form the force execution order.

[0134] Force transmission path models are extracted from the topology reassembly data, and the transmission coefficients of each path in the model are analyzed. A larger transmission coefficient indicates higher force transmission efficiency for that path, and a higher priority for the force transmission of the corresponding mold component during the reset process. All mold components are sorted from largest to smallest according to their corresponding path transmission coefficient values ​​to form the component force execution order. During the reset process, components with higher priority apply force first, or apply a larger force, to ensure that force transmission follows the optimal path, improving reset efficiency and accuracy.

[0135] Step S143: Combine the target position adjustment information and the force execution sequence to plan the initial motion trajectory for each mold reset component. The initial motion trajectory is set according to the shortest path from the current position of the component to the target position, while incorporating the force transmission priority constraint.

[0136] Based on the target position adjustment information and force execution sequence of each mold reset component, the initial motion trajectory is planned. The initial motion trajectory planning is based on the shortest path from the component's current position to the target position, typically using a straight line or circular arc path. Simultaneously, force transmission priority constraints are considered; for high-priority components, their motion trajectories can be appropriately adjusted to avoid interference with other components or to ensure they can reach the target position and function preferentially. During the planning process, kinematic analysis is required to ensure the feasibility and smoothness of the trajectory.

[0137] Step S144: Extract the timestamp information from the morphological force coupling data chain, determine the time node for adjusting the association relationship of each topology node, and use the time node for adjusting the association relationship of each topology node as the segment switching point of the mold motion trajectory.

[0138] The timestamps of all topology node relationship adjustments are extracted from the morphological force coupling data chain. These timestamps correspond to the moments when node relationships change. These timestamps are then used as segment switching points for the mold's motion trajectory. At these switching points, the motion trajectory of the mold components may need to be adjusted according to the new node relationships. By setting these switching points, the entire reset process can be divided into multiple time periods, each corresponding to a stable set of node relationships and motion trajectories.

[0139] Step S145: Divide the initial motion trajectory into multiple continuous trajectory segments according to the segment switching points. Each trajectory segment corresponds to the component motion path within a time interval, so that the relationship between the trajectory segments and the nodes is adjusted and synchronized.

[0140] Based on the determined segmentation switching points, the initial motion trajectory is divided into multiple continuous trajectory segments. Each trajectory segment corresponds to a time interval, within which the node relationships remain stable, and the mold component moves according to that trajectory segment. The division of trajectory segments ensures that the adjustment of the motion trajectory and node relationships is synchronized; when the node relationships change, the motion trajectory switches to the new segment accordingly. Each trajectory segment includes a start point, an end point, a motion path, and corresponding time parameters.

[0141] Step S146: Assign motion speed parameters to each trajectory segment. The motion speed parameters are set according to the force transmission coefficient values ​​within the time interval. The force transmission coefficient values ​​and motion speed values ​​are established in a corresponding relationship according to a preset ratio. When the force transmission coefficient values ​​are in a preset high range, the motion speed values ​​are set to the corresponding range values ​​according to the preset ratio. When the force transmission coefficient values ​​are in a preset low range, the motion speed values ​​are set to the corresponding range values ​​according to the preset ratio.

[0142] Based on the force transmission coefficient value within the time interval corresponding to each trajectory segment, motion speed parameters are assigned to the trajectory segment. A preset proportional relationship between the force transmission coefficient and motion speed is established. When the force transmission coefficient value is in a preset high range, it indicates high force transmission efficiency within that time period, allowing the component to move at a higher speed; the motion speed value is set proportionally to the corresponding high range value. Conversely, when the force transmission coefficient value is in a preset low range, the component's motion speed is correspondingly reduced, set proportionally to the corresponding low range value. This method matches the component's motion speed with the force transmission efficiency, improving the stability and efficiency of the reset process.

[0143] Step S147: Integrate all trajectory segments with corresponding velocity parameters and time intervals to form dynamic motion trajectory data for each mold reset component. The dynamic motion trajectory data for each mold reset component clearly defines the motion path and velocity of the component in different time intervals.

[0144] The dynamic motion trajectory data of each mold reset component is formed by integrating all trajectory segments, corresponding velocity parameters, and time intervals. This data is organized chronologically and details the component's motion path (including start point, end point, and path shape) and velocity within different time intervals. This data allows us to understand the component's motion throughout the entire reset process.

[0145] Step S148: Based on the dynamic motion trajectory data of each mold reset component and the force transmission priority, coordinate the motion timing of each mold component; the coordination of the motion timing of each mold component includes: generating mutually exclusive motion periods of the components in the overlapping area based on the spatial position overlap area data in the dynamic motion trajectory data of adjacent components; generating serialized start instructions for the components based on the force transmission priority.

[0146] After generating the dynamic motion trajectory data of each mold reset component, it is necessary to further coordinate their motion timing to avoid motion interference and ensure the orderly transmission of reset force. First, analyze the motion trajectory of all components to identify areas with overlapping spatial positions and plan mutually exclusive motion periods for these areas; then, determine the starting order of each component according to the force transmission priority, generate serialized start commands, and ensure that high-priority components act first.

[0147] For example, step S1481: Extract the key motion stage data from the dynamic motion trajectory data of each mold reset component, and extract the motion start time data and end time data of each component.

[0148] From the dynamic motion trajectory data of each mold reset component, key motion stages such as acceleration, constant speed, and deceleration are separated. Each stage includes information such as start time, end time, and motion path within the stage. The start time data (start time of the first stage) and end time data (end time of the last stage) of the entire motion process are extracted to determine the total duration and time window of the component's motion.

[0149] Step S1482: Based on the force transmission priority, generate a motion start-up sequence list for each component, and assign a start-up time to each component based on the start-up sequence list.

[0150] Sort the mold reset components according to their force transmission priority from high to low, generating a motion start sequence list. The component with the highest priority is assigned the earliest start time, and the start time of subsequent components is determined based on the start time of the preceding component and a preset time interval. The time interval setting must ensure that subsequent components start only after the preceding component has started moving and established a stable force, avoiding impacts caused by force superposition.

[0151] Step S1483: Calculate the spatial location overlap area data in the dynamic motion trajectory data of adjacent components, and determine the area where the trajectory overlaps and the corresponding time interval.

[0152] Spatial intersection calculations are performed on the dynamic motion trajectory data of all adjacent mold reset components to identify three-dimensional spatial regions where trajectory segments overlap. For each overlapping region, the time interval (entry and exit times) of each component's motion trajectory through that region is extracted. The spatial intersection calculation uses a bounding box collision detection algorithm, firstly screening potentially overlapping trajectory segments through coarse detection, and then performing precise line segment intersection calculations.

[0153] Step S1484: For regions where trajectories overlap, adjust the motion speed parameters or start time of relevant components in that region based on the time staggering principle to generate mutually exclusive passage times.

[0154] When the motion trajectories of two components overlap in a certain area, adjustments are made using a time staggering principle. If component A has a higher priority than component B, the motion speed or start time of component B is adjusted so that component B enters the overlapping area after component A has left. The preferred adjustment method is to change the speed parameter (e.g., reduce the speed of component B before entering the overlapping area). If speed adjustment cannot meet the requirements, the start time of component B is delayed. After adjustment, the motion trajectory of component B must be recalculated to ensure that its total motion time remains within the allowable range.

[0155] Step S1485: Based on the motion start sequence list and the motion start time data and end time data of each component, generate the motion connection logic between components, wherein the motion connection logic includes an instruction to trigger the motion start of the next component based on the completion status of the key motion of the previous component.

[0156] Based on the motion initiation sequence list, trigger conditions for action connection are set for adjacent components. For example, when the preceding component completes the acceleration phase and enters the constant speed phase, the start command for the following component is triggered; or when the preceding component reaches a certain key position point, the following component is allowed to start moving. The connection logic is represented in the form of a state machine, containing three elements: trigger events (such as position reached, speed reached), condition judgments, and action commands, ensuring the orderly connection of component actions.

[0157] Step S1486: Based on the adjusted motion speed parameters, start time, and action connection logic of each component, recalculate the motion stroke data of each component.

[0158] Substitute the adjusted motion velocity parameters (such as velocity values ​​at each stage) and start time into the kinematic equations to recalculate the position coordinates of each component at each time point, generating adjusted motion travel data. The travel data must satisfy the constraints of the motion connection logic, meaning the starting position of the subsequent component must maintain a safe distance from the current position of the preceding component. During the recalculation process, it is necessary to check for any new trajectory overlaps; if so, adjustments must be made again.

[0159] Step S1487: Set motion timing coordination verification points, which are distributed in time sequence throughout the reset process. Each verification point corresponds to a set of component motion status data. Verify the motion status data of each component in the order of motion timing coordination verification points.

[0160] Multiple verification points are evenly distributed along the timeline of the entire reset process, with each verification point corresponding to a specific time value. At each verification point, the current position, velocity, acceleration, and other motion state data of all mold reset components are extracted to check for issues such as overlapping positions, abnormal speeds, or failure to follow the connection logic. The number of verification points is determined based on the total reset duration to ensure coverage of critical motion phases (such as the start-up phase, the overlapping area passage phase, and the stop phase).

[0161] Step S1488: Record the start-up time allocation data, speed parameter adjustment data, action connection logic and verification result data to form motion timing coordination data; integrate the motion timing coordination data with the dynamic motion trajectory data of each mold reset component.

[0162] The start-up time allocation table, speed parameter adjustment records, action connection logic state machine, and verification point result data generated during the motion timing coordination process are summarized to form motion timing coordination data. This motion timing coordination data is integrated with the dynamic motion trajectory data of each component by timestamp, so that the trajectory data includes timing constraint information. The integrated dynamic adaptation reset path data includes both the motion trajectory parameters of the components and the timing parameters that ensure motion coordination.

[0163] Step S149: Extract key position points, motion speed change points, and action switching points from the dynamic motion trajectory data of each mold reset component to form a set of core parameters for dynamically adapting the reset path.

[0164] Key parameters are extracted from the dynamic motion trajectory data of each mold reset component, including key position points (such as start point, end point, and turning point), motion speed change points (such as acceleration points and deceleration points), and action switching points (such as trajectory segment switching points). These parameters together constitute the core parameter set for dynamically adapting the reset path. The core parameter set contains key information describing the reset path and is used for subsequent path optimization and control command generation.

[0165] Step S1410: Integrate the dynamic motion trajectory data, motion timing coordination data and core parameter set of each mold component to generate dynamic adaptation reset path data. The dynamic adaptation reset path data includes the motion trajectory and action timing of each mold component after adjustment according to the topological node relationship.

[0166] The dynamic motion trajectory data, motion timing coordination data (such as mutually exclusive motion periods and serialized start commands), and core parameter sets of each mold component are integrated and organized in a unified format to form dynamic adaptation reset path data. This dynamic adaptation reset path data comprehensively describes the motion trajectory, action timing, speed changes, and other information of each mold component during the reset process, and serves as the direct basis for controlling the mold reset motion.

[0167] Step S150: Input the dynamic adaptation reset path data in reverse into the dynamic topology model of the irregular plate reset, iteratively optimize the topology node association strength and reset force transmission coefficient, and generate the iteratively optimized simulation complete data.

[0168] The generated dynamically adapted reset path data is input back into the dynamic topology model for irregular plate reset. The model simulates the motion process of the mold components based on the path data and calculates the corresponding forces and irregular plate deformation. The simulation results are compared with the actual data in the form-force coupling data chain, and the error between the two is calculated. The association strength of the topology nodes and the reset force transmission coefficient are adjusted according to the error magnitude, and then the simulation and comparison are repeated. Through multiple iterative optimizations, the simulation results of the model are made as close as possible to the actual data, and finally, the iteratively optimized complete simulation data is generated. This complete simulation data includes the optimized topology node parameters, reset force transmission parameters, and mold motion path data, which can accurately reflect the real situation of the irregular plate stamping die reset process.

[0169] Step S151: Extract the motion trajectory coordinate data and action timing data of each component of the mold from the dynamic adaptation reset path data, and convert the motion trajectory coordinate data and action timing data of each component of the mold into node input parameters that can be recognized by the dynamic topology model of the irregular plate reset.

[0170] The motion trajectory coordinate data (such as the position coordinates at each time point) and action timing data (such as start time and motion duration) of each component of the mold are extracted from the dynamic adaptation reset path data. This data is then converted into a node input parameter format that the irregular plate reset dynamic topology model can recognize and process. During the conversion process, coordinate system unification and data format conversion are required to ensure that the input parameters match the model's requirements. Node input parameters include the component's position coordinates, motion time, velocity, and other information, used to drive the model for simulation calculations.

[0171] Step S152: Input the converted node input parameters in reverse into the irregular plate reset dynamic topology model, replacing the original node motion parameters and timing parameters in the irregular plate reset dynamic topology model.

[0172] The transformed node input parameters are input into the dynamic topology model for irregular plate reset, replacing the original node motion parameters (such as position and velocity) and timing parameters (such as start-up time and motion cycle) in the model. By updating these parameters, the model can simulate based on the new dynamically adapted reset path. During the replacement process, the correctness and consistency of the parameters are ensured to avoid deviations in simulation results due to parameter errors.

[0173] Step S153: Extract the initial value of node association strength from the topology reorganization data and compare it with the actual association strength data corresponding to the component motion trajectory in the dynamic adaptation reset path data.

[0174] Initial values ​​for node association strength are extracted from the topology reorganization data; these initial values ​​were determined during the previous dynamic reorganization process. Simultaneously, actual association strength data corresponding to the component's motion trajectory is extracted from the dynamic adaptation reset path data. This actual association strength data is obtained by analyzing the interactions between nodes during component movement. The initial values ​​and actual association strength data are compared, and the differences between them are calculated.

[0175] Step S154: Based on the comparison results, adjust the topology node association strength parameters in the dynamic topology model of the irregular plate reset so that the topology node association strength parameters correspond to the actual association strength data.

[0176] Based on the comparison between the initial node association strength values ​​and the actual data, the topological node association strength parameters in the model are adjusted. If the initial value is lower than the actual association strength data, the association strength parameter is increased; if the initial value is higher than the actual data, the association strength parameter is decreased. The adjustment range is determined by the magnitude of the difference; the greater the difference, the greater the adjustment range. Through adjustment, the association strength parameters in the model can accurately reflect the actual node association situation.

[0177] Step S155: Extract the initial value of the reset force transmission coefficient from the topology reorganization data and compare it with the actual transmission coefficient data corresponding to the force transmission efficiency in the dynamic adaptation reset path data.

[0178] Initial values ​​for the reset force transmission coefficient are extracted from the topology reconfiguration data, and actual transmission coefficient data corresponding to the force transmission efficiency are extracted from the dynamic adaptation reset path data. The actual transmission coefficient data is calculated by measuring the transmission efficiency of the actual force between nodes. The initial values ​​are compared with the actual transmission coefficient data to analyze the deviation between them.

[0179] Step S156: Based on the comparison results, optimize the reset force transmission coefficient in the dynamic topology model of the irregular plate reset so that the reset force transmission coefficient corresponds to the actual transmission coefficient data.

[0180] Based on the comparison between the initial value of the reset force transmission coefficient and the actual data, the reset force transmission coefficient in the model is optimized. If the initial transmission coefficient causes the simulated force transmission efficiency to be lower than the actual transmission coefficient data, the transmission coefficient is increased; otherwise, the transmission coefficient is decreased. Through optimization, the model can accurately simulate the force transmission efficiency between nodes, improving the accuracy of the simulation.

[0181] Step S157: Record the adjustment data of the topology node association strength parameter and the optimization data of the reset force transmission coefficient to form an iterative optimization detail. Input the iterative optimization detail into the irregular plate reset dynamic topology model to drive the irregular plate reset dynamic topology model to recalculate the node association relationship and force transmission path.

[0182] Detailed records are kept of the adjustment data for the topological node association strength parameters (such as the values ​​before and after adjustment, and the adjustment range) and the optimization data for the reset force transmission coefficient, forming an iterative optimization detail. This detail is input into the dynamic topology model for irregular plate reset, and the model recalculates the association relationships and force transmission paths between nodes based on this data. During the recalculation process, the model comprehensively considers the adjusted association strength and transmission coefficient, updating the connection relationships and path parameters of the nodes.

[0183] Step S158: Extract the recalculated node association data, force transmission path data, and component motion trajectory optimization data, and output the iteratively optimized complete simulation data. The iteratively optimized complete simulation data includes the adjusted topology node parameters, reset force transmission parameters, and mold motion path data.

[0184] Node association data (such as the association matrix), force transmission path data (such as the path transmission coefficient), and component motion trajectory optimization data (such as the optimized trajectory coordinates and velocity parameters) are extracted from the recalculated model. These data are then integrated to form the complete simulation data after iterative optimization. This complete simulation data includes topological node parameters, reset force transmission parameters, and mold motion path data optimized through multiple iterations, accurately reflecting the dynamic characteristics of the resetting process of the irregular plate stamping mold.

[0185] In one exemplary embodiment, a simulation system for the resetting process of a shaped plate stamping die, incorporating data modeling, is provided. This simulation system can be a terminal, server, etc., and its internal structure diagram can be as follows: Figure 2As shown, the simulation system for the resetting process of a non-standard sheet metal stamping die, which incorporates data modeling, includes a processor, memory, input / output interface, communication interface, display unit, and input device. The processor, memory, and input / output interface are connected via a system bus, and the communication interface, display unit, and input device are also connected to the system bus via the input / output interface. The processor provides computational and control capabilities. The memory includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores the operating system and computer programs. The internal memory provides the environment for the operation of the operating system and computer programs in the non-volatile storage medium. The input / output interface is used for exchanging information between the processor and external devices. The communication interface is used for wired or wireless communication with external terminals; wireless communication can be achieved through Wi-Fi, mobile cellular networks, near-field communication, or other technologies. When the computer program is executed by the processor, it implements a simulation method for the resetting process of a non-standard sheet metal stamping die, incorporating data modeling. The display unit is used to generate a visually visible image and can be a display screen, projection device, or virtual reality imaging device. The display screen can be an LCD screen or an e-ink screen. The input device can be a touch layer covering the display screen, or a button, trackball, or touchpad set on the shell of the simulation system for resetting the irregular plate stamping die combined with data modeling, or an external keyboard, touchpad, or mouse, etc.

[0186] It should be noted that, in order to simplify the description of the present invention and thus help to understand one or more embodiments of the invention, multiple features may sometimes be grouped into one embodiment, drawing or description thereof in the foregoing description of the embodiments of the present invention.

Claims

1. A method for simulating the resetting process of a stamping die for irregularly shaped plates by combining data modeling, characterized in that, The method includes: Based on the three-dimensional morphological data of the irregular plate and the motion characteristic data of the mold reset component, a dynamic topology model for irregular plate reset is constructed. The dynamic topology model for irregular plate reset is associated with the contact area of ​​the irregular plate and the force transmission path of the mold reset through topology nodes. Collect morphological deformation data and real-time mold reset force data during the resetting process of the irregular plate, and generate a morphological force coupling data chain. The morphological force coupling data chain records the synchronous change information of deformation and corresponding reset force in each contact area of ​​the irregular plate. The morphological force coupling data chain is input into the dynamic topology model of the irregular plate reset, which drives the dynamic topology model of the irregular plate reset to perform dynamic reorganization of topology nodes and generate topology reorganization data. The topology reorganization data reflects the node connection relationship after the deformation of the irregular plate and the reset force are adapted. Based on the topology reorganization data, dynamic adaptation reset path data is generated. The dynamic adaptation reset path data includes the motion trajectory and action sequence of each component of the mold after adjustment according to the topology node relationship. The dynamic adaptation reset path data is input in reverse into the dynamic topology model of the irregular plate reset, and the association strength of the topology nodes and the reset force transmission coefficient are iteratively optimized to generate the complete simulation data after iterative optimization.

2. The simulation method for the resetting process of a stamping die for irregularly shaped plates combined with data modeling according to claim 1, characterized in that, The dynamic topology model for resetting the irregularly shaped plate is constructed based on the three-dimensional morphological data of the irregularly shaped plate and the motion characteristic data of the mold reset component, including: Collect the three-dimensional contour data, thickness distribution data and material elasticity data of the irregularly shaped plate, and generate the three-dimensional morphological data of the irregularly shaped plate through three-dimensional modeling technology. The three-dimensional morphological data of the irregularly shaped plate is used to present the spatial structure and physical parameters of the irregularly shaped plate. Collect motion stroke data, rotation angle data, and component surface contact parameter data of each reset component of the mold to form motion characteristic data of the mold reset component. The motion characteristic data of the mold reset component is used to reflect the motion capability and contact parameter information of the reset component. Extract the data of protruding areas, recessed areas and edge contours from the three-dimensional morphological data of the irregular plate. Determine the surface area of ​​the irregular plate corresponding to the data of protruding areas, recessed areas and edge contours from the three-dimensional morphological data of the irregular plate as the candidate area of ​​the topology node. Each candidate area of ​​the topology node corresponds to the key position where the irregular plate may contact the mold. Extract the component contact end face data, motion direction data, and force output data from the motion characteristic data of the mold reset component. Determine the component action area corresponding to the component contact end face data, motion direction data, and force output data from the motion characteristic data of the mold reset component as the topology node matching area. Each topology node matching area corresponds to a key location where the mold may apply a reset force. Establish the initial association relationship between the candidate region of topological nodes and the matching region of topological nodes, and assign the initial association weight according to the overlap between the spatial position of the irregular plate and the motion trajectory of the mold to form an initial topological association matrix; Based on the initial topology association matrix, each candidate region of the topology node and the corresponding matching region of the topology node are defined as the basic topology node of the irregular plate reset dynamic topology model. Each basic topology node contains spatial coordinate parameters and initial association weight parameters. Construct a force transmission path model between basic topology nodes, and set the path transmission coefficient according to the elastic data of the irregular plate material and the contact parameter data of the mold component. The path transmission coefficient is directly related to the force transmission efficiency of adjacent basic topology nodes. The basic topology nodes, initial topology association matrix and force transmission path model are integrated to form the core structure of the irregular plate reset dynamic topology model. The core structure of the irregular plate reset dynamic topology model is used to define the node composition and association rules. Input the material elastic parameters from the three-dimensional morphological data of the irregular plate and the force output parameters from the motion characteristic data of the mold reset component, fill in the node physical properties and path transmission parameters in the core structure of the dynamic topology model of the irregular plate reset, and form a preliminary dynamic topology model of the irregular plate reset. By using static contact test data between the irregular plate and the mold, the association weights of the basic topology nodes and the force transmission path coefficients in the preliminary dynamic topology model for irregular plate resetting are adjusted to generate the final dynamic topology model for irregular plate resetting.

3. The simulation method for the resetting process of a stamping die for irregularly shaped plates combined with data modeling according to claim 1, characterized in that, The process of collecting morphological deformation data and mold reset force real-time data during the repositioning of the irregularly shaped plate generates a morphological-force coupled data chain, including: Miniature deformation sensors are installed in the candidate regions of topological nodes on the surface of the irregular plate. The sensor acquisition points are fixed according to the spatial distribution of each candidate region of topological nodes, so that each acquisition point corresponds precisely to a candidate region of topological nodes. Pressure sensors are installed in the topological node matching area of ​​the mold, and the sensor detection surfaces are arranged according to the force output direction of each topological node matching area so that the detection surfaces are completely aligned with the contact direction of the irregular plate. The data acquisition programs of the deformation sensor and pressure sensor are started simultaneously, and a unified acquisition time interval is set. The acquisition time interval is set according to the mold reset movement speed, and synchronous data of each movement moment is captured. Real-time deformation data, deformation rate data, and deformation recovery data of each candidate region of the topology node are collected by deformation sensors to form a set of deformation data of the irregular plate shape. The deformation information of each candidate region of the topology node is arranged in the order of collection time in the set of deformation data of the irregular plate shape. Real-time force data, force change rate data, and force transmission duration data of each topological node matching area are collected by pressure sensors to form a real-time mold reset force data set. The real-time mold reset force data set is arranged in the order of collection time, with the force information of each topological node matching area arranged in the order of collection time. Extract deformation data of all candidate regions of topological nodes at the same acquisition time from the deformation data set of the irregular plate shape, and generate instantaneous deformation data frames. Each instantaneous deformation data frame contains the complete deformation parameters of all candidate regions of topological nodes at that acquisition time. Extract the force data of all topological node matching areas at the same acquisition time point from the real-time data set of mold reset force, and generate instantaneous force data frames. Each instantaneous force data frame contains the complete force parameters of all topological node matching areas at that time point. The instantaneous deformation data frame is associated and bound with the corresponding instantaneous force data frame. All associated and bound frame data are connected sequentially according to the acquisition time to form a preliminary form-force coupling data chain. The time stamp information and node corresponding identifier of each frame of data in the preliminary morphological force coupling data chain are supplemented, and the node to which each deformation data and force data belongs and the acquisition time are marked to generate a complete morphological force coupling data chain. The morphological force coupling data chain continuously records the synchronous change information of deformation and corresponding reset force of each contact area of ​​the irregular plate in chronological order. Each data node includes spatial position, deformation parameters, force parameters and time parameters.

4. The simulation method for the resetting process of a stamping die for irregularly shaped plates combined with data modeling according to claim 1, characterized in that, The step of inputting the morphological force coupling data chain into the irregular plate reset dynamic topology model to drive the irregular plate reset dynamic topology model to perform dynamic topology node reorganization and generate topology reorganization data includes: The instantaneous deformation data frame and instantaneous force data frame in the morphological force coupling data chain are analyzed in chronological order to extract the actual deformation of each candidate region of the topology node and the actual force of the corresponding matching region of the topology node. Calculate the deformation force adaptation coefficient for each basic topology node. The deformation force adaptation coefficient is the ratio of the actual deformation of the candidate region of the corresponding topology node to the actual force of the matching region of the topology node, reflecting the degree of adaptation between deformation and force at the node level. Input the deformation force adaptation coefficients of all basic topology nodes into the irregular plate reset dynamic topology model, compare them with the adaptation coefficient threshold initially set in the irregular plate reset dynamic topology model, and filter out abnormal nodes whose adaptation coefficients exceed the threshold range. For abnormal nodes, extract continuous multi-frame deformation and force change data corresponding to the abnormal nodes in the morphological force coupling data chain, analyze the deformation trend and force change trend of the abnormal nodes, and locate the core cause of the node abnormality. Based on the deformation and force change trends of abnormal nodes, the association weights between abnormal nodes and adjacent nodes are adjusted. The association strength of nodes is adjusted according to the degree of fit between deformation and force. When the degree of fit between deformation and force meets the preset standard, the association strength is increased; when it does not meet the preset standard, the association strength is decreased. Reconstruct the force transmission path corresponding to the abnormal node, adjust the path transmission direction according to the deformation trend, so that the force transmission path corresponds to the actual deformation direction of the irregular plate, and optimize the path transmission coefficient to match the actual force transmission efficiency. For normal nodes that do not exceed the adaptation coefficient threshold, the associated weights and path transmission coefficients of the nodes are finely adjusted according to the actual data in the morphological force coupling data chain to make the node parameters fit the actual working conditions. Record all node association weight adjustment data, force transmission path reconstruction data, and parameter fine-tuning data to form node adjustment details. Integrate the node adjustment details to generate the node association matrix and force transmission path model after dynamic reorganization of topology nodes, and label the connection relationship and transmission rules of the reorganized basic topology nodes. The recombined node association matrix, force transmission path model, and all node adjustment parameters are integrated to generate topology recombined data, which is used to reflect the node connection relationship after the deformation of the irregular plate and the reset force are adapted.

5. The simulation method for the resetting process of a stamping die for irregularly shaped plates combined with data modeling according to claim 1, characterized in that, The generation of dynamically adapted reset path data based on topology reorganization data includes: Analyze the node association matrix in the topology reorganization data, extract the spatial coordinate change data of the candidate area of ​​the irregular plate topology node corresponding to the matching area of ​​each mold topology node, and obtain the target position adjustment information of the action area of ​​each component. Extract the force transmission path model from the topology reorganization data, determine the force transmission priority of each mold component, and sort the components according to the path transmission coefficient values ​​to form the force execution order of the components. By combining the target position adjustment information and the force execution sequence, an initial motion trajectory is planned for each mold reset component. The initial motion trajectory is set according to the shortest path from the current position of the component to the target position, while incorporating the force transmission priority constraint. Extract timestamp information from the morphological force coupling data chain, determine the time node for adjusting the association relationship of each topological node, and use the time node for adjusting the association relationship of each topological node as the segment switching point of the mold motion trajectory; The initial motion trajectory is divided into multiple continuous trajectory segments according to the segment switching points. Each trajectory segment corresponds to the component motion path within a time interval, so that the relationship between the trajectory segments and the nodes is adjusted synchronously. Motion speed parameters are assigned to each trajectory segment. The motion speed parameters are set according to the force transmission coefficient values ​​within that time interval. A correspondence is established between the force transmission coefficient values ​​and the motion speed values ​​according to a preset ratio. When the force transmission coefficient values ​​are in a preset high range, the motion speed values ​​are set to the corresponding range values ​​according to the preset ratio. When the force transmission coefficient values ​​are in a preset low range, the motion speed values ​​are set to the corresponding range values ​​according to the preset ratio. By integrating all trajectory segments with their corresponding velocity parameters and time intervals, dynamic motion trajectory data of each mold reset component is formed. The dynamic motion trajectory data of each mold reset component clearly defines the motion path and velocity of the component in different time intervals. Based on the dynamic motion trajectory data of each mold reset component and the force transmission priority, the motion timing of each mold component is coordinated; the coordination of the motion timing of each mold component includes: generating mutually exclusive motion periods of the components in the overlapping area based on the spatial position overlap area data in the dynamic motion trajectory data of adjacent components; and generating serialized start instructions for the components based on the force transmission priority. Extract key position points, motion speed change points, and action switching points from the dynamic motion trajectory data of each mold reset component to form a set of core parameters for dynamically adapting the reset path; By integrating the dynamic motion trajectory data, motion timing coordination data, and core parameter set of each mold component, dynamic adaptation reset path data is generated. The dynamic adaptation reset path data includes the motion trajectory and action timing of each mold component after adjustment according to the topological node relationship.

6. The simulation method for the resetting process of a stamping die for irregularly shaped plates combined with data modeling according to claim 1, characterized in that, The process involves inputting the dynamically adapted reset path data back into the dynamic topology model of the irregular plate reset, iteratively optimizing the topology node association strength and reset force transmission coefficient, generating iteratively optimized complete simulation data, and outputting the iteratively optimized complete simulation data, including: Extract the motion trajectory coordinate data and action timing data of each component of the mold from the dynamic adaptation reset path data, and convert the motion trajectory coordinate data and action timing data of each component of the mold into node input parameters that can be recognized by the dynamic topology model of the irregular plate reset. The converted node input parameters are then input into the irregular plate reset dynamic topology model in reverse order, replacing the original node motion parameters and timing parameters in the irregular plate reset dynamic topology model. Extract the initial value of node association strength from the topology reorganization data and compare it with the actual association strength data corresponding to the component motion trajectory in the dynamic adaptation reset path data; Based on the comparison results, the topology node association strength parameters in the dynamic topology model of the irregular plate reset are adjusted so that the topology node association strength parameters correspond to the actual association strength data. Extract the initial value of the reset force transmission coefficient from the topology reorganization data and compare it with the actual transmission coefficient data corresponding to the force transmission efficiency in the dynamic adaptation reset path data; Based on the comparison results, the reset force transmission coefficient in the dynamic topology model of the irregular plate reset is optimized so that the reset force transmission coefficient corresponds to the actual transmission coefficient data. Record the adjustment data of the topology node association strength parameter and the optimization data of the reset force transmission coefficient to form an iterative optimization detail. Input the iterative optimization detail into the irregular plate reset dynamic topology model to drive the irregular plate reset dynamic topology model to recalculate the node association relationship and force transmission path. Extract the recalculated node association data, force transmission path data, and component motion trajectory optimization data, and output the iteratively optimized complete simulation data. The iteratively optimized complete simulation data includes the adjusted topology node parameters, reset force transmission parameters, and mold motion path data.

7. The simulation method for the resetting process of a stamping die for irregularly shaped plates combined with data modeling according to claim 2, characterized in that, The process involves collecting 3D contour data, thickness distribution data, and material elasticity data of the irregularly shaped plate, and then generating 3D morphological data of the irregularly shaped plate using 3D modeling technology, including: A three-dimensional laser scanning device was used to scan the entire surface of the irregular plate, and the spatial coordinate data of each surface point of the irregular plate were recorded to form the original three-dimensional contour data. The thickness of the irregularly shaped plate is measured in different areas by an ultrasonic thickness gauge, and the thickness values ​​are recorded according to the measurement location to form the original thickness distribution data. The elasticity of the irregularly shaped plate was tested using a materials mechanics testing device. The elastic modulus and Poisson's ratio data of the material were recorded to form the original material elasticity data. The original 3D contour data is denoised to remove abnormal point data generated during the scanning process, resulting in purified 3D contour data. The original thickness distribution data is interpolated to supplement the thickness values ​​of the measurement gap area, thus forming continuous thickness distribution data. The original material elasticity data is averaged to eliminate numerical fluctuations caused by experimental errors and form standard material elasticity data. Import the purified 3D contour data, continuous thickness distribution data, and standard material elasticity data into the 3D modeling software; A three-dimensional solid model of the irregular plate is constructed using 3D modeling software, and thickness parameters and material elasticity parameters are assigned to the corresponding areas in the 3D solid model. The three-dimensional solid model is meshed to form mesh model data containing spatial coordinates and physical parameters; Spatial structure information and physical parameter information are extracted from the mesh model data and integrated to generate three-dimensional morphological data of the irregular plate. The three-dimensional morphological data of the irregular plate fully presents the spatial structure and physical parameters of the irregular plate.

8. The simulation method for the resetting process of a stamping die for irregularly shaped plates combined with data modeling according to claim 3, characterized in that, The process involves collecting real-time deformation data, deformation rate data, and deformation recovery data for each candidate region of a topological node using deformation sensors, forming a set of deformation data for the irregularly shaped plate, including: Set the sampling frequency of the deformation sensor to be consistent with that of the pressure sensor to keep the data acquisition synchronized. The initial morphological data of the candidate region of the topological node is captured by the detection unit of the deformation sensor and used as the deformation comparison benchmark data; During the resetting process of the irregular plate, the real-time morphological data of each candidate region of the topology node is continuously captured, and the difference is calculated with the deformation comparison benchmark data to obtain the real-time deformation data. Based on real-time deformation data collected in two consecutive intervals, the deformation rate data is obtained by calculating the change in deformation per unit time. After the reset action stops, morphological recovery data of the candidate region of the topology node is continuously collected, and process data from the end of the reset to the morphological stabilization is recorded to obtain deformation recovery data; Assign a unique data identifier to each candidate region of the topology node, and associate the data identifier with the corresponding real-time deformation data, deformation rate data, and deformation recovery data; The real-time deformation data, deformation rate data, and deformation recovery data corresponding to each data identifier are sorted according to the acquisition time to form a deformation data sequence for a single node. Collect individual node deformation data sequences from all candidate topological node regions to form a preliminary set of irregular plate morphological deformation data; Supplement the acquisition timestamp and node location information of each data sequence to improve the association information of the irregular plate morphology deformation data set and generate a complete irregular plate morphology deformation data set. The irregular plate morphology deformation data set arranges the deformation information of each topological node candidate region in the order of acquisition time.

9. The simulation method for the resetting process of a stamping die for irregularly shaped plates combined with data modeling according to claim 4, characterized in that, The force transmission path corresponding to the reconstructed abnormal node is adjusted according to the deformation trend to ensure that the force transmission path corresponds to the actual deformation direction of the irregular plate. The path transmission coefficient is optimized to match the actual force transmission efficiency, including: Extract deformation trend data from abnormal nodes to determine the main direction and magnitude of deformation. Based on the data of the main deformation direction, a new force transmission path direction is planned so that the new path direction forms a complementary angle with the main deformation direction; Based on the deformation amplitude data, the number of branches in the force transmission path is set. The deformation amplitude data and the number of branches in the path are allocated according to a preset correspondence, and the number of branches is determined according to the actual deformation requirements. Extract the location data of adjacent normal nodes of abnormal nodes, and use the adjacent normal nodes as the connection nodes of the new force transmission path; Construct path segments between abnormal nodes and adjacent normal nodes to form a preliminary reconstruction force transmission path; Calculate the length and resistance data of the initial reconfiguration force transmission path, and set the initial path transmission coefficient based on the length and resistance data of the initial reconfiguration force transmission path; Extract the actual force transmission efficiency data corresponding to the abnormal nodes in the morpho-force coupling data chain; Adjust the initial path transmission coefficient so that the transmission efficiency corresponding to the adjusted path transmission coefficient corresponds to the actual force transmission efficiency data; Record the direction data, branch data, connection node data, and optimized path transmission coefficient of the reconstructed force transmission path to generate complete force transmission path reconstruction data for abnormal nodes. The force transmission path reconstruction data for abnormal nodes includes the adjusted path transmission direction and the path transmission coefficient that matches the actual force transmission efficiency.

10. A simulation system for the resetting process of a stamping die for irregularly shaped plates, combining data modeling, characterized in that, include: processor; A machine-readable storage medium for storing machine-executable instructions of the processor; The processor is configured to execute the simulation method for resetting irregular plate stamping dies by incorporating data modeling as described in any one of claims 1 to 9 by executing the machine-executable instructions.

Images