林业科学  2018, Vol. 54 Issue (12): 149-156   PDF    
DOI: 10.11707/j.1001-7488.20181217
0

文章信息

张绍群, 颜家雄, 曹雷, 王海涛.
Zhang Shaoqun, Yan Jiaxiong, Cao Lei, Wang Haitao.
基于ADAMS和ABAQUS热磨机研磨木片断裂过程数值模拟
Crack Propagation Simulation of Hot Mill Grinding with Wood Based on ADAMS and ABAQUS
林业科学, 2018, 54(12): 149-156.
Scientia Silvae Sinicae, 2018, 54(12): 149-156.
DOI: 10.11707/j.1001-7488.20181217

文章历史

收稿日期:2017-03-23
修回日期:2017-06-05

作者相关文章

张绍群
颜家雄
曹雷
王海涛

基于ADAMS和ABAQUS热磨机研磨木片断裂过程数值模拟
张绍群, 颜家雄, 曹雷, 王海涛     
东北林业大学机电工程学院 哈尔滨 150040
摘要:【目的】基于ADAMS和ABAQUS软件对热磨机研磨木片的断裂过程进行数值模拟,研究当动、静盘磨齿间的夹角α一定、木片到主轴中心的距离R变化时的场变量输出,从理论上分析磨片破碎区域面积对纤维分离能耗的影响,指导热磨机磨片的优化设计。【方法】结合断裂力学理论,基于数值模拟方法,取动、静盘磨齿间的夹角α=15°,运用ADAMS软件求解木片到主轴中心的距离R分别为370、390、410、430、450 mm 5种情况时木片所受的研磨压力Fm和剪切力Fr,并利用ABAQUS软件将所求得的力加载在木片上进一步模拟木片断裂过程。【结果】1)在热磨机动、静磨齿研磨压力Fm和剪切力Fr的作用下,木片发生断裂破坏,α=15°、R=450 mm时木片的Mises应力变化范围为9.846×10-2~7.493 MPa,第一主应变变化范围为1.235×10-5~2.546×10-2,第三主应变变化范围为-1.231×10-4~-1.853×10-2,应变E23变化范围为-3.481×10-2~3.756×10-3,较大的应力应变主要分布在裂纹附近,木片左右两侧的应力应变逐渐减小。α=15°、R=450 mm时木片节点位移U1的变化范围为-7.605×10-4~1.693×10-3mm,节点位移U3的变化范围为-3.334×10-2~1.966×10-2mm,U3的变化范围大于U1的变化范围,木片主要沿U3方向被撕裂,即木片在被研磨过程中主要在剪切力Fr的作用下沿木片横纹方向被分离。2)在磨片破碎区,当动、静盘磨齿间的夹角α一定时,木片到主轴中心的距离R增大,木片的Mises应力、正应力、Magnitude、节点位移呈增大趋势。当木片到主轴中心的距离R越大,即破碎区面积越大时,木片所受应力应变增大,木片断裂越容易。【结论】1)木片被研磨过程中主要在剪切力Fr的作用下沿木片横纹方向被分离,对磨片进行结构设计和优化分析时应尽可能增加磨齿对木片的剪切力Fr,使磨片破碎区木片的断裂程度更充分,降低粗磨区和精磨区的纤维分离能耗。2)热磨机研磨木片时木片到主轴中心的距离R越大,即破碎区面积越大,木片断裂就越容易,合理设计磨片破碎区域面积可以降低纤维分离能耗。3)运用ABAQUS软件模拟木片到主轴中心的距离R不等时的场变量输出,从理论上分析磨片破碎区域面积对纤维分离能耗的影响,可以指导热磨机磨片的优化设计。
关键词:热磨机    木片断裂    纤维分离    数值模拟    
Crack Propagation Simulation of Hot Mill Grinding with Wood Based on ADAMS and ABAQUS
Zhang Shaoqun, Yan Jiaxiong, Cao Lei, Wang Haitao     
School of Mechanical and Electrical Engineering, Northeast Forestry University Harbin 150040
Abstract: 【Objective】In order to improve the design of the mill grinding disc, the fracture process of grinding the wood chips into fibers by the hot mill was studied in this paper. Based on ADAMS and ABAQUS, the process of grinding the wood chips by the hot mill was numerically simulated, and the process of gradual fracture of the wood chips with the force was explored. When the angle between the grinding and grinding teeth of the moving and stationary discs is constant, the amount of field change was recorded for the distance R from the wood piece to the center of the spindle. The influence of broken area size of the grinding piece on the energy consumption of fiber separation was theoretically analyzed.【Method】Combined with the theory of fracture mechanics, the method of numerical simulation is used. When the angle between the teeth of the moving and stationary discs α is 15°, the numerical simulation method is applied to solve the distance R between the wood chips and the spindle center by using ADAMS software to obtain R=370, 390, 410, 430, 450 mm in five cases, the size of the grinding pressure Fm and the shear force Fr. Then use ABAQUS software to load it on the wood chip and numerically simulate the fracture process.【Result】The ADAMS software was used to simulate the working process of grinding the wood chips by the hot mill, and the simulation result of the grinding pressure Fm and the shearing force Fr of the wood chips α=15° and R=450 mm were obtained. The field variable output result of α=15° and R=450 mm were obtained by using ABAQUS software. The variation trend of Mises stress, normal stress, Magnitude and node displacement of wood chips under different conditions was studied by using Excel scatter plot and regression analysis chart.【Conclusion】1) The damage form of the wood chips during the grinding process is mainly separated by the shear force Fr along the transverse direction of the wood chips. Therefore, during the structural design and optimization analysis of the grinding plates, the grinding teeth should be increased as much as possible. The force of the cutting force Fr makes the fracture degree of the wood chip in the crushing zone more complete, and reduces the energy consumption of fiber separation in the rough grinding zone and the fine grinding zone. 2) When the hot mill grinds the wood chips, the larger the distance R from the wood chips to the center of the main shaft, the larger the area of the crushing area, the easier the fracture of the wood chips is. The reasonable design of the area of the crushing area of the grinding plate can reduce the energy consumption of fiber separation. 3) Using ABAQUS to simulate the field variable output when the distance from the wood chip to the center of the spindle is not equal, the influence of the area of the crushed area on the fiber separation energy consumption is theoretically analyzed, which can guide the optimized design of the grinding piece.
Key words: hot grinding    wood fracture    fiber separation    the numerical simulation    

