基于能量等效的直齿锥齿轮时变啮合刚度计算

针对直齿圆锥齿轮啮合刚度的计算问题,从微元思想出发,将变截面齿廓划分为若干等截面微段齿段,基于能量等效建立微段齿段啮合刚度计算模型,并利用积分方法获得单齿啮合刚度。此外,基于力平衡和变形协调条件进一步提出了齿轮时变啮合刚度计算模型,同时根据几何关系导出了相应传动误差计算式。利用有限元分析对解析计算模型进行了验证,并分析了误差来源。结果表明:利用该模型不仅能够将直齿锥齿轮啮合刚度计算精度保证在2%以内,还达到了快速求解的目的。      

基于预定啮合特性的点啮合齿面设计方法

针对格里森螺旋锥齿轮齿面设计方法无法在整个齿面接触传动过程中有效控制齿面啮合特性的不足,论述一种按预定的啮合特性设计点啮合齿面的理论和方法:啮合齿面的接触迹线上每一点的几何结构按预定的齿面啮合特性要求设计,而与特定的机床结构参数无关.啮合点的二阶接触参数通过弹性齿轮副的载荷-变形效应条件,而不是接触区的位置、大小和形状来设计.实例计算和分析表明本文的点啮合齿面设计方法是格里森齿面设计方法在数控技术条件下的发展和完善.      

一种考虑齿轮副连续啮合过程的接触有限元动力学分析方法

在系统分析齿轮副连续啮合过程不同典型啮合状态特点的基础上,提出了一种可以考虑齿轮副连续啮合过程中啮合齿对变化,受载弹性变形以及滑动摩擦等影响因素的接触有限元分析方法.利用该方法对三种考虑不同因素的齿轮副模型进行了对比分析,并得到了各种因素对齿轮副连续啮合过程动态传递误差及动态接触力特性的影响.研究表明:该方法可以真实模拟齿轮副连续啮合过程中单、双齿对啮合及其交替啮合状态的动力学特性,包括滑动摩擦引起的节点冲击,受载弹性变形引起的啮入、啮出冲击以及时变刚度等激励特性;并可以得到啮入、啮出冲击的大小及作用时间,以及滑动摩擦和齿廓修形对动态啮合特性的影响;动态啮合特性分析结果与齿轮副连续啮合状态的特征完全吻合.      

基于数值的弧齿锥齿轮齿面啮合点位置分析的拟赫兹法

在已知弧齿锥齿轮齿面网格的条件下,分析了单齿啮合时弧齿锥齿轮齿面弹性变形对轮齿啮合点位置的影响。为此,首先形成一整套刚性齿面啮合点的数值计算方法;然后,采用赫兹接触理论计算齿面弹性变形,确定齿轮轮齿的微小转动及由此引起的啮合点位置的变动。将齿轮因齿面变形而产生微小转角,继而进行齿面啮合分析的过程定义为拟赫兹接触分析。结果表明,齿面弹性变形引起的齿轮轮齿的微小转动对啮合点的最终位置有一定的影响。      

螺旋锥齿轮啮合刚度及参数振动稳定性研究

准确计算时变啮合刚度是齿轮系统动力学研究的基础.针对航空高速重载螺旋锥齿轮,基于轮齿接触分析(TCA)和轮齿加载接触分析(LTCA)通过计算瞬时接触点的轮齿变形柔度建立了时变啮合刚度数值模型;将齿轮时变啮合刚度在一个啮合周期内视为逐段线性,基于Floquet理论推导了含时变刚度参数振动系统的状态转换矩阵解析式;通过修正小轮机床调整参数设计三种接触情况,分析了算例齿轮在相同载荷工况下的接触轨迹、传动误差、重合度和时变啮合刚度;采用二自由度齿轮系统动力学模型考察工作转速范围内的周期运动不稳定区间,分析了时变啮合刚度对螺旋锥齿轮系统参数振动稳定性的影响.      

差曲面拓扑的齿轮啮合刚度计算与承载接触分析

基于齿条-齿轮等切共轭产形原理,构建齿面数值模型、ease-off差齿面,对ease-off蕴含的齿面啮合信息进行解析,获得了齿面接触路径、传动误差、接触线瀑布图;综合ease-off拓扑仿真与轮齿刚度非线性单元耦合解析,给出了修形拓扑齿面的啮合刚度、承载传动误差的计算方法。沿接触路径遍历接触线序列,获得了轮齿时变啮合刚度、承载传动误差与载荷分布图;给出了2阶抛物面对称与对角拓扑两种修形形式算例,求出了系列载荷作用下的啮合刚度、承载传动误差、齿面载荷分布。结果显示:随着载荷的增加,轮齿啮合刚度时变效应明显减弱;承载传动误差波动与啮合刚度、修形梯度密切相关;对角修形在啮合刚度、传动误差、载荷分布特性方面好于对称修形。      

