CN-122021197-A - Buffer packaging structure optimization design method and system based on finite element simulation

CN122021197ACN 122021197 ACN122021197 ACN 122021197ACN-122021197-A

Abstract

The invention discloses a buffer packaging structure optimization design method and a buffer packaging structure optimization design system based on finite element simulation, which relate to the technical field of buffer packaging finite element analysis, and the method is used for acquiring preset phase change characteristic data, buffer energy absorption data and simulated working condition boundary data aiming at a transport packaging scene with phase change characteristic contents to construct a fluid-solid coupling finite element model; and obtaining an equivalent stress coefficient envelope value through time-course envelope analysis, correcting the buffering and energy-absorbing data according to the difference between the envelope value and a preset constraint threshold, iterating until convergence conditions are met, and outputting target buffering and energy-absorbing data. The invention can accurately adapt to the phase change characteristics of the content, and effectively improve the protection reliability and the simulation optimization efficiency of the buffer package in the whole transportation period.

Inventors

DONG TIANYI
WANG BAOXIA
Zhu mou
HONG ZETAO

Assignees

安徽农业大学

Dates

Publication Date: 20260512
Application Date: 20260413

Claims (10)

1. The buffer packaging structure optimization design method based on finite element simulation is characterized by comprising the following steps of: s1, acquiring preset phase change characteristic data of the content of a target packaging container, buffering and energy absorbing data of the target packaging container and simulated working condition boundary data of a target transportation line corresponding to the target packaging container; S2, constructing a fluid-solid coupling finite element model according to preset phase change characteristic data and simulated working condition boundary data, and inputting buffer energy absorption data into the fluid-solid coupling finite element model; S3, outputting a first equivalent stress time course curve of a first target area of the target packaging container and a second equivalent stress time course curve of a second target area of the target packaging container in a preset time window through a fluid-solid coupling finite element model, and carrying out time course envelope analysis according to the first equivalent stress time course curve and the second equivalent stress time course curve to obtain a plurality of equivalent stress coefficient envelope values; s4, if any equivalent stress coefficient envelope value exceeds a preset constraint threshold, correcting the buffer energy absorption data at the next preset sampling moment according to the difference value of the two coefficient envelope values, and substituting the corrected buffer energy absorption data into the fluid-solid coupling finite element model; S5, repeatedly executing S3 to S4, and ending the repeated execution if the equivalent stress coefficient envelope values corresponding to all preset sampling moments do not exceed a preset constraint threshold value and the difference value of the buffering energy absorption data corresponding to the two adjacent execution processes does not exceed a preset difference threshold value, and taking the current buffering energy absorption data as target buffering energy absorption data; S6, outputting target buffering energy absorption data.
2. The optimization design method of the buffer packaging structure based on finite element simulation according to claim 1, wherein the steps S3 to S4 are repeatedly executed, and the repeated execution times are accumulated at the same time, and if the accumulated repeated execution times reach a preset maximum execution times and the equivalent stress coefficient envelope value corresponding to the preset sampling time still exceeds a preset constraint threshold, or the difference values of the buffer energy absorption data corresponding to all adjacent preset sampling times exceed a preset difference threshold, the iteration is stopped; And extracting the buffer energy absorption data corresponding to the equivalent stress coefficient envelope value with the minimum absolute difference value with the preset constraint threshold value from the equivalent stress coefficient envelope values corresponding to all preset sampling moments, taking the buffer energy absorption data as target buffer energy absorption data, and entering S6.
3. The buffer packaging structure optimization design method based on finite element simulation, which is disclosed by claim 1, is characterized in that the preset phase change characteristic data is obtained through a preset physical property test and comprises physical property parameters changing in real time along with the phase change process, the phase change process ratio is the ratio of the preset phase change rate of the content to the preset total phase change capacity, the simulated working condition boundary data comprises real-time environment temperature time course data and impact vibration load spectrum of a target transportation line, the buffer energy absorption data is a preset buffer structure unit volume energy absorption threshold, and the physical property parameters changing in real time along with the phase change process ratio at least comprise elastic modulus.
4. The method for optimizing the design of the buffer packaging structure based on finite element simulation according to claim 3, wherein the step of calculating the equivalent stress coefficient envelope value is specifically as follows: the first target area and the second target area comprise a plurality of finite element grid nodes, and a first equivalent stress time-course curve of the finite element grid nodes of the first target area and a second equivalent stress time-course curve of the finite element grid nodes of the second target area are output through a fluid-solid coupling finite element model; For each preset sampling time in a preset time window, carrying out time-course envelope analysis on a first equivalent stress time-course curve corresponding to the preset sampling time, and extracting an envelope peak of the time-course curve as a first stress peak of a plurality of finite element grid nodes of the preset sampling time; Performing time-course envelope analysis on a second equivalent stress time-course curve corresponding to the preset sampling moment, and extracting an envelope peak value of the time-course curve to serve as a second stress peak value of a plurality of finite element grid nodes at the preset sampling moment; Marking common finite element grid nodes under a first target area and a second target area as intersection finite element grid nodes, and taking the ratio of a second stress peak value to a first stress peak value of the same intersection finite element grid node as an equivalent stress coefficient of the intersection finite element grid nodes; And performing envelope fitting by taking a plurality of intersection finite element grid nodes at the preset sampling moment as independent variables to obtain an equivalent stress coefficient envelope curve, and extracting a maximum value in the equivalent stress coefficient envelope curve as an equivalent stress coefficient envelope value at the preset sampling moment.
5. The method for optimizing the design of the buffer packaging structure based on finite element simulation according to claim 4, wherein the step of correcting the buffer energy absorption data is specifically as follows: Extracting equivalent stress coefficient envelope values of the buffering energy absorption data corresponding to preset sampling moments, acquiring preset constraint thresholds, and calculating relative difference values of the buffering energy absorption data and the preset constraint thresholds, wherein a relative difference value calculation formula is as follows: , wherein, The relative difference value is indicated as such, Represents the envelope value of the equivalent stress coefficient, Representing a preset constraint threshold; Simultaneously extracting the phase change process duty ratio corresponding to the preset sampling moment; Acquiring a preset phase change progress matching correction function, wherein the value of the preset phase change progress matching correction function is obtained by weighting a phase change basic value and an elasticity correction value, the phase change basic value is obtained according to a preset mapping rule according to the current phase change progress ratio, and the elasticity correction value is obtained according to the inverse proportion of the elasticity modulus of the content of the target packaging container; Calculating buffering energy absorption data corresponding to the next preset sampling time according to a nonlinear correction formula, wherein the specific nonlinear correction formula is as follows: , wherein, Represents the buffered energy absorption data corresponding to the next preset sampling moment, Represents the buffered energy-absorbing data corresponding to the current preset sampling moment, Representing a preset linear correction basis coefficient, Representing a preset nonlinear correction strengthening coefficient, And representing that the preset phase change process is matched with the correction function.
6. The method for optimizing the design of the cushioning packaging structure based on finite element simulation according to claim 5, wherein the step S5 further comprises: Extracting the maximum value of the equivalent stress coefficient envelope values obtained in the steps S3 to S4 of the two adjacent repeated execution steps, and calculating the relative change rate of the equivalent stress coefficient envelope values; calculating a time window adjustment coefficient and a sampling density adjustment coefficient according to a time window adjustment formula, wherein the time window adjustment formula is as follows: , , wherein, Representing the time window adjustment factor(s), Representing the relative rate of change of the equivalent stress factor envelope value, Indicating a preset time window adjustment sensitivity factor, Representing the sampling density adjustment coefficient, Adjusting a sensitivity coefficient for a preset sampling density; Taking the product of the time window adjusting coefficient and the length of the current time window as the updated time window length, and taking the quotient of the preset sampling time interval and the sampling density adjusting coefficient as the updated preset sampling time interval; substituting the updated time window length and the updated interval of the preset sampling time into the step S3 repeatedly executed in the next round.
7. The method for optimizing the design of the cushioning package structure based on finite element simulation according to claim 6, wherein the step S5 further comprises: if any abrupt change amplitude of impact vibration load spectrum or real-time environment temperature time course data in the simulated working condition boundary data exceeds a preset abrupt change threshold value in the adjacent execution process, immediately assigning a time window adjustment coefficient to be a preset minimum window coefficient, assigning a sampling density adjustment coefficient to be a preset maximum density coefficient, and re-executing the steps S3 to S4 after updating the intervals of the preset time window length and the preset sampling time.
8. The method for optimizing the design of the cushioning packaging structure based on finite element simulation according to claim 7, wherein the step S5 further comprises: If the calculated residual error of the fluid-solid coupled finite element model exceeds a preset residual error threshold value and the difference value of the length of the time window before and after updating exceeds a preset window difference value threshold value, obtaining a difference value of the calculated residual error exceeding the preset residual error threshold value and a difference value of the length of the time window before and after updating exceeding the preset window difference value threshold value, and obtaining a grid adjustment coefficient through weighted summation; Obtaining the finite element grid intervals of the intersection finite element grid nodes, reducing the finite element grid intervals according to grid adjustment coefficients, obtaining new finite element grid intervals, updating the finite element grid intervals of the intersection finite element grid nodes, obtaining new intersection finite element grid nodes, recalculating to obtain new equivalent stress coefficient envelope values, and then re-executing the steps S3 to S4.
9. Buffer packaging structure optimal design system based on finite element simulation, which is characterized by comprising: The data acquisition module is used for acquiring preset phase change characteristic data of the content of the target packaging container, buffering energy absorption data of the target packaging container and simulated working condition boundary data of the target packaging container corresponding to the target transportation line; the finite element model construction module is used for constructing a fluid-solid coupling finite element model according to preset phase change characteristic data and simulated working condition boundary data, and replacing the buffering energy absorption data with the fluid-solid coupling finite element model; The time-course envelope analysis module is used for outputting a first equivalent stress time-course curve of a first target area of the target packaging container and a second equivalent stress time-course curve of a second target area of the target packaging container in a preset time window through the fluid-solid coupling finite element model, and carrying out time-course envelope analysis according to the first equivalent stress time-course curve and the second equivalent stress time-course curve to obtain a plurality of equivalent stress coefficient envelope values; the buffering and energy-absorbing correction module is used for correcting buffering and energy-absorbing data at the next preset sampling moment according to the difference value of any equivalent effect coefficient envelope value and the preset constraint threshold value if the value exceeds the preset constraint threshold value, and substituting the corrected buffering and energy-absorbing data into the fluid-solid coupling finite element model; the iteration module is used for repeatedly executing the time interval envelope analysis module to the buffering and energy absorbing correction module, and if the equivalent stress coefficient envelope values corresponding to all preset sampling moments do not exceed a preset constraint threshold value and the difference value of the buffering and energy absorbing data corresponding to the adjacent two execution processes does not exceed a preset difference threshold value, the repeated execution is ended, and the buffering and energy absorbing data corresponding to the current iteration period is determined to be target buffering and energy absorbing data; and the output module is used for outputting target buffering energy absorption data.
10. A computer readable storage medium storing a computer program, which when executed by a processor performs the method according to any one of claims 1-8.

