Product packaging iteration design system and method based on data analysis

Abstract

This invention discloses a product packaging iterative design system and method based on data analysis, relating to the field of product packaging technology. The invention acquires a simulated packaging process of the product being placed into a packaging shell, marks primary damage risk points, generates a static matching degree map in the static state after simulated packaging, applies environmental load parameters simulating the transportation environment to the static matching degree map, analyzes the relative displacement and internal stress redistribution between the product and the packaging shell, assesses their impact on the deterioration of primary damage risk and static matching degree, identifies points that exacerbate damage risk, adjusts the design scheme of the packaging shell, conducts iterative simulation analysis, and selects the packaging shell with the best overall performance as the final design scheme. This solves common quality problems in cross-regional transportation of cold chain food, such as high-altitude bag swelling and extrusion, and high-temperature softening providing space for extrusion, improving the accuracy and efficiency of packaging scheme determination.

Description

Technical Field

[0001] This invention relates to the field of product packaging technology, and in particular to a product packaging iterative design system and method based on data analysis. Background Technology

[0002] In the food supply chain, with the popularization of fresh food e-commerce, cold chain logistics and cross-regional distribution, the demand for food packaging includes not only isolating external contamination, but also maintaining the sensory quality of food and extending shelf life. Especially for perishable foods such as fresh fruits and vegetables, delicate pastries and soft meat products, their physical characteristics are usually low hardness, high elasticity or easy rheology, which requires the packaging shell to achieve extremely precise mechanical coordination with the food surface.

[0003] Currently, food packaging design mainly relies on the following traditional methods: designing the packaging cavity based on the average external dimensions of the food, allowing for a margin, and then manually filling the food for sealing and drop testing. However, this method completely ignores the softness and compressibility of the food itself. During manual filling, it is often difficult to detect microscopic over-compression, and physical trial and error not only wastes a large number of food samples but also fails to quantify the mechanical damage to the food's cellular structure caused by filling. Furthermore, there is a tendency to enlarge the cavity to avoid damaging the food, which leads to shaking during storage and transportation. For fresh fruits and vegetables, shaking can cause mechanical collisions, creating wounds that can lead to browning and microbial spoilage. For liquid or semi-solid foods, shaking can cause... Traditional methods lack the ability to quantitatively measure the balance between fit and contact pressure. They also only test the overall sealing or impact resistance of the packaging in extreme environments, while ignoring the deformation and damage mechanisms of food packaging. Food often uses modified atmosphere packaging (MAP, filled with nitrogen / carbon dioxide) or has its own respiration (fruits and vegetables produce gas through metabolism). When the cold chain breaks (temperature rises) or is transported by air / high altitude (sudden drop in external air pressure), the gas inside the packaging expands rapidly. Current technology cannot simulate how this sudden change in internal pressure alters the contact state between the food and the outer shell. This often results in contradictory situations where improved pressure resistance leads to more difficult filling, or where reserved expansion space causes food to shake and be damaged at room temperature.

[0004] To address the aforementioned problems, this invention provides a product packaging iterative design system and method based on data analysis. Summary of the Invention

[0005] To achieve the above objectives, the present invention provides the following technical solution: In a first aspect, the present invention provides a product packaging iterative design method based on data analysis, comprising the following specific steps: S1. Obtain the three-dimensional shape and physical properties of the item to be packaged, as well as the three-dimensional shape, physical properties, and hardness parameters of the packaging shell; S2. Simulate the packaging process of the item to be packaged being placed into the outer packaging shell, predict the potential damage locations of the item to be packaged and the outer packaging shell due to contact during the process, and mark them as primary damage risk points. S3. In the static state after the simulated packaging is completed, analyze the contact area, gap distribution and contact pressure distribution of the outer surface of the item to be packaged and the inner surface of the packaging shell at various spatial positions, and generate a static matching degree map. S4. Apply environmental load parameters simulating the transportation environment to the static matching degree map, introduce the internal gas expansion of the packaged object caused by changes in altitude and temperature, analyze the relative displacement and internal stress redistribution of the packaged object and the packaging shell under the action of deformation caused by gas expansion, assess its impact on the primary damage risk and the deterioration of the static matching degree, and identify the points that aggravate the damage risk. S5. Adjust the design of the packaging shell, repeat steps S2 to S4 for iterative simulation analysis, establish an evaluation model that includes damage risk and matching degree, and select the packaging shell with the best overall performance as the final design scheme by comparing the comprehensive scores.

[0006] Preferably, step S1 includes the following specific steps: S11. Extract the three-dimensional surface mesh of the object to be packaged, mark it as the surface of the object to be packaged, obtain the coordinates of the internal volume domain, surface nodes and their external normal vectors, and obtain the physical properties of the object to be packaged, including elastic modulus, Poisson's ratio, density and yield strength. S12. Extract the three-dimensional surface mesh of the packaging shell, mark it as the packaging shell surface, obtain the coordinates of the internal volume domain, surface nodes and their internal normal vectors, and obtain the physical properties and hardness parameters of the packaging shell, including the surface Brinell hardness.

[0007] Preferably, step S2 includes the following specific steps: S21. Fix the outer packaging shell, and displace the object to be packaged along the insertion path with a preset relative velocity vector to construct a penetration depth model. Calculate the shortest distance from the nodes of the object to be packaged to the surface of the outer packaging shell. The penetration depth model includes a penetration depth function, which is expressed as: ,in, This represents the penetration depth at node x at time t. This indicates the position of a node on the surface of the object to be packaged at time t. Indicates the curved surface of the packaging shell. This represents any node on the surface of the packaging shell. Indicates that the surface of the packaging shell is on The normal vector of a point, operators Dot product; S22. Substitute the elastic modulus and Poisson's ratio of the product to be packaged and the packaging shell into the equivalent contact modulus calculation formula to obtain the equivalent contact modulus. The equivalent contact modulus calculation formula is expressed as follows: ,in, Indicates the equivalent contact modulus. This represents the Poisson's ratio of the item to be packaged. This indicates the elastic modulus of the item to be packaged. The Poisson's ratio of the packaging shell. The elastic modulus of the packaging shell is represented by the equivalent contact modulus divided by the characteristic length of the contact area to obtain the equivalent contact stiffness, and the contact pressure is obtained by multiplying the equivalent contact stiffness by the penetration depth. S23. Construct the dynamic equilibrium equations, which are expressed as follows: ,in, Represents the Cauchy stress tensor, which includes normal stress. , , and shear stress , , , Indicates stress divergence, Indicates the density of the item to be packaged. Let represent the acceleration of a particle. This is discretized and solved using the finite element method to obtain the stress tensor field, which is then transformed into equivalent stress. The transformation formula is expressed as: ; S24. Obtain the structural damage index of the object to be packaged by dividing the equivalent stress of the nodes of the object to be packaged by the yield strength of the object to be packaged, and obtain the indentation damage index of the outer shell by dividing the contact pressure by the surface Brinell hardness of the outer shell. S25. During the packaging simulation process where the item to be packaged is placed into the packaging shell, obtain the maximum structural damage index of each node of the item to be packaged, and mark the node of the item to be packaged with the maximum structural damage index greater than the structural damage trigger threshold as the primary damage risk point of the item to be packaged. Obtain the maximum indentation damage index of each node of the packaging shell, and mark the node of the packaging shell with the maximum indentation damage index greater than the indentation damage trigger threshold as the primary damage risk point of the packaging shell.