弧线齿面齿轮应力过程分析

为了给弧线齿面齿轮的齿面接触强度和齿根弯曲强度设计提供理论依据,研究了弧线齿面齿轮的齿面接触应力和齿根弯曲应力随载荷和安装误差的变化规律.在齿面接触分析和承载接触分析的基础上应用弹性理论计算了弧线齿面齿轮副的齿面接触应力和应用有限元应力影响矩阵法计算了该齿轮副的齿根弯曲应力.给出了数字计算实例,计算结果表明:齿面接触强度和齿根弯曲强度在重载时的接触强度和弯曲强度由单齿啮合区的强度决定,轴向安装误差和轴夹角安装误差分别会增加齿面接触应力和齿根弯曲应力,轴夹角安装误差和轴间距安装误差对齿面接触应力影响甚小,而轴向安装误差和轴间距安装误差可以降低齿根弯曲应力,与直齿面齿轮相比,弧线齿面齿轮的接触和弯曲应力明显减小.      

构造拓扑修形齿面的面齿轮传动主动设计

为了改善加工参数较少的面齿轮传动的啮合性能,提出了用给定的啮合性能对面齿轮和小齿轮齿面进行拓扑修形设计的方法.面齿轮用插齿法加工,面齿轮与插齿刀的转角关系由预设的传动误差确定,由此确定面齿轮的拓扑修形齿面.然后用接触路径的位置及其方向和接触椭圆的大小构造小齿轮的拓扑修形齿面,该齿面用共轭点接触法磨齿加工,建立了小齿轮拓扑修形齿面与加工参数之间的线性方程.结果表明:齿面接触分析(TCA)获得的传动误差、接触路径与预设的传动误差、接触路径相同,TCA仿真的椭圆长度与预设椭圆长度相差范围为0.175 2~1.16mm.     

齿轮减速器系统可变固有特性动力学研究

考虑到齿轮传动啮合刚度的波动和传动误差的影响以及轴承支撑刚度的作用,对二级齿轮减速器传动系统进行了理论建模和动态响应分析,并与实验结果进行了比较。结果表明,齿轮传动在单齿啮合区和双齿啮合区之间啮合刚度变化较大;减速器系统的动态特性 (固有频率、固有振型、阻尼等 )随啮合周期而发生变化,呈现出一种可变的动态固有特性。故对于系统进行研究时,可分别按单齿区和双齿区平均啮合刚度进行分析,一般可以满足实际工程要求。     

行星轮系动态啮合应力研究

针对行星传动齿轮系的结构特点,建立了行星齿轮和主动齿轮的非线性接触有限元模型,应用AN-SYS软件进行了动态接触分析的研究。基于研究结果,建立了可以体现行星齿轮系在工作过程的啮合应力变化情况的计算模型,可以直观地观察到行星齿轮系从开始啮合到完全脱开整个过程中的应力分布及变化情况。选择20节点六面体单元,模拟齿轮受载情况,迅速、高效地得出整个齿轮系的强度、刚度分布云图,大大提高了设计效率。     

多激励下某直升机传动系统动载特性

针对由弧齿锥齿轮和行星轮系构成的直升机传动系统,构建了纯扭振动模型,采用集中参数法建立了齿侧间隙非线性动力学方程.通过有限元方法求得了时变啮合刚度,采用4-5阶变步长Runge-Kutta法对动力学方程进行了数值求解,借助动载系数、相图、Poincaré截面图、快速傅里叶变换频谱图等分析手段,研究了传动系统在时变啮合刚度、齿侧间隙、综合传动误差、外载荷等多种激励作用下系统的动载特性.结果表明啮合刚度对传动系统的影响最大,动载系数最大值为1.5;齿侧间隙对系统响应特性的影响是有限的;啮合误差在一定程度上抑制了齿轮系统的振动;外载荷波动对不同速级的影响不同,动载系数最大值发生在并车传动.     

准双曲面齿轮LTCA弯曲变形约束方法

