受力分析提升包装结构稳定性研究

摘 要 为解决塑料食品包装的结构稳定性不足问题，本文从力学角度探究包装结构与受压力时包装最大形变之间的关系，对拱形结构进行受力分析，并使用ANSYS有限元分析系统探究施加相同压力下不同高斯曲率与不同结构的包装的最大形变，设计出可将最大形变减少12.4%的新型拱形包装结构，从理论与数值计算两个方面验证了拱形结构对提升包装结构稳定性的可行性。本实验结果将物理理论与实际应用结合，设计出更稳定的食品包装结构，为包装快递行业提供技术支持与理论参考，以减少因食品过度包装带来的环境污染。

关键词 受力分析;包装结构;高斯曲率;稳定性

固体材料在不同结构形态下的应力响应不同，了解结构与应力之间的关系，是解决建筑、产品、包装等领域设计问题的基础。各种固体和凝聚体材料的结构与应力关系研究，已成为凝聚态物理的一个重要研究方向。随着计算机技术的发展，利用计算方法进行实验设计与数值模拟来解决各种复杂结构形态的应力问题，已成为相关研究的一个核心领域。产品包装是市场广阔的应用领域，但目前国内外较少从结构、受力等方面研究产品包装的稳定性问题，市面上的塑料包装产品大部分采用方形或圆柱形设计，结构单一，不能满足特定的防渗漏、防挤压等目标要求。新的产品设计应具有特定结构形态，且以相关材料的结构形态与应力分析为基础。本项目利用ANSYS有限元分析系统，对塑料材料的拱形结构在不同结构参数下的应力响应进行了细致研究与分析，并在此基础上提出了新的包装产品设计建议方案。

1 包装模型的受力分析

1.1 外卖包装模型的受力分析

本研究中将外卖包装内的液体理想化为静止液体，分析液体作用在平面上的总压力大小。已知重力场中静止压强分布为（取垂直坐标与重力加速度方向一致）

p =p0 +ρgz （1）

其中，p0 为海平面z=0上的大气压强。

柱体容器内水面上部压力分布均匀

p×A =p0×A + ρgz ×A （2）

水面下部压力为p0，因此由液体产生的力为

Fp = ρgz ×A （3）

即作用在水下受力斜面A 上的静水作用力等于以该面积为底面，以该面积面心的深度为高的柱体中的水重，方向垂直于斜面。

对截面为如图1所示倒圆台形的容器，容器形心位置为

对斜面，总压力的作用点为

其中，yc 为平面形心的y 坐标，Icx 为面积对通过形心且平行于ox 轴的轴线的惯性矩。由式（5）可知，压力的中心总位于形心的下方，随着淹没深度的增加，压力中心逐渐靠近形心。本实验中假设液体盛满整个容器，其压力中心即为形心，在分析容器受力时，仅需分析形心位置受力即可。

1.2 塑料包装容器结构设计

设计塑料包装容器时，应保证其力学性能、强度性能、密封性能与稳定性测试均达到合格标准。传统包装结构的侧面是曲率为0的斜面结构，斜面结构较易变形，稳定性较差。

拱形结构的力学特点是能够承载较大的压力。拱的结构由拱券及支座组成，在荷载作用下主要承受轴向压力，有时也承受弯矩。拱的性质使应力可以较为均匀地分布，避免集中在中央导致中央形变过大（见图2）。在受压状态下能把受到的压力分解为向下的压力和向外的推力，受力点将力传递给相邻的区域，外力沿瓶身曲面分散开，使拱形的每个部位都分担一部分力，能有效起到抗挤压的作用。

2 模型设计

本文使用ANSYS有限元软件来建模和模拟计算。在模型设计中，考虑网格划分与曲率改变的难易程度，采用正六面体（见图3）模型近似模拟塑料包装的圆柱体结构，当立方体的边的数目足够多时，可认为立方体趋近于平滑圆柱体，使用正六面体进行静力学结构受力的有限元分析有理论依据与可操作性。在改变模型结构时，调整正六面体每个面的高斯曲率并施加压力，通过观察正六面体模型的最大形变来探究包装设计的拱形结构对应力的响应，验证拱形结构对稳定性的改进情况。

3 研究结果与讨论

通过市场调研，将有限元分析系统中，正六面体模型的材料设置为市面上塑料餐盒常用的材料，即塑性PP材料，查阅市场资料，得到材料各参数如表1所示。

采用高为0.02m，底边边长为0.01m 的PP材料的平行六面体作为包装模型，建立坐标系，坐标系中心位于包装底面正六面体中心位置，z 轴指向包装顶面。仅考虑弹性形变，设包装材料均匀。在有限元分析中，改变包装模型受到的压力，施力点位于包装顶部中心位置，压力方向竖直向下，即沿-z 轴方向，改变压力的数值大小。考虑到塑料包装的塑性形变更易发生在最大形变处，以不同压力下包装模型的最大形变为标准，判断包装的稳定性。将模型进行网格划分后，将模型底面设置为固定支撑，设模型受到的标准重力加速度为9.8066m/s2，进行静力学分析。

