0 引言膜片联轴器是一种高性能金属元件挠性联轴器，由法兰、螺栓、膜片组、衬套和自锁螺母等构成。相较于非金属弹性元件联轴器和刚性联轴器，膜片联轴器可通过金属膜片的变形来补偿传动系统的轴向、角向、径向等不对中安装误差，且具有性能稳定、传递转矩大、使用寿命长以及可靠性高等优点[1]，广泛使用在各类汽轮机组、风力发电机组以及航空传动系统等传动系统中[2-3]。膜片形状主要有圆环型、轮幅型、连杆型、多边型和束腰型。其中，束腰型膜片由于具有更好柔性、去除多余质量、接近等强度、应力水平低等特点，应用更为广泛。膜片联轴器的动力学参数以及安装情况直接影响轴系的动力学特性，进而影响轴系运转过程中的稳定性。因此，研究膜片联轴器的刚度与阻尼特性是设计联轴器结构及分析含膜片联轴器轴系动力学特性的基础。国内外学者对膜片联轴器已经进行了内容丰富的研究。刘子良等[4]推导了四螺栓膜片联轴器的横向刚度，认为双膜片组联轴器刚度呈周期性波动。汪中厚等[5]根据联轴器刚度计算模型，推导了橡胶块、膜片部件对联轴器径向、轴向刚度影响程度的计算公式。刘景明等[6]利用有限元仿真软件分析了结构参数对圆环式膜片联轴器扭转刚度的影响，分析出关键的结构参数。范祥等[7]研究了膜片受载张口对膜片联轴器周向刚度的影响，建立了膜片弓起时扭转刚度的力学模型。Redmond[8]研究了不对中对膜片联轴器动力学特性的影响，认为不对中问题会导致振动问题，需要尽可能减小初始不对中量。Stocco等[9]设计了一种能够承受因轴向力产生的线性及非线性变形的柔性膜片联轴器。Dobre等[10]认为弹性联轴器对补偿轴向偏移有一定作用，轴向和角向位移载荷会对结构强度产生相应的影响。目前，膜片联轴器在国内大型机械中已经得到了广泛应用，但其生产制造主要依赖于经验，对其性能及关键参数的研究还不够充分。膜片结构参数对联轴器刚度的影响规律以及刚度、阻尼的建模还鲜见于文献，有较大的研究需求和研究空间。本文利用有限元分析软件，对膜片联轴器刚度、阻尼进行了仿真分析，并进行了试验验证。1 膜片联轴器静刚度仿真计算1.1　有限元模型建立膜片联轴器主要由膜片组、紧固件及输入、输出法兰构成。本文基于Ansys有限元仿真软件，建立膜片联轴器整体仿真模型。为得到较好的网格质量、提高仿真效率，简化了螺栓、螺母的建模，将螺栓、螺母部件六角头简化为圆柱头。采用多区网格划分方法划分法兰座、膜片、螺栓以及凸缘衬套为六面体单元，并通过面网格细化部件间的接触圆柱表面；采用扫掠网格划分法划分螺母、球面垫圈和衬套等标准圆柱形特征部件为六面体单元；为较好仿真模拟膜片的弯曲变形，膜片单元中每个膜片使用两层单元。为了获得更好的收敛性，对膜片间的接触对进行了相同尺寸的网格划分。有限元建模、网格划分及膜片结构参数示意分别如图1（a）、图1（b）所示。膜片联轴器基本参数如表1所示。10.16578/j.issn.1004.2539.2024.04.001.F001图1膜片联轴器结构、网格划分Fig. 1Diaphragm coupling structure and mesh division10.16578/j.issn.1004.2539.2024.04.001.T001表1联轴器基本参数Tab. 1Basic parameters of the diaphragm coupling膜片材料法兰材料螺栓材料螺母材料SUS30142CrMo30Ni4CrMoA30CrMoA其他材料分布圆直径D/mm膜片厚度h/mm膜片数量Z30CrMnSiA1050.4513螺栓孔数n螺栓孔径Dk/mm最小宽度Wb/mm膜片外径Dm/mm61012130膜片联轴器接触设置如图2所示。螺母与衬套表面因预紧力作用而保持相对位置不变；螺母与螺栓建模不考虑螺纹影响，且不发生相互滑移；垫片与凸缘衬套实际装配关系为过盈配合。因此，将上述接触表面间的接触关系设为绑定接触。其余部件间的接触大多为圆柱面接触，需考虑接触变形问题，无法简化为线性接触，因此，设定为有摩擦或无摩擦非线性接触。金属间界面摩擦参考文献[11]，无润滑不锈钢之间的摩擦因数为0.15~0.2，有润滑不锈钢材质之间摩擦因数为0.1~0.12。膜片采用表面覆盖有低摩擦因数的涂层来减小接触摩擦力，膜片间的摩擦因数取0.12，其他接触表面之间的摩擦因数取0.15。膜片承受载荷作用主要以弯曲变形为主，初始设置膜片接触面之间的接触刚度系数为0.1，进一步调节膜片间节点穿透系数、网格质量等，使膜片应力、支反力等结果达到收敛，则此时相关参数设置较为合理。10.16578/j.issn.1004.2539.2024.04.001.F002图2膜片联轴器接触设置Fig. 2Contact setting of the diaphragm coupling1.2　静刚度仿真结果膜片联轴器各向刚度仿真的约束及边界条件加载点如图3所示。全局坐标位于膜片组垂直于z轴方向上的中心面处。设立远程控制点1和控制点2，两点均与全局坐标点重合。从动法兰内圆柱面空间自由度与控制点1相耦合，主动法兰圆柱面与控制点2相耦合，且圆柱面设定为刚体，忽略其变形影响。10.16578/j.issn.1004.2539.2024.04.001.F003图3膜片联轴器静刚度仿真约束及边界条件Fig. 3Simulation constraints and boundary conditions of the static stiffness of the diaphragm coupling扭转、轴向、角向及径向刚度仿真时，对远程控制点1施加点约束，约束其空间6自由度。对远程控制点2分别施加转矩、轴向位移、角向弯曲角度以及径向位移载荷，约束其他方向自由度。施加外载荷前，对螺栓施加给定预紧力。由于膜片联轴器两侧交错安装，旋转过程中膜片段所在位置发生变化。联轴器角向偏角载荷产生的支反力矩在旋转过程中是波动的。图4给出了角向刚度在旋转周期内取得最大值以及最小值的两个特殊位置，Rot x和Rot y分别表示膜片联轴器绕x轴偏转和绕y轴偏转。同时，径向刚度在1个旋转周期内也具有交变特性，考虑大变形效应，径向刚度值在y、-y轴位移方向也存在偏差。10.16578/j.issn.1004.2539.2024.04.001.F004图4膜片联轴器角向偏转特殊位置Fig. 4Special positions of the angular deflection of the diaphragm coupling（a）绕x轴偏转 （b）绕y轴偏转为尽可能逼近真实静刚度值，用每个载荷间隔内的平均刚度值等效替代当前对应总累计载荷下的静刚度值。其中，转矩增量为100 N·m、轴向位移增量为0.2 mm、角向偏角增量为0.05°、径向位移增量为0.02 mm。得到各向静刚度仿真结果如图5所示。10.16578/j.issn.1004.2539.2024.04.001.F005图5膜片联轴器静刚度仿真结果Fig. 5Simulation results of the static stiffness of the diaphragm coupling根据图5得出结论：①转矩较小时，螺栓预紧力提供的摩擦力矩平衡转矩，转角较小，膜片联轴器的扭转刚度随着传递转矩呈非线性增长；随着传递转矩的增大，膜片和螺栓同时发生变形，以抵抗转矩作用，联轴器扭转静刚度呈非线性下降趋势；之后由于传递转矩过大，受压段膜片发生屈曲失稳，导致膜片联轴器扭转刚度快速下降直至平稳。②膜片联轴器的轴向刚度随轴向位移的增大呈“指数型”增长趋势。轴向位移较小时，联轴器受拉刚度大于受压刚度；位移逐渐增大时，联轴器受压刚度大于受拉刚度，且两者差距逐渐增大。③膜片联轴器的角向刚度随角向不对量的增大呈非线性先减后增趋势。两个特殊位置的膜片联轴器角向刚度表现为：绕图4中y轴的刚度值较大，绕-y轴的刚度较小。④膜片联轴器的径向刚度随径向位移增大呈非线性下降趋势。随着径向位移增大，受压段膜片易出现屈曲变形，联轴器径向刚度快速下降。其中，沿y轴正方向偏移时，膜片最容易发生屈曲失稳，导致径向刚度下降最快。1.3　膜片结构参数对静刚度影响基于上述有限元模型，分别调整膜片最小宽度（10~16 mm）、螺栓孔分布圆直径（90~110 mm）、膜片厚度（0.35~0.55 mm）以及膜片螺栓孔径（9~13 mm）这4个膜片主要结构尺寸参数，分析膜片结构尺寸参数对膜片各向静刚度的影响。考虑到膜片不同结构尺寸值差距较大，为更好地显示不同结构参数变化对联轴器单向静刚度的影响，采用平均刚度表征联轴器静刚度，并对膜片结构参数进行归一化处理，得到膜片结构参数对膜片静刚度的影响，如图6所示。10.16578/j.issn.1004.2539.2024.04.001.F006图6膜片结构参数对静刚度的影响Fig. 6Influence of diaphragm structural parameters on the static stiffness根据图6得出结论：①增大膜片螺栓孔分布圆直径，膜片联轴器扭转刚度提高，而轴向、角向和径向刚度均降低。这表明增大膜片螺栓孔分布圆直径可提高联轴器承扭能力，同时提高其补偿能力。②增大膜片厚度，膜片扭转和径向刚度先增后减。因为膜片联轴器整体刚度还需考虑螺栓等部件变形，膜片厚度过大，考虑大变形效应，螺栓剪切变形严重，所以联轴器整体刚度反而下降。而轴向以及角向刚度随着膜片厚度增大而非线性增大，表明增大膜片厚度，联轴器不对中补偿能力降低。③增大膜片最小宽度值，膜片联轴器扭转、角向以及径向静刚度均增大，轴向刚度呈先增后减趋势，但整体变化较小。④增大膜片螺栓孔径，联轴器扭转、角向以及径向刚度均增大。而联轴器轴向刚度随膜片螺栓孔径的增大而减小，整体变化也较小。总体上看，改变膜片分布圆直径以及膜片厚度对联轴器静刚度影响较大。2 动态刚度及等效阻尼系数仿真计算2.1　扭转及轴向瞬态动力学仿真结果将转矩和轴向位移载荷条件分布设置为M(t)=MAsin(ωt)和D(t)=DAsin(ωt)的正弦激励力，得到不同幅值、频率下转矩M(t)与转角φ(t)、轴向支反力F(t)与轴向位移D(t)的迟滞回线，分别如图7、图8所示。10.16578/j.issn.1004.2539.2024.04.001.F007图7不同激励幅值下迟滞回线Fig. 7Hysteresis loops under different excitation amplitudes10.16578/j.issn.1004.2539.2024.04.001.F008图8不同激励频率下迟滞回线Fig. 8Hysteresis loop under different excitation frequencies由图7可以看出，膜片联轴器在扭转正弦振动下，迟滞回线形状类似斜椭圆形态。在未发生较大屈曲变形情况下，联轴器扭转刚度变化较小，曲线形状变化较小；在轴向正弦振动下，膜片联轴器轴向刚度随轴向位移呈指数型变化趋势，迟滞回线形状由近似斜椭圆形态转变为S形态，且越靠近轴向位移零点，阻尼耗能越小。由图8可以看出，在扭转及轴向正弦激励加载下，改变膜片联轴器激励频率，在单循环周期内的摩擦耗能基本不变，迟滞回线力或位移响应极值点变化较小。这说明膜片联轴器扭转和轴向振动1周期内，结构阻尼耗能与激励频率无关；而小激励频率情况下，结构动刚度变化较小。2.2　动态刚度以及等效阻尼系数计算首先，基于迟滞回线求解等效动态刚度及阻尼系数的过程，通过函数拟合得到振动周期内的非线性弹性恢复力函数以及迟滞非线性阻尼力函数；其次，将两者叠加，得到1个周期的非线性迟滞恢复力函数，基于拟合优度R2值判定拟合程度；最后，基于拟合得到振动周期内的非线性弹性恢复力函数以及迟滞非线性阻尼力函数，进一步计算等效动刚度及阻尼系数。2.2.1 非线性迟滞恢复力与迟滞回线非线性迟滞恢复力是位移x、速度x˙、振幅A和频率f的函数，即[12]2996[13]8Ftx，x˙=Fk+Fc=∑i=1mk2i-1（A，f）x2i-1+   C（A，f）x˙nA，fsgn（x˙）x=Asinωt=Asin2πftx˙=Aωcos ωt=2πfA2-x2sgnx˙ （1）式中，Fk为非线性弹性恢复力；Fc为迟滞非线性阻尼力；k2i-1A，f为刚度函数；CA，f为阻尼函数；nA，f为阻尼成分函数。当n=0时，仅有干摩擦阻尼；当0n≤1时，既有干摩擦阻尼又有黏性阻尼；当n1时，表示有高阶阻尼[13]7。膜片联轴器摩擦主要为干摩擦阻尼，但联轴器受载后，膜片间相对接触压力及相对滑动量呈现非线性变化。因此，将联轴器干摩擦阻尼力等效为黏性阻尼力，通过求解1个周期内阻尼耗能并以此求解出等效黏滞阻尼系数来表征膜片联轴器的阻尼特性。由图9所示膜片联轴器扭转动态迟滞回线可知，该回线为对称的类椭圆曲线。因此，将阻尼力以标准椭圆函数拟合，即n取1，则阻尼力满足下列公式Fc=c（A，f）A2-x2sgnx˙c（A，f）=2πfC（A，f）=b(A，f)/a （2）式中，c（A，f）为椭圆短长半轴之比。而椭圆长半轴a即为外加激励振幅A。10.16578/j.issn.1004.2539.2024.04.001.F009图9迟滞回线Fig. 9Hysteresis loop迟滞回线是以原点为中心的倾斜椭圆，其面积S为1个周期内消耗的能量。图9中，Fm为最大位移时所受的激励力；A为最大位移；F0为位移为0时所受激励力；a为长半轴的长度；b为短半轴的长度。轴向振动动态迟滞回线越靠近位移零点，摩擦作用越小，摩擦阻尼耗能越少。因此，将阻尼力以标准椭圆函数拟合与实际情况不符。为此，在式（1）迟滞非线性阻尼力项中引入高阶阻尼的影响，则迟滞非线性阻尼力变为Fc=∑i=1mC2i-1A，fx˙2i-1sgnx˙ （3）2.2.2 动刚度与等效阻尼系数在联轴器工程实际中，为方便计算及对比，常将动刚度曲线等效为频率范围内的某个具体值，即等效动刚度。等效动刚度计算式为Ke=2∫0AFkdxA2 （4）每个周期消耗的能量计算公式[12]2996为ΔU=2∫-AAFcdx，x˙≥0 （5）等效阻尼系数为c=ΔU∫0Tx˙2dt （6）若迟滞非线性阻尼力曲线为标准椭圆曲线，则单周期消耗的能量及等效阻尼系数公式可分别简化为ΔU=πabc=baω （7）2.2.3 扭转动刚度与扭转等效阻尼系数图10所示为基于式（1）~式（3）对膜片联轴器扭转和轴向正弦激励仿真迟滞回线进行非线性最小二乘拟合的曲线，拟合优度采取R2系数表征。结果表明，拟合曲线与仿真计算曲线吻合良好且拟合优度均达到0.99以上，说明拟合模型在求解联轴器非线性迟滞恢复力时具有良好的适用性。10.16578/j.issn.1004.2539.2024.04.001.F010图10扭转及轴向迟滞回线拟合结果Fig. 10Fitting results of torsional and axial hysteresis loops基于式（4）~式（7），计算膜片联轴器的扭转等效动刚度Ket和等效阻尼系数ct、轴向等效动刚度Kea和等效阻尼系数ca随激励幅值及频率变化情况，结果分别如表2、表3所示。10.16578/j.issn.1004.2539.2024.04.001.T002表2不同载荷幅值下扭转和轴向等效动刚度、等效阻尼系数Tab. 2Equivalent dynamic stiffness and damping coefficient of torsion and axial directions under different load amplitudes类型幅值A/（N·m）频率f/Hz等效动刚度Ket/（N·m/rad）等效阻尼系数ct/ （N·m·s/rad）扭转200103.55×10556.32400103.79×105104.1500103.76×105318.391 000103.59×105511.71 500103.48×105379.1类型幅值A/mm频率f/Hz等效动刚度Kea/（N/m）等效阻尼系数ca/ （N·s/m）轴向0.555.75×10581.321.056.57×10587.631.558.06×105148.272.051.09×106244.132.551.71×106398.0810.16578/j.issn.1004.2539.2024.04.001.T003表3不同载荷频率下扭转和轴向等效动刚度、等效阻尼系数Tab. 3Equivalent dynamic stiffness and damping coefficient of torsion and axial directions under different load frequencies类型幅值A/（N·m）频率f/Hz等效动刚度Ket/（N·m/rad）等效阻尼系数ct/ （N·m·s/rad）扭转1 000103.59×105511.71 000203.587×105253.91 000303.582×105181.31 000403.577×105120.41 000503.565×105101.8类型幅值A/mm频率f/Hz等效动刚度Kea/（N/m）等效阻尼系数ca/ （N·s/m）轴向1106.57×10543.71206.56×10521.951306.54×10514.121406.52×10510.451506.49×1058.23由表3可知：1）振动频率远低于共振频率时，膜片联轴器扭转等效动刚度与扭转静刚度一致；随着转矩的增大，呈先非线性增长后再减小的趋势。扭转等效阻尼系数随转矩变化情况同等效扭转刚度一致；而轴向等效动刚度也同轴向静刚度一致，呈非线性指数型增长趋势。轴向等效阻尼系数随轴向位移的增大，呈非线性递增趋势。2）振动频率远低于共振频率时，膜片联轴器的扭转和轴向等效动刚度基本不变，而等效阻尼系数随着激励频率的增大而逐渐减小，且下降趋势变缓。该变化结果也符合式（6）、式（7）中等效阻尼系数与激励频率成反比的特性。3 膜片联轴器刚度、阻尼试验分析3.1　试验样件与设备为了验证仿真结果的合理性，对膜片联轴器开展了静态扭转、拉伸以及动态扭转、拉伸试验。试验设备与试验件分别如图11、图12所示。10.16578/j.issn.1004.2539.2024.04.001.F011图11电液伺服扭转疲劳试验机（左）和动静万能试验机（右）Fig. 11Electro-hydraulic servo torsion fatigue testing machine （left） and dynamic and static universal testing machine （right）10.16578/j.issn.1004.2539.2024.04.001.F012图12试验件Fig. 12Test piece试验装置采用PNW-5型电液伺服扭转疲劳试验机及PA-100型电液伺服动静万能试验机。扭转试验台最大转矩为±5 000 N∙m，转矩测量精度为示值的±1%，最大扭转角为±50°（分辨率为±0.02°），频率范围为0.001~100 Hz。轴向拉压试验台试验力为±100 kN，测量误差为示值的1%，频率范围为0.1~100 Hz（分辨率为0.01 Hz）。3.2　静刚度试验3.2.1　扭转静刚度试验基于扭转及拉伸试验台，开展联轴器静态扭转及拉伸试验。由于根据试验台安装夹具设计出的法兰与标准膜片联轴器法兰结构不同，为验证仿真模型及结构的正确性，基于试验件形状建立新的膜片联轴器仿真模型，不改变仿真接触条件，仅改变法兰结构外形，如图13所示。仿真载荷及约束条件为在输入、输出法兰盘端面的面中心添加远程控制点，仅释放其转动自由度。输出法兰与试验台连接的螺栓孔设为固定约束。输入法兰与试验台连接的螺栓孔耦合到面中心远程控制点，在输入法兰面中心远程控制点施加转矩。10.16578/j.issn.1004.2539.2024.04.001.F013图13试验件仿真模型Fig. 13Simulation model of the test piece利用扭转试验台进行扭转静刚度试验验证。在图12所示扭转主动法兰上缓慢施加转矩到2 200 N·m，利用角位移传感器采集加载过程中主动法兰盘的转角。绘制的转角及扭转刚度随转矩变化的曲线如图14所示。10.16578/j.issn.1004.2539.2024.04.001.F014图14静态扭转试验及仿真曲线Fig. 14Static torsion test and simulation curves由于在试验件与试验台安装的夹具之间沿周向存在多处间隙，初始阶段转矩增加，试验转角曲线φ试验转角增量大于仿真转角曲线φ仿真，反映到刚度值上，就是试验扭转刚度出现较大幅度下降的情况。随着转矩进一步增大，联轴器各结构件接触变形产生相应转矩，抵消输入转矩，试验与仿真转角增长幅度基本相同，转角曲线基本平行。扭转刚度曲线随转矩增大呈基本靠拢趋势。在不同转矩下，忽略试验初始阶段扭转刚度，计算对比两者等转矩间隔内的平均扭转刚度，得到仿真与试验扭转静刚度相对误差，如表4所示。结果表明，膜片联轴器仿真模型及参数设置较为合理。10.16578/j.issn.1004.2539.2024.04.001.T004表4扭转静刚度仿真与试验结果Tab. 4Simulation and test results of the torsional static stiffness转矩/（N·m）扭转静刚度/（105 N·m/rad）误差/%仿真试验400~6001.541.603.75600~1 0001.491.462.051 000~1 4001.471.432.801 400~1 8001.461.413.553.2.2　轴向静刚度试验利用轴向位移试验台进行轴向静刚度试验验证。从上端连接盘缓慢施加轴向位移载荷到2 mm，利用传感器采集加载过程中连接盘受到的支反力，绘制支反力及轴向刚度随转矩变化的曲线，如图15所示。膜片静态轴向拉压试验与仿真轴向支反力-轴向位移曲线基本重合；拉压轴向刚度曲线呈指数型增长趋势，且受压刚度大于拉伸刚度，与仿真结果一致。10.16578/j.issn.1004.2539.2024.04.001.F015图15静态轴向拉压试验及仿真曲线Fig. 15Static axial tension and compression test and simulation curves计算对比两者等轴向位移间隔内的平均刚度，得到仿真与试验轴向静刚度的相对误差，如表5所示。结果表明，膜片联轴器仿真模型及参数设置准确性较合理。10.16578/j.issn.1004.2539.2024.04.001.T005表5轴向静刚度仿真与试验结果Tab. 5Simulation and test results of the axial static stiffness位移/mm轴向静刚度/（105 N/m）误差/%仿真试验-2~-1.219.8020.352.70-1.2~-0.49.239.735.14-0.4~04.74.82.080~0.45.335.635.330.4~1.28.698.983.231.2~217.7618.463.793.3　动刚度、阻尼试验3.3.1　扭转动刚度、阻尼试验转矩以正弦波的形式加载，得到不同激励幅值及频率下膜片联轴器的转矩及转角值，经低通滤波降噪处理，绘制扭转振动迟滞回线，如图16所示。10.16578/j.issn.1004.2539.2024.04.001.F016图16扭转振动迟滞回线随幅值与频率变化的趋势Fig. 16Trend chart of torsional vibration hysteresis loops changing with amplitude and frequency联轴器扭转简谐振动迟滞回线同仿真结果一致，扭转激励幅值越大，迟滞回线包围面积越大，摩擦阻尼耗能随振动幅值的增大而增大。曲线位移幅值端点连线斜率降低，表明联轴器扭转动刚度随扭振幅值增大而减小；整体上迟滞回线随扭振激励频率的变化较小，表明低频率下阻尼耗能及刚度基本上与激励频率无关。试验获得的膜片联轴器扭转等效动刚度及阻尼系数随转矩幅值及频率变化的趋势如图17所示。10.16578/j.issn.1004.2539.2024.04.001.F017图17扭振等效刚度及阻尼系数随幅值及频率变化的趋势Fig. 17Trend chart of the torsional equivalent stiffness and the damping coefficient changing with amplitude and frequency扭转等效动刚度随激励幅值增大呈先增后减的非线性趋势。小激励频率范围内，扭转等效动刚度随频率增大而略微下降；扭转等效阻尼系数随激励幅值增大呈先增后减的非线性趋势，随频率增大呈非线性下降趋势。3.3.2　轴向动刚度、阻尼试验轴向以正弦波的形式加载，得到不同激励幅值及频率下膜片联轴器的支反力及位移值，经低通滤波降噪处理，绘制轴向振动迟滞回线，如图18所示。10.16578/j.issn.1004.2539.2024.04.001.F018图18轴向振动迟滞回线随幅值与频率变化的趋势Fig. 18Trend chart of the hysteresis loop of axial vibration changing with amplitude and frequency联轴器轴向简谐振动迟滞回线呈S形，随振幅增大，迟滞回线包围面积增大，且靠近位移零点，曲线紧缩越明显。曲线位移幅值端点连线斜率增大，表明联轴器平均轴向动刚度随轴向幅值增大而非线性增大；整体上轴向振动迟滞回线包含面积随激励频率的变化较小。由于联轴器轴向刚度较小，轴向振动固有频率较低，迟滞回线支反力极值点下降明显，表明轴向等效动刚度随频率增大而减小。试验获得的膜片联轴器轴向等效动刚度及阻尼系数随轴向位移幅值及频率变化趋势如图19所示。10.16578/j.issn.1004.2539.2024.04.001.F019图19轴向振动等效刚度及阻尼系数随幅值及频率变化的趋势Fig. 19Trend chart of the equivalent stiffness and the damping coefficient of axial vibration changing with amplitude and frequency轴向等效动刚度随激励幅值增大呈非线性增长趋势，而随频率增大呈非线性下降趋势；轴向等效阻尼系数随激励幅值增大呈非线性增长，增长趋势逐渐变缓，随激励频率增大呈非线性下降趋势。对比分析膜片联轴器动态扭转及轴向等效动刚度与阻尼系数的仿真与试验结果，两者变化趋势一致。但由于联轴器法兰件的仿真模型与试验件法兰存在尺寸差异，且由于试验台采用液压驱动方式，试验阻尼耗能引入了试验台黏性阻尼耗能，以及仿真各零件之间接触的摩擦因数为查表所得，所以，仿真与试验的动刚度和阻尼在数值上存在误差。4 结论运用有限元仿真法，研究了膜片联轴器静刚度特性以及结构参数对其影响的规律，研究了激励幅值、激励频率对等效动刚度及阻尼系数的影响，并开展扭转与轴向拉伸试验，验证了仿真结果。主要结论与展望如下：1）试验与仿真结果的对比表明，仿真模型及相关接触、边界条件设置基本合理。2）考虑非线性接触以及大变形效应，膜片联轴器4类静刚度中角向刚度和轴向刚度较小，膜片联轴器具有良好的角向和轴向误差补偿作用。由于膜片的薄片特征，承受转矩以及径向载荷时，受压段膜片易发生屈曲失稳现象，导致扭转刚度及径向刚度随载荷增大而减小，且变化幅度较大。3）膜片螺栓孔分布圆直径及膜片厚度对膜片刚度影响较大。可通过增大膜片螺栓孔分布圆直径来提高联轴器承扭及补偿能力；可适当减小膜片厚度来提高其补偿能力。4）类椭圆函数可以较好地构建膜片联轴器动力学模型，解算非线性刚度和阻尼，更好地描述膜片联轴器的迟滞阻尼特性。5）在简谐转矩以及轴向位移载荷激励下，激励幅值的变化对联轴器扭转及轴向等效动刚度和阻尼系数的影响较大；而激励频率的变化对联轴器扭转与轴向振动1周期内的阻尼耗能影响较小，等效阻尼系数与频率成反比关系。