[论文解读] Elastic properties of cellular polypropylene films: Finite element simulations and their comparison with experiments

Enis Tuncer, Michael Wegener|arXiv (Cornell University)|Apr 7, 2004

Cellular and Composite Structures参考文献 20被引用 5

一句话总结

本研究建立了一个有限元模型，用于模拟二维桁架状多孔聚丙烯薄膜的有效杨氏模量，提出了一种依赖于固体体积分数和微观结构几何形状的幂律关系。该模型与实验数据高度吻合，表明胞壁厚度和长宽比（a/b）对刚度有关键影响，最优刚度出现在 a/b ≈ 2 且 t/a ≈ 0.075 时。

ABSTRACT

The Young's modulus of a two-dimensional truss-like structure is simulated by using the finite element method. A power-law expression is proposed for the effective Young's modulus of the system. The obtained numerical results are compared with the experimental data of the {\em anisotropic thin cellular polypropylene films}. At high solid volume fractions ($>0.4$), the average shape of the cells are lateral, and their dimensions have around one-to-five ratio. As the samples are inflated further, volume fraction of the solid is decreased, the average shape approach a diamond-like structure with one-to-two ratio. In addition the effective Young's modulus of the system increases. It is concluded that valuable structural information can be obtained by analyzing the experimental data and the numerical simulations, which take into account the material's micro-structural information, simultaneously.

研究动机与目标

• 使用有限元分析模拟薄型各向异性多孔聚丙烯薄膜的有效弹性性能。
• 开发一种均质化框架，以捕捉微观结构几何形状对宏观刚度的影响。
• 将数值预测结果与不同固体体积分数下杨氏模量的实验测量值进行比较。
• 通过基于仿真的逆分析方法，从实验数据中提取微观结构参数（如胞元长宽比、壁厚）。
• 理解在膨胀过程中胞壁拉伸与结构演化对力学响应的作用。

提出的方法

• 在具有周期性边界条件的二维桁架状结构代表性单元胞上进行有限元（FE）模拟。
• 通过施加单轴载荷并测量模拟域内的应变响应，计算有效杨氏模量。
• 提出幂律表达式：$ E^{ ext{e}} = ilde{C} q^n E_s $，其中 $ q $ 为固体体积分数，$ E_s $ 为固体模量，$ n $ 取决于微观几何形状。
• 通过改变胞元长宽比（a/b）和归一化壁厚（t/a），生成等厚度线以与实验结果对比。
• 将薄型聚丙烯薄膜的实验数据（1–10 MPa，$ E_s acksimeq 1 $ GPa）进行归一化，并拟合为二次对数表达式。
• 通过对比模拟的刚度趋势与实验数据在膨胀诱导的微观结构变化下的表现，验证模型。

实验结果

研究问题

• RQ1 多孔聚丙烯薄膜的有效杨氏模量如何依赖于固体体积分数和微观结构几何形状？
• RQ2 有限元模型能否准确预测具有演化微观结构的各向异性、薄型多孔薄膜的弹性响应？
• RQ3 哪些微观结构参数（如胞元长宽比、壁厚）控制了膨胀过程中观察到的刚度转变？
• RQ4 在高固体分数下，胞壁拉伸对力学性能的影响程度如何？
• RQ5 能否通过逆向建模，利用实验刚度数据推断潜在的微观结构特征？

主要发现

• 多孔聚丙烯薄膜的有效杨氏模量对固体体积分数 $ q $ 呈幂律依赖，对数-对数二次拟合结果为：$ \text{log}(E^{ ext{e}}/E_s) = 10x^2 - 8.9x - 0.84 $，其中 $ x = q $。
• 当 $ q > 0.4 $ 时，微观结构以横向胞元为主，长宽比 $ a/b acksimeq 5 $，且胞壁拉伸显著影响刚度。
• 当 $ q < 0.4 $ 时，数据呈现弯曲趋势且 $ dE^{ ext{e}}/dq < 0 $，表明向刚度更高的菱形结构转变，此时 $ a/b acksimeq 2 $，壁厚拉伸效应减小。
• 模型与实验数据高度一致：在 $ t/a = 0.075 $ 时的模拟等厚度线与实验的二次拟合结果极为接近。
• 刚度转变的临界壁厚估计为 $ t/a acksimeq 0.05 $，低于此值时，尽管固体分数较低，材料仍表现出更高刚度。
• 实验数据分布在 $ 2 acksimeq a/b acksimeq 6 $ 和 $ 0.03 acksimeq t/a acksimeq 0.12 $ 范围内，表明制造过程会选择特定的微观结构参数。