[0008] Preferably, step S3 includes the following specific steps: S31. Extract the final displacement vectors of all nodes of the object to be packaged and the outer shell when the simulated packaging is completed. Superimpose the initial node coordinates with the final displacement vectors to obtain the spatial coordinates of the object to be packaged and the outer shell under static equilibrium. Update the outer normal vector of the outer surface of the object to be packaged and the inner normal vector of the inner surface of the outer shell according to the spatial coordinates under static equilibrium. S32. Construct a static gap model to calculate the gap between the item to be packaged and the outer packaging shell. The static gap model includes a static gap calculation formula, which is expressed as follows: ,in, This represents the gap distance at node x. This represents the spatial coordinates of the object to be packaged under static equilibrium. express For the nearest point of the packaging shell, This represents the inner normal vector of the updated inner surface of the packaging shell. When the gap distance is greater than 0, it is divided into a gap area. When the gap distance is less than or equal to 0, it is divided into a contact area. At this time, the equivalent static penetration depth is the negative value of the gap distance. The residual contact pressure is obtained by multiplying the equivalent contact stiffness by the equivalent static penetration depth. The bonding area is obtained by integrating the area of ​​the contact area. The global bonding rate is obtained by dividing the bonding area by the initial total surface area of ​​the packaged item. S33. For the gap area, obtain the value 1 minus the ratio of the gap distance to the tolerance threshold, compare it with 0, and take the larger value as the local matching degree of the gap area. For the contact area, obtain the absolute value of the difference between the residual contact pressure and the ideal contact pressure, divide it by the maximum allowable pressure, and then subtract this value from the value 1 to obtain the local matching degree of the contact area. S34. Based on the static surface mesh of the object to be packaged, input the gap value, contact pressure and local matching degree scalar attribute at each mesh node, and establish the mapping relationship between local matching degree and color. Attach the mapped color to the three-dimensional static surface to generate a static matching degree map.

[0009] Preferably, step S4 includes the following specific steps: S41. Obtain environmental load parameters of the transportation environment, including altitude and ambient temperature. Substitute the altitude into the air pressure calculation formula to obtain the external atmospheric pressure. The air pressure calculation formula is expressed as: ,in, Indicates external atmospheric pressure. This indicates the initial absolute pressure inside the item to be packaged. Indicates the rate of temperature lapse. Indicates altitude, This indicates the initial temperature inside the item to be packaged. Represents gravitational acceleration. Indicates the molar mass of air. The ideal gas constant is represented by the internal absolute pressure obtained by dividing the ambient temperature by the initial temperature inside the packaged item and then multiplying by the initial absolute pressure inside the packaged item. The net expansion pressure difference is obtained by subtracting the external atmospheric pressure from the internal absolute pressure. S42. On the inner surface mesh of the packaged object, the net expansion pressure difference is converted into a uniformly distributed surface force load by multiplying the net expansion pressure difference by the outer normal vector of the inner surface of the packaged object. The isotropic thermal strain tensor caused by temperature change is obtained by multiplying the environmentally sensitive deformation index of the packaged object, the difference between the ambient temperature and the initial temperature inside the packaged object, and the scalar product of the second-order unit tensor. S43. Construct the static mechanical equilibrium equations and calculate the additional displacement vector. The static mechanical equilibrium equations are expressed as follows: ,in, Represents the gradient operator, This represents the additional displacement vector caused by environmental loads. This represents the elastic stiffness tensor. Represents the thermal strain tensor. Indicates to First, calculate the gradient, then transpose the gradient result. Subtract the additional displacement vector from the gap distance and multiply it by the inner normal vector of the inner surface of the packaging shell to obtain the updated gap distance. Obtain the updated gap distance within the gap area. Divide the area where the updated gap distance is less than or equal to 0 into the new contact area. Update the equivalent static penetration depth of the contact area and the new contact area. The updated equivalent static penetration depth is the negative value of the updated gap distance. Multiply the updated gap distance by the equivalent contact stiffness to obtain the updated residual contact pressure. S44. The amount of structural damage deterioration of the packaged object is obtained by dividing the difference between the updated equivalent stress and the equivalent stress of the nodes of the packaged object by the yield strength of the packaged object. The amount of indentation damage deterioration of the packaged shell is obtained by dividing the difference between the updated residual contact pressure and the residual contact pressure by the surface Brinell hardness of the packaged shell. S45. For the updated gap area, obtain the value 1 minus the ratio of the updated gap distance to the tolerance threshold, compare it with 0, and take the larger value as the local matching degree of the updated gap area. For the updated contact area, obtain the absolute value of the difference between the updated residual contact pressure and the ideal contact pressure, divide it by the maximum allowable pressure, and then subtract this value from the value 1 to obtain the local matching degree of the updated contact area. Obtain the matching degree deterioration by subtracting the local matching degree from the updated local matching degree. S46. When the maximum structural damage index is less than or equal to the structural damage trigger threshold and the amount of structural damage deterioration of the packaged item is greater than the structural damage deterioration threshold, or the maximum indentation damage index is less than or equal to the indentation damage trigger threshold and the amount of indentation damage deterioration of the packaging shell is greater than the indentation damage deterioration threshold, or the item has been marked as a primary damage risk point and the amount of matching degree deterioration is greater than the maximum matching degree attenuation, the corresponding node will be marked as an aggravated damage risk point.

[0010] Preferably, step S5 includes the following specific steps: S51. Adjust the design scheme of the packaging shell. Based on the original packaging shell mesh, introduce a shape modification vector to obtain the coordinates of the surface nodes of the packaging shell after iteration. Adjust the physical properties and hardness parameters of the original packaging shell into a variable sequence that changes with iteration. S52. Substitute the iterated packaging shell design into steps S2 to S4 for iterative simulation, and extract the iterated static matching features, environmental degradation features, and damage risk features. The static matching features include the iterated global average matching degree, which is obtained by integrating the updated local matching degree on the initial total surface area of ​​the object to be packaged and then dividing it by the initial total surface area of ​​the object to be packaged. The environmental degradation features include the iterated matching degree degradation amount. The damage risk features include the iterated maximum structural damage index, the iterated maximum indentation damage index, and the total number of aggravated damage risk points after iteration. S53. Subtract the ratio of the maximum structural damage index after iteration to the structural damage triggering threshold from the numerical value 1, compare its magnitude with 0, and take the larger value as the structural damage evaluation score. Subtract the ratio of the maximum indentation damage index after iteration to the indentation damage triggering threshold from the numerical value 1, compare its magnitude with 0, and take the larger value as the shell indentation evaluation score. Subtract the ratio of the total number of aggravated damage risk points after iteration to the maximum tolerable risk points from the numerical value 1, compare its magnitude with 0, and take the larger value as the risk penalty evaluation score. Subtract the ratio of the amount of matching degree deterioration after iteration to the amount of maximum matching degree decay from the numerical value 1, compare its magnitude with 0, and take the larger value as the stability evaluation score. Divide the global average matching degree after iteration by the largest average matching degree among all iterative schemes to obtain the matching degree evaluation score. S54. Establish an evaluation model that includes damage risk and matching degree. The evaluation model obtains a comprehensive score by weighted summation of structural damage evaluation score, shell indentation evaluation score, risk penalty evaluation score, stability evaluation score and matching degree evaluation score, and selects the packaging shell design scheme corresponding to the highest comprehensive score as the final design scheme.