20世纪60年代中期,中密度纤维板问世,随后取得了迅速发展。纤维是中密度纤维板的质量基础,研磨纤维是中密度纤维板制造的重要工序之一。热磨机研磨纤维的过程为蒸汽软化后的木片送入研磨室磨盘,使木片受压缩、拉伸、剪切、扭转、冲击、摩擦和水解等多次重复的外力作用,最终导致纤维分离,其实质是木片原料的断裂和解离(陈光伟,2012徐大鹏,2012)。磨片是热磨机进行研磨分离的关键部件,直接参与研磨工作,磨片的结构参数设计是否合理,会直接影响纤维形态质量、热磨机能量耗散和磨片平均使用寿命(陈光伟等,2010)。多年来,国内外学者从不同角度研究热磨机磨片结构影响因素,并取得了不少研究成果(Berg, 2001; Chen et al., 2008; Resch et al., 2012; Fan et al., 1999; Hua et al., 2008; Li et al., 2007; Pelletier et al., 2013; Qin et al., 1999; Santos et al., 2009; Shida et al., 2010)。Olson等(2003)选用5种磨片采用建模法与试验法相结合对纤维分离强度和比能耗之间的关系进行研究,得出磨片齿刃比负荷与选择函数间具有显著相关性;Roux等(2009)运用建模法分析精磨机磨片磨齿对纤维的切断作用,得出纸浆研磨参数可以预测磨齿切断纤维的情况;Holmberg等(2009)通过试验对磨齿研磨木片和切断的过程进行研究,得出木材纤维的断裂特性,并分析了纤维分离中磨齿的作用;May(1973)研究磨片结构参数对木材研磨强度的影响,从理论角度分析了纤维长度的形成机制;沈立新(1999)对磨片齿形结构和选择理论依据进行研究,分析磨片的几何形状,并以扇形分区磨片和圆环分区直齿磨片的设计方法为例进行了详细阐述;李艳丽等(2008)设计扇块分区磨片的整体结构,通过仿真分析方法研究磨片在静载荷作用下的应力应变和热分布情况,并以降低磨片质量为优化目标对磨片结构进行了优化设计;邵贤林等(2013)研究磨片间隙变化量对制浆质量和磨片使用寿命的影响,并对常用磨片间隙调节方式的特点和应用进行了阐述。

