第45卷第5期 西 安交通 大 学 学报 Vo1.45 No.5 Mav 2011 2011年5月 JOURNAL OF XI AN JIAOTONG UNIVERSITY 平面结合面切向接触阻尼分形模型及其仿真 张学良,温淑花,兰国生,丁红钦,张宗阳,王晓伟,刘志恒 (太原科技大学机械工程学院,030024,太原) 摘要:基于接触分形理论和结合面接触阻尼耗能机理,以及球体与平面接触时的阻尼耗能理论,建 立了平面结合面切向接触阻尼的分形模型,通过数值仿真直观揭示了平面结合面切向接触阻尼耗 能与结合面法向载荷以及结合面表面粗糙轮廓分形维数之间的非线性关系.仿真结果表明:平面结 合面切向接触阻尼的耗能随着结合面法向载荷的增大而减小;当结合面表面粗糙轮廓分形维数小 于等于1.2时,平面结合面切向接触阻尼耗能随结合面表面粗糙轮廓分形维数的增大而增大;当结 合面表面粗糙轮廓分形维数大于1.2时,平面结合面切向接触阻尼耗能随结合面表面粗糙轮廓分 形维数的增大而减小.仿真结果验证了模型的有效性. 关键词:平面结合面;接触阻尼;分形模型 中图分类号:TH1]3.1 文献标志码:A文章编号:0253—987X(2011)05—0074—04 Fractal Model for Tangential Contact Damping of Plane Joint Interfaces with Simulation ZHANG Xueliang,WEN Shuhua,LAN Guosheng,DING Hongqin, ZHANG Zongyang,WANG Xiaowei,LIU Zhiheng (Mechanical Engineering College,Taiyuan University of Science and Technology,Taiyuan 030024,China) Abstract:Based on the contact fractal theory and mechanism of contact damping dissipating ener— gy of joint interfaces as well as damping dissipating energy theory between a sphere and a plane, a fractal model for tangential contact damping of joint interfaces is proposed.Numerical simula— tion reveals the nonlinear relationship between tangential contact damping dissipating energy and the normal load on joint interfaces.It is concluded that the tangential contact damping dissipating energy of joint interfaces decreases with the normal load on joint interfaces.The tangential con— tact damping dissipating energy of j oint interfaces increases with the fractal dimension D when D ≤1.2,but when D>1.2 it decreases with D.The effectiveness of the model is verified by the experiment. Keywords:plane joint interfaces;tangential contact damping;fractal model 结合面在机械结构中大量存在,由于存在着接 触刚度和接触阻尼,因此对机械和机械结构的静动 态特性产生了重要影响.研究表明,机械结构的阻尼 约有90 9/6是由其中的结合面所产生的,因此人们对 结合面接触刚度和接触阻尼做了许多理论与实验研 究工作__】],并提出了结合面接触刚度的分形模 了总结,认为在结合面间不存在油介质的情况下,法 向交变力使结合面之间耗能很小,而切向交变力则 耗损振动能量形成阻尼作用.文献E43在对结合面阻 尼机理综合分析研究的基础上,认为,结合面间的微 观滑移是固定结合面接触阻尼产生的主要原因.文 型l2 ].文献E7]对结合面阻尼机理的研究成果进行 献E8—10]对结合面切向动态特性的实验研究结果表 明,结合面切向动态特性具有迟滞非线性,并且这种 收稿日期:2010—07—28. 作者简介:张学良(1964一),男,教授,博士生导师. 基金项目:国家自然科学基金资助项目 (50775153). 第5期 张学良,等:平面结合面切向接触阻尼分形模型及其仿真 迟滞非线性是结合面阻尼产生的主要原因,且每个 (2c 一( /2)2/(D-1) ’ k—H/a一 = L,振动周期的阻尼耗能与振动频率无关.文献[10]认 为,结合面间的微观滑移是结合面切向接触阻尼产 生的原因,但由于结合面接触阻尼的复杂性,因此对 其所做的理论研究工作依然很少.本文基于接触分 形理论和结合面的阻尼耗能机理,以及球体与平面 一d /E ( 式中:以 为临界接触面积;n 为最大接触点的接触 面积; (以)为微接触点的(面积)分布函数 ;D为 结合面表面粗糙轮廓分形维数;G为结合面表面粗 糙轮廓分形特征长度;H表示较软材料的硬度; 、, 表示较软材料的屈服强度. 接触时的阻尼耗能理论,研究建立了结合面接触阻 尼分形模型,并通过数值仿真研究揭示了结合面的 有关参数对接触阻尼的影响规律. 将式(4)、式(6)和式(7)代人式(5),得 (2一v)丌I/2Ta(2一一D)0 。 1平面结合面接触阻尼分形模型 平面结合面实质上是由2个粗糙表面组成的, 可以将其简化为一个粗糙表面与一个理想的刚性平 面的接触问题.对于粗糙表面上的单个微凸体,可以 将其近似等效为球体,其等效曲率半径为R.单个微 凸体与平面接触时,一个简谐振动周期内的接触阻 尼耗能可以表示为l1 e一 一—36G'— , r/ ̄p’ r一(以/兀)一 以 兀 (1)nf “ 一。 一 一以:。一 ] l (8) ]+ 。 4 ̄1—3 3 2D(可) /2E DG (D-1)ap/Z[盘 一 k荷a DaD/2缸 1州,。≠1.5 (9) P ===兀 G n( )+ (1o)3胁 va 以,1、 , D 一1. 5 式中:G 、 和 分别为结合面两接触材料的当量剪 式中:E为结合面两接触材料的弹性模量.由于 和 均是口1的函数,因此由式(8)和式(9)、式(10) 就建立了Ea和P 之间的非线性隐函数关系,从而 构成了平面结合面切向接触阻尼的分形模型. 切弹性模量、泊松比和摩擦系数; 为作用于微凸体 上的法向预加载荷;t为作用于微凸体上的切向动 态载荷的幅值;r为微接触区域的半径;a为接触 面积. 本文基于如下3个基本假设E3-4,6J:①粗糙表面 的微观形貌各向同性;②粗糙表面上各微凸体之间 2平面结合面切向接触阻尼分形模型 的数值仿真 由式(8)、式(9)和式(10)可知,给定固定的a,, 的相互作用可以忽略不计;③各微凸体所受的力与 其接触面积的大小成正比.作用于a上的t可以表 示为 t=aT/a (2) 便可以计算得到P..和 .根据文献[2—4],k一1.0, —1.0,G一1.19×10~1TI,E一2.1×10儿Pa,G 一 8.0×10∞Pa, 一3.53×10 Pa,v一0.3, 一0.2, 式中:T为作用于整个结合面上的切向动态载荷的 T一2.0×10一N,而D一1.1~1.9.图1、图2给出 幅值; 为结合面上的真实接触面积. 作用于a上的法向面压力为 P—ap /a (2一v)7【 /。Taa。/2 了数值仿真计算结果,且图中所给数据均满足结合 面间不发生宏观相对滑移的条件. (3) ,^、 数值仿真计算结果表明:①在平面结合面切向 简谐动载荷幅值一定,且结合面问不发生宏观相对 滑移的情况下(即为固定平面结合面), 随着 . .的增大而减小,二者之间具有很强的非线性关系,说 明在作用于平面固定结合面上的切向动载荷的幅值 一36G a: 式中:P..为作用于整个结合面上的法向预加载 荷 ]. 结合结合面接触阻尼耗能机理,平面结合面的 切向接触阻尼耗能为 Ed一 (以)da (5) (6) 定的情况下,平面结合面法向载荷增大时,平面结 合面的切向接触阻尼耗能则减小,这与平面结合面法 向和切向接触刚度的情况相反;②Ed随着D的变化 情况较为复杂,当D一1.1~1.2时,EI随着D的增大 而增大,当D>I.2时, 随着D的增大而减小,这与 平面结合面法向和切向接触刚度的情况也是相反的. http://ww、vI jdxb.CYI ( )一百D 西安交通大学学报 第45卷 拿 Z 莘 2 (a)D一1.1 (a)D在1.1~1.2时, 随D的增大而增大 蛊 Z 是 0 岫 (b)D一1.5 4 6 8 10 pul ̄lq Ⅲ.zNz_0_【 5 4 3 2 1 0 (b)D在1.2~1.4时, 随D的增大而减小 ‘0f o.8} 翠z叫 f 0-2} 【 2 4 6 8 1O (c)D一1.6 0_0O pJpN (c)D在1.3~1.7时,Ed随D的增大而减小 图2结合面表面粗糙轮廓分形维数 对接触阻尼耗能的影响 阻尼损耗因子 随着结合面法向面压力的增大而减 小,这与上述平面结合面切向接触阻尼分形模型的 仿真结论是一致的.图3c说明平面结合面切向接触 Pu/PN 阻尼损耗因子随着结合面粗糙度R 的减小而减小, (d)D=1.8 图1结合面接触阻尼耗能与法向载荷之间的关系 而结合面粗糙度的减小说明D在增大,因此平面结 合面切向接触阻尼损耗因子随着D的增大而减小, 这与上述仿真结论相一致. 3平面结合面切向接触阻尼实验 为了验证平面结合面切向接触阻尼分形模型的 有效性,本文对平面结合面切向接触阻尼进行了实 4结论 (1)平面固定结合面切向接触阻尼与作用在结 合面上的法向载荷之间具有很强的非线性关系,且 验研究,实验装置与方法见文献[12],实验结果如图 3所示.图3a、b和C均表明:平面结合面切向接触 随着法向载荷的增大,平面固定结合面切向接触阻 http://wⅥ jdxb.cn 第5期 张学良,等:平面结合面切向接触阻尼分形模型及其仿真 77 参考文献: [1]张学良,温淑花.机械结合面静动态特性研究回顾与 展望[J].太原重型机械学院学报,2002,23(3):276— 281. ZHANG Xueliang,WEN Shuhua.Review and pros— pect of the research on the static and dynamic charac— teristics of machine joint surfaces[J].Journal of Taiyuan Heavy Machinery Institute,2002,23(3):276— 281. (a)钢/铸铁磨削 E2]张学良,黄玉美,温淑花.结合面接触刚度分形模型研 究[J].农业机械学报,2000,31(4):89—91. ZHANG Xueliang,HUANG Yumei,WEN Shuhua. Fractal model of contact stiffness of joint surfaces I-J]. Transactions of the Chinese Society of Agricultural Machinery,2000,31(4):89—91. [3]张学良,温淑花.基于接触分形理论的结合面切向接 触刚度分形模型[J].农业机械学报,2002,33(2): 91—97. ZHANG Xueliang,WEN Shuhua.A fractal model of tangential contact stiffness of joint surfaces based on (b)铸铁/铸铁磨削 the contact fractal theory[J].Transactions of the Chi— nese Society of Agricultural Machinery,2002,33(2): 91—97. [4] 张学良.机械结合面动态特性及应用[M].北京:中 国科学技术出版社,2002:9. E5]温淑花,张学良,文晓光,等.结合面法向接触刚度分 形模型与仿真r-J-1.农业机械学报,2009,40(11):197— 202. WEN Shuhua,ZHANG Xueliang,WU Meixian,et a1.Fractal model and simulation of normal contact (c)钢/钢磨、铣削 stiffness of joint interfaces and its simulation[J]. Transactions of the Chinese Society for Agricultural 图3结合面切向接触阻尼损耗因子的实验结果 Machinery,2009,40(11):197—202. 尼减小,因此对于干结合面,要想通过调节结合面的 法向载荷(面压力)来同时提高结合面的接触刚度和 接触阻尼是不能实现的. [6]温淑花,张学良,文晓光,等.结合面切向接触刚度分 形模型及其仿真[J].农业机械学报,2009,40(12): 223—227. (2)结合面表面粗糙轮廓分形维数的大小在一 定程度上反映了其粗糙程度,较小的分形维数对应 着较大的表面粗糙度.对于干结合面,当结合面较为 粗糙(如D一1.1~1.2)时,平面结合面切向接触阻 尼随分形维数的增大而增大;当结合面的表面粗糙 度较小(如D>1.2)时,结合面切向接触阻尼随分形 WEN Shuhua,ZHANG Xueliang,WEN Xiaoguang, et a1.Fractal model of tangential contact stiffness of joint interfaces and its simulation I-J].Transactions of the Chinese Society for Agricultural Machinery,2009, 40(12):223—227. [7]戴德沛.阻尼减振降噪技术I-M].西安:西安交通大 学出版社,1986:6. 维数的增大而减小.因此,减小结合面表面粗糙度不 利于改善平面结合面的接触阻尼,但却能提高平面 结合面的接触刚度. (3)实验结果验证了本文所提出的平面结合面 切向接触阻尼分形模型的有效性. [8]ROGERS P F,BOOTHROYD G.Damping at metallic interfaces subjected tO oscillating tangential loads[J]. ASME J Eng for Industry,1975(8):1087—1093. (下转第136页) http://www.jdxb.cn 136 西安交通大学学报 第45卷 表4功能模块与结构模块映射表 功能模块 X向支撑 X向防护 y向支撑 向防护 [J].Annals of the CIRP,1996,45(1):1-6. 1-3]LEEWB,LAY H,LIU Z Z,et a1.A fuzzy analytic hierarchy process approach in modular product design 结构模块 X床身 X床身护板 立柱加滑座 立柱护板 口].Expert Systems,2001,18(1):32—41. [4] GU P,S0SALE S Product modularization for life cy— cle engineering[J].Robotics and Practice,1999(15): 387—401. [5]刘小鹏,张卫国,钟毅芳.机床模块化设计中的模块创 建及应用[J].华中理工大学学报,2000,28(5):16—17. LIU Xiaopeng,ZHANG Weiguo,ZHONG Yifang. The creation of machine tool module in modularization Z向支撑 Z向防护 B向支撑 B向进给功能 w向支撑 W向进给功能 X向进给 y向进给 Z向进给 X向位置检测 y向床身 Z床身护板 工作台 主轴单元 [J].Journal of Huazhong University of Science and Technology,2000,28(5):16—17. [63冯超.基于虚拟样机的虚拟轴机床设计研究[D].南 京:南京理工大学,2002. 丝杠螺母 丝杠螺母 [7]齐尔麦,徐燕申,谢艳.模块化数控机床概念设计的研 究[J].组合机床与数控化加工技术,2003(1):15—18. QI Ermai,XU Yanshen,XIE Yah.Research on the 齿轮齿条 conceptual design of modular NC machine tool[J]. Modular Machine Tool&Automatic Manufacturing y向位置检测 Z向位置检测 B向位置检测 光栅尺 Technique,2003(1):15—18. w向位置检测 冷却 润滑 排屑 [8]宗鸣摘,刘旭东,李湘嫒.产品模块化设计中模块划分 冷却系统 润滑系统 排屑系统 刀库、自动换刀装置 的多角度、分级特性讨论VJ].CAD/CMA与制造业信 息化,2003(1):50-52. ZONG Mingdi,LIU Xudong,LI Xiangyuam A multi— perspective hierarchical division approach for a product 自动换刀 数控 in its modular design[J].Digital Manufacturing Indus— try,2003(1):5O一52. 数控操作系统 [9]姜慧,徐燕中,谢艳.机械产品模块划分方法的研究 [J].制造技术与机床,1999(3):7-9. JIANG Hui,XU Yanzhong,XIE Yah.Research on 结果表明本文所提出的方法是有效可行的. 参考文献: [1]童时中.模块化原理设计方法及应用[M].北京:中国 标准出版社,2000:10—11. modular dividing method of machinery products[J]. Manufacturing Technology&Machine Tool,1999(3): 7-9. [1O]王莲芬.层次分析法[M].北京:中国人民大学出版 社,1982:18—24. [2]ENIXON G,VON YXKULL A,AMSTROM A.Mod— ularity:the basis for product and factory reengineering (编辑管咏梅) (上接第77页) [9]DEKONINCK c Deformation properties of metallic contact surfaces of joints under the influence of dy— El1]YAMADA A,KAKUBARI T.Prediction of dynamic characteristics of beam containing junction[J].Jour— nal of the Japan Society of Precision Engineering, 1986,52(12):2051-2057. namic tangential loads[J].International Journal of Machine Design and Research,1972(12):193—199. [12]ZHANG Xueliang,WEN Shuhua.Experiment re— search on tangential dynamic characteristics of ma— [10]堤正臣.机床螺栓结合面的阻尼机理[J].机床译丛, 1980(2):26—29. TSUTSUMI M.The damping mechanism for machine chined plane joint interfaces[J].Advanced Materials Research,2011,145:584—589. tool bolted joints[J].Machine Tools Translations, 1980(2):26—29. (编辑http:∥www.jdxb.cn 管咏梅)