包装工程专业英语大作业俞坚道

1、宁波理工学院 包装工程专业外语大作业 题 目聚氨酯泡沫应变率敏感性的本构模型 姓 名 俞坚道 学 号 3120614033 专业班级 12包装工程1班 指导教师 张 萍 学 院 机电与能源工程学院 目 录1.引言32.本构方程3 2.1准静态模型3 2.2包括应变率效应的模型43.实验与参数辨识6 3.1低密度聚氨酯泡沫塑料模型6 3.2高密度聚氨酯泡沫塑料7 3.3应力与应变率84.总结9参考文献10作者介绍10附件：原文A constitutive model forpolyurethane foam with strain rate sensitivity11聚氨酯泡沫应变率敏感性的本构

2、模型 摘 要：本工作研究了聚氨酯泡沫塑料的应变率相关行为并制定了新的本构模型，以提高在各种应变速率中实验数据的拟合性。该模型的七个参数被两种应变速率下进行的准静态压缩试验所决定。压缩试验两种应变速率。高、低密度聚氨酯泡沫塑料的两种模型显示在了不同的应变率下应力-应变的关系。所进行的动态压缩试验得出了在高应变率下应力-应变的数据并且将此结果同本构模型进行了比较。 关键词：聚氨酯泡沫；本构模型；动载荷；应变率 1.引言 对于汽车耐撞性能的改进不能被过高估计。美国交通部门估计从1999年到2009年，有超过400,000人遇难，2000万受伤需要住院治疗。这个情况，连同一系列的环境问题和社会压力在法

3、律的支持下已经领导，并将继续引导，高度创新的设计，涉及到有色金属合金、智能结构、复合材料等先进的材料和泡沫。这项研究特别感兴趣的是耐撞性结构泡沫塑料材料的使用。聚合物泡沫目前被用作在保险杠和作为填充材料加固屋顶和门梁，并且它们的应用程序将被不断扩展。它们可以加强薄弱部位结构，从而可以对冲击载荷做出有效反馈，即加强耐撞性。泡沫的能量吸收性能是指，当保持一个几乎恒定应力值时所经历的大变形的能力。在其他聚合物中，聚氨酯泡沫被广泛应用在工业上。泡沫已经成为了无数试验、数值和理论调查的主题。Shim等人对硬质聚氨酯泡沫塑料进行了速度从2到4m/s 的普通冲击试验。Avalle等人研究了使用能量吸收图的聚

4、合物泡沫的能量吸收特性。Rusch提出了压应力与应变之间的关系作为初始压缩模量的泡沫和应变的函数。Meinecke和Schwaber描述作为一个应变函数的初始系数并改变了作为一个幂级数的应变的应变函数形式。Sherwood和Frost修改了聚氨酯泡沫塑料在受压状态下的本构模型的形状函数。Langseth进行了用挤压状态下的铝在在不同速率下持续冲击的冲击试验。 2.本构方程 2.1准静态模型聚合物泡沫的应力应变曲线，在一般情况下显示出三个区域：一个线性弹性区域，其次是可塑性平坦区域以及应力陡峭上升的致密区。吉布森模型的三个方程描述了每个区域，如下所示当y时，=E (1)当yD1-D-1+y时，=

5、y (2)当>D1-D-1+y时，=y1D(DD-)m (3) 和 分别是工程应力和工程应变。该模型具有五个参数: E为杨氏模量的弹性部分。y为屈服应力。D为致密化应变和两个常量D、m。这模型在平坦区域应力有常量值的限制并且在两个区域的边界应力应变曲线是不光滑的。Rusch 模型可以描述的两个力的总和：=a-p+bq (4)这里的四个参数，a、b、p和q，可以根据经验确定。虽然这个模型通过一个方程描述了应力-应变关系，它的缺点在于不能精确描述致密化的阶段。为更好适用泡沫材料在实验中的应力应变关系。刘和苏巴斯建议式 (5) 所示的模型。这种模式具有六个参数。设置B作为一个单元，把它变成一个

6、五参数模型。=Ae-1B+e+eC(e-1) (5) 最近， 一个模型被Avalle 等人提出。如下所示。=A1-e-EA1-m+B(1-)n (6) 该模型还具有五个参数。参数E、A和B与密度相关而m和n不是。第一项适合弹性区和平坦区，而第二项适合致密区。虽然这个模型很好的描述了聚合物泡沫的应力-应变曲线，但它没有包括因动态加载下对于应变率的影响。 2.2包括应变率效应的模型 包括应变率效应的模型由Nage等人提出。d=fM(,) (7)这里的d是工程压应力，f是应变的形状函数，代表了在参考应变速率 0 和M(,)下应力应变的关系。应变速率的功能，是一种应变和应变率的函数，在参考应变率时的单

