中密度纤维板切削加工机理的研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
中密度纤维板(Medium Density Fiberboard,缩写为MDF)已经在人造板生产和应用中,占据了令人瞩目的地位。但MDF的二次加工,如切削加工,存在加工表面质量差、粉尘污染严重等诸多问题。本文作为国家自然科学基金资助项目“木材及木质复合材料高速切削过程中切屑流控制机理的研究”的重要部分,主要侧重于MDF切削过程中的切屑形成机理、切削力、铣削粉尘特性及切屑流特性进行分析研究。
     在MDF切屑形成机理研究方面,采用高速摄像技术和材料微观分析技术,进行了MDF正交切削试验,分析了切屑形成过程的基本特征,建立了MDF切削机理模型;进而分析了切削厚度、刀具前角及切削速度等不同切削参数对切削变形的影响,分析了切屑形态演变过程;建立了切削参数与切削加工表面破碎率的关系。研究结果表明:在MDF切削过程中,随着切削厚度的增大,MDF切屑由以粉尘状切屑为主首先逐步向以颗粒状切屑为主转变,进而向以准连续状切屑为主转变;在较大的切削厚度条件下,降低刀具前角或增大切削速度有助于MDF切屑由以准连续状切屑为主向以颗粒状切屑为主转变;在较小的切削厚度条件下,增大刀具前角或减小切削速度有助于MDF切屑由以粉尘状切屑为主向以颗粒状切屑为主转变。采用图像处理技术计算表面破碎率,不仅可以作为MDF已加工表面质量的评价指标,而且可以作为系统研究MDF切削机理、切削力等的重要手段。
     在MDF切削力研究方面,分析了不同切屑形成过程中切向力时域波形特征,分析了切削参数对MDF切向力和径向力的影响趋势,并且建立MDF切向力和径向力经验计算公式。研究结果表明:MDF准连续状切屑生成过程中,切削力‐时间曲线上的切向力呈一定周期性波动,切向力变化较大。颗粒状切屑形成时,切削力‐切削时间曲线上开始出现比较有规律性的锯齿状波形,这些高高低低的锯齿波峰值点对应着大大小小的MDF挤压颗粒状切屑的形成过程。粉尘状切屑形成时,平均切削力较小,围绕在42N左右变化,切削力‐时间曲线波动平稳。切削厚度对切削力的影响显著,随着切削厚度增大,切向力急剧增大,而径向力逐渐减小;刀具前角对切削力的影响较显著,随着刀具前角的增大,切向力减小较快,而径向力减小趋势较慢;切削速度对切削力的影响较显著,随着切削速度增大,切向力是逐渐减小,而径向力趋势表现并不明显。
     在MDF铣削加工产生的粉尘特性研究方面,分析了平均铣削厚度对MDF铣削粉尘粒度分布、典型粒度形貌及MDF已加工表面粗糙度的影响,并且分析了粉尘粒度对甲醛释放量的影响。研究结果表明:MDF铣削粉尘粒度分布主要集中在0.5-450μm区间,铣削MDF粉尘颗粒基本上为狭长型的纤维束、胶束或者短小型纤维束。随着平均铣削厚度增加,MDF铣削粉尘的平均粒度{D(4,3),D(3,2)}及中位径D(50)呈逐渐增大的趋势,而铣削速度影响较进给速度影响显著。随着平均铣削厚度的增加,MDF铣削表面质量降低,而铣削速度影响较进给速度影响显著。因此,可以认为,适当地增加铣刀转速有助于提高MDF铣削表面质量。此外,相同质量的MDF铣削粉尘,随着平均粒径减小,其甲醛释放量趋于增大。
     在MDF切屑流研究方面,基于高速摄像技术和图像识别与处理技术,采集木质复合材料铣削过程中切屑流的瞬态图像,分析MDF铣削切屑流形成机理、流场速度及扩散角度的变化规律。研究结果表明:切屑流形成过程分成两个阶段,第一个阶段为切削层材料离开工件形成切屑开始,到与齿槽接触;第二个阶段为切屑离开齿槽之后向外飞射,最后扩散在大气中。切屑流流场速度变化也分为两个阶段,第一个阶段为切屑流从形成时,具备一定流场速度后,开始趋于减小,第二个阶段为切屑流流场速度突然达到峰值,然后速度随之趋于减小。木质复合材料高速铣削时,随着刀齿切入工件,铣刀转速会趋于降低;而随着刀齿离开工件瞬间,刀齿线速度迅速恢复到额定速度并保持稳定。在切屑流流场内部的不同位置,切屑的速度存在差异。在同一时刻,切屑流中部的切屑颗粒速度最大,靠近刀齿部分速度次之,远离刀齿部分速度最小。随着铣刀转速的增加,切屑流的扩散角减小,在不同的木质复合材料铣削加工中,MDF铣削过程中的扩散角最大,变化范围为59.1‐66.5°。
     本文研究成果具有一定理论及应用价值,填补了国内外关于MDF切削机理研究的空白,为MDF切削过程中产生粉尘污染的控制,提供了可靠的理论依据。