Description

Buffer packaging structure optimization design method and system based on finite element simulation Technical Field The application relates to the technical field of finite element analysis of buffer packaging, in particular to a buffer packaging structure optimization design method and system based on finite element simulation. Background In the long-distance transportation packaging field with the phase-change characteristic contents, the integrity and quality stability of the contents in the transportation process are directly determined by the protective performance of the buffer packaging structure, and the design core of the buffer packaging is to accurately match the mechanical load characteristic in the transportation process so as to provide continuous and effective protection for the contents. The existing buffer package simulation design technology mostly builds a simulation analysis model based on physical parameters in a fixed state of the content, does not fully consider a continuous phase change process of the content caused by environmental temperature fluctuation in a transportation process, cannot accurately represent a continuous dynamic change rule of the physical parameters of the content in the phase change process, and cannot accurately restore a continuously evolving fluid-solid coupling mechanical action relationship between the content and a package structure in the phase change process. The static design mode based on the fixed physical parameters cannot effectively capture the most unfavorable working condition of the stress of the buffer structure in the whole transportation period, the condition that the local stress exceeds a protection threshold value in the phase change process is easy to occur, and meanwhile, the reasonable balance of the buffer protection performance and the light weight of the package is difficult to realize, or the unnecessary consumption of the packaging material is caused by the protection redundancy, or the damage of the content occurs in the transportation process due to the insufficient protection capability. To sum up, the core defect of the prior art is that in the design scene of the long-distance transportation buffer package with the content with the phase change characteristic, the dynamic physical property change of the continuous phase change of the content cannot be adapted, so that the reliability of the full-period protection of the buffer package is insufficient. Disclosure of Invention In order to solve the technical problems, the buffer packaging structure optimization design method and system based on finite element simulation are provided, and the technical scheme solves the problems in the background technology. According to the first aspect, the embodiment of the application provides an optimization design method of a buffer packaging structure based on finite element simulation, which comprises the following steps of S1, obtaining preset phase change characteristic data of contents of a target packaging container, buffer energy absorption data of the target packaging container and simulated working condition boundary data corresponding to a target transportation line of the target packaging container, S2, constructing a fluid-solid coupling finite element model according to the preset phase change characteristic data and the simulated working condition boundary data, inputting the buffer energy absorption data into the fluid-solid coupling finite element model, S3, outputting a first equivalent stress time course curve of a first target area of the target packaging container and a second equivalent stress time course curve of a second target area of the target packaging container through the fluid-solid coupling finite element model in a preset time window, performing time course envelope analysis according to the first equivalent stress time course curve and the second equivalent stress time course curve, obtaining a plurality of equivalent stress coefficient envelope values, S4, correcting the buffer energy absorption data at the next preset sampling time according to the difference value of the two equivalent stress coefficient envelope values, replacing the buffer energy absorption data after correction with the buffer energy absorption data at the next preset sampling time with the fixed constraint threshold value, S5, repeatedly executing S3 to S4, and repeatedly executing the buffer energy absorption data corresponding to the preset stress coefficient envelope values at the preset time coefficient envelope values which are not exceeding the preset constraint threshold value, and repeatedly executing the buffer coefficient data when the buffer coefficient envelope values are not corresponding to the preset buffer coefficient of the preset value is not corresponding to the preset, and executing the buffer data is not corresponding to the preset and the preset when the difference value is executed. In a second aspect, the embo