摩擦式离合器滑摩过程及温度场分析

来源 :计算机辅助工程 | 被引量 : 0次 | 上传用户:mario2
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  摘 要:为更真实地反映摩擦式离合器滑摩过程中的温度场,提出1种利用动力学软件与有限元软件相互辅助分析的方法.在试验测得的负载扭矩的基础上,建立摩擦片和内鼓轮及气胎的模型;利用MSC Adams对摩擦式离合器进行动力学仿真,通过对摩擦片加压计算出滑摩时间和主、从动端的相对转速;将该相对转速用于有限元热机耦合场分析,通过摩擦片和内鼓轮的相对转动形成滑摩、生成热量,得离合器各时刻的温度场.结果表明:离合器在接触后3.064 s时主、从动端接合完毕,且摩擦片瞬态最高温度只有65.83 ℃.该方法对摩擦式离合器滑摩过程及温度场分析具有一定的实际意义.关键词:摩擦式离合器; 滑摩时间; 热机耦合; 温度场中图分类号:TH133.4; TH117.1; TB115 文献标志码:
  A
  Analysis on frictional sliding process andtemperature field of friction clutch
  CHEN Huiliang
  1a, WU Xiugen
  1b, FU Kunkun2, SHEN Mingxing2
  (1. a. School of Aerospace Eng. & Mechanics; b. School of Civil Eng., Tongji Univ., Shanghai 200092, China;2. 704th Research Institute, China Shipbuilding Industry Co., Shanghai 200031, China)
  Abstract: To reflect the temperature field of friction clutch in the frictional sliding process more realistically, a method using multibody dynamics software and finite element software which are used to perform mutual aidedanalysis is proposed. Based on the load torque obtained from experiments, the models of friction plate, inner drum and pneumatic tyre are built; multibody dynamics simulation of friction clutch is performed by MSC Adams to get frictional sliding time and relative rotational speed between master side and slave side; the relative rotational speed is used for finite element analysis on thermalstructural coupling, the frictional sliding and heat are generated by relative motion between friction plate and inner drum, and the temperature filed is obtained at each time step. The results show that the master side and slave side complete bonding engagement at 3.064 s after clutch contacts, and the highest temperature is 65.83 ℃ only. The method is practical for the analysis on frictional sliding process and temperature field of friction clutch.Key words: friction clutch; frictional sliding time; structuralthermal coupling; temperature field
  收稿日期:2010[KG*9〗03[KG*9〗29 修回日期:2010[KG*9〗09[KG*9〗02作者简介: 陈惠亮(1986—),男,上海人,硕士研究生,研究方向为计算固体力学,(Email)chl041037@sina.com0 引 言
  摩擦式离合器是机械工程中重要的传动装置,其主要功能为实现主、从动轴在短时间内方便、高效、平稳地接合和分离.在摩擦过程中,摩擦表面间的相互滑摩会产生大量热量,使离合器温度升高,而温度过高会引起离合器性能改变,甚至导致事故发生,有必要分析离合器的滑摩过程及其温度场.
  离合器热负荷的计算主要采用经验公式,利用滑摩功能简单估算出摩擦物体的温度
  [1],但不够准确.目前,对于摩擦片的温度场分析主要采用有限元法,通过计算摩擦表面的热流密度
  [24],考虑热传导、热交换来计算温度场.但每个时刻摩擦表面都相对运动,未接触表面和接触表面的传热和换热边界条件不一致,因此温度分布也不相同.若不考虑转动,则不能体现实际动态接触产生不一致的温度场.这种缺点可通过设定摩擦表面按实际情况运转来克服:BABU等
  [5]提出考虑旋转的动态接触耦合分析温度场的方法,但只计算了二维情况,且其负载扭矩通过编程利用简单公式得出,滑摩转矩影响滑摩时间,对温度场影响很大.
  本文利用动力学软件与有限元软件相互辅助分析的方法,通过动力学分析考虑实际的负载扭矩情况,计算得到滑摩时间和主、从动端的相对转速.考虑离合器三维模型,利用有限元热机耦合分析计算离合器的瞬态温度场,并考虑接触传热系数,计算出稳态温度场.1 离合器接合过程分析
  摩擦式离合器接合过程分为4个阶段
  [6]:第1阶段,随着气胎的逐渐加压,气胎逐渐膨胀,消除摩擦表面的间隙;第2阶段,当间隙消除后,随着气胎压力的进一步增大,摩擦力逐渐增加,即摩擦转矩逐渐增大,但此阶段滑摩转矩不足以克服负载扭矩,从动端静止不动;第3阶段,随着气胎压力的进一步增大,摩擦力也进一步增加,当滑摩转矩大于负载转矩时,在滑摩转矩的作用下,从动端开始加速,直到主从动端角速度相等,达到稳定状态;第4阶段,从动端和主动端的角速度达到稳定,气胎压力继续增加到最大,此时滑摩转矩最大.2 动力学仿真
   图 1 摩擦片和内鼓轮及气胎的模型通过试验测得负载扭矩,建立摩擦片和内鼓轮及气胎的模型(见图1),应用MSC Adams进行动力学仿真,分析中转动惯量影响扭矩,本文通过改变密度设定主、从动端的转动惯量.将各摩擦片初始位置设定为即将与内鼓轮接触的位置,即摩擦片与内鼓轮贴合无间隙.通过压力控制其接触.从气胎逐渐加压消除间隙到刚与内鼓轮接触需要0.8 s,动力学仿真从刚接触瞬间开始分析,即从第2阶段开始分析.动摩擦因数为0.3,静摩擦因数为0.35.设以恒定的25 r/min转速模拟主动端电机的持续运转;设定摩擦片和气胎为径向移动副,并在摩擦片上设定沿径向力.设定的约束关系和运动关系见表1.表 1 设定的约束关系和运动关系气胎摩擦片内鼓轮大地气胎无径向滑动副,并有径向力作用无旋转副,并设定扭矩摩擦片径向滑动副,并有径向力作用无接触(考虑摩擦)无内鼓轮无接触(考虑摩擦)无旋转副,并设定恒定转速大地旋转副,并设定扭矩无旋转副,并设定恒定转速无负载扭矩对确定滑摩时间起重要作用,本文由试验测得滑摩扭矩,输入到MSC Adams中,负载扭矩时程曲线见图2.在开始阶段由于摩擦力矩小,负载扭矩随着摩擦力矩的增大而增大,此时内鼓轮速度不变;随着摩擦力矩的增加,负载扭矩也随之增加;最后达到稳定.与上文分析一致.分析起始时刻为摩擦片与内鼓轮刚接触时的时刻,分析时长为
  9.4 s.内鼓轮与气胎的角速度时程曲线见图3,在摩擦片与内鼓轮接触到3.064 s时主、从动端转速达到稳定、速度相等,即此时接合完成.图 2 负载扭矩时程曲线图 3 内鼓轮与气胎的角速度时程曲线3 温度场分析
  MSC Marc采用弱耦合方式进行热机耦合分析,即先进行热传导分析,然后进行应力分析.温度影响变形,而变形不影响温度.控制方程为[WTHX]M[WTBX]u¨+[WTHX]D[WTBX]u·+[WTHX]K[WTBX](T,t)[WTHX]u[WTBX]=[WTHX]F[WTBX][JY](1)[WTHX]c[WTBX](T)T·+KT(T)[WTHX]T[WTBX]=[WTHX]Q[WTBX]+[WTHX]Q[WTBX]F[JY](2)式中:[WTHX]M[WTBX]为质量矩阵;[WTHX]D[WTBX]为阻尼矩阵;[WTHX]K[WTBX](T,t)为与温度、时间相关的刚度矩阵;F为外力矢量;[WTHX]u[WTBX]为节点位移矢量;[WTHX]T[WTBX]为节点温度矢量;Q为热通量向量;QF为摩擦热量;c[WTBX]为比热容矩阵;KT为热传导矩阵.
  在摩擦片上施加压力边界条件模拟气胎加压,并约束摩擦片最外层节点的轴向及周向自由度,使摩擦片只能沿径向移动,设定热机耦合所需要的分析材料参数见表2,取内鼓轮为刚性传热体,并设定角速度(可由上文动力学分析得到),摩擦片用三维热机耦合全积分单元分析.当稳定接合时摩擦力不做功,不生成热量.热机耦合分析时间为5 s,室温为20 ℃,不考虑接触换热,进行瞬态热机耦合分析.表 2 分析材料参数参数摩擦片内鼓轮弹性模量/GPa3206泊松比0.20.3密度/(kg/m3)2 0007 860比热容/(J/(kg·K))1 100420热导率/(W/(m·K))155接触传热系数/(W/(m2·K))9090由于气胎消除间隙需0.8 s,在进行热分析时,这段时间不产生热量.之后,摩擦片温度升高较快,在滑摩结束且主、从动端角速度相同即3.864 s时,摩擦片温度局部达65.83 ℃,见图4;在
  3.864~5 s内,离合器表面高温部分向低温部分传热,在5 s时最高温度只有52.74 ℃.取表面3个节点分析发现,在刚接合时节点温度迅速升高;但在接合未结束时节点温度有所下降,这主要是由于温度升高过快,与周围温差大,热传导作用大于摩擦生热作用,温度下降;到滑摩结束时,由于不产生热量,温度依然下降,摩擦片节点温度时程曲线见图5.图 4 3.864 s时摩擦片的温度分布云图 图 5 摩擦片节点温度时程曲线内鼓轮的温度升高较慢,在滑摩结束即3.864 s时,摩擦片的温度局部最高只有22.67 ℃,见图6;在3.864~5 s时,离合器表面高温部分向低温部分传热,在5 s时最高温度只有21.91 ℃.取表面3个节点分析发现,在刚接合时各节点温度迅速升高,但在接合未结束时温度就有所下降,这主要是由于温度升高过快,与周围温差大,热传导作用大于摩擦生热作用,温度下降,到滑摩结束时由于不产生热量,温度依然下降,见图7.图 6 3.864 s时内鼓轮温度图 7 内鼓轮节点温度时程曲线图 8 离合器接合过程中摩擦力做功时程曲线离合器接合过程中摩擦力做功时程曲线见图8,摩擦力做功转换成热量,在稳定接合后共做功
  148 339.6 J.为分析离合器换热过程,将稳定接合即3.864 s时各节点的温度作为初始温度条件,分析离合器传热,分析时间为1 800 s,时间较长,可近似为稳态过程.不考虑与环境的换热,考虑摩擦片和内鼓轮的接触换热,接触传热系数见表2.1 800 s时离合器温度场分布云图见图9,可知温度梯度从内鼓轮无侧面端向有侧面端逐渐变小.图 9 1 800 s时离合器温度分布云图
  摩擦片上节点温度时程曲线见图10,发现摩擦片表面温度下降非常快,这时主要为摩擦片的热传导和与内鼓轮的换热作用,在30 s内温度已下降到30 ℃以下.
  (a)时程曲线(b)局部放大图 10 摩擦片上节点温度时程曲线
  内鼓轮节点温度时程曲线见图11,发现在开始阶段温度下降也非常明显,在50 s后有一定上升,这主要是由于摩擦片的热交换作用大于表面节点的热传导作用;随着摩擦片表面温度的迅速下降,温度逐渐趋于稳定.(a)时程曲线(b)局部放大图 11 内鼓轮节点温度时程曲线
  4 结 论
  通过动力学软件与有限元软件相互辅助分析,利用动力学软件在处理约束关系强和有限元软件在处理变形和温度场强的特点,对离合器的滑摩过程和温度场进行分析,得离合器从摩擦片和内鼓轮刚接触到稳定接合时间为3.864 s,在不考虑换热时摩擦片瞬态最高温度为65.83 ℃.该方法具有一定的实际意义.参考文献:[1] 张嗣伟. 石油钻机的气动摩擦离合器[M]. 北京: 石油工业出版社, 1980: 3034.
  [2] BALZS C, KROLY V, ALBERT A, et al. Fe thermal analysis of a ceramic clutch[J]. Tribology Int, 2009, 42(5): 714723.[3] 徐晓科, 莫易敏, 汤春球, 等. 微型汽车离合器摩擦副温度场研究[J]. 现代制造工程, 2009(1): 7477.[4] 林腾蛟, 李润方, 杨成云, 等. 湿式摩擦离合器瞬态热传导过程数值仿真[J]. 机械科学与技术, 2003, 22(1): 3941.
  [5] BABU A S, PRASAD N S. Coupled field finite element analysis of railway block brakes[J]. Proc Inst Mech Engineer, Part F: J Rail & Rapid Transit, 2009, 223(4): 345352.[6] 杨勇强. 矿用防爆机车湿式摩擦离合器的优化与仿真研究[D]. 西安: 西安科技大学, 2005.(编辑 于 杰)
  HyperWorks在日立300 t大型矿用卡车研发中的应用
  由日立重卡制造公司设计并在其加拿大安大略省圭尔夫市工厂生产的大型矿用卡车,体宽9.144 m
  (30英尺),能够运载重达316 t的有效负载.如何降低材料成本,达到ISO标准并获得客户满意的有效负载性能,是日立重卡制造公司的设计工程师在设计过程中面临的重大挑战.
  2005年,时任日立重卡制造公司技术分析经理HyperWorks专家Harry TEMPELMAN把HyperWorks带到日立.他把工作重点集中在使设计方案满足ISO 3471防翻滚保护系统规范(ROPS)和防落物保护系统规范(FOPS)的同时,降低驾驶室焊接结构钢的重量.
  通过很短时间内的多次设计和验证,日立重卡制造公司在计算过程中获得1个性能超出ISO标准的卡车结构.物理测试表明,HyperWorks能够准确预测此结构所有的非线性行为.
其他文献
摘要:为对构音异常声道实行辅助治疗,使用Matlab对声道进行动态仿真。采用极点形式在牛顿插值的基础上进行迭代计算出共振峰,对声道进行非线性分段,使用Matlab实现分段后声道形状显示;使用参数Corana改进算法,优化9个声道参数。在进行大量迭代计算获得合适的声道参数后,显示优化后的声道形状。该形状可以为构音治疗师的治疗提供参考。使用声门激励模型产生用于激励的声音,能够合成用于反馈的声音。实验证
期刊
LMS国际公司振动噪声、疲劳试验和仿真解决方案技术交流会在上海召开    2009年4月22日,LMS国际公司振动噪声、疲劳试验和仿真解决方案技术交流会在上海召开,来自泛亚汽车技术中心,上海大众、上海汇众、克莱斯勒、联合汽车电子、汉高、上海海马、上海微小卫星、上柴、常柴、中船等企业、和高校的100多位科研人员参会,《计算机辅助工程》作为支持媒体应邀参会。  会上,石银明博士首先对LMS国际公司及其
期刊
主题:引领自主创新,促进产业升级  时间:2009年7月28-29日  地点:甘肃兰州  主办单位:中国机械工程学会机械工业自动化分会、中国力学学会及其促进工程应用与产业结合工作委员会、中国软件行业协会数学软件分会、陕西省国防科技工业信息化协会  会务承办:北京诺维特机械科学技术发展中心  本届会议适逢CAE年会创办5周年,将由来自科技主管部门、科研院所、制造企业、高等院校的领导、两院院士、知名专
期刊
一年一度的CAE业界盛会——中国CAE工程分析技术年会于2009年7月28—29日在兰州成功举办. 本届年会由中国机械工程学会机械工业自动化分会、中国力学学会、中国力学学会促进工程应用与产业结合工作委员会、中国软件行业协会数学软件分会以及陕西省国防科技工业信息化协会主办,北京诺维特机械科学技术发展中心承办. 来自航空、航天、兵器、船舶、汽车、机械、土木、材料等领域的350多人参会. 《计算机辅助工
期刊
摘要:为提高飞机燃油系统管路设计效率,基于CAA二次开发技术以及CATIA管路模型几何与拓扑信息,开发具有3层架构的飞机燃油系统管路设计仿真软件. 该软件通过中间模型关联CAD和CAE软件,解决它们之间反复的数据传递和重复建模等问题;利用CATIA V5 Tubing Design工作台开发管流分析功能模块,实现飞机燃油系统管路建模时的在线流动分析.实例表明,该软件可提高燃油系统管路设计效率,具有
期刊
摘要:为解决传统二维光路设计和传输仿真中三维可视化不足的问题,通过增强光路设计的三维交互能力,提出基于Open Inventor的可视化三维激光光路布局设计和传输仿真. 借助Pro/ENGINEER完成光学元件和光路模型的建模,以Open Inventor为三维可视化仿真开发平台进行光学元件和光路模型及模型定位功能的管理,实现光学元件和光路模型的三维可视化布局,建立传输仿真场景数据库,并采用基于粒
期刊
摘要:为研究飞艇蒙皮PBO织物的编织单胞并优化编织方案,假设单胞为宏观均匀的复合材料,基于MSC Patran,借助PCL建立参数化2D平纹织物模型,并分析编织纤维间隙、编织纤维截面夹角和树脂层厚度等参数与面内刚度的关系.结果表明,降低编织纤维间隙、减小编织纤维截面夹角、减薄树脂层厚度,可以提高织物面内刚度. 该模型及其模块可以拓展MSC Patran的功能,所开发的程序可为该类材料的编织设计提供
期刊
摘要:为提高结构强度计算效率,增强飞机等结构强度计算工作的协同性,促进强度设计知识和经验的传承,基于CAESAM开发环境构建StrenBox结构强度自动化计算软件平台. 就该平台的基本理论和方法进行阐释,并以1个飞机结构的典型舱段结构为例进行实际应用的演示和说明. 最后,根据飞机等结构强度设计和研究的当前状态,对结构强度自动化计算软件平台的进一步工作进行展望.   关键词:结构强度; 软件平台;
期刊
摘要:为扩展B样条曲线,提出1种均匀二重混合双曲多项式B样条曲线. 该样条曲线在span{sinh t,cosh t,tsinh t,tcosh t,1,t,…,tk-6,tk-5}空间上均匀产生,其中k是大于等于5的整数. 证明k阶二重混合双曲多项式B样条基的性质和二重混合双曲多项式B样条曲线的性质. 二重混合双曲多项式B样条曲线精确地包含双曲多项式B样条曲线. 给出这种新样条曲线的细分公式并证
期刊
摘要:针对粗糙集只能处理量化数据,容错和推广能力较差的缺点以及BP神经网络的维数灾难问题,提出1种基于信息熵的粗糙集属性离散化方法. 该方法利用粗糙集对属性进行约简,解决BP神经网络的维数灾难问题,并将BP神经网络用于模式分类补偿粗糙集属性约简用于模式分类时的不足. 实例分析表明该方法具有较好的故障诊断效果.   关键词:粗糙集; 信息熵; 离散化; 神经网络; 故障诊断  中图分类号:TP212
期刊