[0011] Secondly, the present invention provides a product packaging iterative design system based on data analysis, comprising: The parameter acquisition module is used to acquire the three-dimensional shape and physical properties of the item to be packaged, as well as the three-dimensional shape, physical properties, and hardness parameters of the packaging shell. The primary damage point marking module is used to simulate the packaging process of the product to be packaged being placed into the packaging shell, predict the potential damage locations of the product to be packaged and the packaging shell caused by contact during the process, and mark them as primary damage risk points. The matching degree generation module is used to analyze the contact area, gap distribution and contact pressure distribution of the outer surface of the packaged item and the inner surface of the packaging shell at various spatial locations in a static state after the simulated packaging is completed, and generate a static matching degree map. The aggravated damage point marking module is used to apply environmental load parameters simulating the transportation environment to the static matching degree map, introduce the internal gas expansion of the packaged object caused by changes in altitude and temperature, analyze the relative displacement and internal stress redistribution of the packaged object and the packaging shell under the action of deformation caused by gas expansion, assess its impact on the primary damage risk and the deterioration of the static matching degree, and identify the aggravated damage risk points. The design scheme evaluation module is used to adjust the design scheme of the packaging shell, conduct iterative simulation analysis, establish an evaluation model that includes damage risk and matching degree, and select the packaging shell with the best overall performance as the final design scheme by comparing comprehensive scores.

[0012] Thirdly, the present invention provides a storage medium comprising stored instructions, wherein, when the instructions are executed, the device in which the storage medium is located executes the product packaging iterative design method based on data analysis as described above.

[0013] Fourthly, the present invention provides an electronic device including a memory and one or more instructions, wherein one or more instructions are stored in the memory and configured to be executed by one or more processors as described above in the data analysis-based product packaging iterative design method.

[0014] Compared with existing technologies, the beneficial effects of this invention are: by dynamically simulating the filling process, it accurately captures the contact pressure distribution between the packaging shell and the food surface at the moment of assembly, enabling early identification of primary damage risk points that are highly likely to cause food surface cracking or shape collapse. By calculating the bonding area and contact pressure distribution, it identifies which parts of the interference fit have approached the yield limit of the food, leading to juice leakage or tissue necrosis. It finds the optimal geometric cavity solution between anti-shaking packaging and anti-static pressure loss. By introducing the external air pressure difference caused by altitude reduction and the gas expansion force caused by temperature difference, it simulates the internal gas exerting pressure on the packaged item (food) like a balloon. By analyzing the secondary deformation of food under this inward-outward thrust, the deterioration process of food being re-squeezed against the inner wall of the packaging shell is assessed. This effectively solves common quality problems in cross-regional transportation of cold chain food, such as high-altitude bag swelling and deterioration, and high-temperature softening providing space for compression. It accurately identifies the risk points of aggravated damage under normal conditions but deterioration under extreme temperature and pressure. By adjusting the local cavity curvature and material hardness of the packaging shell, it performs automated closed-loop iteration. Through comprehensive evaluation and scoring, it supports multiple rounds of virtual simulation improvement, significantly reducing the number of physical prototype tests, shortening the food packaging R&D cycle, reducing development costs, and improving the accuracy and efficiency of packaging solution determination. Attached Figure Description

[0015] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0016] Figure 1 A schematic diagram of the product packaging iterative design method based on data analysis provided in an embodiment of the present invention; Figure 2 This is a schematic diagram of the S2 process of the product packaging iterative design method based on data analysis provided in an embodiment of the present invention. Figure 3 This is a schematic diagram of the S4 process of the product packaging iterative design method based on data analysis provided in an embodiment of the present invention. Figure 4 A schematic diagram of the structure of a product packaging iterative design system based on data analysis provided in an embodiment of the present invention. Detailed Implementation

[0017] The technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative effort are within the scope of protection of the present invention.

[0018] In this invention, the terms “comprising,” “including,” or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, or apparatus. Without further limitation, an element defined by the phrase “comprising one…” does not exclude the presence of other identical elements in the process, method, or apparatus that includes said element.

[0019] Please see Figure 1 This invention provides a product packaging iterative design method based on data analysis, including the following specific steps: S1. Obtain the three-dimensional shape and physical properties of the item to be packaged, as well as the three-dimensional shape, physical properties, and hardness parameters of the packaging shell; In this embodiment, S1 includes the following specific steps: S11. Extract the three-dimensional surface mesh of the object to be packaged, mark it as the surface of the object to be packaged, obtain the coordinates of the internal volume domain, arbitrary nodes on the surface and their external normal vectors. The external normal vectors point to the outside of the object to be packaged. Obtain the physical properties of the object to be packaged, including elastic modulus, Poisson's ratio, density and yield strength. S12. Extract the 3D surface mesh of the packaging shell and mark it as the packaging shell surface. Obtain the coordinates of the internal volume domain, arbitrary nodes on the surface, and their internal normal vectors. The internal normal vectors point to the inside of the shell (cavity direction). Obtain the physical properties and hardness parameters of the packaging shell. The hardness parameters include the surface Brinell hardness. The physical properties of the object to be packaged and the packaging shell are obtained from the material handbook. The hardness parameters are obtained from the Brinell hardness tester. The internal volume domain is obtained from the CAD model. Use a 3D scanner to scan the object to be packaged to extract the 3D geometry. Then use finite element preprocessing software to divide the mesh and determine the space it occupies.

