基于传递函数的多缸内燃机轴系扭振数值计算方法：理论、实践与优化

基于传递函数的多缸内燃机轴系扭振数值计算方法：理论、实践与优化一、引言1.1研究背景与意义多缸内燃机作为现代动力系统的核心部件，在汽车、船舶、工程机械、发电设备等众多领域中都发挥着不可替代的关键作用。在汽车领域，多缸内燃机为车辆提供强劲动力，确保其能在不同路况下稳定行驶，满足人们日常出行和货物运输需求；船舶领域，它是船舶航行的动力源泉，支撑着远洋运输和内河航运的高效运作；工程机械依靠多缸内燃机的动力驱动，实现各类复杂工程作业；发电设备中，多缸内燃机则是备用电源和分布式发电的重要选择，保障电力供应的稳定性和可靠性。多缸内燃机工作时，轴系会受到周期性变化的气体压力、往复惯性力和旋转惯性力等的作用，这些力会产生周期性变化的激振力矩，从而使轴系产生扭转振动。轴系扭振会导致一系列严重问题，如轴系零件磨损加剧，降低零件使用寿命，增加维修成本；噪声剧增，影响工作环境和操作人员身心健康；功率下降，降低设备工作效率；甚至可能导致断轴等严重事故，危及设备安全和人员生命。在轴系扭振数值计算中，传递函数是一种极为重要的工具。传递函数能够清晰地描述系统输入与输出之间的关系，将系统的动态特性以数学形式展现出来。通过传递函数，可深入分析系统的固有频率、振型等重要参数，准确预测轴系在不同工况下的扭振响应。在实际应用中，传递函数能为轴系的优化设计提供科学依据，指导工程师合理调整轴系结构和参数，有效降低扭振幅值，提高轴系的稳定性和可靠性。研究基于传递函数的多缸内燃机轴系扭振数值计算方法具有重大的理论和实际意义。理论层面，能进一步完善多缸内燃机动力学理论体系，加深对轴系扭振复杂现象的认识和理解，为后续研究提供坚实的理论基础。实际应用中，该研究成果可直接应用于多缸内燃机的设计研发过程，帮助工程师提前发现轴系扭振问题并采取有效措施加以解决，从而缩短研发周期，降低研发成本，提高产品质量和性能。对于运行中的多缸内燃机，基于传递函数的轴系扭振数值计算方法可用于状态监测和故障诊断，及时发现轴系扭振异常，预测故障发生，为设备维护和检修提供准确依据，保障设备安全稳定运行，减少因故障导致的停机损失和生产事故。1.2国内外研究现状在多缸内燃机轴系扭振数值计算方法领域，国内外学者已开展了大量研究，取得了一系列丰硕成果。国外方面，早期就有学者关注到内燃机轴系扭振问题，并开展了相关研究。随着计算机技术的飞速发展，数值计算方法逐渐成为研究轴系扭振的重要手段。在基于传递函数的轴系扭振数值计算方法研究中，国外学者提出了多种理论和算法。文献[具体文献1]通过建立详细的轴系动力学模型，深入探讨了传递函数在轴系扭振分析中的应用，推导出了计算轴系固有频率和振型的传递函数表达式，为后续研究奠定了坚实的理论基础。该研究利用传递函数对轴系进行模态分析，清晰地揭示了轴系各阶固有频率和对应的振型特征，使研究者能够直观地了解轴系的振动特性。在实际应用方面，文献[具体文献2]将基于传递函数的计算方法应用于某大型船舶多缸内燃机轴系扭振分析中，通过精确计算轴系在不同工况下的扭振响应，成功预测出轴系可能出现共振的转速范围，为船舶动力系统的安全运行提供了有力保障。研究过程中，通过对船舶实际运行数据的采集和分析，验证了基于传递函数的计算方法的准确性和可靠性。国内在该领域的研究起步相对较晚，但发展迅速。众多高校和科研机构投入大量资源开展相关研究，取得了显著进展。在理论研究方面，文献[具体文献3]针对国内常见的多缸内燃机结构特点，对传递函数的计算方法进行了优化和改进，提高了计算精度和效率。该研究通过引入新的数学模型和算法，有效解决了传统传递函数计算方法中存在的精度不足和计算效率低下的问题，使得基于传递函数的轴系扭振数值计算方法更加适用于国内多缸内燃机的实际情况。在实验研究方面，文献[具体文献4]搭建了多缸内燃机轴系扭振实验台，通过实验测量轴系的扭振响应，并与基于传递函数的数值计算结果进行对比分析，进一步验证了计算方法的正确性和可靠性。实验过程中，对实验数据进行了详细的处理和分析，通过对比实验结果和数值计算结果，深入探讨了实验误差的来源和影响因素，为进一步改进实验方法和提高计算精度提供了重要依据。尽管国内外在基于传递函数的多缸内燃机轴系扭振数值计算方法研究方面已取得诸多成果，但仍存在一些不足之处。部分研究在建立轴系模型时，对一些复杂因素的考虑不够全面，如轴系的非线性特性、部件之间的接触刚度变化等，这些因素可能会对轴系扭振响应产生较大影响，导致计算结果与实际情况存在一定偏差。在计算效率方面，对于一些大型复杂的多缸内燃机轴系，现有的计算方法可能需要耗费大量的计算时间和计算资源，难以满足实际工程快速分析的需求。在实验研究方面，实验设备和实验方法仍有待进一步完善，以获取更准确、更全面的实验数据，为理论研究和数值计算提供更有力的支持。目前的研究在不同工况下轴系扭振的瞬态响应研究还不够深入，对于轴系在启动、加速、减速等动态过程中的扭振特性了解有限，这对于全面评估轴系的可靠性和安全性具有一定的局限性。1.3研究目标与内容本研究旨在深入探究基于传递函数的多缸内燃机轴系扭振数值计算方法，通过理论分析、实例计算和算法优化等多方面研究，完善该计算方法，提高其计算精度和效率，为多缸内燃机轴系的设计、优化和故障诊断提供更为准确可靠的理论依据和技术支持。在理论分析方面，将全面深入研究传递函数在多缸内燃机轴系扭振分析中的理论基础。详细推导传递函数的表达式，深入剖析其物理意义，清晰阐述传递函数与轴系扭振特性之间的内在联系。全面考虑轴系的各种复杂因素，如轴系的非线性特性、部件之间的接触刚度变化、阻尼特性等，建立更为精确的轴系动力学模型。在该模型中，充分考虑轴系各部件的实际结构和工作状态，对轴系的弹性、惯性、阻尼等参数进行准确描述，确保模型能够真实反映轴系的扭振特性。通过理论分析，深入研究这些复杂因素对传递函数和轴系扭振响应的影响规律，为后续的数值计算提供坚实的理论支撑。基于建立的轴系动力学模型，运用数值计算方法求解轴系的扭振响应。选择合适的数值计算方法，如有限元法、模态叠加法等，并对这些方法进行详细的比较和分析，确定最适合本研究的数值计算方法。利用选定的数值计算方法，对不同类型和工况下的多缸内燃机轴系进行扭振响应计算。在计算过程中，严格控制计算参数，确保计算结果的准确性和可靠性。将计算结果与实验数据或实际运行数据进行对比验证，详细分析计算结果与实际情况之间的差异，找出产生差异的原因，提出相应的改进措施，以提高计算方法的准确性。为了进一步提高基于传递函数的多缸内燃机轴系扭振数值计算方法的计算效率，对现有的计算算法进行优化。针对大型复杂的多缸内燃机轴系，研究采用并行计算技术、降阶模型等方法，有效减少计算时间和计算资源的消耗。在并行计算技术方面，充分利用多核处理器的优势，将计算任务合理分配到各个处理器核心上，实现并行计算，提高计算效率。在降阶模型方面，通过合理简化轴系模型，去除对扭振响应影响较小的部件和因素，建立降阶模型，在保证计算精度的前提下，大幅减少计算量。通过优化算法，使计算方法能够更好地满足实际工程快速分析的需求，为多缸内燃机的设计和运行提供更高效的技术支持。通过对多缸内燃机轴系扭振的研究，为轴系的优化设计提供科学依据。根据扭振计算结果，深入分析轴系的薄弱环节和潜在问题，提出针对性的优化措施。例如，通过调整轴系的结构参数，如轴径、轴长、惯量分布等，改变轴系的固有频率，避免共振的发生；优化轴系的材料选择和制造工艺，提高轴系的强度和刚度，降低扭振应力；合理设计轴系的减振装置，如采用阻尼减振器、硅油减振器等，有效抑制扭振振幅，提高轴系的稳定性和可靠性。通过优化设计，降低轴系扭振对内燃机性能和可靠性的影响，提高多缸内燃机的整体性能和使用寿命。1.4研究方法与技术路线本研究采用理论分析、数值模拟和实验验证相结合的研究方法，从多个角度深入探究基于传递函数的多缸内燃机轴系扭振数值计算方法。在理论分析方面，深入剖析传递函数在多缸内燃机轴系扭振分析中的理论根基。全面考虑轴系的各种复杂因素，如轴系的非线性特性、部件之间的接触刚度变化、阻尼特性等，推导传递函数的表达式，明确其物理意义，阐释传递函数与轴系扭振特性的内在联系，为后续研究筑牢理论基础。通过理论分析，揭示轴系扭振的基本规律，为数值计算和实验研究提供理论指导。运用数值模拟方法，基于建立的轴系动力学模型，求解轴系的扭振响应。选择有限元法、模态叠加法等合适的数值计算方法，并对其进行细致的比较和分析，确定最契合本研究的方法。利用选定方法，对不同类型和工况下的多缸内燃机轴系进行扭振响应计算，将计算结果与实验数据或实际运行数据对比验证，分析差异原因，提出改进措施，提高计算方法的准确性。数值模拟能够快速、高效地获取轴系扭振的相关数据，为理论分析和实验研究提供有力支持。为了验证理论分析和数值模拟的结果，开展实验研究。搭建多缸内燃机轴系扭振实验台，模拟实际工况，测量轴系的扭振响应。将实验结果与理论计算和数值模拟结果进行对比分析，进一步验证基于传递函数的多缸内燃机轴系扭振数值计算方法的正确性和可靠性。实验研究能够获取真实的轴系扭振数据，为理论和数值研究提供验证依据，同时也能发现理论和数值研究中未考虑到的问题，推动研究的深入发展。本研究的技术路线如图1所示：确定研究目标与内容：明确基于传递函数的多缸内燃机轴系扭振数值计算方法的研究目标，确定理论分析、数值计算和实验验证等具体研究内容。理论分析：研究传递函数理论，考虑轴系复杂因素，推导传递函数表达式，分析其与扭振特性的联系。数值计算方法选择与模型建立：对比有限元法、模态叠加法等数值计算方法，确定合适方法，建立轴系动力学模型。数值计算与结果分析：运用选定方法计算不同工况下轴系扭振响应，与实验或实际数据对比，分析差异并改进。实验验证：搭建实验台，测量轴系扭振响应，对比实验与理论、数值计算结果，验证方法的正确性和可靠性。结果总结与应用：总结研究成果，提出轴系优化设计建议，将研究成果应用于实际工程。[此处插入技术路线图][此处插入技术路线图]通过上述研究方法和技术路线，本研究将全面深入地探究基于传递函数的多缸内燃机轴系扭振数值计算方法，为多缸内燃机轴系的设计、优化和故障诊断提供坚实的理论依据和技术支持。二、多缸内燃机轴系扭振及传递函数理论基础2.1多缸内燃机轴系扭振概述2.1.1扭振的产生原因多缸内燃机工作时，轴系扭振的产生是多种复杂因素共同作用的结果，其中气体压力、往复惯性力和旋转惯性力起着关键作用。在每个工作循环中，气缸内的气体压力会呈现出周期性的剧烈变化。当活塞处于压缩冲程末期时，气缸内气体被压缩，压力急剧升高；在做功冲程，气体迅速膨胀，推动活塞下行，此时气体压力达到峰值，对活塞产生强大的推力，并通过连杆传递给曲轴，形成周期性变化的激振力矩。对于四冲程内燃机，每完成一个工作循环，曲轴旋转两圈，每个气缸的气体压力变化会产生特定频率的激振力，多个气缸的激振力相互叠加，使得轴系受到复杂的周期性激励，从而引发扭振。内燃机的活塞、连杆等部件在工作过程中做往复直线运动，由于其速度和加速度的不断变化，会产生往复惯性力。当活塞从气缸一端运动到另一端时，速度先从零逐渐增大，然后又逐渐减小到零，加速度方向也相应改变，这就导致往复惯性力的大小和方向呈周期性变化。这种周期性变化的往复惯性力通过连杆作用于曲轴，成为轴系扭振的又一重要激振源。旋转部件，如曲轴、飞轮等，由于质量分布不均匀或制造误差，在旋转过程中会产生旋转惯性力。即使部件质量分布均匀，在高速旋转时，由于离心力的作用，也会对轴系产生附加的动态载荷。这些旋转惯性力同样会引起轴系的扭振。除了上述主要因素外，轴系自身的结构和特性也对扭振的产生有重要影响。轴系的弹性使其在受到外力作用时能够发生弹性变形，为扭振的产生提供了条件。当轴系受到周期性的激振力矩时，轴段会发生扭转弹性变形，储存弹性势能；当激振力消失或方向改变时，轴段又会依靠弹性恢复力回到原来的位置，在这个过程中，弹性势能和动能相互转化，形成扭振。轴系的阻尼特性则对扭振起着抑制作用，它能够消耗扭振过程中的能量，使扭振振幅逐渐减小。实际轴系的阻尼来源较为复杂，包括轴系各部件之间的摩擦阻尼、润滑油的粘性阻尼以及材料内部的阻尼等。不同类型和工况下的内燃机，其轴系的阻尼特性会有所差异，从而对扭振的影响也不尽相同。2.1.2扭振的危害轴系扭振会对多缸内燃机产生诸多不良影响，严重威胁其零部件的疲劳寿命、传动系统的稳定性以及整机的性能和可靠性。轴系扭振会导致轴系零部件承受交变应力，加速其疲劳磨损。当轴系发生扭振时，轴段会产生周期性的扭转变形，使得轴表面的材料受到拉伸和压缩交替作用，从而产生交变应力。这种交变应力会在轴的表面和内部形成微观裂纹，随着裂纹的逐渐扩展，最终导致零部件疲劳断裂。对于曲轴、凸轮轴等关键部件，疲劳断裂可能会引发严重的故障，导致内燃机停机甚至损坏。扭振还会加剧轴与轴承之间的摩擦和磨损，降低轴承的使用寿命，增加维修成本。在一些高速运转的内燃机中，轴系扭振引起的轴承磨损问题尤为突出，需要频繁更换轴承，影响设备的正常运行。轴系扭振会对内燃机的传动系统稳定性造成严重影响。当扭振发生时，轴系的转速会出现波动，导致传动系统中的齿轮、链条等部件受到冲击载荷。这会使齿轮之间的啮合不稳定，产生噪声和振动，加速齿轮的磨损，甚至导致齿轮断齿。扭振还会影响到内燃机与其他设备之间的动力传递，降低传动效率，影响整个系统的工作性能。在汽车发动机中，轴系扭振可能会导致变速器换挡困难、离合器片磨损加剧等问题，影响车辆的驾驶舒适性和操控稳定性。轴系扭振会降低内燃机的整机性能和可靠性。扭振会消耗一部分能量，导致内燃机输出功率下降，燃油经济性变差。扭振还会引起内燃机的振动和噪声增大，影响工作环境和操作人员的身心健康。在一些对噪声要求较高的场合，如船舶机舱、发电机组等，轴系扭振产生的噪声问题需要采取额外的降噪措施来解决。严重的扭振还可能导致内燃机的工作循环不稳定，影响燃烧过程，增加污染物排放，降低内燃机的环保性能。如果扭振问题得不到及时解决，长期积累可能会引发更严重的故障，如内燃机失控、爆炸等，危及设备安全和人员生命。2.1.3常见的扭振分析方法在多缸内燃机轴系扭振研究中，Holzer表法、系统矩阵法、传递矩阵法等是常用的分析方法，每种方法都有其独特的优缺点。Holzer表法是一种经典的轴系扭振分析方法，具有概念清晰、计算过程直观的优点。它基于轴系的动力学基本方程，通过列表计算的方式，逐步求解轴系各节点的扭振参数。在计算过程中，需要根据轴系的结构和参数，列出各节点的运动方程和力平衡方程，然后通过迭代计算，得到轴系的固有频率和振型。这种方法适用于简单轴系的扭振分析，对于初学者来说，易于理解和掌握。当轴系结构较为复杂，节点数量较多时，Holzer表法的计算过程会变得繁琐，计算工作量大幅增加，而且容易出现计算错误。系统矩阵法将轴系视为一个整体系统，通过建立系统的动力学方程，利用矩阵运算求解轴系的扭振特性。该方法能够全面考虑轴系各部件之间的相互作用，适用于复杂轴系的分析。在建立系统矩阵时，需要考虑轴系的刚度矩阵、质量矩阵和阻尼矩阵等，通过求解系统的特征方程，可以得到轴系的固有频率和振型。系统矩阵法的计算精度较高，但对计算资源的要求也较高，需要使用专业的数值计算软件进行求解。由于系统矩阵的建立较为复杂，需要对轴系的结构和参数有深入的了解，因此对于一些复杂的轴系，应用系统矩阵法进行分析具有一定的难度。传递矩阵法是一种基于能量守恒原理的轴系扭振分析方法，它将轴系划分为若干个单元，通过建立单元之间的传递关系，求解轴系的扭振响应。传递矩阵法的计算效率较高，能够快速得到轴系的固有频率和振型。在计算过程中，只需要知道轴系各单元的参数，就可以通过传递矩阵的运算，得到整个轴系的扭振特性。这种方法适用于各种类型的轴系，尤其是对于大型复杂轴系的分析，具有明显的优势。传递矩阵法在处理非线性问题时存在一定的局限性，对于一些含有非线性元件的轴系，需要进行特殊处理才能应用传递矩阵法进行分析。这些常见的扭振分析方法在多缸内燃机轴系扭振研究中都发挥着重要作用。在实际应用中，需要根据轴系的具体结构和分析要求，选择合适的分析方法，以准确、高效地求解轴系的扭振特性。2.2传递函数基本原理2.2.1传递函数的定义与数学表达传递函数是描述线性定常系统动态特性的重要数学工具，在零初始条件下，其定义为系统输出量的拉普拉斯变换与引起该输出的输入量的拉普拉斯变换之比。对于一个线性定常系统，若输入量为x(t)，输出量为y(t)，则传递函数G(s)可表示为：G(s)=\frac{Y(s)}{X(s)}其中，Y(s)=L[y(t)]是输出量y(t)的拉普拉斯变换，X(s)=L[x(t)]是输入量x(t)的拉普拉斯变换，s为复变量。以一个简单的二阶线性定常系统为例，其动力学方程为：m\ddot{y}(t)+c\dot{y}(t)+ky(t)=f(t)式中，m为质量，c为阻尼系数，k为刚度系数，f(t)为输入力，y(t)为输出位移。在零初始条件下，对上述方程两边进行拉普拉斯变换，根据拉普拉斯变换的微分性质：L[\dot{y}(t)]=sY(s)，L[\ddot{y}(t)]=s^{2}Y(s)，可得：ms^{2}Y(s)+csY(s)+kY(s)=F(s)整理后得到该系统的传递函数为：G(s)=\frac{Y(s)}{F(s)}=\frac{1}{ms^{2}+cs+k}从这个例子可以看出，传递函数将系统的输入输出关系以复频域的形式简洁地表达出来，分母中的各项系数m、c、k取决于系统本身的固有特性，而分子通常为1（在这个简单的力学系统中），它反映了系统与外界激励的关系。通过传递函数，我们可以方便地分析系统对不同输入信号的响应特性，为系统的设计、分析和优化提供有力的工具。2.2.2传递函数在振动系统中的应用在振动系统中，传递函数具有重要的应用价值，它能够深入描述系统的动态特性，为系统的分析和设计提供关键依据。传递函数可以用于求解系统的响应。当已知振动系统的传递函数G(s)和输入激励的拉普拉斯变换X(s)时，通过公式Y(s)=G(s)X(s)，即可求得系统输出响应的拉普拉斯变换Y(s)，再对Y(s)进行拉普拉斯反变换，就能得到系统在时域中的响应y(t)。对于一个受到简谐激励力F(t)=F_0\sin(\omegat)作用的单自由度振动系统，其传递函数为G(s)=\frac{1}{ms^{2}+cs+k}，先将激励力F(t)进行拉普拉斯变换得到F(s)，然后通过上述公式计算Y(s)，再反变换得到系统的位移响应y(t)，从而清晰地了解系统在该激励下的振动情况。传递函数在分析系统的频率响应方面发挥着关键作用。频率响应是指系统对不同频率输入信号的稳态响应特性。将传递函数G(s)中的复变量s用j\omega（j为虚数单位，\omega为角频率）代替，得到频率响应函数G(j\omega)，其模\vertG(j\omega)\vert表示系统在频率\omega处的幅值放大倍数，相位\angleG(j\omega)表示输出与输入之间的相位差。通过绘制频率响应曲线，即\vertG(j\omega)\vert-\omega和\angleG(j\omega)-\omega曲线，可以直观地了解系统对不同频率激励的响应特性，确定系统的共振频率、带宽等重要参数。当\vertG(j\omega)\vert达到最大值时，对应的频率即为共振频率，此时系统的振动最为剧烈。传递函数还广泛应用于振动系统的模态分析。模态分析是研究系统固有频率和振型的重要方法。通过求解系统的特征方程\vertsI-A\vert=0（对于状态空间模型，A为系统矩阵，I为单位矩阵），得到系统的固有频率\omega_n，而振型则可以通过求解相应的特征向量得到。传递函数与系统的模态密切相关，其极点（即传递函数分母为零的根）对应着系统的固有频率。通过分析传递函数的极点和零点分布，可以深入了解系统的模态特性，判断系统的稳定性。若传递函数的极点具有负实部，则系统是稳定的；若存在正实部的极点，则系统不稳定。2.2.3传递函数与多缸内燃机轴系扭振的联系将传递函数应用于多缸内燃机轴系扭振分析，能够建立起轴系扭振系统的传递函数模型，从而深入研究轴系的扭振特性。在多缸内燃机轴系中，把轴系视为一个线性系统，轴系所受的周期性激振力矩作为输入，轴系各节点的扭角或扭矩作为输出。通过对轴系的动力学分析，可建立其传递函数模型。考虑一个简单的多缸内燃机轴系集中质量模型，将轴系划分为若干个集中质量和弹性轴段，每个集中质量具有转动惯量J_i，相邻集中质量之间的弹性轴段具有扭转刚度k_i和阻尼系数c_i。根据牛顿第二定律和轴系的力学关系，可列出轴系的运动微分方程。对于第i个集中质量，其运动方程为：J_i\ddot{\theta}_i(t)+c_i(\dot{\theta}_i(t)-\dot{\theta}_{i-1}(t))+k_i(\theta_i(t)-\theta_{i-1}(t))-c_{i+1}(\dot{\theta}_{i+1}(t)-\dot{\theta}_i(t))-k_{i+1}(\theta_{i+1}(t)-\theta_i(t))=T_i(t)其中，\theta_i(t)为第i个集中质量的扭角，T_i(t)为作用在第i个集中质量上的激振力矩。在零初始条件下，对上述运动微分方程进行拉普拉斯变换，经过一系列的数学推导和整理，可以得到轴系扭振系统的传递函数矩阵[G(s)]。传递函数矩阵中的元素G_{ij}(s)表示第j个输入（激振力矩）对第i个输出（扭角或扭矩）的传递函数关系。通过这个传递函数模型，能够计算轴系在不同工况下的扭振响应，分析轴系的固有频率和振型。当给定轴系的激振力矩频谱[T(s)]时，通过矩阵运算[\theta(s)]=[G(s)][T(s)]，可得到轴系各节点扭角的拉普拉斯变换[\theta(s)]，再进行拉普拉斯反变换，就能得到轴系各节点在时域中的扭角响应\theta(t)。通过分析传递函数的极点，可确定轴系的固有频率，进而判断轴系在运行过程中是否会发生共振现象，为轴系的设计和优化提供重要依据。三、基于传递函数的多缸内燃机轴系扭振数值计算模型构建3.1轴系当量系统的建立3.1.1当量系统的概念与转化原则在多缸内燃机轴系扭振分析中，由于实际轴系结构复杂，包含众多形状不规则、材料特性各异的部件，直接对其进行精确的动力学分析难度极大。为简化分析过程，提高计算效率，引入轴系当量系统的概念。轴系当量系统是将实际复杂的轴系简化为一个由若干集中质量、弹性轴段和阻尼元件组成的等效系统，该系统在动力学特性上与实际轴系基本相同，能够反映实际轴系的扭振特性。将实际轴系转化为当量系统时，需遵循严格的等效原则，以确保当量系统能准确替代实际轴系进行扭振分析。转动惯量等效是重要原则之一，要求当量系统中各集中质量的转动惯量之和等于实际轴系中相应部件的转动惯量之和。对于内燃机的曲轴，其结构复杂，包含多个曲拐和轴颈，在转化为当量系统时，需将每个曲拐和轴颈的转动惯量合理分配到对应的集中质量上，使集中质量的转动惯量总和与实际曲轴的转动惯量相等。通过精确计算和合理分配转动惯量，保证当量系统在惯性特性上与实际轴系一致，从而准确反映轴系在扭振过程中的惯性作用。刚度等效原则确保当量系统中弹性轴段的扭转刚度与实际轴系中相应轴段的扭转刚度相等。轴系的扭转刚度决定了轴在受到扭矩作用时的扭转变形程度，对于多缸内燃机轴系，不同轴段的材料、尺寸和结构差异会导致其扭转刚度各不相同。在建立当量系统时，需根据实际轴段的材料特性（如弹性模量、剪切模量）、几何尺寸（如轴径、轴长），运用材料力学和弹性力学的相关理论，准确计算各轴段的扭转刚度，并将其等效为当量系统中弹性轴段的刚度。对于连接曲轴和飞轮的轴段，通过计算其实际的扭转刚度，在当量系统中设置具有相同刚度的弹性轴段，以保证当量系统能正确模拟实际轴系的弹性变形特性。阻尼等效也是必不可少的原则。实际轴系中的阻尼来源广泛，包括部件之间的摩擦阻尼、润滑油的粘性阻尼以及材料内部的阻尼等。在当量系统中，需综合考虑这些阻尼因素，通过合理的方法确定阻尼元件的阻尼系数，使当量系统的阻尼特性与实际轴系一致。可通过实验测量实际轴系的阻尼特性，或参考相关的经验公式和数据，确定阻尼元件的阻尼系数，确保当量系统在能量耗散方面与实际轴系相符，从而准确模拟轴系扭振过程中的能量衰减现象。通过遵循转动惯量、刚度和阻尼等效原则，建立的轴系当量系统能够在动力学特性上与实际轴系高度相似，为基于传递函数的多缸内燃机轴系扭振数值计算提供准确可靠的模型基础，使后续的计算结果更接近实际情况，为轴系的设计、优化和故障诊断提供有力支持。3.1.2轴系各部件的当量参数计算轴系各部件的当量参数计算是建立轴系当量系统的关键环节，下面分别对曲轴、飞轮、联轴器、齿轮等主要部件的当量转动惯量、当量刚度和当量阻尼等参数进行详细计算。曲轴作为多缸内燃机轴系的核心部件，其当量参数的准确计算至关重要。曲轴的结构复杂，由多个曲拐、主轴颈和连杆轴颈等组成。计算曲轴的当量转动惯量时，可将曲轴简化为由若干个集中质量组成的系统，每个集中质量的转动惯量根据其对应的曲拐和轴颈的实际结构和质量分布进行计算。对于每个曲拐，可将其质量等效集中到曲拐的中心，根据转动惯量的计算公式J=\sum_{i=1}^{n}m_ir_i^2（其中m_i为第i个质量微元的质量，r_i为该质量微元到转动轴的距离），计算出每个曲拐的转动惯量，然后将所有曲拐和轴颈的转动惯量相加，得到曲轴的当量转动惯量。曲轴的当量刚度主要取决于轴段的材料特性和几何尺寸。根据材料力学中的扭转刚度公式k=\frac{GJ_p}{l}（其中G为材料的剪切弹性模量，J_p为轴的极惯性矩，l为轴段长度），对于不同的轴段，分别计算其极惯性矩J_p=\frac{\pi}{32}d^4（d为轴径），再结合材料的剪切弹性模量和轴段长度，可得到各轴段的扭转刚度。将相邻集中质量之间的轴段刚度作为当量系统中相应弹性轴段的刚度，从而确定曲轴的当量刚度。曲轴的当量阻尼计算较为复杂，涉及多种阻尼因素。可通过实验测量或参考相关经验公式来确定。一种常用的方法是将曲轴的阻尼等效为粘性阻尼，根据实验数据或经验公式确定阻尼系数c，使当量系统中的阻尼力能够近似反映实际曲轴在扭振过程中的能量耗散情况。飞轮通常被视为一个质量集中的刚体，其当量转动惯量可根据飞轮的实际质量m和回转半径r，利用公式J=mr^2进行计算。对于常见的圆盘状飞轮，若已知其外径D和内径d，以及材料密度\rho，则其质量m=\frac{\pi}{4}\rho(D^2-d^2)h（h为飞轮厚度），回转半径r=\sqrt{\frac{D^2+d^2}{8}}，进而计算出飞轮的当量转动惯量。飞轮与曲轴之间的连接轴段的刚度决定了飞轮与曲轴之间的相对扭转刚度，可按照与曲轴轴段刚度计算相同的方法，根据连接轴段的材料特性、几何尺寸计算其当量刚度。飞轮的阻尼相对较小，在一些简化计算中，可忽略不计；若需考虑，可根据实际情况，参考类似结构的阻尼数据，确定一个合理的当量阻尼系数。联轴器在轴系中起到连接和传递扭矩的作用，其当量参数的计算与联轴器的类型密切相关。对于弹性联轴器，通常将其主、从动部分分别视为两个集中质量，其转动惯量可根据联轴器的结构尺寸和材料密度进行计算。弹性联轴器的弹性元件提供了连接部分的扭转刚度，其当量刚度可通过实验测量或根据弹性元件的材料特性和结构参数，利用相关的力学模型进行计算。例如，对于橡胶弹性联轴器，可根据橡胶的弹性模量、橡胶元件的形状和尺寸，计算其扭转刚度。弹性联轴器的阻尼主要来源于弹性元件的内摩擦，可通过实验或经验公式确定其当量阻尼系数。齿轮在轴系中参与动力传递，其当量转动惯量计算时，需考虑齿轮的质量、齿数、模数等因素。对于单个齿轮，可将其质量等效集中到齿圈中心，根据转动惯量计算公式计算其转动惯量。在计算齿轮系统的当量转动惯量时，还需考虑齿轮之间的传动比，通过传动比将从动齿轮的转动惯量折算到主动轴上。齿轮的啮合刚度很大，在一般的轴系扭振计算中，可近似取系统中最大刚度的1000倍作为齿轮啮合的当量刚度。齿轮在啮合过程中会产生一定的阻尼，主要来源于齿面之间的摩擦和润滑油的粘性阻尼。可通过实验测量或参考相关文献中的经验公式，确定齿轮的当量阻尼系数。在一些研究中，根据齿轮的工作条件和润滑情况，给出了计算齿轮阻尼系数的经验公式，可用于估算齿轮的当量阻尼。通过准确计算轴系各部件的当量参数，为建立完整、准确的轴系当量系统提供了坚实的数据基础，确保基于传递函数的轴系扭振数值计算能够更精确地反映实际轴系的扭振特性。3.1.3建立完整的轴系当量系统模型在完成轴系各部件当量参数的计算后，按照实际轴系的结构和连接方式，将各部件的当量模型连接起来，构建出完整的多缸内燃机轴系当量系统模型。以某直列四缸内燃机轴系为例，从自由端开始，依次连接硅油减振器、曲轴的四个曲拐、飞轮等部件的当量模型。硅油减振器可简化为一个由壳体惯量和惯性轮惯量组成的当量惯量，与曲轴前端的第一个曲拐通过具有一定当量刚度和当量阻尼的弹性轴段相连。曲轴的每个曲拐都作为一个集中质量，相邻曲拐之间通过当量刚度和当量阻尼的轴段连接，以模拟曲轴的实际弹性和阻尼特性。飞轮作为一个集中质量，连接在曲轴的后端，其与最后一个曲拐之间同样通过当量刚度和当量阻尼的轴段相连。在连接各部件当量模型时，需严格遵循实际轴系的结构和连接关系，确保当量系统模型能够准确反映实际轴系的动力学特性。对于存在齿轮传动的轴系，需根据齿轮的传动比，将从动齿轮的当量转动惯量折算到主动轴上，并正确设置齿轮啮合的当量刚度和当量阻尼。若轴系中还包含联轴器、离合器等部件，也需按照其实际连接位置和工作特性，将其当量模型准确地融入到整个轴系当量系统中。构建完成的轴系当量系统模型可以用数学表达式进行描述，通常采用动力学方程的形式。对于由n个集中质量组成的轴系当量系统，其动力学方程可表示为矩阵形式：[M]\{\ddot{\theta}\}+[C]\{\dot{\theta}\}+[K]\{\theta\}=\{T\}其中，[M]为质量矩阵，由各集中质量的当量转动惯量组成；[C]为阻尼矩阵，包含各弹性轴段和阻尼元件的当量阻尼系数；[K]为刚度矩阵，由各弹性轴段的当量刚度构成；\{\theta\}为各集中质量的扭角向量；\{T\}为作用在各集中质量上的激振力矩向量。通过建立这样的数学模型，可以方便地运用数值计算方法求解轴系的扭振响应。在实际计算中，可根据具体的数值计算方法（如有限元法、模态叠加法等），对上述动力学方程进行离散化处理，将连续的系统转化为离散的数值模型，以便在计算机上进行求解。利用有限元软件，将轴系当量系统划分为若干个有限元单元，通过对每个单元的动力学方程进行组装，得到整个轴系的有限元方程，然后运用数值算法求解该方程，得到轴系各节点的扭角、扭矩等扭振响应参数。完整的轴系当量系统模型的建立，为基于传递函数的多缸内燃机轴系扭振数值计算提供了核心的计算模型，使得后续能够通过数学分析和数值计算，深入研究轴系在各种工况下的扭振特性，为轴系的优化设计和故障诊断提供有力的工具。3.2传递函数模型的推导3.2.1基于动力学方程的传递函数推导在建立多缸内燃机轴系扭振的传递函数模型时，从轴系的动力学基本方程出发是关键的第一步。对于前文构建的轴系当量系统模型，其动力学方程通常可表示为二阶常微分方程组。以一个由n个集中质量组成的轴系为例，其动力学方程的矩阵形式为：[M]\{\ddot{\theta}\}+[C]\{\dot{\theta}\}+[K]\{\theta\}=\{T\}其中，[M]为质量矩阵，其对角元素为各集中质量的当量转动惯量J_i（i=1,2,\cdots,n），非对角元素为0，即[M]=diag(J_1,J_2,\cdots,J_n)；[C]为阻尼矩阵，其元素C_{ij}反映了第i个和第j个集中质量之间的阻尼关系；[K]为刚度矩阵，元素K_{ij}表示第i个和第j个集中质量之间的扭转刚度；\{\theta\}为各集中质量的扭角向量，\{\theta\}=\begin{bmatrix}\theta_1&\theta_2&\cdots&\theta_n\end{bmatrix}^T；\{T\}为作用在各集中质量上的激振力矩向量，\{T\}=\begin{bmatrix}T_1&T_2&\cdots&T_n\end{bmatrix}^T。为了推导传递函数，在零初始条件下，对上述动力学方程两边进行拉普拉斯变换。根据拉普拉斯变换的微分性质，L[\dot{\theta}_i(t)]=s\Theta_i(s)，L[\ddot{\theta}_i(t)]=s^{2}\Theta_i(s)（其中\Theta_i(s)为\theta_i(t)的拉普拉斯变换），可得：[M]s^{2}\{\Theta(s)\}+[C]s\{\Theta(s)\}+[K]\{\Theta(s)\}=\{T(s)\}进一步整理得到：\left([M]s^{2}+[C]s+[K]\right)\{\Theta(s)\}=\{T(s)\}则轴系扭振系统的传递函数矩阵[G(s)]可表示为：[G(s)]=\left([M]s^{2}+[C]s+[K]\right)^{-1}传递函数矩阵[G(s)]中的元素G_{ij}(s)表示第j个输入（激振力矩T_j(s)）对第i个输出（扭角\Theta_i(s)）的传递函数关系，即：\Theta_i(s)=G_{ij}(s)T_j(s)通过这样的推导，建立了多缸内燃机轴系扭振的传递函数模型，它将轴系的输入激振力矩与输出扭角在复频域中联系起来，为后续分析轴系的扭振特性提供了有力的数学工具。通过对传递函数矩阵的分析，可以求解轴系的固有频率、振型以及在不同激励下的扭振响应等重要参数，深入了解轴系的动态特性。3.2.2考虑非线性因素的传递函数修正在实际的多缸内燃机轴系中，存在多种非线性因素，如间隙、摩擦、非线性刚度等，这些因素会对轴系的扭振特性产生显著影响，因此需要对基于线性假设推导的传递函数进行修正，以更准确地描述轴系扭振特性。轴系中各部件之间的连接可能存在间隙，如轴承间隙、齿轮啮合间隙等。间隙的存在使得轴系在受力时，会出现一段无作用力的区域，只有当扭角超过一定值后，力才开始作用，这就导致了轴系的动力学特性呈现非线性。对于存在间隙的情况，可以采用分段线性化的方法对传递函数进行修正。将间隙区域和正常受力区域分别进行分析，建立相应的动力学方程。在间隙区域，由于没有力的作用，轴系的运动方程相对简单；而在正常受力区域，按照常规的线性动力学方程进行处理。然后通过适当的边界条件，将这两个区域的方程进行衔接，从而得到考虑间隙非线性的传递函数。轴系中的摩擦也是不可忽视的非线性因素，它主要包括轴与轴承之间的滑动摩擦、齿轮啮合时的摩擦等。摩擦的存在会消耗轴系扭振的能量，影响扭振的振幅和相位。常用的库仑摩擦模型可以表示为：F_f=\muN\text{sgn}(\dot{\theta})其中，F_f为摩擦力，\mu为摩擦系数，N为正压力，\text{sgn}(\dot{\theta})为符号函数，根据扭角速度\dot{\theta}的正负取值。在考虑摩擦非线性时，将摩擦力项加入到轴系的动力学方程中，然后对修正后的动力学方程进行拉普拉斯变换，得到考虑摩擦非线性的传递函数。在建立动力学方程时，需要根据实际情况准确确定摩擦力的方向和大小，以确保修正后的传递函数能够准确反映摩擦对轴系扭振的影响。一些轴系部件的刚度可能呈现非线性特性，如橡胶减振器等弹性元件，其刚度会随着变形的增大而发生变化。对于非线性刚度，通常采用非线性弹簧模型来描述，如幂函数形式的弹簧力-位移关系：F=kx^n其中，F为弹簧力，k为弹簧刚度系数，x为位移，n为非线性指数，n\neq1时表示非线性刚度。在考虑非线性刚度时，将非线性弹簧模型代入轴系的动力学方程，然后进行拉普拉斯变换，得到考虑非线性刚度的传递函数。在实际应用中，需要通过实验或其他方法准确确定非线性刚度的参数，以保证传递函数的准确性。通过对这些非线性因素的分析和相应的传递函数修正，可以使建立的传递函数模型更加符合多缸内燃机轴系的实际工作情况，提高轴系扭振数值计算的精度，为轴系的设计、优化和故障诊断提供更可靠的依据。在实际计算中，需要根据轴系的具体结构和工作条件，合理考虑各种非线性因素的影响，选择合适的修正方法，确保计算结果能够准确反映轴系的扭振特性。3.3数值计算方法与步骤3.3.1离散化处理与求解算法选择在基于传递函数的多缸内燃机轴系扭振数值计算中，首先需要对连续的轴系扭振系统进行离散化处理，将其转化为适合数值计算的离散模型。常见的离散化方法有有限元法、集中质量法等。有限元法是一种广泛应用的离散化方法，它将连续的轴系分割成有限个单元，通过对每个单元的分析和组装，得到整个轴系的数值模型。在轴系扭振分析中，通常将轴段划分为梁单元或杆单元，每个单元具有相应的节点和自由度。对于一个轴系，可根据其结构特点和精度要求，将轴段划分为若干个梁单元，每个梁单元的两端为节点，节点具有扭角和扭矩等自由度。通过对每个梁单元的动力学方程进行推导和组装，得到整个轴系的有限元方程，从而实现对轴系扭振的数值求解。有限元法的优点是能够精确地模拟轴系的复杂结构和边界条件，适用于各种形状和材料的轴系。其计算精度较高，能够得到轴系各节点的详细扭振信息，如扭角、扭矩、应力等。在分析具有变截面、阶梯轴等复杂结构的轴系时，有限元法能够准确地考虑轴段的几何形状和材料特性对扭振的影响。集中质量法是将轴系中的质量集中在若干个离散点上，这些离散点称为集中质量点，相邻集中质量点之间通过无质量的弹性轴段和阻尼元件连接，形成离散的轴系模型。在建立集中质量模型时，需要根据轴系各部件的质量分布和转动惯量，合理确定集中质量点的位置和大小。对于曲轴，可将每个曲拐的质量集中到曲拐的中心，作为一个集中质量点；将轴段的质量平均分配到相邻的集中质量点上。集中质量法的优点是计算简单、计算效率高，适用于对计算精度要求不是特别高的初步分析。在早期的轴系扭振研究中，集中质量法得到了广泛应用，能够快速地计算出轴系的固有频率和振型等基本扭振特性。在完成离散化处理后，需要选择合适的数值求解算法来求解轴系扭振的动力学方程。常用的数值求解算法有Runge-Kutta法、Newmark法等。Runge-Kutta法是一种基于数值积分的求解算法，具有较高的精度和稳定性。它通过在不同的时间步长上对微分方程进行近似求解，逐步得到系统的响应。对于轴系扭振的动力学方程，可将其转化为一阶微分方程组，然后应用Runge-Kutta法进行求解。Runge-Kutta法的优点是对不同类型的微分方程都具有较好的适应性，能够处理非线性、时变等复杂情况。在轴系扭振计算中，Runge-Kutta法能够准确地计算出轴系在不同工况下的扭振响应，特别是在处理瞬态响应时，表现出较好的性能。Newmark法是一种隐式的数值求解算法，主要用于求解结构动力学中的二阶微分方程，在轴系扭振分析中也有广泛应用。该方法通过对加速度和速度进行线性插值，建立时间步长内的运动方程，从而求解系统的响应。Newmark法的优点是计算精度较高，能够有效地处理高频振动问题，并且在一定程度上具有无条件稳定性，即对于较大的时间步长，仍然能够保证计算结果的稳定性。在多缸内燃机轴系扭振计算中，当轴系存在较高频率的振动成分时，Newmark法能够准确地捕捉到这些高频振动，得到较为精确的扭振响应结果。在实际应用中，需要根据轴系的具体特点和计算要求，综合考虑离散化方法和求解算法的选择。对于复杂结构的轴系，有限元法结合Runge-Kutta法或Newmark法能够提供高精度的计算结果，但计算量较大；对于简单轴系或初步分析，集中质量法结合较为简单的求解算法，如四阶Runge-Kutta法，能够在保证一定精度的前提下，提高计算效率。3.3.2计算流程与参数设置基于传递函数的轴系扭振数值计算流程包括输入参数设置、计算过程迭代和结果输出等关键环节。在输入参数设置环节，需要准确输入轴系的各种参数。轴系的几何参数是基础，包括轴径、轴长等，这些参数直接影响轴系的转动惯量和刚度。轴径的大小决定了轴的抗扭能力，轴长则影响轴系的整体刚度和振动特性。轴系各部件的材料参数也至关重要，如弹性模量、密度、泊松比等。弹性模量反映了材料抵抗弹性变形的能力，对轴系的刚度有直接影响；密度决定了部件的质量，进而影响转动惯量；泊松比则在计算轴系的应力和应变时发挥作用。轴系的阻尼系数也是不可忽视的参数，它反映了轴系在振动过程中能量耗散的特性。阻尼系数的大小会影响扭振的振幅和响应的衰减速度，准确确定阻尼系数对于得到准确的计算结果至关重要。除了这些物理参数，还需输入激励力矩的频谱特性，包括各阶谐波的幅值和相位。激励力矩是轴系扭振的驱动力，其频谱特性决定了轴系在不同频率下的响应情况。在计算过程迭代环节，首先根据输入的参数建立轴系的动力学方程。利用离散化方法将连续的轴系转化为离散模型，再结合传递函数理论，建立起描述轴系扭振的动力学方程。然后选择合适的数值求解算法对动力学方程进行求解。在求解过程中，需要设置合理的时间步长。时间步长的大小直接影响计算精度和计算效率。较小的时间步长可以提高计算精度，但会增加计算量和计算时间；较大的时间步长虽然可以提高计算效率，但可能会导致计算精度下降，甚至出现数值不稳定的情况。因此，需要根据轴系的振动特性和计算要求，通过试算或参考相关经验，确定合适的时间步长。在迭代计算过程中，还需进行收敛性判断。通过比较相邻迭代步的计算结果，如扭角、扭矩等参数的变化量，判断计算是否收敛。若变化量小于设定的收敛精度，则认为计算收敛，停止迭代；否则，继续进行迭代计算，直到满足收敛条件为止。在结果输出环节，将计算得到的轴系扭振响应结果进行整理和输出。输出的结果包括轴系各节点的扭角、扭矩随时间的变化曲线，这些曲线能够直观地展示轴系在不同时刻的扭振情况。通过分析扭角和扭矩的变化曲线，可以了解轴系的振动幅值、频率等信息，判断轴系是否存在共振等异常情况。还可以输出轴系的固有频率和振型等模态参数。固有频率是轴系的重要特征参数，它反映了轴系自身的振动特性，与激励频率的匹配情况会影响轴系的扭振响应。振型则描述了轴系在振动时各节点的相对位移关系，对于分析轴系的振动形态和薄弱环节具有重要意义。将这些结果以图表、数据文件等形式输出，以便后续对轴系扭振特性进行深入分析和研究。可将扭角和扭矩随时间的变化曲线绘制在坐标系中，清晰地展示其变化趋势；将固有频率和振型等参数整理成表格形式，方便对比和分析。四、案例分析与结果验证4.1某型号多缸内燃机轴系扭振实例计算4.1.1内燃机基本参数与工况设定本研究选取某型号直列六缸柴油机作为研究对象，该柴油机在船舶动力系统中广泛应用，其结构紧凑、功率强劲，能够满足船舶在不同航行工况下的动力需求。其主要结构参数如表1所示：[此处插入表格1：某型号直列六缸柴油机主要结构参数][此处插入表格1：某型号直列六缸柴油机主要结构参数]参数名称数值气缸数6缸径120mm行程130mm压缩比17.5额定转速1800r/min额定功率200kW曲轴材料40Cr弹性模量2.1×10¹¹Pa剪切模量8.1×10¹⁰Pa该柴油机采用四冲程工作循环，在一个工作循环中，曲轴旋转两圈，每个气缸依次完成进气、压缩、做功和排气四个冲程。进气冲程中，活塞下行，进气门打开，新鲜空气被吸入气缸；压缩冲程中，活塞上行，进气门和排气门关闭，空气被压缩；做功冲程中，喷油器向气缸内喷射燃油，燃油与高温高压的空气混合燃烧，产生强大的气体压力，推动活塞下行做功；排气冲程中，活塞上行，排气门打开，燃烧后的废气被排出气缸。在实际运行中，该柴油机存在多种工况，为全面研究其轴系扭振特性，本研究设定了如表2所示的三种典型工况：[此处插入表格2：某型号直列六缸柴油机典型工况设定][此处插入表格2：某型号直列六缸柴油机典型工况设定]工况转速(r/min)负荷(%)工况1120050工况2150075工况31800100在工况1下，柴油机处于中低转速、部分负荷运行状态，此时船舶可能在低速航行或进行一些轻载作业。工况2代表了柴油机在中等转速、中等负荷下的运行情况，是船舶正常航行时常见的工况。工况3为柴油机的额定工况，在该工况下，柴油机以额定转速和满负荷运行，输出最大功率，通常用于船舶加速、顶风航行或重载作业等情况。通过对这三种典型工况的研究，可以较为全面地了解该柴油机轴系在不同运行条件下的扭振特性。4.1.2基于传递函数的数值计算过程根据前文建立的轴系当量系统模型和传递函数模型，对选定的某型号直列六缸柴油机轴系在三种典型工况下的扭振进行数值计算。首先，根据表1中的结构参数，计算轴系各部件的当量参数。对于曲轴，将每个曲拐的质量等效集中到曲拐中心，根据转动惯量计算公式，计算出每个曲拐的当量转动惯量。已知曲拐的质量为m_i，质心到曲轴轴线的距离为r_i，则曲拐的当量转动惯量J_{crank,i}=m_ir_i^2。对于该六缸柴油机，六个曲拐的当量转动惯量分别为J_{crank,1}=0.12kg·m²、J_{crank,2}=0.12kg·m²、J_{crank,3}=0.12kg·m²、J_{crank,4}=0.12kg·m²、J_{crank,5}=0.12kg·m²、J_{crank,6}=0.12kg·m²。相邻曲拐之间的轴段当量刚度根据材料力学中的扭转刚度公式k=\frac{GJ_p}{l}计算，其中G为材料的剪切模量，J_p为轴的极惯性矩，l为轴段长度。经计算，相邻曲拐间的当量刚度k_{i,i+1}=3.5×10^7N·m/rad（i=1,2,3,4,5）。飞轮的当量转动惯量根据其实际质量和回转半径计算，已知飞轮质量m_{flywheel}=50kg，回转半径r_{flywheel}=0.3m，则飞轮的当量转动惯量J_{flywheel}=m_{flywheel}r_{flywheel}^2=50×0.3²=4.5kg·m²。飞轮与曲轴后端连接轴段的当量刚度同样根据上述扭转刚度公式计算，为k_{6,flywheel}=4.0×10^7N·m/rad。在计算阻尼系数时，考虑到轴系中存在多种阻尼因素，通过参考相关经验数据和实验测量，确定轴系的当量阻尼系数c_{i,i+1}=500N·m·s/rad（i=1,2,3,4,5），c_{6,flywheel}=600N·m·s/rad。根据轴系各部件的当量参数，建立轴系的动力学方程：[M]\{\ddot{\theta}\}+[C]\{\dot{\theta}\}+[K]\{\theta\}=\{T\}其中，质量矩阵[M]为：[M]=\begin{bmatrix}J_{crank,1}&0&0&0&0&0&0\\0&J_{crank,2}&0&0&0&0&0\\0&0&J_{crank,3}&0&0&0&0\\0&0&0&J_{crank,4}&0&0&0\\0&0&0&0&J_{crank,5}&0&0\\0&0&0&0&0&J_{crank,6}&0\\0&0&0&0&0&0&J_{flywheel}\end{bmatrix}阻尼矩阵[C]为：[C]=\begin{bmatrix}c_{1,2}&-c_{1,2}&0&0&0&0&0\\-c_{1,2}&c_{1,2}+c_{2,3}&-c_{2,3}&0&0&0&0\\0&-c_{2,3}&c_{2,3}+c_{3,4}&-c_{3,4}&0&0&0\\0&0&-c_{3,4}&c_{3,4}+c_{4,5}&-c_{4,5}&0&0\\0&0&0&-c_{4,5}&c_{4,5}+c_{5,6}&-c_{5,6}&0\\0&0&0&0&-c_{5,6}&c_{5,6}+c_{6,flywheel}&-c_{6,flywheel}\\0&0&0&0&0&-c_{6,flywheel}&c_{6,flywheel}\end{bmatrix}刚度矩阵[K]为：[K]=\begin{bmatrix}k_{1,2}&-k_{1,2}&0&0&0&0&0\\-k_{1,2}&k_{1,2}+k_{2,3}&-k_{2,3}&0&0&0&0\\0&-k_{2,3}&k_{2,3}+k_{3,4}&-k_{3,4}&0&0&0\\0&0&-k_{3,4}&k_{3,4}+k_{4,5}&-k_{4,5}&0&0\\0&0&0&-k_{4,5}&k_{4,5}+k_{5,6}&-k_{5,6}&0\\0&0&0&0&-k_{5,6}&k_{5,6}+k_{6,flywheel}&-k_{6,flywheel}\\0&0&0&0&0&-k_{6,flywheel}&k_{6,flywheel}\end{bmatrix}\{\theta\}为各集中质量的扭角向量，\{T\}为作用在各集中质量上的激振力矩向量。在零初始条件下，对动力学方程两边进行拉普拉斯变换，得到：\left([M]s^{2}+[C]s+[K]\right)\{\Theta(s)\}=\{T(s)\}则轴系扭振系统的传递函数矩阵[G(s)]为：[G(s)]=\left([M]s^{2}+[C]s+[K]\right)^{-1}对于不同工况，根据柴油机的工作循环和负荷情况，确定作用在轴系上的激振力矩。在四冲程柴油机中，每个气缸的做功冲程产生的激振力矩是轴系扭振的主要激励源。激振力矩的大小和相位与气缸的工作顺序、负荷以及转速有关。以工况1（转速1200r/min，负荷50%）为例，通过查阅相关资料和理论计算，得到各气缸在一个工作循环内的激振力矩随曲轴转角的变化规律，将其进行傅里叶分解，得到各阶谐波激振力矩的幅值和相位。将激振力矩的拉普拉斯变换\{T(s)\}与传递函数矩阵[G(s)]相乘，得到轴系各节点扭角的拉普拉斯变换\{\Theta(s)\}：\{\Theta(s)\}=[G(s)]\{T(s)\}再对\{\Theta(s)\}进行拉普拉斯反变换，得到轴系各节点在时域中的扭角响应\theta(t)。在数值计算过程中，采用四阶Runge-Kutta法进行求解，设置时间步长\Deltat=1×10^{-4}s，以确保计算精度和稳定性。通过迭代计算，得到各工况下轴系各节点的扭角随时间的变化曲线，以及轴系的固有频率和振型等关键数据。4.2实验测试与数据采集4.2.1实验台架搭建与测试方案设计为了验证基于传递函数的多缸内燃机轴系扭振数值计算方法的准确性，搭建了专门的实验台架，其结构示意图如图2所示。实验台架主要由某型号直列六缸柴油机、测功机、扭矩传感器、转速传感器、扭振传感器、信号采集系统和计算机等部分组成。[此处插入实验台架结构示意图][此处插入实验台架结构示意图]柴油机作为实验的核心部件，为轴系提供动力，使其产生扭振。测功机用于模拟不同的负载工况，通过调节测功机的加载电阻，改变负载大小，从而实现对柴油机不同工况的模拟。在本实验中，根据前文设定的三种典型工况，通过测功机分别加载相应的负载，使柴油机在1200r/min、50%负荷，1500r/min、75%负荷和1800r/min、100%负荷工况下稳定运行。扭矩传感器安装在柴油机与测功机之间的连接轴上，用于实时测量轴系传递的扭矩，其测量精度为±0.5%FS，能够准确捕捉扭矩的变化。转速传感器采用磁电式传感器，安装在曲轴自由端附近，通过测量曲轴的转速，为实验提供转速数据，其测量精度为±1r/min。扭振传感器是实验的关键部件之一，选用高精度的应变片式扭振传感器。在轴系的关键位置，如曲轴自由端、各曲拐之间的轴段以及飞轮处，共布置了7个扭振传感器。曲轴自由端的传感器用于测量轴系自由端的扭振情况，反映轴系扭振的起始状态；各曲拐之间轴段的传感器能够监测不同轴段的扭振响应，分析扭振在轴系中的传播特性；飞轮处的传感器则用于测量飞轮的扭振情况，了解飞轮对轴系扭振的影响。每个扭振传感器的测量精度为±0.01°，能够满足实验对扭振测量精度的要求。信号采集系统负责采集扭矩传感器、转速传感器和扭振传感器输出的信号。采用多通道数据采集卡，其采样频率最高可达100kHz，能够满足对高频扭振信号的采集需求。采集卡将模拟信号转换为数字信号，并传输至计算机进行后续处理。在测试方案中，针对不同工况，分别进行多次测量，每次测量持续时间为30s，以获取稳定的实验数据。在每个工况下，测量前先让柴油机在该工况下稳定运行5min，使轴系达到稳定的扭振状态，然后开始采集数据。在数据采集过程中，同时记录柴油机的转速、扭矩、油温、水温等参数，以便对实验结果进行全面分析。4.2.2实验数据采集与处理在完成实验台架搭建和测试方案设计后，按照设定的三种典型工况，即工况1（1200r/min、50%负荷）、工况2（1500r/min、75%负荷）和工况3（1800r/min、100%负荷），依次进行实验。在每个工况下，启动柴油机，待其稳定运行后，通过信号采集系统同步采集轴系各测点的扭振信号、扭矩信号和转速信号。实验过程中，由于受到各种干扰因素的影响，采集到的原始信号中往往包含噪声和其他无用信息，因此需要对原始信号进行处理，以提高信号质量，为后续分析提供准确的数据。采用低通滤波器对原始信号进行滤波处理，去除高频噪声。根据轴系扭振的频率特性，设置低通滤波器的截止频率为1000Hz，该频率能够有效滤除高频噪声，同时保留轴系扭振的主要频率成分。利用小波降噪算法对信号进行进一步降噪处理。小波降噪算法能够根据信号的局部特征，自适应地调整降噪参数，有效去除信号中的噪声，同时保留信号的细节信息。通过多次实验对比，选择合适的小波基函数和分解层数，以达到最佳的降噪效果。在本实验中，选择db4小波基函数，分解层数为5，经过小波降噪处理后，信号的信噪比得到显著提高，能够更清晰地反映轴系扭振的真实情况。经过滤波和降噪处理后，对信号进行特征提取。计算轴系各测点的扭振幅值、频率和相位等特征参数。扭振幅值反映了轴系扭振的剧烈程度，通过对各测点扭振幅值的分析，可以确定轴系扭振的主要部位和最大扭振幅值出现的位置。频率是轴系扭振的重要特征参数之一，通过傅里叶变换等方法，分析扭振信号的频率成分，确定轴系的固有频率和主要激励频率，判断轴系是否存在共振现象。相位信息则反映了各测点扭振的相对位置关系，对于分析轴系的振型和振动传播特性具有重要意义。以工况2为例，经过处理后得到的曲轴自由端扭振幅值随时间的变化曲线如图3所示。从图中可以看出，在该工况下，曲轴自由端的扭振幅值在一定范围内波动，最大值约为0.3°，最小值约为0.1°，平均扭振幅值为0.2°左右。通过对扭振信号的频谱分析，得到其主要频率成分集中在100Hz-300Hz之间，其中200Hz左右的频率成分幅值较大，与该工况下柴油机的工作频率相匹配，说明该频率成分是由柴油机的工作循环引起的主要激励频率。[此处插入工况2下曲轴自由端扭振幅值随时间变化曲线][此处插入工况2下曲轴自由端扭振幅值随时间变化曲线]通过对不同工况下实验数据的采集和处理，得到了轴系在各种工况下的扭振特性，为后续与数值计算结果的对比分析提供了实验依据，有助于验证基于传递函数的多缸内燃机轴系扭振数值计算方法的准确性和可靠性。4.3计算结果与实验结果对比分析4.3.1对比分析方法与指标选择为了全面、准确地评估基于传递函数的多缸内燃机轴系扭振数值计算方法的准确性和可靠性，选择了误差分析和频谱分析作为主要的对比分析方法，并确定了扭振振幅、频率等关键对比指标。误差分析是评估计算结果与实验结果偏差程度的重要手段。通过计算扭振振幅和频率的绝对误差和相对误差，能够直观地了解计算值与实验值之间的差异大小。绝对误差直接反映了计算值与实验值之间的差值，计算公式为：E_{abs}=|x_{cal}-x_{exp}|，其中E_{abs}为绝对误差，x_{cal}为计算值，x_{exp}为实验值。相对误差则是绝对误差与实验值的比值，以百分比的形式表示，更能体现误差的相对大小，计算公式为：E_{rel}=\frac{|x_{cal}-x_{exp}|}{x_{exp}}\times100\%。在对比轴系各节点的扭振振幅时，通过计算绝对误差和相对误差，可以清晰地判断计算结果在振幅大小上与实验结果的偏离程度，从而评估计算方法的精度。频谱分析是从频率域的角度对计算结果和实验结果进行对比。通过对扭振信号进行傅里叶变换，将时域信号转换为频域信号，得到扭振信号的频谱图。在频谱图中，能够直观地观察到扭振信号的主要频率成分及其幅值分布。对比计算结果和实验结果的频谱图，可以判断计算方法是否能够准确捕捉到轴系扭振的主要频率成分，以及各频率成分的幅值是否与实验结果相符。若计算结果的频谱图中，主要频率成分的位置和幅值与实验结果的频谱图一致，则说明计算方法能够准确反映轴系扭振的频率特性；反之，则说明计算方法存在一定的误差。扭振振幅是衡量轴系扭振剧烈程度的重要指标，直接影响轴系的安全性和可靠性。在对比分析中，重点关注轴系关键部位，如曲轴自由端、各曲拐之间的轴段以及飞轮处的扭振幅值。通过对比计算得到的扭振幅值与实验测量值，可以评估计算方法在预测轴系扭振剧烈程度方面的准确性。若计算得到的扭振幅值与实验值接近，则说明计算方法能够较好地预测轴系的扭振情况；若两者差异较大，则需要进一步分析原因，改进计算方法。频率是轴系扭振的另一个关键指标，它反映了轴系扭振的周期性和振动特性。在对比分析中，主要对比轴系的固有频率和主要激励频率。固有频率是轴系自身的特性参数，与轴系的结构、质量和刚度等因素有关。通过对比计算得到的固有频率与实验测量值，可以验证计算方法对轴系固有特性的描述是否准确。主要激励频率则与内燃机的工作循环和工况有关，对比计算结果和实验结果中主要激励频率的一致性，能够判断计算方法是否能够准确反映内燃机工作过程中对轴系的激励作用。4.3.2结果对比与误差讨论将某型号直列六缸柴油机轴系扭振的数值计算结果与实验测试结果进行详细对比，并对两者之间的差异和误差来源进行深入分析，以评估基于传递函数的数值计算方法的准确性和可靠性。在扭振振幅方面，以工况2（1500r/min、75%负荷）为例，实验测量得到的曲轴自由端最大扭振幅值为0.32°，而数值计算得到的结果为0.30°。通过计算绝对误差E_{abs}=|0.30-0.32|=0.02°，相对误差E_{rel}=\frac{|0.30-0.32|}{0.32}\times100\%\approx6.25\%。在其他工况下，各测点的扭振幅值计算结果与实验结果也存在一定差异，具体数据如表3所示：[此处插入表格3：不同工况下各测点扭振幅值对比（单位：°）][此处插入表格3：不同工况下各测点扭振幅值对比（单位：°）]工况测点实验值计算值绝对误差相对误差（%）工况1曲轴自由端0.210.190.029.52工况1曲拐1-曲拐2间轴段0.150.130.0213.33工况2曲拐3-曲拐4间轴段0.250.230.028.00工况3飞轮处0.400.370.037.50从表中数据可以看出，在不同工况下，数值计算得到的扭振幅值与实验值的绝对误差在0.02°-0.03°之间，相对误差在6.25%-13.33%之间。扭振幅值的差异可能由多种因素导致。在建立轴系当量系统模型时，虽然遵循了转动惯量、刚度和阻尼等效原则，但实际轴系的结构和参数存在一定的不确定性，如轴系各部件的加工精度、装配间隙等，这些因素可能导致当量参数与实际值存在偏差，从而影响扭振计算结果。实验测量过程中也存在一定的误差，包括传感器的测量误差、信号传输过程中的干扰以及数据采集系统的精度限制等，这些误差也会使实验测量值与真实值存在差异。在频率方面，通过频谱分析得到实验结果和数值计算结果的主要频率成分。以工况3（1800r/min、100%负荷）为例，实验测得的轴系主要激励频率为250Hz，对应柴油机的工作频率；数值计算得到的主要激励频率为248Hz。两者频率差值为\Deltaf=250-248=2Hz，相对误差为\frac{2}{250}\times100\%=0.8\%。在固有频率方面，实验测量得到轴系的一阶固有频率为120Hz，数值计算结果为118Hz，频率差值为\Deltaf=120-118=2Hz，相对误差为\frac{2}{120}\times100\%\approx1.67\%。不同工况下轴系频率对比数据如表4所示：[此处插入表格4：不同工况下轴系频率对比（单位：Hz）][此处插入表格4：不同工况下轴系频率对比（单位：Hz）]工况主要激励频率（实验值）主要激励频率（计算值）频率差值相对误差（%）一阶固有频率（实验值）一阶固有频率（计算值）频率差值相对误差（%）工况117517321.14908822.22工况221020820.9510510321.90工况325024820.8012011821.67从频率对比结果来看，数值计算得到的主要激励频率和一阶固有频率与实验值较为接近，频率差值在2Hz左右，相对误差在0.80%-2.22%之间。频率计算结果与实验值存在差异的原因主要包括：在建立传递函数模型时，虽然考虑了轴系的主要结构和参数，但对于一些细微的结构变化和材料特性的不均匀性等因素可能考虑不足，这些因素会对轴系的固有频率产生一定影响；实验测量过程中，由于受到外界干扰和测量设备精度的限制，可能导致测量得到的频率存在一定误差。综合扭振振幅和频率的对比结果，基于传递函数的多缸内燃机轴系扭振数值计算方法能够较好地预测轴系的扭振特性，计算结果与实验结果在趋势上基本一致，主要参数的误差在可接受范围内。该计算方法仍存在一定的改进空间，需要进一步优化轴系当量系统模型和传递函数模型，考虑更多的实际因素，以提高计算精度。在实验方面，也需要进一步提高测量设备的精度，减少测量误差，为数值计算提供更准确的实验验证数据。五、计算方法的优化与改进5.1算法优化策略5.1.1提高计算精度的方法在基于传递函数的多缸内燃机轴系扭振数值计算中，提高计算精度是关键目标之一，可通过增加计算步长、改进数值算法、优化离散化方案等多种方法来实现。适当减小计算步长能够显著提高计算精度。在数值计算中，计算步长的大小直接影响着计算结果的准确性。以四阶Runge-Kutta法为例，其计算过程基于泰勒级数展开，步长越小，泰勒级数展开的截断误差就越小，计算结果就越接近真实值。对于一些对精度要求极高的轴系扭振计算，如航空发动机轴系，由于其工作转速高、振动特性复杂，采用较小的计算步长能够更精确地捕捉轴系扭振的细节变化，避免因步长过大而导致的计算误差积累，从而提高计算精度。减小计算步长会显著增加计算量和计算时间，因此需要在精度和效率之间进行权衡。可根据轴系的具体振动特性和计算要求，通过试算确定合适的计算步长。在计算初期，可采用较大的步长进行初步计算，了解轴系扭振的大致趋势，然后逐步减小步长，对关键部位和重要频率段进行精细计算，以在保证精度的前提下提高计算效率。改进数值算法是提高计算精度的重要途径。传统的数值算法在处理复杂的轴系扭振问题时，可能存在精度不足的问题。例如，对于具有强非线性特性的轴系，传统的线性数值算法可能无法准确描述轴系的动态行为。此时，可采用自适应变步长算法，如Runge-Kutta-