总体来看,以往研究主要采用固体力学和流体力学等方法对磨片与木片原料之间的相互作用关系进行分析和建模,很难直接描述木片原料受力后逐步断裂的过程机制。鉴于此,本研究基于ADAMS和ABAQUS软件对热磨机研磨木片的断裂过程进行数值模拟,从理论上分析磨片破碎区域面积对纤维分离能耗的影响,以指导热磨机磨片的优化设计。

1 热磨机作用下木片受力分析

热磨机研磨木片时,木片断裂主要在研磨室内进行,磨片破碎区受力状况最为复杂,因此,本研究以热磨机研磨木片时的破碎区为例进行分析,如图 1所示。设热磨机工作时的输出功率为P(kW),转速为n(r·min-1),木片尺寸为L×a×h(mm)(长×宽×高),磨片内径为R0(mm),磨片外径为R1(mm),木片到主轴中心的距离为Ri(mm);动磨盘随主轴高速旋转时,动、静盘磨齿间的夹角为α(锐角),动盘磨齿对木片的剪切力为Fr(N),方向垂直于动盘磨齿齿面指向木片一侧;静磨盘磨齿对木片的反作用力为Fs(N),方向垂直于静盘磨齿齿面指向木片一侧;木片随着动磨盘高速旋转产生的离心力为Fc(N),方向沿径向指向磨盘外侧;动磨盘对木片的研磨压力为Fm(N),方向垂直于动磨盘面指向精磨盘;静磨盘对木片的反向支撑力为Fn(N),方向与研磨压力方向相反;静盘磨齿对木片运动的摩擦阻力为Ff(N),方向与木片运动方向相反,木片受到的蒸煮压力为S

图 1 热磨机作用下木片的受力分析 Fig. 1 Mechanical analysis of wood chips under the action of hot mill

分析木片受力情况可知,使木片断裂的力主要是动盘磨齿对木片的剪切力Fr沿x轴方向的分力Fr1和动磨盘对木片的研磨压力Fm,而离心力Fc、蒸煮压力S和摩擦阻力Ff等主要与纤维的运动状态相关,促使木片向磨片外部运动,影响纤维在热磨区的停留时间,对木片断裂的影响很小。

2 基于ADAMS热磨机研磨木片的动力学仿真分析 2.1 仿真工作环境分析

热磨机研磨木片,是在动、静磨齿对木片的研磨压力Fm和剪切力Fr共同作用下,木片沿纤维结合力最弱的方向断裂,受离心力Fc、蒸煮压力S的影响很小。实际生产中,由于木片尺寸大小不一、研磨室内木片剪切面随机性很大以及研磨室内木片运动轨迹很难预测,对热磨机研磨木片的仿真分析存在一定难度;随动磨盘的高速旋转,木片在动、静磨片的挤压和剪切作用下被逐渐破碎和切断,由于木材横纹强度远小于顺纹强度,软化后更弱,因此木片沿横纹方向被剪切的可能性最大。本研究对单片木片沿横纹方向被剪切的情况进行分析。同时为了简化仿真模型,作以下假设:1)研磨室内被研磨木片与动磨盘磨齿接触瞬间木片处于静止状态;2)待研磨木片尺寸均匀,其磨片破碎区木片长度方向与磨齿侧面平行,宽度方向沿齿面的法向方向。