[0020] S2. Simulate the packaging process of the item to be packaged being placed into the outer packaging shell, predict the potential damage locations of the item to be packaged and the outer packaging shell due to contact during the process, and mark them as primary damage risk points. Please see Figure 2 In this embodiment, S2 includes the following specific steps: S21. Fix the outer packaging shell in the global coordinate system, and displace the object to be packaged along the insertion path with a preset relative velocity vector to simulate the motion state during actual insertion. The preset relative velocity vector can be set to 0.1 m / s. Construct a penetration depth model and calculate the shortest distance from the node of the object to be packaged to the surface of the outer packaging shell. The penetration depth model includes a penetration depth function, which is expressed as: ,in, This represents the penetration depth at node x at time t. This indicates that the node of the item to be packaged has come into contact with the surface of the outer packaging. This indicates that the node of the item to be packaged is not in contact with the surface of the outer packaging. This indicates the position of a node on the surface of the object to be packaged at time t. Indicates the curved surface of the packaging shell. This represents any node on the surface of the packaging shell. Indicates that the surface of the packaging shell is on The normal vector of a point takes a negative value because the surface of the packaging shell is... The normal vector of a point points outwards; to obtain the distance from the interior to the surface, we need to invert it. (Operator) This represents the dot product. This indicates traversing all points on the outer casing. Take the minimum distance; S22. Substitute the elastic modulus and Poisson's ratio of the product to be packaged and the packaging shell into the equivalent contact modulus calculation formula to obtain the equivalent contact modulus. The equivalent contact modulus calculation formula is expressed as: ,in, The equivalent contact modulus reflects the stiffness of the contact area; the smaller the equivalent contact modulus, the easier the contact area is to deform. This represents the Poisson's ratio of the item to be packaged. This indicates the elastic modulus of the item to be packaged. The Poisson's ratio of the packaging shell. The elastic modulus of the packaging shell is represented by the equivalent contact modulus, which is divided by the characteristic length of the contact area to obtain the equivalent contact stiffness. The smaller the characteristic length of the contact area, the more concentrated the contact area, and the higher the stiffness. The characteristic length of the contact area is taken as the contact radius. , Indicates the contact radius. The contact force is measured by a pressure sensor. The equivalent radius of curvature is represented by the 3D scanner. For regular geometric shapes, the radius of the object is directly measured by a 3D scanner. For irregular shapes, the geometric model of the object is obtained by 3D scanning. The radius of curvature of the contact area is then extracted by the curvature analysis tool of the CAD software. The contact pressure is obtained by multiplying the equivalent contact stiffness by the penetration depth. This pressure is used as a distributed load and acts on the contact surfaces of both objects along the opposite direction of the outer normal vector of the object to be packaged and the direction of the inner normal vector of the packaging shell. S23. Construct the dynamic equilibrium equations, which are expressed as follows: ,in, This represents the Cauchy stress tensor, a 3×3 matrix used to describe the stress state at any point inside an object, including normal stress. , , and shear stress , , , Represents stress divergence, describing the rate of spatial change of stress. Indicates the density of the item to be packaged. This represents the acceleration of a particle, obtained by differentiating the relative velocity. The discretization and solution are performed using the finite element method. The specific steps are as follows: the continuous volume domain of the object to be packaged is divided into a finite number of elements (e.g., tetrahedral elements), and the number of nodes in each element is obtained (e.g., a tetrahedral element has 4 nodes). The displacement of the nodes within an element can be expressed as: ,in, This represents the nodal displacement within the element. Let represent a linear shape function. The fundamental properties of a shape function are that it takes the value 1 at its own node and 0 at other nodes, and the sum of all shape functions at all nodes is always 1. Let represent the displacement vector of the i-th node at time t. To represent the total number of nodes within an element, establish the constitutive relation: ,in, Represents the strain tensor. , This represents the displacement gradient tensor, used to describe how displacement changes with position. Indicates transpose. express The transpose of the tensor (exchanging rows and columns). It is a symmetric tensor. The elastic stiffness tensor is represented by Young's modulus and Poisson's ratio measured through tensile tests. These values ​​are then substituted into the elastic constitutive relation for isotropic materials to construct the elastic stiffness tensor. The elastic stiffness tensor is a 6×6 matrix (because the stress / strain tensor has 6 independent components). Its elements are determined by Young's modulus and Poisson's ratio (e.g., diagonal elements reflect the relationship between normal stress and normal strain, while off-diagonal elements reflect the relationship between shear stress and shear strain). It is used to reflect the ability of food materials to resist deformation. Representing the thermal strain tensor, and Substituting into the dynamic equilibrium equation, we get Using the Galerkin method, the variational function of the shape function is applied to both sides of the dynamic equilibrium equation. Integrating, we get: , This represents the volume domain of the object to be packaged, i.e., the internal space of the object. The operator `:` represents a double dot product, which is the sum of corresponding element-wise multiplication. Represents a volume element. Represents the gradient operator, It represents a tiny change in displacement. Indicates the boundary. Describe the boundary infinitesimal element. This represents surface forces at the boundary, such as contact pressure. , It is the boundary normal vector, i.e., the surface of the object to be packaged. The unit integral is transformed into a system of dynamic algebraic equations: , The mass matrix represents the mass distribution of an object, and is composed of the element mass matrix. Assembled, that is, all units Assembled into the global quality matrix according to node number. In this example, the item to be packaged has 10 nodes, and each unit has 4 nodes. It is a 4×4 matrix. It is a 10×10 matrix. During assembly, Add the elements to The corresponding node positions (e.g., nodes 1-4 of element 1) Nodes 1-4 of element 1 and 3-6 of element 2 correspond to... (Nodes 3-6), ultimately yielding the global mass matrix. , , Represents a shape function matrix, express The transpose of the matrix, Represents the volume domain of a unit cell. Indicated in unit Integrating over the infinitesimal volume element. This represents the global nodal acceleration vector (the second derivative of displacement with respect to time). ), The stiffness matrix represents the stiffness distribution of an object and is derived from the element stiffness matrix. Assembled , The strain-displacement matrix represents the transformation of nodal displacements into strains. The elements of the strain-displacement matrix are the partial derivatives of the shape functions with respect to coordinates. The dimension depends on the element type, such as the tetrahedral element. It is a 6×12 matrix, corresponding to 6 strain components and 12 nodal displacement components. express The transpose of the matrix, Represents the global node displacement vector. This represents the global load vector. During assembly, the load vectors of all elements are... By arranging the nodes according to their numbers and adding them to the global array, the load distribution of the entire object is obtained. The Newmark method is then used to solve the system of algebraic equations. Then, the strain was calculated. ) and stress ( We obtain the stress tensor fields in the spatial and temporal domains and transform them into equivalent stresses, where the transformation formula is expressed as: ; S24. The structural damage index of the packaged object is obtained by dividing the equivalent stress of the node of the packaged object by the yield strength of the packaged object. A structural damage index greater than 1 indicates that yield damage has occurred, such as cracking. The indentation damage index of the packaging shell is obtained by dividing the contact pressure by the surface Brinell hardness of the packaging shell. An indentation damage index greater than 1 indicates that surface plastic indentation has occurred, such as scratches. During the placement process, the contact force and stress fluctuate with time, and the damage index will also change dynamically. S25. During the packaging simulation process where the product to be packaged is placed into the outer packaging shell, the maximum structural damage index of each node of the product to be packaged is obtained. The nodes of the product to be packaged with a maximum structural damage index greater than the structural damage trigger threshold are marked as primary damage risk points of the product to be packaged, indicating the location where dents, cracks or internal structural failures may occur. The structural damage trigger threshold is the critical stress value at which structural damage (such as plastic deformation or fracture) occurs, usually taken as the yield strength of the food material. The maximum indentation damage index of each node of the outer packaging shell is obtained. The nodes of the outer packaging shell with a maximum indentation damage index greater than the indentation damage trigger threshold are marked as primary damage risk points of the outer packaging shell, indicating the location where surface scratches, indentations or peeling may occur. The indentation damage trigger threshold is the critical contact pressure or penetration depth at which indentation damage (such as surface dents or material yielding) occurs, usually taken as the surface Brinell hardness or the maximum indentation depth.