The production and application of Medium Density Fiberboard (MDF) were developing sofast that it has already taken a significant position in the man-made board industry. But thefurther processing of the MDF boards, cutting for example, causes a series of problems such asbad processing quality, serious dust pollution and so on. This paper was an important part of theproject named “Study on the dust flow control of the wood and wood based materials in a highspeed cutting condition”, which was subsidized by Natural Science Foundation of China(NSFC). And this paper placed an emphasis on the study of chip forming mechanism, cuttingforce, dust characteristics of milling and the dust flow control.
     The MDF chip forming mechanism was investigated by testing the MDF orthogonalcutting with the technology of high speed camera and micro-analysis of the materials. The basiccharacteristics of the chip formation process and the influence of the cutting parameters such ascutting thickness, rake angle of the tool and cutting speed to the chip formation process havebeen studied. And the relationship between cutting parameters and the ratio of the surfacedamage had been analyzed. The result showed that: with the increasing of the cutting thickness,the dominated style of MDF chips changed from dust chip to granular chip and thenquasi-continuous chip; by increasing the rake angle or the cutting speed under the condition oflarge cutting thickness, the main chip style of the MDF changed from quasi-continuous chip togranular chip. While the cutting thickness was small, the leading chip type of MDF changedfrom granular chip to dust chip with the decreasing of the rake angle or the increasing of thecutting speed. The calculated ratio of surface damage based on the image processing technologycan not only be used as the evaluation index of the quality of processed surface, but also used asthe main approach, which was reasonably used in this paper, to study the cutting mechanismand cutting force of MDF systematicly.
     The characters of the time domain wave of cutting force during the chip forming processand the influence of the cutting parameters on the variation of MDF tangential force and radialforce were studied in the aspect of MDF cutting force research, which were leaded to anestablishment of experience formula for calculation tangential force and radial force of MDFcutting. It was showed that: approximate periodic fluctuations of Cutting Force-Time curve ofthe quasi-continuous chip forming process are found, and the variation of the tangential forceand radial force were obvious. During the formation of the granular chip formation, the Cutting Force-Time curve showed a regular serrated waveform, which corresponding to the differentsizes of MDF granular chip caused by the extrusion. When the cutting thickness was small, theaverage tangential force was comparatively small, more or less than42N, and the wave of theCutting Force-Time curve is comparatively smooth. The cutting thickness had obviousinfluence on the cutting force, and with the increasing of the cutting thickness the tangentialforce increased sharply, while the radial force decreased slowly. The rake angle of the tool alsohad some influence on the cutting force, with the increasing of the rake angle the tangentialforce decreased quickly while the radial force decreased slowly. The cutting speed also hadsome influence on the cutting force, with the increasing of the cutting speed the tangential forcedecreased gradually while the trend of variation of radial force had no significant impact.
     By not only testing the particle size distribution of the dust, typical particle sizemicrograph and the processed surface roughness with different average milling thickness butanalyzing the influence of different particle sizes on the formaldehyde emission, thecharacteristics of the dust produced during the MDF milling had been researched. The studyresult showed that: The distribution of the particle size was mainly ranging from0.5to450mand the particles of the dust chip can be classified into long-narrow fiber bundle, micelle orsmall fiber bundle. And with the increasing of the average milling thickness, the averageparticle size {D (4,3), D (3,2)} and the median-particle-size D (50) all tended to increase. Andthe cutting speed had more significant influence than the feed rate. However, with theincreasing of the cutting thickness, the processed surface quality decreased, in this case, thecutting speed also had more significant influence than the feed rate. Based on this result, properincreasing of the cutting speed can be considered as a helpful method to improve the processedsurface quality.
     Based on the high speed video technology and the image identifying and processingtechnology, the transient images of the chip flow during the wood-based composites’ millingprocess were recorded. The formation mechanism and the variation law of the velocity of flowfield and the diffusion angle were studied as a main part of the research on MDF chip flow. Thestudy result showed that: The forming process of the chip flow can be classified into two stages;the first stage was ranging from the uncut chip leaving away from the work piece to the chipcontacting with the gullet, and the second stage was that the chips left away from the gullet anddiffused into the air in the end. The velocity of the flow field can also be classified into twostages; the first one was the velocity of the flow field, got at the beginning of the chip flowformation, started to decrease, and the second one was the velocity of flow field increasedsuddenly to a peak and then decreased. While milling the wood-based composites with high speed, when the cutter tooth started to cut into the work piece the milling speed tended todecrease and at the moment of the cutter tooth leaving the work piece the linear velocity of thecutter tooth recovered quickly to the rated speed and then kept steady. In the same time, thevelocity of the chip in the middle of the chip flow was the fastest, the second location was theplace near the cutter tooth and the location away from the cutter tooth had the lowest velocity.With the increasing of the milling speed, the diffusion angle of the chip flow decreased. Amongthe wood-based composites, the MDF had the largest diffusion angle during the milling process,ranging from53.3to66.5°.
     Some significant theoretical and application value were showed in this paper which notonly filled in a gap on the study of the mechanism of MDF cutting at home and abroad, but alsoa theoretical basis for the dust pollution control during the MDF cutting process was provided.
