Enterprise production information management sharing method and system based on digital twinning

Abstract

The invention relates to the technical field of production information processing, in particular to an enterprise production information management sharing method and system based on digital twinning. The method comprises the following steps: collecting a multi-view three-dimensional point cloud in a to-be-wired product to construct a digital twin model of the product; executing three-dimensional path optimization within constraints according to cable connection starting and ending coordinates in the digital twin model of the product, and then obtaining a digital cable guide path limited by physical bending of the cable; and in the wiring operation process, a wiring conformity index is calculated and shared through a digital cable guiding path, so that dynamic guiding of wiring and management sharing of production information are realized. According to the method, the digital twin guide path is constructed and the wiring conformity is calculated in real time, so that the wiring quality consistency is improved and the production information sharing is realized.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of production information processing, and in particular to an enterprise production information management sharing method and system based on digital twinning. BACKGROUND

[0002] The existing computer case wiring technology mainly relies on static two-dimensional wiring diagrams and textual instructions, which cannot intuitively present the actual layout and hierarchical relationship of cables in three-dimensional space. It is difficult for operators to accurately understand complex wiring requirements, especially when dealing with the spatial layout and mutual relationship of multiple types of cables (such as power lines, data lines, control lines, etc.), which is prone to understanding deviations. Moreover, after the operators complete the wiring, the quality inspection process mainly relies on manual visual inspection, lacking objective and quantitative evaluation standards and tools, making it difficult to find wiring defects in hidden areas, such as excessive bending, tight or loose bundling of cables, etc., which may cause product quality problems. In particular, in a mass production environment, the wiring quality of different operators varies significantly, making it difficult to ensure product consistency, resulting in insufficient wiring quality and efficiency, and failing to achieve case wiring assembly guidance and wiring quality information sharing. SUMMARY

[0003] Therefore, the present application provides an enterprise production information management sharing method and system based on digital twinning to solve at least one of the above technical problems.

[0004] To achieve the above purpose, an enterprise production information management sharing method based on digital twinning is applied to computer case wiring operations, which comprises the following steps:

[0005] Step S1: Collecting multi-view three-dimensional point clouds inside the product to be wired to construct a product digital twinning model;

[0006] Step S2: Mapping the physical constraints in the product digital twinning model into heat dissipation duct wiring constraints and rotating area wiring constraints, and performing three-dimensional path optimization according to the cable connection start and end coordinates in the product digital twinning model within the constraints, and marking them as ideal cable paths;

[0007] Step S3: Identifying and correcting path segments that have geometric interference with the ideal cable path and the preset cable bending limit constraint, thereby generating a digital cable guide path with cable physical bending limitation;

[0008] Step S4: During the wiring operation, collecting three-dimensional data in the wiring operation area to obtain the actual wiring path, comparing the actual wiring path with the digital cable guide path, calculating and sharing the wiring compliance index to realize dynamic guidance of wiring and management sharing of production information.

[0009] The application also provides an enterprise production information management sharing system based on digital twinning, which executes the enterprise production information management sharing method based on digital twinning as described above, and comprises:

[0010] A digital twinning modeling module is configured to collect multi-view three-dimensional point clouds inside a product to be wired to construct a product digital twinning model.

[0011] A wiring path optimization module is configured to map physical constraints in the product digital twinning model into heat dissipation duct wiring constraints and rotating area wiring constraints, and perform three-dimensional path optimization in the constraints according to cable connection start and end coordinates in the product digital twinning model, and mark the ideal cable path.

[0012] A guide path generation module is configured to perform cable bending limit constraint on the ideal cable path, correct all geometric interference path segments in the ideal cable path in a full range of the ideal cable path, and obtain a digital cable guide path with cable bending limit.

[0013] A work information sharing module is configured to collect three-dimensional data in a wiring work area to obtain an actual wiring path during a wiring work process, compare the actual wiring path with the digital cable guide path, calculate and share a wiring compliance index, and realize dynamic guidance of wiring and management sharing of production information.

[0014] The application has the following beneficial effects:

[0015] On the one hand, the digital twinning model is constructed by collecting multi-view three-dimensional point clouds inside a product to be wired, and a digital cable guide path is generated, so that the traditional static and abstract two-dimensional wiring diagram and text description are upgraded to dynamic and intuitive three-dimensional visual work guidance. This process is based on accurate reproduction of a real three-dimensional space of a product, can clearly present actual movements, spatial levels and avoidance relationships of various cables in a complex structure, and enables an operator to accurately understand wiring requirements, significantly reduces wiring errors caused by understanding deviation, and effectively improves the accuracy of first wiring, especially when various cables such as power lines and data lines are arranged in an interleaved manner.

[0016] In another aspect, the physical constraints such as heat dissipation air duct, rotating area, etc. are mapped as wiring constraints in the digital twin model, and three-dimensional path optimization is performed between the start and end coordinates of the cable connection to ensure that the generated ideal cable path is not only the optimal connection path, but also complies with the core design requirements of the product (such as heat dissipation performance, structure interference avoidance, etc.) from the source; In addition, by identifying and correcting the path segments that geometrically interfere with the physical bending limit of the cable, it can proactively avoid quality risks caused by excessive bending and stress concentration of the cable, accurately identify and correct subtle path deviations caused by tight or loose cable bundling, etc. Hidden defects ensure the high consistency and reliability of product wiring quality between different batches and different operators. BRIEF DESCRIPTION OF DRAWINGS

[0017] Fig. 1 Step flowchart of the enterprise production information management sharing method based on digital twin of the present application;

[0018] Fig. 2 Module schematic diagram of the enterprise production information management sharing system based on digital twin of the present application;

[0019] Fig. 3 Case wiring example diagram of the present application;

[0020] The implementation, functional features and advantages of the present application will be further described with reference to the embodiments and the accompanying drawings. DETAILED DESCRIPTION

[0021] The technical method of the present application will be described clearly and completely below in combination with the accompanying drawings. Obviously, the described embodiments are part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present application.

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

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

[0024] To achieve the above object, please refer to Figs. 1 to 3 The application provides an enterprise production information management sharing method based on digital twinning, comprising the following steps:

[0025] Step S1: Collecting multi-view three-dimensional point clouds inside the product to be wired to construct a product digital twinning model;

[0026] In the embodiment of the application, a 2U rack-mounted server case is taken as the product to be wired. Before the collection starts, the empty case is fixed on a three-dimensional automatic rotating table. An industrial-grade structured light three-dimensional scanner with a scanning accuracy of 0.05 mm is used to collect multi-angle and high-density point cloud data of the case. The scanning process covers the top, bottom, front, back, left and right six main viewing angles of the case, and in the interior of the case, at least six different inclined angle local fine scanning is additionally performed for the key wiring areas such as the power supply compartment, the mainboard mounting area and the hard disk backplane, and a total of not less than 12 viewing angle scanning operations are completed.

