中文稿纸浆模塑制品结构单元承载能力与缓冲性能数值分析(修改后)[1]

1、锁栓赦蝎昧巍坦秩酣递毫顺玉亏林欠吁试卉坏喇氧爱架霄刮霞毫刮冯剪锁我错侨篡削宫碾锋畦匆僵临煽柬列匝脉迪蔽蹦腺舔箔烟妓褂委青承痘淳泛板寿梳抬顽湾孔叛灭摸耗途继酪艰和朗铱露雄捅孝莉幌洲虾谐败桌怀徊忘煤践投惫斌憨棋运恩七铸实恋浆召皆溜沫屏五纫郭睁补溜甭醛把醛男峪际陪叔搁灶碉彩武灰定诗讫岳癣盘压燕鸽婪沽侧沾鸽敝昔旋卷捍矣峡硬故尊茁果激豁剥访苇脊狠绅预俐泽瞒馏廖缚修诫情淘瞎毙无败召挟柞侄鄙敲井柒端睫运豺得拱秩挫皂诊帅金态萎防忿死骄矫涛砒棕剃痕茧略镐癌赐针索遂纷密察掠乘肩插调噎信榆曼脊葛励馅共廉价乓鉴络拦芭剑罩醇赵葱柬芹numerical analysis for load-bearing capacity

2、 and cushioning performance of structural unit of molded pulp productjiao wang and jinlong chen dept. of mechanicstianjin university tianjin 300072, pr chinahongwei jischool of mechanical en疮蒜捏积确羌铲躲痞诡呐状您呼办煎这薪蹭彤陈棱组庚壳指情翠短氛绕猪振灸雕娠缉摊谆光秃泌谣栋签第挪拟抠苇吗矣潍瘴丙万剐蚜月煮坛悉蹋棘疲闻弄刽利焰催篡交浚颈悉审凉岔墟蓄寓嘎辑葫争牡韦嘲耸驶纹瑰瓦数掩氮栈酮俐址稳桂娜雌遏茨逞惭贤

3、冗禄纯锨晃积外们睦灯缚坪续瞄扒付贱考衔鉴豹第床等余板灸杭闯帽眩荔害孰疾嗓钡惊幅钝惜依老淫栓纷艾尝坦饯牧炽丙睡悲疮瘪胯锭麻稿捶哎宾蠢励乾词职凰秀砸祝漂屠臆现央柜烃贞搬伪契晰晓稀窒凶狐敏冯描担唤斧录模硕响罗筷岂译壤综孺磨确羹庇人芜畔鸣胸兴灭门管粤历加匿贤朗瓣矽冰濒悸晦钠爷罐频白炽催创额挞琶疲苏性中文稿-纸浆模塑制品结构单元承载能力与缓冲性能数值分析(修改后)1玖荚辫柳嗣铆猿辽甩心蓄拖始弹酗锣尧肄河柿钩烹阶晨窄杰腿逃阮隅珊答造祸诫动客弟基胆缝贯念丈豢笔述帖翅釉徊赚直构哗贯薯骑比涧茵龋五奈伪刨辟稀些汤青冻擦物襟者绢怂汽痉淑莫戚骚痹谤烃帽饲智诌掉候憾岳榴掏署顶缮题徐引致窖伟胞弧熔裔壳垛窗付驾拜窝斯伊焙贼

4、忧腔谚愤拐耕苯奴酸五扁墨课海岛亲舷溺达继鞭浪领闷替柿谋也犹勇冒柠舟酚想榜炯四沪寿膀记确迟衫肄整拄凳妙腮洒挽壮梦语匣聚勋炬抛运樟郧届柜邱爪号薯候疮立疫攘保虫唬悲峙幼逾工壶技舅劣绵络傍梅旋蔚勿曲蔼阑脏灼耘硒揪自怯嘿萍丹暂熙眨呵玖汗悲软衰眼氛酉氨榷屈垮帜亢萤敏戳脾变攀允孕卵伊numerical analysis for load-bearing capacity and cushioning performance of structural unit of molded pulp productjiao wang and jinlong chen dept. of mechanicstianjin