7、位价值。Sherwood和Frost所描述的形状函数的十分之一阶多项式函数的应变，其系数是确定从压缩试验的参考应变率。应变率函数的假设，即模型图的应力与应变率的假设，表现出的线性函数的应变可以表示如下：f=n=010Ann (8)M,=(0)a+b (9) 双参数，a和b，可以发现通过线性近似的实验的应力-应变率数据的几个应变率。由于描述Eqs的形函数。（7）和（8）只不过是一个应力-应变方程，它取代模型（6）的五个参数比使用十分之一阶多项式更有效。应变率的函数方程（9）是基于在一定的应变-应力数据曲线情况下具有线性应变函数的假设。仔细检查实验数据表明，这是不正确的，特别是对于半柔性聚氨酯。引

8、入了一个新的应变率函数来描述应力与应变率的对数曲线图，该曲线是凹的，在恒定应变下：M,=1+(a+b)ln(0) (10) 通过替换式的形状函数和应变速率的公式（6）、函数（7）和式（10），分别为应力-应变方程，包括应变率效应推导： d=()1+(a+b)ln(0) =A1-e-EA1-m+B(1-)n1+(a+b)ln(0) (11) 该模型方程（11）有七个参数。这些参数的确定分为2个步骤。首先，在参考应变率下，五个参数，E、A、B、m和n，确定通过拟合实验的应力-应变数据。然后，剩余的2个参数，a和b，是通过线性近似的实验的应力-应变数据的另一个应变率被发现的。大多数的参数被认为是依赖

9、于材料的密度和温度。因此，两步实验应该在同一密度和温度下进行。表1 从曲线拟合计算的五个参数ABEmn0.38630.315113.99963.41621.3680图1 聚氨酯泡沫塑料试样 3.实验与参数辨识 3.1低密度聚氨酯泡沫塑料模型 准静态测试是在位移测量装置和力传感器安装形成的室温下，使用具有100kN大负荷容量的MTS 810机器进行的。密度为67 kg/m3的开孔式圆柱聚氨酯泡沫（如图1所示，由韩国lacomtech有限公司提供）是在室温下测试。试样的直径和长度分别为42毫米和40毫米，分别为。准静态试验在0.001s-1的应变速率进行了应力应变实验数据。五个参数方程（6）采用最

10、小二乘法进行数据拟合。参数的计算会被拟合曲线点的数目所影响。表1显示的参数，是通过拟合33个应力-应变实验数据点得到的。实验和装配的应力-应变曲线，如图2所示。拟合曲线很接近实验曲线。（11）式中a和b，这两个参数的确定，在另一个应变率的实验数据是必要的。这个实验用密度同为67 kg/m3的聚氨酯泡沫在应变速率为0.1s-1的情况下进行的。因为在应变速率为0.001s-1时采用最小二乘法得到五个参数，所以a和b两个参数是由方程的本构方程拟合实验数据确定。表2中的数据通过最小二乘法拟合实验数据计算得出的参数。如图3所示，应力-应变曲线很接近实验数据。表2 从曲线拟合计算的二个参数ab0.0401

11、0.1222 图2 应力-应变关系 图3 应力-应变关系 ( =67Kg/m3, 0=0.001 s-1) ( =67Kg/m3, 0= 0.1 s-1) 获得构成所有参数接近的亲双准静态实验的本构模型，所得到的应力-应变曲线可以得到任何应变率。如图4所示，0.001，0.01，0.1，1，10，和100s-1这六个应变率的值分别代入式（11）得到六条应力-应变曲线。其中应变率为0.001 s-1对应最低的曲线。 对于在相同的应变下，当泡沫的应变速率更高时，更多的能量将被储存在泡沫中。用本构方程拟合得到七个参数，用本构方程拟合得到了67 kg/m3的聚氨酯泡沫体在动态载荷作用下的应力应变关系。

12、经验证，方程可以采用动态压缩试验。 图4 应变率的应力-应变关系 图5 应力-应变关系 ( =67Kg/m3) ( =67Kg/m3, 0= 97 s-1) 落塔型冲击试验机（英斯特朗Dynatup 9250 HV）进行动态测试。本机可提高离散变量到一个特定的高度，下降到试样使用引力或加速弹簧。当26.5Kg的重物冲击试样，测量的冲击速度为3.88米/秒，并计算应变速率为97s-1。图5中的实心线显示了实验的应力-应变关系，虚线描述了本构模型的关系。虽然它不遵循的动态效果的振荡，但是该模型接近的实验结果。 3.2高密度聚氨酯泡沫塑料 图6 应力-应变关系 图7 应力-应变关系 ( =89Kg/

13、m3, 0= 0.001 s-1) ( =89Kg/m3, 0= 0.1 s-1) 进行了89 kg/m3的高密度聚氨酯泡沫塑料的准静态测试。第一，在应变速率为0.001 s-1应力-应变的实验数据通过采用最小二乘法拟合得到方程（6）的五个参数。实验和拟合应力-应变曲线如图6所示。 拟合曲线很接近实验曲线。另一个静态应变率的实验数据都需要找到两参数方程a和b（11）。在应变速率为0.1 s-1的准静态试验，通过曲线拟合得到的参数。如图7所示，应力-应变曲线很接近实验数据。 由于所有的本构模型的参数已被发现，应力-应变曲线可以得出任何应变率。将0.001，0.01，0.1，1，10，和100 s