3.1 包装侧面曲率为0时

改变包装受到的压力，探究不同压力下包装的最大形变情况如图4所示。

图4说明，仅考虑弹性形变，包装侧面曲率为0时，包装的最大形变与受到的压力正相关，最大形变与压力关系的直线拟合曲线的比例系数为0.2984，当压力增加1N时，最大形变增加0.2984×10-7m。

3.2 改变包装侧面为高斯曲率时

在有限元分析中采用高斯曲率作为曲率的度量，在微分几何中，曲面上一点的高斯曲率即该点主曲率κ1 与κ2 的乘积。高斯曲率的值仅依赖于曲面上的距离如何测量，更能反映曲面在两个方向上的弯曲情况。分别对包装模型的顶部相同位置施加不同大小的压力，记录改变包装侧面曲率（如图5所示）时，相同压力下包装的最大形变的变化。

由图6可知，高斯曲率越大，相同压力下的包装的最大形变越小，包装的最大形变随压力线性变化。控制模型受到的压力大小与压力的受力点不变，改变模型侧面的高斯曲率。当压力为4N时，研究不同曲率包装模型的最大形变变化。

如图7所示，对数据进行对数函数拟合后发现，增加包装侧面的高斯曲率可以减小相同压力下包装的最大形变，提升包装的稳定性。压力为4N，当侧面高斯曲率为0.11/mm2 时，包装的侧面最大形变减少了10%，当侧面高斯曲率为0.21/mm2 时，包装侧面最大形变减少了11.5%。由数据拟合发现，包装侧面的高斯曲率与受到压力下包装最大形变的关系呈对数关系，即随着包装侧面的高斯曲率的增加，最大形变的减小速度变慢，且高斯曲率越大，包装的开模难度越大，因此不能通过一直增加包装的高斯曲率提升包装的稳定性。

3.3 侧面两个方向的曲率的配合

当仅改变包装的每个面的高斯曲率时，包装的稳定性提升有限，随着每个面的高斯曲率的增加，包装的开模难度增大，为了进一步提升包装的稳定性，结合拱形结构的优势，在已有的曲率下增加小面积的大曲率部分作为包装加强筋。经过图7分析，认为在包装曲率为0.21/mm2 时，包装的最大形变较小，且开模难度较小，因此以该曲率模型为基础，增加小面积的曲率，分析模型受到压力时的最大形变。

在原有曲率上进一步改变曲率部分面积的曲率，得到结果如表2所示。当小面积拱形曲率为0.221/mm2 时，同一压力下模型的最大形变相较单一拱形模型减小了1.03%。且在已有拱形上增加拱形结构时，后加的拱形结构的高斯曲率的大小对减小包装最大形变的影响不大。

综合以上有限元分析结果可知，将塑料包装的侧面结构改为有曲率的拱形结构可以有效减少受到压力时的包装最大形变，从而提升模型的稳定性，在拱形结构的基础上增加进一步的拱形结构可以进一步提升模型的稳定性，但效果不大。

4 模型设计

使用3D建模软件Rhino犀牛，将通过理论分析设计的外卖包装模型可视化，得到的模型结构图如图8所示。将侧面改进为拱形结构，为了提升模型的稳定性，进一步将模型的底面进行改进，使用了较为稳定的花瓣形结构作为包装底部。

设计的模型在主体上采用碗身加盖的结构。其中模型侧壁采用对称的拱形设计，整个包装的拱形结构高斯曲率为0.21/mm2，凸起部分高斯曲率为0.221/mm2，经过有限元分析，该结构可将最大形变减少12.4%。为了保证包装的稳定88性，底部采用花瓣形结构，通过凹陷设计和拉长筋的结构把容器内的压力在容器底部均衡分散。为方便运输，模型顶部为加装有把手的盖子，其中盖子底部设计有简单的螺纹，更好地增加了外卖盒的密闭性。盖子上加装的简易把手能够使得外卖包装在运输过程中更加便于携带，减少塑料袋在运输过程中的使用。盒子的四壁与盒子的底部转折处设计成了大曲率半径过度，防止因转折处半径过小而使得吹塑厚度不足与应力集中导致容器受挤压与跌落时的凹陷与破裂。

5 结论