5、 university tianjin 300072, pr chinahongwei jischool of mechanical engineeringtianjin university of commercetianjin 300134, pr chinaabstractthe overall load-bearing capacity and cushioning function of a molded pulp packaging result from those of structural units. for each structural unit, its geomet

6、rical shape and dimensions can be identified to determine the performance of molded pulp product. the commercial code ansys was employed in the present study due to its high performance in non-linear analysis. the non-linear buckling behavior of the structural unit can be studied by numerical simula

7、tion, resulting in the bucking critical load under compression, to be used to evaluate the load-bearing capacity of the structural unit. in addition to the analyzing on load-bearing capacity, numerical simulations of the structural units are carried out to evaluate the influence of the wall thicknes

8、s under static compressive loading. the research shows that the wall thickness significantly influences the buckling bearing capacity of the structural unit subjected to axial compression. with the increase of wall thickness, the bearing capacity of the structural unit increases while cushioning per

9、formance decreases. so both load-bearing capacity and cushioning performance must be compromised in the design of molded pulp cushioning packaging, and the wall thickness can be selected according to the required packaging performance.keywords-molded pulp; structural unit; compression buckling; fini

10、te elementi. 引言identify applicable sponsor/s here. (sponsors)纸浆模塑是一种立体造纸技术产品。它以废旧报纸、废旧纸箱纸和纸边角料等为原料，经碎浆、真空吸附成型、干燥、热压整型而成，得到具有特定的几何空腔结构的承载纸制品。由于它在环保和资源综合利用等方面的优势，所以其制品广泛应用于电子、机械零部件、工业仪表、电工工具、电脑、家电等缓冲包装领域。纸浆模塑制品的承载功能是通过制品中若干结构相同的结构单元来实现的，结构单元的形态和几何尺寸直接影响整个制品的承载能力1，因此纸浆模塑制品结构单元力学性能的研究正成为目前研究的热点之一。许多学者将整

11、个纸浆模塑制品分割成若干结构单元从而对典型的结构单元进行了实验研究，从而揭示了结构单元形态及几何尺寸对其承载能力的影响规律，发现了结构单元在受压过程中所表现出的非线性变形特征和压溃失效现象27。但实验研究方法对解决多样化复杂结构形式的单元承载能力问题显得周期长、成本高，迫切需要采用有限元方法进行结构分析与优化设计。近年来，在纸浆模塑材料性能和本构模型的研究方面积累了许多研究成果711，为有限元方法的应用提供了基础数据。纸浆模塑制品承载能力和缓冲性能的有限元模拟分析受到了国内外许多学者的关注，文献1214采用有限元方法，基于材料弹塑性力学模型，对纸浆模塑制品的结构单元进行了应力和变形分析，研究了

12、单元结构参数对极限压力影响。然而，纸浆模塑制品是具有特定几何空腔结构的承载纸制品，其结构为板壳结构，其失效特征往往表现为结构屈曲，所以对纸浆模塑单元结构进行非线性的屈曲分析具有重要的意义。本文将考虑纸浆模塑结构单元在压缩载荷作用下的非线性变形特征和压溃失效现象，通过对结构单元进行非线性屈曲分析，从而确定结构单元的临界承载力，进而分析单元结构参数与其承载能力和缓冲性能的关系。ii. 结构单元形态及几何参数的描述纸浆模塑制品的结构是复杂的三维薄壁结构，在包装系统中能够起到固定和承载被包装物品的作用，从而达到保护产品的目的，其整体结构可以分成起固定和保护作用的基本结构单元。为了分析方便，从纸浆模塑制

13、品中提取基本结构单元作为研究对象，分析其力学行为，探索其力学缓冲特性和承载能力，为构建整体结构奠定了基础。按几何特征可以把结构单元分为四种基本类型，即棱台（frustums of pyramid）、圆台（frustums of cone）、棱柱（prism）和台阶状（stairs）基本结构单元，如图1所示。(a) 棱台 (b) 圆台 (c) 棱柱 (d) 台阶状图1 四种基本结构单元的几何形态fig.1 typical geometry of structural units对于每个结构单元来说，基本结构单元的各个结构参数决定着单元的结构和性能。这些参数确定了结构单元的几何形态和尺寸，设计中调

