V8发动机曲轴的动态特性研究与结构优化（可编辑）

V8发动机曲轴的动态特性研究与结构优化（可编辑） 密级: 分类号: 编号:: 工学硕士学位论文 发动机曲轴的动态特性研究与结构优化 硕士研究生:王艳玲 :马星国教授 指导教师 学科、专业:机械及理论 ’ 、、, 沈阳理工犬字 厂? ?? 年月。 ?,?鑫噼,? ;,,?;;.、?。。 分类号: 密级:: 编号: 工学硕士学位论文 发动机曲轴的动态特性研究与结构优化 硕士研究生: 王艳玲 : 指导教师 马星国教授 : 学位级别 。二硕卜 : 学科、专‘业 机械设计及王单论 : 所在单位 机械:程学院 论文提交划: 年 , 年月沦文答辩期: 学位授予单位: 沈阳理工大学:‘ ’ ..:.:. ’? ’’: ’ 硒:: ’。 : ‘, :; 【.\.: ?硕士学位论文原创性声明 本人郑重声明:本论文的所有工作,是在导师的指导下,由作者本 人 独立完成的。有关观点、、数据和文献的引用已在文中指出,并与参 考文献相对应。除文中已注明引用的内容外,本论文不包含任何其他个人 或集体己经公开发表的作品成果。对本文的研究做出重要贡献的个人和集 体,均己在文中以明确方式标明。本人完全意识到本声明的法律结果由本 人承担。 作者签字: 埘新盼 日 日期:沙碑弓月 学位论文版权使用授权书 本学位论文作者完全了解沈阳理工大学有关保留、使用学位论文的规 定,即:沈阳理工大学有权保留并向国家有关部门或机构送交学位论文的 复印件和磁盘,允许论文被查阅和借阅。本人授权沈阳理工大学可以将学 位论文的全部或部分内容编入有关数据库进行检索,可以采用影 印、缩印 或其它复制手段保存、汇编学位论文。 保密的学位论文在解密后适用本授权书 一 签 学位论文作者签名: 教 币 钐论。 指日 匕日 日 期: 埘沙 名期 “ 弘伽\,\,、:,:,;,::::::;;:.?;:.,:,;,目录 目录 第章 绪论? .课的目的与意义??一 .发动机曲轴的研究和应用概况一 ..发动机曲轴的研究概况..发动机曲轴国内外研究概况 .论文的主要研究内容?一 ..研究目标..研究内容..关键性问题? ..技术路线第章多体动力学基础理论一 .多柔体系统动力学的基础理论. .. 多柔体系统中的坐标系? ..相对变形场的描述 ..变形体上任一点的位移、速度和加速度??.. ..多柔体系统动力学中广义力?一 .多柔体系统动力学方程 ..运动约束方程?一 .. 系统方程一 ..多柔体系统动力学方程的求解? 第章 曲轴轴系动力学仿真模型的建立.. .曲轴系动力学模型机构简介? .曲柄连杆机构的运动学分析? ..活塞和连杆的运动计算?.沈阳理工大学硕士学位论文 ??????????????????????????????????????????????????????? ?????????二????????一一 一 .曲柄连杆机构动力学分析..曲柄连杆机构的质量分布. ..曲柄连杆机构中的作用力. .有限元模型的建立..几何模型的导入?. ..单元类型的选择? ..定义材料属性..有限元网格的划分 .有限元模型设定 ..设定主自由度和主自由度节点?. ..曲轴有限元模型验证??. .曲轴系动力学模型的建立.. 建模实现 .施加主动载荷及约束? .曲轴系仿真结果 .本章小结 第章曲轴有限元分析及结构优化 .载荷状况确定? .边界条件的处理 ..轴颈分布载荷计算 ..曲轴旋转惯性力和重力? ..相邻曲拐传递的扭矩..曲轴约束条件处理 .基于动力学仿真的曲轴动态应力计算 ..应力分析? .曲轴强度校核? ..曲轴疲劳强度安全系数校核.曲轴优化 ..优化结果分析??..本章小结~ 结论与展望?一? 参考文献??一“ 攻读硕士学位期间发表的论文和获得的科研成果一, 致 谢?一~ 、沈阳理工大学硕士学位论文 摘 要 曲轴是柴油机中最重要的零部件之一,它承受复杂、交变的冲击载荷,从而使曲 轴有着复杂的动态特性,其动态特性的好坏直接影响到柴油机的工作可靠性、经济性 及整机寿命.因此对曲轴的动态特性进行理论和实际的研究分析具有重要的现实意 义。传统分析方法的简化难以满足实际的需要,本文运用有限元法和多体动力学仿真 相结合的方法,建立柴油机曲轴轴系在实际工况下的多体动力学模型,进行曲轴 的动力学仿真计算,研究曲轴在即定工况下的实时动态响应,确定危险部位及极限应 力,并进行结构优化。 本文应用 软件和多体动力学仿真软件,有限元分析软件, 构造了一个柴油机曲轴轴系曲柄连杆机构的虚拟样机模型。基于这个虚拟样机模 型,用软件进行了运动学、动力学分析,计算出各构件在任意时刻、任意位置 的运动规律位移,速度,加速度等和构件问的作用力,得到了发动机工作过程中作 用在曲轴上的动态载荷,运用有限元分析软件对曲轴进行较为精确的应力、变形分析, 从而得到曲轴在整个工作行程中的响应,并给出了动态仿真过程。从而可以找到曲轴 工作时的危险部位,通过强度校验对曲轴进行优化。 研究结果表明,利用虚拟样机技术,结合有限元分析手段可以完成发动机曲轴的 动力学响应分析工作,获得了很好的仿真结果,同时对于发动机曲轴的改进设计,提 高发动机设计水平及提高发动机整机性能有着重要意义,是既经济又有效的科学化手 段。 关键词:曲轴,多体动力学,动力学响应,仿真 ,, . ., 科 . .?. , .. . ?, .,, , ,, . .. . ., .第章绪论 绪论 第章 .课题的目的与意义 曲轴作为发动机的重要部件之一,它的受力情况较为复杂,是发动机设计的重点 和难点。传统的设计方法的处理难以满足实际的需要,而通过运用三维实体造型软件 以及有限元软件建立模型,则可以精确分析曲轴的承载。模态分析用于确定设计机构 或机器部件的振动特性,即结构的固有频率与振型,它们正是承受动态载荷结构设计 中的重要参数。通过对曲轴的分析研究,从而对发动机的动态特性有了深入的了解, 掌握动态工作载荷的作用规律,对系统进行响应分析和评估,可以从中找出结构在动 态性能上所存在的问题。如果把其动态特性的提高作为重要目标进行结构的动念设计 和优化,则可将发动机设计由经验、类比和静态设计方法改变为动态、优化的设计方 法,从而提高产品的设计水平。由此可见,对于该课题的研究是具有重要的现实意义 的。 随着科技的发展和社会的进步,柴油机正朝着轻质、高速、大马力的方向发展, 曲轴是柴油机的关键部件,曲轴的质量特性直接影响了柴油机的使用寿命和安全可靠 性。对曲轴进行强度分析和有限元分析就成为新产品开发设计中必不可少的环节,因 此,对柴油机进行多柔体动力学分析,对新产品的开发设计具有一定的指导意义,对 提高柴油机的使用性能有十分重要的意义。也为新产品设计人员能够设计出高质量、 高性能的柴油机提供了理论依据。 曲轴是柴油机的主要部件之~,其尺寸参数在很大程度上影响并决定着柴油机的 整体尺寸和重量,而且曲轴是承担柴油机中热能向动能转换的关键部件,同时柴油机 产生的功率由它导出,因此在柴油机的设计与改进中曲轴占有极其重要的地位。曲轴 还是柴油机中受力最复杂的零件之~,它承受着气缸内气体作用力、往复运动质量惯 性力及旋转运动质量惯性力引起的周期性变化载荷并对外输出扭矩,因此曲轴承受交 变的拉伸、压缩、弯曲和扭转的复合应力,同时还有弯曲振动和扭转振动。在周期性 交变载荷作用下,曲轴系统可能在柴油机工作频率内产生共振、附加动念应力,致使沈阳理:【大学硕士学位论文 曲轴过早的出现弯曲疲劳破坏和扭转疲劳破坏。因此研究曲轴的动态性能,掌握曲轴 在工作过程中的振动规律以及载荷的变化规律,对提高曲轴的可靠性和寿命具有十分 重要的意义。 传统的机械设计效率低的主要原因就是,只有等到根据产品设计制造的样品 经过验证后,才能确定产品设计方案是否合理可行。这就大大延长了产品设计周期以 及造成了更多不必要的浪费。现阶段随着计算机的发展,利用计算机对其进行动态仿 真,可以进行虚拟设计仿真,这样可以缩短设计周期、提高设计质量。本课题运用多 体动动力学理论和有限元分析技术结合,建立多柔体动力学模型,通过模拟接近实际 工况的边界条件下进行动力学仿真,获取模拟一个工作周期的仿真结果,将其作为边 界条件在有限元软件中进行瞬态响应分析,从而获取曲轴的接近 实况下的动力学响应 和力学特性。为后续的曲轴结构优化提供了指导依据。曲轴是一异形长轴类零件,具 有轴线不连续,长径比大,形状和结构比较复杂等特点,决定其实际结构的几何参数 很多,且由于其受到冲击、扭转及弯曲应力和扭转应力等多种载荷的综合作用,所以 在实际工作中各个部位的不同时刻必然呈现出复杂的运动形态。因此构建合理的多柔 体动力学模型,进行精确的多柔体动力学分析以及更加精确的瞬态分析,为柴油机曲 轴的设计和优化提供理论参考和现实依据,具有十分重要的现实意义。本课题试图运 用结构动力学、模态分析理论、有限元理论及多体动力学分析理论,并有效利用专业 计算机软件研究曲轴的动态特性。 综上所述,鉴于曲轴作为柴油机最重要的运动部件,本课题结合工程实际项目对 曲轴连杆系机构进行多体动力学分析,并对曲轴进行了实时的动力学响应分析。利 用这种研究方法可以高效的对构件进行分析,并实时的获得想要的结果和数据,对缩 短设计周期,提高分析精度、设计水平和市场竞争力都具有重要的现实意义。因此, 本课题的研究具有纵深的理论研究背景和广阔的应用前景,具有十分重要的理论和实 际工程意义。 .发动机曲轴的研究和应用概况 ..发动机曲轴的研究概况 目前,在对发动机曲轴进行动态响应分析过程中,由于实际工况中曲轴承受 活塞、连杆传递的爆发压力的交变载荷作用,受力情况及其复杂,采用传统的单纯有第章绪论 限元分析方法很难完成对曲轴运行过程中动态变化的边界条件描述。一般采用局限于 某个或某几个特定时刻的,过度简化的曲轴加载边界条件,其计算结果往往与实测 偏差较大,不能真实地反映曲轴在实际运行工况下的力学特性。同时本文中所分析的 发动机曲轴为引进机型,设计经验和现有计算公式严重缺乏,单纯依靠传统方法完 成曲轴实时动态响应的分析工作是不可能的。 ..发动机曲轴国内外研究概况 在年代末年代初至今,国内外在应用有限元分析技术对曲轴进行 分析方面 取得了大量的成果。。,且利用有限元法进行曲轴强度分析的工作已经取得了许多有用 的成果,如文献?用有限元程序对曲轴进行了强度校核,计算了疲劳安全系数;文献 ?运用有限元分析技术对柴油机曲轴进行了强度分析,从而探索了一条连续梁 法和有限元法相结合的新方法;文献口们研究了曲轴的弯曲应力受其结构参数变化引起 影响;文献?蜘对运用有限元分析方法进行了曲轴圆角形状的优化;文献?”对曲轴用 有限元法进行了模态分析,大学的..对曲轴做了有限元分析, ..等研究了曲轴的动态响应引;德 并实现了声振耦合分析“; 国大众公司.和.对某直列四缸水冷柴油机机体作了静刚度和模态 分析引,以上这些研究成果对于应用有限元分析技术对曲轴进行强度分析方面都起到 了极大的推动作用。但也都明显一些不足之处,他们对于边界条件的处理过于简化。 曲轴是一个复杂的构件,在每一个部位的不同时刻都在承受着不同的载荷,这些因素 在以往的研究工作中并没有充分的考虑,所以分析结果也必然产生很大的误差。但把 单拐三维模型的计算结果与曲轴的空间刚架模型相结合,对提高后者的计算精度还是 大有裨益的。 我国在应用有限元分析技术对发动机研究方面也取得了大量的成果。特别是在大 型柴油机领域,有限元分析技术已经在很多方面得以应用?。。旧。在曲柄连杆机构的其 它部件如曲轴、主轴承座、连杆、高压汽缸等在设计开发阶段,都有运用有限元分析 技术进行结构强度及疲劳强度分析的实际经验。其中也不乏对零部件动态情况下的有 限元分析。南京理工大学的王良囤、胡德波通过建立曲轴的三维几何模型基础上, 应用有限元分析技术进行了交变载荷作用下的疲劳强度分析,并通过试验对比验证了 其可靠性 。西安交通大学的王迪生、张勇运用有限元分析软件,通过编制计算机程 序,对武汉理工大学硕士学位论文压气缸的静态及动态应力进行分析,为探讨高压气沈阳理工大学硕士学位论文 缸断裂机理和容限寿命受动态应力分析影响方面提供了理论基 础和研究方法【;河南 工程学院的李立琳、吴素珍建立了柴油机油底壳三维模型,运用有限元分析 技术进行了自由模态分析,为后续进行的动力学响应分析奠定了基础,同时的有限元 分析模型,对其进行了自由模态实验,为柴油机下一步的动力响应分析提供一条较为 精确的方法,同时对于整机的动态特性分析和结构的优化设计莫定了理论基础,也为 柴油机实际结构的改进提供了有效的途径乜铂:沈阳工业学院李松波等人也对曲轴进行 了模态分析,并找出了曲轴振动中的危险区域乜;太原理工大学应用有限元分析技术 对型柴油机曲轴进行过振动响应分析乜引;南通柴油机股份有限公司姜峰等 人对柴油机曲轴进行了模态分析实验,并在此基础上进行了模拟特定工况下的 受迫响应分析和计算,得出了曲轴动态特性参数和特定工况下的响应变形,据此对 柴油机曲轴进行了修改乜。李梅乜采用整体曲轴进行应力计算,并建立了主轴 承座的有限元模型,在主轴颈处定义接触,在连杆轴颈处施加面 载荷。胡金蕊乜旷删等 综合运用有限元分析技术与多柔体动力学分析技术,根据模拟仿真一个工作周期后得 获得的动态载荷作为边界条件,在有限元软件中将其加在曲轴表面节点上,通过分析 计算得到了曲轴的动态应力分析状况。其计算结果更加接近实际。但是仍然没有考虑 曲轴一轴承系统摩擦、刚度和强度耦合问题对曲轴应力的影响。 综上所述,运用多体动力学方法将是曲轴动态强度研究的主要方向,通过对曲轴 进行的模态分析,以获得该曲轴的固有频率和振型,分析曲轴的结构参数对曲轴动态 特性的影响。特别是对曲轴进行结构瞬态响应分析以及由此引起的整个曲轴的应力变 化。 .论文的主要研究内容 ..研究目标 建立考虑弹性支撑作用下的曲轴在实际工况下的动力学模型,研究曲轴在即定工 况下的实时动态响应,确定危险部位及极限应力,并进行结构优化。 ..研究内容 建立包括曲轴、连杆、活塞在内的曲柄连杆机构的三维几何模型。 建立多刚体动力学模型并仿真。 建立曲轴和机体的有限元模型。并生成软件所需要的柔性体模态中性文件。第章绪论 进行曲轴的模态分析,分析曲轴的固有频率和振型,研究曲轴的结构参数对曲轴 动态特性的影响。 进行多柔体动力学分析,并输出载荷文件。研究曲轴的不同部位瞬态受力情况, 为曲轴的进一步瞬态响应分析提供边界条件。 对曲轴进行结构瞬态响应分析及结构优化。 ..关键性问题 简化数字模型。 正确建立曲柄连杆系统的柔性体模型。 正确施加复杂的载荷。 ..技术路线 本文综合采用虚拟样机技术、柔性多体动力学理论和有限元分析方法对曲轴进行 动态分析,更加精确地进行了曲轴在实际工况下的实时动力学响应分析。 应用软件,建立包括曲轴、连杆、活塞在内的曲柄连杆机构的三维 几何模型,对模型进行必要的简化,忽略对动态特性分析影响不大的特征。 将建立好的曲轴和轴承三维几何模型分别导入软件中,设置好相应的 单元类型和材料属性,进行网格化分。 将处理好的有限元模型导入软件,建立主自由度节点,并建立刚性层。 并生成/所需要的模态中性文件。 进行模态对比,确保引入后,与原来有限元模型的动态特性一致。 应用/软件建立曲柄连杆机构的多刚体动力学模型后,将曲轴和轴承的 刚性体部件替换柔性体部件,完成多柔体动力学建模。 将根据示功图拟合出的燃气爆发压力分另:在个活塞上。 进行多柔体动力学仿真。仿真时问设置一个周期.,步,每隔度提取 一次数据。 将多柔体动力学仿真后所产生的载荷文件按术.格式导出。 应用软件将已经处理好的有限元模型进行应力恢复后,将多柔体仿真所产 生的术.载荷文件以命令流的形式加载,进行曲轴的瞬态分析。 根据应力分析结果进行结构优化。图?建模软件、有限元软件和机械系统动力学分析软件的数据图第章多体动力学基础理论 第章多体动力学基础理论 虚拟样机技术的核心理论是多体系统动力学,多体系统是指由多个物体通过运动 副连接的复杂机械系统。多体系统动力学是由多刚体系统动力学与多柔体系统动力学 组成的。多体系统动力学的根本目的是应用计算机技术进行复杂机械系统的动力学分 析与仿真。研究这些系统的动力学需要建立非线性运动方程、能量表达式、运动学表 达式以及其它一些量的公式。多柔体系统动力学的研究对象是由大量刚体和柔体组成 的系统。 。多柔体系统动力学的基础理论 .. 多柔体系统中的坐标系 研究多柔体系统,首先要选择合适的坐标系。图.中的五而五为惯性坐标系, 该坐标系不随时问而变化。为了描述柔体的运动,在柔体上建立一个动坐标系葺:玛。 该动坐标系可以相对惯性坐标系进行有限的移动和转动。描述该坐标系移动和转动的 坐标称参考坐标。 幽.单个柔体的坐标系和坐标描述 图.柔体的浮动坐标系 对于刚体而言,只要把动坐标系固定在刚体上的某一点,那么整个刚体的空问位 置就可以由动坐标系来确定,因为刚体上的任意一点相对动坐标系的位置都不发生变 化。与刚体不同,柔性体是变形体,体内各点的相对位置时时刻刻都在变化,只靠动沈阳理工大学硕士学位论文 坐标系不能准确描述该柔性体在惯性坐标系中的位置,因此,引入弹性坐标来描述柔 性体上各点相对动坐标系统的变形。这样柔性体上任一点的运动就是动坐标系的“刚 性”运动与弹性变形的合成运动。由于柔体上各点之问有相对运动,所以动坐标系的 选择不是采用连体坐标系,而需要采用随着柔性体形变而变化的坐标系,即“浮动坐 标系”。图.为柔体采用浮动坐标系的例子。其中动坐标系五,.,始终沿杆的切 线方向,动坐标系,.沿杆的割线方向。 在研究多柔体系统时,合适的坐标系是非常重要的。在确定浮动坐标系时有两点 准则:、便于方程建立求解;、柔性体刚体运动与变形运动的耦合尽量小。目前常 见的浮动坐标系大致有如下种,局部附着框架、中心惯性主轴框 架、蒂斯拉德框架、 巴克凯恩斯框架以及刚体模态框架。采用何种需因实际情况而 定。图.所示为多柔 体系统与多刚体系统的区别与联系及其坐标的选择。 图.多柔体系统与多刚体系统的关系 对于柔体上的某一点,相对惯性系的位置矢量可用下式表示略去 各矢量的上 脚标 ,月么万 ? 式中一动坐标系原点位置矢量; 万一点相对动坐标系的位置矢量; 一坐标转换矩阵。 对于变形体,点可表示为第章多体动力学基础理论 式中瓦一点没有变形时的相对位置矢量; 巧一点相对变形矢量。 所有各点的变形构成一相对变形场。 ..相对变形场的描述 对于变形体杆、梁、板等,变形场可用一组无穷级数来表示,即 啄。?吼五,六五,恐,屯 弓:?包反,。鼠五,艺,恐 七 巧,?,纹。,而,屯 式中订,、历,、订,一构件上一点,恐,毛的变形分量; %、反、一依赖时间的系数,也称“坐标”; 六,。,‰一基函数,它们应满足变形体的边界条件。 如果上述级数收敛,那么?式就代表实际的变形场。无穷级数给计 算带来了 困难,为了减少自由度,可用有限级数来代替?式,这类似于模态 分析中用头 几阶振型代替无穷级数的做法。这种离散化的方法也称为瑞利?? 罩兹法,即 , 乃。圭?吼五 七 』 圭 /\ 瓯鼠 一吩 圭 纹 一吁 。?料。?料 从理论上可以证明,只要,,足够大,方程?完全可以得到足够精度 的近 似。公式可更简洁的写成: 历,, 式中 历广嘛万门万,;缈一弹性坐标矢量,元素为%,包和。变形体上任 一点的位移、速度和加速度 由?、?及?式可得变形体上任~点的位移 名尺彳瓦, ? 对上式求导可得速度表达式 易月彳玩彳‘:影 公式?中第二项可进行下述处理。 ?琢:矽 矽为矢量,其元素为动坐标系旋转坐标对时间的导数: 为矩阵,定义为: ~ 融时?私 ,,为旋转坐标的个数。 这样?式可写成更简洁的形式 哆 ,曰 ~ 彳,厂‘雾三亘 式中,是×单位矩阵,口矗矽 口,,三:,彳。 对?式进~步求导可得加速度矢量 毒硅 一 ?。多柔体系统动力学中广义力 ’ 。,等量柔体系统动力学分析中需考虑的广义力有弹性变形广义 力、外力引起的广义 力和约束的广义力。 弹性广义力 矗性力引起变形体的虚功可表 ~由弹性力学的理论可知,对于各向同性的材料, 不成式中仃、为应力和应变矢量。 应变可表示变形矢量巧的函数 一 面 工巳阵 乘子矧 为微分 厂 式 旦甄 一 一% 一甄 一 一% 一 一%旦甄. 一阢 一帆一畹, 由公式?可得 占, ? 为弹性系数对称矩阵。所以有 ? 把?和?式代入?式经整理得 ? 形,万,,, ? 刀,,, 为了后面推导方便,?式也可写成 ? 万形:一臼 ;,黾曼暑蓍 广义外力 作用在多体系统某一构件上的外力所做的功为:噬 ? 上式还可写成 孵鳞,嘲 约束广义力 .弹簧~减振器~力发生器单元 帅烹黧件?用;单拜减振器~力发生器静‖点之问连接,,/分另 为弹簧刚度、阻尼和激振力。%为弹簧原长。尸和点的相对位置 矢量可表示为 ,尺彳‘万?’一彳万。 ~ 。、. , 式中, 乞,其长度为 ,历蕊可以证明长度对时问的导数为 ? 式中声为沿尸,的单位矢量,,为 矿矗矽彳酱/一天』一/矽/一彳口; ? 或 白一 ? 虚功方程为 一。 ? 式中为弹簧~减振器~力发生器的作用力。 尼,一%户 ? 虚位移万,为第章多体动力学基础理论 ? 已哇一《 所以虚功表达式可进一步变成 巧形一矿‘’ 。丑 式中 盯一只声 ;‖一只声’。 .多柔体系统动力学方程 ..运动约束方程 在多体系统中并不是所有广义坐标都是独立的,各构件之间由各 种铰链联接在一 起。例如悬架中常用的柱铰、球铰、万向节铰、滑移铰等。在列 写约束方程时要考虑 弹性变形,这是与多刚体系统的主要区别。最后的形式为一组非 线性方程组 ? , 式中玎??矿是系统全部广义坐标,为广义坐标个数; 。??。、,是约束方程的个数。各约束方程问要求线性独立。 对方程?求导,可得由广义坐标表示的速度并;速度方程 ? 口 ? 牙一,曹,牙,幻 式中下标记号为偏导数符号。 .. 系统方程 对第个柔体或刚体,拉格朗日方程形式为 蕊 。,. ,?、 协, 瓦四一四呼以?/ 硇‘/ 式中:广义力包括单元弹性变形引起的和外加载荷包括单元问载 荷引起的 广义力两部分。兄为拉氏乘子。由?和式得 一 鼍 由公式~不难得&::沈阳理工大学硕士学位论文 础, 丢雾等穹砧虿一毒?丁虿,’ 令 沪岍将矿四丁 该项也称为速度二次项,其中包含陀螺力和哥氏力引起的力学效 应。把上面结果 代入?式得 ? 江,,??,% ‘百刀?笋允砭残, 一? 垅,艘 《 吻嘛 《 略嘞略 ~~,彩 酶 .尼伊酊 ? ?,????,??? ? ?................................... ????????? 一 五憾批孙撕。, 一旧一旧旧一 ~ 臣,中 一其 都 体 以上方程与约束方程?一起构成多体系统动力学方程,它是由微分方程组 和非线性几何方程组组成的。全系统方程由上述各方程组装而成。要想得到上述方程 组的封闭解往往是十分困难的,一般只能采取数值解法。 ..多柔体系统动力学方程的求解 对于由微分方程组和约束方程组构成的多体系统动力学方程的求解问题已有一些 研究成果。本文采用的方法是:首先利用运动学关系消除非独立变量,使微分方程只 含有独立变量,然后直接积分求解。求出独立变量的位移、速度和加速度后再通过运 动学关系求出非独立变量的位移、速度和加速度。积分过程采用?法。具 体做法如下 对形如下式的方程 《五 , , 【 首先对约束方程求二次导 ‘互一,一乞,一‘季。第章多体动力学基础理论 令 一一,亘一口。口 称为加速度方程的右端项。 , 孑《孑蔓妻一万 卜驯簪野 式中 。百 ;。『冬《;户曼一物『每 由上式可求出历和兄。按此方法求出的加速度,约束方程是自动满足的,因为方 尽管对方程可进行直接积分求解。但由于误差的积累,计算结果的误差 较大,而且很容易使计算过程发散。为了改进计算,提高精度,可采用下述方法。 首先对独立和非独立变量利用运动约束先进行分离,然后只对独立变量进行积分, 先把广义坐标写成『门的形式,其中为独立坐标,。为非独立坐标。 . 等价于 , ? 即 ? 式中 ~?‖????????????????塑三奎堂堡主堂垡笙壅 另外,对约束方程两次求导后得 电 、. 血唾一 计算时只对状态变量少。,.,求解,非独立变量及速度由约束方程 和? 式求出在动力学方程中,去掉速度和加速度项,得到静力分析基 本公式 ‰第 【 /令?,方程?变成 五一疋 ? 对?式取变分得 嘭五一兄国 ~ 应用变量分离方法可得 套 、? 由上式可得 岍引瑚吼 式中吃八 把?式代入?式,并简化注意到吃为零得 豆 ? 式中哀一醚 方程?即为静力平衡方程,由它可解出给定力下的平衡位置时的独立坐标 数值。对该非线性方程组的求解,可采用?方法。 动力学部分的计算如图.所示。图.计算框图 一 .沈阳理工大学硕士学位论文 第章 曲轴轴系动力学仿真模型的建立 为了进行曲轴系统多体动力学分析,曲轴轴系的动力学仿真模型是必需的,本章 详细分析研究了曲轴轴系仿真模型建立的方法及选择的理由。首先,利用在软 件中建立曲柄连杆机构的多刚体动力学模型。在软件建立曲轴和轴承座的 三维实体模型,根据分析需要选择合适的单元类型和单元尺寸在 软件进行 结构离散:然后在软件设定主自由度节点并生成模态中性文件。导入中 建立了曲轴轴系多柔体动力学仿真模型,为动力学仿真计算打下基础。 .曲轴系动力学模型机构简介 曲轴系动力学模型机构主要为曲柄连杆机构,包括活塞连杆组、曲柄飞轮组等内 燃机的主要运行件。其功用是将活塞的往复运动转化为曲轴的旋转运动,且将作用在 活塞上的燃气压力转化为扭矩,借助飞轮向外输出,从而实现热能向机械能的转化, 是内燃机传递运动和动力的机构。 曲轴 曲轴作为发动机的中枢,其通过连杆的上下往复运动变成循环旋转运动。 是发动机上的一个非常重要的机件,其材料是由碳素结构钢或球墨铸铁制成的。有两 个重要的部位:主轴颈、连杆轴颈。主轴颈被安装在机体上,连杆轴颈与连杆大头孔 连接,连杆小头孔与汽缸活塞连接,是一个典型的曲柄滑块机构。发动机工作就是, 活塞经过混合压缩气的爆炸,推动活塞做直线运动,并通过连杆将力传给曲轴,使曲 轴将直线运动变成旋转运动。曲轴的旋转是发动机的动力源,也是整个机械系统的源 动力。 曲轴是复杂的弹性连续体,其动态特性对发动机整体的动态特性有重要影响。曲 轴在工作过程中,要承受周期性变化的气体爆发压力、活塞连杆组往复惯性力、旋转 惯性力、附件不规则阻力矩及外界反作用力等多种载荷共同作用。在多种载荷及约束 的综合作用下,曲轴表现出复杂的运动形态。 在中曲轴被分为单个曲柄之后,整个曲轴活塞连杆机构就演变为个如图第章曲轴轴系动力学仿真模型的建立 .所示的曲柄滑块机构的组合。 图.曲柄滑块机构 其中:为曲柄,为连杆,为活塞。曲柄保留绕轴旋转自由度; 曲柄与连杆保留相对旋转自由度;连杆与活塞之间保留相对旋转自由度;活塞与机体 之间仅保留沿方向移动自由度。 个曲柄构件之间以哑元通过胡克铰链弹簧力和万向铰链连接。胡克副和万向副 约束两个构件之问的个平动自由度和个旋转自由度,两个构件之间有个旋 转自由度。哑元是一种特殊的构件,无几何形状、无转动惯量和质量。根据多体动力 学理论可以建立曲轴的动力学方程。 连杆 连杆是曲柄连杆机构中两端分别与主动和从动构件铰接以传递运动和力的杆 件。它是发动机中的重要零件,它连接着活塞和曲轴,其作用是将活塞的往复运动转 变为曲轴的旋转运动,并把作用在活塞上的力传给曲轴以输出功率。工作时承受三个 方面的力:活塞上的气体压力、活塞和连杆小头的往复惯性力以及连杆本身绕活塞销 做变速摆动产生的横向惯性力。由于上述力均是周期性变化的,且带有冲击性,因而 连杆是在拉、压和弯曲等交变应力下工作的。因此要求连杆具有足够的刚度、疲劳强 度和较小的质量。 活塞 活塞的构造可分为三部分,活塞顶部:它承受气体压力,它是燃烧室的组成 部分;活塞头部:作用是承受气体压力;安装活塞环,与活塞环一起密封气缸; 将热量通过活塞环传给气缸壁。活塞裙部:作用是对活塞在气缸内的往复运动起 导向作用,并承受侧压力。 活塞作为发动机的重要运动件之一,它所处的工作条件相当严酷,高温、高压、 高速、润滑不良的条件下工作。因此就要要求活塞要有足够的刚度和强度,传力可靠; 要有足够的刚度和强度,传力可靠;导热性好,耐高压、耐高温、耐磨损;质量小,沈阳理工大学硕士学位论文 重量轻,尽可能减小往复惯性力;热膨胀系数小,摩擦系数小等性能要求。 活塞与气缸盖、气缸共同构成燃烧室,是发动机做功的关键零件,通过活塞的往 复运动,完成发动机的进气、压缩、膨胀做功和排气过程。而在这些过程中,活塞顶 部气压的不断地循环改变,从而导致活塞的上下来回运动。这样内燃机燃料的热能就 转化成为活塞的上下运动的动能,而活塞就是充当这种转化的介质。活塞的功能包括 组成燃烧室,承受燃气作用力,并把它传给连杆:密封气缸,防止 燃气泄漏及润滑油 窜入燃烧室:将热量传给气缸壁:将连杆的侧压力传给气缸壁等。 活塞是发动机的“心脏”,承受交变的机械负荷和热负荷,是发动机中工作条件 最恶劣的关键零部件之一。活塞的功能是承受燃气压力和惯性力,并将燃气压力通过 活塞销传给连杆,从而驱使曲轴旋转对外做功。 .曲柄连杆机构的运动学分析 发动机的工作机构主要采用中心式曲柄连杆机构,其气缸中心线通过曲轴的旋转 中心,结构简图如图.所示: 图.中心曲柄连杆机构示意图 图中,一连杆长度,指连杆大、小头孔中心的距离: ,一曲柄半径,指曲柄销中心与曲轴旋转中心的距离: 口一曲轴转角,指曲轴偏离气缸中心线的角度: ‖一连杆摆角,指连杆中心线在其摆动平面内偏离气缸中心线的角度: 缈~曲轴旋转角速度: 一活塞位移,指活塞由上止点女向下比点运动的距离,上、下止点分别指活塞在 一,一第章曲轴轴系动力学仿真模型的建立 气缸中运动所达到的距离曲轴旋转中心最远、最近的位置。活塞 在气体爆发压力作用 下作往复运动,并通过连杆将往复运动转化为曲轴的旋转运动。 在曲柄连杆机构进行 运动学、动力学分析时,通常近似的认为曲轴作匀速旋转运动, 并将机构的各种运动 学和动力学参数表示为曲轴转角口的函数。 ..活塞和连杆的运动计算 由图.可得,活塞的位移 ,.一 /、 式中一一; 五一。 活塞速度 ? ,:,缈?? 活塞加速度 ‖ ? 缈 名 ?一 生出 竺吣 浠 口一‖ 式中可按牛顿二项式定理展开成如下级数: ‖:一口一五口~ 口一??. ’ 实际中为了计算方便起见,一般取前两项: 。‖?一三五凹 ? ? 式中. 口,称为活塞的一级加速度 ,允 口,称为活塞的二级加速度。 ? ‖ 一/?。‘。 连杆的角速翰刊器刊?蕊‘可沈阳理工大学硕士学位论文 。‖ 一九。 ‘口? 连杆的角加速如一矿加以嵩一矿篙黑器 .曲柄连杆机构动力学分析 ..曲柄连杆机构的质量分布 曲柄连杆机构的所有运动零件可按运动性质不同分为活塞组、曲 拐和连杆组三部 分,其质量分别用。、%和%表示。 通过对曲柄连杆机构的运动分析可知,活塞组作往复运动,曲拐 作回转运动,连 杆组作往复运动和回转运动的复合运动。连杆是作复杂平面运动的零件,为便于计算, 将连杆组用换算质量聊。和%来代换。。是集中作用在连杆小头中心处,并只做往复 直线运动的质量;,是集中作用在连杆大头中心处,并只做沿着圆周做旋转运动的质 量。据文献?。,工程上通常近似取: %.?.。,聊:.一.。 往复直线运动部分的质量//,可表示为活塞质量聊。与换算到连杆小头中心的质量 之和,即,% 不平衡回转质量,.可表示为曲拐质量///。与换算到大头中一““.”圭杆质量//:之和, 即聊,.%// ..曲柄连杆机构中的作用力 气体作用力 作用在活塞上的气体作用力。等于活塞上、下两面的空间气体压力差与活塞顶面 积的乘积,即 旷等?’ 四 式中:。??活塞上的气体作用力,; ??气缸内的气体压力,:第章曲轴轴系动力学仿真模型的建立 ’??曲轴箱内的气体压力,,一般取一.; 卜活塞直径,; 往复惯性力 在曲柄连杆机构中,往复惯性力为 ? 一%以一,缈 往复惯性力作用在气缸中心线上,它的变化规律和活塞加速度相 同,负号表示 的方向与活塞加速度的方向相反。 旋转惯性力 各种曲柄连杆机构的旋转惯性力均为 ? ,一, 式中,??旋转运动质量,,当?等于常数时,旋转惯性力,.的大小 是不变的, 其方向总是沿着曲柄半径方向向外。 作用在活塞上的总作用力 在活塞销中心处,同时作用着气体作用力;和往复惯性力乞,由于 作用力的方向 都沿着气缸中心线,故只需代数相加,即可求得合力合为 ? 合 式中合??活塞上的总作用力,; ??往复惯性力,; ;??活塞上的气体作用力,。 活塞上的总作用力合分解与传递 合分解成两个力:沿连杆轴线作用的力,和把活塞压向气缸壁的 侧向力,如 图所示,其中沿着连杆的作用力为: 脚合南尘 四七