[0027] In one implementation manner of the embodiment of the application, it is assumed that the size of the 2U server case is 482 mm x 88 mm x 700 mm. The case is placed on the rotating table, and first, global scanning is performed, the rotating table is rotated by 30 degrees each time, the scanner collects data once, and a total of 12 groups of global point clouds are obtained. Subsequently, the scanner probe is inserted into the interior of the case and is aligned with the CPU power supply port of the mainboard area, the dense area of the SATA interface and the like, and scanning is performed from-45°, 0° and +45° three pitch angles. After all the scanning is completed, a total of about 800 million original three-dimensional coordinate points are collected. The point cloud processing software (such as Geomagic Wrap) is used to automatically register and fuse the point cloud data under the multi-view, and redundant and noise points are removed, and finally a highly accurate unified three-dimensional point cloud model containing all structural features (such as screw holes, cable tie fixing points, sheet metal folded edges, component outlines and the like) in the interior of the case is generated. Further, the point cloud model is converted into a polygon mesh model (Mesh Model) with a topological structure through a surface reconstruction technology, and finally a digital twinning model corresponding to the physical case 1:1 is generated, and the model is stored in a STEP format file, thereby providing an accurate geometric basis for subsequent path planning.

[0028] Step S2: Map the physical constraints in the product digital twin model into the heat dissipation air duct wiring constraints and the rotating area wiring constraints, and perform three-dimensional path optimization within the constraints according to the cable connection start and end coordinates in the product digital twin model, and mark it as an ideal cable path;

[0029] In the embodiment of the application, based on the server chassis digital twin model, the key physical constraints in the model are defined digitally in the three-dimensional design software environment (such as SOLIDWORKS or CATIA). First, three 80mm heat dissipation fans in the chassis are identified, and the space swept by the rotation of the blades is defined as a cylinder with a diameter of 82mm (1mm safety margin is reserved) and a height of 25mm, which is the rotating area wiring constraint, and any cable path cannot enter this area. Secondly, according to the server design cooling scheme, a main air flow channel from the front panel hard drive cage to the rear CPU heat sink area is defined, which is abstracted as a cuboid space with a size of 400mm×150mm×20mm, which is the heat dissipation air duct wiring constraint, and the cable should be avoided to cross this core air duct to ensure the cooling efficiency.

[0030] Step S3: Identify and correct the path segment of the ideal cable path that geometrically interferes with the pre-set cable bending limit constraint, thereby generating a digital cable guide path with physical bending limit of the cable;

[0031] In the embodiment of the application, taking the example of laying a 24-Pin ATX main power cable from the power supply unit (PSU) to the mainboard, the physical properties of the cable are defined in the process database in advance, wherein the cable bending limit constraint is set to be not less than 25mm (determined based on the cable diameter and internal conductor structure) The minimum bending radius of the cable is analyzed segment by segment along the coordinate point sequence of the ideal cable path, and whether there is excessive bending is identified by calculating the curvature radius of each point on the path.

[0032] In an implementation manner of the embodiment of the present application, it is assumed that the ideal cable path is detected to make a sharp turn when bypassing a chassis sheet metal column, and the three continuous coordinate points P10 (180, 50, 410), P11 (185, 50, 405), and P12 (190, 55, 405) constitute the sharp turn. Through geometric calculation, it is determined that the curvature radius of the circular arc passing through the three points is 18 mm. Since 18 mm is less than the preset bending limit constraint of 25 mm, the path segment is determined as a “geometrically interfered path segment”. Then, the path correction program is triggered: taking P10 and P12 as control points, a new coordinate point is automatically inserted in the path segment to generate a smooth curve (circular arc segment) that satisfies the minimum bending radius limit to replace the original P10-P11-P12 polyline segment. The curvature radius of any point on the corrected new path segment is greater than or equal to 25 mm. After checking and correcting all such geometrically interfered points on the entire ideal cable path, the final path obtained is the digital cable guide path that takes into account the space obstacle avoidance and the physical limitations of the cable and can be followed by the operator in practice.

[0033] Step S4: During the wiring operation, three-dimensional data in the wiring operation area is collected to obtain an actual wiring path, and by comparing the actual wiring path with the digital cable guide path, a wiring compliance index is calculated and shared to realize dynamic guidance of wiring and management and sharing of production information.

[0034] In an implementation manner of the embodiment of the present application, when the operator completes the preliminary fixation of the 24-Pin cable, the actual wiring path obtained in real time is compared and analyzed with the digital cable guide path displayed in the AR glasses. The iterative closest point (ICP) algorithm is used to calculate the average distance deviation between the two paths. It is assumed that the average deviation calculated is 4.5 mm, and the preset maximum allowed deviation is 10 mm. According to the formula: wiring compliance index = (1- average deviation / maximum allowed deviation) x 100%, the current compliance is calculated as (1-4.5 / 10) x 100% = 55%. Since the value is lower than the qualified threshold of 95%, the cable segment with a deviation of more than 5 mm is highlighted in red in the AR glasses, and a dripping prompt sound is emitted, dynamically guiding the operator to adjust the position of the cable segment. After the operator adjusts, the compliance is recalculated to be 98.2%. At this time, the cable path in the AR becomes green, indicating that the wiring is qualified. Finally, the wiring compliance index of 98.2%, together with the operator ID, wiring time, product serial number, and other information, is packaged into a piece of production quality data, which is uploaded to the manufacturing execution system (MES) of the enterprise in real time, realizing the management and sharing of production information among the design, assembly, and quality inspection links.

[0035] Preferably, before mapping the physical constraints in the product digital twin model into the heat dissipation air duct wiring constraints and the rotating area wiring constraints in step S2, it further comprises the following steps:

[0036] Extract the contour boundary coordinates of the mainboard, power supply and hard disk bracket in the product digital twin model, identify the CPU cooler top plane and the chassis rear fan center point coordinates, and calculate the line connecting the two points as the main air duct axis direction vector;

[0037] Determine the air duct safety envelope according to the main air duct axis direction vector;

[0038] Detect the installation position and blade radius size of all fans in the chassis, draw a circular exclusion area with the fan center as the center and the blade radius + 15mm as the safety radius, and mark it as a rotating envelope.

[0039] In the embodiment of the application, before performing three-dimensional path optimization, the constructed server chassis digital twin model is subjected to automatic physical constraint identification and definition, the semantic information of the digital twin model (such as a STEP AP242 format file, which retains the assembly information and metadata of the components) is queried through an API interface. The component marked as "CPU_Cooler" in the model is automatically identified and extracted. The three-dimensional coordinates of the geometric center of the top surface of the smallest bounding box of the component are determined, which are marked as the CPU cooler top plane center point P_cpu. Subsequently, the rear exhaust fan marked as "Fan_Rear_1" is identified in the same way, and the geometric center point coordinates are obtained, which are marked as the chassis rear fan center point P_fan.

[0040] In one implementation manner of the embodiment of the application, it is assumed that the coordinate system origin of the server chassis is located at the lower left corner of the front end, and after analyzing the digital twin model, it is obtained that the coordinates of the CPU cooler top plane center point P_cpu are (210, 85, 300) mm; the coordinates of the chassis rear 120mm exhaust fan center point P_fan are (210, 85, 690) mm, and the connecting vector between the two points is calculated as the main air duct axis direction vector V_airflow: V_airflow = P_fan - P_cpu = (0, 0, 390), which indicates that the main air duct is a straight air flow channel along the depth direction (Z axis positive direction) of the chassis, straightly pointing to the rear exhaust fan from the CPU area. Next, the air duct safety envelope is determined according to the main air duct axis direction vector. The envelope is a cuboid space with a preset width and height, with the line connecting P_cpu and P_fan as the central axis. The width and height values are preset according to the thermal simulation report of the chassis design, aiming to protect the core air flow area from cable interference.