2.2 研磨压力Fm和剪切力Fr求解

选取单片木片作为研究对象,热磨机研磨木片的状态简化为图 2所示。

图 2 热磨机研磨木片的仿真模型 Fig. 2 The simulation model of grinding wood chip in hot mill 1.静磨盘Static grinding disc; 2.木片Wood chip;3.动磨盘Dynamic grinding disc.

运用ADAMS软件模拟热磨机研磨本片的工作过程,首先在ADAMS软件中建立热磨机研磨木片的仿真模型,构建仿真工程环境,然后以热磨机实际工况对仿真模型添加约束条件、设置各运动副、定义动磨盘转速n=1 500 r·min-1,进行仿真分析。α=15°(动、静盘磨齿间的夹角)、R=450 mm (木片到主轴中心的距离)时热磨机研磨木片的研磨压力Fm和剪切力Fr的仿真结果如图 3所示,可以看出由于动磨盘高速旋转,磨齿对单片木片的作用时间很短,剪切力Fr和研磨压力Fm呈先增大后减小的变化趋势。

图 3 α=15°、R=450 mm时木片所受的研磨压力Fm和剪切力Fr仿真结果 Fig. 3 The simulation results of grinding pressure(Fm) and shear stress (Fr) of wood chip at α=15°, R=450 mm a.木片所受剪切力Fr沿x方向的分力The shear force (Fr) along the x direction of wood chip; b.木片所受剪切力Fr沿y方向的分力The shear force (Fr) along the y direction of wood chip; c.木片所受研磨压力The grinding pressure (Fm) of wood chip.
3 木片断裂裂纹扩展数值模拟 3.1 有限元模型

木片进入研磨室后排列方向是随机的,很难预测,为了方便,本仿真限定在以下工况下进行:1)仿真针对磨片破碎区木片长度方向与磨齿侧面相平行、宽度方向沿齿面法向方向的情况进行研究;2)以单片木片主要沿横纹方向被剪切的情况进行研究。

图 4所示,以磨片上磨齿的Fr沿木片横纹方向被剪切、Fm沿木片顺纹方向被挤压的情况进行仿真分析。

图 4 木片仿真模型的弹性主轴方向 Fig. 4 The direction of elastic principal axis of wood chip simulation model
3.2 材料属性定义及网格划分

采用ABAQUS软件建立热磨机研磨木片的有限元模型,用于仿真的木材为杨木,其基本参数为:ρ=0.386 g·cm3EL= 5 352 MPa,ER=1 204.7 MPa,ET=1 039 MPa,GRL= 621 MPa,GTL=466.3 MPa,GRT=94.3 MPa,μRT=0.452,μLR=0.64,μLT=0.66。进行网格划分,整个三维有限元模型均采用C3D8R,即八结点线性减缩积分六面体单元,结构化网格划分后得到57 702个单元83 753个节点。

3.3 边界条件和载荷定义

图 5所示,热磨机研磨木片时木片侧面与磨片底面接触,在ABAQUS软件中定义木片侧面边界条件为U2=0、U3=0、UR1=0;磨齿侧面与木片顶面的一部分接触,在ABAQUS软件中定义木片侧面边界条件为U2=0、U3=0、UR1=0、UR2=0。热磨机研磨木片时,木片断裂主要是由于动磨盘随主轴沿逆时针方向高速旋转对木片施加剪切力和挤压力作用,对木片添加载荷如图 6所示。

图 5 定义边界条件 Fig. 5 Defining boundary conditions
图 6 添加载荷 Fig. 6 Loading addition
3.4 仿真结果 3.4.1 裂纹扩展过程

仿真木片长度(25 mm)方向为径向(R),宽度(10 mm)方向为纵向(L),高度(5 mm)方向为弦向(T),热磨机研磨木片(10 mm×5 mm)时的断裂过程如图 7所示。木片被研磨过程中主要在剪切力Fr的作用下沿木片横纹方向被分离,随时间增加木片达到塑性破坏,裂纹附近区域的应力分布呈非对称现象,沿木片左右两侧依次减小。

