基于ANSYS的秸秆类生物质冷成型仿真分析.pdf
《基于ANSYS的秸秆类生物质冷成型仿真分析.pdf》由会员分享,可在线阅读,更多相关《基于ANSYS的秸秆类生物质冷成型仿真分析.pdf(7页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、2009 年1 2 月农业机械学报第4 0 卷第1 2 期基于A N S Y S 的秸秆类生物质冷成型仿真分析孙启新1张仁俭1董玉平2(1 淮海工学院机械工程学院,连云港2 2 2 0 0 5;2 高效洁净机械制造教育部重点实验室(山东大学),济南2 5 0 0 6 1)【摘要】运用A N S Y S 有限元理论对秸秆类生物质成型徐变过程进行模拟分析,根据其生物机理和化学成分建立生物质徐变有限元模型。生成接触对,通过非线性大变形分析得到成型徐变规律和成型压力一变形关系。结果表明生物质冷成型徐变过程包括松散、压紧、固化3 个阶段,是塑性变形和粘性变形的结合。关键词:生物质冷成型A N S Y S
2、 模拟分析中图分类号:T K l 6文献标识码:AS i m u l a t i o nA n a l y s i so fC o m p r e s s i n gM o l d i n gu n d e rG e n e r a lC o n d i t i o nf o rS t r a wB i o m a s sB a s e do nA N S Y SS u nQ i x i n lZ h a n gR e n j i a n lD o n gY u p i n g z(1 C o l l e g eo f M e c h a n i c a lE n g i n e e r i
3、n g,H u a i h a iI n s t i t u t eo fT e c h n o l o g y,L i a n y u n g a n g2 2 2 0 0 5,C h i n a2 K e yL a b o r a t o r y0,H i g hE f f i c i e n c ya n dC l e a nM e c h a n i c a lM a n u f a c t u r e(S h a n d o n gU n i v e r s i t y),M i n i s t r yo fE d u c a t i o n,J i n a n2 5 0 0 6 1
4、,C h i n a)A b s t r a c tS i m u l a t i o na n a l y s i so fm o l d i n gp r o c e s sf o rs t r a wb i o m a s sw a sc a r r i e db a s e do nA N S Y S T h em e c h a n i c sm o l dw a se s t a b l i s h e da c c o r d i n gt oi t sb i o l o g ym e c h a n i s ma n dc h e m i c a lc o m p o n e
5、n t,a l o n gw i t hc o n t a c tp a i r s T h e nm o l d i n gr u l e sa n dr e l a t i o n sb e t w e e np r e s s u r ea n dd e f o r m a t i o nw e r eo b t a i n e dt h r o u g hn o n l i n e a r i t ya n a l y s i s T h ea n a l y s i sr e s u l t ss h o w e dt h a tb i o m a s sm o l d i n gu
6、 n d e rg e n e r a lc o n d i t i o np r o c e s si n c l u d e dt h r e es t a g e s:l o o s e,c o m p a c t i o n,s o l i d i f i c a t i o n,w h i c hc o m b i n e dt h ep l a s t i cd e f o r m a t i o nw i t ht h ev i s c o u sd e f o r m a t i o n T h eo b t a i n e dr e s u l t sa r ew e l l
7、i na g r e e m e n tw i t he x p e r i m e n td a t a K e yw o r d sB i o m a s s,M o l d i n gu n d e rg e n e r a lc o n d i t i o n,A N S Y S,S i m u l a t i o n,A n a l y s i s引言生物质压缩成型是将生物质废弃物用机械加压的方法,使原来松散、无定形的原料压缩成具有一定形状、密度较大的固体成型燃料。生物质在经过压缩成型之后,其密度、强度和燃烧性能都有了质的改善,大大提高了生物质作为燃料的品质和性能,可以为日常生活、工
8、业生产提供高效清洁的能源。我国生物质能资源非常丰富。但是,作为一种散抛型低容积密度的能源存在形式,生物质能源具有资源分散、能量密度低、容积密度小、储运不方便等缺点。特别是我国地域广阔,地形复杂。严重制约了生物质能的大规模应用。所以生物质高品位转换技术的研究便成为开发利用生物质能的重点。而研制开发成本低、产能高、操作简单、体积小、便于田间小户型作业的成型设备成为这项技术推广应用的关键所在。近年来在生物质压缩成型技术的改进创新发展基础上,本文利用有限元理论对生物质的压缩成型徐变过程和力学特性进行模拟分析,对成型过程中压力和变形的变化规律作深入分析,为研究其最佳成型条件和工艺流程寻求理论依据。l 生
9、物质致密成型技术1 1 基本理论从生物学角度。生物质是由纤维素、半纤维和木质素构成。这3 种成分构成植物的支持骨架。其中,收稿日期:2 0 0 9 0 5 1 0 修回日期:2 0 0 9 0 7 1 1作者简介:孙启新,讲师主要从事机械设计及生物质能源研究。E m a i l:s u n q i x i n s d 1 6 3 o D m通讯作者:董玉平。教授。博士生导师。主要从事机械设计理论及生物质能源研究,E m a i l:d o n g y p s d u e d u 万方数据第1 2 期孙启新等:基于A N S Y S 的秸秆类生物质冷成型仿真分析1 3 1纤维素组成微细纤维,构成
10、纤维细胞壁的网状骨架,赋予植物茎秆弹性和机械强度;而半纤维素和木质素则是填充在纤维之间和微细胞之间的“粘结剂”和“填充剂”。在一般植物中,这3 种成分的质量占总质量的8 0 9 0。是生物质的主要组成成分。生物质典型结构如图1 所示【1 q J。2)图1 生物质典型结构F i g 1T y p i c a ls t r u c t u r eo fb i o m a s s1 木质素2 纤维素3 半纤维素根据国内外学者对粘结成型机理的研究,可以将成型块内部粘结状况及结合力分成5 种类型:固体桥;非自由移动粘结剂下的粘附力和内聚力;自由移动液体的表面张力和毛细管作用力;固体颗粒之间的吸引力(分子
11、间的范德华力或静电力);颗粒间的交错粘结。除了在特殊的成型过程中,物料压缩成型利用了部分或全部的粘接机理以外,通常成型块在聚合过程中只有一种或两种粘接机理起着主要的作用L3|。一般情况,成型块的粘结强度随着成型压力增加而增大。在不添加粘结剂的成型过程中,秸秆颗粒在外部压缩力的作用下相互滑动,颗粒间的孔隙减小,颗粒在压力作用下发生塑性变形,并达到粘结成型的目的。对大颗粒而言,颗粒之间交错粘接为主;对于很小的颗粒而言(粉粒状),颗粒之间以吸引力(分子间的范德华力或静电力)粘结为主。因此生物质徐变应包括粘性变形和塑性变形。1 2 生物质徐变有限元模型基于以上分析,如何真实地模拟纤维的分裂软化,木质素
12、的固化成型是有限元分析的关键。同时颗粒物料应属于非连续离散性介质,而非连续介质力学的基本理论还很不完善,使其在工程上的应用受到一定的限制。为研究问题的方便而又能充分反映生物质受压变形情况,本文借用钢筋混凝土梁的粘塑性分析理论,将纤维等效为梁的纵横向筋,木质素等效为混凝土,并采用离散式模型。单元应力一位移关系为【4 J蠢=B 五(1)式中考单元应变率矢量五单元节点位移矢量县单元位移一应变转换矩阵,取决于选用的形式函数式中,上圆点表示对时间的导数,下同。结构物平衡方程为I。,B T d d V=素4(2)式中者单元应力率矢量R4 节点荷载率矢量单元内应力一应变关系,应反映出非弹性材料的特点,表示为
13、=D(壶一句)(3)式中D 单元材料本构矩阵面徐变应变率矢量徐变应变率是应力、时间和温度的函数。一般表达式为 5】由=f(口,t,T)(4)式中厂-单元徐变应变函数z 单元徐变时间T 单元徐变温度将式(3)代入式(2)得I。B T D(主一而)d V=蠢4(5)将式(1)代入式(5),移项整理得(Jy B T DB d V)u=】i 4+Jy B T D 而d V(6)式(6)可简化为K 五=R4+R _(7)其中K=l,罗T D B d V(8)素k 1,B T D d V(9)式中K 组装后的总体刚度矩阵R,初始徐变率矢量R4 节点荷载率矢量式(7)即为考虑徐变的有限元基本方程。对于非时效
14、材料K 由常量组成,对徐变计算的主要难度在于如何近似获得R 4,因为它是某时间段内变化着的应力函数。在有限元分析中采用L I N K 单元模拟纤维,c o n c r e t e 单元模拟木质素。同时为了较真实地虚拟成型过程,采用直径3 0m m、高6 0m m、壁厚1m m 的圆筒状模具,并使用P L A N E4 2 单元模拟成型模具。由于细胞构造、排列等方面的原因,植物的材料性质多为各向异性。如木材沿树干纵向(习惯称顺 万方数据1 3 2农业机械学报2 009 焦纹)强度较横向(横纹)大得多。因此,有限元计算时采用正交各向异性材料模型。由于秸秆类生物质力学性能参数的测定至今没有得到一个公
15、认的准确数值,且与杉木相比两者的主要组成成分都是纤维素,因此在计算分析中所需材料模型的弹性模量、泊松比等参数参考了文献 3 报道的杉木力学数据,应力一应变曲线参考了文献 6 9 的实验数据。2A N S Y S 分析计算由于成型挤压过程中被挤压物料和成型模具之间存在挤压和滑动摩擦现象,进行有限元分析时需建立接触副,用以计算两个物体之间的这种运动。2 1 网格划分首先为不同几何模型选择对应单元类型和材料模型;然后使用M e s h T o o l 进行智能网格划分。生物质模型和成型模具模型的单元尺寸分别为2 级和3 级。2 2 建立接触副A N S Y S 支持3 种接触方式 1 0-l z】:
16、点一点、点一面、面一面,每种接触方式使用的接触单元适用于某一类问题。对于刚体一柔体的面一面的接触单元,刚性面被当作“目标”面,分别用T a r g e l 6 9 和T a r g e l 7 0来模拟2 D 和3 一D 的“目标”面,柔性体的表面被当作“接触”面,用C o n t a l 7 1,C o n t a l 7 2,C o n t a l 7 3,C o n t a l 7 4 来模拟。本文采用目标单元1 6 9 和接触单元1 7 1 进行面一面接触分析计算。分别选取生物质物料作为柔性体,成型模具为刚性体,“目标”面分别选取两部分的边界线。2 3 建立分析条件生物质物料的左边采用
17、对称约束,在其顶部施加压力作为载荷。接触分析是典型的非线性问题,同时生物质也是非线性材料,所以在求解前。进行计算结果收敛性设定。首先,在S o l nc o n t r o l 中设定L a r g ed i s p l a c e m e n ts t a t i c,T i m ea te n do fl o a d s t e p 设定为1 0,N u m b e ro fs u n s t e p 设定为2 0;然后,在A n a l y s i so p t i o n 中设L a r g ed e f o r me f f e c t 为O n;最后打开非线性计算收敛预测器。2 4
18、 分析计算点击C u r r e n tL S 进行求解。2 4 1 挤压成型徐变过程图2、图3 是生物质徐变初始阶段和结束阶段的变形图,从图中可清晰看出其成型过程中位移、体积变化情况。物料变形具有一定的层次,从上到下其变形越来越小,顶端最大。底部变形较小。分析可得各部分受力、变形不均匀,采用模具成型得到的制成品各部分致密度不相同,一端较致密,另一端容易散落。另外,生物质边缘变形比中部变形有一定迟缓,这是由磨擦所造成的。所以选择受力、变形均匀的挤压方法是提高制品品质的关键因素。图2f=0 5s 时变形F i g 2D e f o r m a t i o no ft=0 58图3t=1 0s 时
19、变形F i g 3D e f o r m a t i o no ff=1 0s2 4 2 应力分布图4 为第三强度定律应力图,图5 为等效塑性应力图。图4t=1 0s 时第三强度定律应力图F i g 43 r dp r i n c i p a ls t r e s so ft=1 03从图4、图5 所示应力分布得出,生物质上部应力较大,下部较小,这与上部变形大于下部相一致。与材料模型相比较得出,物料上部徐变以粘性变形为主,而下部刚进入粘性变形,以塑性变形为主。对于成型模具其内壁直角处应力较大,其主要原因是这一部位即承受横向挤压力又承受上部挤压力和摩 万方数据第1 2 期孙启新等:基于A N S
20、 Y S 的秸秆类生物质冷成型仿真分析1 3 3图5t=1 0s 时等效塑性应力图F i g 5P l a s t i ce q u i v a l e n ts t 硝so ft=1 0s擦力在此产生的弯曲力偶。如何消除此集中应力是提高模具寿命的关键。2 4 3 压力和变形关系从图6 中得出,物料挤压初始阶段在较小压力作用下产生了较大变形,随着时间推移变形越来越缓慢,而且在成型后阶段变形会出现徐变停滞现象。根据图6、图7 可将生物质徐变过程分为3 个阶段:松散阶段、压紧阶段、固化阶段。燹彤m m图6 生物质压力一变形曲线F i g 6R e l a t i o nb e t w e e nl
21、 o a da n dd e f o r m a t i o no fb i o m a s s在松散阶段(图6、7 中曲线从0 点到A 点),初始物料较松散,在很小的压力下,产生较大变形。表明在此阶段,变形以克服物料之间的空隙为主,物料的体积减小较快。在压紧阶段(图6、7 中曲线从A 点到B 点),变形的增大率逐渐减弱。表明在此阶段,秸秆物料的变形从克服物料之间的空隙为主转向以克服物料自身内部的空隙为主。物料颗粒发生破裂,纤维开始分裂软化,徐变以塑性变形为主。在固化阶段(图6、7 中曲线从B 点到C 点),物董奋压缩嚣l U k g m-图7 生物质压力一压缩密度曲线F i g 7R e l
22、 a t i o nb e t w e e nl o a da n dd e n s i t yo fb i o m a s s料变形与压力的关系近似呈线性关系。随着压力的增加,变形的增加非常缓慢。表明在此阶段,秸秆物料颗粒之间的空隙和秸秆自身的空隙己经被完全克服,纤维进一步软化,木质素开始固化,徐变以粘性变形为主。在压紧和固化阶段,变形出现停滞,其原因一方面与物料中存在的大量气泡和水分有关,另一方面与大量纤维的突然软化分裂有关。在整个成型过程中。除去水分外物料质量损失较少,所以密度变化与位移变化曲线相似。从图6、7 可以得出:物料成型与其含水率有一定关系。本文有限元分析数据与文献 1 3 1
23、 4 中经过试验得出的秸秆类生物质压力一变形数据很相似。所以,本文所采用方法较为合理,分析计算准确。3 结论(1)准确建立生物质的徐变模型,模拟其生物机理是有限元分析合理性的关键。由生物质徐变图可得,对物料从各个方向施加均匀的受力是保证成型品品质的关键。(2)根据压力一变形图,生物质成型徐变要经过松散、压紧、固化3 个阶段。(3)由于挤压力和摩擦力的存在,成型模具产生较大应力集中,这是柱塞式成型设备的主要缺陷。(4)运用A N S Y S 对生物质成型徐变过程进行模拟分析,其变形过程及其压力与变形的关系与已有实验数据对比,其数据可靠,为生物成型机理的进一步研究提供了较好的依据和方法。参考文献1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 ANSYS 秸秆 生物 成型 仿真 分析
限制150内