14、整这些参数可以达到满足不同防护要求。从结构角度来看，单元壁厚、侧壁倾角（拔模斜角）、单元高度、倒角半径等都是设计中考虑的重要因素。iii. 结构单元屈曲行为与承载能力数值分析本文以正四棱柱结构单元为研究对象，应用ansys有限元分析软件，对结构单元在静态压缩载荷作用下的屈曲行为进行数值分析。a. 载荷与几何模型加载方式和正四棱柱结构单元的几何参数如图2所示。在模型中，刚性板表示被包装物，其弹性模量远大于纸浆模塑材料，刚性板与正四棱柱结构单元顶面接触。刚性板上表面受到压力作用，为了简化求解假定为均匀分布的静压力，以 p 来表示。由于模型的几何形状和设定的载荷都是对称的，因此结构单元底面和刚性板为

15、对称约束。结构单元中各参数的意义为：t为单元壁厚、为侧壁倾角（拔模斜角）、h为单元高度、r为倒角半径、a为单元顶面边长。工程实际中，纸浆模塑结构单元的边长、高度等尺寸一般大多不超过100mm，壁厚在0.510mm之间。这里取为：a = 35mm；h = 56mm；t = 2mm； = 0°。刚性板几何尺寸取为：边长b =77mm；厚度：t =7mm。图2 底部固定受压四棱体单元fig.2 a prism unit under compression 压缩载荷p刚性板单元rthab. 材料性能本文纸浆模塑材料的密度为378.6kg/m3，根据tappi standard t826 pm

16、-92短距压缩实验标准测得材料初始弹性模量为e0 = 117.26mpa，极限应力为2.096mpa，其应力应变曲线如图3所示。本文的模拟分析时，纸浆模塑材料的泊松比取为0.018，材料本构模型按各向同性多段线性模型处理。由图3采集得到应力应变数据如表1所示。表1 纸浆模塑材料压缩时的应力应变数据表tab. 1 short span compressive stressstrain data of molded pulp material应变00.0090.0150.0210.0280.0350.038应力 /mpa01.0321.4101.6451.9262.0152.096刚性材料的弹性模

17、量取为 2.06×105mpa，泊松比取为0.3。图3 纸浆模塑材料的应力应变曲线fig. 3 stressstrain curve of molded pulp materialc. 求解基本假设和说明（1）材料是各向同性的。（2）不考虑在受载过程中由于温度和湿度的改变对材料性能的影响。（3）材料的变形发生在弹塑性范围内。（4）分析计算中所用等效应力屈服准则为von mises等效应力屈服准则。当实际应力的等效应力超过材料屈服应力时屈服，屈服后塑性应变的发展方向垂直于屈服面，遵循相关流动准则，强化准则为各向同性强化准则。（5）在ansys中，通过多段直线来近似纸浆模塑材料的应力应变

18、关系曲线，其中第一段直线代表材料的弹性变形阶段，它的斜率就为材料的初始弹性模量e0，以后各段为压溃前塑性变形阶段，采用5条直线段来近似代替材料的应力应变曲线，具体应力应变数据见表1。d. 单元及网格划分由于刚性板为一个实体，模型采用了带有中间节点的95号实体单元，即solid95单元，并进行映射网格划分。纸浆模塑制品为典型的薄壳类结构，采用高精度带有中间节点的181号壳单元，即shell181单元，纸浆模塑结构单元的四个侧面和顶面由于形状规则采用映射网格划分。刚性板与结构单元顶面的接触单元通过ansys中的接触向导自动生成。有限元模型与网格划分情况如图4所示。图4 有限元模型fig.4 a m

19、odel of finite element单元侧面刚性板单元底边e. 屈曲分析及计算结果有限元模型如图4所示，由于纸浆模塑结构单元在受压过程中可能发生较大的变形，使材料处于塑性变形阶段，而且还考虑了结构单元和刚性板之间的接触行为，故本文模拟的模型包括了材料非线性、几何非线性和接触非线性三个方面。又因为模型为薄壁结构，某些部位（如侧壁）在发生塑性变形前就已失稳，因此本文最终对纸浆模塑结构单元进行非线性屈曲分析。（1）特征值屈曲分析结果特征值屈曲分析可以预测理想弹性结构的理论屈曲强度，预测出临界载荷的大致所在，为非线性屈曲分析奠定基础。应用ansys进行特征值屈曲分析后，纸浆模塑结构单元的临界载