[0041] In another implementation form of the embodiment of the application, based on the industry heat dissipation design specification, the preset air duct safety width of the 2U server case is 180 mm, and the safety height is 50 mm. An oriented bounding box with a length of 390 mm, a width of 180 mm, and a height of 50 mm is constructed with the vector V_airflow as the orientation. The geometric center of the cuboid coincides with the midpoint of the line connecting P_cpu and P_fan, and the longest side of the cuboid is parallel to V_airflow. This cuboid space is marked as “heat dissipation air duct wiring constraint” and will be regarded as a high-cost or impassable area in subsequent path optimization.

[0042] It should be noted that the oriented envelope body constructed based on the main air duct axis direction vector can more accurately match the actual air duct layout that is inclined or non-orthogonal, thereby releasing more available space for wiring in non-critical areas while ensuring heat dissipation efficiency.

[0043] At the same time, start to automatically detect all fan components in the case to define the rotating envelope, traverse all components in the digital twin model, and filter out the components with the metadata label “Fan”. For each identified fan component, read its assembly position (i.e., the center point coordinates) and the “Blade_Radius” (blade radius) parameter pre-stored in its attributes.

[0044] In one implementation form of the embodiment of the application, it is assumed that a total of 3 fan components are detected in the case: “Fan_Rear_1” located at the rear and “Fan_Front_1” and “Fan_Front_2” located at the front end. The parameters are read as follows: Fan_Rear_1: center point (210, 85, 690), blade radius 58 mm (corresponding to a 120 mm fan); Fan_Front_1: center point (140, 60, 25), blade radius 38 mm (corresponding to an 80 mm fan); and Fan_Front_2: center point (280, 60, 25), blade radius 38 mm (corresponding to an 80 mm fan). Take the center of each fan as the center of the circle, and take “blade radius + 15 mm” as the safety radius to generate a cylindrical exclusion area, i.e., the rotating envelope. The height of the cylinder is equal to the thickness of the fan (for example, 25 mm). The safety radius of Fan_Rear_1 is calculated as 73 mm. The safety radius of Fan_Front_1 and Fan_Front_2 is calculated as 38 mm + 15 mm = 53 mm. Three cylindrical exclusion areas are generated in the model, and they are uniformly marked as rotating area wiring constraints.

[0045] It is important to note that the safety radius of 15mm is a critical engineering preset, which aims to: 1) provide a buffer for the slight deformation or axial movement of the fan blades when rotating at high speed; 2) prevent the cable from being sucked into the rotating blades due to the air flow negative pressure effect (Venturi effect); 3) provide the necessary safety operation distance for the fixation of the cable in this area (such as using a cable tie).

[0046] Especially important is that determining the air duct safety envelope according to the main air duct axis direction vector includes:

[0047] Taking the top plane profile of the CPU cooler as the starting section and the fan profile at the rear of the case as the terminal section, stretch along the main air duct axis direction vector to generate an initial air duct corridor;

[0048] Perform a spatial Boolean operation in the initial air duct corridor to subtract all solid entities in the product digital twin model except the cooler and the fan to obtain an effective air duct space;

[0049] Calculate the minimum cross-sectional area of the effective air duct space along the axis direction, and construct the air duct safety envelope with the effective air duct space and the minimum cross-sectional area.

[0050] In the embodiment of the application, two key cross-sectional profiles are accurately identified and extracted in the digital twin model: the starting section is positioned to the CPU cooler component, and the two-dimensional profile of the topmost cooling fin is extracted, which is a rectangle. The terminal section is positioned to the 120mm exhaust fan at the rear of the case, and the two-dimensional profile of the circular air vent is extracted.

[0051] In one implementation manner of the embodiment of the application, by analyzing the digital twin model, the starting section is obtained as a rectangular profile with a center at (210, 85, 300) and a size of 100mm x 80mm; the terminal section is obtained as a circular profile with a center at (210, 85, 690) and a radius of 60mm. Call the lofting (Loft) function of the three-dimensional modeling engine, take the two sections as input, and stretch along the previously calculated main air duct axis direction vector V_airflow = (0, 0, 390) to generate a smooth transition three-dimensional entity. The entity smoothly transitions from a rectangular base to a circular top, and its shape accurately covers the theoretical core airflow area from the CPU cooler to the rear exhaust fan, which is the initial air duct corridor.

[0052] In another implementation manner of the embodiment of the application, in the space covered by the initial air duct corridor, other solid entities are identified, such as four DDR5 memory sticks and a large capacitor on the motherboard. Perform a large-scale spatial Boolean operation, specifically a subtraction operation, which can be abstracted as:

[0053] S_effective = S_initial - (∑S_component);

[0054] Where: S_effective represents the final effective air duct space; S_initial represents the three-dimensional entity of the initial air duct corridor; ∑S_component represents the geometric union of all other solid entities (such as memory sticks, capacitors, motherboard heat sinks, etc.) located inside the initial air duct corridor.

[0055] It should be noted that performing this Boolean operation is crucial because it converts a regular, theoretical air duct model into an irregular, but physically completely accurate, available airflow space model. For example, a memory stick will form four regular "grooves" under the side of the effective air duct space. Any subsequent cable path planning that ignores these "grooves" will incorrectly plan the cable path into the physically occupied space by the memory stick.

[0056] In an implementation manner of the embodiment of the present application, the effective air duct space is virtually "sliced" along the main air duct axis direction (Z axis) from Z = 300 mm to Z = 690 mm with a step size of 0.5 mm, and the cross-sectional area of each slice is calculated. At Z = 310 mm (just leaving the CPU heat sink), the cross-sectional area is 100 mm x 80 mm = 8000 mm 2 At Z = 365 mm, the slice just passes through the highest point of the memory stick, and the cross-sectional shape at this time is the original rectangle minus the cross-sectional area of the four memory sticks. Assuming that the size of each memory stick at this cross-section is 30 mm x 8 mm, the cross-sectional area is 8000 mm 2 - 4 x (30 mm x 8 mm) = 7040 mm 2 At Z = 450 mm, the slice passes through a tall motherboard capacitor, and the projected area of the capacitor at this cross-section is π x (5 mm) 2 ≈ 78.5 mm 2 At this time, the cross-sectional area is about 7921.5 mm 2 . By traversing all the slices, it is determined that the minimum cross-sectional area A_min occurs at the location of the memory stick, and the value is 7040 mm 2 . Finally, the air duct safety envelope is defined as a composite constraint body: its geometric shape is the effective air duct space entity obtained after Boolean operation, and it is accompanied by the minimum cross-sectional area A_min = 7040 mm 2 .

[0057] Preferably, step S2 comprises the following steps:

[0058] Step S21: determining the heat dissipation air duct wiring constraint according to the air duct safety envelope;

[0059] Step S22: determining a rotation region routing constraint according to the rotation envelope;

[0060] Step S23: determining a start point and an end point of the cable to be routed according to a preset socket position of the cable connector in the product digital twin model;

[0061] Step S24: performing three-dimensional path optimization according to the start point and the end point of the cable to be routed, and then performing repulsion value constraint through the heat dissipation air duct routing constraint and the rotation region routing constraint, so as to lock a three-dimensional space trajectory with the minimum cumulative repulsion value as a path center line;

[0062] Step S25: generating a three-dimensional tubular entity along the path center line in the product digital twin model according to the cable diameter of the cable to be routed as a cross-sectional diameter, so as to form an ideal cable path.

[0063] In the embodiment of the application, taking routing of a CPU 8-Pin power cable as an example, the connection information of the CPU 8-Pin cable is queried by accessing the component database in the digital twin model. The database records two components connected by the cable: a power supply unit (PSU) and a motherboard (Motherboard).

[0064] In one implementation manner of the embodiment of the application, by analyzing the geometric model of the PSU component, the connector marked as “CPU_POWER_OUT” is found, and the geometric center coordinates of the 8 pins thereof are extracted, and the average value thereof is calculated as the start point P_start. Similarly, the socket marked as “CPU_POWER_IN” is found on the motherboard model, and the average center coordinates thereof are calculated as the end point P_end. Assuming that P_start (power supply side) = (35.0, 70.0, 685.0) mm and P_end (motherboard side) = (180.0, 105.0, 280.0) mm are calculated, the two precise three-dimensional coordinate points will be input to the path optimization algorithm.

[0065] In the embodiment of the application, the improved rapid extended random tree (RRT) algorithm is adopted to perform three-dimensional path optimization. The algorithm randomly samples and extends the path in the free space of the case until the path connecting P_start and P_end is found. The key is that the cost function of the algorithm is designed to not only include the geometric length of the path, but also introduce the repulsion value from the routing constraint. Specifically, for any point P in the space, the repulsion value C_repulsion(P) is constructed by the following function model:

[0066] c_repulsion(P) = w_rot x C_rot(P) + w_air x C_air(P);

[0067] where w rot and w air are weight coefficients to adjust the importance of different constraints. In this embodiment, w rot is set to ∞ (to represent absolute forbidden), and w air is set to 1.5. C rot (P) is the repulsive force from the rotation region routing constraint. If point P is inside any of the rotation envelopes, C rot (P) = ∞; otherwise, C rot (P) = 0. C air (P) is the repulsive force from the air duct routing constraint, which is calculated as:

[0068] C air (P) = (1 / d(P, S air )) x (A max / A min );

[0069] where d(P, S air ) is the shortest Euclidean distance from point P to the surface of the air duct safety envelope (S air ); A max is the maximum cross-sectional area of the air duct envelope (e.g., 8000 mm 2 ), and A min is the minimum cross-sectional area (7040 mm 2 ). The ratio (A max / A min ) ≈ 1.136, as an amplification factor, means that the repulsive force will be stronger near the narrowest region of the air duct.

[0070] In one implementation of the embodiment, after the RRT algorithm generates a candidate path (composed of a series of points {P1, P2,..., Pn}), it calculates the cumulative repulsive force value C total :