[0021] S3. In the static state after the simulated packaging is completed, analyze the contact area, gap distribution and contact pressure distribution of the outer surface of the item to be packaged and the inner surface of the packaging shell at various spatial positions, and generate a static matching degree map. In this embodiment, S3 includes the following specific steps: S31. Extract the final displacement vectors of all nodes of the object to be packaged and the outer shell when the simulated packaging is completed, that is, the total displacement of the nodes during the placement process. Superimpose the initial node coordinates with the final displacement vectors to obtain the spatial coordinates of the object to be packaged and the outer shell under static equilibrium. Update the outer normal vector of the outer surface of the object to be packaged and the inner normal vector of the inner surface of the outer shell according to the spatial coordinates under static equilibrium. During the placement process, the object to be packaged and the outer shell will undergo elastic deformation and plastic residue. Only by obtaining the final geometry can the gap distribution and contact state between the two be accurately analyzed. S32. Construct a static gap model to calculate the gap between the item to be packaged and the packaging shell. The static gap model includes the static gap calculation formula, which is expressed as follows: ,in, This represents the gap distance at node x. This represents the spatial coordinates of the object to be packaged under static equilibrium. express To find the nearest point on the outer shell of the packaging, take any point on the outer surface of the product to be packaged, and find the nearest point on the inner surface of the outer shell along its normal vector direction or using a minimum distance search algorithm. This represents the inner normal vector of the updated packaging shell's inner surface. When the gap distance is greater than 0, it indicates that the two are not in contact and there is a gap, which is defined as the gap area. When the gap distance is less than or equal to 0, it indicates that the two are interfering or adhering, which is defined as the contact area. At this time, the equivalent static penetration depth is the negative value of the gap distance. The residual contact pressure is obtained by multiplying the equivalent contact stiffness by the equivalent static penetration depth. The area integration of the contact area is used to obtain the adhering area. The global adhering rate is obtained by dividing the adhering area by the initial total surface area of ​​the object to be packaged. The initial total surface area of ​​the object to be packaged is extracted by scanning the object to be packaged with a 3D scanner. After placement, the object to be packaged stops moving, the inertial force disappears, and the contact pressure is determined by static deformation (such as residual extrusion after elastic recovery). If the residual pressure exceeds the maximum allowable pressure of the material, it will lead to overpressure damage. S33. For the gap area, obtain the value 1 minus the ratio of the gap distance to the tolerance threshold. The tolerance threshold is the maximum allowable penetration depth of the contact area, representing the gap tolerance between the packaging shell and the packaged item. It is determined according to the packaging design specifications (such as the ISTA2A standard) (e.g., to avoid direct contact between the packaging and the packaged item, and to prevent friction damage). Compare its value with 0, and take the larger value as the local matching degree of the gap area. For the contact area, obtain the absolute value of the difference between the residual contact pressure and the ideal contact pressure and divide it by the maximum allowable pressure. The ideal contact pressure difference is usually taken as half of the maximum allowable pressure. The maximum allowable pressure is the maximum contact pressure that the material can withstand. The yield pressure of the material is measured by compression test (when the contact pressure exceeds the yield pressure, the material undergoes plastic deformation). Then, subtract this value from the value 1 to obtain the local matching degree of the contact area. S34. Based on the static surface mesh of the object to be packaged, input the gap value, contact pressure, and local matching degree scalar attributes at each mesh node, and establish a mapping relationship between local matching degree and color. For example, when the local matching degree is close to 1, it is mapped to the ideal fit color (such as dark green, indicating small gap and appropriate pressure). When the local matching degree decreases and is in the gap area (large gap), it is mapped to the gap tolerance color (such as a gradient from blue to red, the gradient is determined by the gap distance, the larger the gap distance, the redder the color). When the local matching degree decreases and is in the contact area (abnormal pressure), it is mapped to the overpressure risk color (such as orange-yellow, the brightness is determined by the residual contact pressure, the larger the residual contact pressure, the brighter the color). Attach the mapped colors to the three-dimensional static surface to generate a static matching degree map. In the map, cool colors (blue / red) indicate the presence of a suspended gap (excessive gap), green indicates effective packaging (small gap and appropriate pressure), and warm colors (orange-yellow) indicate the presence of overpressure (pressure exceeds the allowable value).