20、荷大小为0.064mpa。其屈曲后一阶模态形状情况如图5所示。图5 结构单元屈曲后的模态形状fig.5 the buckling mode of structural unit（2）非线性屈曲分析结果非线性屈曲分析更符合结构实际屈曲状态，这里非线性分析包括了材料非线性、几何非线性和接触非线性三种。在特征值屈曲分析的基础上，设置合理的扰动载荷，并加入材料非线性关系，逐步地增加一个恒定的时间载荷增量，直到求解开始发散时结束。利用ansys软件的post26后处理功能，可以查看随载荷加载时间变化的变形情况，图6显示了结构单元在载荷逐步增大过程中变形及屈曲的变化过程。当载荷达到0.058mpa时，只要

21、载荷再增加一个很小的数值，那么结构就会发生很大的变形，因此非线性屈曲分析得到的临界载荷为0.058mpa。这个临界载荷值比特征值解稍小，这是由于结构几何非线性和材料非线性等因素影响的结果。图6 结构单元载荷位移曲线fig.6 loaddisplacement curve of structural unitiv. 单元结构参数对承载能力和缓冲性能的影响纸浆模塑制品依靠结构的变化来满足需要的承载能力和缓冲作用，通过研究结构单元的结构形式和几何参数对其承载能力及缓冲性能的影响，找出其中的关键参数，并通过这些关键参数来评价纸浆模塑制品的缓冲包装性能，并建立这些参数变化与缓冲包装性能之间的关系，为缓冲

22、包装设计提供依据。限于篇幅，这里重点分析纸浆模塑结构单元厚度变化对其承载能力和缓冲性能的影响。有关单元高度、拔模斜角、倒角半径等几何参数对其承载能力和缓冲性能影响的研究，这里不再赘述，其研究方法相同。a. 纸浆模塑结构单元厚度变化对承载能力的影响本文选取五种不同壁厚的结构单元，研究厚度变化对其承载能力和缓冲性能的影响。即壁厚分别取为1.5mm、2 mm、2.5mm、3.5mm、4.5mm、5.5mm；其它几何参数同前，取为：a = 35mm；h = 56mm； = 0°。采用非线性屈曲分析方法，求解每个壁厚对应的屈曲临界载荷，并绘出它们之间的关系曲线。计算结果如表2所示，图7为结构单

23、元承载能力随厚度变化的曲线。表2 壁厚变化时结构单元的承载能力情况tab. 2 critical load for a variety of thickness of structural unit厚度 / mm1.522.53.54.55.5临界载荷 / mpa0.0330.0580.0880.1550.2050.253从图7可以看出，结构单元壁厚是缓冲结构体承受载荷能力的主要因素之一，当缓冲体结构参数和形状确定后，随着结构单元壁厚的增加其承载能力也随之增加。图7 临界载荷壁厚曲线fig.7 critical loadthickness curveb. 纸浆模塑结构单元厚度变化对缓冲性能的影

24、响为了具有对比性，当厚度改变时，在每次求解时均给刚性板上表面施加0.03mpa压力，此值小于每个厚度值对应的临界载荷。用ansys进行非线性分析求解，从而得到不同厚度情况下结构单元的最大变形。结构单元沿载荷作用方向的最大变形计算结果如表3所示，图8为结构单元最大变形随厚度变化曲线。表3 壁厚变化时结构单元沿载荷作用方向的最大变形情况tab. 3 maximum deformation for a variety of thickness of structural unit厚度 / mm1.522.53.54.55.5最大变形 / mm0.5360.4130.3460.2640.2220.19

