机械毕业设计(论文)-机械零件自由模态数字试验误差分析【全套设计】.doc_第1页
机械毕业设计(论文)-机械零件自由模态数字试验误差分析【全套设计】.doc_第2页
机械毕业设计(论文)-机械零件自由模态数字试验误差分析【全套设计】.doc_第3页
机械毕业设计(论文)-机械零件自由模态数字试验误差分析【全套设计】.doc_第4页
机械毕业设计(论文)-机械零件自由模态数字试验误差分析【全套设计】.doc_第5页
已阅读5页,还剩40页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

机械零件自由模态数字试验误差分析扬州大学广陵学院本科生毕业论文 毕业论文题目 机械零件自由模态数字试验误差分析 学 生 姓 名 专 业 机械设计制造及其自动化 班 级 机械81001(广陵) 指 导 教 师 完 成 日 期 2014 年 6 月 2 日 摘要 近年来,随着计算机技术的普及和计算速度的不断提高,ANSYS有限元分析在工程设计和分析中得到越来越广泛的应用,已经成为解决复杂的工程分析计算问题的有效途径。本课题选用了平面结构作为研究的模型,利用ANSYS进行建模、确定边界条件,施加载荷,划分单元网格,进行自由模态分析模态分析。且改变单元结构、网格大小等因素,讨论其对模态分析的影响。一般情况下,零件都是具有约束的,所以求解具有约束情况下的模态更具有实际意义。选用不同约束的梁作为研究对象,求其固有频率,且讨论了约束刚度对固有频率的影响。最后对悬挂法测固有频率进行分析,讨论了吊装形式对分析结果的影响,研究其自由模态与实际情况下的试验误差。关键词:有限元,ansys,模态分析,约束刚度全套设计,加153893706 Abstract In recent years, with the popularization of computer technology and computational speed is increased, the ANSYS finite element analysis has been widely used in engineering design and analysis, has become the effective way to solve complex engineering analysis of computational problems. The project uses the plane structure as the research model, use ANSYS for modeling, boundary conditions, loading, dividing grid unit, the free modal analysis modal analysis. And the change of cell structure, the size of the grid and other factors, discuss its influence on modal analysis. In general, the parts are constrained, so the solving of modal constraint condition has more practical significance. Choose different restrained beam as the research object, for its natural frequency, and discusses the constraint stiffness influence on the natural frequency. At the end of the suspension method to measure the natural frequency was analyzed, the influence of hoisting form discussing the analysis results, the test error of the free modal and the actual situation.Key word:Finite element ,ANSYS,Modal analysis,Restraint stiffness目录摘要IAbstractII目录III第一章 绪论11.1 引言11.2 本课题研究内容31.3本章小结3第二章 振动研究的理论基础42.1 振动基础理论42.2 连续系统自由振动52.3本章小结7第三章利用Ansys软件进行模态分析83.1 ansys有限元软件简介83.2 ansys自由模态分析的基本操作93.3约束下梁的模态分析133.4本章小结19 第四章 利用ansys进行悬挂法分析仿真204.1无约束下箱体自由模态分析204.2悬挂法测试箱体固有频率244.3本章小结29致谢30总结31 参考文献32附录133附录234附录336附录438IV1 沈酝籍 机械零件自由模态数字试验误差分析第一章 绪论1.1 引言在数学中,有限元法(FEM,Finite Element Method)是一种为求得偏微分方程边值问题近似解的数值技术。它通过变分方法,使得误差函数达到最小值并产生稳定解。类比于连接多段微小直线逼近圆的思想,有限元法包含了一切可能的方法,这些方法将许多被称为有限元的小区域上的简单方程联系起来,并用其去估计更大区域上的复杂方程。它将求解域看成是由许多称为有限元的小的互连子域组成,对每一单元假定一个合适的(较简单的)近似解,然后推导求解这个域总的满足条件(如结构的平衡条件),从而得到问题的解。这个解不是准确解,而是近似解,因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解,而有限元不仅计算精度高,而且能适应各种复杂形状,因而成为行之有效的工程分析手段。早期的有限元主要关注于某个专业领域,比如应力或疲劳,但是,一般来说,物理现象都不是单独存在的。例如,只要运动就会产生热,而热反过来又影响一些材料属性,如电导率、化学反应速率、流体的粘性等等。这种物理系统的耦合就是我们所说的多物理场,分析起来比我们单独去分析一个物理场要复杂得多。很明显,我们需要一个多物理场分析工具。在上个世纪90年代以前,由于计算机资源的缺乏,多物理场模拟仅仅停留在理论阶段,有限元建模也局限于对单个物理场的模拟,最常见的也就是对力学、传热、流体以及电磁场的模拟。看起来有限元仿真的命运好像也就是对单个物理场的模拟。这种情况已经开始改变。经过数十年的努力,计算科学的发展为我们提供了更灵巧简洁而又快速的算法,更强劲的硬件配置,使得对多物理场的有限元模拟成为可能。新兴的有限元方法为多物理场分析提供了一个新的机遇,满足了工程师对真实物理系统的求解需要。有限元的未来在于多物理场求解。40 沈酝籍 机械零件自由模态数字试验误差分析近年来,随着计算机技术的普及和计算速度的不断提高,ANSYS有限元分析在工程设计和分析中得到越来越广泛的应用,已经成为解决复杂的工程分析计算问题的有效途径。现在越来越多的设计制造都离不开ANSYS有限元分析计算,例如机械制造、材料加工、航天航空、汽车、土木建筑、电子电器、国防军工船舶、铁道、石化等。因此掌握和研究ANSYS有限元分析对我们工作中的设计制造具有重要的意义。此外,在ANSYS有限元分析的应用中,建立的模型向大型化,计算结果向精确化发展,所以ANSYS有限元的进阶技术整体法与局部法的应用也越来越普遍化,也成为ANSYS有限元应用不可缺少的部分。本课题的目的在于巩固和拓展我们在校期间所学的基本知识和专业知识,训练我们综合运用所学知识,提高分析和解决问题的能力。灵活运用各种知识,把使用有限元分析软件ANSYS同巩固和提高自己已有知识统一起来,把所掌握的新软件同解决实际问题统一起来,全面提高我们的能力。通过运用有限元分析软件各种轴的应力分析,使我掌握了有限元分析软件ANSYS的基本知识及其基本操作,会用程序设计语言精确建立模型,确定边界条件,划分单元网格,施加载荷,以及对模型进行应力应变分析。同时也要学会如何使用子模型法对模型进行应力分析。课题的研究内容所涉及到有限元模态技术的发展,借助梁振动的微分方程求解并得到了悬臂梁各阶固有频率,计算不同约束梁的固有频率随约束刚度的变化等。在龙英,滕召金,赵福水所著有限元模态分析与发展趋势;陈忠所著滚动轴承及其支承刚度计算;李东旭所著高等结构部动力学等著作中均有提及或实例。1.2 本课题研究内容 本课题选用了平面作为研究的模型,利用ANSYS进行建模、确定边界条件,施加载荷,划分单元网格,进行自由模态分析模态分析。且改变单元结构、网格大小等因素,讨论其对模态分析的影响。一般情况下,零件都是具有约束的,所以求解具有约束情况下的模态更具有实际意义。选用不同约束的梁作为研究对象,求其固有频率,且讨论了约束刚度对固有频率的影响。最后对悬挂法测固有频率进行分析,讨论了吊装形式对分析结果的影响,研究其自由模态与实际情况下的试验误差。因为本课题进行模态分析,要不断改变单元尺寸、单元实常数等,这样在ANSYS中直接建模就比较麻烦。所以本课题中均采用ANSYS参数化设计语言(ANSYS Paramelric Design Language,简称APDL),用智能分析的手段进行模型的建立、加载、求解和数据后处理。这样建立的APDL命令流文件不仅便于保存和交流,而且有利于多次修改,多次重复分析。1.3本章小结了解了有限元的发展历史,应用领域,发展趋势。简单介绍本课题的研究内容以及分析方法。确定对有约束情况下的零件进行模态分析的实际意义。第二章 振动研究的理论基础2.1 振动基础理论 在机器设备的运转过程中,不可避免地会产生振动,而在振动信号中,含有包括机器运行状态的大量的信息。一旦机器运转异常或者发生故障时,振动信号将会发生一些变化,这些变化具体表现为频率成分、相位差别、幅值大小和能量分布状况的改变等等。振动信号的特征和性质与机器故障、系统固有特性这两大因素均有关联。但即使是同一种故障,发生在机器的不同部位,振动信号所反映出的振动特征和响应可能会有非常大的差异。所以,故障类型和振动信号反映出的振动特征,并不是相互对应的确定关系,这就是基于振动的故障识别的难点所在。所以,研究振动的识别对机器设备故障诊断技术的完善有着极为重要的意义。本质上来说,任何一个振动系统都是一个动力系统,对于大多数情况而言,振动系统所受到的激励和响应都随着时间变化,且系统的响应大都依赖初始条件以及外部激励。绝大部分实际振动系统非常复杂,因此,想要在数学建模分析的时候把所有的细节情况都考虑进来是绝对不可能的。若想预测在确定的激励下振动系统的响应情况,一般需要对振动模型进行简化和抽象,只考虑一些比较重要的因素。大多数实际系统都是连续的,具有无限多个自由度。一些连续系统的振动特性可用偏微分方程描述,求解偏微分方程,是十分困难的。并且,许多偏微分方程并不存在解析解。另一方面,建立多自由度系统的振动方程只要求解一组常微分方程,这相对来说要简单的多。因此,为了分析的简化,经常将连续系统近似为多自由度系统。对于一个物体来说,固有频率是它的重要参数,且固有频率与物体的结构特征有着密不可分的联系,找寻它们之间的关系是研究机械结构振动时常需要面对的一个重要问题。由于机油冷却器的缺陷会引起其固有频率的变化,且固有频率比较容易通过实验的方法测得,故本文通过分析缺陷对机油冷却器固有频率的影响情况,寻找判定机油冷却器产品的好坏的途径和方法。机械系统可视为由质量、刚度和阻尼各元素以一定形式组成而成。略去实际的机械结构的阻尼不计,把它简化为由若干个无弹性的质量和无质量的弹性元件所组成的力学模型,称为弹簧质量系统。对于一个有无限多个自由度的系统,有无限多个固有频率,每个频率对应一个闻有振型,可通过系统的特征方程来确定其固有频率。2.2 连续系统自由振动 当没有外激励作用,即梁弯曲自由振动时,等截面弯曲横向自由振动的运动微分方程中包含四阶微分导数和二阶时间导数 (s1) 设方程(xx)的解对时间和空间是分离的,令y(x,t)=Y(x)F(t) (s2) 将式(s2)带入方程(s1),可得 (s3) (0xL) (s4) 方程(s3)飞通解为简谐函数 (s5) 式中,A和B为积分常数,由两个初始条件确定。通过解方程(s4)可以得到振型函数的一般表达式,这里振型函数Y(x)必须满足相应的边界条件。其中,铰支端的边界条件为Y(x)=0, ,(x=0或x=L) (s6) 令,则方程(S4)可化简为 (s7) 方程(s7)是一个四阶常系数线性微分方程,其解为 (s8) 这就是粱振动的振型函数,其中 为积分常数,可以用四个边界条件来确定其中三个积分常数及导出特征方程,从而确定粱弯曲振动的固有频率和振型函数Y(x)。粱支架的边界条件为Y(x)=0,Y(L)=0, (s9) 将第一组边界条件带入式(s8)及其二阶导函数,得 (s10) 将第二组边界条件带入式(s8)及其二阶导函数,得 (s11) 因为当时,故得,于是,特征方程为 (s12) 它的根为 (s13) (s14) 与此相应固有频率为 (s15) 相应振型函数为 (s16) 因为振型只确定系统中个点振幅的相对值,不能唯一的确定振幅的大小,故其表达式无需再带常数因子,其振型函数为 (s17) 现选用半径为0.05m,长度为1m的粱,材料属性:弹性模量E为21011Pa,密度为7800kg/m3,求其模态。由公式(s16)可知 (s18)由公式(S18)算出悬臂梁前三阶非零固有频率为:一阶为70.92Hz,二阶为451.31Hz,三阶为1243.12Hz。【21】2.3本章小结 简单了解了振动的简化与抽象,求解常微分方程,得到各阶的固有频率求解公式为。且带入数据求解了悬臂梁的前三阶非零固有频率,得到具有实际意义的固有频率参数。第三章利用Ansys软件进行模态分析3.1 ansys有限元软件简介 ANSYS是目前世界顶端的有限元商业应用程序,是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。由世界上最大的有限元分析软件公司之一的美国ANSYS开发,它能与多数CAD软件接口,实现数据的共享和交换,如Pro/Engineer,NASTRAN,IDEAS,AutoCAD等,是现代产品设计中的高级CAD工具之一。美国JohnSwanson博士于1970年创建ANSYS公司后,便开发出了该应用程序,以此用计算机模拟工程结构分析,历经30多年的不断完善和修改,现成为全球最受欢迎的应用程序。 ANSYS是一种广泛的商业套装工程分析软件。所谓工程分析软件,主要是在机械结构系统受到外力负载所出现的反应,例如应力、位移、温度等,根据该反应可知道机械结构系统受到外力负载后的状态,进而判断是否符合设计要求。一般机械结构系统的几何结构相当复杂,受的负载也相当多,理论分析往往无法进行。想要解答,必须先简化结构,采用数值模拟方法分析。由于计算机行业的发展,相应的软件也应运而生,ANSYS软件在工程上应用相当广泛,在机械、电机、土木、电子及航空等领域的使用,都能达到某种程度的可信度,颇获各界好评。使用该软件,能够降低设计成本,缩短设计时间。到20世纪80年代初期,国际上较大型的面向工程的有限元通用软件主要有:ANSYS,NASTRAN,ASKA,ADINA,SAP等。以ANSYS为代表的工程数值模拟软件,是一个多用途的有限元法分析软件,它从1971年的2.0版本与今天的8.0版本已有很大的不同,起初它仅提供结构线性分析和热分析,现在可用来求结构、流体、电力、电磁场及碰撞等问题的解答。它包含了前置处理、解题程序以及后置处理,将有限元分析、计算机图形学和优化技术相结合,已成为现代工程学问题必不可少的有力工具。ANSYS软件是第一个通过ISO9001质量认证的大型分析设计类软件,是美国机械工程师协会(ASME)、美国核安全局(NQA)及近20种专业技术协会认证的标准分析软件。在国内第一个通过了中国压力容器标准化技术委员会认证并在国务院17个部委推广使用。3.2 ansys自由模态分析的基本操作 以21m的平面(材料属性为:弹性模量E为21011Pa,密度为7800kg/m3)为例,分析其自由模态,且讨论了结构单元、网格单元对ansys模态分析结果的影响。结构单元分别为PLANE42(PLANE42被称为2D结构单元,用于模拟平面实体结构。该单元由4个节点定义,每个节点有2个自由度,即沿节点坐标系x和y方向的平动位移,单元模型如图3-1所示)和PLANE82(PLANE82称为2D8节点结构实体单元,是4节点PLANE42的高阶单元,对四边形和三角形的混合网格具有较高的精度,即使是不规则形状,其精度降低也很小,该单元采用协调的位移差值函数,因此能够很好的适应曲线边界。该单元由8个节点定义,每个节点有2个自由度,即沿节点坐标系x和y方向的平动位移,单元模型如图3-2所示);取网格单元大小分别为0.5、0.25、0.125、0.0625(如图3-3所示)。由于要改变结构单元和网格单元,求得不同情况下的自由模态,所以对其进行参数化处理,使用ansys参数化语言APDL比较方便。求得平面自由模态的结果(由于是自由模态,所以前三阶为0)如表3-1所示(使用APDL语言见附录1)。图 3-1 PLANE42单元几何图 3-2 PLANE82单元几何表3-1不同情况下平面的固有频率值单元类型单元长度自由模态各阶频率值45678PLAN420.5889.18128616132417.32431.60.25843.441259.71442.521272136.30.125829.9612531395.92030.42050.80.0625826.441251.3138420052028.4PLAN820.5829.191251.41398.52058.72076.60.25825.51250.71381.22001.12024.60.125825.261250.71380.11996.82021.10.0625825.251250.71380.11996.52020.9 0.5 0.25 0.125 0.0625图3-3不同单元长度的平面有限元模型 由图3可知,随着单元长度的减小,固有频率的值趋于稳定,即此时的模态分析结果更为精确。且PLANE82单元的变化比PLANE42单元小,即PALNE82单元对网格单元大小的要求比PLANE42单元更小。所以进行模态分析时,对结构单元、单元大小的选择十分重要。由表F1、图3-4、图3-5,可知,不同阶数的固有频率值随单元大小的变化也是有所差异的,高阶固有频率对网格单元长度的要求更高,即要使高阶固有频率的值更加精确,则网格单元长度要取的更小。图3-4不同单元4阶固有频率随单元长度的变化图3-5不同单元8阶固有频率随单元长度的变化 由于模型比较简单,计算量不大,各种模态分析方法的比较也没有什么区别。常用方法有Block Lanczos法,PCG Lanczos法,子空间法,缩聚法。其比较见图3-6【23】3-6 针对对称矩阵的特征值求解方法比较3.3约束下梁的模态分析 一般情况下,零件都是具有约束的,所以求解具有约束情况下的模态更具有实际意义。在求解具有约束情况下的模态时,一般将约束当成刚性约束。但是大部分约束都是具有弹性的,例如轴承约束轴。以半径为0.05m,长度为1m的粱为例,材料属性:弹性模量E为21011Pa,密度为7800kg/m3,求其在不同约束刚度下的模态。图3-7约束梁的几何模型 首先选择结构单元beam3粱单元,设置单元属性:半径r为0.05,长度l为1,设置材料属性:弹性模量E为2e11,泊松比为0.3,密度7800。用combin14单元来模拟弹性约束,通过对单元关键选项KEYOPT(2)的选择,设置三种弹簧单元。分别为:X方向轴向拉伸弹簧,Y方向轴向拉伸弹簧,Z方向扭转弹簧。且分别设置单元实常数为K1、K2、K3、K4、K5、K6(对应图3-7所示)图3-8粱有限元模型选择选择分析类型为模态分析Modal,设置扩展模态为3、提取的振型数为3,频率范围默认全部,选择分析方法Block Lanczos Method。通过对K1到K6的不同赋值模拟不同约束的梁,如对K1、K2、K3取值而使K4、K5、K6恒为零,则为悬臂梁。取K1、K2、K3为1020N/m,可将其当作刚性约束,模态分析结果为:1阶72.391Hz,2阶449.79Hz,3阶1242.5Hz,与上章节理论计算结果基本相同,说明对梁的结构单元、单元大小等的选择是符合计算精度的。一般轴承的刚度为109N/m左右,所以对K1、K2、K3的取值范围定为1061012N/m。计算时因为轴有不止一个方向的约束,所以采用控制变量法,即分别对单个方向的约束刚度进行取值进行模态分析,使其他方向的约束刚度不变且取足够大,本文取为1020N/m。其模态如表3-2所示(使用APDL语言见附录2)。 由3-9图可知约束刚度对梁的一阶固有频率影响较大,且随着约束刚度的增大固有频率的值趋于稳定,即刚性约束下的固有频率。则只有当约束刚度达到使固有频率趋于稳定时的值,才能将约束当成刚性约束。且不同方向的约束刚度对一阶固有频率的影响也是有所差异的,其中Y轴向约束刚度至109N/m趋于稳定, 而X轴向约束刚度和Z方向扭转约束刚度至108N/m就已经趋于稳定,即Y轴向约束刚度对梁的一阶固有频率影响最大。由图3-10、3-11可知,随着约束刚度的表3-2悬臂粱的固有频率固有频率阶数123约束刚度lgk1620.30672.391449.79764.15572.391449.79872.391201.05449.79972.391449.79584.031072.391449.791114.71172.391449.791242.51272.391449.791242.5lgk2619.82117.56614.98750.824143.2621.06869.613290.89688.34972.111428.751074.21072.363447.721227.31172.389449.5812411272.391449.771242.3lgk3631.686331.771026.3760.846394.491117.9870.944441.171219.6972.243448.8812401072.377449.71242.21172.39449.781242.51272.391449.791242.5增大,固有频率值的变化有所不同,2阶及以上的固有频率在刚度较小的区域有一段保持不变,从图3-12中可以更加直观的看出,且固有频率不变时的值正好为上阶刚性约束时的值。而各方向约束刚度对二阶、三阶固有频率的影响大致相同,不过随着阶数的升高,X轴向约束刚度对固有频率影响越大。由图3-12可知,单方向约束刚度对不同阶固有频率的影响也不一样,其中X轴向约束刚度至108N/m一阶固有频率趋于稳定;至109N/m二阶固有频率趋于稳定;至1011N/m三阶固有频率趋于稳定,即高阶固有频率对约束刚度的要求更高。 例如选择双列向心短圆柱棍子轴承3182120(d=100mm,D=150mm,b=37mm,i=2,z=30,=11mm,L=11mm,r=0.8mm)。轴承的预紧量g=-5,外圈与箱体的配合过盈量,F=4900N。则支承的刚度为1.707109N/m文献编号 陈忠,即K2=1.707109N/m,取K1、K3为1020N/m,则悬臂梁前三阶固有频率:一阶为72.227Hz、二阶为437.53Hz、三阶为1145.4Hz。若将其当成刚性约束则误差分别为:一阶0.227%;二阶为2.73%;三阶为7.81%。图 3-9不同约束刚度下悬臂梁的一阶固有频率图 3-10不同约束刚度下悬臂梁的二阶固有频率图 3-11不同约束刚度下悬臂梁的三阶固有频率图 3-12不同X轴向约束刚度下悬臂梁的固有频率值 对K1、K2、K5赋值,其余刚度为零,则为简支梁。对简支梁进行模态分析,记录不同约束刚度下前三阶固有频率,如表3-3所示(使用APDL语言见附录3)。由图(3-13、3-14、3-15)可知,各阶固有频率随不同约束刚度的变化,简支梁与悬臂梁相同。一阶固有频率随约束刚度的增大而变大,且趋向于稳定值,二阶及以上固有频率在约束刚度较低时有一段保持不变。对于不同方向刚度,Y轴向对固有频率的影响较大,且随着阶数的增多,X轴向刚度对固有频率的影响随之增大。表3-3简支粱的固有频率固有频率(Hz)阶数123约束刚度(N/m)lgk1620.306202.88804.14764.155202.88804.148201.05202.88804.149202.88584.03804.1410202.88804.141114.711202.88804.141274.712202.88804.141293.5lgk2628.4949.503457.86784.065154.94489.298169.57441.88744.779198.88743.051295.710202.47797.911295.711202.84803.521295.712202.87804.081295.7图3-13 不同X轴向约束刚度下简支梁各阶固有频率图3-14 不同约束刚度下简支梁的一阶固有频率图3-15不同约束刚度下简支梁的二阶固有频率3.4本章小结了解有限元软件ANSYS的诞生背景以及应用领域。对同一平面取不同的结构单元和网格单元对比其对各阶固有频率的影响。随着单元长度的减小,固有频率的值趋于稳定,即此时的模态分析结果更为精确,高阶固有频率对网格单元长度的要求更高。且使用高阶单元结构令分析结果更为精确。计算了不同约束梁的固有频率随约束刚度的变化。约束刚度对梁的一阶固有频率影响较大,且随着约束刚度的增大固有频率的值趋于稳定,即刚性约束下的固有频率。其中Y轴向约束刚度对梁的一阶固有频率影响最大。且随着阶数的增多,X轴向刚度对固有频率的影响随之增大。 第4章 利用ansys进行悬挂法分析仿真4.1无约束下箱体自由模态分析 实际测量零件的模态时,通常选择悬挂法测量零件的自由模态,即利用绳索将零件悬挂起来,使零件自由悬浮在空中。然后用力锤敲击零件,通过传感器测试零件上各点的振动。下面选择箱体为测试零件(长宽高均为1m,壁厚为0.1m的箱体),由上章节可知,选择高阶单元有利于得到更加精确的结果,且对网格大小的要求较低,所以此处选用SOLID95单元(SOLID95称为3D20节点结构实体单元,是SOLID45的高阶单元,对不规则形状也具有较好的精度;由于采用协调的位移差值函数,可很好的适应曲线边界。该单元由20个节点定义,每个节点有3个自由度,即沿节点坐标系x、y和z方向的平动位移,单元模型如图4-16所示)。材料属性:弹性模量E为2e11,泊松比为0.3,密度7800。为了更好的控制网格单元大小,选择6面体单元对零件进行网格划分,如图4-17所示。首先进行自由模态分析选择,选择分析类型为模态分析Modal,设置扩展模态为10、提取的振型数为10,频率范围默认全部,选择分析方法Block Lanczos Method。结果如表4-4所示(由于是自由模态,前6阶为零),振型如图4-18所示表 4-4箱体自由模态各阶固有频率阶数5678910固有频率(Hz)4.54E-045.72E-04195.87349.89591.62591.62图4-16 SOLID95单元几何图 4-17箱体有限元模型第7阶第8阶第9阶第10阶图 4-18箱体自由模态各阶振型4.2悬挂法测试箱体固有频率 采用悬挂法测出的模态是否为自由模态,且悬挂绳索的刚度对模态分析是否具有影响。利用弹簧单元COMBIN14单元表示悬挂的绳索,如图4-19所示。且通过对单元关键选项KEYOPT(2)的选择,使其为一维纵向弹簧单元,自由度UY。通过对单元实常数K来控制弹簧刚度,不同弹簧刚度(106N/m、109N/m、1012N/m)对应的各阶固有频率如表F5所示(使用APDL语言见附录4)图 4-19悬挂法有限元模型表 4-5不同弹簧刚度对应的各阶固有频率固有频率(Hz)阶数5678910lgk(N/m)65.72E-042.7669195.87349.89591.62591.6295.80E-0472.273195.87349.89591.62591.62125.97E-04127.44195.87349.89591.62591.62 由表4-5可知不同弹簧刚度对应的各阶固有频率基本相同,除第6阶固有频率。所以利用悬挂法测出的固有频率即为自由模态频率。第6阶固有频率的不同是因为利用弹簧单元来模拟绳索具有系统误差,因为绳索只有单向刚度,即Y负方向,而弹簧单元具有双向刚度,限制了箱体上节点的Y位移,振型如图D20所示。其它振型如图D21所示图 4-20 悬挂法模态分析第6阶振型第7阶第8阶第9阶第10阶图 4-21 悬挂法模态分析各阶振型 上述为单根绳索悬挂情况,但有时零件利用单根绳索不易悬挂,进而要使用2根或多根绳索悬挂,由于固有频率与弹簧刚度变化的基本无关,所以此处弹簧刚度K选定为1109N/m。2根绳索时的有限元模型如图4-22所示,各阶固有频率随弹簧刚度的变化如表4-6所示。其固有频率与自由模态相比略有误差,除了第5、6阶模态,此处与单根绳索相同,还是由于绳索只有单向刚度,即Y负方向,而弹簧单元具有双向刚度,且2根绳索的影响更大。图 4-22悬挂法(2根绳索)有限元模型表 4-6根绳索时 各阶固有频率随弹簧刚度的变化阶数5678910固有频率(Hz)103.89132.78208.31349.89591.62593.884.3本章小结 利用ansys进行悬挂法分析仿真,利用弹簧单元COMBIN14单元表示悬挂的绳索,弹簧单元选取不同刚度时箱体的固有频率值,并且与自由模态分析进行对比。得到结论为悬挂绳索的刚度对模态分析没有影响,且固有频率值恒为自由模态时的值。最后使用悬挂法(两根)进一步验证结果,得到结论相同。所以利用悬挂法测出的固有频率即为自由模态频率。总结 根据梁振动的微分方程求解并得到了悬臂梁各阶固有频率。简单了解了振动的简化与抽象,求解常微分方程,得到各阶的固有频率求解公式。且带入数据求解了悬臂梁的前三阶非零固有频率。 本课题选用了平面结构作为研究的模型,利用ANSYS进行建模、确定边界条件,施加载荷,划分单元网格,进行自由模态分析模态分析。且改变单元结构、网格大小等因素,讨论其对模态分析的影响。随着单元长度的减小,固有频率的值趋于稳定,即此时的模态分析结果更为精确,高阶固有频率对网格单元长度的要求更高。且使用高阶单元结构令分析结果更为精确。计算了不同约束梁的固有频率随约束刚度的变化。约束刚度对梁的一阶固有频率影响较大,且随着约束刚度的增大固有频率的值趋于稳定,即刚性约束下的固有频率。其中Y轴向约束刚度对梁的一阶固有频率影响最大。且随着阶数的增多,X轴向刚度对固有频率的影响随之增大。 利用ansys进行悬挂法分析仿真,利用弹簧单元COMBIN14单元表示悬挂的绳索,弹簧单元选取不同刚度时箱体的固有频率值,并且与自由模态分析进行对比。得到结论为悬挂绳索的刚度对模态分析没有影响,且固有频率值恒为自由模态时的值。最后使用悬挂法(两根)进一步验证结果,得到结论相同。所以利用悬挂法测出的固有频率即为自由模态频率。致谢首先衷心地感谢我的导师吴志学教授,感谢恩师在学业上的悉心指导。恩师朴实的为人、渊博的学识、务实的作风以及对事业的执着,总是让学生仰慕不已。在导师的精心栽培下,学生在文献收集与整理、新知识的吸收与消化、创新思维方法等方面均获得了相当的研究经验。再次向吴老师致以最崇高的敬意和最诚挚的感谢!特别感谢父母及亲人,感谢他们多年来的支持、关心和理解,使作者的大学学业得以顺利完成。向所有关心和帮助作者的各位师长、同学和朋友们致谢! 参考文献1 ANSYS高级技术指南2 蒋玉川,李章政, 弹性力学与有限元方法 , 化学工业出版社 ,2010;3 夏建芳,叶南海, 有限元法原理与ANSYS应用;4 任学平,高耀东, 弹性力学基础及有限单元法, 华中科技大学出版社, 2006;5 孙泽世,王锋,子模型技术在桥梁分析上的应用;6 宋志安,于涛,李红艳,等,机械结构有限元分析ANSYS与ANSYS Workbench 工程应用,国防工业出版社,2010;7 廖日东,有限元法原理简明教程,北京理工大学出版社,2009;8 鲁建霞,苟惠芳,有限元法的基本思想与发展过程;9 宁连旺,ANSYS有限元分析理论与发展;10 李黎明,ANSYS有限元分析实用教程 ,清华大学出版社 ,2005;11 张红松,胡仁喜,康士廷等,ANSYS 13.0 有限元分析;12 许文本,焦群英,机械振动与模态分析基础. 北京:机械工业出版社. 1998;13 师汉民,机械振动系统:分析测试建模对策(上、下册)(第二版). 华中科技大学出版社. 200414 刘国庆,杨庆东 ANSYS工程应用教程机械篇 2003 15 杨帅,王太勇.,竖直方向弹性约束悬臂梁的固有频率分析 天津大学报 2011(44): 18-2216 陈忠,滚动轴承及其支承刚度计算 煤矿机械 2006(3)387-388;17 杨康,韩涛,ANSYS在模态分析中的应用 佳木斯大学学报 2005;18 袁安富,陈俊 ANSYS在模态分析中的应用;19 梁军,赵登峰 模态分析方法综述;20 龙英,滕召金,赵福水 有限元模态分析与发展趋势;21 李东旭,高等结构部动力学(P319-P329) 科学出版社 2010;22 龚曙光,谢桂兰,黄云清,ANSYS参数化编程 机械工业出版社2009;23 曾攀,雷丽萍,方刚 ,结构的建模与分析 2010;24 王新敏,李义强,许宏伟,ANSYS结构分析单元与应用 2011附录1/PREP7 A=0.5 ET,1,PLANE42 MPTEMP,1,0 MPDATA,EX,1,2e11 MPDATA,PRXY,1,0.3 MPTEMP,1,0 MPDATA,DENS,1,7800 RECTNG,-1,1,0,1, ESIZE,A,0,MSHAPE,0,2D MSHKEY,0 AMESH,ALL FINISH ANTYPE,2 MODOPT,LANB,6 EQSLV,SPAR MXPAND,6, , ,0 LUMPM,0 PSTRES,0MODOPT,LANB,6,0,0, ,OFF /STATUS,SOLUSOLVE 附录2/PREP7NUM=30PI=3.1415926R=0.05A=PI*R*RI=0.25*PI*R*4L=1ET,1,BEAM3R,1,A,IMP,EX,1,2E11MP,PRXY,1,0.3 MP,DENS,1,7800SECTYPE,1,BEAM,CSOLID,0SECOFFSET, CENTSECDATA,R N,1,0,0N,NUM,L,0FILL,1,NUM,NUM-2TYPE,1 MAT,1REAL,1*DO,N,1,NUM-1E,N,N+1*ENDDO/SOLANTYPE,MODALMODOPT,LANB,6MXPAND,6SOLVE附录3/PREP7 ET

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论