纸板包装的模拟是业界一个有趣且重要的课题。然而,对此类材料的行为进行建模并非易事。它需要物理测试数据和合适的材料模型。迄今为止,标准方法涉及使用希尔势函数定义正交各向异性线弹性和塑性。许多在 Abaqus 环境中描述这种方法的文章都是与波兹南生命科学大学和波兹南理工大学合作开发的。该研究团队的一篇值得关注的论文是 T. Garbowski、T. Gajewski 和 A. Knitter-Piątkowska 撰写的 《模拟和数字折痕线对瓦楞纸板和包装力学参数的影响》 。
在最新版本的 Abaqus 中,引入了一种新的材料模型,这可能会显著改变纸板包装模拟的方法。这就是纸板塑性模型,专门用于模拟纸板的永久变形。它源自剑桥麻省理工学院 QS Xia 的博士论文 《纸板非弹性变形和分层力学》 。
纸板塑性模型适用于 Abaqus/Standard 和 Abaqus/Explicit 求解器(用于平面应力元素和 3D 连续体元素),并考虑纸板在三个主要方向上的行为:
ZD – 贯穿厚度,
MD(机器方向) ——沿着纤维,
CD(横向) —— 横跨纤维。
该模型必须与正交各向异性线弹性一起使用,其中 v_13=v_23=0。它描述了平面内、厚度方向以及面外剪切的解耦行为。
夏庆生提出的面内响应塑性准则引入了塑性“泊松比”的概念。面内行为的屈服面由六个子面组成:
MD 方向张力
CD 方向张力
MD-CD 方向正剪切
MD 方向压缩
CD 方向压缩
MD-CD 方向的负剪切
屈服面:
σ_xy=0 的屈服面:
使用以下关键字完成定义:
*PAPERBOARD PLASTICITYν_1p, ν_2p, ν_4p, ν_5p
其中:ν_1p - 方向 1 拉伸的塑性“泊松比”,ν_2p - 方向 2 拉伸的塑性“泊松比”,v_4p - 方向 1 压缩的塑性“泊松比”,v_5p - 方向 2 压缩的塑性“泊松比”。
该模型还考虑了平面内硬化,由以下关键字定义:
*PAPERBOARD HARDENINGσ_01, A_1, B_1, C_1, σ_02, A_2, B_2, C_2σ_03, A_3, B_3, C_3, σ_04, A_4, B_4, C_4σ_05, A_5, B_5, C_5
硬化参数根据相应的屈服子表面进行编号。
厚度方向 (ZD) 的响应在拉伸时为线弹性响应,在压缩时为弹塑性响应。弹性部分使用以下关键字定义:
*PAPERBOARD THICKNESS COMPRESSION ELASTICµ, σ_tZ
其中 µ 和 σ_tZ 是定义 ZD 方向弹性行为的参数。
该方向的硬化定义为:
*PAPERBOARD THICKNESS COMPRESSION PLASTICITYA_σ, B_σ, C_σ
其中:A_σ、B_σ、C_σ——硬化参数。
最后考虑的行为是平面外剪切,使用以下关键字定义:
*PAPERBOARD TRANSVERSE SHEAR PLASTICITYA_τ, B_τ, C_τ
其中 A_τ、B_τ、C_τ 为硬化参数。
*Material, name=paper*Density1e-03,*elastic,type=engineering constants5600., 2000., 100., 0.37, 0., 0., 4100., 400.400.*paperboard plasticity,exponent=40.5, 0.133,0.5, 0.133,*paperboard hardening12, 19.,260.,800.,6.5,7.4,160.,160.6.,6.,375.,200.,7.3,6.,160.,300.6.3,9.,310.,225.,*paperboard thickness compression elastic, void ratio = 1.120.026,0.36*paperboard thickness compression plasticity0.4378, -1., 6.5,*paperboard transverse shear plasticity5., 3., 2.,
纸板塑性模型全面描述了纸板在各个方向上的行为。与通用塑性模型(通常针对金属设计)不同,该模型专为模拟此类材料而设计,非常适合纸板包装分析。
与任何新的材料模型一样,建议先测试单个单元或基本样本,以了解其行为。