图 7 裂纹扩展过程 Fig. 7 Crack propagation process
3.4.2 不同条件下场变量输出结果

为了研究木片在磨片破碎区内不同位置时对木片断裂的影响,运用ADAMS求解木片到主轴中心的距离R分别为370、390、410、430、450 mm 5种情况时木片所变的研磨压力Fm和剪切力Fr,并将求得的研磨压力Fm和剪切力Fr分别加载在木片上,运用ABAQUS对热磨机作用下木片的裂纹扩展进行数值模拟。图 8所示为α=15°、R=450 mm时的场变量输出结果,不同条件下热磨机研磨木片时木片的Mises应力、正应力、Magnitude和节点位移的变化趋势利用Excel散点图作线性回归分析,如图 9所示。

图 8 α=15°、R=450 mm时场变量输出结果 Fig. 8 Field variables output results at α=15°, R=450 mm a.Mises应力云图Mises stress cloud image;b.第一主应变云图First principal strain cloud image; c.第三主应变云图Third principal strain cloud image;d.应变E23云图Strain E23 cloud image;e.节点位移U1云图Node displacement U1 cloud image;f.节点位移U3云图Node displacement U3 cloud image.
图 9 木片到主轴中心的距离R不相等时场变量输出线性回归 Fig. 9 Linear regression of field variables output when the distance between wood chip and spindle center is not equal
4 结论

1) 在热磨机动、静磨齿研磨压力Fm和剪切力Fr作用下,木片发生断裂破坏,α=15°、R=450 mm时木片的Mises应力变化范围为9.846×10-2 ~7.493 MPa,第一主应变变化范围为1.235×10-5~ 2.546×10-2,第三主应变变化范围为-1.231×10-4~-1.853×10-2,应变E23变化范围为-3.481×10-2~3.756×10-3,较大的应力应变主要分布在裂纹附近,木片左右两侧的应力应变逐渐减小。α=15°、R=450 mm时木片节点位移U1的变化范围为-7.605×10-4~1.693×10-3mm,节点位移U3的变化范围为-3.334×10-2~1.966×10-2mm,U3的变化范围大于U1的变化范围,木片主要沿U3方向被撕裂,即木片在被研磨过程中主要在剪切力Fr的作用下沿木片横纹方向被分离。因此,在对磨片进行结构设计和优化分析时,应尽可能增加磨齿对木片的剪切力Fr,使磨片破碎区木片的断裂程度更充分,降低粗磨区和精磨区的纤维分离能耗。

2) 在磨片破碎区,当动、静盘磨齿间的夹角α一定时,木片到主轴中心的距离R增大,木片的Mises应力、正应力、Magnitude、节点位移呈增大趋势。当木片到主轴中心的距离R越大,即破碎区面积越大时,木片所受应力应变增大, 木片断裂越容易。因此,合理设计磨片破碎区域面积对纤维分离质量和能耗具有重要影响。