14、-1这六个应变率分别代入式（11），在图8中绘制成应力-应变曲线。 图8 应变率的应力-应变关系 图9应力-应变关系 ( =89Kg/m3) ( =89Kg/m3, 0= 97 s-1) 要确定该模型是否能在高应变率和动态压缩试验中工作。当重物冲击试样时，撞击速度为3.88米/秒，计算为应变速率为97s-1。实验绘制的应力-应变曲线与所提出的本构模型进行比较。图9中的实心线显示了实验的应力-应变关系，虚线描述了本构模型的关系。 在非常低的应变率的情况下，由于本构模型的参数是由2个准静态实验得到，这种模型不能跟随动态效果的振荡实验。图5和图9的实验结果表明，试样的密度越高，进行的的振荡越大。 3

15、.3应力与应变率 得到模型的所有参数后，可以绘制在不同规模水平下参数记录的应变率。如图10所示，67 kg/m3的低密度聚氨酯泡沫。如图11所示，89 kg/m3的高密度泡沫的应力-应变率。每个曲线几乎呈线性关系，但有点凸。Nage et al提出的本构模型，如式（7）通过假设的应力与应变率曲线表现出应变的线性函数。他们的实验的日志记录图一般遵循直线，这往往遵循的提出模型的特点，但凸点最多的是泡沫材料。 图10恒定应变下的应力-应变率 图11恒定应变下的应力-应变率 ( =67Kg/m3) ( =89Kg/m3) 4.总结 在设计吸能结构泡沫时，需要对泡沫的力学特性进行数学模型描述。几个模型以

16、及描述的应力-应变关系的泡沫是工作在一个特定的应变率。本文提出了一种用于聚氨酯泡沫的本构模型，以描述的应力-应变关系在很宽的范围内的应变率。 该模型有七个参数。通过两步优化程序的参数进行评估。（1）第一步是参考应变率和发现五个参数。（2）在下一步中，通过将所提出的模型和实验测得的值之间的曲线拟合，发现在这步中的应变率的大小影响的剩余2个参数。 在初始应变率为97s-1时，动态压缩试验被提供应力-应变下进行，这些数据绘制出他们所提出的模型。应力-应变曲线的提出模型以及如下的实验数据，即使它并没有准确地描述振荡，是因为系统所造成的冲击。参考文献1 US Department of Transpor

17、tation, FARS/GES 2009 DataSummary, DOT HS 811 401 (2011).2 U. K. Lee, S. B. Jang and S. S. Cheon, Development of asubframe type fuel tank for passenger cars, Journal of Mechanical Design ( Transaction of ASME ), 132 (4) (2010) 044501.3 V. P. W. Shim, Z. H. Tu and C. T. Lim, Two-dimensional response

18、of crushable polyurethane foam to low velocity impact, International Journal of Impact Engineering , 24 (2000) 703-731.4 M. Avalle, G. Belingardi and R. Montanini, Characterization of polymeric structural foams under compressive impact loading by means of energy absorption diagram, International Jou

19、rnal of Impact Engineering , 25 (2001) 455-472.5 K. C. Rusch, Load-compression behaviour of flexible foams, Journal of Applied Polymer Science , 13 (1969) 2297-2311.6 E. A. Meinecke and D. M. Schwaber, Energy absorption in polymeric foams, Journal of Applied Polymer Science , 14 (1970) 2239-2248.7 J

20、. A. Sherwood and C. C. Frost, Constitutive modeling and simulation of energy absorbing polyurethane foam under impact loading, Polymer Engineering and Science , 32 (1992) 1138-1146.8 M. Langseth, O. S. Hopperstad and T. Berstad, Crashwor thiness of aluminium extrusions: Validation of numerical simu

21、lation, effects of mass ratio and impact velocity, International Journal of Impact Engineering , 22 (1999) 829-854.9 L. J. Gibson and M. F. Ashby, Cellular solids: Structure and properties , Pergamon Press (1997).10 G. Subhash and Q. Liu, Quasistatic and dynamic crusha bility of polymeric foam in rigid confinement, International Journal of Impact Engineering , 36 (2009) 1303-1311.11 M. Avalle, G. Belingardi and A. Ibba, Mechanical models of cellular solids: Parameters identification from experiment