本项目为北京师范大学珠海校区2022年大学生创新创业训练计划校级推选优秀项目（已答辩结项）。本文关注日常生活中常使用的塑料食品包装的结构稳定性问题，利用ANSYS有限元分析系统对塑料包装模型进行了结构改进，探究了不同曲率下包装受到压力后的最大形变的变化。通过有限元分析得出拱形结构可以提升包装的稳定性，随着高斯曲率的增大，包装的最大形变减小。本文创新性地在单一拱形结构的基础上进一步增加小面积的大曲率拱形结构，经过有限元分析发现该结构对提升包装稳定性有一定效果，且在实际生产过程中开模难度较小，最终得到当包装的拱形结构高斯曲率为0.21/mm2，凸起部分高斯曲率为0.221/mm2 时，可将最大形变减少12.4%，从结构上提升了塑料食品包装的稳定性，从源头上减少因包装稳定性不足导致的过度包装以提升稳定性造成的环境污染。改进的拱形塑料包装有助于节省塑料包装的原料，减少塑料污染。

本项目利用ANSYS有限元分析系统，对塑料材料的拱形结构在不同结构参数下的应力响应进行了细致研究与分析，将物理理论与实际应用结合，改进的塑料包装为包装快递行业提供技术支持与理论参考，减少食品过度包装带来的环境污染，提升人民生活质量。

参 考 文 献

[1] 熊跃龙，叶海明，杜文康，等. 物态方程曲面上临界点高斯曲率的物理意义[J]. 物理与工程，2022，32（4）：5-7，23.

XIONG Y L， YE H M， DU W K， et al. Physical significanceof Gaussian curvature at critical points of equation ofstate surfaces[J]. Physics and Engineering， 2022， 32（4）：5-7， 23. （in Chinese）

[2] 刘利. 拱结构稳定研究的传递矩阵法[D]. 长沙：湖南大学，2019.

[3] 薛璟，张晓霞. 餐饮外卖包装盒创新设计研究[J]. 艺术科技， 2019， 32（1）： 184-185.

XUE J， ZHANG X X. Research on Innovative Design ofFood Takeaway Box[J]. Art Science and Technology，2019， 32（1）： 184-185. （in Chinese）

[4] 王安霞，侯贺维. 以用户体验为中心的外卖食品包装创新设计研究[J]. 包装工程， 2020， 41（10）：186-192.

WANG A X， HOU H W. Research on innovative design oftake-away food packaging based on user experience[J].Packaging Engineering， 2020， 41（10）： 186-192. （in Chinese）

[5] 王超颖. 压拱型预应力型钢混凝土框架试验研究及有限元分析[D].济南：济南大学， 2019

[6] 严国平，彭震奥，钟飞，等. 精密包装弹性支撑缓冲系统振动特性修正模型及分析[J]. 振动与冲击， 2021， 40（22）： 251-258.

YAN G P， PENG Z A， ZHONG F， et al. Precision packagingcorrection model and elastic buffer system vibrationcharacteristics analysis[J]. Journal of Vibration and Shock，2021， 40（22）： 251-258. （in Chinese）

[7] 苏潇阳，朱国敬，康厚军. 具有弹性转动约束的斜拉桥多索浅拱模型及其面内自由振动分析[J]. 地震工程与工程振动，2021， 41（3）：75-85.

SU X Y， ZHU G J， KANG H J. Multi-cable-shallow archmodel and in-plane free vibration analysis of cable-stayed bridgewith elastic rotation constraint[J]. Earthquake Engineering andEngineering Vibration， 2021， 41 （3）： 75-85. （in Chinese）

[8] 雷家艳，崔怡观，周志成，等.刚性吊杆拱桥主要受力构件冲击效应差异性[J].厦门大学学报（自然科学版），2022，61（4）：600-607.

LEI J Y， CUI Y G， ZHOU Z C， et al. Difference of impacteffects of Main stressed members of rigid derrick arch bridge[J]. Journal of Xiamen University （Natural Science Edition），2022， 61（4）： 600-607. （in Chinese）

[9] 胡雅芬，赵增辉，付靖宇.从一道力学竞赛题解析扭曲变形效应对梁应力分布的影响[J].力学与实践，2023，45（1）：200-205.

HU Y F， ZHAO Z H， FU J Y. Analysis of the influence ofdistortion effect on stress distribution of beam from a mechanicalcompetition problem[J]. Mechanics and Practice，2023，45（1）： 200-205. （in Chinese）

[10] 张伟为，康元顺，崔哲华，等.基于正交试验方法的大型有面外支撑杆X撑结构的屈曲分析和优化设计[J].工程力学，2022，39（S1）：261-271.

ZHANG W W， KANG Y S， CUI Z H， et al. Bucklinganalysis and optimization design of large X-braced structurewith out-of-plane support based on orthogonal test method[J]. Engineering Mechanics， 2022，39（S1）：261-271. （in Chinese）

基金项目： 北京师范大学校级大学生创新创业训练计划项目（X202119027057）。