发泡聚苯乙烯缓冲性能分析方法本科毕业论文1

毕业设计（论文） 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 目 发泡聚苯乙烯缓冲性能分析方法 姓 名 俞坚道 学 号 3120614033 　 专业班级 12包装工程1班 指导教师 卢富德 学 院 机电与能源工程学院 完成日期 2016年4月30日 毕业设计（论文）任务书 1、 毕业设计（论文）题目 发泡聚苯乙烯缓冲性能分析方法 2、 学生姓名 俞坚道 学号 3120614033 专业班级 12包装（1）班 任务书发放日期 2015年10月 3、 指导教师对毕业设计（论文）的进度安排及任务要求 1． 主要任务与目标 （1）建立MATLAB/GUI界面，求解EPS最大加速度-静应力曲线。 （2）建立EPS的缓冲本构模型，求解易损件系统的冲击响应。 2． 主要内容与基本要求 （1）在产品-EPS衬垫模型，首先获取所发表文献中的有关EPS的 数据，并汇总不同密度的 曲线，构建EPS材料库。 （2）在易损件-质量主体-EPS衬垫模型，首先得到EPS的本构模型，然后得到EPS传递给产品质量主体的加速度离散脉冲，另外根据三次样条插值计算得到，脉冲激励下对易损件系统的冲击响应。 基本要求： （1）查找和翻译与课题相关的英文文献； （2）文献综述≥300 0字 个人自传范文3000字为中华之崛起而读书的故事100字新时代好少年事迹1500字绑架的故事5000字个人自传范文2000字 （后附所查阅的文献总数≥20篇，其中外文≥3篇）； （3）外文翻译≥3000字； （4）开题报告≥3500字（后附文献总数≥20篇，其中外文≥3篇）； （5）界面设计和相关程序； （6）毕业论文≥1万字。 3． 计划 项目进度计划表范例计划下载计划下载计划下载课程教学计划下载 进度 （1）2015年10月至2016年1月：阅读文献，建立MATLAB/GUI界面，求解EPS最大加速度-静应力曲线。 （2）2016年1月至2016年3月：建立EPS的缓冲本构模型。 （3）2016年3月至2016年5月：撰写论文，准备答辩。 4． 主要参考文献 [1] Liu Q, Subhash G. A phenomenological constitutive model for foams under large deformation[J]. Polymer Engineering and Science, 2004, 44 (5): 463~473. [2] Liu Q, Toole B O, Behavior pattern and parametric characterization for low density crushable foams[J]. Journal of Materials Processing Technology. 2007, 191( 1-3)：73~76. [3] Subhash G, Liu Q L. Quasistatic and dynamic crushability of polymeric foams in rigid confinement[J]. International Journal of Impact Engineering, 2009, 36(10-11): 1303~1311. [4] Avalle M, Belingardi G, Ibb A. Mechanical models of cellular solids: Parameters identification from experimental tests[J]. International Journal of Impact Engineering, 2007, 34(1): 3~27. [5] Sherwood J A, Frost C C. Constitutive modeling and simulation of energy absorbing polyurethane foam under impact loading[J]. Polymer Engineering and Science, 1992, 32(16): 1138~1146. [6] Jeong K Y， Chon S S， Munshi M B. A constitutive model for polyurethane foam with strain rate sensitivity[J]. Journal of mechanical Science and Technology, 2012, 26(7): 2033~2038. [7] Zhang J, Kikuchi N, Li V, etal. Constitutive modeling of polymeric foam material subjected to dynamic crash loading[J]. International Journal of Impact Engineering 1998, 21(5): 369~386. 起讫日期： 2015年10月20日 至 2016年6月3日 指导教师签名 职称 年 月 日 4、 分院审核意见 负责人签名 年 月 日 毕业设计(论文)开题报告 (含文献综述、外文翻译) 题 目发泡聚苯乙烯缓冲性能分析方法 姓 名 俞坚道 学 号 3120614033 　 专业班级 12包装工程1班 指导教师 卢富德 学 院 机电与能源工程学院 完成日期 2016年3月5日 目 录 I 目 录 1文献综述 11. 引言 12. 聚苯乙烯泡沫(EPS)的定义、性能 12.1 聚苯乙烯泡沫(EPS)的定义 12.2 聚苯乙烯泡沫(EPS)的性能 22.2.1 聚苯乙烯泡沫的化学性能 22.2.2 聚苯乙烯泡沫的物理性能 33. 聚苯乙烯泡沫(EPS)的发展和应用 33.1 聚苯乙烯泡沫(EPS)的发展 33.2 聚苯乙烯泡沫(EPS)的应用 33.2.1 聚苯乙烯泡沫（EPS）的具体应用 43.2.2 苯乙烯泡沫（EPS）产品废弃物的具体应用 44. 聚苯乙烯泡沫(EPS)的国内外的研究现状 44.1 聚苯乙烯泡沫(EPS)的国内研究现状 64.2 聚苯乙烯泡沫(EPS)的国外研究现状 75. 总结 8参考文献 10开题报告 101. 研究意义 102. 研究背景 143. 主要研究工作 143.1 研究内容 143.2 研究重点与难点 143.2.1 编辑求解EPS缓冲曲线程序 143.2.2 编辑求解EPS单自由度缓冲作用下衬垫设计 153.2.3 编辑求解EPS二自由度缓冲作用下易损件的加速度响应程序 154. 研究方法与实施进度计划 16参考文献 18外文翻译和原稿 文献综述 发泡聚苯乙烯缓冲性能分析方法 1. 引言 聚苯乙烯泡沫是一种轻型高分子聚合物。它是采用聚苯乙烯树脂加入发泡剂，同时加热进行软化，产生气体，形成一种硬质闭孔结构的泡沫塑料。EPS泡沫是一种热塑性材料，白色珠状颗粒，每立方米体积内含有300-600万个独立密闭气泡，内含空气的体积为98%以上,相对密度1.05，由于空气的热传导性很小，且又被封闭于泡沫塑料中而不能对流，所以EPS是一种隔热保温性能非常优良的材料。这种均匀封闭的空腔结构，同时也使EPS泡沫具有吸水性小，质量轻及较高的机械强度等特点。EPS是当前最轻的包装材料，在冰箱、冷柜、冷藏库等方面广为应用，土木建筑中也用它来作保温隔热材料。 2. 聚苯乙烯泡沫(EPS)的定义、性能 2.1 聚苯乙烯泡沫(EPS)的定义 聚苯乙烯(英文polystyrene，简称PS)是由含有不饱和键的苯乙烯(英文styrene)单体在引发剂(或催化剂)存在下进行聚合反应生成。化学反应方程式如下: 聚苯乙烯泡沫(英文Expanded polystyrene，简称EPS)是由聚苯乙烯(PS)经加热发泡后形成的具有微细闭孔结构的泡沫塑料[1] 。 2.2 聚苯乙烯泡沫(EPS)的性能 EPS由聚苯乙烯颗粒发泡而成，根据发泡的方式分为:模型中发泡和挤出法发泡，在EPS成型过程中，聚苯乙烯颗粒中的戊烷受热汽化，在颗粒中膨胀形成许多封闭的空腔，正是这种均匀的封闭空腔结构决定了EPS具有许多材料所没有的特性。 2.2.1 聚苯乙烯泡沫的化学性能 EPS的化学特性从其本质上说与聚苯乙烯树脂相同，对于一般的酸、碱、动植物油、盐类等有较好的抗化学性，而对于芳香族碳化氢、卤族碳化氢、酮类等矿油系药品具有易溶解的性质。因此，要注意防止与这些物质接触。另外，EPS具有耐久性，在自然环境下具有耐霉变、不受白蚁的影响 [2] 。 2.2.2 聚苯乙烯泡沫的物理性能 （1）密度。EPS的密度由成型阶段聚苯乙烯颗粒的膨胀倍数决定，一般介于10～45kg/m3之间，作为工程中使用的EPS其表观密度一般在15～30kg/m3。密度是EPS的一个重要指标，其各项力学性能几乎与它的密度成正比关系 [3] 。 （2）保温隔热性能。EPS具有良好的保温性能是因为它由完全封闭的多面体形蜂窝构成，蜂窝的直径0.08～0.15mm，蜂窝壁厚为0.001mm。EPS由约98%的空气和2%的聚苯乙烯组成。截留在蜂窝内的空气是一种热的不良导体，因而对EPS的热绝缘（保温）性能起决定性的作用。与含有其他气体的泡沫塑料不同，空气长期留在蜂窝内，所以保温效果稳定不变。 （3）力学性能。EPS泡沫塑料的变形量是随负载时间的增大而加大。只有在生产两天后才能达到它的承载能力。刚刚出模的制件对压力很敏感，因为泡孔中蒸汽和残留发泡剂的冷凝会导致部分真空，这要等到泡孔吸入空气后才能达到压力平衡。EPS泡沫塑料属于硬性泡沫塑料类。在荷载情况下，它呈现粘弹性，这是一种脆硬性材料所具有的特性。弯曲强度、压缩强度和抗拉强度都与EPS泡沫塑料的表观密度成正比。 （4）吸水性能。与其他许多泡沫塑料不一样，EPS泡沫塑料是不吸湿的，即使将它浸没在水中，也仅仅吸很少量的水，由于蜂窝壁是不透水的，水仅能从熔融的蜂窝之间微小通道透入泡沫塑料，不言而喻，吸收的水量（透入的水量）取决于EPS泡沫塑料原材料在加工时的性能和加工条件（特别在发泡时）。随着EPS泡沫制件密度的提高，水的吸收和水蒸汽的透过率下降，水蒸汽的扩散阻力系数上升。 （5）导电性能。由于EPS泡沫塑料具有较低的导电性能，因此它易产生自身带电现象，特别是其生产过程中难免要遇到不同程度的摩擦。这样，带静电的珠粒会吸附在管道中，造成管道的堵塞；而制成的包装产品带上静电，既容易吸尘，又会对某些电器产品造成影响。为此，完全有必要对一些EPS包装增添防静电保护剂，即增加EPS泡沫塑料自身的导电性能。 3. 聚苯乙烯泡沫(EPS)的发展和应用 3.1 聚苯乙烯泡沫(EPS)的发展 自上世纪50年代由德国BASF公司开发EPS珠粒生产工艺后，EPS泡沫塑料由于成型工艺简单及设备简易可行，并可制成各种形状、不同密度的产品，因而发展迅速。现在EPS已成为苯乙烯树脂三大产品(GPPS、HIPS、EPS)中重要的产品之一。我国EPS工业从1958年自行研制的悬浮聚苯乙烯塑料上市至21世纪的今天，EPS得到空前的发展，著名的厂家有龙王、兴达、台达等 [4] 。 3.2 聚苯乙烯泡沫(EPS)的应用 EPS泡沫塑料为轻质半硬性多孔材料，它由无数闭孔组成，含97%以上的空气，具有导热系数低、冲击强度高、吸水性小、隔音和耐老化等特点。改性后具有抗静电、阻燃等性能，可广泛应用于缓冲包装、漂浮器具、保温隔热和建筑材料等。 3.2.1 聚苯乙烯泡沫（EPS）的具体应用 (1)房屋建筑、冷热管道、汽车、机车、船舶、冷藏库、电器、军用兵器、炮弹等的外包装，以达到保温、隔热、隔音、防震的目的。 (2)用于救生设备、浮标以及海上围油栏以减少油船卸油时对水域的污染，是优良的漂浮材料。 (3)代替铸造机械部件中木模，用于实型铸造和精密铸造。 (4)与其他多种材料进行复合，制成建筑物中内外墙体材料，代替红砖使用。 (5)用于美术艺术雕刻造型，舞台布景，电影、电视道具，制作动物标本及工艺美术品的内衬等。 (6)由于其微波穿透损耗率极低，可用作雷达和微波通讯部件。 (7)其泡沫珠粒还是一种新型的过滤材料，可用于种子发芽和水面育秧等。 3.2.2 苯乙烯泡沫（EPS）产品废弃物的具体应用 (1)捣碎后作土壤改良剂，用于改进粘性太大的土壤。还可与粘土混在一起用作制砖或直接与水泥等混合，浇成轻质水泥地面、屋顶或预制板。这种材料具有阻燃性并能防霜。 (2)用作路基底层修筑路面，能使修建在承载能力很差的土壤上的道路具有优良的性能。 (3)常温下裂解制成PS树脂，可替代传统的醇酸漆、酚醛漆等产品。 (4)用解溶剂型解质法消泡工艺进行消泡回收处理，消泡后的物料运往有关塑料制品厂和粘合剂厂使用，效果良好。 (5)利用固体废弃物再生资源处理设备及技术将废EPS和废木屑混合，可生产木屑含量超过50%的超级木材，这种再生木材制成的器具具有天然林木的全部功能 [5] 。 4. 聚苯乙烯泡沫(EPS)的国内外的研究现状 4.1 聚苯乙烯泡沫(EPS)的国内研究现状 佟富强等[6]等通过对聚苯乙烯泡沫材料进行不同形变速率的压缩实验，经过数学处理得到不同形变速率条件下材料的缓冲系数—最大应力曲线，找出压缩速度对材料缓冲性能影响的变化规律，可用以指导实际应用。 阳以本[7]介绍了我国现阶段普遍采用的螺旋传动，蒸汽加热发泡工艺生产的可发性聚苯乙烯泡沫材料密度与压缩强度(压缩50%)的关系，分析成果应用的可行性给其摩托车生产工厂带来好的经济效益(EPS应用于摩托车成品包装)。 杜骋等[8]在《聚苯乙烯泡沫(EPS)的特性及应用分析》文中也指出聚苯乙烯泡沫是一种性能优良的路基轻质填料，具有轻质、高强、较强的化学稳定性和水稳定性、良好的力学性能且施工方便简单等优点，在国外道路工程中有较为广泛的应用。我国对EPS的研究和应用较少，文章对EPS的物理化学性能、力学性能、EPS作为路基轻质填料的结构设计方法、EPS在道路工程中的应用等方而作了较为全而的介绍和分析，对我国使用EPS有借鉴作用。 毛快[9]主要研究了EPS的压缩蠕变性能。EPS材料作为一种热塑性材料，其特点决定其在使用过程中与其他黏弹性材料一样会发生蠕变和松弛现象，因此要考虑EPS材料的性能随时间的变化规律。因为EPS作承重构件时通常受压，长期蠕变则需要几天或更长的时间，且发现仅当应力超过阀值水平时才产生。针对不同密度的EPS进行了短期压缩蠕变实验研究，讨论应力水平和密度对蠕变的影响。 朱向荣[10]采用三轴试验研究EPS在不同加载频率、不同循环次数和不同围压的下强度和模量等的变化规律，为EPS工程应用提供了一些理论依据。 卢富德等[11]针对泡沫结构的压缩性能，介绍了泡沫结构的试验压缩响应与数值模拟方面的研究进展。结果表明：相对密度、环境温度、压缩应变率及微观结构等参数对泡沫的缓冲性能影响显著。相对密度越大，泡沫结构的屈服应力越大，吸收能量的能力越大；泡沫的吸收能量能力一般随环境温度的增加而减小；由于泡沫的基体材料表现率相关性，应变率增加导致吸收能量能力的增加；当微观结构不同时，泡沫的细胞分布导致结构的缓冲性能差异显著。 高德等[12]根据植物秸秆纤维聚苯乙烯缓冲材料静态压缩实验数据，研究了材料的缓冲性能，并考虑植物秸秆纤维的影响，对Sherwood- Frost本构关系框架进行了扩充，建立了非线性本构关系模型，并利用实验数据成功识别模型参数此种描述植物纤维类材料非线性力学行为的方法，为进一步研究和开发植物纤维聚苯乙烯材料提供了理论基础。 都学飞等[13]比较分析了４种厚度EPS缓冲包装材料的压缩变形回复性、外力－位移曲线等性能。结果表明，EPS缓冲材料的回复性、永久变形与材料的厚度有很大关系，而它的外力－位移曲线的走向大致相同，厚度越大吸收的能量越多，缓冲性能越好。 龙志坚[14]制备了发泡量相同的微发泡聚苯乙烯, 采用扫描电镜(SEM)和Image-pro图像处理软件对微发泡聚苯乙烯的微孔尺寸进行了观察和统计;建立微球模型, 分析了微孔尺寸大小对微发泡聚苯乙烯力学性能的影响。结果表明:细小而均匀的泡孔对微发泡聚苯乙烯力学性能的提高有较明显的促进作用, 微球模型的计算结果与宏观力学性能的影响规律有很好的重现性。 4.2 聚苯乙烯泡沫(EPS)的国外研究现状 Horvath[15]用边长为5cm的EPS立方体试件在应变速率10 mm/min的条件下，采用应变控制形式进行了无侧限单轴压缩试验，得到压缩应力-应变曲线并对该曲线进行了分析。但没有对多种密度与多种加载速率的情况进行试验比较和分析。 Duskov[16]采用直径为10 cm、高为20 cm的圆柱体EPS试件在20 kPa的作用下进行蠕变研究，得出蠕变曲线，从蠕变曲线的分析中可以看出：EPS材料的蠕变主要发生在加载初期，随着加载时间的增长，蠕变的速率趋于稳定，在加载1年以后蠕变的速率几乎接近常数。 Kwang Young Jeong[17]研究了聚氨酯泡沫塑料的应变率相关行为并制定了新的本构模型，以提高在各种应变速率中实验数据的拟合性。该模型的七个参数被两种应变速率下进行的准静态压缩试验所决定。压缩试验两种应变速率。高、低密度聚氨酯泡沫塑料的两种模型显示在了不同的应变率下应力-应变的关系。所进行的动态压缩试验得出了在高应变率下应力-应变的数据并且将此结果同本构模型进行了比较。 Qunli Liu[18]认为泡沫在受到大的变形的情况下，加压力负载时，可在五参数模型完全确定应力-应变响应的三个基本特征，即线性、可塑性状应力平台和致密化的阶段。此外，根据不同的屈服强度和硬化状或软化的约束条件，该模型的参数可以系统变化，以确定泡沫的初始密度。 James M. Gibert[19]探讨了压盘影响开孔泡沫缓冲材料动力学的定性认识。一个超弹性材料模型是用来描述阻尼和滞回特性为线性粘弹性的非线性泡沫的应力–应变关系。使用一个简单的非线性不连续模型的跌落试验以及数值模拟，研究探讨对其的物理影响。数值研究表明，该模型能够提供预测的冲击脉冲的形状、持续时间和幅度，在不同的静态应力和下降高度。模型所产生的动态缓冲曲线保留凹向上的“槽”的实验曲线的形状特征。此外，该模型表明，给定的下降条件下，冲击吸收的最佳振幅取决于泡沫的厚度和横截面面积。 Giampiero Pampolini[20]利用一个模型耦合的非线性弹性和粘度的概念来描述一种聚氨酯泡沫进行单轴循环压缩响应。非弹性效应是由于泡沫的粘度特性，而应变局部化和滞后归因于应变能密度的非凸性。但该模型并没有重现第一次装卸周期的响应曲线。用现象学的损伤法来描述在第一个加载周期期间发生的损害。开发了一个识别材料常数的程序，并在两个单调和循环变形过程中粘度和损害相互作用进行了讨论。 5. 总结 实践证明，EPS有着隔热性好、自重轻、自立性强等优良性能，基于EPS的性能和优点，它已经被大量的推广和使用。虽然它具有非常多的特点以及可人的前景，但是我国对EPS本身的性能缺乏深入全面的研究,在使用过程中尚有许多需要解决的问题，EPS的回收等环保因素仍然值得人们考虑。在现代社会中，不论是国内还是国外，“白色污染”是人们共同关注的话题。一些国家已经在大力开展3R运动：要求做到废塑料的减量化、再利用、再循环。目前，在德、日、美等国家，对包装材料的回收处理十分重视，因此减少和合理使用EPS泡沫塑料是迫在眉睫的。因此本次毕业设计选择了对EPS泡沫塑料进行研究。EPS泡沫塑料的缓冲性能研究方法分析是本次毕业设计的研究方向。有理由相信，在相关科研人员的努力下，通过理论计算和实验验证等方法，对包装衬垫进行不断地优化，从而使EPS泡沫塑料的使用量逐步减少。 参考文献 [1] 熊志远. 聚苯乙烯泡沫（EPS）力学行为的实验研究[D]. 湘潭：湘潭大学，2007. [2] 颜志平. 泡沫聚苯乙烯(EPS)力学性能的室内试验研究[J]. 中南公路工程，2005，04：42～45. [3] 张卫兵. 聚苯乙烯泡沫(EPS)的特性及其在道路工程中的应用[J]. 公路，2004，05：146～149. [4] 李忠，郭丽. 聚苯乙烯泡沫(EPS)综述[J]. 四川建材，2012，05：10～11. [5] 汪萍. 可发性聚苯乙烯泡沫塑料板材的生产技术和应用[J]. 现代塑料加工应用，2001，02：31～34 [6] 佟富强，刘宝璋，杜金来. 形变速率对聚苯乙烯泡沫材料缓冲性能的影响[J]. 塑料包装，1996，04：32～34. [7] 阳以本. 可发性聚苯乙烯泡沫材料密度与压缩强度的关系及其应用[J]. 包装工程，1995，02：27～32. [8] 杜骋，杨军. 聚苯乙烯泡沫(EPS)的特性及应用分析[J]. 东南大学学报(自然科学版)，2001，03：138～142. [9] 毛快，张俊彦. 聚苯乙烯泡沫(EPS)压缩蠕变实验研究[J]. 材料导报，2007，S3：468～470. [10] 朱向荣，EPS动力特性试验研究[J]. 岩土工程学报. 2005，27（11）：1253～1256 [11] 卢富德，刘雄建，高德. 多胞缓冲材料本构模型与应用进展[J]. 浙江大学学报(工学版)，2014，07：1336～1344． [12] 高德，刘壮，孙智慧. 秸秆纤维EPS缓冲包装材料性能的研究[J]. 农业机械学报，2006，08：201～204. [13] 都学飞，欧阳效卓，张汪年. EPS缓冲材料的静态压缩性能的试验研究[J]. 包装工程，2012，03：39～41+45. [14] 龙志坚，龚维，肖开折，李宏，何颖. 微孔发泡聚苯乙烯力学性能的分析与探讨[J]. 现代塑料加工应用，2011，04：23～26. [15] Horvath J S. Expanded polystyrene( EPS) geofoam：an introduction to materialbehavior[J]. Geotextiles and Geomembrances，1994 (13)：26～280. [16] Dnskov M. Materials research on EPS20 and EPS15 under representative conditions inpavement Structures[J]. Geotextiles and Geomembrances，1997，(15)：17～18 [17] Kwang Young Jeong，Seong Sik Cheon，Mahbubul Basit Munshi. A constitutive model for polyurethane foam with strain rate sensitivity[J]. Journal of Mechanical Science and Technology，2012，267：. [18] Qunli Liu，Ghatu Subhash. A phenomenological constitutive model for foams under large deformations[J]. Polym Eng Sci，2004，443：. [19] James M. Gibert,Gregory S. Batt. Impact Oscillator Model for the Prediction of Dynamic Cushion Curves of Open Cell Foams[J]. Packag. Technol. Sci.，2015，283：. [20] Giampiero Pampolini，Michel Raous. Nonlinear elasticity，viscosity and damage in open-cell polymeric foams[J]. Archive of Applied Mechanics，2014，8412：. 开题报告 发泡聚苯乙烯缓冲性能分析方法 1. 研究意义 随着国民素质的提高，环保问题已经深入人心，“白色污染”也成为了人们共同关注的热门话题。发泡聚苯乙烯缓冲材料的质量比较轻省、具有良好的缓冲和保温性能，基于发泡聚苯乙烯缓冲材料的性能和优点，它已经被人们大量的推广和频繁的使用，但是发泡聚苯乙烯缓冲材料的回收问题值得人们认真考虑。滥用发泡聚苯乙烯缓冲材料，不仅是对泡沫原材料的浪费，更会造成危害极大的“白色污染”。目前，在美国等发达国家和其他欧洲城市，已经对废弃的包装泡沫塑料不断进行回收、循环利用，杜绝出现过度包装的现象，减少和合理使用发泡聚苯乙烯缓冲材料是一个迫在眉睫的问题。因此本次毕业设计选择了对发泡聚苯乙烯缓冲材料的缓冲性能进行研究，基于MATLAB软件，对发泡聚苯乙烯泡沫的缓冲性能进行理论分析，以此优化发泡聚苯乙烯泡沫缓冲衬垫的包装方案，减少缓冲衬垫的使用量，从而在一定程度上对“白色”污染起到了缓解作用。 2. 研究背景 程志胜等（2000）[1]等通过动态和静态压缩试验，分析了聚苯乙烯泡沫塑料衬垫的缓冲性能的基本特征，并应用弹塑性理论，建立了既能反映静态应力应变规律，又能反映其缓冲性能基本特征的非线性数学模型，同时根据试验数据识别了模型参数。 胡俊等（2015）[2]研究了四种不同密度的聚苯乙烯EPS泡沫材料单轴压缩下的应力-应变关系。在Gibson模型和Rusch模型基础上建立了EPS泡沫单轴压缩下应力-应变关系模型，并对建立的模型中各参数分别进行了定义;通过对能量吸收图、吸能效率图、理想吸能效率图的分析表明:应力在0.687MPa～1.038MPa范围内，密度为55kg/m3的EPS泡沫吸收的能量最大，当应力>1.038MPa时，吸收的能量随密度的增加而增加;四种密度的泡沫(28kg/m3, 40kg/m3, 55kg/m3, 70kg/m3)在吸能能力最佳时的应力分别为0.396MPa、0.565MPa、0.866MPa、1.222MPa;密度为55kg/m3的EPS泡沫最接近于理想吸能材料。 Changfeng Ge（2015）[3]研究了工业上使用的端盖、边缘或角落泡沫和试验样本的平面泡沫之间的应力-应变关系差异和动态影响。一个滞后周期压缩试验，在数字图像相关技术的结合，和击缓冲试验进行评估的静态应力应变、能量耗散和扁角泡沫最大加速度。该研究发现，在压缩试验期间角落泡沫发生的静态应力是平面泡沫的23%，这取决于压缩速度。并研究建议，传统的缓冲曲线向右偏移水平约23–35%为了包装设计师运用角泡沫。 卢富德等（2015）[4]为研究易损件在泡沫缓冲作用下冲击响应，建立了二自由度非线性产品包装系统模型，分析了易损件与泡沫黏弹性耦合后的跌落加速度特性，讨论破损边界上方区域内衬垫面积、厚度及其比值对易损件加速度响应的影响。结果表明，满足约束条件的易损件的最大加速度响应与缓冲衬垫面积和厚度关系成条状特性。随着缓冲衬垫面积与厚度增大，易损件加速度先减小后增大。易损件最大响应加速度与缓冲衬垫体积关系呈现勾状，最后引入系数表征缓冲衬垫的经济性和易损件安全性。 郭勇等[5]等基于缓冲材料宏观本构模型,建立不同缓冲材料组成的串联力学模型。介绍求解串联缓冲系统动力学响应的虚拟质量方法及包装结构优化设计方法。利用MATLAB/GUI编制用于缓冲系统结构设计的界面程序,并脱离MATLAB平台。为缓冲包装动力学设计提供直接方法。 佟富强等[6]等通过对聚苯乙烯泡沫材料进行不同形变速率的压缩实验，经过数学处理得到不同形变速率条件下材料的缓冲系数-最大应力曲线，找出压缩速度对材料缓冲性能影响的变化规律，可用以指导实际应用。 阳以本[7]介绍了我国现阶段普遍采用的螺旋传动，蒸汽加热发泡工艺生产的可发性聚苯乙烯泡沫材料密度与压缩强度(压缩50%)的关系，分析成果应用的可行性给其摩托车生产工厂带来好的经济效益(EPS应用于摩托车成品包装)。 杜骋等[8]在《聚苯乙烯泡沫(EPS)的特性及应用分析》文中也指出聚苯乙烯泡沫是一种性能优良的路基轻质填料，具有轻质、高强、较强的化学稳定性和水稳定性、良好的力学性能且施工方便简单等优点，在国外道路工程中有较为广泛的应用。我国对EPS的研究和应用较少，文章对EPS的物理化学性能、力学性能、EPS作为路基轻质填料的结构设计方法、EPS在道路工程中的应用等方而作了较为全而的介绍和分析，对我国使用EPS有借鉴作用。 毛快[9]主要研究了EPS的压缩蠕变性能。EPS材料作为一种热塑性材料，其特点决定其在使用过程中与其他黏弹性材料一样会发生蠕变和松弛现象，因此要考虑EPS材料的性能随时间的变化规律。因为EPS作承重构件时通常受压，长期蠕变则需要几天或更长的时间，且发现仅当应力超过阀值水平时才产生。针对不同密度的EPS进行了短期压缩蠕变实验研究，讨论应力水平和密度对蠕变的影响。 朱向荣[10]采用三轴试验研究EPS在不同加载频率、不同循环次数和不同围压的下强度和模量等的变化规律，为EPS工程应用提供了一些理论依据。 卢富德等[11]针对泡沫结构的压缩性能，介绍了泡沫结构的试验压缩响应与数值模拟方面的研究进展。结果表明：相对密度、环境温度、压缩应变率及微观结构等参数对泡沫的缓冲性能影响显著。相对密度越大，泡沫结构的屈服应力越大，吸收能量的能力越大；泡沫的吸收能量能力一般随环境温度的增加而减小；由于泡沫的基体材料表现率相关性，应变率增加导致吸收能量能力的增加；当微观结构不同时，泡沫的细胞分布导致结构的缓冲性能差异显著。 高德等[12]根据植物秸秆纤维聚苯乙烯缓冲材料静态压缩实验数据，研究了材料的缓冲性能，并考虑植物秸秆纤维的影响，对Sherwood- Frost本构关系框架进行了扩充，建立了非线性本构关系模型，并利用实验数据成功识别模型参数此种描述植物纤维类材料非线性力学行为的方法，为进一步研究和开发植物纤维聚苯乙烯材料提供了理论基础。 都学飞等[13]比较分析了４种厚度EPS缓冲包装材料的压缩变形回复性、外力－位移曲线等性能。结果表明，EPS缓冲材料的回复性、永久变形与材料的厚度有很大关系，而它的外力－位移曲线的走向大致相同，厚度越大吸收的能量越多，缓冲性能越好。 龙志坚[14]制备了发泡量相同的微发泡聚苯乙烯, 采用扫描电镜(SEM)和Image-pro图像处理软件对微发泡聚苯乙烯的微孔尺寸进行了观察和统计;建立微球模型, 分析了微孔尺寸大小对微发泡聚苯乙烯力学性能的影响。结果表明:细小而均匀的泡孔对微发泡聚苯乙烯力学性能的提高有较明显的促进作用, 微球模型的计算结果与宏观力学性能的影响规律有很好的重现性。 Horvath[15]用边长为5cm的EPS立方体试件在应变速率10 mm/min的条件下，采用应变控制形式进行了无侧限单轴压缩试验，得到压缩应力-应变曲线并对该曲线进行了分析。但没有对多种密度与多种加载速率的情况进行试验比较和分析。 Duskov[16]采用直径为10 cm、高为20 cm的圆柱体EPS试件在20 kPa的作用下进行蠕变研究，得出蠕变曲线，从蠕变曲线的分析中可以看出：EPS材料的蠕变主要发生在加载初期，随着加载龄期的增长，蠕变的速率趋于稳定，在加载1年以后蠕变的速率几乎接近常数。 Kwang Young Jeong[17]研究了聚氨酯泡沫塑料的应变率相关行为并制定了新的本构模型，以提高在各种应变速率中实验数据的拟合性。该模型的七个参数被两种应变速率下进行的准静态压缩试验所决定。压缩试验两种应变速率。高、低密度聚氨酯泡沫塑料的两种模型显示在了不同的应变率下应力-应变的关系。所进行的动态压缩试验得出了在高应变率下应力-应变的数据并且将此结果同本构模型进行了比较。 Qunli Liu[18]认为受到大的变形建议。当加压力负载时，可在五参数模型完全确定应力-应变响应的三个基本特征，即线性、可塑性状应力平台和致密化的阶段。此外，根据不同的屈服强度和硬化状或软化的约束条件，该模型的参数可以系统变化，以确定泡沫的初始密度。 James M. Gibert[19]探讨了压盘影响开孔泡沫缓冲材料动力学的定性认识。一个超弹性材料模型是用来描述阻尼和滞回特性为线性粘弹性的非线性泡沫的应力–应变关系。使用一个简单的非线性不连续模型的跌落试验以及数值模拟，研究探讨对其的物理影响。数值研究表明，该模型能够提供预测的冲击脉冲的形状、持续时间和幅度，在不同的静态应力和下降高度。模型所产生的动态缓冲曲线保留凹向上的“槽”的实验曲线的形状特征。此外，该模型表明，给定的下降条件下，冲击吸收的最佳振幅取决于泡沫的厚度和横截面面积。 Giampiero Pampolini[20]利用一个模型耦合的非线性弹性和粘度的概念来描述一种聚氨酯泡沫进行单轴循环压缩响应。非弹性效应是由于泡沫的粘度特性，而应变局部化和滞后归因于应变能密度的非凸性。但该模型并没有重现第一次装卸周期的响应曲线。用现象学的损伤法来描述在第一个加载周期期间发生的损害。开发了一个准确的识别材料常数的程序，并在两个单调和循环变形过程中粘度和损害的相互作用进行了讨论。 3. 主要研究工作 3.1 研究内容 本文围绕发泡聚苯乙烯泡沫（EPS）的缓冲性能进行系统研究，计算在不同的缓冲面积、不同的跌落高度、不同的衬垫厚度等情况下，经过聚苯乙烯泡沫缓冲作用下，易损件的加速度响应。所有研究内容将在理论分析的前提和试验的基础下撰写。本论文的研究内容大致为： 第一章：绪论：介绍发泡聚苯乙烯泡沫（EPS）的缓冲性能、国内外研究背景、意义以及思路等。 第二章：利用MATLAB/GUI界面建立EPS缓冲数据库求解缓冲曲线：建立MATLAB/GUI界面，基于图片的像素值，求解EPS最大加速度-静应力等缓冲曲线 第三章：利用MATLAB/GUI求解EPS缓冲作用下单自由度缓冲性能：主要通过包装件地脆值、包装件的质量、跌落高度，以及EPS的应力-应变曲线来设计缓冲衬垫的缓冲面积和衬垫厚度。 第四章：利用MATLAB/GUI求解二自由度EPS缓冲性能：建立EPS的缓冲本构模型，计算绘制在不同的发泡聚苯乙烯泡沫（EPS）的缓冲作用下，易损件的加速度响应曲线。 第五章：总结与展望 3.2 研究重点与难点 3.2.1 编辑求解EPS缓冲曲线程序 利用MATLAB程序在曲线本身坐标轴和EPS材料应力-应变曲线照片的像素坐标（以图片左下角像素点，水平方向为x轴，竖直方向为y轴）建立一系列的转换关系，通过曲线点的拾取，可获得精确的坐标值，并通过相关公式对其加以处理，从而获得缓冲系数-最大应力曲线和最大加速度-静应力曲线。 3.2.2 编辑求解EPS单自由度缓冲作用下衬垫设计 在计算缓冲衬垫的面积、衬垫厚度等设计程序的条件下，明确求解缓冲衬垫的面积、衬垫厚度等计算步骤。将这些求解的步骤编辑成程序，用MATLAB/GUI制作成程序界面。将编辑完成的程序和界面生成相应的软件包，用于设计缓冲衬垫的最优解。只要选择EPS材料的密度，设定包装件地脆值、包装件的质量、跌落高度，便可求解缓冲衬垫的面积和厚度。 3.2.3 编辑求解EPS二自由度缓冲作用下易损件的加速度响应程序 编辑程序利用MATLAB/GUI建立工作界面，只要输入发泡聚苯乙烯泡沫（EPS）衬垫的厚度、缓冲面积，包装件的质量、跌落高度，便可以计算并绘制出在制定的EPS缓冲衬垫作用下易损件的加速度响应曲线。 4. 研究方法与实施进度计划 对本课题的研究主要采用理论研究和实际相结合的方法，在指导老师的指导下对课题进行研究。 在理论研究阶段，通过各种途径收集中、外文献资料，对前人的研究成果进行阅读、总结、归纳、和分析，为文章的撰写提供大量客观可靠的事实依据。 在实际解决问题阶段，总结前人的研究结果，利用MATLAB/GUI软件进行编程，生成应用文件，对EPS泡沫缓冲性能进行计算，并通过试验进行验证。 最后掌握理论和实际资料，结合大学期间所学的专业知识进行论文的写作。 为了有序的完成论文各部分的写作，笔者制订了如下的进度计划： 第一阶段：2015年10月——2015年11月，阅读并收集中、外文献资料。 第二阶段：2015年11月——2015年12月，完成外文翻译，撰写文献综述和开题报告初稿。 第三阶段：2015年12月——2016年1月，对开题报告和文献综述初稿进行修改完善。 第四阶段：2016年1月——2016年2月，建立MATLAB/GUI界面，建立EPS缓冲数据库求解缓冲曲线。 第五阶段：2016年2月——2016年3月，建立MATLAB/GUI界面，求解EPS缓冲作用下单自由度缓冲性能和二自由度EPS缓冲性能。 第六阶段：2016年3月——2016年4月，撰写论文。 第七阶段：2016年4月——2016年5月，修改论文，准备答辩。 参考文献 [1] 熊志远. 聚苯乙烯泡沫（EPS）力学行为的实验研究[D]. 湘潭：湘潭大学，2007. [2] 颜志平. 泡沫聚苯乙烯(EPS)力学性能的室内试验研究[J]. 中南公路工程，2005，04：42～45. [3] 张卫兵. 聚苯乙烯泡沫(EPS)的特性及其在道路工程中的应用[J]. 公路，2004，05：146～149. [4] 李忠，郭丽. 聚苯乙烯泡沫(EPS)综述[J]. 四川建材，2012，05：10～11. [5] 汪萍. 可发性聚苯乙烯泡沫塑料板材的生产技术和应用[J]. 现代塑料加工应用，2001，02：31～34 [6] 佟富强，刘宝璋，杜金来. 形变速率对聚苯乙烯泡沫材料缓冲性能的影响[J]. 塑料包装，1996，04：32～34. [7] 阳以本. 可发性聚苯乙烯泡沫材料密度与压缩强度的关系及其应用[J]. 包装工程，1995，02：27～32. [8] 杜骋，杨军. 聚苯乙烯泡沫(EPS)的特性及应用分析[J]. 东南大学学报(自然科学版)，2001，03：138～142. [9] 毛快，张俊彦. 聚苯乙烯泡沫(EPS)压缩蠕变实验研究[J]. 材料导报，2007，S3：468～470. [10] 朱向荣，EPS动力特性试验研究[J]. 岩土工程学报. 2005，27（11）：1253～1256 [11] 卢富德，刘雄建，高德. 多胞缓冲材料本构模型与应用进展[J]. 浙江大学学报(工学版)，2014，07：1336～1344． [12] 高德，刘壮，孙智慧. 秸秆纤维EPS缓冲包装材料性能的研究[J]. 农业机械学报，2006，08：201～204. [13] 都学飞，欧阳效卓，张汪年. EPS缓冲材料的静态压缩性能的试验研究[J]. 包装工程，2012，03：39～41+45. [14] 龙志坚，龚维，肖开折，李宏，何颖. 微孔发泡聚苯乙烯力学性能的分析与探讨[J]. 现代塑料加工应用，2011，04：23～26. [15] Horvath J S. Expanded polystyrene( EPS) geofoam：an introduction to materialbehavior[J]. Geotextiles and Geomembrances，1994 (13)：26～280. [16] Dnskov M. Materials research on EPS20 and EPS15 under representative conditions inpavement Structures[J]. Geotextiles and Geomembrances，1997，(15)：17～18 [17] Kwang Young Jeong，Seong Sik Cheon，Mahbubul Basit Munshi. A constitutive model for polyurethane foam with strain rate sensitivity[J]. Journal of Mechanical Science and Technology，2012，267：. [18] Qunli Liu，Ghatu Subhash. A phenomenological constitutive model for foams under large deformations[J]. Polym Eng Sci，2004，443：. [19] James M. Gibert,Gregory S. Batt. Impact Oscillator Model for the Prediction of Dynamic Cushion Curves of Open Cell Foams[J]. Packag. Technol. Sci.，2015，283：. [20] Giampiero Pampolini，Michel Raous. Nonlinear elasticity，viscosity and damage in open-cell polymeric foams[J]. Archive of Applied Mechanics，2014，8412：. 外文翻译和原稿 聚氨酯泡沫应变率敏感性的本构模型 A constitutive model for polyurethane foam with strain rate sensitivity 摘 要： 本工作研究了聚氨酯泡沫塑料的应变率相关行为并制定了新的本构模型，以提高在各种应变速率中实验数据的拟合性。该模型的七个参数被两种应变速率下进行的准静态压缩试验所决定。压缩试验两种应变速率。高、低密度聚氨酯泡沫塑料的两种模型显示在了不同的应变率下应力-应变的关系。所进行的动态压缩试验得出了在高应变率下应力-应变的数据并且将此结果同本构模型进行了比较。 关键词：聚氨酯泡沫；本构模型；动载荷；应变率 1．引言 对于汽车耐撞性能的改进不能被过高估计。美国交通部门估计从1999年到2009年，有超过400,000人遇难，2000万受伤需要住院治疗。这个情况，连同一系列的环境问题和社会压力在法律的支持下已经领导，并将继续引导，高度创新的设计，涉及到有色金属合金、智能结构、复合材料等先进的材料和泡沫。这项研究特别感兴趣的是耐撞性结构泡沫塑料材料的使用。聚合物泡沫目前被用作在保险杠和作为填充材料加固屋顶和门梁，并且它们的应用程序将被不断扩展。它们可加强薄弱部位，从而可以对冲击载荷做出有效反馈，即加强耐撞性。泡沫的能量吸收性能是当保持一个几乎恒定应力值时所经历的变形能力。在其他聚合物中，聚氨酯泡沫被广泛用在工业上。 泡沫已经成为了无数试验、数值和理论调查的主题。Shim等人对硬质聚氨酯泡沫塑料进行了速度从2—4m/s 的普通冲击试验。Avalle等人研究了使用能量吸收图的聚合物泡沫的能量吸收特性。Rusch提出了压应力与应变之间的关系作为初始压缩模量的泡沫和应变的函数。Meinecke和Schwaber描述作为一个应变函数的初始系数并改变了作为一个幂级数的应变的应变函数形式。Sherwood和Frost修改了聚氨酯泡沫塑料在受压状态下的本构模型的形状函数。Langseth进行了用挤压状态下的铝在在不同速率下持续冲击的冲击试验。 2. 本构方程 2.1 准静态模型 聚合物泡沫的应力应变曲线，在一般情况下显示出三个区域：一个线性弹性区域，其次是可塑性平坦区域以及应力陡峭上升的致密区。吉布森模型的三个方程描述了每个区域，如下所示 当 时， QUOTE (2.1) 当 QUOTE (2.2) QUOTE 时， 当 QUOTE QUOTE 时， (2.3) 和 分别是工程应力和工程应变。该模型具有五个参数: E为杨氏模量的弹性部分。 QUOTE QUOTE 为屈服应力。为致密化应变和两个常量D、m。这模型在平坦区域应力有常量值的限制并且在两个区域的边界应力应变曲线是不光滑的。Rusch 模型可以描述的两个力的总和： QUOTE (2.4) 这里的四个参数，a、b、p和q，可以根据经验确定。虽然这个模型通过一个方程描述了应力-应变关系，它的缺点在于不能精确描述致密化的阶段。为更好适用泡沫材料在实验中的应力应变关系。刘和苏巴斯建议式 (2.5) 所示的模型。这种模式具有六个参数。设置B作为一个单元，把它变成一个五参数模型。 QUOTE (2.5) 最近， 一个模型被Avalle 等人提出。如下所示。 QUOTE (2.6) 该模型还具有五个参数。参数E、A和B与密度相关而m和n不是。第一项适合弹性区和平坦区，而第二项适合致密区。虽然这个模型很好的描述了聚合物泡沫的应力-应变曲线，但它没有包括因动态加载下对于应变率的影响。 2.2 包括应变率效应的模型 包括应变率效应的模型由Nage等人提出。 (2.7) 这里的 是工程压应力， 是应变的形状函数，代表了在参考应变速率和 下应力-应变的关系。应变速率的功能，是一种应变和应变率的函数，在参考应变率时的单位价值。Sherwood和Frost所描述的形状函数的十分之一阶多项式函数的应变，其系数是确定从压缩试验的参考应变率。应变率函数的假设，即模型图的应力与应变率的假设，表现出的线性函数的应变可以表示如下： QUOTE (2.8) 双参数，a和b，可以发现通过线性近似的实验的应力-应变率数据的几个应变率。由于描述Eqs的形函数。（2.7）和（2.8）只不过是一个应力-应变方程，它取代模型（6）的五个参数比使用十分之一阶多项式更有效。应变率的函数方程（2.9）是基于在一定的应变-应力数据曲线情况下具有线性应变函数的假设。仔细检查实验数据表明，这是不正确的，特别是对于半柔性聚氨酯。引入了一个新的应变率函数来描述应力与应变率的对数曲线图，该曲线是凹的，在恒定应变下： (2.10) 通过替换式的形状函数和应变速率的公式（2.6）、函数（2.7）和式（2.10），分别为应力-应变方程，包括应变率效应推导： QUOTE QUOTE (2.11) 该模型方程（2.11）有七个参数。这些参数的确定分为2个步骤。首先，在参考应变率下，五个参数，E、A、B、m和n，确定通过拟合实验的应力-应变数据。然后，剩余的2个参数，a和b，是通过线性近似的实验的应力-应变数据的另一个应变率被发现的。大多数的参数被认为是依赖于材料的密度和温度。因此，两步实验应该在同一密度和温度下进行。 3．实验与参数辨识 3.1 低密度聚氨酯泡沫塑料模型 准静态测试是在位移测量装置和力传感器安装形成的室温下，使用具有100kN大负荷容量的MTS 810机器进行的。密度为67kg/m3的开孔式圆柱聚氨酯泡沫（如图3.1所示，由韩国lacomtech有限公司提供）是在室温下测试。 图3.1 聚氨酯泡沫塑料试样 试样的直径和长度分别为42毫米和40毫米，分别为。准静态试验在0.001s-1的应变速率进行了应力应变实验数据。五个参数方程（2.6）采用最小二乘法进行数据拟合。参数的计算会被拟合曲线点的数目所影响。表1显示的参数，是通过拟合33个应力-应变实验数据点得到的。实验和装配的应力-应变曲线，如图2所示。拟合曲线很接近实验曲线。（2.11）式中a和b，这两个参数的确定，在另一个应变率的实验数据是必要的。这个实验用密度同为67kg/m3的聚氨酯泡沫在应变速率为0.1s-1的情况下进行的。因为在应变速率为0.001s-1时采用最小二乘法得到五个参数，所以a和b两个参数是由方程的本构方程拟合实验数据确定。表3.2中的数据通过最小二乘法拟合实验数据计算得出的参数。如图3.3所示，应力-应变曲线很接近实验数据。 表3.1 从曲线拟合计算的五个参数 A B E m n 0.3863 0.3151 13.9996 3.4162 1.3680 表3.2 从曲线拟合计算的二个参数 a b 0.0401 0.1222 图3.2 应力-应变关系 图3.3 应力-应变关系 ( =67kg/m3， =0.001s-1 ) ( =67kg/m3， =0.1s-1 ) 获得构成所有参数接近的亲双准静态实验的本构模型，所得到的应力-应变曲线可以得到任何应变率。如图3.4所示，0.001、0.01、0.1、1、10和100s-1 这六个应变率的值分别代入式（2.11）得到六条应力-应变曲线。其中应变率为0.001s-1 对应最低的曲线。 对于在相同的应变下，当泡沫的应变速率更高时，更多的能量将被储存在泡沫中。用本构方程拟合得到七个参数，用本构方程拟合得到了67kg/m3的聚氨酯泡沫体在动态载荷作用下的应力应变关系。经验证，方程可以采用动态压缩试验。 图3.4 应变率的应力-应变关系 图3.5 应力-应变关系 ( =67kg/m3 ) ( =67kg/m3， =97s-1 ) 落塔型冲击试验机（英斯特朗Dynatup 9250 HV）进行动态测试。本机可提高离散变量到一个特定的高度，下降到试样使用引力或加速弹簧。当26.5kg的重物冲击试样，测量的冲击速度为3.88m/s，并计算应变速率为97s-1。图3.5中的实心线显示了实验的应力-应变关系，虚线描述了本构模型的关系。虽然它不遵循的动态效果的振荡，但是该模型接近的实验结果。 3.2 高密度聚氨酯泡沫塑料 图3.6 应力-应变关系 图3.7 应力-应变关系 ( =89kg/m3， =0.001s-1 ) ( =89kg/m3， =0.1s-1 ) 进行了89kg/m3的高密度聚氨酯泡沫塑料的准静态测试。第一，在应变速率为0.001s-1应力-应变的实验数据通过采用最小二乘法拟合得到方程（2.6）的五个参数。实验和拟合应力-应变曲线如图3.6所示。 拟合曲线很接近实验曲线。另一个静态应变率的实验数据都需要找到两参数方程a和b（2.11）。在应变速率为0.1s-1 QUOTE 的准静态试验，通过曲线拟合得到的参数。如图3.7所示，应力-应变曲线很接近实验数据。 由于所有的本构模型的参数已被发现，应力-应变曲线可以得出任何应变率。将0.001、0.01、0.1、1、10和100s-1 QUOTE 这六个应变率分别代入式（2.11），在图3.8中绘制成应力-应变曲线。 要确定该模型是否能在高应变率和动态压缩试验中工作。当重物冲击试样时，撞击速度为3.88m/s，计算为应变速率为97s-1。实验绘制的应力-应变曲线与所提出的本构模型进行比较。图3.9中的实心线显示了实验的应力-应变关系，虚线描述了本构模型的关系。 图3.8 应变率的应力-应变关系 图3.9应力-应变关系 ( =89kg/m3) ( =89kg/m3， =97s-1 ) 在非常低的应变率的情况下，由于本构模型的参数是由2个准静态实验得到，这种模型不能跟随动态效果的振荡实验。图3.5和图3.9的实验结果表明，试样的密度越高，进行的的振荡越大。 3.3 应力与应变率 得到模型的所有参数后，可以绘制在不同规模水平下参数记录的应变率。如图3.10所示，67 kg/m3的低密度聚氨酯泡沫。如图3.11所示，89 kg/m3的高密度泡沫的应力-应变率。每个曲线几乎呈线性关系，但有点凸。Nage et al提出的本构模型，如式（2.7）通过假设的应力与应变率曲线表现出应变的线性函数。他们的实验的日志记录图一般遵循直线，这往往遵循的提出模型的特点，但凸点最多的是泡沫材料。 图3.10恒定应变下的应力-应变率 图3.11恒定应变下的应力-应变率 ( =67kg/m3) ( =89kg/m3) 4. 总结 在设计吸能结构泡沫时，需要对泡沫的力学特性进行数学模型描述。几个模型以及描述的应力-应变关系的泡沫是工作在一个特定的应变率。本文提出了一种用于聚氨酯泡沫的本构模型，以描述的应力-应变关系在很宽的范围内的应变率。 该模型有七个参数。通过两步优化程序的参数进行评估。（1）第一步是参考应变率和发现五个参数。（2）在下一步中，通过将所提出的模型和实验测得的值之间的曲线拟合，发现在这步中的应变率的大小影响的剩余2个参数。 在初始应变率为97s-1时，动态压缩试验被提供应力-应变下进行，这些数据绘制出他们所提出的模型。应力-应变曲线的提出模型以及如下的实验数据，即使它并没有准确地描述振荡，是因为系统所造成的冲击。 参考文献 [1] US Department of Transportation, 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 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 Journal 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. 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 simulation, 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 experimental tests, International Journal of Impact Engineering , 34 (2007) 3-27. [12] Q. Liu, G. Subhash and X. L. Gao, A parametric study on crushability of open-cell structural polymeric foams, Journal of Porous Materials , 12 (2005) 233-248. [13] Q. Liu and G. Subhash, A phenomenological constitutive model for foams under large deformation, Polymer Engineering and Science , 44 (5) (2004) 463-473. [14] A. Nagy, W. L. Ko and U. S. Lindholm, Mechanical behavior of foamed materials under dynamic compression, Journal of Cellular Plastics , 10 (1974) 127-134. 摘 要 缓冲系统由产品-缓冲衬垫组成。典型破损认为只要产品加速度达到极限，即为产品失效，因此设计缓冲材料常用的 和 。但在实际分析中，产品的破坏不但与 有关，还与脉宽及脉冲的形状有关，即需要考虑产品的易损件，这时的系统应为易损件-质量主体-缓冲衬垫组成。由于EPS具有良好的缓冲性能，本文利用MATLAB/GUI，对两类形式的缓冲系统进行响应求解与设计。本课题研究内容为： （1）在产品-EPS衬垫模型，首先获取所发表文献中的有关EPS的 数据，并汇总不同密度的 曲线，构建EPS材料库。 （2）在易损件-质量主体-EPS衬垫模型，首先得到EPS的本构模型，然后得到EPS传递给产品质量主体的加速度离散脉冲，另外根据三次样条插值计算得到，脉冲激励下对易损件系统的冲击响应。 本课题基于MATLAB/GUI，设计了一系列简单实用的应用模块，使用者只需要在MATLAB/GUI界面中输入和选择所需要的参数，便能完成实现EPS缓冲性能、缓冲衬垫设计的输出。 关键词： EPS；MATLAB／GUI缓冲衬垫；最大加速度-静应力曲线；易损件 Abstract The buffer system is composed of product quality - buffer pad. The typical damage is as long as the acceleration of products to reach the limit that is product failure. Therefore, characterization of cushioning materials commonly used cushion coefficient-maximum stress curve, and the maximum acceleration static stress curve. But in actual analysis, the destruction of products not only with a maximum acceleration of products related, but also pulse width and pulse shape, that is, to consider vulnerable parts inside the product, then the system should be composed of three parts as vulnerable parts - quality subject - buffer pad. Because EPS has good cushioning property, this paper uses the MATLAB/GUI tool to solve the two kinds of the buffer system. The content of this research is: (1) In the product quality of EPS cushion model, we first get the published literature related to EPS stress-strain data, summarize the stress-strain curves of different densities, construct EPS materials library. (2) in the wearing parts: main quality of EPS cushion model, EPS has firstly obtained the constitutive model, then get EPS transfer to the quality of products the main acceleration discrete pulses, also based on the cubic spline interpolation calculated pulse excitation the damageable parts of shock system response. This paper based on the MATLAB software, a series of simple and practical application module, users only need to in Matlab / GUI interface to input and select the required parameters, can realize the output of EPS cushion performance, the cushioning design. Keywords： EPS；MATLAB / GUI buffer pad；Maximum acceleration static stress curve；Vulnerable parts 目 录 I 摘 要 IIAbstract 11. 绪 论 11.1 研究背景和意义 11.2 发泡聚苯乙烯泡沫的介绍 11.2.1 发泡聚苯乙烯泡沫的定义 11.2.2 发泡聚苯乙烯泡沫的特性 21.3 国内外研究现状 21.3.1 国内研究现状 31.3.2 国外研究现状 41.4 研究内容 52 利用MATLAB/GUI界面建立EPS缓冲数据库求解缓冲曲线 52.1 建立EPS缓冲性能数据库 52.1.1 读取所需处理的应力-应变曲线 52.1.2 寻找图像坐标刻度与图片的像素关系 62.1.3 将图像曲线上选取的点转化为像素坐标系上的点坐标值 62.1.4 建立新图像坐标，将选取的点的像素坐标值转化为真实的试验数值 72.2 利用应力-应变数据计算并绘制缓冲曲线 72.2.1 读取MATLAB的工作文件夹中的文本文件 2.2.2 建立新图形窗口，绘制EPS的曲线图 7 2.2.3 绘制成曲线图 7 2.2.4 绘制曲线图 8 92.3 EPS缓冲数据库界面程序制作 102.4示 例 113 利用MATLAB/GUI求解EPS缓冲作用下单自由度缓冲性能 113.1 求解EPS缓冲作用下单自由度缓冲性能程序设计 113.1.1 设计输入单自由度的缓冲环境的程序 123.1.2 计算缓冲衬垫的最小缓冲面积 123.1.3 计算缓冲衬垫的厚度 133.2 EPS缓冲作用下单自由度缓冲性能界面程序制作 133.4 示例 154 利用MATLAB/GUI求解二自由度EPS缓冲性能 154.1 包装产品易损件的EPS本构模型 174.2 缓冲过程分析以及运行程序设计 174.2.1 EPS衬垫压缩阶段缓冲分析 184.2.2 EPS衬垫回弹阶段缓冲分析 194.2.3 包装产品易损件加速度响应分析 204.3 EPS缓冲作用下易损件的加速度响应曲线界面程序制作 204.4 示例 225. 总结与展望 23参考文献 25附 录 25附录1 建立EPS缓冲性能数据库程序 27附录2 利用应力-应变数据计算并绘制缓冲曲线程序 30附录3 求解单自由度EPS缓冲衬垫设计程序 34附录4 求解二自由度EPS缓冲作用下缓冲性能程序 39致 谢 1. 绪 论 1.1 研究背景和意义 随着国民素质的提高，环保问题已经深入人心，“白色污染”也成为了人们共同关注的热门话题。EPS的质量比较轻省、具有良好的缓冲和保温性能，基于EPS材料的性能和优点，它已经被人们大量的推广和频繁的使用，但是EPS材料的回收问题值得人们认真考虑[1]。滥用EPS材料，不仅是对泡沫原材料的浪费，更会造成危害极大的“白色污染”。目前，在美国等发达国家和其他欧洲城市，已经对废弃的包装泡沫塑料进行回收、循环利用，杜绝出现过度包装的现象。减少和合理使用发泡聚苯乙烯缓冲材料是一个迫在眉睫的问题，因此本次毕业设计选择了对发泡聚苯乙烯缓冲材料的缓冲性能进行研究，基于MATLAB软件，对发泡聚苯乙烯泡沫的缓冲性能进行理论分析，以此优化发泡聚苯乙烯泡沫缓冲衬垫的包装方案，减少缓冲衬垫的使用量，从而在一定程度上对“白色”污染起到了缓解作用。 1.2 发泡聚苯乙烯泡沫的介绍 1.2.1 发泡聚苯乙烯泡沫的定义 聚苯乙烯是在催化剂的作用下，使苯乙烯单体发生了加聚反应而生成的一种高分子聚合物[2] REF _Ref450387863 \n \h \* MERGEFORMAT 。 聚苯乙烯泡沫(英文Expanded polystyrene，简称为EPS)是由聚苯乙烯(简称PS)颗粒在加热的作用下，使内部的戊烷发生汽化，膨胀发泡后形成具有密闭空腔结构的硬质泡沫塑料[3]。根据泡沫生产发泡的过程分类，一般有模具内部发泡和挤出发泡两种方式。 1.2.2 发泡聚苯乙烯泡沫的特性 （1）密度和力学特性 EPS的各项力学特性都是随着它的密度的增大而提高, 因此EPS的密度是一个十分重要的指标，表观的密度取决于PS颗粒的发泡膨胀的倍数大小，一般情况下，是每升体积在10至45克之间[4] REF _Ref450387913 \n \h \* MERGEFORMAT 。 （2）保温性能 EPS材料中的每一个小珠粒是单独密闭的一个中空结构，使它具有许多其它材料所不具备的保温特性。空气本身就是一种热的不良导体，如果将空气留在EPS小珠粒的空腔结构之中，在内部的空气就不容易甚至不能发生对流，从而使这种材料具有特别显著的保温隔热效果。 （3）容易产生静电 EPS材料的导电性能较差，然而在EPS泡沫的生产和运输的过程中，材料和材料之间会相互发生不同程度的摩擦，容易使EPS泡沫自身产生静电。所以，EPS泡沫生产厂家会在制造的过程中加入一些具有防静电效果的药剂，用来提高聚苯乙烯泡沫的导电能力。 1.3 国内外研究现状 1.3.1 国内研究现状 佟富强[5]等通过前期实验和后期数据处理对EPS缓冲衬垫进行了研究和分析，拟合出了在不同的压缩速率下EPS缓冲衬垫的 QUOTE 曲线，明确了压缩的速率对EPS缓冲性能的影响。 毛快[6]等主要研究了不同密度的发泡聚苯乙烯缓冲材料在较短时间内的压缩状态下，产生蠕变和应力松弛的影响。并且发泡聚苯乙烯发生的弹性蠕变和应力松弛也与实验时间有一定的关系。 朱向荣[7]等研究了不同的受压缩速率和加载次数对发泡聚苯乙烯材料缓冲性能的影响，如表1-1所示。 表1-1 不同的加载次数对EPS的缓冲性能的影响 循环次数/次 屈服应变 屈服应力/kPa 弹性模量/MPa 10 3.97 70.77 1.68 100 2.47 104.68 4.33 500 2.28 81.67 3.27 5000 3.59 69.77 2.12 10000 2.10 67.62 3.11 卢富德[8]等根据发泡聚苯乙烯的动态压缩实验和数据模拟等方法，对发泡聚苯乙烯的物理压缩性能进行了分析，如图1.1所示。实验和模拟得出结论：发泡聚苯乙烯的缓冲性能主要由它的密度、内部结构、所处环境等其他方面决定的。发泡聚苯乙烯的缓冲性能与相对密度呈正比，但与周围环境温度呈反比。在不同内部结构的情况下，发泡聚苯乙烯的缓冲性能也截然不同。 图1.1 典型泡沫的应力－应变曲线 高德[9]等研究了秸秆纤维和EPS复合材料，建立了该材料的力学模型，通过实验和数据拟合较完整地表述了这种衬垫的非线性的力学行为，为进一步开发这种缓冲材料提供了切实的理论依据。 都学飞[10]等研究了不同厚度的EPS缓冲材料的缓冲力学性能。研究显示：EPS缓冲材料的力学性能与材料的厚度有很大关系:缓冲衬垫的厚度与它具有的缓冲能力呈正比，如图1.2所示。 图1.2 厚度对EPS永久压缩率的影响 1.3.2 国外研究现状 Horvath[11]用长、宽、高均为5厘米的发泡聚苯乙烯作为试验样件，在应变速率为0.167 mm/s的状态下进行了压缩泡沫的测试，得到了压缩 QUOTE 曲线并对该曲线进行了分析。但没有对多种密度与多种加载速率的缓冲情况进行试验比较和分析。 Duskov[12]采用直径为0.1m、高为0.2m的圆柱发泡聚苯乙烯作为试验样件，在20 kPa的作用下进行了蠕变的测试实验，得出该材料相应的测试曲线，根据曲线得出了结论：聚苯乙烯缓冲材料的蠕变一般发生在压力加载的初期，而蠕变的速率随着压力加载时间的增长，逐渐稳定下来，甚至在加载1年以后蠕变速率几乎接近于某个常数。 Qunli Liu[13]认为泡沫在受到较大变形的情况下，增加压力负载时，可用五参数模型完全确定 响应的三个基本特征，即线性阶段、可塑性状屈服平台阶段和致密化的阶段。此外，根据不同的屈服强度和硬化或软化的约束条件，该模型的参数可以系统变化，以确定泡沫的初始密度。 1.4 研究内容 发泡聚苯乙烯泡沫材料在包装设计以及工程等其他领域起着十分重要的作用，但是由于它难以降解等其它缺点，对人类的生存环境造成了严重的影响。所以发泡聚苯乙烯缓冲材料的包装优化设计备受学术界的关注。在相同缓冲效果的情况下，保证包装产品不受损坏，如何减少发泡聚苯乙烯的使用量已经成为了一个研究热点，这不但减少了发泡聚苯乙烯泡沫材料的浪费，并解决了部分过度包装问题，更在一定程度上缓解了废弃泡沫对环境污染。 本课题基于MATLAB/GUI软件编程，从使用者的角度出发，设计了一系列简单而又实用的应用模块，主要应用于：不同密度EPS的缓冲性能的收集、单自由度EPS缓冲衬垫的设计、二自由度EPS缓冲作用下包装产品易损件加速度响应的计算。使用者只需要在MATLAB/GUI界面中输入和选择所需要的参数，设定所需要的缓冲环境，而不需要修改源程序的代码，便能完成实现发泡聚苯乙烯缓冲性能的获取、单自由度缓冲衬垫的设计、二自由度包装产品易损件加速度响应计算的输出。 2 利用MATLAB/GUI界面建立EPS缓冲数据库求解缓冲曲线 在包装行业使用发泡聚苯乙烯泡沫对产品进行包装设计时，需要提供发泡聚苯乙烯泡沫的缓冲曲线和数据。而在学术文献里获取的发泡聚苯乙烯泡沫的 曲线，若使用普通方法，很难获取精确的坐标数值，导致无法将该曲线转化成其他缓冲曲线，如 [14] QUOTE 曲线，建立不同密度的EPS缓冲性能数据库 曲线和 QUOTE 曲线等，这无疑浪费了学术文献中发泡聚苯乙烯泡沫的研究资源。如果利用MATLAB程序在曲线本身坐标轴和图片像素坐标（以图片左下角的像素点，竖直为y轴，水平为x轴）建立一系列的转换关系，通过曲线的识别，可获得精确的坐标值，并通过相关公式对其加以处理，从而获得。这对文献上发泡聚苯乙烯泡沫 曲线的研究提供了一种便捷而又精确的方法，更能有效地避免对该材料进行重复实验。此外，也可以将实验获得的缓冲材料应力和应变的数据以文本的格式，利用MATLAB程序导入，自动转化成EPS的 、 和 曲线。 2.1 建立EPS缓冲性能数据库 2.1.1 读取所需处理的应力-应变曲线 在文献中选择并拍取所需要的某个密度的发泡聚苯乙烯缓冲材料的 曲线，以.jpg、.bmp等格式放入MATLAB的工作文件中。通过uigetfie等函数，使MATLAB以对话框的形式，读取工作文件夹中的图片文件，然后导入工作区。 2.1.2 寻找图像坐标刻度与图片的像素关系 用text(x,y，'string')函数提示用户选取在x轴上点取相邻两个坐标刻度点。执行程序ginput（2），可为使用者提供十字光标，能更加精确地在图片坐标x轴上点取相邻两个坐标刻度点，MATLAB会自动记录这两个刻度点的像素坐标，分别为 和 。这样在照片坐标x轴上一单位刻度对应的数值与像素坐标上 个像素值建立起联系，即图片x坐标轴一单位刻度包含 个照片像素。同理，在函数的提示下，再一次利用ginput（2）函数在图片坐标y轴上点取相邻两个坐标刻度点,分别为 和 ，图片x坐标轴一单位刻度包含 个像素。 假设聚苯乙烯泡沫应力-应变曲线图片上x轴和y轴一单位刻度数值分别为A和B，所以可以通过程序算出：在x轴方向上，一个像素在图片坐标中的值 ；在y轴方向上，一个像素在图片坐标中的值 。 2.1.3 将图像曲线上选取的点转化为像素坐标系上的点坐标值 假设选取点个数为C，通过编辑程序，在图片界面上显示出一段提示语“请在应力-应变曲线上选取C个点 ”。同时执行程序ginput（C），使用者可以用出现的十字光标，在应力-应变曲线上精确的选取C个点，以矩阵的形式，返回这些点的像素坐标值。一般情况下，应力-应变曲线从坐标原点开始，所以选取点从坐标原点开始，在水平方向上间距尽量要相等，且从左向右拾取，不可颠倒次序。在应力-应变曲线上选取点的个数C，可根据曲线而决定，只要能将应力-应变曲线的形状和具体位置准确无误地选取出来。 2.1.4 建立新图像坐标，将选取的点的像素坐标值转化为真实的试验数值 假设在曲线上选取的点分别是 、 、 、...、 。便能建立起以 为原点，向右每增加一个像素，x轴数值便增加 QUOTE ，向上每增加一个像素，y轴数值便增加 。选取点在新图像坐标轴中的横向数值为该点的水平像素数值与 之差，再乘以 ；纵向数值为该点的竖直像素数值与 之差，再乘以 。例如：在EPS应力-应变曲线上选取的点 ，它在新图像坐标的x轴方向上的值 ；y轴方向上的值 。如此，将选取的点的像素坐标值转化为真实的试验数值 、 、 、...、 。 将这些用上述关系转化出来的真实试验数值，形成一个矩阵 ，赋值于xn。最后，使用MATLAB函数dlmwrite('eps test.txt',xn,'newline','pc')，便可在MATLAB的工作文件夹中出现一个被命名为“eps test”的文本文件，即该材料应力-应变的数值。 将不同密度的EPS材料应力-应变曲线导入MATLAB/GUI界面，通过上述方式转化成各自的应力-应变数据，以文本文件的格式储存在一起，从而建立起一个不同密度EPS材料缓冲性能的数据库。 2.2 利用应力-应变数据计算并绘制缓冲曲线 2.2.1 读取MATLAB的工作文件夹中的文本文件 将经过图片像素关系转化或者实验获得的数据，以文本文件的格式储存于MATLAB的工作文件夹中。通过uigetfie等函数，使MATLAB以对话框的形式，读取工作文件夹中的数据文本文件，然后导入工作区。 2.2.2 建立新图形窗口，绘制EPS的 曲线图 若读取的文本文件里的数据是EPS材料发生应变 所产生对应的应力 。文件的数据中应变 和应力 是成对出现，而形成的一个矩阵。利用figure函数命令，来创建一个新图形窗口。在着新的图形窗口中，通过plot( , )函数，便可将这个矩阵中的坐标点绘制成二维图，点和点之间呈线性连接。如此，便将图片像素关系转化或者实验获得的坐标数据绘制成了EPS的应力-应变曲线图。 2.2.3 绘制成 曲线图 使用uigetfie等函数，将选取点转化后的试验数值 、 、 、...、 导入工作区。根据缓冲系数公式（2.1）进行计算，绘制EPS的 曲线。 （2.1） 式中： ——EPS的缓冲系数； ——不同应变所对应最大应力，N/m2； ——不同应力所对应的形变能，N m。 当应变为 时，该材料的 等于其应力-应变曲线函数对应变在0到 上求定积分，即应力在0到 内，应力-应变曲线与横坐标轴所围成的面积。在MATLAB程序中，使用循环语句进行面积叠加的方法，近似计算该材料的 的值。然后代入缓冲系数的公式，计算其值，如式（2.2）、式（2.3）、式（2.4）所示。 当应力为 QUOTE =0时， QUOTE ； （2.2） 当应力为 时， QUOTE ； （2.3） 当应力为 QUOTE 。 （2.4） QUOTE ， 时， 将 QUOTE 、、 、...、 QUOTE 形成一个矩阵, 通过运用函数plot ，可在新的图形窗口中，绘制成EPS的 曲线图。 2.2.4 绘制 QUOTE 曲线图 通过缓冲系数公式（2.1） QUOTE ，计算缓冲系数-最大应力值，形成的矩阵。设定跌落高度H = 0.6m，在EPS缓冲衬垫的厚度为2cm、3cm、4cm、5cm、6cm、7cm和8cm这7种情况下，将缓冲系数-最大应力 代入静应力公式（2.2），算出静应力 QUOTE ，得到矩阵。 QUOTE （2.5） 式中： ——EPS产生的静应力，Pa； ——EPS衬垫的厚度，cm； ——应变所产生的最大应力，Pa； ——EPS的缓冲系数； ——实验跌落高度，m。 得到矩阵 后，将矩阵代入最大加速度公式（2.3），计算最大加速度 ,得到不同应力 所对应的最大加速度 ，即矩阵 。 式中： ——最大无量纲加速度。 通过绘图函数，绘制成EPS的 2.3 EPS缓冲数据库界面程序制作 在计算缓冲系数、最大加速度等程序的条件下，明确求解缓冲系数、最大加速度等步骤，如图2.1和图2.2所示。将这些求解的步骤编辑成程序，用MATLAB/GUI制作成程序界面，如图2.3和图2.4所示。将编辑完成的程序和界面生成相应的软件包，用于绘制EPS的 、 和 曲线[5]。只要拍取所需要的EPS材料应力-应变曲线照片，将这些照片导入到制作好的应用程序之中，同时设定好跌落高度H，便可求解并绘制EPS材料的 、 和 等曲线[15]。 图2.1 建立EPS缓冲性能数据库的流程图 图2.2 计算并绘制缓冲曲线的流程图 图2.3 将曲线图片转化点坐标数据的设计界面 图2.4 计算并绘制缓冲曲线的设计界面 2.4示 例 将EPS材料应力-应变曲线图片导入工作区，利用“epschuli1”程序对图像进行处理后，所得到的曲线坐标值会生成一个文本文件，如图2.6。并将这个文本文件保存在相应的工作区。然后再利用“epschuli2”程序，将该文本文件中的数据，计算并绘制成EPS的 、 和 曲线，如图2.7所示。 图2.5 导入程序的EPS应力-应变曲线图片 图2.6 图片曲线转化成坐标的文本文件 图2.7 程序绘制EPS的 、 和 曲线 基于拍摄所得的EPS缓冲材料应力-应变曲线图片，利用MATLAB/GUI开发的程序，建立不同密度的EPS缓冲性能的数据库，求解并绘制材料的应力-应变曲线等缓冲曲线，并可脱离MATLAB程序。该软件方便、易操作，导入应力-应变曲线图片，可计算获得EPS的 、 和 数据。从而便于包装缓冲结构设计，促进了实验数据共享，更好地用于实践。 3 利用MATLAB/GUI求解EPS缓冲作用下单自由度缓冲性能 3.1 求解EPS缓冲作用下单自由度缓冲性能程序设计 EPS缓冲衬垫，主要是为了防止包装件不触底的情况下，增加包装件受冲击的时间，来减少冲击载荷；使振动系统的频率远离包装件的固有频率，防止发生共振，避免包装件在运输中所受的冲击和振动作用下造成的损坏[16]。 主要通过包装件的脆值Gc、包装件的质量m、跌落高度H，以及EPS的应力-应变曲线来设计缓冲衬垫的缓冲面积A和衬垫厚度h。 3.1.1 设计输入单自由度的缓冲环境的程序 在学术文献中收集不同密度的应力-应变曲线的照片。用第二章所述的方式，基于MATLAB/GUI，将照片中的曲线转化成的应力-应变数据，并以文本的形式储存在MATLAB的工作文件夹中，建立起数据库，如：25kg/m3 、80 、60 QUOTE 、50等密度。设置一个popupmenu下拉菜单，打开属性框，在string属性中填写不同的密度，每一个不同密度值用回车分开，反映在程序中。当使用者选中某个密度时，该密度的发泡聚苯乙烯缓冲材料对应的文本文件中的应力-应变数据可以直接导入程序中。求该密度的缓冲材料发生应变时，产生的应变能，等于其应力-应变曲线函数在0到对应应变 内， QUOTE 上求定积分，即应力在0到曲线与横坐标轴围成的面积[17]。在程序中，使用循环语句进行面积叠加的方法，近似计算该材料发生所应变对应产生的应变能。然后使用缓冲系数的公式，计算不同应变对应的缓冲系数值。不一样的应变所对应的 和缓冲系数，可以以此组成 数组。利用min(i)函数在这组数列里选择出缓冲系数的最小值，所对应的为 值。在制作包装衬垫时，应该尽量选在C最小值的附近，因为这样比较节省包装衬垫的材料。通过这种方法可以得到不同密度的C以及对应的 值[18]。 再将包装件的脆值Gc、包装件的质量m和跌落高度H，通过str2num函数输入的字符串转化成数值导入到运行的程序之中。 3.1.2 计算缓冲衬垫的最小缓冲面积 对于密度较大的包装件，需要将设计较大的缓冲面积，而对于密度较小的产品，则需要减小受压接触面积，计算最适合缓冲面积。 将设定的包装件的脆值Gc、包装件的质量m、模拟环境的高度 H和程序根据所选密度材料的 数据计算获得的 ，代入公式（3.1）中，便可获得所选密度发泡聚苯乙烯缓冲材料最适合的大小。 QUOTE （3.1） 式中， ——包装衬垫的最小缓冲面积，cm2； ——包装件的脆值，g； ——包装件的质量，kg； ——重力加速度，m/s2； ——应变所产生的最大应力，kPa。 3.1.3 计算缓冲衬垫的厚度 将设定的包装件的脆值 、跌落高度H和程序根据所选密度材料的 数据获得的缓冲系数C，代入公式（3.2）中，可获取所选密度材料的厚度h。 QUOTE （3.2） 式中， ——包装衬垫的最佳厚度，cm； ——包装件的脆值，g； ——缓冲衬垫的最小缓冲系数； ——包装件的跌落高度，m。 由此可计算得到缓冲衬垫的大小和最佳厚度，需要根据公式（3.3）进行校核。如果满足此式，则可以采用输出的缓冲衬垫的大小和最佳厚度进行制作；否则，需要重新计算缓冲衬垫。 QUOTE （3.3） 式中， ——包装衬垫的最小缓冲面积，cm2； ——包装衬垫的最佳厚度，cm。 3.2 EPS缓冲作用下单自由度缓冲性能界面程序制作 在计算缓冲衬垫的面积、衬垫厚度等设计程序的条件下，明确求解缓冲衬垫的面积、衬垫厚度等计算步骤，如图3.1所示。将这些求解的步骤编辑成程序，用MATLAB/GUI制作成程序界面，如图3.2所示。将编辑完成的程序和界面生成相应的软件包，用于设计缓冲衬垫的最优解。只要选择EPS材料的密度，设定包装件地脆值Gc（g）、包装件的质量m（kg）、跌落高度H（m），便可求解缓冲衬垫的面积和厚度。 图3.1 单自由度EPS缓冲衬垫设计流程图 图3.2 单自由度EPS缓冲衬垫的设计界面 3.4 示例 打开“untitied”应用文件，得到一个可以输入的应用程序的界面，如图3.3所示。首先选择EPS衬垫材料的密度：有25、50、60、80kg/m3cm2 QUOTE 等密度可供选择，若无所需的发泡聚苯乙烯缓冲材料的密度。可在MATLAB工作文件中加入所需材料的应力-应变实验数据。然后在应用程序的界面上输入跌落高度H（m）、包装件脆值Gc（g）、包装件的质量m（kg）来设定跌落环境，点击“计算缓冲系数最优解”可以输出该环境下发泡聚苯乙烯缓冲衬垫的设计最优参数，即缓冲面积A（）和衬垫厚度h（cm）。例如，在界面中选择包装产品所需的发泡聚苯乙烯缓冲材料的密度为50kg/m3，并输入“包装件地脆值Gc（g）=12、包装件的质量m（kg）=1、跌落高度H（m）=0.5”，然后按下“计算缓冲系数最优解”按钮，可输出该环境下EPS衬垫的设计最优参数，缓冲面积A（cm2）=13.44和衬垫厚度h（cm）=8.53608，如图3.3所示。 图3.3 EPS缓冲作用下单自由度缓冲衬垫设计图 基于MATLAB/GUI开发出的应用程序用于求解不同环境下EPS衬垫的设计最优参数，并可脱离MATLAB程序运行，为单自由度缓冲作用下的缓冲包装结构设计提供理论数据。 4 利用MATLAB/GUI求解二自由度EPS缓冲性能 在目前的学术研究中，已经有通过有限元的方法获取某种衬垫的应力-应变曲线，来计算并设计不同跌落环境中这种材料所需要的衬垫形状。在实验模拟和理论推导的基础上，通过建立相关衬垫的力学模型，来进行衬垫设计。高德[14]等通过缓冲测试建立了一些缓冲衬垫的多个自由度的本构模型，计算对整个产品的保护性能，但未计算易损件的加速度相应。包装件的质量出现问题，一般都集中在某几个易损部位，包装质量的好坏基本上取决于衬垫对易损件的保护。本次毕业设计，是基于MATLAB/GUI程序，通过输入相关参数。来设置跌落和缓冲环境，根据三次样条插值的方法，来计算二自由度的发泡聚苯乙烯缓冲性能，并输出包装易损件的加速度相应图和包装件冲击脉冲图。 4.1 包装产品易损件的EPS本构模型 本次毕业设计研究的是具有线弹性的包装易损件在发泡聚苯乙烯泡沫塑料的缓冲保护作用下，和包装件共同发生跌落冲击作用的整个力学过程。整个过程的简图，如图4.1所示。 图4.1 EPS缓冲衬垫的跌落模型 在考虑易损件质量远小于包装件的质量的情况下，则从泡沫压缩开始计算，可以列出这个缓冲跌落模型的初始动力学方程，由式（4.1）、（4.2）、（4.3）、（4.4）可见： QUOTE （4.1） QUOTE （4.2） QUOTE （4.3） QUOTE （4.4） 式中： ——包装件的质量，kg； ——包装易损件的质量，g； ——包装易损件的速度，m/s2； ——包装易损件的速度，m/s； ——包装件的相对位移，m； ——包装易损件的弹性系数； ——环境设定的跌落高度，m； ——缓冲衬垫的弹性应力，N/m2 ——有效的缓冲面积，m2。 经过查阅文献、实验研究和理论论证，得到发泡聚苯乙烯属于非线性弹簧，它的缓冲本构模型[19]，由式（4.5）可见： （4.5） 式中： ——包装件对缓冲衬垫产生的应力，N/m2； = 0.02MPa/s； = 0.12MPa/s； = 0.008MPa/s； =13. 2MPa/s； QUOTE MPa/s。 = 0.0037MPa/s； 由此，可得到发泡聚苯乙烯的弹性力F，由式（4.6）可见 （4.6） 将发泡聚苯乙烯的弹性力F代入公式（4.2），忽略易损件的影响，可将该式转化为式（4.7）： （4.7） 根据公式（4.1）和公式（4.7）可得到该缓冲跌落模型的最终动力学方程。 4.2 缓冲过程分析以及运行程序设计 本次设计主要分析的是两大阶段：第一阶段为包装件的缓冲衬垫发生压缩开始到压缩到衬垫厚度的最小值；第二阶段为包装衬垫开始回弹到恢复到原型。利用str2num函数将输入包装件的质量M、跌落高度H、缓冲面积A、缓冲衬垫的厚度h和包装易损件的质量m的字符串转化成数值导入到运行的程序之中。通过对参数的输入，来设计跌落和缓冲的环境。 4.2.1 EPS衬垫压缩阶段缓冲分析 用format函数将程序输出格式设置成为双精度。设置求解器的参数，用 odeset函数来对变量options进行赋值，用来检测发泡剂苯乙烯泡沫压缩的过程。用ode45函数作为求解器来求解发泡聚苯乙烯的动力学方程，求解前定义速度[0；- ]为初始值列向量，在时间为0至0.1秒之间，通过对程序输入包装件质量、产品易损件质量、有效缓冲面积、衬垫厚度等，来计算包装件在不同时刻所发生的位移和速度。 根据公式（4.2）转化，计算出包装件的无量纲加速度，如公式（4.8）所示。 （4.8） 式中： 计算包装件将发泡聚苯乙烯衬垫压缩至最小值时，即压缩的距离 最大时 QUOTE (4.9) QUOTE (4.10) odeset函数对发泡剂苯乙烯泡沫压缩过程的检测设置，需要设置聚苯乙烯泡沫压缩至厚度最小值的临界点。定义一个events函数，作为事件的检测函数，进行过零点检测。编辑程序direction = -1，设置检测过零点的方向为由正到负。当由正到负的方向检测到零点这一指定条件时，通过设置isterminal = 1，来终止ode45函数的运行。 用axes和plot函数将衬垫前期压缩时间和对应包装件的加速度倍数进行绘图，即泡沫压缩阶段的冲击脉冲图。 4.2.2 EPS衬垫回弹阶段缓冲分析 设置回弹缓冲求解器的参数，用来识别缓冲衬垫回弹至恢复的过程。用ode45函数作为求解器来求解发泡聚苯乙烯回弹方程，求解前定义位移[ ；0]为初始值列向量，在时间为0至0.1秒之间，通过输入包装件质量、产品易损件质量、有效缓冲面积、衬垫厚度以及衬垫压缩至最小值的应力和应变，计算包装件在回弹过程中发生的位移和所对应的速度。 定义一个发泡聚苯乙烯泡沫回弹的函数，命名为rebound函数。设置发泡聚苯乙烯缓冲衬垫的弹性模量为 Pa。通过公式（4.11）来计算在回弹的过程中包装件的加速度。 (4.11) 根据公式（4.11），计算出包装件的加速度倍数，如公式（4.12）所示。 (4.12) odeset函数对发泡剂苯乙烯泡沫压缩过程的检测设置，需要设置聚苯乙烯泡沫压缩至厚度最小值的临界点。定义一个events2函数，设置由正到负的方向过零点的识别。当满足识别的条件时，终止ode45函数的运行。 用绘图函数绘制发泡聚苯乙烯泡沫回弹阶段的冲击脉冲图。 4.2.3 包装产品易损件加速度响应分析 设置包装易损件的固有频率 = 80Hz QUOTE 以及弹性阻尼 = 0.1，由于包装产品易损件的加速度响应是离散的，需要定义一个离散型的微分函数。根据三次样条插值的方法，以三次方的形式在缓冲衬垫开始压缩至回弹离开这一过程中，对时间和所对应的最大加速度进行内插，来计算包装产品易损件的加速度响应。 计算包装产品易损件的加速度响应方程，如式（4.13）所示。 QUOTE （4.13） 在式（4.13）中同时减去 ，可将其转化为式（4.14）。 QUOTE （4.14） 令 ，则将式（4.14）转化为式（4.15）。 QUOTE （4.15） 式中： ——包装易损件的相对位移，m； ——包装易损件的相对速度，m/s； ——包装易损件的相对加速度，m/s2。 在式（4.15）两边同除以m，得到式（4.16）。 QUOTE （4.16） 整理公式（4.16），可得包装产品易损件的相对加速度，如式（4.17）所示。 QUOTE （4.17） 式中： ——易损件衬垫的固有频率，Hz； ——易损件衬垫的阻尼系数。 则包装产品易损件的加速度，可以通过式（4.18）计算得到。 QUOTE （4.18） 在程序最后设置“清零”按钮，将程序显示冲击脉冲图和易损件加速度响应图清空，便于再一次使用。 4.3 EPS缓冲作用下易损件的加速度响应曲线界面程序制作 在计算包装件的冲击脉冲和产品易损件的加速度响应等设计程序的基础上，使用MATLAB/GUI，通过句柄、设置提示文本、参数输入框、绘图等函数,将求解包装件的冲击和易损件加速度等步骤做成软件包，如图4.2所示。该软件包用于绘制包装件冲击脉冲图和易损件加速度响应图。只要在程序的界面上输入“包装产品质量M、有效缓冲面积A、产品易损件质量m、缓冲衬垫的厚度h和产品跌落高度H”来定义冲击和缓冲环境，就可以求解包装件的冲击脉冲和产品易损件加速度响应[20]。 图4.2二自由度EPS缓冲作用下易损件加速度响应计算界面 4.4 示例 打开“yisunjian”应用文件，得到一个可以输入的应用程序的界面，如图4.3所示。使用者可以在这个界面上输入“包装件质量M、缓冲面积A、易损件质量m、缓冲衬垫厚度h、跌落高度H ”来定义振动环境得到“冲击脉冲图”和“易损件加速度响应图”。例如，在界面中输入“包装件质量M = 6 kg、缓冲面积A = 0.012m2、易损件质量m = 0.01kg、缓冲衬垫厚度h = 0.02m、跌落高度H = 0.5m”，然后按下“确定计算”按钮，可在界面的图像显示区域看到该统计环境下的“冲击脉冲图”和“易损件加速度响应图”。 图4.3二自由度EPS缓冲作用下易损件加速度响应计算结果 基于MATLAB/GUI设计开发出的应用程序，用于求解并绘制二自由度EPS缓冲作用下包装件的冲击和易损件加速度，并可脱离MATLAB程序。生成的应用软件简单实用，操作方便，为设计二自由度EPS缓冲衬垫的设计提供了理论基础，更好地用来保护包装产品及其易损件。 5. 总结与展望 当今世界，不分国内还是国外，可持续发展治理环境的呼声亦是愈来愈高。中国已经把环境保护作为一项基本的国策。在欧洲其他发达的城市，已经实行将废弃的缓冲泡沫做到减少用量、循环使用。实践证明：EPS有着缓冲性能好、隔热性佳等优良性能，所以它被大量的推广和使用，但是发泡聚苯乙烯缓冲衬垫易发生“白色污染”等环保问题需要解决。 本次毕业设计选择了对发泡聚苯乙烯缓冲材料的缓冲吸能进行研究。基MATLAB软件，从实际使用者的角度出发，编辑使用程序，设计了四个简单实用的应用程序来研究发泡聚苯乙烯缓冲材料的缓冲性能，并相应生成了操作方便的应用软件。 一、“epschuli1” 程序界面：将不同密度的EPS材料应力-应变曲线的照片导入工作区，通过点坐标的选取将曲线图片转化成各自的应力-应变数据，以文本文件的格式储存在工作文件夹之中，以此建立不同密度的发泡聚苯乙烯衬垫材料的缓冲性能数据库。 二、“epschuli2”程序界面：将转化后的不同密度的发泡聚苯乙烯衬垫材料应力-应变数据，通过界面逻辑函数等程序的计算，获得缓冲衬垫材料的 和 数据，利用绘图软件绘制成EPS的 、 和 曲线。 三、“untitied”程序界面：首先选择EPS衬垫材料的密度，然后在应用程序的界面上输入跌落高度H、包装件地脆值Gc、包装件的质量m来设定跌落环境，点击“计算缓冲系数最优解”可以输出该环境下EPS衬垫的设计最优参数，即缓冲面积A和衬垫厚度h。 四、“yisunjian”程序界面：首先输入“包装件质量M、缓冲面积A、易损件质量m、缓冲衬垫厚度h、跌落高度H”来定义振动环境, 才有三次样条插值的方法，得到“包装件冲击脉冲图”和“易损件加速度响应图”。 有理由相信，在相关科研人员的不懈努力下，通过理论计算和实验验证等多种方法，对包装衬垫的设计进行不断地优化，从而使发泡聚苯乙烯泡沫塑料的使用量逐步减少，使环境的污染得到缓解。 参考文献 [1] 刘晔，祁书艳，史锃莹，任家仪. 竹纤维发泡缓冲包装材料的研究开发[J]. 中国包装，2006，03：26～28. [2] 何辉，王守伏，张继尹. 聚苯乙烯泡沫塑料的回收与利用[J]. 高分子材料科学与工程，2009，v. 2512：162-165. [3] 熊志远. 聚苯乙烯泡沫（EPS）力学行为的实验研究[D]. 湘潭：湘潭大学，2007. [4] 张卫兵. 聚苯乙烯泡沫(EPS)的特性及其在道路工程中的应用[J]. 公路，2004，05：146～149. [5] 佟富强，刘宝璋，杜金来. 形变速率对聚苯乙烯泡沫材料缓冲性能的影响[J]. 塑料包装，1996，04：32～34. [6] 毛快，张俊彦. 聚苯乙烯泡沫(EPS)压缩蠕变实验研究[J]. 材料导报，2007，S3：468～470. [7] 朱向荣，EPS动力特性试验研究[J]. 岩土工程学报. 2005，27（11）：1253～1256. [8] 卢富德，刘雄建，高德. 多胞缓冲材料本构模型与应用进展[J]. 浙江大学学报(工学版)，2014，07：1336～1344. [9] 高德，刘壮，孙智慧. 秸秆纤维EPS缓冲包装材料性能的研究[J]. 农业机械学报，2006，08：201～204. [10] 都学飞，欧阳效卓，张汪年. EPS缓冲材料的静态压缩性能的试验研究[J]. 包装工程，2012，03：39～41+45. [11] Horvath J S. Expanded polystyrene( EPS) geofoam:an introduction to materialbehavior[J]. Geotextiles and Geomembrances，1994(13)：26～280. [12] Dnskov M. Materials research on EPS20 and EPS15 under representative conditions inpavement Structures[J]. Geotextiles and Geomembrances，1997，(15)：17～18 [13] Qunli Liu，Ghatu Subhash. A phenomenological constitutive model for foams under large deformations[J]. Polym Eng Sci，2004，443：. [14] 薛飞. 跌落冲击测试系统设计与EPE缓冲隔振性能研究[D]. 太原：太原科技大学，2014. [15] 郭勇，卢富德，高德，王振林. 基于MATLAB/GUI缓冲包装动力学优化设计[J]. 振动与冲击，2014，02：81～83. [16] 叶柏彰,张涛. 家电产品要打造“绿色包装”纸蜂窝包装材料“走红家电”[J]. 包装世界，2004，02：67～70. [17] 沈训乐. 发泡缓冲材料本构模型研究及其CAD应用[D]. 咸阳：陕西科技大学,2012. [18] 郭彦峰，赵杰，许文才. 组合衬垫缓冲系数的分析与应用[J]. 中国包装工业，2000，06：39～41. [19] Fude Lu, Weiming Tao, De Gao. Virtual Mass Method for Solution of Dynamic Response of Composite Cushion Packaging System[J]. Packag. Technol. Sci., 2013, 26：. [20] 童忠钫，胡强，李培玉. 采用泡沫衬垫的产品缓冲包装优化设计的研究[J]. 浙江大学学报(自然科学版)，1994，01：41～48. 附 录 附录1 建立EPS缓冲性能数据库程序 function varargout = epschuli1(varargin) gui_Singleton = 1; gui_State = struct('gui_Name', mfilename, ... 'gui_Singleton', gui_Singleton, ... 'gui_OpeningFcn', @epschuli1_OpeningFcn, ... 'gui_OutputFcn', @epschuli1_OutputFcn, ... 'gui_LayoutFcn', [] , ... 'gui_Callback', []); if nargin && ischar(varargin{1}) gui_State.gui_Callback = str2func(varargin{1}); end if nargout [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); else gui_mainfcn(gui_State, varargin{:}); End function epschuli1_OpeningFcn(hObject, eventdata, handles, varargin) handles.output = hObject; guidata(hObject, handles); function varargout = epschuli1_OutputFcn(hObject, eventdata, handles) varargout{1} = handles.output; function pushbutton1_Callback(hObject, eventdata, handles) [filename, pathname] = ... uigetfile({'*.jpg';'*.bmp';'*.tif';'*.*'},'File Selector'); a=importdata(filename); figure image(a) text(10,250,'1.请在X轴上量取2点','color','b','Fontsize',10) xb1=ginput(2); xb=0.2/(xb1(2,1)-xb1(1,1)); text(10,200,'2.请在y轴上量取2点','color','b','Fontsize',10) yb1=ginput(2); yb=0.05/(yb1(2,2)-yb1(1,2)); text(10,150,'3.请在应力应变曲线上取30点','color','b','Fontsize',10) x=ginput(30); xx=(x(:,1)-x(1,1))*xb; yy=(x(:,2)-x(1,2))*yb; xn=[xx yy]; dlmwrite('eps test.txt',xn,'newline','pc') 附录2 利用应力-应变数据计算并绘制缓冲曲线程序 function varargout = epschuli2(varargin) gui_Singleton = 1; gui_State = struct('gui_Name', mfilename, ... 'gui_Singleton', gui_Singleton, ... 'gui_OpeningFcn', @epschuli2_OpeningFcn, ... 'gui_OutputFcn', @epschuli2_OutputFcn, ... 'gui_LayoutFcn', [] , ... 'gui_Callback', []); if nargin && ischar(varargin{1}) gui_State.gui_Callback = str2func(varargin{1}); end if nargout [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); else gui_mainfcn(gui_State, varargin{:}); End function epschuli2_OpeningFcn(hObject, eventdata, handles, varargin) handles.output = hObject; guidata(hObject, handles); function varargout = epschuli2_OutputFcn(hObject, eventdata, handles) varargout{1} = handles.output; function pushbutton1_Callback(hObject, eventdata, handles) [filename, pathname] = ... uigetfile({'*.txt';'*.mdl';'*.mat';'*.*'},'File Selector'); a=importdata(filename); xx=a(:,1);yy=a(:,2); figure plot(xx,yy) xlabel('应变') ylabel('应力/Mpa') function pushbutton2_Callback(hObject, eventdata, handles) [filename, pathname] = ... uigetfile({'*.txt';'*.mdl';'*.mat';'*.*'},'File Selector'); a=importdata(filename); xx=a(:,1);yy=a(:,2); s(1)=0; for i=1:length(xx)-1 s(i+1)=s(i)+(yy(i+1)+yy(i))*(xx(i+1)-xx(i))/2; c(i)=yy(i)/s(i+1); end aa=length(xx)/5; figure plot(yy(aa:length(xx)-1),c(aa:length(xx)-1)) xlabel('最大应力/Mpa') ylabel('缓冲系数/C') function pushbutton3_Callback(hObject, eventdata, handles) [filename, pathname] = ... uigetfile({'*.txt';'*.mdl';'*.mat';'*.*'},'File Selector'); a=importdata(filename); xx=a(:,1);yy=a(:,2); s(1)=0; for i=1:length(xx)-1 s(i+1)=s(i)+(yy(i+1)+yy(i))*(xx(i+1)-xx(i))/2; c(i)=yy(i)/s(i+1); end aa=length(xx)/5; yyy=yy(aa:length(xx)-1); ccc=c(aa:length(xx)-1); figure for j=1:7 h(j)=0.01*j+0.01;H=0.6;st=(yyy*h(j))./ccc'/H; Gm=yyy./st;Gm=yyy./st;semilogx(st(aa:length(st))*1000,Gm(aa:length(st))); hold on end xlabel('静应力/Kpa') ylabel('最大加速度/g') 附录3 求解单自由度EPS缓冲衬垫设计程序 function varargout = untitled(varargin) gui_Singleton = 1; gui_State = struct('gui_Name', mfilename, ... 'gui_Singleton', gui_Singleton, ... 'gui_OpeningFcn', @untitled_OpeningFcn, ... 'gui_OutputFcn', @untitled_OutputFcn, ... 'gui_LayoutFcn', [] , ... 'gui_Callback', []); if nargin && ischar(varargin{1}) gui_State.gui_Callback = str2func(varargin{1}); end if nargout [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); else gui_mainfcn(gui_State, varargin{:}); end function untitled_OpeningFcn(hObject, eventdata, handles, varargin) handles.output = hObject; guidata(hObject, handles); function varargout = untitled_OutputFcn(hObject, eventdata, handles) varargout{1} = handles.output; function midu_popupmenu_Callback(hObject, eventdata, handles) global a1; val=get(hObject ,'value'); switch (val) case 1 a1=importdata('12Kg eps.txt'); case 2 a1=importdata('12Kg eps.txt'); case 3 a1=importdata('20Kg eps.txt'); case 4 a1=importdata('30Kg eps.txt'); case 5 a1=importdata('40Kg eps.txt'); case 6 a1=importdata('55Kg eps.txt'); case 7 a1=importdata('70Kg eps.txt'); end function midu_popupmenu_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function jisuan_togglebutton_Callback(hObject, eventdata, handles) global a1; gc=str2num(get(handles.gc_edit,'string')); m=str2num(get(handles.m_edit,'string')); h=str2num(get(handles.h_edit,'string')); global a1; xx=a1(:,1);yy=a1(:,2); s(1)=0; for i=1:length(xx)-1 s(i+1)=s(i)+(yy(i+1)+yy(i))*(xx(i+1)-xx(i))/2; cc(i)=yy(i)/s(i+1); end aa=length(xx)/5; yyy=yy(aa:length(xx)-1); ccc=cc(aa:length(xx)-1); [min_of_cc i]=min(ccc); sm=yyy(i); a=(gc*m*9.8/100)/sm; h1=min_of_cc*h*100./gc; set(handles.h1_edit,'string',h1) set(handles.a_edit,'string',a) function gc_edit_Callback(hObject, eventdata, handles) function gc_edit_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function m_edit_Callback(hObject, eventdata, handles) function m_edit_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function h_edit_Callback(hObject, eventdata, handles) function h_edit_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function a_edit_Callback(hObject, eventdata, handles) function a_edit_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function h1_edit_Callback(hObject, eventdata, handles) function h1_edit_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end 附录4 求解二自由度EPS缓冲作用下缓冲性能程序 （1）求解二自由度EPS缓冲性能主程序 function varargout = yisunjian(varargin) gui_Singleton = 1; gui_State = struct('gui_Name', mfilename, ... 'gui_Singleton', gui_Singleton, ... 'gui_OpeningFcn', @yisunjian_OpeningFcn, ... 'gui_OutputFcn', @yisunjian_OutputFcn, ... 'gui_LayoutFcn', [] , ... 'gui_Callback', []); if nargin && ischar(varargin{1}) gui_State.gui_Callback = str2func(varargin{1}); end if nargout [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); else gui_mainfcn(gui_State, varargin{:}); End function yisunjian_OpeningFcn(hObject, eventdata, handles, varargin) guidata(hObject, handles); function varargout = yisunjian_OutputFcn(hObject, eventdata, handles) varargout{1} = handles.output; function m_edit_Callback(hObject, eventdata, handles) function m_edit_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function m1_edit_Callback(hObject, eventdata, handles) function m1_edit_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function H_edit_Callback(hObject, eventdata, handles) function H_edit_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function A_edit_Callback(hObject, eventdata, handles) function A_edit_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function h2_edit_Callback(hObject, eventdata, handles) function h2_edit_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function jisuan_togglebutton_Callback(hObject, eventdata, handles) b1=0.02; b2=0.12; b3=13.2; b4=0.008; b5=0.0037*0; b6=11.3e-6*0; m=str2num(get(handles.m_edit,'string')); H=str2num(get(handles.H_edit,'string')); A=str2num(get(handles.A_edit,'string')); h2=str2num(get(handles.h2_edit,'string')); m1=str2num(get(handles.m1_edit,'string')); E=6e6;%把这些参数反映在界面上 format long options=odeset('events',@events1); [t,y]=ode45(@EPS,[0,0.1],[0; sqrt(2*9.8*H)],options,m,m1,A,h2); aa= 1/m*A*1e6*((b1*((y(:,1))/h2)+b2*tanh(b3*y(:,1)/h2)+b4*y(:,1)./(0.85-y(:,1)/h2)/h2).*(1+b5*y(:,2)/h2)+b6*y(:,2)/h2)/9.8; nnn=length(t); s0=aa(nnn)*9.8*m/A; e0=y(nnn,1)/h2; axes(handles.axes2); plot(t,aa,'r--.');%压缩泡沫的过程 hold on options=odeset('events',@events2); [t1,y1]=ode45(@Rebound,[0,0.1],[e0*h2; 0],options,m,m1,A,h2,e0,s0); aa1=1/m*A*(E*(y1(:,1)/h2-e0)+s0)/9.8; nnn1=length(t1); tttt=[t;t1+t(nnn)+0.000001];yyyy=[aa;aa1]*9.8; x=[tttt yyyy]; axes(handles.axes2); plot(x(:,1),x(:,2)/9.8,'r-.');%回弹的过程 w=80;xi=0.1; options=odeset('events',@events22); [tt,yy]=ode45(@DiscreteDiffFun,[0,0.1],[0;0],options,w,xi,x); axes(handles.axes1); plot(tt,yy(:,1)/9.8) function qingling_togglebutton_Callback(hObject, eventdata, handles) cla reset; axes(handles.axes2); cla reset; （2）求解二自由度EPS缓冲性能EPS子程序 function dy=EPS(t,y,m,m1,A,h2) b1=0.02; b2=0.12; b3=13.2; b4=0.008; b5=0.0037*0; b6=11.3e-6*0; dy=[y(2);-1/m*A*1e6*((b1*((y(1))/h2)+b2*tanh(b3*y(1)/h2)+b4*y(1)/h2/(0.85-y(1)/h2))*(1+b5*y(2)/h2)+b6*y(2)/h2)]; （3）求解二自由度EPS缓冲性能DiscretediffFun子程序 function dy=DiscreteDiffFun(t,y,w,xi,x) aa=x; xx=aa(:,1)-aa(1,1);yy=aa(:,2);nn=length(xx)*10; tt=max(xx);tt1=tt+0.000001; xxx1=linspace(tt1,10*tt,nn)'; xxx=[xx;xxx1];yyy=[yy;zeros(nn,1)]; zz=interp1(xxx,yyy,t,'cubic'); dy=[y(2);-w^2*y(1)-2*xi*w*y(2)+w^2*zz]; （4）求解二自由度EPS缓冲性能Rebound子程序 function dy=Rebound(t,y,m,m1,A,h2,e0,s0) E=6e6; dy=[y(2);-1/m*A*(E*(y(1)/h2-e0)+s0)]; （5）求解二自由度EPS缓冲性能events1子程序 function [value,isterminal,direction]=events1(t,y,m,m1,A,h2) value=y(2); isterminal=1; direction=-1; （6）求解二自由度EPS缓冲性能events22子程序 function [value,isterminal,direction]=events22(t,y,w,xi,x) E=6e6; value=0; isterminal=1; direction=-1; 致 谢 “岁月如斯夫，不舍昼夜”，我在大学这座“象牙塔”里的美好时光，也从憧憬逐渐成为了回忆。 本次毕业论文的结稿，离不了同学们的关心，离不了亲友们的支持，更离不了老师们的指导。在此之际，我首先要感谢卢富德老师，在冗杂的教学工作之余，在论文的前期准备、中期实验设计以及后期论文写作等方面，给予了我悉心的关怀和耐心的指导。卢老师严谨的教学态度、敬业的工作精神和高尚的学术情怀都激励着我在学术的高山上不断地求索、攀登。 感谢学院的各位领导和老师们，特别是包装工程研究所的各位恩师。在我四年的大学生涯中，你们认真严谨的工作态度，奔竞不息的进取精神，诲人不倦的治学作风，缜密活泼的指导思维，让我获益良多，这将对我今后的学习、生活和工作等方面产生巨大而又深远的影响。同时还感谢和我一起学习、共同成长的包装工程专业的同学们。回想四度寒暑我们共同度过的一点一滴，回想大家在工作、学习和生活上对我的无私帮助。 感谢我的父亲和母亲，在四年大学期间，父母亲对我的信任和支持是我努力学习奋勇拼搏的无形动力。 大学时光一刹那，刹那之间成永恒。感谢你们：我的老师、我的同学、我的理工，因为有你，我的大学生活多彩缤纷；因为有你，我的大学生涯充满意义；因为有你，我的大学时光与众不同。最后真挚向你们道一声 “谢谢”，祝愿同学们将今日之所学，开拓明日成功之道路。 俞坚道 2016年4月30日 毕业设计（论文）原创性声明和使用授权说明 原创性声明 本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得 及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。 作 者 签 名：　　 　 　　日　 期：　　 　 　　 ​​​​​​​​​​​​ 指导教师签名：　　 　 　　 日　　期：　　 　 　　 使用授权说明 本人完全了解 大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。 作者签名：　　 　 　　 日　 期：　　 　 　　 ​​​​​​​​​​​​ 学位论文原创性声明 本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。 作者签名： 日期： 年 月 日 学位论文版权使用授权书 本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权　　 　 　大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。 涉密论文按学校规定处理。 作者签名： 日期： 年 月 日 导师签名： 日期： 年 月 日 指导教师评阅书 指导教师评价： 一、撰写（设计）过程 1、学生在论文（设计）过程中的治学态度、工作精神 □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、学生掌握专业知识、技能的扎实程度 □ 优 □ 良 □ 中 □ 及格 □ 不及格 3、学生综合运用所学知识和专业技能分析和解决问题的能力 □ 优 □ 良 □ 中 □ 及格 □ 不及格 4、研究方法的科学性；技术线路的可行性； 设计方案 关于薪酬设计方案通用技术作品设计方案停车场设计方案多媒体教室设计方案农贸市场设计方案 的合理性 □ 优 □ 良 □ 中 □ 及格 □ 不及格 5、完成毕业论文（设计）期间的出勤情况 □ 优 □ 良 □ 中 □ 及格 □ 不及格 二、论文（设计）质量 1、论文（设计）的整体结构是否符合撰写规范？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、是否完成指定的论文（设计）任务（包括装订及附件）？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 三、论文（设计）水平 1、论文（设计）的理论意义或对解决实际问题的指导意义 □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、论文的观念是否有新意？设计是否有创意？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 3、论文（设计说明书）所体现的整体水平 □ 优 □ 良 □ 中 □ 及格 □ 不及格 建议成绩：□ 优 □ 良 □ 中 □ 及格 □ 不及格 （在所选等级前的□内画“√”） 指导教师： （签名） 单位： （盖章） 年 月 日 评阅教师评阅书 评阅教师评价： 一、论文（设计）质量 1、论文（设计）的整体结构是否符合撰写规范？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、是否完成指定的论文（设计）任务（包括装订及附件）？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 二、论文（设计）水平 1、论文（设计）的理论意义或对解决实际问题的指导意义 □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、论文的观念是否有新意？设计是否有创意？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 3、论文（设计说明书）所体现的整体水平 □ 优 □ 良 □ 中 □ 及格 □ 不及格 建议成绩：□ 优 □ 良 □ 中 □ 及格 □ 不及格 （在所选等级前的□内画“√”） 评阅教师： （签名） 单位： （盖章） 年 月 日 教研室（或答辩小组）及教学系意见 教研室（或答辩小组）评价： 一、答辩过程 1、毕业论文（设计）的基本要点和见解的叙述情况 □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、对答辩问题的反应、理解、表达情况 □ 优 □ 良 □ 中 □ 及格 □ 不及格 3、学生答辩过程中的精神状态 □ 优 □ 良 □ 中 □ 及格 □ 不及格 二、论文（设计）质量 1、论文（设计）的整体结构是否符合撰写规范？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、是否完成指定的论文（设计）任务（包括装订及附件）？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 三、论文（设计）水平 1、论文（设计）的理论意义或对解决实际问题的指导意义 □ 优 □ 良 □ 中 □ 及格 □ 不及格 2、论文的观念是否有新意？设计是否有创意？ □ 优 □ 良 □ 中 □ 及格 □ 不及格 3、论文（设计说明书）所体现的整体水平 □ 优 □ 良 □ 中 □ 及格 □ 不及格 评定成绩：□ 优 □ 良 □ 中 □ 及格 □ 不及格 （在所选等级前的□内画“√”） 教研室主任（或答辩小组组长）： （签名） 年 月 日 教学系意见： 系主任： （签名） 年 月 日 学位论文原创性声明 本人郑重声明：所呈交的学位论文，是本人在导师的指导下进行的研究工作所取得的成果。尽我所知，除文中已经特别注明引用的内容和致谢的地方外，本论文不包含任何其他个人或集体已经发表或撰写过的研究成果。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式注明并表示感谢。本人完全意识到本声明的法律结果由本人承担。 学位论文作者（本人签名）： 年 月 日 学位论文出版授权书 本人及导师完全同意《中国博士学位论文全文数据库出版章程》、《中国优秀硕士学位论文全文数据库出版章程》(以下简称“章程”)，愿意将本人的学位论文提交“中国学术期刊（光盘版）电子杂志社”在《中国博士学位论文全文数据库》、《中国优秀硕士学位论文全文数据库》中全文发表和以电子、网络形式公开出版，并同意编入****《中国知识资源总库》，在《中国博硕士学位论文评价数据库》中使用和在互联网上传播，同意按“章程”规定享受相关权益。 论文密级： □公开 □保密（___年__月至__年__月）(保密的学位论文在解密后应遵守此协议) 作者签名：_______ 导师签名：_______ _______年_____月_____日 _______年_____月_____日 独 创 声 明 本人郑重声明：所呈交的毕业设计(论文)，是本人在指导老师的指导下，独立进行研究工作所取得的成果，成果不存在知识产权争议。尽我所知，除文中已经注明引用的内容外，本设计（论文）不含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体均已在文中以明确方式标明。 本声明的法律后果由本人承担。   作者签名: 二〇一〇年九月二十日   毕业设计（论文）使用授权声明 本人完全了解**学院关于收集、保存、使用毕业设计（论文）的规定。 本人愿意按照学校要求提交学位论文的印刷本和电子版，同意学校保存学位论文的印刷本和电子版，或采用影印、数字化或其它复制手段保存设计（论文）；同意学校在不以营利为目的的前提下，建立目录检索与阅览服务系统，公布设计（论文）的部分或全部内容，允许他人依法合理使用。 （保密论文在解密后遵守此规定）   作者签名: 二〇一〇年九月二十日 致 谢 时间飞逝，大学的学习生活很快就要过去，在这四年的学习生活中，收获了很多，而这些成绩的取得是和一直关心帮助我的人分不开的。 首先非常感谢学校开设这个课题，为本人日后从事计算机方面的工作提供了经验，奠定了基础。本次毕业设计大概持续了半年，现在终于到结尾了。本次毕业设计是对我大学四年学习下来最好的检验。经过这次毕业设计，我的能力有了很大的提高，比如操作能力、分析问题的能力、合作精神、严谨的工作作风等方方面面都有很大的进步。这期间凝聚了很多人的心血，在此我表示由衷的感谢。没有他们的帮助，我将无法顺利完成这次设计。 首先，我要特别感谢我的知道***老师对我的悉心指导，在我的论文书写及设计过程中给了我大量的帮助和指导，为我理清了设计思路和操作方法，并对我所做的课题提出了有效的改进方案。***老师渊博的知识、严谨的作风和诲人不倦的态度给我留下了深刻的印象。从他身上，我学到了许多能受益终生的东西。再次对周巍老师表示衷心的感谢。 其次，我要感谢大学四年中所有的任课老师和辅导员在学习期间对我的严格要求，感谢他们对我学习上和生活上的帮助，使我了解了许多专业知识和为人的道理，能够在今后的生活道路上有继续奋斗的力量。 另外，我还要感谢大学四年和我一起走过的同学朋友对我的关心与支持，与他们一起学习、生活，让我在大学期间生活的很充实，给我留下了很多难忘的回忆。 最后，我要感谢我的父母对我的关系和理解，如果没有他们在我的学习生涯中的无私奉献和默默支持，我将无法顺利完成今天的学业。 四年的大学生活就快走入尾声，我们的校园生活就要划上句号，心中是无尽的难舍与眷恋。从这里走出，对我的人生来说，将是踏上一个新的征程，要把所学的知识应用到实际工作中去。 回首四年，取得了些许成绩，生活中有快乐也有艰辛。感谢老师四年来对我孜孜不倦的教诲，对我成长的关心和爱护。 学友情深，情同兄妹。四年的风风雨雨，我们一同走过，充满着关爱，给我留下了值得珍藏的最美好的记忆。 在我的十几年求学历程里，离不开父母的鼓励和支持，是他们辛勤的劳作，无私的付出，为我创造良好的学习条件，我才能顺利完成完成学业，感激他们一直以来对我的抚养与培育。 最后，我要特别感谢我的导师***老师、和研究生助教***老师。是他们在我毕业的最后关头给了我们巨大的帮助与鼓励，给了我很多解决问题的思路，在此表示衷心的感激。老师们认真负责的工作态度，严谨的治学精神和深厚的理论水平都使我收益匪浅。他无论在理论上还是在实践中，都给与我很大的帮助，使我得到不少的提高这对于我以后的工作和学习都有一种巨大的帮助，感谢他耐心的辅导。在论文的撰写过程中老师们给予我很大的帮助，帮助解决了不少的难点，使得论文能够及时完成，这里一并表示真诚的感谢。 致 谢 这次论文的完成，不止是我自己的努力，同时也有老师的指导，同学的帮助，以及那些无私奉献的前辈，正所谓你知道的越多的时候你才发现你知道的越少，通过这次论文，我想我成长了很多，不只是磨练了我的知识厚度，也使我更加确定了我今后的目标：为今后的计算机事业奋斗。在此我要感谢我的指导老师——***老师，感谢您的指导，才让我有了今天这篇论文，您不仅是我的论文导师，也是我人生的导师，谢谢您！我还要感谢我的同学，四年的相处，虽然我未必记得住每分每秒，但是我记得每一个有你们的精彩瞬间，我相信通过大学的历练，我们都已经长大，变成一个有担当，有能力的新时代青年，感谢你们的陪伴，感谢有你们，这篇论文也有你们的功劳，我想毕业不是我们的相处的结束，它是我们更好相处的开头，祝福你们！我也要感谢父母，这是他们给我的，所有的一切；感谢母校，尽管您不以我为荣，但我一直会以我是一名农大人为荣。 通过这次毕业设计，我学习了很多新知识，也对很多以前的东西有了更深的记忆与理解。漫漫求学路，过程很快乐。我要感谢信息与管理科学学院的老师，我从他们那里学到了许多珍贵的知识和做人处事的道理，以及科学严谨的学术态度，令我受益良多。同时还要感谢学院给了我一个可以认真学习，天天向上的学习环境和机会。 即将结束*大学习生活，我感谢****大学提供了一次在**大接受教育的机会，感谢院校老师的无私教导。感谢各位老师审阅我的论文。 毕业设计（论文）原创性声明和使用授权说明 原创性声明 本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得 及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。 作 者 签 名：　　 　 　　日　 期：　　 　 　　 ​​​​​​​​​​​​ 指导教师签名：　　 　 　　 日　　期：　　 　 　　 使用授权说明 本人完全了解 大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。 作者签名：　　 　 　　 日　 期：　　 　 　　 ​​​​​​​​​​​​ 学位论文原创性声明 本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。 作者签名： 日期： 年 月 日 学位论文版权使用授权书 本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权　　 　 　大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。 涉密论文按学校规定处理。 作者签名： 日期： 年 月 日 导师签名： 日期： 年 月 日 独 创 声 明 本人郑重声明：所呈交的毕业设计(论文)，是本人在指导老师的指导下，独立进行研究工作所取得的成果，成果不存在知识产权争议。尽我所知，除文中已经注明引用的内容外，本设计（论文）不含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体均已在文中以明确方式标明。 本声明的法律后果由本人承担。   作者签名: 年 月 日   毕业设计（论文）使用授权声明 本人完全了解**学院关于收集、保存、使用毕业设计（论文）的规定。 本人愿意按照学校要求提交学位论文的印刷本和电子版，同意学校保存学位论文的印刷本和电子版，或采用影印、数字化或其它复制手段保存设计（论文）；同意学校在不以营利为目的的前提下，建立目录检索与阅览服务系统，公布设计（论文）的部分或全部内容，允许他人依法合理使用。 （保密论文在解密后遵守此规定）   作者签名: 年 月 日 基本要求：写毕业论文主要目的是培养学生综合运用所学知识和技能，理论联系实际，独立分析，解决实际问题的能力，使学生得到从事本专业工作和进行相关的基本训练。毕业论文应反映出作者能够准确地掌握所学的专业基础知识，基本学会综合运用所学知识进行科学研究的方法，对所研究的题目有一定的心得体会，论文题目的范围不宜过宽，一般选择本学科某一重要问题的一个侧面。 毕业论文的基本教学要求是： 1、培养学生综合运用、巩固与扩展所学的基础理论和专业知识，培养学生独立分析、解决实际问题能力、培养学生处理数据和信息的能力。2、培养学生正确的理论联系实际的工作作风，严肃认真的科学态度。3、培养学生进行社会调查研究；文献资料收集、阅读和整理、使用；提出论点、综合论证、总结写作等基本技能。 毕业论文是毕业生总结性的独立作业，是学生运用在校学习的基本知识和基础理论，去分析、解决一两个实际问题的实践锻炼过程，也是学生在校学习期间学习成果的综合性总结，是整个教学活动中不可缺少的重要环节。撰写毕业论文对于培养学生初步的科学研究能力，提高其综合运用所学知识分析问题、解决问题能力有着重要意义。 毕业论文在进行编写的过程中，需要经过开题报告、论文编写、论文上交评定、论文答辩以及论文评分五个过程，其中开题报告是论文进行的最重要的一个过程，也是论文能否进行的一个重要指标。 撰写意义：1.撰写毕业论文是检验学生在校学习成果的重要措施，也是提高教学质量的重要环节。大学生在毕业前都必须完成毕业论文的撰写任务。申请学位必须提交相应的学位论文，经答辩通过后，方可取得学位。可以这么说，毕业论文是结束大学学习生活走向社会的一个中介和桥梁。毕业论文是大学生才华的第一次显露，是向祖国和人民所交的一份有份量的答卷，是投身社会主义现代化建设事业的报到书。一篇毕业论文虽然不能全面地反映出一个人的才华，也不一定能对社会直接带来巨大的效益，对专业产生开拓性的影响。但是，实践证明，撰写毕业论文是提高教学质量的重要环节，是保证出好人才的重要措施。 2.通过撰写毕业论文，提高写作水平是干部队伍“四化”建设的需要。党中央要求，为了适应现代化建设的需要，领导班子成员应当逐步实现“革命化、年轻化、知识化、专业化”。这个“四化”的要求，也包含了对干部写作能力和写作水平的要求。 3.提高大学生的写作水平是社会主义物质文明和精神文明建设的需要。在新的历史时期，无论是提高全族的科学文化水平，掌握现代科技知识和科学管理方法，还是培养社会主义新人，都要求我们的干部具有较高的写作能力。在经济建设中，作为领导人员和机关的办事人员，要写指示、通知、总结、 调查报告 行政管理关于调查报告关于XX公司的财务调查报告关于学校食堂的调查报告关于大米市场调查报告关于水资源调查报告 等应用文；要写说明书、广告、解说词等说明文；还要写科学论文、经济评论等议论文。在当今信息社会中，信息对于加快经济发展速度，取得良好的经济效益发挥着愈来愈大的作用。写作是以语言文字为信号，是传达信息的方式。信息的来源、信息的收集、信息的储存、整理、传播等等都离不开写作。 论文种类：毕业论文是学术论文的一种形式，为了进一步探讨和掌握毕业论文的写作规律和特点，需要对毕业论文进行分类。由于毕业论文本身的内容和性质不同，研究领域、对象、方法、表现方式不同，因此，毕业论文就有不同的分类方法。 按内容性质和研究方法的不同可以把毕业论文分为理论性论文、实验性论文、描述性论文和设计性论文。后三种论文主要是理工科大学生可以选择的论文形式，这里不作介绍。文科大学生一般写的是理论性论文。理论性论文具体又可分成两种：一种是以纯粹的抽象理论为研究对象，研究方法是严密的理论推导和数学运算，有的也涉及实验与观测，用以验证论点的正确性。另一种是以对客观事物和现象的调查、考察所得观测资料以及有关文献资料数据为研究对象，研究方法是对有关资料进行分析、综合、概括、抽象，通过归纳、演绎、类比，提出某种新的理论和新的见解。 按议论的性质不同可以把毕业论文分为立论文和驳论文。立论性的毕业论文是指从正面阐述论证自己的观点和主张。一篇论文侧重于以立论为主，就属于立论性论文。立论文要求论点鲜明，论据充分，论证严密，以理和事实服人。驳论性毕业论文是指通过反驳别人的论点来树立自己的论点和主张。如果毕业论文侧重于以驳论为主，批驳某些错误的观点、见解、理论，就属于驳论性毕业论文。驳论文除按立论文对论点、论据、论证的要求以外，还要求针锋相对，据理力争。 按研究问题的大小不同可以把毕业论文分为宏观论文和微观论文。凡届国家全局性、带有普遍性并对局部工作有一定指导意义的论文，称为宏观论文。它研究的面比较宽广，具有较大范围的影响。反之，研究局部性、具体问题的论文，是微观论文。它对具体工作有指导意义，影响的面窄一些。 另外还有一种综合型的分类方法，即把毕业论文分为专题型、论辩型、综述型和综合型四大类： 1．专题型论文。这是分析前人研究成果的基础上，以直接论述的形式发表见解，从正面提出某学科中某一学术问题的一种论文。如本书第十二章例文中的《浅析领导者突出工作重点的方法与艺术》一文，从正面论述了突出重点的工作方法的意义、方法和原则，它表明了作者对突出工作重点方法的肯定和理解。2．论辩型论文。这是针对他人在某学科中某一学术问题的见解，凭借充分的论据，着重揭露其不足或错误之处，通过论辩形式来发表见解的一种论文。3．综述型论文。这是在归纳、总结前人或今人对某学科中某一学术问题已有研究成果的基础上，加以介绍或评论，从而发表自己见解的一种论文。4．综合型论文。这是一种将综述型和论辩型两种形式有机结合起来写成的一种论文。如《关于中国民族关系史上的几个问题》一文既介绍了研究民族关系史的现状，又提出了几个值得研究的问题。因此，它是一篇综合型的论文。 写作步骤：毕业论文是高等教育自学考试本科专业应考者完成本科阶段学业的最后一个环节，它是应考者的 总结 性独立作业，目的在于总结学习专业的成果，培养综合运用所学知识解决实际 问题 的能力。从文体而言，它也是对某一专业领域的现实问题或 理论 问题进行 科学 研究 探索的具有一定意义的论说文。完成毕业论文的撰写可以分两个步骤，即选择课题和研究课题。 首先是选择课题。选题是论文撰写成败的关键。因为，选题是毕业论文撰写的第一步，它实际上就是确定“写什么”的问题，亦即确定科学研究的方向。如果“写什么”不明确，“怎么写”就无从谈起。 教育部自学考试办公室有关对毕业论文选题的途径和要求是“为鼓励理论与工作实践结合，应考者可结合本单位或本人从事的工作提出论文题目，报主考学校审查同意后确立。也可由主考学校公布论文题目，由应考者选择。毕业论文的总体要求应与普通全日制高等学校相一致，做到通过论文写作和答辩考核，检验应考者综合运用专业知识的能力”。但不管考生是自己任意选择课题，还是在主考院校公布的指定课题中选择课题，都要坚持选择有科学价值和现实意义的、切实可行的课题。选好课题是毕业论文成功的一半。 第一、要坚持选择有科学价值和现实意义的课题。科学研究的目的是为了更好地认识世界、改造世界，以推动社会的不断进步和发展 。因此，毕业论文的选题，必须紧密结合社会主义物质文明和精神文明建设的需要，以促进科学事业发展和解决现实存在问题作为出发点和落脚点。选题要符合科学研究的正确方向，要具有新颖性，有创新、有理论价值和现实的指导意义或推动作用，一项毫无意义的研究，即使花很大的精力，表达再完善，也将没有丝毫价值。具体地说，考生可从以下三个方面来选题。首先，要从现实的弊端中选题，学习了专业知识，不能仅停留在书本上和理论上，还要下一番功夫，理论联系实际，用已掌握的专业知识，去寻找和解决工作实践中急待解决的问题。其次，要从寻找科学研究的空白处和边缘领域中选题，科学研究。还有许多没有被开垦的处女地，还有许多缺陷和空白，这些都需要填补。应考者应有独特的眼光和超前的意识去思索，去发现，去研究。最后，要从寻找前人研究的不足处和错误处选题，在前人已提出来的研究课题中，许多虽已有初步的研究成果，但随着社会的不断发展，还有待于丰富、完整和发展，这种补充性或纠正性的研究课题，也是有科学价值和现实指导意义的。 第二、要根据自己的能力选择切实可行的课题。毕业论文的写作是一种创造性劳动，不但要有考生个人的见解和主张，同时还需要具备一定的客观条件。由于考生个人的主观、客观条件都是各不相同的，因此在选题时，还应结合自己的特长、兴趣及所具备的客观条件来选题。具体地说，考生可从以下三个方面来综合考虑。首先，要有充足的资料来源。“巧妇难为无米之炊”，在缺少资料的情况下，是很难写出高质量的论文的。选择一个具有丰富资料来源的课题，对课题深入研究与开展很有帮助。其次，要有浓厚的研究兴趣，选择自己感兴趣的课题，可以激发自己研究的热情，调动自己的主动性和积极性，能够以专心、细心、恒心和耐心的积极心态去完成。最后，要能结合发挥自己的业务专长，每个考生无论能力水平高低，工作岗位如何，都有自己的业务专长，选择那些能结合自己工作、发挥自己业务专长的课题，对顺利完成课题的研究大有益处。 致 谢 这次论文的完成，不止是我自己的努力，同时也有老师的指导，同学的帮助，以及那些无私奉献的前辈，正所谓你知道的越多的时候你才发现你知道的越少，通过这次论文，我想我成长了很多，不只是磨练了我的知识厚度，也使我更加确定了我今后的目标：为今后的计算机事业奋斗。在此我要感谢我的指导老师——***老师，感谢您的指导，才让我有了今天这篇论文，您不仅是我的论文导师，也是我人生的导师，谢谢您！我还要感谢我的同学，四年的相处，虽然我未必记得住每分每秒，但是我记得每一个有你们的精彩瞬间，我相信通过大学的历练，我们都已经长大，变成一个有担当，有能力的新时代青年，感谢你们的陪伴，感谢有你们，这篇论文也有你们的功劳，我想毕业不是我们的相处的结束，它是我们更好相处的开头，祝福你们！我也要感谢父母，这是他们给我的，所有的一切；感谢母校，尽管您不以我为荣，但我一直会以我是一名农大人为荣。 通过这次毕业设计，我学习了很多新知识，也对很多以前的东西有了更深的记忆与理解。漫漫求学路，过程很快乐。我要感谢信息与管理科学学院的老师，我从他们那里学到了许多珍贵的知识和做人处事的道理，以及科学严谨的学术态度，令我受益良多。同时还要感谢学院给了我一个可以认真学习，天天向上的学习环境和机会。 即将结束*大学习生活，我感谢****大学提供了一次在**大接受教育的机会，感谢院校老师的无私教导。感谢各位老师审阅我的论文。 � EMBED Word.Picture.8 \* MERGEFORMAT ��� 宁波理工学院 � EMBED Word.Picture.8 \* MERGEFORMAT ��� 宁波理工学院 PAGE _1234567955.unknown _1234568019.unknown _1234568051.unknown _1234568067.unknown _1234568075.unknown _1234568083.unknown _1234568087.unknown _1234568089.unknown _1234568091.unknown _1234568093.unknown _1234568094.unknown _1234568092.unknown _1234568090.unknown _1234568088.unknown _1234568085.unknown _1234568086.unknown _1234568084.unknown _1234568079.unknown _1234568081.unknown _1234568082.unknown _1234568080.unknown _1234568077.unknown _1234568078.unknown _1234568076.unknown _1234568071.unknown _1234568073.unknown _1234568074.unknown _1234568072.unknown _1234568069.unknown _1234568070.unknown _1234568068.unknown _1234568059.unknown _1234568063.unknown _1234568065.unknown _1234568066.unknown _1234568064.unknown _1234568061.unknown _1234568062.unknown _1234568060.unknown _1234568055.unknown _1234568057.unknown _1234568058.unknown _1234568056.unknown _1234568053.unknown _1234568054.unknown _1234568052.unknown _1234568035.unknown _1234568043.unknown _1234568047.unknown _1234568049.unknown _1234568050.unknown _1234568048.unknown _1234568045.unknown _1234568046.unknown _1234568044.unknown _1234568039.unknown _1234568041.unknown _1234568042.unknown _1234568040.unknown _1234568037.unknown _1234568038.unknown _1234568036.unknown _1234568027.unknown _1234568031.unknown _1234568033.unknown _1234568034.unknown _1234568032.unknown _1234568029.unknown _1234568030.unknown _1234568028.unknown _1234568023.unknown _1234568025.unknown _1234568026.unknown _1234568024.unknown _1234568021.unknown _1234568022.unknown _1234568020.unknown _1234567987.unknown _1234568003.unknown _1234568011.unknown _1234568015.unknown _1234568017.unknown _1234568018.unknown _1234568016.unknown _1234568013.unknown _1234568014.unknown _1234568012.unknown _1234568007.unknown _1234568009.unknown _1234568010.unknown _1234568008.unknown _1234568005.unknown _1234568006.unknown _1234568004.unknown _1234567995.unknown _1234567999.unknown _1234568001.unknown _1234568002.unknown _1234568000.unknown _1234567997.unknown _1234567998.unknown _1234567996.unknown _1234567991.unknown _1234567993.unknown _1234567994.unknown _1234567992.unknown _1234567989.unknown _1234567990.unknown _1234567988.unknown _1234567971.unknown _1234567979.unknown _1234567983.unknown _1234567985.unknown _1234567986.unknown _1234567984.unknown _1234567981.unknown _1234567982.unknown _1234567980.unknown _1234567975.unknown _1234567977.unknown _1234567978.unknown _1234567976.unknown _1234567973.unknown _1234567974.unknown _1234567972.unknown _1234567963.unknown _1234567967.unknown _1234567969.unknown _1234567970.unknown _1234567968.unknown _1234567965.unknown _1234567966.unknown _1234567964.unknown _1234567959.unknown _1234567961.unknown _1234567962.unknown _1234567960.unknown _1234567957.unknown _1234567958.unknown _1234567956.unknown _1234567923.unknown _1234567939.unknown _1234567947.unknown _1234567951.unknown _1234567953.unknown _1234567954.unknown _1234567952.unknown _1234567949.unknown _1234567950.unknown _1234567948.unknown _1234567943.unknown _1234567945.unknown _1234567946.unknown _1234567944.unknown _1234567941.unknown _1234567942.unknown _1234567940.unknown _1234567931.unknown _1234567935.unknown _1234567937.unknown _1234567938.unknown _1234567936.unknown _1234567933.unknown _1234567934.unknown _1234567932.unknown _1234567927.unknown _1234567929.unknown _1234567930.unknown _1234567928.unknown _1234567925.unknown _1234567926.unknown _1234567924.unknown _1234567907.unknown _1234567915.unknown _1234567919.unknown _1234567921.unknown _1234567922.unknown _1234567920.unknown _1234567917.unknown _1234567918.unknown _1234567916.unknown _1234567911.unknown _1234567913.unknown _1234567914.unknown _1234567912.unknown _1234567909.unknown _1234567910.unknown _1234567908.unknown _1234567899.unknown _1234567903.unknown _1234567905.unknown _1234567906.unknown _1234567904.unknown _1234567901.unknown _1234567902.unknown _1234567900.unknown _1234567895.unknown _1234567897.unknown _1234567898.unknown _1234567896.unknown _1234567893.unknown _1234567894.unknown _1234567892.unknown _1234567891.unknown _1234567890.doc