参考文献(References)
陈光伟, 花军, 纪伟, 等. 2010. 磨片结构对纤维分离过程中能量转换机理的影响. 东北林业大学学报, 38(8): 109-110, 114.
(Chen G W, Hua J, Ji W, et al. 2010. Effects of abrasive disc- structure on energy transformation during fiber separation. Journal of Northeast Forestry University, 38(8): 109-110, 114. DOI:10.3969/j.issn.1000-5382.2010.08.033 [in Chinese])
陈光伟. 2012.热磨法磨片纤维分离机理的模型分析与实验研究.哈尔滨: 东北林业大学博士学位论文, 22-30.
(Chen G W. 2012. Model analysis and experimental study on mechanism of thermomechanical refining plate separating fiber.Harbin: PhD thesis of Northeast Forestry University, 22-30. [in Chinese])
李艳丽, 高秀华, 陈淑清. 2008. 扇块分区热磨机磨片的热应力分析与优化设计. 林业机械与木工设备, 36(1): 29-31.
(Li Y L, Gao X H, Chen S Q. 2008. Thermal stress analysis and optimization design of fan block to zone grinding mill. Forestry Machinery & Woodworking Equipment, 36(1): 29-31. DOI:10.3969/j.issn.2095-2953.2008.01.008 [in Chinese])
邵贤林, 王平. 2013. 盘磨机磨盘间隙的调整方式与分析. 纸和造纸, 32(3): 21-25.
(Shao X L, Wang P. 2013. Adjustment and analysis of disc gap in disc refiner. Paper and Paper Making, 32(3): 21-25. [in Chinese])
沈立新. 1999. 浅谈盘磨磨片齿型的设计与选择. 上海造纸, (1): 27-30, 33.
(Shen L X. 1999. Discussion on the design and selection of disc grinding die. Shanghai Paper Making,, (1): 27-30, 33. [in Chinese])
徐大鹏. 2012. MDF纤维分离机理及磨片齿形结构参数优化设计的研究.哈尔滨: 东北林业大学硕士学位论文, 4-15.
(Xu D P. 2012. Study on the fiber separation mechanism of MDF and the optimization design of plate-bar structural parameters. Harbin: MS thesis of Northeast Forestry University, 4-15.[in Chinese])
Berg J E. 2001. Effect of impact velocity on the fracture of wood as related to the mechanical pulping process. Wood Science and Technology, 35(4): 343-351. DOI:10.1007/s002260100100
Chen G W, Hua J. 2008. Test analysis and deduction of factors affecting fiber separating quality and mechanism refiner-disc working. International Association of Wood Products Societies: 227-228.
Fan K, Hatzikifiakos S G, Avramidis S. 1999. Determination of the surface fractal dimension from sorption isotherms of five softwoods. Wood Science and Technology, 33(2): 139-149. DOI:10.1007/s002260050105
Hua J, Chen G W. 2008. Theoretic deduction and computation of material movement track in parting fiber. International Association of Wood Products Societies(IAWPS2008): 235-236.
Holmberg S, Persson K, Petersson H. 2009. Nonlinear mechanical behavior and analysis of wood and fibre materials. Computers & Structures, 72(4/5): 459-480.
Li J G, Pang S S, Scharpf E W. 2007. Modeling of thermal energy demand in MDF production. Forest Products Journal, 57(9): 97-104.
Olson J A, Drozdiak J, Martinez M, et al. 2013. Characterizing fiber shortening in low-consistency refining using a comminution model. Powder Technology, 129(1/3): 122-129.
May W D. 1973. A theory of chip refining the origin of fibre length. Pulp and Paper Research Institude of Canada, 24: 245-253.
Pelletier A L, Zhao Y, Lei X, et al. 2013. Improved fiber separation and energy reduction in thermomechanical pulp refining using enzyme-pretreated wood. BioResources, 8(3): 3385-3398.
Qin D C, Guan N, Jiang X M. 1999. Morphology of wood failure in relation to the variation in tensile strength parallel to grain of three hard pines. Journal of the Institute of Wood Science, 15(1): 1-5.
Resch E, Kaliske M. 2012. Numerical analysis and design of double-shear dowel-type connections of wood. Engineering Structures, 41: 234-241. DOI:10.1016/j.engstruct.2012.03.047
Roux J C, Mayade T L. 2009. Modeling of the particle breakage kinetics in the wet mills for the paper industry. Powder Technology, 105(1/3): 237-242.
Santos C L, De Jesus A M P, Morais J J L, et al. 2009. Quasi-static mechanical behaviour of a double-shear single dowel wood connection. Construction and Building Materials, 23(1): 171-182. DOI:10.1016/j.conbuildmat.2008.01.005
Shida S, Hiziroglu S. 2010. Evaluation of shear strength of Japanese wood species as a function of surface roughness. Forest Products Journal, 60(4): 400-404. DOI:10.13073/0015-7473-60.4.400