[0022] S4. Apply environmental load parameters simulating the transportation environment to the static matching degree map, introduce the internal gas expansion of the packaged object caused by changes in altitude and temperature, analyze the relative displacement and internal stress redistribution of the packaged object and the packaging shell under the action of deformation caused by gas expansion, assess its impact on the primary damage risk and the deterioration of the static matching degree, and identify the points that aggravate the damage risk. Please see Figure 3 In this embodiment, S4 includes the following specific steps: S41. Obtain environmental load parameters of the transportation environment, including altitude and ambient temperature. Altitude is obtained through GPS positioning, and ambient temperature is obtained through a temperature sensor. Substitute the altitude into the air pressure calculation formula to obtain the external atmospheric pressure. The air pressure calculation formula is expressed as: ,in, Indicates external atmospheric pressure. This indicates the initial absolute pressure inside the item to be packaged, typically in the form of standard atmospheres. This represents the temperature lapse rate, describing the rate at which temperature decreases for every 1 meter increase in altitude, typically taken as 0.0065 K / m. Indicates altitude, This indicates the initial temperature inside the item to be packaged, which is collected by a temperature sensor. This represents the acceleration due to gravity, taken as 9.8 m / s². Express the molar mass of air, taken as 0.029 kg / mol. The ideal gas constant is 8.314 J / (mol·K). The internal absolute pressure is obtained by dividing the ambient temperature by the initial temperature inside the package and then multiplying it by the initial absolute pressure inside the package. The effective net expansion pressure difference between the internal gas and the inner wall of the package is obtained by subtracting the external atmospheric pressure from the internal absolute pressure. The transportation environment will cause changes in the internal gas pressure. For example, as altitude increases, the external atmospheric pressure decreases, and the internal pressure increases relatively. As temperature increases, the internal pressure also increases, which in turn causes gas expansion and puts pressure on the packaging structure. S42. On the inner surface grid of the product to be packaged, the net expansion pressure difference is converted into a uniformly distributed surface force load by multiplying the net expansion pressure difference by the outer normal vector of the inner surface of the product (pointing outward, i.e., the direction of gas expansion). The isotropic thermal strain tensor caused by temperature change is obtained by multiplying the environmentally sensitive deformation index of the product, the difference between the ambient temperature and the initial temperature inside the product, and the scalar product of the second-order unit tensor. The second-order unit tensor is a 3×3 matrix, with 1s on the diagonal and 0s on the off-diagonal, representing isotropic strain. The steps for obtaining the environmentally sensitive deformation index are as follows: Select a representative food sample in the target packaging state, and test it under standard atmospheric pressure (e.g., 101.325 kPa) and... At room temperature (e.g., 20℃), the initial geometric shape of the food is measured and recorded (initial volume and surface area are obtained through 3D scanning), along with the initial moisture content, free gas content, and other physical states of the food. A baseline zero point is established. The food sample is placed in a sealed, transparent, and visible environmental simulation chamber. This test chamber can precisely programmatically control the temperature inside (simulating fluctuations from cold chain to room temperature) and connect to a vacuum pump or high-pressure gas source to precisely adjust the air pressure inside the chamber (simulating the rapid pressure drop process from sea level to high-altitude air transport). Inside the test chamber, a series of discrete environmental nodes are set to form a two-dimensional temperature-pressure matrix. For example, the temperature is set to 0℃, 10℃, 20℃, etc. At 30℃, with air pressures set at 101kPa, 80kPa, and 60kPa (corresponding to altitudes of approximately 0 meters, 2000 meters, and 4000 meters), the system automatically traverses these temperature and pressure combinations. Sufficient time is allowed at each node to ensure steady-state heat conduction and pressure equilibrium within the food. After steady-state conditions are achieved at each temperature and pressure node, a 3D optical scanner placed outside the environmental chamber scans the food in situ through high-transparency glass. Using a point cloud comparison algorithm, the real-time 3D surface area and volume of the food under that specific temperature and pressure combination are directly extracted. Compared to traditional methods like water displacement or geometric caliper measurements, this method does not contact the food and does not disrupt its original expansion state under pressure. In a free state, food may experience irreversible rupture or collapse. However, in practical applications, food is wrapped in packaging. When food is placed in a standard-shaped constrained package made of a flexible but inextensible film, the volume expands in the thermobaric chamber and has nowhere to be released. Instead, it exerts pressure on the constrained package. By measuring the tension change on the surface of the constrained package, the equivalent volume expansion thrust and anisotropic deformation distribution of the food under the constrained state can be deduced. Discrete data points (input quantities are temperature and air pressure, output quantities are volume change rate and anisotropic deformation ratio) are imported into data analysis software, and a continuous mathematical surface equation is fitted using a multivariate nonlinear regression algorithm. S43. Construct the static mechanical equilibrium equations and calculate the additional displacement vector. The static mechanical equilibrium equations are expressed as follows: ,in, Represents the gradient operator, This represents the additional displacement vector caused by environmental loads. This represents the elastic stiffness tensor. Represents the thermal strain tensor. Indicates to First, calculate the gradient, then take the transpose of the gradient calculation result. The specific steps for solving the static equilibrium equations are as follows: Divide the volume domain of the object to be packaged into a finite number of elements, obtain the number of nodes in each element, and the displacement of the nodes within the element can be expressed as: ,in, Let the displacement vector of the i-th node be represented, and the constitutive relation be established: ,at this time Using the Galerkin method, multiplying both sides of the static equilibrium equations by the variation of the shape function and integrating, we obtain: Transform the unit integral into a global system of algebraic equations: , Represents the global node displacement vector. This represents the global load vector, which is composed of the equivalent loads and boundary forces caused by thermal strain. Solving the system of algebraic equations yields The updated gap distance is obtained by subtracting the additional displacement vector dot product of the inner normal vector of the inner surface of the packaging shell from the gap distance. The updated gap distance is obtained by dividing the gap area into gap areas in step S32. The area where the updated gap distance is less than or equal to 0 is divided into the new contact area. The equivalent static penetration depth of the updated contact area and the new contact area is updated. The updated equivalent static penetration depth is the negative value of the updated gap distance. The updated residual contact pressure is obtained by multiplying the updated gap distance by the equivalent contact stiffness. S44. The amount of structural damage deterioration of the packaged object is obtained by dividing the difference between the updated equivalent stress and the equivalent stress of the node of the packaged object by the yield strength of the packaged object. The equivalent stress is updated in step S23. The amount of indentation damage deterioration of the packaged shell is obtained by dividing the difference between the updated residual contact pressure and the residual contact pressure by the surface Brinell hardness of the packaged shell. S45. For the updated gap area (the newly added contact area has been deleted), obtain the value 1 minus the ratio of the updated gap distance to the tolerance threshold, compare it with 0, and take the larger value as the local matching degree of the updated gap area. For the updated contact area (the newly added contact area has been added), obtain the absolute value of the difference between the updated residual contact pressure and the ideal contact pressure, divide it by the maximum allowable pressure, and then subtract this value from the value 1 to obtain the local matching degree of the updated contact area. Subtract the corresponding local matching degree calculated in step S32 from the updated local matching degree to obtain the matching degree deterioration. S46. When the maximum structural damage index is less than or equal to the structural damage trigger threshold and the amount of structural damage deterioration of the packaged object is greater than the structural damage deterioration threshold, or the maximum indentation damage index is less than or equal to the indentation damage trigger threshold and the amount of indentation damage deterioration of the packaging shell is greater than the indentation damage deterioration threshold, or when it has been marked as a primary damage risk point and the amount of matching degree deterioration is greater than the maximum matching degree attenuation, the corresponding node is marked as an aggravated damage risk point. The amount of matching degree deterioration is superimposed and rendered on the original static matching degree map. For example, the area with the most severe decrease in matching degree is highlighted in red, that is, the spatial location where the gap area deteriorates into a high-risk interference area. Finally, the environmental load deterioration matching degree map is output.

[0023] S5. Adjust the design of the packaging shell, repeat steps S2 to S4 for iterative simulation analysis, establish an evaluation model that includes damage risk and matching degree, and select the packaging shell with the best overall performance as the final design scheme by comparing the comprehensive scores.

[0024] In this embodiment, S5 includes the following specific steps: S51. Adjust the design scheme of the packaging shell. Based on the original packaging shell mesh, introduce a shape modification vector. The shape modification vector reflects the design adjustments such as cavity expansion and local wall thickness increase leading to inner surface normal offset. Obtain the surface node coordinates of the packaging shell after iteration. Adjust the physical properties and hardness parameters of the original packaging shell into a variable sequence that changes with iteration. S52. Substitute the iterated packaging shell design into steps S2 to S4 for iterative simulation, and extract the static matching features, environmental degradation features, and damage risk features after iteration. The static matching features include the global average matching degree after iteration. The global average matching degree is obtained by integrating the updated local matching degree on the initial total surface area of ​​the object to be packaged and then dividing it by the initial total surface area of ​​the object to be packaged. The environmental degradation features include the amount of matching degree degradation after iteration. The damage risk features include the maximum structural damage index after iteration, the maximum indentation damage index after iteration, and the total number of aggravated damage risk points after iteration. S53. Subtract the ratio of the maximum structural damage index after iteration to the structural damage trigger threshold from the numerical value of 1, compare its magnitude with 0, and take the larger value as the structural damage evaluation score. Subtract the ratio of the maximum indentation damage index after iteration to the indentation damage trigger threshold from the numerical value of 1, compare its magnitude with 0, and take the larger value as the shell indentation evaluation score. Subtract the ratio of the total number of aggravated damage risk points after iteration to the maximum tolerable risk points from the numerical value of 1, compare its magnitude with 0, and take the larger value as the risk penalty evaluation score. Subtract the ratio of the amount of matching degree deterioration after iteration to the amount of maximum matching degree decay from the numerical value of 1, compare its magnitude with 0, and take the larger value as the stability evaluation score. Divide the global average matching degree after iteration by the maximum average matching degree among all iterative schemes to obtain the matching degree evaluation score. The maximum tolerable risk points and the maximum matching degree decay are obtained through engineering settings. The maximum tolerable risk points are usually 5 in history, and the maximum matching degree decay is usually 0.1 in history. S54. Establish an evaluation model that includes damage risk and matching degree. The evaluation model obtains a comprehensive score by weighted summation of structural damage evaluation score, shell indentation evaluation score, risk penalty evaluation score, stability evaluation score and matching degree evaluation score. The packaging shell design scheme corresponding to the highest comprehensive score is selected as the final design scheme.