[0071] C total =∑[length(Pi, Pi+1) + β x (C repulsion (P i + C repulsion (P i +1)) / 2];

[0072] where β is a tuning parameter to balance the length and repulsive force, and length(Pi, Pi+1) represents the straight-line geometric distance between two consecutive path points Pi and Pi+1 in the candidate path. The goal of the algorithm is to find a path that minimizes C total.

[0073] In one implementation of the embodiment, the database records that the cable has a diameter of 3.5 mm, and the Sweep function of the three-dimensional modeling engine is invoked to sweep along the path centerline with a circular cross-section of diameter 3.5 mm. This operation generates a smooth three-dimensional tubular entity in the digital twin model that represents the actual footprint of the CPU 8-Pin cable. It is noted that this three-dimensional tubular entity is the final ideal cable path.

[0074] Preferably, step S3 comprises the following steps:

[0075] Step S31: sliding window processing is performed on the sequence of three-dimensional coordinate points along the ideal cable path, and local curvatures within each window are obtained;

[0076] Step S32: a cable bending limit constraint is set based on the local curvatures within each window;

[0077] Step S33: all path segments in the ideal cable path that geometrically interfere with the cable bending limit constraint are taken as to-be-corrected path segments;

[0078] Step S34: three-dimensional path optimization is re-performed between the start point and the end point of each to-be-corrected path segment to determine a corrected wiring path segment;

[0079] Step S35: the corrected wiring path segment is replaced with the to-be-corrected path segment, and is spliced with the path segments in the ideal cable path that do not interfere, thereby obtaining a digital cable guide path.

[0080] In one implementation manner of the embodiment of the present application, the coordinate point sequence {P1, P2,..., Pn} of the ideal cable path (i.e., the path center line) of the generated CPU 8-Pin power cable is taken as input, and a fixed-size sliding window is moved along the sequence, which covers three continuous coordinate points (Pi-1, Pi, Pi+1) each time. For example, the size of the sliding window is set to 3 points. When the window slides to a position centered on point Pi, a circle is uniquely determined in the three-dimensional space through the three points Pi-1, Pi, and Pi+1, and the radius R_i of the circle is calculated. The local curvature K_i is defined as the reciprocal of the radius. By traversing the entire coordinate point sequence from i=2 to i=n-1, a local curvature curve corresponding to the ideal cable path is obtained, which accurately describes the bending degree of each point on the path.

[0081] In one implementation manner of the embodiment of the present application, the database specifies that the minimum allowable bending radius of the CPU 8-Pin cable is 6 times the cable diameter (3.5 mm), i.e., R_min=6×3.5mm=21mm. Therefore, the maximum allowable curvature K_max is calculated as follows: K_max=1 / R_min=1 / 21mm≈0.0476mm -1 This K_max value is set as the global cable bending limit constraint.

[0082] In one implementation manner of the embodiment of the present application, when the curvature curve is traversed, it is found that from the 35th coordinate point P35 to the 42nd coordinate point P42, the corresponding local curvatures K35 to K42 are all greater than K_max=0.0476mm -1For example, at P38, the calculated bending radius is only 15mm, and its curvature K38 = 1 / 15 = 0.0667mm, far beyond the limit value. The entire line segment from the previous point P34 (the start point of interference) to the next point P43 (the end point of interference) of P35 is marked as a to-be-modified path segment. It should be noted that when marking the to-be-modified path segment, one or several points are usually extended to both ends to ensure that the new path after modification can be smoothly connected with the original path. -1

[0083] In another implementation manner of the embodiment of the present application, in this local space, in addition to the original chassis firmware, the repulsive force source of the rotation area, and the like, a curvature-related "internal stress potential field" is introduced. The "energy" of the potential field is proportional to the curvature of the path. The potential energy function model can be set as:

[0084] U_curvature(P) = c x (max(0, K(P) - K_max)) 2 ;

[0085] Wherein, K(P) is the instantaneous curvature of the path at point P; K_max is the maximum allowable curvature; and C is a penalty coefficient. The planner starts from P34 and searches for a path to P43 with the lowest total potential energy (including the repulsive potential energy and the internal stress potential energy). Due to the existence of the potential energy function, any path that attempts to bend beyond the limit will be subject to a huge "energy penalty", so the algorithm will naturally converge to a smooth path with the curvature of all points being less than or equal to K_max. The newly generated path is the modified wiring path segment.

[0086] In one implementation manner of the embodiment of the present application, the to-be-modified part {P35,..., P42} in the coordinate point sequence {P1,..., P34, P35,..., P42, P43,..., Pn} of the original ideal cable path is deleted, and then the coordinate point sequence {P'35, P'36,..., P'm} of the generated modified wiring path segment is inserted between P34 and P43. The spliced new coordinate point sequence is {P1,..., P34, P'35,..., P'm, P43,..., Pn}. The splicing points (P34, P'35) and (P'm, P43) are smoothed (such as using spline interpolation) to ensure that the first and second derivatives of the path at the connection are continuous, that is, the path itself and the curvature change are smooth. After all the to-be-modified path segments on the ideal cable path are replaced and spliced, the final path obtained is the digital cable guide path that meets the space obstacle avoidance constraint and the physical bending limit constraint and can be directly used in the production site.

[0087] Preferably, the step S4 comprises the following steps:​

[0088] Step S41: periodically synchronously scan the real-time work area point cloud data of the wiring work area;

[0089] Step S42: extract the cable skeleton line from the real-time work area point cloud data as the actual wiring center line;

[0090] Step S43: calculate the path deviation amount based on the actual wiring center line and the digital cable guide path, wherein the path deviation amount includes the main path compliance deviation and the bending compliance deviation;

[0091] Step S44: perform weighted processing according to the main path compliance deviation and the bending compliance deviation to obtain the wiring compliance index.

[0092] In one implementation of the embodiment of the application, at the assembly station, an AzureKinect DK depth camera is fixed on the gantry above the operation table, and its field of view completely covers the wiring work area of the 2U server case. When the operator starts to lay the CPU 8-Pin power cable, the depth camera is triggered to perform synchronous scanning at a preset period (for example, 2Hz, i.e. 2 times per second), and real-time work area point cloud data containing the case, installed components, the operator's hand, and the cable being laid are captured.

[0093] In another implementation of the embodiment of the application, a semantic segmentation model based on a PointNet++ deep learning network is used. The model is trained on a large amount of labeled case wiring point cloud data and can accurately segment the point cloud clusters belonging to "cable" from dynamic point cloud. After segmentation, a three-dimensional point cloud cluster representing only the CPU 8-Pin cable is obtained. Then, the "L1-median skeletonization" algorithm is applied to the cable point cloud to calculate the center axis of the point cloud cluster, even if the point cloud has noise or partial loss. Finally, a curve composed of discrete three-dimensional coordinate points is output, which is determined as the actual wiring center line.

[0094] In one implementation of the embodiment of the application, discrete Frechet distance is used to measure the similarity of two paths in space position. This distance value is not sensitive to the local stretching, compression and slight translation of the path, and can well reflect the overall shape matching degree of the two paths. Assuming that the calculated Frechet distance is D_frecjet=8.2mm, the deviation value is obtained by normalizing the distance:

[0095] D_pos=D_frechet / L_guide;

[0096] Where L_guide is the total length of the digital cable guide path. Assuming L_guide = 450mm, then D_pos = 8.2 / 450≈0.0182.

[0097] In one implementation of the embodiment of the application, the local curvatures of each corresponding point on the actual wiring centerline and the digital cable guide path are calculated respectively to obtain two curvature curves K_actual(s) and K_guide(s), where s is the arc length parameter along the path length, and then the root mean square error (RMSE) between the two curvature curves is calculated as the bending compliance deviation:

[0098] D_curv = sqrt[∫(K_actual(s)―K_guide(s)) 2 ds / L_guide];

[0099] This value reflects the difference between the bending smoothness of the actual wiring and the ideal state. Assuming D_curv = 0.005mm -1 is calculated, the application believes that the positional accuracy of the main path is more important than the fine imitation of the bending, so the weights w_pos = 0.7 and w_curv = 0.3 are set, and the wiring compliance index Q_score is calculated by the following model:

[0100] Q_score = (1―(w_pos×D_pos+w_curv×D_curv))×100%;

[0101] It should be noted that before substituting the formula, the original deviation value needs to be normalized (D_pos, D_curv) to map it to the [0, 1] interval. Normalization is usually based on the maximum allowed deviation value. Assuming the maximum allowed value of D_pos is 0.05 and the maximum allowed value of D_curv is 0.02mm -1 .

[0102] Preferably, determining the heat dissipation air duct wiring constraint according to the air duct safety envelope comprises:

[0103] Performing a start cross-section scan on the air duct safety envelope along the main air duct axis direction vector, locking the geometric plane with the smallest cross-sectional area in the air duct safety envelope as the core air duct cross-section;

[0104] Three-dimensional voxelization of the air duct safety envelope into an air duct voxel grid;

[0105] For each voxel in the air duct voxel grid, calculate the repulsive force value according to the reciprocal of the area of the cross-section it is in and the preset air duct weight coefficient;

[0106] The repulsion value of the voxel corresponding to the intersection of the core air duct cross section in the air duct voxel grid is multiplied by a preset penalty factor to construct the heat dissipation air duct wiring constraint.

[0107] In the embodiment of the application, the main air duct axis direction vector V_airflow calculated in the above embodiment is the scanning direction, and the cross section analysis is performed on the air duct safety envelope obtained by the Boolean operation. Starting from Z=300mm (the top of the CPU radiator), a virtual scanning plane is moved along the positive direction of the Z axis at a step length of 0.5mm until Z=690mm (the rear fan). The area of the two-dimensional section formed by the intersection of the scanning plane and the air duct safety envelope is calculated at each step length position, and by comparing the area values of all sections, it is found that at Z=365mm, due to the occupation of the memory stick, the cross-sectional area reaches a minimum value of 7040mm 2 The geometric plane at Z=365mm is locked and marked as the core air duct cross section.

[0108] In an implementation manner of the embodiment of the application, the size of the voxel is set to 2mm*2mm*2mm. The air duct safety envelope is completely filled with such cubic units to form a three-dimensional grid. Any voxel with a center point located inside the envelope is retained to form an air duct voxel grid.

[0109] In another implementation manner of the embodiment of the application, the calculation model of the repulsion value is as follows:

[0110] Repulsion(V)=W_airway×(1 / Area(S_v));

[0111] Wherein, Repulsion(V) is the repulsion value of the voxel V, W_airway is a preset air duct weight coefficient for adjusting the influence weight of the air duct repulsion in the whole path optimization cost function. Here, W_airway is set to 1000 according to the process experience, and Area(S_v) is the cross-sectional area of the position of the voxel V. For example, for a voxel V_wide located at Z=310mm (the cross-sectional area is 8000mm 2 ), the repulsion value is Repulsion(V_wide)=1000×(1 / 8000)=0.125. For a voxel V_narrow located at Z=365mm (the cross-sectional area is 7040mm 2 ), the repulsion value is Repulsion(V_narrow)=1000×(1 / 7040)≈0.142.

[0112] In one implementation of the embodiment of the application, the penalty factor PF = 5.0 is set. For the previously calculated voxel V_narrow located at the core air duct cross section, its final repulsion value Repulsion_final(V_narrow) will be updated as: Repulsion_final(V_narrow) = Repulsion(V_narrow) x PF = 0.142 x 5.0 = 0.71.

[0113] Preferably, determining the rotation region routing constraint according to the rotation envelope comprises:

[0114] marking the internal three-dimensional space of the rotation envelope as an absolute forbidden region;

[0115] performing a start equidistant outward expansion operation based on the geometric outer surface of the absolute forbidden region to obtain a repulsion influence zone;

[0116] setting a three-dimensional repulsion field according to the absolute forbidden region and the repulsion influence zone, assigning any point in the field located in the absolute forbidden region with a maximum repulsion value, and assigning any point located in the repulsion influence zone according to an inverse function of its distance to the surface of the absolute forbidden region, to form the rotation region routing constraint; wherein the maximum repulsion value is an infinite numerical value.

[0117] In the embodiment of the application, the internal space of the rotation envelope previously generated for each fan (for example, the cylinder with a radius of 73 mm and a height of 25 mm corresponding to the rear fan Fan_Rear_1) is directly marked as an absolute forbidden region in the calculation environment of path optimization. In the voxel-based path planning environment, the Cost attribute of any voxel whose center point falls within the cylinder is directly set to a special numerical value representing infinity, which means that any path search algorithm (such as A or RRT) will never expand path nodes into these voxels when expanding its search tree, thus fundamentally eliminating the possibility of the cable path crossing the fan rotation area.

[0118] In one implementation of the embodiment of the application, an equidistant outward expansion (Offset) operation is performed based on the geometric outer surface of the absolute forbidden region (cylinder with a radius of 73 mm) of Fan_Rear_1. The outward expansion distance is pre-set according to process safety specifications, which is set to 20 mm in this case. This operation generates a new coaxial cylinder with a radius of 73 mm + 20 mm = 93 mm, and the repulsion influence zone is defined as the annular space between the two coaxial cylinders. It should be noted that this 20 mm distance is a safety buffer distance evaluated based on the diameter of the cable, the fixing method (such as the size of the cable tie head), and the potential air flow disturbance effect, to ensure that even with a small deviation in cable layout, it will not enter the dangerous area.

[0119] Preferably, setting the cable bending limit constraint based on the local curvature within each window comprises:

[0120] Obtaining a material minimum bending radius of the cable to be routed;

[0121] Traversing the local curvatures of all windows based on the material minimum bending radius to screen a minimum allowable bending radius of the cable to be routed; wherein the minimum allowable bending radius is not less than its material minimum bending radius;

[0122] Searching for a region in the ideal cable path where the radius of curvature is less than the minimum allowable bending radius to obtain a routing sharp corner point;

[0123] Constructing an interference sphere with each of the routing sharp corner points as the sphere center and the minimum allowable bending radius as the sphere radius;

[0124] Performing a three-dimensional solid Boolean union operation on the interference spheres to form the cable bending limit constraint.

[0125] In an implementation manner of the embodiment of the present application, the material minimum bending radius attribute of the cable is recorded in the database. The value is determined by the cable supplier according to the industry standard IPC / WHMA-A-620 and the like through repeated bending fatigue tests on the cable sample, and is a physical limit to ensure that the internal conductor of the cable does not break and the insulating layer does not crack. The material minimum bending radius R_material_min of the cable is obtained by querying, which is 15 mm.

[0126] In an implementation manner of the embodiment of the present application, in order to facilitate manual operation and ensure long-term reliability, the bending radius during actual routing should be at least 6 times the cable diameter (3.5 mm), that is, 6*3.5 mm = 21 mm. The process requirement value (R_process_req) is compared with the material limit value (R_material_min), and the larger one of the two is taken as the final minimum allowable bending radius 21 mm.

[0127] In an implementation manner of the embodiment of the present application, a continuous coordinate point sequence {P36, P37, P38, P39, P40} is identified on the ideal cable path, and the corresponding local bending radii are {19 mm, 16 mm, 15 mm, 17 mm, 20 mm} respectively, all of which are less than 21 mm. The point with the minimum bending radius in this sequence, that is, P38 (bending radius is 15 mm), is marked as a routing sharp corner point. Multiple such routing sharp corner points can be found by traversing the entire path.

[0128] In one implementation of the embodiment, a virtual sphere is constructed in the three-dimensional space, with the three-dimensional coordinates of the found corner point P38 as the center of the sphere, and with the minimum allowed bend radius R_allow_min = 21 mm as the radius of the sphere. This sphere is referred to as the interference sphere. It is noted that the physical meaning of this interference sphere is that any path that attempts to form a bend radius smaller than R_allow_min will have its center of curvature inside this sphere. Therefore, as long as the modified path does not interfere with this sphere, it is guaranteed that the bend radius of the path in this region meets the requirement.

[0129] In another implementation of the embodiment, it is assumed that three independent interference spheres A1, A2, A3 are generated. In order to form a unified constraint region, a three-dimensional Boolean union operation is performed on these spheres: A1∪A2∪A3; the result of the operation is one or more (if the spheres do not intersect) irregular, smoothly connected geometric entities, which are determined to be the cable bend limit constraints for this cable.

[0130] Preferably, calculating the path deviation amount based on the actual wiring centerline and the digital cable guide path comprises:

[0131] calculating the maximum path deviation distance between the actual wiring centerline and the digital cable guide path, denoted as the stem path compliance deviation;

[0132] calculating the local curvatures of the actual wiring centerline and the digital cable guide path along their lengths, respectively, and then performing an integral difference operation based on the two local curvatures, denoted as the bend compliance deviation.

[0133] In the embodiment, each coordinate point P_actual_i on the actual wiring centerline is traversed. For each point, the closest point P_guide_proj on the digital cable guide path is found. Then, the Euclidean distance d_i between the two points is calculated, and after traversing all the points on the actual wiring centerline, a distance sequence {d1, d2,..., dn} is obtained. The stem path compliance deviation D_pos is defined as the maximum value in this sequence: D_pos = max(d1, d2,..., dn).

[0134] In one implementation of the embodiment, assume that in real-time evaluation of the wiring of the CPU 8-Pin cable, most of the distance values in the calculated distance sequence are between 3mm and 5mm, but at a corner close to the side wall of the case, a calculated distance value is d_peak = 9.8mm. Therefore, the main path compliance deviation D_pos of this evaluation is determined to be 9.8mm. It should be noted that using the maximum path deviation distance instead of the average distance as an indicator has higher quality control sensitivity. The average distance can mask local serious deviations (for example, the cable touches a heating element at a certain point), and the maximum deviation distance can accurately capture the most non-compliant position in the entire wiring process.

[0135] In the embodiment, in order to evaluate whether the smoothness and bending shape of the operator wiring meet the design requirements, the bending characteristics of the two paths are quantitatively compared. Specifically, the actual wiring center line and the digital cable guide path are parameterized, and their lengths are normalized to the interval [0, L], where L is the total length of the guide path. Then the local curvature functions of the two paths at any arc length s are calculated, denoted as K_actual(s) and K_guide(s) respectively. The bending compliance deviation D_curv is defined as the integral of the absolute value of the difference between the two curvature functions over the entire path length.

[0136] In another implementation of the embodiment, the mathematical expression of the integral difference operation is:

[0137] D_curv = ∫|K_actual(s) - K_guide(s)|ds);

[0138] Since the path is composed of discrete points in actual calculation, numerical integration (such as the trapezoidal rule) is used to approximate the integral value:

[0139] D_curv ≈ Σ[|K_actual_i - K_guide_i| × Δs_i];

[0140] Where: K_actual_i and K_guide_i are the average curvatures of the two paths at the i-th corresponding segment; Δs_i is the length of the segment.

[0141] Especially important is that before realizing the dynamic guidance of wiring and the management sharing of production information, the following steps are further included:

[0142] Obtain the wiring standard compliance set by the current case production batch;

[0143] According to the wiring standard compliance and the wiring compliance index, real-time single-point compliance calculation is performed to generate single-point compliance data;

[0144] mapping the single-point compliance data to the coordinate points on the digital cable guide path to generate a spatial compliance distribution map;

[0145] dividing the spatial compliance distribution map into regions by a preset compliance threshold to obtain a high-compliance region, a medium-compliance region and a low-compliance region respectively;

[0146] performing standard guide strength calculation on the medium-compliance region by the wiring standard compliance to generate a medium-region guide instruction;

[0147] performing weighted guide strength calculation on the high-compliance region and the low-compliance region respectively by the preset guide weight based on the wiring standard compliance to obtain a high-region confirmation signal and a low-region warning instruction respectively;

[0148] performing visual coding processing on the medium-region guide instruction, the high-region confirmation signal and the low-region warning instruction to realize dynamic guidance of wiring.

[0149] In the embodiment of the application, before starting a new batch of production tasks, a wiring standard compliance (SCL) is set for the batch of products through a process parameter configuration interface of a manufacturing execution system (MES), and the value is a percentage, reflecting the quality requirement level of the batch of products. For a high-reliability server batch used in a data center, the SCL is set to 95%. And for a standard PC batch used in an ordinary office environment, the SCL is set to 85%. The SCL value will be pushed to the system of each assembly station. Assuming that the current batch is a high-reliability server, SCL = 95%.

[0150] In an implementation manner of the embodiment of the application, for any point P_actual_i on the actual wiring center line, a corresponding point P_guide_proj on the digital cable guide path is found, and a local deviation d_i (including a comprehensive deviation of position and curvature) between them is calculated. The calculation model of single-point compliance C_point(i) is as follows:

[0151] C_point(i) = 1 - (d_i / D_max_local);

[0152] Wherein D_max_local is a preset maximum allowed local deviation value, for example, 5 mm. If the local deviation d_j of a point P_actual_j is calculated as 1.5 mm, then the single-point compliance C_point(j) of the point is 1 - (1.5 / 5) = 0.7, that is, 70%. The calculation is performed along the entire actual wiring path to generate a single-point compliance data sequence {C_point(1), C_point(2),...} corresponding to the path coordinate points one by one.

[0153] In an implementation manner of the embodiment of the present application, the single-point compliance is mapped from the range of 0% to 100% to a gradient color spectrum from red to green. For example, 100% compliance corresponds to pure green, 50% corresponds to yellow, and 0% corresponds to pure red. When the operator looks at the cabinet, the guide path displayed in the AR glasses is no longer a single color, but a color curve that intuitively shows the distribution of the wiring quality in space. The guide path with color coding is the spatial compliance distribution map. The spatial compliance distribution map is regionally divided by a preset compliance threshold, and high-compliance, medium-compliance and low-compliance regions are obtained respectively.

[0154] In an implementation manner of the embodiment of the present application, the calculation model of the midpoint threshold of the standard line and the full score line is: T_high=SCL+(100-SCL) / 2; T_low=SCL-(100-SCL) / 2 (the symmetric point of the standard line and the full score line about the standard line, but not lower than a certain absolute lower limit, such as 70%); for SCL=95%, T_high=95+(100-95) / 2=97.5%; T_low=95-(100-95) / 2=92.5%.

[0155] Accordingly, the path segment with single-point compliance >97.5% on the spatial compliance distribution map is divided into a high-compliance region; the path segment with compliance between 92.5% and 97.5% is divided into a medium-compliance region; and the path segment with compliance <92.5% is divided into a low-compliance region. Finally, the medium-region guide instruction, the high-region confirmation signal and the low-region warning instruction are visually coded to realize dynamic guidance of wiring.

[0156] The present application also provides an enterprise production information management sharing system based on digital twinning, which executes the enterprise production information management sharing method based on digital twinning as described above, and comprises:

[0157] A digital twinning modeling module S101 is configured to collect multi-view three-dimensional point clouds inside a product to be wired to construct a product digital twinning model;

[0158] A wiring path optimization module S102 is configured to map physical constraints in the product digital twinning model into heat dissipation duct wiring constraints and rotating area wiring constraints, and perform three-dimensional path optimization in the constraints according to cable connection start and end coordinates in the product digital twinning model, and mark the ideal cable path;

[0159] A guide path generation module S103 is configured to perform cable bending limit constraint on the ideal cable path, and correct all path segments with geometric interference in the whole duration of the ideal cable path to obtain a digital cable guide path with cable bending limit constraint;

[0160] The job information sharing module S104 is configured to collect three-dimensional data in the cabling job area to obtain an actual cabling path during the cabling job, compare the actual cabling path with the digital cable guide path, calculate and share a cabling compliance index, so as to realize dynamic guidance of cabling and management sharing of production information.

[0161] Therefore, from any viewpoint, the embodiments should be considered as being exemplary and not limiting, the scope of the application being defined by the appended claims and not by the above description, and all the changes falling within the meaning and the scope of the equivalent elements of the patent file are therefore intended to be comprised in the application.

[0162] The foregoing is considered as a specific implementation of the application, enabling those skilled in the art to understand or implement the application. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the application. Therefore, the application will not be limited to these embodiments shown herein, but will conform to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for enterprise production information management and sharing based on digital twins, characterized in that, Applied to computer chassis cabling operations, the method includes the following steps: Step S1: Collect 3D point clouds from multiple perspectives inside the product to be wired in order to construct a digital twin model of the product; Step S2: Map the physical constraints in the product digital twin model to heat dissipation airflow wiring constraints and rotation area wiring constraints, and perform three-dimensional path optimization within the constraints according to the start and end coordinates of the cable connection in the product digital twin model, and mark it as the ideal cable path. Step S3: Identify and correct the path segments where the ideal cable path and the preset cable bending limit constraint cause geometric interference, thereby generating a digital cable guide path with physical bending limit of the cable; Step S4: During the cabling operation, collect three-dimensional data in the cabling operation area to obtain the actual cabling path. By comparing the actual cabling path with the digital cable guidance path, calculate and share the cabling compliance index to achieve dynamic guidance of cabling and management and sharing of production information.

2. The enterprise production information management and sharing method based on digital twins according to claim 1, characterized in that, Before mapping the physical constraints within the product's digital twin model to heat dissipation airflow wiring constraints and rotation area wiring constraints in step S2, the following steps are also included: Extract the outline boundary coordinates of the motherboard, power supply, and hard drive bracket from the product's digital twin model, identify the coordinates of the top plane of the CPU cooler and the center point of the rear fan of the chassis, and calculate the line connecting the two points as the direction vector of the main airflow axis. Determine the safety envelope of the air duct based on the direction vector of the main air duct axis; Check the installation position and blade radius of all fans inside the chassis. Draw a circular restricted area with the fan center as the center and the blade radius + 15mm as the safety radius, and mark it as the rotating envelope.

3. The enterprise production information management and sharing method based on digital twins according to claim 2, characterized in that, Step S2 includes the following steps: Step S21: Determine the wiring constraints of the heat dissipation air duct based on the air duct safety envelope; Step S22: Determine the wiring constraints of the rotation region based on the rotation envelope; Step S23: Determine the start and end points of the cable to be laid based on the preset socket positions of the cable connectors in the product digital twin model; Step S24: Perform three-dimensional path optimization based on the starting and ending points of the cable to be laid, and then use heat dissipation airflow routing constraints and rotation area routing constraints to constrain the repulsion value, so as to lock the three-dimensional spatial trajectory with the minimum cumulative repulsion value as the path centerline. Step S25: Based on the diameter of the cable to be laid as the cross-sectional diameter, generate a three-dimensional tubular entity in the product digital twin model along the center line of the path to form an ideal cable path.

4. The enterprise production information management and sharing method based on digital twins according to claim 1, characterized in that, Step S3 includes the following steps: Step S31: Perform sliding window processing on the sequence of three-dimensional coordinate points along the ideal cable path to obtain the local curvature within each window; Step S32: Set cable bending limit constraints based on the local curvature within each window; Step S33: Select all path segments in the ideal cable path that geometrically interfere with the cable bending limit constraint as path segments to be corrected; Step S34: Between the start and end points of each path segment to be corrected, re-execute 3D path optimization to determine the corrected routing path segment; Step S35: Replace the path segment to be corrected with the corrected wiring path segment, and splice it with the path segment in the ideal cable path that does not interfere with each other, so as to obtain the digital cable guide path.

5. The enterprise production information management and sharing method based on digital twins according to claim 1, characterized in that, Step S4 includes the following steps: Step S41: Periodically and synchronously scan the real-time point cloud data of the cabling operation area; Step S42: Extract the cable skeleton line based on the real-time work area point cloud data, and use it as the actual cabling center line; Step S43: Calculate the path deviation based on the actual cabling centerline and the digital cable guide path, wherein the path deviation includes the trunk path compliance deviation and the bend compliance deviation. Step S44: Weight the routing compliance deviation and the bending compliance deviation to obtain the routing compliance index.

6. The enterprise production information management and sharing method based on digital twins according to claim 3, characterized in that, The wiring constraints for heat dissipation air ducts are determined based on the air duct safety envelope, including: A starting section scan is performed on the safety envelope of the air duct along the direction vector of the main air duct axis. Within the safety envelope of the air duct, a geometric plane with the smallest cross-sectional area is locked as the core air duct section. The three-dimensional voxelization of the air duct safety envelope is converted into an air duct voxel grid; For each voxel in the duct voxel grid, the repulsive force value is calculated based on the reciprocal of the area of ​​its cross section and the preset duct weight coefficient. Multiply the repulsion value of the voxel corresponding to the core airflow section in the airflow voxel grid by a preset penalty factor to construct the heat dissipation airflow wiring constraint.

7. The enterprise production information management and sharing method based on digital twins according to claim 3, characterized in that, Determining the routing constraints for the rotating region based on the rotating envelope includes: The internal three-dimensional space of the rotating envelope is defined as an absolutely forbidden region; Using the geometric outer surface of the absolutely forbidden zone as a reference, perform an equidistant outward expansion operation to obtain the repulsive force influence zone; A three-dimensional repulsive field is set up based on the absolute forbidden zone and the repulsive influence zone. Any point in the field located in the absolute forbidden zone is assigned the maximum repulsive force value, and any point in the repulsive influence zone is assigned a value based on its distance from the surface of the absolute forbidden zone, so as to form a wiring constraint for the rotating region; where the maximum repulsive force value is an infinite value.

8. The enterprise production information management and sharing method based on digital twins according to claim 4, characterized in that, Setting cable bending limit constraints based on the local curvature within each window includes: Obtain the minimum bending radius of the material of the cable to be laid; The minimum allowable bending radius of the cable to be laid is selected by traversing the local curvature of all windows based on the minimum bending radius of the material; wherein the minimum allowable bending radius is not less than the minimum bending radius of the material. Search for areas in the ideal cable path where the radius of curvature is less than the minimum allowable bending radius to obtain the wiring sharp corners; Construct an interference sphere with each of the wiring corner points as the center and the minimum allowable bending radius as the radius; Perform a three-dimensional solid Boolean union operation on the interference sphere to form the cable bending limit constraint.

9. The enterprise production information management and sharing method based on digital twins according to claim 5, characterized in that, The path deviation is calculated based on the actual cabling centerline and the digital cable guide path, including: Calculate the maximum path deviation distance between the actual cabling centerline and the digital cable guide path, and record it as the backbone path conformity deviation. Calculate the local curvature of the actual cabling centerline and the digital cable guide path along their lengths, and then perform an integral difference calculation based on the two local curvatures, which is recorded as the curvature compliance deviation.

10. A digital twin-based enterprise production information management and sharing system, characterized in that, For executing the enterprise production information management and sharing method based on digital twins as described in claim 1, the enterprise production information management and sharing system based on digital twins includes: The digital twin modeling module is used to collect 3D point clouds from multiple perspectives inside the product to be wired, in order to build a digital twin model of the product. The cabling path optimization module is used to map the physical constraints in the product digital twin model to the heat dissipation airflow cabling constraints and the rotation area cabling constraints, and to perform three-dimensional path optimization within the constraints based on the start and end coordinates of the cable connections in the product digital twin model, and mark the module as the ideal cable path. The guide path generation module is used to apply cable bending limit constraints to the ideal cable path and simultaneously correct all geometrically interfering path segments throughout the entire duration of the ideal cable path to obtain a digital cable guide path with cable bending constraints. The job information sharing module is used to collect three-dimensional data of the cabling operation area during the cabling operation to obtain the actual cabling path. By comparing the actual cabling path with the digital cable guidance path, it calculates and shares the cabling compliance index to realize dynamic guidance of cabling and management and sharing of production information.

Images