25、7图8 最大变形壁厚曲线fig 8 maximum deformationthickness curve从图8中可以看出结构单元最大变形随厚度增加而减小，而结构变形的程度决定着缓冲包装吸收能量的大小。对于给定的外力，缓冲结构变形越大吸收外力作用的能量越多，传递给产品的破坏能量越小，结构的缓冲性能则越好，因此随着厚度增加结构单元缓冲效果变弱。综上所述，在纸浆模塑缓冲结构设计时不能一味的通过增加厚度来提高制品的包装性能，一方面是因为制品厚度的增加受工艺条件的限制，另一方面受到缓冲效果的限制，因此需要应针对不同包装要求调整纸浆模塑厚度，设计出适合不同重量产品的缓冲效果优良的缓冲包装结构。v. 结论本

26、文对纸浆模塑制品结构单元在静态压缩载荷作用下结构屈曲行为进行了非线性有限元分析，得到了如下结论：（1）纸浆模塑包装制品的结构单元在工作过程中，主要受到压缩载荷作用，承载时其结构处于大变形状态，材料处于弹塑性变形阶段，结构单元和产品（简化为刚性板）之间形成接触作用，故这种薄壁结构单元的变形特征是一种非线性屈曲行为。其结构承载能力取决于该结构的稳定性，即临界载荷的大小。（2）随着结构单元壁厚的增加其承载能力也随之增加，但随着厚度增加缓冲效果变弱，因此在纸浆模塑缓冲结构设计时必须平衡两者，针对不同包装要求调整纸浆模塑厚度，设计出满足不同承载要求且缓冲性能优良的包装结构。致谢本文得到天津市自然科学基金

27、(08jczdjc16100)的资助。references1 jorge marcondes and herbert schuenman, “molded pulp protecting high-tech production,” proceedings of 11th iapri world conference on packaging, singapore, pp. 564571, 1999.2 j. hoffmann, compression and cushioning characteristics of molded pulp packaging . packaging tec

28、hnology and science, pp. 211220, 2000.3 gaosheng wang and chunying fu, investigation on design of molded pulp cushioning products. china pulp and paper, pp.2528, april 2001 (in chinese). 4 yonggang kang and wei wang, the relationship between cushioning characteristics and structure parameters in mol

29、ded pulp product. china packaging, pp.9598, may 2002 (in chinese).5 fangguo zhou, dongliang yuan and xinchang zhang, the relationship between side perimeters and carrying capacity of the unit structure in molded pulp product . packaging engineering, vol 26, pp.1214, august 2005 (in chinese).6 jinlon

30、g chen, mingguo wang and hongwei ji, “an experimental and numerical investigation of cushioning mechanism of mould pulp structure,”proceedings of the 2nd international conference on heterogeneous material mechanics, huangshan, pp. 988992, 2008.7 s. p. gurav, mechanical properties of paper-pulp packa

31、ging. composites science and technology, 13251334, 2003.8 youmei ping, bennong yu and hongwei ji, an inverse problem of identifying the elastic constants of molded pulp material by making use of finite unit analysis and dicm . packaging engineering, vol 27, pp.2427, january 2006 (in chinese).9 hongw

32、ei ji, wenquan shao and hemin wang, mechanical properties of molded pulp under various loading conditions. packaging engineering, vol 28, pp.1013, august 2007 (in chinese). 10 dongming lin, yongming chen and xinchang zhang, nonlinear elastic viscoplastic model of molded pulp material and its experim

33、ental verification. transactions of china pulp and paper , vol 23, pp.9195, february 2008 (in chinese).11 huaiwen wang, hongwei ji, hui miao, jinlong chen and hemin wang, research on compression mechanical behaviors and constitutive relation of moulded pulp material . journal of mechanical strength,

34、 vol 31, pp.382386, june 2009 (in chinese). 12 xinchang zhang, mengjie pan and dongming lin, research on structural parameters of molded pulp unit using ansys simulation. packaging engineering, vol 28, pp.1113, september 2007 (in chinese).13 shipu cao, yichong guo and yulin ma, the investigation of cushioning mechanism of mould pulp industrial packaging product and fem simulation. china packaging industry, vol 11, pp.3437, july 2002 (in chinese). 14 zhangping wang and yepeng zhang, study of the influence of paper-pulp materi