[0025] The steps for obtaining the above weights are as follows: obtain historical parameters of the product to be packaged and the packaging shell, as well as the historical final selected design scheme; import the obtained historical parameters of the product to be packaged and the packaging shell into each step of this embodiment to obtain the packaging shell with the best overall performance as the final design scheme; import the historical final selected design scheme and the final design scheme obtained in this step into MATLAB for fitting, and obtain the set of weight values ​​with the highest judgment accuracy.

[0026] Please see Figure 4 This invention also provides a product packaging iterative design system based on data analysis, including: The parameter acquisition module is used to acquire the three-dimensional shape and physical properties of the item to be packaged, as well as the three-dimensional shape, physical properties, and hardness parameters of the packaging shell. The primary damage point marking module is used to simulate the packaging process of the product to be packaged being placed into the packaging shell, predict the potential damage locations of the product to be packaged and the packaging shell caused by contact during the process, and mark them as primary damage risk points. The matching degree generation module is used to analyze the contact area, gap distribution and contact pressure distribution of the outer surface of the packaged item and the inner surface of the packaging shell at various spatial locations in a static state after the simulated packaging is completed, and generate a static matching degree map. The aggravated damage point marking module is used to apply environmental load parameters simulating the transportation environment to the static matching degree map, introduce the internal gas expansion of the packaged object caused by changes in altitude and temperature, analyze the relative displacement and internal stress redistribution of the packaged object and the packaging shell under the action of deformation caused by gas expansion, assess its impact on the primary damage risk and the deterioration of the static matching degree, and identify the aggravated damage risk points. The design scheme evaluation module is used to adjust the design scheme of the packaging shell, conduct iterative simulation analysis, establish an evaluation model that includes damage risk and matching degree, and select the packaging shell with the best overall performance as the final design scheme by comparing comprehensive scores.

[0027] This invention also provides a storage medium, which includes stored instructions, wherein, when the instructions are executed, the device where the storage medium is located is controlled to execute the product packaging iterative design method based on data analysis as described above.

[0028] This invention also provides an electronic device, specifically including a memory and one or more instructions, wherein one or more instructions are stored in the memory and configured to be executed by one or more processors as described above in the data analysis-based product packaging iterative design method.

[0029] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between embodiments can be referred to mutually. Each embodiment focuses on describing the differences from other embodiments. In particular, for system or system embodiments, since they are basically similar to method embodiments, the description is relatively simple. For relevant parts, refer to the description of the method embodiments. The systems and system embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without creative effort.

[0030] Those skilled in the art will further recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the present invention.

[0031] The above description of the disclosed embodiments enables those skilled in the art to make or use the invention. Various modifications to these embodiments will be readily apparent to those skilled in the art. The general principles defined in this invention may be implemented in other embodiments without departing from the spirit or scope of the invention. Therefore, the invention is not to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A product packaging iterative design method based on data analysis, characterized in that, The specific steps include the following: S1. Obtain the three-dimensional shape and physical properties of the item to be packaged, as well as the three-dimensional shape, physical properties, and hardness parameters of the packaging shell; S2. Simulate the packaging process of the item to be packaged being placed into the outer packaging shell, predict the potential damage locations of the item to be packaged and the outer packaging shell due to contact during the process, and mark them as primary damage risk points. S3. In the static state after the simulated packaging is completed, analyze the contact area, gap distribution and contact pressure distribution of the outer surface of the item to be packaged and the inner surface of the packaging shell at various spatial positions, and generate a static matching degree map. The specific steps include: S31. Extract the final displacement vectors of all nodes of the object to be packaged and the outer shell when the simulated packaging is completed. Superimpose the initial node coordinates with the final displacement vectors to obtain the spatial coordinates of the object to be packaged and the outer shell under static equilibrium. Update the outer normal vector of the outer surface of the object to be packaged and the inner normal vector of the inner surface of the outer shell according to the spatial coordinates under static equilibrium. S32. Calculate the gap between the item to be packaged and the outer packaging shell. When the gap distance is greater than 0, it is divided into a gap area. When the gap distance is less than or equal to 0, it is divided into a contact area. At this time, the equivalent static penetration depth is the negative value of the gap distance. Obtain the residual contact pressure through the equivalent contact stiffness and the equivalent static penetration depth. Integrate the area of ​​the contact area to obtain the bonding area. Obtain the global bonding rate through the bonding area and the initial total surface area of ​​the item to be packaged. S33. For the gap area, obtain the value 1 minus the ratio of the gap distance to the tolerance threshold, compare it with 0, and take the larger value as the local matching degree of the gap area. For the contact area, obtain the absolute value of the difference between the residual contact pressure and the ideal contact pressure, divide it by the maximum allowable pressure, and then subtract this value from the value 1 to obtain the local matching degree of the contact area. S34. Based on the static surface mesh of the object to be packaged, input the gap value, contact pressure and local matching degree scalar attribute at each mesh node, and establish the mapping relationship between local matching degree and color. Attach the mapped color to the three-dimensional static surface to generate a static matching degree map. S4. Apply environmental load parameters simulating the transportation environment to the static matching degree map, introduce the internal gas expansion of the packaged goods caused by changes in altitude and temperature, analyze the relative displacement and internal stress redistribution of the packaged goods and the packaging shell under the action of deformation caused by gas expansion, assess its impact on the primary damage risk and the deterioration of the static matching degree, and identify points that aggravate the damage risk; the specific steps include: S41. Obtain environmental load parameters of the transportation environment, including altitude and ambient temperature, obtain external atmospheric pressure, obtain internal absolute pressure through ambient temperature, initial temperature inside the packaged item and initial absolute pressure inside the packaged item, and obtain net expansion pressure difference through internal absolute pressure and external atmospheric pressure. S42. On the inner surface grid of the packaged item, the net expansion pressure difference is converted into a uniformly distributed surface force load. The isotropic thermal strain tensor caused by temperature change is obtained by using the environmentally sensitive deformation index of the packaged item, the ambient temperature, the initial temperature inside the packaged item, and the second-order unit tensor. S43. Calculate the additional displacement vector. Obtain the updated gap distance by using the gap distance, the additional displacement vector, and the inner normal vector of the inner surface of the packaging shell. Obtain the updated gap distance within the gap area. Divide the area where the updated gap distance is less than or equal to 0 into the new contact area. Update the equivalent static penetration depth of the contact area and the new contact area. The updated equivalent static penetration depth is the negative value of the updated gap distance. Obtain the updated residual contact pressure by using the updated gap distance and the equivalent contact stiffness. S44. Obtain the amount of structural damage deterioration of the packaged object by means of the updated equivalent stress, the equivalent stress of the nodes of the packaged object and the yield strength of the packaged object, and obtain the amount of indentation damage deterioration of the packaged shell by means of the updated residual contact pressure, the residual contact pressure and the surface Brinell hardness of the packaged shell. S45. For the updated gap area, obtain the value 1 minus the ratio of the updated gap distance to the tolerance threshold, compare the difference with 0, and take the larger value as the local matching degree of the updated gap area. For the updated contact area, obtain the local matching degree of the updated contact area, and obtain the matching degree deterioration amount through the updated local matching degree and the local matching degree. S46. When the maximum structural damage index is less than or equal to the structural damage trigger threshold and the amount of structural damage deterioration of the packaged item is greater than the structural damage deterioration threshold, or the maximum indentation damage index is less than or equal to the indentation damage trigger threshold and the amount of indentation damage deterioration of the packaging shell is greater than the indentation damage deterioration threshold, or when a primary damage risk point has been marked and the amount of matching degree deterioration is greater than the maximum matching degree attenuation, the corresponding node will be marked as an aggravated damage risk point. S5. Adjust the design of the packaging shell, repeat steps S2 to S4 for iterative simulation analysis, establish an evaluation model that includes damage risk and matching degree, and select the packaging shell with the best overall performance as the final design scheme by comparing the comprehensive scores.