针对现有准双曲面齿轮加载接触分析（LTCA）受载后得到的弯曲变形中包含由于柔性引起的变形和加载点与轮坯之间剪切变形的问题，提出了轮齿有限元模型约束方法用于消除柔性引起的变形和剪切变形。该约束方法在原有约束方法的基础上引入了约束轮坯内圈和轮齿中间平面的有限元模型，并给出了该约束方法对应的齿面法向柔度矩阵的计算过程。为验证该约束方法的正确性，将原约束方法、该约束方法及实验结果进行对比。结果表明:约束方法主要影响齿面加载接触印迹的范围，最大载荷时原有约束方法加载接触印迹到达齿面小端，而该约束方法的加载接触印迹在最大载荷时距离小端仍有3~4mm与实验结果吻合更好。      

计及套圈变形的薄壁角接触球轴承性能分析

为了准确计算柔性支承轴承负荷分布,以外圈孤立点支承薄壁角接触球轴承为例,建立了三维有限元分析模型,考虑外圈受载后弯曲变形对滚动体负荷分布的影响,分析不同工况下轴承负荷分布,寿命及刚度等性能的变化规律,并将有限元法的计算结果与拟动力学法及实验结果进行分析比较.结果表明,有限元方法计算结果与实验结果更为接近.      

弹性环刚度特性研究

本文针对两种弹性环刚度的计算方法——解析法和有限元法进行了对比,分析了各自的适用条件,并基于有限元法分析了涡动角、凸台参数和配合间隙对弹性环刚度的影响,为弹性环的结构设计提供了参考。     

高效非结构网格变形与流场插值方法

非结构动网格在含动边界的非定常流动模拟中应用十分广泛。为了提高网格变形的计算效率和对复杂外形的适用性,发展了高效的K近邻-径向基函数（KNN-RBF）动网格方法。并采用高效的交界面节点一致性方法实现了网格变形的并行化。由于边界的运动和变形形式多种多样,在自主开发的计算流体力学（CFD）程序中实现了动网格的用户自定义函数（UDF）功能模块,可以灵活地设置此类问题中的边界运动。在模拟大变形问题时,采用了全局网格重构的策略对网格质量较差的区域进行处理。提出了快速查找新旧网格最近点的两级KD（K-Dimensional）树方法,并根据最近点在新旧网格间进行流场插值。从而建立了高效健壮的并行流场插值方法。通过对若干算例进行模拟,验证了所发展方法的正确性和健壮性。     

阻尼器切向接触刚度的有限元分析及应用

提出了一种利用有限单元法求解干摩擦阻尼器接触刚度的方法.利用该方法研究了接触面初始正压力、摩擦系数、阻尼器材料属性和阻尼器接触面半径等对接触刚度的影响规律,发现接触刚度随着这些参数的增大而变大.并通过曲线拟合得到了一定条件下接触刚度随接触面正压力变化的函数关系,为求解带阻尼器叶片非线性响应时实时准确考虑接触刚度随正压力的变化提供了方便.最后,利用所得刚度函数研究了接触刚度对带凸肩叶片非线性振动响应的影响,结果表明:随着接触刚度的增加,叶片峰值频率逐渐升高,振动峰值逐渐减小,且各自都逐渐趋于一个定值.      

轮盘端面齿连接的周向拉杆转子刚度研究

针对轮盘端面齿连接的周向拉杆转子模型,理论分析得到转子各组成部分刚度计算方法.应用商业有限元软件ANSYS,基于三维有限元非线性接触算法计算预紧力下轮盘端面齿连接段的弯曲刚度和扭转刚度,得到刚度修正系数和刚度无量纲系数.分析周向拉杆刚度和轮盘刚度的关系,提出拉杆刚度比例系数.计算结果表明:刚度系数无量纲系数小于1时,刚度修正系数保持不变;刚度无量纲系数大于1时,此时预紧力不足,齿面发生滑移,刚度修正系数会突然降低,拉杆弯曲刚度比例系数会增加.刚度无量纲系数具有明确的物理意义,可为带端面齿的燃气轮机转子设计和故障诊断提供参考.      

非结构动网格生成的弹簧-插值联合方法

 提出了基于弹簧法的两点改进,以解决边界发生大变形时的非结构网格变形问题。为了使边界运动引起的网格变形能更好地由边界传递到内部网格中,提出了一种新的基于Delaunay网格插值的弹簧倔强系数逐层改进方法。为了提高弹簧法的计算效率,引入背景网格和直接插值方法,提出了弹簧-插值法。弹簧-插值法首先生成计算域的背景网格(粗网格),然后由弹簧法求解边界运动引起的背景网格变形,最后利用变形后的背景网格直接插值得到计算网格的变形。算例结果表明：改进后的方法一方面有效地提高了动网格的变形能力和变形后的网格质量;另一方面通过降低弹簧法的求解规模,显著地提高了动网格的变形效率。