引文
[1]科尔曼,库恩齐,施塔姆.木材学与木材工艺学原理—人造板[M].北京:中国林业出版社,1984.
    [2]李坚.走向21世纪的木质复合材料[J].世界林业研究,1995,3:34-40.
    [3]朱元鼎.我国人造板工业现状与发展趋势[J].林产工业,1998,(4):1-3.
    [4]刘邦庆.我国人造板工业发展中的问题与建议[J].林产工业,1998,(1):1-4.
    [5]张绍明.木材加工工艺[M].北京:高等教育出版社,2001.
    [6]徐咏兰.中密度纤维板制造[M].北京:中国林业出版社,1995.
    [7]许秀雯.纤维板生产工艺与技术[M].哈尔滨:东北林业大学出版社,1988.
    [8]王凯,袁东岩.积极利用农业剩余物大力发展人造板生产[J].木材工业,1990,(2):35-40.
    [9] Clarke P A. Panel products-past, present and future development[J].J Inst Wood Sci,1991,12:233-241.
    [10] Irle M, Loxton C. Manufacture and use of panel products in the UK[J].J Inst Wood sci,1996,14:21-26.
    [11] Chapman K M. Composite panels in Australia and New Zealand: a focus on MDF.In: Proceeding of the32ndInternational Particleboard/Composite Materials Symposium, Pullman, WA, USA,1998,75-82.
    [12] Sun X, Hammett A L, West C D. Hardwood use in China’s furniture industry [J]. For Prod J,1999,49:51-59.
    [13] Sun X, Hammett A L. Chinese furniture industry: its development and wood use [J]. For Prod J,1999,49:31-35.
    [14] Fuller B. Opportunities versus expectations: world market prospects for OSB and MDF.In: Proceedingsof the35thInternational Particleboard/Composite Materials Symposium, Pullman, WA, USA,2001:75-83.
    [15] Wu Q, Vlosky R P. Panel products: a perspective from furniture and cabinet manufacturers in thesouthern United States [J]. For Prod J,2001,50:45-58.
    [16] Olah D, Smith R, Hansen B. Wood material use in the U.S. cabinet industry1999to2001[J]. For Prod J,2003,53:25-31.
    [17] Duncanson J. Made in China-China becomes#1in MDF production[J]. Pap Tree Lett,2004,4:1-2.
    [18]南京林业大学.木材切削原理与刀具[M].北京:中国林业出版社,1983.
    [19]肖正福.木材切削刀具学[M].哈尔滨:东北林业大学出版社,1992.
    [20]李黎.木材切削原理与刀具[M].北京:中国林业出版社,2005.
    [21] Penman D, Olsson O J. Bowman C C. Automatic inspection of reconstituted wood panels for surfacedefects [J]. Proc of SPIE,1993,1823:284-292.
    [22] Salim H, Pantadee K,Kosonkorn, et al. Evaluation of surface roughness of Thai medium densityfiberboard(MDF)[J]. Inst J Building Environ,2006,41:527–533.
    [23] Chung K Y K, Cuthbert R J, Revell G S,et al. A study on dust emission, particle size distribution andformaldehyde concentration during machining of medium density fibreboard[J]. Ann Occupation HygieneSoc,2000,44:455–466.
    [24] McKenzie W M. Fundamental Analysis of the Wood Cutting Process[D].Ann Arbor: University ofMichigan,1961.
    [25] McKenzie W M. Fundamental Aspects of the Wood Cutting Process[J].Forest Products Journal,1961,10(9):447~456.
    [26] McKenzie W M. The Relationship Between the Cutting Properties of Wood and its PhysicalandMechanical Properties[J].Forest Products Journal,1962,12(6):287-294.
    [27] Franz N C. An Analysis of the Wood Cutting Process [D]. Ann Arbor: University of Michigan,1958.
    [28] Kivimaa, Eero. Cutting Force in Woodworking[J].The State Inst For Tech Res,1950,18:25-27.
    [29] Kivimaa Eero. Investigating Rotary Veneer Cutting with the Aid of a Tension Test[J]. Forest ProductJournal,1956,6(7):251-255.
    [30] Leney Lawrence. Mechanism of Veneer Formation at the Cellular Level[D]. Missouri:University ofMissouri,1960.
    [31] Leney Lawrence. A Photographic Study of Veneer Formation[J]. Forest Product Journal,1960,10(3):133-139.
    [32] McMuillin Charles W. The Relation of Mechanical properties of Wood and Nosebar Pressure in theProduction of Veneer[J].Forest Product Journal,1958,8(1):23-32.
    [33] Koch P. Utilization of Hardwoods Growing on Southern Pine Sites[M]. Washington D C: US Dept ofAgriculture, Forest Service,1985.
    [34] Koch P. Wood Machining Process[M].New York: The Ronald Press Company,1964.
    [35] Stewart H A. Analysis of Orthogonal Woodcutting Across the Grain[J].Wood Science,1971,12(1):38-45.
    [36]黄彦三等. Basic Analysis of Mechanism in Wood-Cutting Analysis of Cutting Energy in OrthogonalCutting Parallel to Grain [J].木材学会志,1974,20(2):77-82.
    [37] Pascal Triboulot. An Application of Fracture Mechanics to the wood-Cutting Process [J].木材学会志,1983,29(2):111-117.
    [38]大矢智等. Analysis of Wood Slicing Mechanism along the Grain I. Previous check occurrenceconditions [J].木材学会志,1992,38(12):1098-1104.
    [39]大矢智等. Analysis of Wood Slicing Mechanism along the Grain II. The transitive conditions offracture type from advanced check to bending rupture [J].木材学会志,1993,39(5):529-534.
    [40] Barney E Klamecki. Discontinuous Chip Formation in Woodcutting-A Catastrophe TheoryDescription[J]. Wood Science,1979,12(1):32-37.
    [41]王毅.竹节部切屑形成机理的研究[D].南京,南京林业大学,2005.
    [42]杨永福.竹材切削加工性能研究[D].北京,北京林业大学,2005.
    [43]王毅,曹平祥.竹材切削变形的初步研究[J].木材加工机械,2005,03:7-12.
    [44] Wen Ching Su,Yiren Wang, Nanfen Zhu, et al. Effect of tool angles on the chips generated duringmilling of wood by straight router-bits[J].Wood Sci,2003,49:271-274.
    [45] Su W C, Wang Y. Effect of the helix angle of router-bits on chip formation and power consumption inmilling of solid wood[J]. Wood Sci,2002,48:126-131.
    [46] Zhao C, Ando K, Hatori N, et al. Effects of conditions on concentration of respirable wood dust underbelt sanding(in Japanese)[J].Mokuzai Gakaishi,1999,45:306-312.
    [47]曹平祥,郝宁仲,王瑾.中密度纤维板切削力的研究[J].南京林业大学学报,1996,20(2):75-79.
    [48] Dippon J, Ren H, Engin F, et al. Cutting mechanics of meduium density fiberboard machining [J]. TechRept. The National Res. Council Canada Innovation Centre, Vancouver, B.C., Canada,1999.
    [49] Dippon J, Ren H, Amara F B, et al. Orthogonal cutting mechanics of medium densityfiberboards[J].Forest Products Journal,2000,50(7/8):25–30.
    [50] Engin S, Altintas Y, Amara F B. Mechanics of routing medium density fiberboard[J]. For Product Journal,2000,50:65–69.
    [51] Lin J T, Bhattacharyya D, Ferguson W G. Chip formation in the machining of SiC-particle-reinforcedaluminum-matrix composites[J].Composites Sci Technol,1998,58:258-291.
    [52] Lin J T, Van Houts J, Bhattacharyya D. Machinability investigation of medium density fiberboard [J].Holzforschung,2006,60:71–77.
    [53] Aguilera A, Meausoone P J, Martin P. Wood material influence in routing operations: the MDF case[J].Holz Roh.undWerkstoff,2000,58:278–283.
    [54] Davim J P, Clemente V C, Silva S. Surface roughness aspects in milling MDF(medium densityfiberboard)[J].The International Journal of Advanced Manufacturing Technology,2009,40(1-2):49–55.
    [55] Salje E, Druckhammer J. Qualit¨atskontrolle bei der Kantenbearbeitung[J]. Holz Roh Werkst,1984,42:187–192.
    [56] Feld H, Steward J, Horn T, et al. Carbide selection for woodworking tooling[J]. Proc of IWMS,2005,17:26–28.
    [57] Grzegorz Kowaluk, Waldemar Szymanski, Bartosz Palubicki, et al Examination of tools of differentmaterials edge geometry for MDF milling[J]. European Journal of wood and wood products.2009,67(2)173-176.
    [58] Grzegorz Kowaluk, Waldemar Szymanski,Adam Gilewicz, et al. Surface roughness aspects in machinecutting of medium density fibreboards with modified cutters on a CNC woodworking machine[J]. Forestryand Wood Technology.2011,75:202-2-9.
    [59] Davim J P, Reis P, Conceicao A C. Experimental study of drilling glass fiber reinforced plasticsmanufactured by hand lay-up[J]. Journal Composite Sci Technol.2004,64:287-297.
    [60] Davim J P, Clemente V, Silva S. Evaluation of delamination in drilling medium densityfibreboard.Proceedings of the Institution of Mechanical Engineers,Part B[J].Journal of EngineeringManufacture,2007,221(4):655-658.
    [61] Davim J P, Campos R J, Abrao A M. Delamination assessment after drilling medium density fiberboardusing digital image analysis. Holzforschung.2007,61(3):294-300.
    [62] Davim J P, Gaitonde V N, Karnik S R. An investigative study on drilling of medium-density fiberboard(MDF) using response surface models[J].The International Journal of Advanced Manufacturing Technology,2008,37(1/2):49-57.
    [63] Davim J P, Clemente V C, Silva S. Drilling investigation of MDF[J]. Int J Mater Processing Technol.2008,203(1/3):537-541.
    [64] Godfrey C, Onwubolu, Shivendra Kumar. Response surface methodology-based approach to CNCdrilling operations[J].Journal of Materials Processing Technology,2006,171(1):41-47.
    [65] Gaitonde V N, Karnik S R, Davim J P. Prediction and minimization of delamination in drilling ofmedium density fiberboard using response surface methodology and taguchi design[J]. Material ManufactureProcess.2008,23(4):377-384.
    [66] Gaitonde V N, Karnik S R, Davim J P. Taguchi multiple-performance characteristics optimization indrilling of medium density fiberboard to minimize delamination using utility concept[J]. Journal MaterProcess Technol.2008,196:73-78.
    [67] Prakash S, Palanikumar K, Manoharan N. Optimization of delamination factor in drilling mediumdensity fiberboards (MDF) using desirability based approach[J].The International Journal of AdvancedManufacturing Technology,2009,45:370–381.
    [68] Fischer R, Gottlober C, Rehm K, et al. A milling cutter as a screw-cutting instead ofhacking[C].In:Proceeding of the17thInternational wood Machining Seminar Rosenheim, Rosenheim,Germany,2005.
    [69] Gottlober C, Oertel M, Petrak A, et al. Development of milling tools with extreme inclination angle forplanning of wood[C].In:Proceeding of the19thInternational wood Machining Seminar,Nanjing,2009.
    [70] Leitz.Leitz Handbook for Woodworking Machine Tools[M].Nanjing.
    [71] Leuco.Leuco Handbook for Machining Wood,Wooden Materials and Plastics[M].Taicang.
    [72]陈日曜.金属切削原理[M].北京:机械工业出版社,2002.
    [73]陆剑中,孙家宁.金属切削原理与刀具[M].北京:机械工业出版社,2005.
    [74] Trent E M, Wright Paul. Metal Cutting[M].Woburn:Butterworth Heinemann Press,2000.
    [75]艾兴.高速切削加工技术[M].北京:国防工业出版社,2003.
    [76]杨合,詹梅.材料加工过程实验建模方法[M].西安:西北工业大学出版社,2008.
    [77]梅长林,范金城.数据分析方法[M].北京:高等教育出版社,2006.
    [78]陈杰. MATLAB宝典[M].北京:电子工业出版社,2010.
    [79]任彩亮,郭江伟.机械加工中切削力预报分析[J].机械研究与应用,2005,18(6):5-6.
    [80]龙震海,王西彬,刘志兵.高速铣削高强度钢时切削力影响因素析因与偏回归研究[J].航空制造技术,2005,9:77-80.
    [81]赵立.木材工业环境保护[M].北京:中国林业出版社,1985.
    [82]那斌.木地板厂的粉尘控制研究[J].江苏林业科技.2002,29(2):24-25.
    [83] Schlunssen V, Vinzents P S, Mikkelsen A B, et al. Wood dust exposure in the Danish furniture industryusing conventional and passive monitors[J].Annals of Occupational Hygiene,2001,45(2):157-164.
    [84] Harper M, Muller B S, Bartolucci A. Determining particle size distributions in the inhalable size rangefor wood dust collected by air samplers[J].Journal of Environmental Monitoring,2002,4(5):642-647.
    [85] Nylander L A, Dement J M. Carcinogenic effects of wood dust: Review and discussion[J]. AmericanJournal of Industrial Medicine,1993,24(5):619-647.
    [86]周捍东.我国木材工业粉尘污染的控制现状与进展[J].林产工业,2002,29(6):10-12.
    [87]周捍东.我国木材工业粉尘污染的控制现状与进展(续)[J].林产工业.2003,30(1):19-21.
    [88] Palmqvist J, Gustafsson S I.Emission of dust in planing and milling of wood [J]. Holz als Roh-undWerkstoff.1999,57(3):164-170.
    [89]赵建光,金维洙.木材加工企业的粉尘污染及防治措施[J].林业机械与木工设备,2008,36(12):52-54.
    [90] Koch P. Wood Machining Processes [M].New York: The Ronald Press Company,1964.
    [91] GB/T17657-1999.人造板及饰面人造板理化性能试验方法[S].
    [92]郭晓磊,曹平祥,江应.平均铣削厚度对MDF铣削粉尘粒度分布及表面粗糙度的影响[J].安徽农业科学,2010,38(19):10078-10080.
    [93]金维洙.减少木材加工企业粉尘污染的切削技术[J].中国安全科学学报,2006,16(10):100-104.
    [94]曹平祥.新型木工套装铣刀[J].林产工业,2009,36(5):66-69.
    [95] Palmqvist J, Gustafsson S I. Emission of dust in planning and milling of wood [J]. Holz als Roh-undWerkstoff,1999,57(3):164-170.
    [96] Rautio S, Hynynen P, Welling I, et al. Modelling of airborne dust emissions in CNC MDF milling[J].Holz Roh Werkst,2007,65(5):335–341.
    [97]曹平祥,周波,郭晓磊,等. MDF铣削过程中切屑流边界及扩散角的视频识别[J].林业科学,2010,46(1):112-116.
    [98]张东东,许宏庆,何枫,等.气固两相射流瞬时速度场和浓度场PIV研究[J].清华大学学报(自然科学版),2003,3(11):734-737.
    [99]李维礼.木材工业气力输送及场内运输机械[M].北京:中国林业出版社,1990.
    [100]董志勇.射流力学[M].北京:科学出版社,2005.
    [101]罗惕乾.流体力学[M].北京:机械工业出版社,2007.
    [102]车得福,李会雄.多相流及其应用[M].西安:西安交通大学出版社.2007,344-349.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700