2. The product packaging iterative design method based on data analysis according to claim 1, characterized in that, S2 includes the following specific steps: S21. Fix the outer packaging shell, and move the object to be packaged along the insertion path with a preset relative velocity vector, and calculate the shortest distance from the node of the object to be packaged to the surface of the outer packaging shell. S22. Obtain the equivalent contact modulus based on the elastic modulus and Poisson's ratio of the packaged item and the packaging shell. Obtain the equivalent contact stiffness through the equivalent contact modulus and the characteristic length of the contact area. Obtain the contact pressure through the equivalent contact stiffness and the penetration depth. S23. Discretize and solve the stress tensor field using the finite element method, and then convert it into equivalent stress. S24. Obtain the structural damage index of the packaged object by the equivalent stress of the nodes of the packaged object and the yield strength of the packaged object, and obtain the indentation damage index of the packaged shell by the contact pressure and the surface Brinell hardness of the packaged shell. S25. During the packaging simulation process where the item to be packaged is placed into the packaging shell, obtain the maximum structural damage index of each node of the item to be packaged, and mark the node of the item to be packaged with the maximum structural damage index greater than the structural damage trigger threshold as the primary damage risk point of the item to be packaged. Obtain the maximum indentation damage index of each node of the packaging shell, and mark the node of the packaging shell with the maximum indentation damage index greater than the indentation damage trigger threshold as the primary damage risk point of the packaging shell.

3. The product packaging iterative design method based on data analysis according to claim 2, characterized in that, S5 includes the following specific steps: S51. Adjust the design scheme of the packaging shell. Based on the original packaging shell mesh, introduce a shape modification vector to obtain the coordinates of the surface nodes of the packaging shell after iteration. Adjust the physical properties and hardness parameters of the original packaging shell into a variable sequence that changes with iteration. S52. Substitute the iterated packaging shell design into steps S2 to S4 for iterative simulation, and extract the iterated static matching features, environmental degradation features, and damage risk features. The static matching features include the iterated global average matching degree, which is obtained by integrating the updated local matching degree on the initial total surface area of ​​the object to be packaged and the initial total surface area of ​​the object to be packaged. The environmental degradation features include the iterated matching degree degradation amount. The damage risk features include the iterated maximum structural damage index, the iterated maximum indentation damage index, and the total number of aggravated damage risk points after iteration. S53. Subtract the ratio of the maximum structural damage index after iteration to the structural damage triggering threshold from the numerical value 1, compare its magnitude with 0, and take the larger value as the structural damage evaluation score. Subtract the ratio of the maximum indentation damage index after iteration to the indentation damage triggering threshold from the numerical value 1, compare its magnitude with 0, and take the larger value as the shell indentation evaluation score. Subtract the ratio of the total number of aggravated damage risk points after iteration to the maximum tolerable risk points from the numerical value 1, compare its magnitude with 0, and take the larger value as the risk penalty evaluation score. Subtract the ratio of the amount of matching degree deterioration after iteration to the amount of maximum matching degree decay from the numerical value 1, compare its magnitude with 0, and take the larger value as the stability evaluation score. Divide the global average matching degree after iteration by the largest average matching degree among all iterative schemes to obtain the matching degree evaluation score. S54. Establish an evaluation model that includes damage risk and matching degree. The evaluation model obtains a comprehensive score by weighted summation of structural damage evaluation score, shell indentation evaluation score, risk penalty evaluation score, stability evaluation score and matching degree evaluation score, and selects the packaging shell design scheme corresponding to the highest comprehensive score as the final design scheme.

4. The product packaging iterative design method based on data analysis according to claim 3, characterized in that, S1 includes the following specific steps: S11. Extract the three-dimensional surface mesh of the object to be packaged, mark it as the surface of the object to be packaged, obtain the coordinates of the internal volume domain, surface nodes and their external normal vectors, and obtain the physical properties of the object to be packaged, including elastic modulus, Poisson's ratio, density and yield strength. S12. Extract the three-dimensional surface mesh of the packaging shell, mark it as the packaging shell surface, obtain the coordinates of the internal volume domain, surface nodes and their internal normal vectors, and obtain the physical properties and hardness parameters of the packaging shell, including the surface Brinell hardness.

5. A product packaging iterative design system based on data analysis, used to implement the product packaging iterative design method based on data analysis as described in any one of claims 1-4, characterized in that, include: The parameter acquisition module is used to acquire the three-dimensional shape and physical properties of the item to be packaged, as well as the three-dimensional shape, physical properties, and hardness parameters of the packaging shell. The primary damage point marking module is used to simulate the packaging process of the product to be packaged being placed into the packaging shell, predict the potential damage locations of the product to be packaged and the packaging shell caused by contact during the process, and mark them as primary damage risk points. The matching degree generation module is used to analyze the contact area, gap distribution and contact pressure distribution of the outer surface of the packaged item and the inner surface of the packaging shell at various spatial locations in a static state after the simulated packaging is completed, and generate a static matching degree map. The aggravated damage point marking module is used to apply environmental load parameters simulating the transportation environment to the static matching degree map, introduce the internal gas expansion of the packaged object caused by changes in altitude and temperature, analyze the relative displacement and internal stress redistribution of the packaged object and the packaging shell under the action of deformation caused by gas expansion, assess its impact on the primary damage risk and the deterioration of the static matching degree, and identify the aggravated damage risk points. The design scheme evaluation module is used to adjust the design scheme of the packaging shell, conduct iterative simulation analysis, establish an evaluation model that includes damage risk and matching degree, and select the packaging shell with the best overall performance as the final design scheme by comparing comprehensive scores.

6. A storage medium, characterized in that, The storage medium includes stored instructions, wherein, when the instructions are executed, the device containing the storage medium is controlled to perform the data analysis-based iterative design method for product packaging as described in any one of claims 1-4.

7. An electronic device, characterized in that, It includes a memory, and one or more instructions, wherein one or more instructions are stored in the memory and configured to be executed by one or more processors as described in any one of claims 1-4.

Images