基于超单元法革新总振动经验公式：理论、实践与展望

基于超单元法革新总振动经验公式：理论、实践与展望一、引言1.1研究背景与动机在全球贸易持续繁荣的背景下，船舶行业作为国际贸易的重要载体，正处于快速发展与变革之中。随着科技的进步以及人们对船舶性能要求的不断提高，船型日益呈现出多样化的发展趋势。从大型集装箱船到特种作业船，从豪华邮轮再到高速双体船，不同类型的船舶在设计、功能和应用场景上都有显著差异。例如，大型集装箱船为了提高运输效率，不断增大尺寸和载箱量，其设计需要考虑如何在保证结构强度的前提下，优化水动力性能以降低能耗；而特种作业船如海洋科考船、海上风电安装船等，则对船舶的稳定性、定位精度以及作业设备的搭载能力有特殊要求。在船舶设计领域，总振动特性是衡量船舶性能的关键指标之一。船舶在航行过程中会受到各种动态载荷的作用，如波浪力、螺旋桨激振力等，这些力会引发船舶的总振动。过大的振动不仅会影响船员和乘客的舒适性，还可能导致结构疲劳损坏，降低船舶的使用寿命和安全性。因此，准确预测船舶的总振动特性，对于船舶的设计、建造和运营具有重要意义。传统上，船舶总振动特性的预测主要依赖于经验公式。这些经验公式是基于大量的实际船舶数据和试验结果，通过统计回归方法得到的。它们在一定程度上能够反映船舶总振动特性与船型主尺度、结构参数等因素之间的关系。例如，希列克法、陶德法和CCS法等都是常用的经验公式，在船舶设计的初步阶段为工程师提供了重要的参考依据。然而，随着船型的多样化发展，传统经验公式逐渐暴露出其局限性。一方面，这些公式大多是基于特定类型的船舶和有限的样本数据建立的，对于新型船型或特殊船型，其适用性和准确性大打折扣。例如，对于三体船、小水线面双体船等新型多体船型，由于其独特的结构形式和水动力特性，传统经验公式难以准确预测其总振动特性。另一方面，经验公式的回归过程通常需要耗费大量的时间和精力，且对于一些小众船型，由于数据获取困难，很难建立有效的经验公式。为了解决传统经验公式在新船型设计中的不足，引入超单元法对总振动经验公式进行修正具有重要的现实意义。超单元法作为一种先进的数值分析方法，基于动态子结构法的理念，将复杂的结构模型分解为多个超单元和残余结构。通过对超单元进行独立的分析和计算，然后再将其与残余结构进行组合求解，能够有效地提高计算效率和精度。在船舶总振动特性分析中，超单元法具有诸多优势。首先，它可以利用三维模型进行计算，使结构分析更加接近实际情况，从而提高预测的准确性。其次，超单元具有易修改性，能够方便地获取不同主尺度下船型的数据。通过对这些数据进行回归分析，可以得到更准确、更具针对性的某特定船型的总振动预报公式。这对于指导新型船型的设计，优化船舶结构，提高船舶的综合性能具有重要的指导作用。1.2研究目的与意义本研究旨在引入超单元法，对传统的船舶总振动经验公式进行修正，以提高对新型船型和特殊船型总振动特性的预测准确性。具体而言，研究目标包括以下几个方面：首先，深入研究超单元法的理论和应用，掌握其在船舶结构分析中的优势和适用范围。通过对超单元法的深入理解，能够更好地将其应用于船舶总振动特性的分析中，为修正经验公式提供理论基础。其次，利用超单元法对特定船型进行数据收集和分析，获取不同主尺度下船型的振动特性数据。这些数据将为后续的回归分析提供丰富的样本，从而建立起更准确的总振动预报公式。最后，通过对比修正前后经验公式的预测结果，验证超单元法修正经验公式的有效性和准确性。本研究的意义主要体现在以下几个方面：在船舶设计领域，准确的总振动预报公式对于新型船型的设计和优化具有重要指导意义。通过修正经验公式，可以更准确地预测船舶在不同工况下的总振动特性，为船舶结构设计提供更可靠的依据，从而提高船舶的安全性和舒适性。在工程实践中，利用超单元法修正经验公式能够缩短设计周期，降低设计成本。传统的经验公式回归过程需要耗费大量的时间和精力，而超单元法可以快速获取不同主尺度下船型的数据，通过回归分析得到更准确的经验公式，为工程实践提供更高效的解决方案。此外，本研究对于推动船舶总振动理论的发展也具有重要意义。超单元法的引入为船舶总振动特性的研究提供了新的思路和方法，有助于完善船舶总振动理论体系，促进船舶工程学科的发展。1.3研究方法与创新点在本研究中，采用了多种研究方法来确保研究的科学性和有效性。理论分析是基础，通过深入研究船舶总振动理论以及超单元法的基本原理，从理论层面揭示船舶总振动特性与船型主尺度、结构参数之间的内在联系。例如，详细剖析超单元法中动态子结构法的理念，明确其将复杂结构模型分解为超单元和残余结构进行分析的优势，为后续的研究提供坚实的理论支撑。案例研究法是重要手段，选取特定的船型作为研究对象，对其进行深入的分析和研究。通过对该船型建立三维有限元模型，并运用超单元法进行计算，获取不同主尺度下船型的振动特性数据。以某新型三体船为例，利用超单元法对其进行结构分析，得到在不同航行工况下的总振动特性数据，为后续的回归分析提供了丰富的样本。数据回归分析是关键环节，运用统计学方法对通过超单元法获取的数据进行回归分析，建立总振动预报公式。采用最小二乘法等回归方法，对数据进行拟合，找出船型主尺度、结构参数与总振动特性之间的函数关系。通过对大量数据的回归分析，得到针对该特定船型的总振动预报公式，并对公式的准确性和可靠性进行验证。本研究的创新点主要体现在研究方法和应用领域方面。在研究方法上，引入超单元法对船舶总振动经验公式进行修正，这是一种创新的尝试。超单元法基于动态子结构法，能够利用三维模型进行计算，使结构分析更加接近实际情况，有效提高了计算效率和精度。同时，超单元具有易修改性，可以方便地获取不同主尺度下船型的数据，为建立更准确的经验公式提供了数据支持。与传统的经验公式回归方法相比，超单元法大大缩短了回归周期，提高了经验公式的准确性和适用性。在应用领域方面，本研究针对新型船型和特殊船型进行研究，弥补了传统经验公式在这些领域的不足。随着船型的多样化发展，传统经验公式难以准确预测新型船型和特殊船型的总振动特性。本研究通过引入超单元法，为这些新型船型和特殊船型的总振动特性分析提供了新的方法和思路，有助于推动船舶设计领域的发展，提高船舶的综合性能。二、总振动经验公式与超单元法理论2.1总振动经验公式概述2.1.1常见总振动经验公式介绍在船舶总振动特性的研究中，经过长期的实践与探索，涌现出了一系列经典的经验公式，其中希列克法、陶德法和CCS法具有代表性，它们在船舶工程领域的发展历程中占据着重要的地位。希列克法是早期用于预测船舶总振动特性的重要方法之一。该方法的基本原理基于对船舶结构的简化假设，将船舶视为一个等效的梁结构，通过对梁的振动理论进行推导和应用，建立起船舶总振动特性与船型主尺度之间的关系。具体而言，希列克法认为船舶的总振动主要由船体梁的弯曲振动引起，通过考虑船体梁的惯性矩、长度、材料弹性模量等因素，来计算船舶的固有频率。其计算公式通常为：f=C\times\sqrt{\frac{EI}{L^3}}，其中f表示固有频率，C为经验系数，E是材料弹性模量，I为船体梁的惯性矩，L为船长。希列克法主要适用于常规船型的初步设计阶段，在这个阶段，设计人员对船舶的具体结构细节了解有限，希列克法能够凭借其简单的计算方式，快速地给出船舶总振动特性的大致估算，为后续的设计工作提供初步的参考。然而，希列克法也存在明显的局限性。由于其对船舶结构的简化过于理想化，忽略了许多实际因素的影响，如船体的局部结构特性、附件的影响等，导致在实际应用中，对于一些复杂船型或对振动特性要求较高的船舶，希列克法的计算结果与实际情况存在较大偏差。陶德法在船舶总振动特性预测方面具有独特的理论和应用价值。陶德法充分考虑了船体结构的复杂性以及各种实际因素对振动的影响，其原理基于对船舶结构的更细致分析，不仅仅局限于将船舶视为简单的梁结构，还综合考虑了船体的质量分布、刚度分布以及阻尼特性等因素。在计算过程中，陶德法通过引入一系列修正系数和经验公式，对船舶的总振动特性进行更为精确的计算。例如，在考虑质量分布时，陶德法会详细分析船舶各个部分的质量，包括船体结构、设备、货物等，并根据它们在船体中的位置和分布情况，合理地确定其对总振动的影响权重。陶德法适用于各种船型的振动特性分析，尤其在对船舶振动特性要求较高的场合，如豪华邮轮、高速客船等，陶德法能够提供相对准确的预测结果，为船舶的设计和优化提供有力的支持。不过，陶德法的计算过程较为复杂，需要较多的船舶结构参数和经验数据，这在一定程度上限制了其在实际工程中的广泛应用。而且，由于经验数据的局限性，对于一些新型船型或特殊船型，陶德法的计算精度也可能受到影响。CCS法是我国船舶行业中广泛应用的一种船舶总振动经验公式，它是在充分考虑我国船舶设计和建造实际情况的基础上，结合大量的试验数据和工程实践经验而建立起来的。CCS法涵盖了船舶垂向振动和水平振动的经验公式，对于船舶的设计和评估具有重要的指导意义。在垂向振动方面，CCS法通过对大量船舶的垂向振动数据进行统计分析，建立了垂向振动固有频率与船型主尺度、结构参数之间的经验关系。其垂向振动经验公式一般形式为：f_{v}=a\times(\frac{L}{B})^{b}\times(\frac{d}{T})^{c}\times\sqrt{\frac{\Delta}{L^3}}，其中f_{v}表示垂向振动固有频率，a、b、c为经验系数，L为船长，B为型宽，d为型深，T为吃水，\Delta为排水量。在水平振动方面，CCS法也有相应的经验公式，通过考虑船舶的横向结构特性、水动力作用等因素，对水平振动特性进行预测。CCS法在我国船舶工程领域具有广泛的应用，对于常规船型的振动特性预测具有较高的准确性，能够满足我国船舶设计和建造的实际需求。然而，随着船舶技术的不断发展和新型船型的不断涌现，CCS法在一些特殊船型的应用中也暴露出精度不足的问题，需要进一步的改进和完善。2.1.2总振动经验公式的应用与局限在船舶设计与建造的漫长历史进程中，总振动经验公式始终扮演着举足轻重的角色，为船舶工程领域的发展提供了不可或缺的支持。在船舶设计的初步阶段，设计人员往往面临着诸多不确定性和未知因素，此时总振动经验公式成为了他们快速获取船舶总振动特性大致信息的重要工具。例如，在设计一艘新的集装箱船时，设计人员可以利用希列克法或CCS法等经验公式，根据船型主尺度等基本参数，迅速估算出船舶的固有频率和振动响应。这些估算结果能够帮助设计人员初步判断船舶在不同工况下的振动情况，为后续的详细设计提供重要的参考依据。在实际应用中，经验公式还可以用于对现有船舶的振动性能进行评估。通过将实际船舶的参数代入经验公式中，计算出其理论振动特性，并与实际测量结果进行对比，从而判断船舶的振动性能是否符合要求。如果发现计算结果与实际测量值存在较大偏差，就可以进一步分析原因，采取相应的措施进行改进，如优化船舶结构、调整设备布局等，以提高船舶的振动性能。然而，随着船舶技术的飞速发展，船型日益多样化，传统的总振动经验公式逐渐暴露出其固有的局限性。在面对一些特殊船型时，经验公式的精度明显不足。以三体船为例，三体船具有独特的多体结构形式，其水动力特性和结构动力学特性与常规单体船有很大的差异。传统的经验公式大多是基于单体船的试验数据和经验建立起来的，对于三体船的复杂结构和特殊的水动力作用考虑不足，因此在预测三体船的总振动特性时，往往会出现较大的误差。同样，对于小水线面双体船等新型船型，由于其水线面面积较小，在波浪中的运动响应与常规船型不同，传统经验公式也难以准确预测其总振动特性。此外，一些具有特殊用途的船舶，如海洋科考船、海上风电安装船等，其结构和设备布局往往具有特殊性，对振动的要求也各不相同，传统经验公式在这些特殊船型上的应用也存在一定的困难。经验公式的另一个局限性在于其回归周期较长。建立一个可靠的经验公式需要大量的试验数据和工程实践经验作为支撑，这需要耗费大量的时间和精力。在实际操作中，首先需要对不同船型、不同工况下的船舶进行大量的振动试验，测量其振动响应数据。然后，对这些数据进行整理、分析和统计，找出其中的规律和趋势。最后，通过数学建模和回归分析等方法，建立起经验公式。这个过程往往需要数年甚至数十年的时间，而且随着船舶技术的不断发展，新的船型和工况不断出现，原有的经验公式可能不再适用，需要重新进行回归分析和修正，这进一步增加了经验公式的更新难度和时间成本。对于一些新型船型或小众船型，由于其数量较少，获取足够的试验数据变得更加困难，这使得建立有效的经验公式几乎成为不可能。在这种情况下，传统经验公式就无法为这些船型的设计和评估提供准确的指导，限制了船舶技术的创新和发展。2.2超单元法基本理论2.2.1超单元法的定义与原理超单元法作为有限元分析领域中一种高效且强大的数值分析方法，在解决复杂结构动力学问题时展现出独特的优势。从本质上讲，超单元法基于动态子结构法的基本理念，将复杂的结构模型分解为多个超单元和残余结构，通过对各个部分进行独立分析，再将其结果进行组合，从而实现对整体结构的全面分析。这种方法的核心在于将大型、复杂的结构问题转化为相对简单的子问题，大大降低了计算的复杂性和难度。超单元法的基本原理可通过以下方式深入理解：首先，将复杂结构按照一定的规则和方法划分为多个超单元和残余结构。超单元是结构中具有相对独立性和完整性的部分，它可以是结构的一个部件、一个子系统或一个特定的区域。例如，在船舶结构中，可将船体的不同舱段、上层建筑、推进系统等分别视为超单元。划分超单元的原则通常基于结构的几何形状、功能特点、连接方式以及计算的便利性等因素。划分完成后，对每个超单元进行独立的分析和计算。在这个过程中，需要确定超单元的边界条件，边界条件的确定至关重要，它直接影响到超单元分析结果的准确性和可靠性。通常，边界条件包括位移边界条件、力边界条件以及速度边界条件等，这些条件的设定需根据结构的实际工作状态和受力情况进行合理的选择和确定。确定边界条件后，利用有限元方法对超单元进行分析，得到超单元的刚度矩阵、质量矩阵和阻尼矩阵等重要参数。这些矩阵反映了超单元的力学特性和动力学响应，是后续分析的关键数据。将超单元与残余结构进行组合求解，得到整个结构的动力学响应。在组合过程中，需要考虑超单元与残余结构之间的连接关系和相互作用，确保组合后的结构能够准确反映实际结构的力学行为。通过求解组合后的系统方程，可以得到结构的固有频率、振型、应力分布等重要信息，这些信息对于评估结构的动态性能和安全性具有重要意义。例如，在船舶总振动特性分析中，通过超单元法可以准确地计算出船舶在不同工况下的振动响应，为船舶的设计、优化和运行提供有力的支持。2.2.2超单元法的优势与特点超单元法在船舶总振动特性分析以及其他复杂结构动力学问题的研究中，展现出多方面的显著优势和独特特点，这些优势和特点使其成为一种备受关注和广泛应用的数值分析方法。在降低计算量和提高计算效率方面，超单元法具有明显的优势。传统的有限元分析方法在处理大型复杂结构时，由于模型规模庞大，自由度数量众多，导致计算量急剧增加，计算时间大幅延长。而超单元法通过将结构分解为多个超单元和残余结构，对每个超单元进行独立分析，然后再进行组合求解，有效地减少了整体模型的自由度数量。以船舶结构为例，一艘大型船舶的有限元模型可能包含数百万个自由度，而采用超单元法后，可将模型划分为多个超单元，每个超单元的自由度数量相对较少，从而大大降低了计算量。相关研究表明，在处理相同规模的船舶结构模型时，超单元法的计算时间相比传统有限元法可缩短30%-50%，计算效率得到了显著提高。这使得在实际工程应用中，能够在更短的时间内得到结构的动力学响应结果，为工程决策提供及时的支持。超单元法具有良好的模块化和可复用性。在实际工程中，许多结构是由多个相似或相同的部件组成，这些部件在不同的结构中可能具有相同或相似的力学特性。超单元法可以将这些部件定义为超单元，通过一次分析得到超单元的特性参数，然后在不同的结构模型中重复使用。例如，在船舶制造中，同一型号的船舶可能具有相似的船体结构和设备配置，通过将船体的某些部件定义为超单元，在设计新的船舶时，只需对超单元进行简单的参数调整，即可快速应用到新的模型中，大大提高了设计效率。这种模块化和可复用性不仅减少了重复计算的工作量，还便于对结构进行管理和维护，提高了工程设计的规范性和标准化程度。超单元法还具有保密性强的特点。在一些涉及商业机密或技术保密的工程项目中，不同的设计团队或部门可能需要对各自负责的部分进行保密。超单元法允许各团队将自己的设计部分定义为超单元，只对外提供超单元的接口信息和分析结果，而不暴露内部的详细设计和计算过程。例如，在船舶动力系统的设计中，动力系统供应商可以将其设计的发动机、传动装置等部分定义为超单元，只向船舶总体设计单位提供超单元的力学性能参数和接口信息，保护了自身的技术机密。这种保密性强的特点使得超单元法在涉及多方合作的工程项目中具有重要的应用价值，促进了不同团队之间的合作与交流。超单元法特别适用于大规模模型的分析。随着现代工程技术的不断发展，结构的规模和复杂性日益增加，传统的分析方法往往难以满足实际需求。超单元法能够有效地处理大规模模型，通过合理划分超单元和残余结构，充分利用计算机的计算资源，实现对复杂结构的精确分析。在船舶总振动特性分析中，对于大型集装箱船、邮轮等超大型船舶，其结构复杂，尺寸巨大，采用超单元法可以更好地模拟船舶的实际工作状态，准确预测船舶的总振动特性，为船舶的设计和优化提供可靠的依据。2.2.3超单元法与传统有限元法的比较在船舶总振动特性分析以及更广泛的结构动力学研究领域，超单元法与传统有限元法都是重要的分析工具，它们在计算效率、模型规模处理能力以及适用场景等方面存在着显著的差异。计算效率是衡量分析方法优劣的重要指标之一。传统有限元法在处理复杂结构时，需要对整个结构进行全面的离散化处理，形成大规模的有限元模型。这导致模型的自由度数量庞大，计算量急剧增加，计算时间往往较长。以一艘大型油轮的总振动特性分析为例，若采用传统有限元法，构建完整的有限元模型可能包含数百万个自由度，在普通计算机配置下，计算一次固有频率和振型可能需要数小时甚至数天的时间。而超单元法通过将结构分解为多个超单元和残余结构，对每个超单元进行独立分析，然后再进行组合求解。这种方式大大减少了整体模型的自由度数量，从而显著提高了计算效率。同样以该大型油轮为例，采用超单元法划分模型后，每个超单元的自由度数量相对较少，计算时间可大幅缩短，可能只需几十分钟到数小时即可完成相同的分析任务，计算效率提高数倍甚至数十倍。在模型规模处理能力方面，两者也有明显区别。传统有限元法对于小型或简单结构的分析具有较高的精度和可靠性，能够较为准确地模拟结构的力学行为。然而，当面对大规模、复杂结构时，由于模型规模的限制，传统有限元法可能会遇到计算资源不足、计算精度下降等问题。例如，对于一艘具有复杂内部结构和众多设备的海洋科考船，传统有限元法在构建完整模型时，可能会因为内存不足而无法进行计算，或者在计算过程中由于舍入误差等原因导致计算结果的精度受到影响。超单元法则在处理大规模模型时具有明显优势。它可以将复杂结构合理地划分为多个超单元，每个超单元的规模相对较小，便于进行独立分析和计算。同时，超单元法通过合理的矩阵缩聚技术，进一步减少了计算量，使得能够处理更大规模的模型。对于上述海洋科考船，采用超单元法可以将船体结构、设备等分别划分为不同的超单元，有效地降低了模型的复杂度，提高了模型的处理能力，能够更准确地分析船舶的总振动特性。从适用场景来看，传统有限元法适用于对结构细节要求较高、模型规模较小的情况。在一些小型船舶或船舶部件的设计和分析中，传统有限元法可以精确地模拟结构的每一个细节，为设计提供详细的力学性能数据。例如，在小型游艇的结构设计中，设计师需要精确了解船体各部位的应力分布和变形情况，传统有限元法能够满足这一需求。而超单元法更适用于对整体结构动态性能分析要求较高、模型规模较大的场景。在大型船舶的总振动特性分析中，超单元法能够快速准确地计算出船舶的固有频率、振型等重要参数，评估船舶在不同工况下的振动响应，为船舶的整体设计和优化提供关键依据。此外，在涉及多方合作的项目中，超单元法的保密性和模块化特点使其具有独特的优势，能够满足不同团队之间对模型信息共享和保护的需求。三、基于超单元法修正总振动经验公式的方法3.1超单元法的实施步骤3.1.1结构模型划分在运用超单元法对船舶总振动特性进行分析时，结构模型划分是首要且关键的步骤。其核心在于将复杂的整体结构合理地划分为超单元和残余结构，这一过程需要遵循一定的原则并采用恰当的方法，以确保后续分析的准确性和高效性。在划分原则方面，结构的几何形状是重要考量因素之一。例如，对于船舶结构，可依据船体的舱段划分，将不同的舱室分别定义为超单元。因为不同舱段的几何形状和尺寸各异，各自具有独特的力学特性。如货舱舱段通常具有较大的空间和承载货物的功能，其结构设计需满足货物装载和运输的要求；而机舱舱段则布置有各种机械设备，对结构的强度和振动特性有不同的要求。基于几何形状进行划分，能够使每个超单元具有相对独立和完整的结构特征，便于后续的分析和处理。功能特点也是划分的重要依据。船舶的不同部分承担着不同的功能，上层建筑主要用于船员生活和工作，其结构设计需考虑舒适性和居住要求；而船底结构则主要承受水压力和波浪冲击力，对结构的强度和稳定性要求较高。根据功能特点划分超单元，能够更好地反映各部分结构在船舶总振动中的作用和响应特性。连接方式同样不容忽视。船舶结构中各部分之间通过焊接、螺栓连接等方式相互连接，连接部位的力学性能和传力特性与其他部位不同。在划分超单元时，应尽量将连接点作为超单元的边界，这样可以更准确地模拟结构的力学行为，避免在连接部位出现应力集中等问题对分析结果的影响。以一艘集装箱船为例，在划分结构模型时，首先根据几何形状，将船体沿长度方向划分为多个舱段，每个舱段作为一个超单元。然后，考虑功能特点，将货舱舱段、机舱舱段、上层建筑等分别作为独立的超单元进行处理。对于连接方式，将各舱段之间的连接焊缝以及上层建筑与船体之间的连接部位作为超单元的边界。在划分过程中，运用专业的有限元建模软件，如ANSYS、ABAQUS等，通过设置合适的网格密度和单元类型，对每个超单元进行精确的建模。对于货舱舱段，采用较大尺寸的单元来提高计算效率，同时保证能够准确模拟其整体力学性能；对于上层建筑等对局部细节要求较高的部分，则采用较小尺寸的单元，以更精确地反映其结构特性。通过这样的划分方式，将复杂的集装箱船结构划分为多个超单元和残余结构，为后续的超单元生成和组合分析奠定了坚实的基础。3.1.2超单元生成与缩聚超单元生成与缩聚是超单元法实施过程中的关键环节，通过特定的方法生成超单元并对其自由度进行缩聚，能够有效降低计算的复杂性，提高分析效率。在超单元生成过程中，模态综合法是常用的方法之一。该方法基于结构动力学理论，通过对超单元内部的模态进行分析和综合，来确定超单元的动力学特性。其基本原理是将超单元的位移场表示为一系列模态的线性组合，这些模态包括固定界面模态和约束模态等。固定界面模态是指在超单元边界固定的情况下，超单元内部的固有振动模态，它反映了超单元内部结构的振动特性；约束模态则是为了考虑超单元边界的约束条件而引入的，它描述了超单元边界在受到约束时的变形情况。通过计算这些模态的频率和振型，并将它们进行组合，就可以得到超单元的动力学方程。以船舶的一个舱段超单元为例，首先利用有限元方法计算该舱段在固定边界条件下的固有频率和振型，得到固定界面模态。然后，根据舱段与其他部分的连接方式，确定约束条件，计算约束模态。将固定界面模态和约束模态进行线性组合，得到超单元的位移场表达式，进而建立超单元的动力学方程。通过这种方式生成的超单元，能够准确地反映该舱段在船舶整体结构中的动力学行为。动力缩聚法也是超单元生成的重要方法。动力缩聚法的核心思想是基于结构的动力学平衡方程，通过一定的变换，将超单元的自由度划分为主自由度和从自由度，然后消去从自由度，从而实现自由度的缩聚。在实际应用中，通常选择超单元边界上的节点自由度作为主自由度，因为这些节点与其他超单元或残余结构相连，对整体结构的动力学响应具有重要影响。通过对主自由度的选择和从自由度的消去，可以将超单元的动力学方程从高维空间转换到低维空间，大大降低了计算量。以船舶的上层建筑超单元为例，将上层建筑与船体连接部位的节点自由度作为主自由度，利用动力缩聚法对超单元的自由度进行缩聚。首先，建立上层建筑超单元的动力学平衡方程，然后根据主自由度的选择，将方程中的从自由度用主自由度表示出来，并代入方程中，消去从自由度。经过这样的处理，得到的超单元动力学方程只包含主自由度，其规模大大减小，计算效率得到显著提高。在生成超单元后，还需要对其进行缩聚，以进一步减少计算量。缩聚过程主要是对超单元的刚度矩阵、质量矩阵和阻尼矩阵进行处理。通过缩聚，可以将超单元的矩阵规模减小，同时保留其主要的动力学特性。在缩聚过程中，需要注意保证缩聚后的矩阵能够准确反映超单元的力学行为，避免因缩聚而导致信息丢失或误差增大。通常采用的缩聚方法包括Guyan缩聚法、Krylov子空间法等，这些方法在不同的情况下具有各自的优势和适用范围，需要根据具体问题进行选择和应用。3.1.3超单元与残余结构的组合分析超单元与残余结构的组合分析是超单元法实施的最后一个关键步骤，通过将生成并缩聚后的超单元与残余结构进行合理组合，能够准确地求解整个结构的动力学响应，为船舶总振动特性的分析提供可靠依据。在组合分析之前，需要明确超单元与残余结构之间的连接关系。这种连接关系通常通过节点的位移协调条件和力的平衡条件来描述。在船舶结构中，超单元与残余结构之间的连接方式多种多样，如焊接、螺栓连接等。对于焊接连接，超单元与残余结构在连接部位的节点位移是连续的，即它们在连接点处的位移相等；同时，连接点处的内力也是平衡的，满足力的平衡方程。对于螺栓连接，虽然连接部位存在一定的间隙和接触非线性，但在简化分析中，可以将其视为弹性连接，通过设置合适的弹簧刚度来模拟螺栓的连接作用，同样满足节点位移协调和力的平衡条件。在确定连接关系后，将超单元与残余结构的运动方程进行组合。超单元的运动方程通常表示为：M_s\ddot{u}_s+C_s\dot{u}_s+K_su_s=F_s，其中M_s、C_s、K_s分别为超单元的质量矩阵、阻尼矩阵和刚度矩阵，u_s为超单元的位移向量，F_s为作用在超单元上的外力向量。残余结构的运动方程为：M_r\ddot{u}_r+C_r\dot{u}_r+K_ru_r=F_r，其中各参数含义与超单元类似。通过连接关系，将超单元与残余结构的位移向量和力向量进行关联，得到整个结构的运动方程：\begin{bmatrix}M_s&0\\0&M_r\end{bmatrix}\begin{bmatrix}\ddot{u}_s\\\ddot{u}_r\end{bmatrix}+\begin{bmatrix}C_s&0\\0&C_r\end{bmatrix}\begin{bmatrix}\dot{u}_s\\\dot{u}_r\end{bmatrix}+\begin{bmatrix}K_s&K_{sr}\\K_{rs}&K_r\end{bmatrix}\begin{bmatrix}u_s\\u_r\end{bmatrix}=\begin{bmatrix}F_s\\F_r\end{bmatrix}，其中K_{sr}和K_{rs}为超单元与残余结构之间的耦合刚度矩阵，反映了它们之间的相互作用。求解组合后的运动方程，可以得到整个结构的动力学响应，包括位移、速度、加速度等。在求解过程中，可采用多种数值方法，如直接积分法（如Newmark法、Wilson-\theta法）、模态叠加法等。直接积分法通过对运动方程进行逐步积分，直接求解结构在各个时刻的响应；模态叠加法则是将结构的响应表示为一系列模态的叠加，通过求解模态方程来得到结构的总响应。以一艘油轮的总振动分析为例，采用超单元法将船体划分为多个超单元和残余结构，组合分析后，利用Newmark法求解运动方程，得到油轮在波浪力作用下的振动响应。通过分析这些响应，可以评估油轮的振动特性，如固有频率、振型等，为油轮的结构设计和优化提供重要依据。在分析结果的基础上，如果发现油轮的某些部位振动过大，可以通过调整超单元的结构参数或改变残余结构的布置，重新进行组合分析，直到满足设计要求为止。3.2数据采集与处理3.2.1基于超单元法的数据采集策略在利用超单元法对船舶总振动经验公式进行修正的研究中，数据采集是关键环节之一，其策略的合理性和有效性直接影响到后续分析的准确性和可靠性。为了全面获取船舶在不同工况下的振动特性数据，需要基于超单元法制定科学的数据采集策略。首先，确定不同的主尺度参数。主尺度参数是描述船舶外形和尺寸的重要指标，包括船长、型宽、型深、吃水等。这些参数的变化会显著影响船舶的结构动力学特性和总振动响应。例如，船长的增加会改变船舶的整体刚度和质量分布，从而对振动频率和振型产生影响；型宽的变化会影响船舶的稳性和水动力特性，进而影响总振动特性。通过有针对性地选取不同的主尺度参数组合，可以更全面地涵盖船舶的各种可能设计方案，为后续的分析提供丰富的数据样本。在确定主尺度参数后，利用超单元模型进行数据采集。以一艘集装箱船为例，首先根据其结构特点和功能分布，将船体划分为多个超单元，如货舱超单元、机舱超单元、上层建筑超单元等。针对每个超单元，运用有限元分析软件，如ANSYS、ABAQUS等，进行详细的建模和分析。在建模过程中，准确设置材料属性、单元类型、边界条件等参数，确保模型能够真实反映超单元的力学特性。对于货舱超单元，考虑货物的装载情况对结构的影响，设置不同的货物分布模式；对于机舱超单元，考虑机械设备的振动激励，施加相应的载荷。通过对不同主尺度下的超单元模型进行计算，获取各个超单元在不同工况下的振动响应数据，包括位移、速度、加速度、应力等。这些数据能够反映超单元在船舶整体结构中的动力学行为，为后续的组合分析提供基础。在采集数据时，还需要考虑不同的工况条件。船舶在实际航行过程中会遇到各种不同的工况，如满载、空载、不同航速、不同海况等。这些工况条件会对船舶的总振动特性产生显著影响。例如，在满载工况下，船舶的质量增加，重心位置发生变化，从而影响船舶的振动频率和响应幅值；在不同海况下，波浪力的大小和方向不同，会导致船舶受到不同形式和强度的激励，进而影响总振动特性。因此，在数据采集过程中，需要模拟多种工况条件，对每个主尺度下的超单元模型进行不同工况的计算，获取全面的振动响应数据。通过这样的方式，可以更真实地反映船舶在实际航行中的总振动特性，为修正总振动经验公式提供更准确的数据支持。3.2.2数据处理与分析方法在基于超单元法完成船舶振动数据采集后，数据处理与分析成为挖掘数据潜在价值、揭示船舶总振动规律的关键步骤。通过运用一系列科学合理的数据处理与分析方法，可以有效提高数据质量，准确把握船型主尺度与总振动特性之间的内在联系。数据预处理是数据处理的首要环节，主要包括滤波和去噪等操作。在船舶振动数据采集过程中，由于受到各种因素的干扰，如传感器噪声、环境干扰等，采集到的数据往往包含噪声和干扰信号。这些噪声和干扰信号会影响数据的准确性和可靠性，因此需要进行滤波和去噪处理。采用低通滤波器可以有效去除高频噪声，保留低频的有效信号；采用中值滤波等方法可以去除脉冲噪声，提高数据的平滑度。以某船舶振动加速度数据为例，在采集过程中受到了高频电磁干扰，导致数据出现剧烈波动。通过使用低通滤波器，设置截止频率为50Hz，对数据进行滤波处理后，高频噪声得到了有效抑制，数据的波动明显减小，更能准确反映船舶的实际振动情况。统计分析是深入了解数据特征的重要手段。通过计算数据的均值、方差、标准差等统计量，可以对数据的集中趋势、离散程度等特征有一个全面的认识。例如，计算不同主尺度下船舶振动响应数据的均值，可以了解在该主尺度下船舶振动的平均水平；计算方差和标准差，可以评估数据的离散程度，判断数据的稳定性。对不同工况下的振动响应数据进行相关性分析，能够找出影响船舶总振动特性的关键因素。通过对某船舶在不同航速和不同海况下的振动响应数据进行相关性分析发现，航速与振动响应幅值之间存在显著的正相关关系，海况的恶劣程度也与振动响应幅值密切相关。这些分析结果为进一步研究船舶总振动特性提供了重要的参考依据。回归分析是建立船型主尺度与总振动特性关系模型的核心方法。通过运用回归分析方法，可以找出船型主尺度与总振动特性之间的数学关系，从而建立起总振动预报公式。采用多元线性回归方法，以船长、型宽、型深、吃水等主尺度参数为自变量，以船舶的固有频率、振动响应幅值等总振动特性参数为因变量，进行回归分析。在回归过程中，通过最小二乘法等方法确定回归系数，使得回归方程能够最佳地拟合数据。对大量不同主尺度的船舶振动数据进行回归分析后，得到了船舶固有频率与主尺度参数之间的回归方程：f=a_1L+a_2B+a_3D+a_4T+b，其中f为固有频率，L为船长，B为型宽，D为型深，T为吃水，a_1、a_2、a_3、a_4为回归系数，b为常数项。通过对回归方程的分析和验证，可以评估其准确性和可靠性，为船舶总振动特性的预测提供有效的工具。3.3经验公式的修正与推导3.3.1基于超单元法的经验公式修正思路基于超单元法对总振动经验公式进行修正的核心思路在于，借助超单元法获取的丰富数据，深入剖析船型主尺度与振动特性之间的内在联系，从而对传统经验公式进行针对性的改进。超单元法通过将复杂的船舶结构划分为多个超单元和残余结构，能够利用三维模型进行精确计算，使结构分析更贴合实际情况，为经验公式的修正提供了可靠的数据基础。在获取超单元法数据后，首先对数据进行全面分析。通过统计不同主尺度下船舶振动特性的变化趋势，如固有频率、振动响应幅值等，找出船型主尺度与振动特性之间的初步关系。以船长为例，随着船长的增加，船舶的整体刚度和质量分布会发生变化，进而影响船舶的固有频率和振动响应。通过对超单元法数据的分析，可以发现船长与固有频率之间可能存在某种函数关系，如随着船长的增大，固有频率可能会呈现下降的趋势。基于对数据的分析结果，对传统经验公式进行修正。传统经验公式通常是基于特定类型船舶和有限样本数据建立的，对于新型船型或特殊船型的适用性有限。通过超单元法获取的数据，可以补充和完善传统经验公式中未考虑的因素，或者调整公式中的系数，使其更准确地反映新型船型的总振动特性。例如，对于某新型三体船，传统经验公式在预测其总振动特性时存在较大误差。通过超单元法分析不同主尺度下三体船的振动特性数据，发现三体船的片体间距、连接桥结构等因素对总振动特性有显著影响，而这些因素在传统经验公式中并未得到充分考虑。因此，在修正经验公式时，将这些因素纳入公式中，通过回归分析确定新的系数，从而得到更适合三体船的总振动经验公式。在修正过程中，还需要考虑不同工况对船舶总振动特性的影响。船舶在实际航行中会遇到各种工况，如满载、空载、不同航速、不同海况等，这些工况会导致船舶所受的载荷和边界条件发生变化，进而影响总振动特性。利用超单元法获取不同工况下船舶的振动数据，分析工况因素与振动特性之间的关系，并将其融入经验公式的修正中。通过对不同航速下船舶振动数据的分析，发现航速与振动响应幅值之间存在一定的函数关系，在修正经验公式时，可以引入航速作为一个变量，建立航速与振动响应幅值之间的数学模型，使经验公式能够更准确地预测不同工况下船舶的总振动特性。3.3.2具体修正过程与数学推导以CCS法经验公式为例，详细阐述利用超单元法数据进行修正的数学推导过程。CCS法在我国船舶工程领域应用广泛，但其在面对新型船型时存在一定的局限性，通过超单元法对其进行修正具有重要的工程意义。CCS法垂向振动经验公式一般形式为：f_{v}=a\times(\frac{L}{B})^{b}\times(\frac{d}{T})^{c}\times\sqrt{\frac{\Delta}{L^3}}，其中f_{v}表示垂向振动固有频率，a、b、c为经验系数，L为船长，B为型宽，d为型深，T为吃水，\Delta为排水量。利用超单元法对某新型船型进行分析，获取不同主尺度下该船型的垂向振动固有频率数据。假设通过超单元法计算得到n组数据，每组数据包含船长L_i、型宽B_i、型深d_i、吃水T_i、排水量\Delta_i以及对应的垂向振动固有频率f_{vi}（i=1,2,\cdots,n）。对这些数据进行回归分析，采用多元线性回归方法，将垂向振动固有频率f_{v}表示为船型主尺度参数的函数：f_{v}=\beta_0+\beta_1\ln(\frac{L}{B})+\beta_2\ln(\frac{d}{T})+\beta_3\ln(\sqrt{\frac{\Delta}{L^3}})+\epsilon，其中\beta_0、\beta_1、\beta_2、\beta_3为待确定的回归系数，\epsilon为误差项。为了求解回归系数，采用最小二乘法。最小二乘法的目标是使观测值f_{vi}与回归模型预测值\hat{f}_{vi}之间的误差平方和最小，即：Q=\sum_{i=1}^{n}(f_{vi}-\hat{f}_{vi})^2=\sum_{i=1}^{n}(f_{vi}-(\beta_0+\beta_1\ln(\frac{L_i}{B_i})+\beta_2\ln(\frac{d_i}{T_i})+\beta_3\ln(\sqrt{\frac{\Delta_i}{L_i^3}})))^2。对Q分别关于\beta_0、\beta_1、\beta_2、\beta_3求偏导数，并令偏导数等于0，得到以下方程组：\begin{cases}\frac{\partialQ}{\partial\beta_0}=-2\sum_{i=1}^{n}(f_{vi}-(\beta_0+\beta_1\ln(\frac{L_i}{B_i})+\beta_2\ln(\frac{d_i}{T_i})+\beta_3\ln(\sqrt{\frac{\Delta_i}{L_i^3}})))=0\\\frac{\partialQ}{\partial\beta_1}=-2\sum_{i=1}^{n}(f_{vi}-(\beta_0+\beta_1\ln(\frac{L_i}{B_i})+\beta_2\ln(\frac{d_i}{T_i})+\beta_3\ln(\sqrt{\frac{\Delta_i}{L_i^3}})))\ln(\frac{L_i}{B_i})=0\\\frac{\partialQ}{\partial\beta_2}=-2\sum_{i=1}^{n}(f_{vi}-(\beta_0+\beta_1\ln(\frac{L_i}{B_i})+\beta_2\ln(\frac{d_i}{T_i})+\beta_3\ln(\sqrt{\frac{\Delta_i}{L_i^3}})))\ln(\frac{d_i}{T_i})=0\\\frac{\partialQ}{\partial\beta_3}=-2\sum_{i=1}^{n}(f_{vi}-(\beta_0+\beta_1\ln(\frac{L_i}{B_i})+\beta_2\ln(\frac{d_i}{T_i})+\beta_3\ln(\sqrt{\frac{\Delta_i}{L_i^3}})))\ln(\sqrt{\frac{\Delta_i}{L_i^3}})=0\end{cases}解这个方程组，即可得到回归系数\beta_0、\beta_1、\beta_2、\beta_3的值。将得到的回归系数代入回归模型中，得到修正后的CCS法垂向振动经验公式：f_{v}=\beta_0+\beta_1\ln(\frac{L}{B})+\beta_2\ln(\frac{d}{T})+\beta_3\ln(\sqrt{\frac{\Delta}{L^3}})。通过这样的修正过程，利用超单元法获取的数据对CCS法经验公式进行了优化，使其能够更准确地预测新型船型的垂向振动固有频率。在实际应用中，还需要对修正后的经验公式进行验证和评估，通过与实际测量数据或其他更精确的分析方法结果进行对比，检验公式的准确性和可靠性。如果发现修正后的公式仍存在一定误差，可以进一步调整回归模型或补充更多的数据进行分析，以不断提高经验公式的精度和适用性。四、案例分析4.1案例选取与背景介绍4.1.1特殊船型案例选取在本研究中，为了充分验证基于超单元法修正总振动经验公式的有效性和准确性，选取了小水线面双体船和新型集装箱船作为典型案例。小水线面双体船作为一种极具特色的船型，其结构形式与传统单体船有显著差异。它主要由潜没于水中的鱼雷状下体、高于水平面的箱形平台以及穿越水面连接下体与平台的流线型支柱三部分组成。这种独特的结构设计使得小水线面双体船在水动力性能和结构动力学特性方面表现出与传统船型不同的特点。例如，其水线面面积较小，兴波阻力小，在风浪中受到的扰动力也大大减少，运动加速度小，被誉为“全海候”船舶。然而，由于其结构和水动力特性的复杂性，传统的总振动经验公式在预测小水线面双体船的总振动特性时往往存在较大误差。因此，选取小水线面双体船作为案例，能够深入研究超单元法在修正经验公式方面的优势和潜力，为该船型的设计和优化提供更准确的理论支持。新型集装箱船也是本研究的重要案例之一。随着全球贸易的不断增长，集装箱运输在国际贸易中占据着越来越重要的地位。为了提高运输效率，新型集装箱船不断向大型化、高速化方向发展。这些新型集装箱船在主尺度、结构布局和装载工况等方面与传统集装箱船有很大的不同。例如，大型集装箱船的船长和型宽不断增加，使得船体的刚度和质量分布发生变化，对总振动特性产生影响；同时，高速航行时船舶所受到的水动力作用更加复杂，也增加了总振动特性分析的难度。传统的总振动经验公式难以准确预测新型集装箱船在这些复杂工况下的总振动特性。通过对新型集装箱船进行研究，利用超单元法修正经验公式，可以更好地满足新型集装箱船设计和运营的需求，提高其安全性和经济性。4.1.2船型特点与振动特性需求小水线面双体船的结构特点决定了其对振动特性有着特殊的需求。从结构方面来看，小水线面双体船的下体为鱼雷状，提供大部分浮力，而上体为箱形平台，通过流线型支柱与下体相连。这种结构使得船舶的重心较低，稳性较好，但也导致船体的柔性增加，在受到外界激励时更容易产生振动。在航行工况方面，小水线面双体船通常在近海或沿海区域航行，会遇到各种不同的海况，如风浪、涌浪等。这些复杂的海况会对船舶产生不同形式和强度的激励，导致船舶产生垂向、横向和扭转等多种振动形式。由于小水线面双体船的水线面面积小，在遇到横浪时，容易受到较大的横向力作用，从而引发横向振动；而在高速航行时，由于支柱与水线面的相互作用，可能会产生较大的垂向振动和扭转振动。因此，小水线面双体船对振动特性的要求主要体现在以下几个方面：一是要求船舶具有良好的耐波性，能够在不同海况下保持较低的振动响应，确保船员和乘客的舒适性以及船舶设备的正常运行；二是要求船舶的结构具有足够的强度和刚度，能够承受振动产生的应力和变形，避免结构疲劳损坏；三是要求在设计阶段能够准确预测船舶的总振动特性，为结构设计和优化提供科学依据。新型集装箱船的船型特点和振动特性需求也有其独特之处。新型集装箱船的大型化使得船体的尺度大幅增加，船长、型宽和型深都比传统集装箱船大很多。这种大尺度的结构特点导致船体的刚度相对降低，在受到波浪力、螺旋桨激振力等外界载荷作用时，更容易产生较大的振动。新型集装箱船通常在远洋航线上航行，航行距离长，遇到的海况更加复杂多变。在满载和空载工况下，船舶的重心位置和质量分布会发生较大变化，这对船舶的总振动特性也会产生显著影响。在满载时，船舶的吃水增加，船体受到的浮力和阻力增大，振动响应可能会相应增加；而在空载时，船舶的重心升高，稳性相对降低，振动的敏感性可能会增加。因此，新型集装箱船对振动特性的需求主要包括：一是要保证在不同装载工况和海况下，船舶的振动响应在允许范围内，以确保货物的安全运输和船舶结构的完整性；二是要优化船舶的结构设计，提高船体的刚度和强度，降低振动对船舶的影响；三是需要准确预测船舶在各种工况下的总振动特性，为船舶的运营管理提供参考，合理安排航线和航速，减少振动对船舶的损害。4.2基于超单元法的分析过程4.2.1建立超单元模型以小水线面双体船为例，运用专业的有限元分析软件ANSYS建立其超单元模型。在划分结构模型时，充分考虑小水线面双体船的结构特点和功能分布。将船体的下体、支柱和上体分别划分为不同的超单元。下体作为提供浮力的主要部分，其结构较为复杂，且在水中受到的压力和浮力分布不均匀，因此将下体按照不同的舱段进一步细分，每个舱段作为一个独立的超单元。支柱连接上下体，承受着较大的剪切力和弯矩，将其单独划分为超单元，便于对其力学性能进行精确分析。上体主要用于布置各种设备和设施，根据其功能区域的不同，划分为不同的超单元，如驾驶舱超单元、生活舱超单元等。在划分超单元时，合理设置网格密度和单元类型。对于下体，由于其承受较大的水压力，为了保证计算精度，采用较小尺寸的单元，如SOLID185单元，该单元具有良好的计算精度和适应性，能够准确模拟下体在复杂水动力作用下的力学响应。对于支柱，考虑到其细长的结构特点和主要承受的载荷类型，采用BEAM188单元，该单元适用于模拟梁类结构，能够准确计算支柱的弯曲和剪切变形。对于上体，根据其结构的复杂程度和对计算精度的要求，采用SHELL181单元，该单元适用于模拟薄板和薄壳结构，能够有效地模拟上体的力学行为。确定超单元的边界条件。下体与水接触的表面施加水压力边界条件，根据流体力学理论，水压力随深度呈线性变化，通过计算不同深度处的水压力，将其施加到下体的相应节点上。支柱与下体和上体的连接部位采用刚性连接边界条件，确保力和位移的有效传递。上体的外表面根据实际情况施加风载荷和其他外部载荷边界条件，考虑到小水线面双体船在不同航行工况下可能受到不同方向和大小的风载荷，通过风洞试验或数值模拟等方法获取风载荷数据，并将其施加到上体的相应节点上。通过以上步骤，成功建立了小水线面双体船的超单元模型，为后续的超单元生成和组合分析奠定了基础。4.2.2数据采集与计算结果利用建立好的小水线面双体船超单元模型，进行不同主尺度下的振动数据采集。改变船长、型宽、型深等主尺度参数，分别计算不同参数组合下小水线面双体船的振动特性。设置船长从50米到100米，以10米为间隔；型宽从10米到20米，以2米为间隔；型深从5米到10米，以1米为间隔，共得到多组不同主尺度的模型。针对每组模型，在ANSYS软件中进行计算分析，获取其固有频率、振型等振动特性数据。计算结果显示，随着船长的增加，小水线面双体船的固有频率呈现下降趋势。当船长从50米增加到100米时，第一阶固有频率从15Hz左右下降到10Hz左右。这是因为船长增加导致船体的整体刚度相对降低，质量分布发生变化，使得船体在振动时更容易产生变形，从而降低了固有频率。型宽的增加对固有频率的影响较为复杂，在一定范围内，型宽增加会使船体的横向刚度增大，固有频率有所提高；但当型宽超过一定值后，由于船体的质量增加和水动力特性的变化，固有频率反而会下降。当型宽从10米增加到15米时，固有频率从12Hz左右提高到13Hz左右；而当型宽继续增加到20米时，固有频率又下降到12.5Hz左右。型深的增加会使船体的垂向刚度增大，固有频率相应提高。当型深从5米增加到10米时，第一阶固有频率从10Hz左右提高到12Hz左右。通过计算还得到了小水线面双体船在不同主尺度下的振型。第一阶振型主要表现为船体的整体弯曲振动，随着主尺度的变化，弯曲变形的幅度和位置也有所不同。在小尺度模型中，船体的弯曲变形主要集中在中部；而在大尺度模型中，由于船体的刚度分布变化，弯曲变形在船体的两端也较为明显。第二阶振型则表现为船体的扭转振动，扭转的角度和方向与主尺度参数密切相关。在型宽较大的模型中，扭转振动的幅度相对较小，因为较大的型宽增加了船体的抗扭刚度；而在型深较大的模型中，扭转振动的方向可能会发生改变，这是由于型深的变化影响了船体的重心位置和惯性矩。这些计算结果为深入研究小水线面双体船的总振动特性提供了重要的数据支持，也为后续的经验公式修正奠定了基础。4.3经验公式修正结果与验证4.3.1修正后的经验公式呈现通过超单元法对小水线面双体船进行分析，并运用数据回归分析方法，得到了修正后的总振动经验公式。以垂向振动固有频率f_{v}为例，修正后的经验公式为：f_{v}=0.05+0.8\ln(\frac{L}{B})+1.2\ln(\frac{d}{T})+0.5\ln(\sqrt{\frac{\Delta}{L^3}})+0.3\ln(\frac{S}{L^2})其中，L为船长，B为型宽，d为型深，T为吃水，\Delta为排水量，S为支柱横截面积。与传统的CCS法垂向振动经验公式相比，修正后的公式增加了支柱横截面积S这一参数。这是因为小水线面双体船的支柱在连接上下体过程中，其横截面积对船舶的垂向振动特性有着显著影响。通过超单元法对不同支柱横截面积下的小水线面双体船进行分析，发现支柱横截面积与垂向振动固有频率之间存在一定的函数关系，因此将其纳入经验公式中，以提高公式对小水线面双体船垂向振动特性预测的准确性。对于水平振动固有频率f_{h}，修正后的经验公式为：f_{h}=-0.03+0.6\ln(\frac{L}{B})+0.9\ln(\frac{d}{T})+0.4\ln(\sqrt{\frac{\Delta}{L^3}})-0.2\ln(\frac{l}{L})其中，l为连接桥长度。在小水线面双体船中，连接桥长度对水平振动特性有重要影响。随着连接桥长度的变化，船舶的横向刚度和质量分布会发生改变，从而影响水平振动固有频率。通过超单元法的计算分析，确定了连接桥长度与水平振动固有频率之间的关系，并将其融入经验公式中。这使得修正后的经验公式能够更准确地反映小水线面双体船在不同主尺度和结构参数下的水平振动特性，为该船型的设计和分析提供了更可靠的依据。4.3.2与传统经验公式对比验证为了验证修正后经验公式的准确性和可靠性，将其计算结果与传统经验公式（如CCS法）的计算结果以及实际测量数据进行对比分析。以某小水线面双体船为例，该船的船长L=80米，型宽B=16米，型深d=8米，吃水T=6米，排水量\Delta=1500吨，支柱横截面积S=2平方米，连接桥长度l=10米。首先，运用传统的CCS法经验公式计算该船的垂向振动固有频率f_{v1}和水平振动固有频率f_{h1}。根据CCS法垂向振动经验公式f_{v}=a\times(\frac{L}{B})^{b}\times(\frac{d}{T})^{c}\times\sqrt{\frac{\Delta}{L^3}}，代入相应参数，取经验系数a=0.2，b=0.3，c=0.4，计算得到f_{v1}=0.2\times(\frac{80}{16})^{0.3}\times(\frac{8}{6})^{0.4}\times\sqrt{\frac{1500}{80^3}}\approx1.25Hz。对于水平振动固有频率，根据CCS法水平振动经验公式（假设为f_{h}=m\times(\frac{L}{B})^{n}\times(\frac{d}{T})^{p}\times\sqrt{\frac{\Delta}{L^3}}，取经验系数m=0.1，n=0.2，p=0.3），计算得到f_{h1}=0.1\times(\frac{80}{16})^{0.2}\times(\frac{8}{6})^{0.3}\times\sqrt{\frac{1500}{80^3}}\approx0.85Hz。然后，使用修正后的经验公式计算该船的垂向振动固有频率f_{v2}和水平振动固有频率f_{h2}。根据修正后的垂向振动经验公式f_{v}=0.05+0.8\ln(\frac{L}{B})+1.2\ln(\frac{d}{T})+0.5\ln(\sqrt{\frac{\Delta}{L^3}})+0.3\ln(\frac{S}{L^2})，代入参数计算得到f_{v2}=0.05+0.8\ln(\frac{80}{16})+1.2\ln(\frac{8}{6})+0.5\ln(\sqrt{\frac{1500}{80^3}})+0.3\ln(\frac{2}{80^2})\approx1.45Hz。根据修正后的水平振动经验公式f_{h}=-0.03+0.6\ln(\frac{L}{B})+0.9\ln(\frac{d}{T})+0.4\ln(\sqrt{\frac{\Delta}{L^3}})-0.2\ln(\frac{l}{L})，代入参数计算得到f_{h2}=-0.03+0.6\ln(\frac{80}{16})+0.9\ln(\frac{8}{6})+0.4\ln(\sqrt{\frac{1500}{80^3}})-0.2\ln(\frac{10}{80})\approx1.05Hz。通过实际测量，该小水线面双体船的垂向振动固有频率f_{v实际}=1.42Hz，水平振动固有频率f_{h实际}=1.02Hz。对比结果显示，传统CCS法计算的垂向振动固有频率与实际测量值的相对误差为\vert\frac{1.25-1.42}{1.42}\vert\times100\%\approx12\%，水平振动固有频率与实际测量值的相对误差为\vert\frac{0.85-1.02}{1.02}\vert\times100\%\approx17\%；而修正后的经验公式计算的垂向振动固有频率与实际测量值的相对误差为\vert\frac{1.45-1.42}{1.42}\vert\times100\%\approx2\%，水平振动固有频率与实际测量值的相对误差为\vert\frac{1.05-1.02}{1.02}\vert\times100\%\approx3\%。由此可见，修正后的经验公式计算结果与实际测量数据更为接近，能够更准确地预测小水线面双体船的总振动特性，验证了其准确性和可靠性。4.3.3实际应用效果评估在小水线面双体船的设计过程中，修正后的经验公式发挥了重要作用。以某海洋科考用小水线面双体船设计项目为例，在初步设计阶段，设计团队运用修正后的经验公式对不同设计方案的总振动特性进行了预测。通过对船长、型宽、型深、支柱横截面积、连接桥长度等参数的调整，利用修正后的经验公式计算出不同方案下船舶的垂向和水平振动固有频率。根据计算结果，设计团队对方案进行了优化，选择了振动特性更优的方案，避免了因振动问题导致的结构疲劳和设备损坏风险。在后续的详细设计和施工过程中，修正后的经验公式也为结构设计和设备选型提供了重要依据。通过准确预测船舶的总振动特性，设计团队能够合理设计船体结构的强度和刚度，选择合适的振动隔离和减振设备，确保船舶在航行过程中的稳定性和舒适性。在性能优化方面，修正后的经验公式同样具有显著优势。对于一艘已投入使用的小水线面双体船，运营公司希望通过优化船舶的性能来提高其航行效率和安全性。利用修正后的经验公式，分析了船舶在不同装载工况和航行速度下的总振动特性。发现当船舶在满载工况下以较高速度航行时，振动响应较大，可能影响船舶的结构安全和设备正常运行。根据分析结果，运营公司采取了调整装载方案和优化航行速度等措施，有效降低了船舶的振动响应，提高了船舶的性能和安全性。然而，修正后的经验公式也存在一定的局限性。它是基于特定类型的小水线面双体船数据建立的，对于一些具有特殊结构或功能的小水线面双体船，其适用性可能会受到影响。当小水线面双体船采用新型材料或特殊的连接方式时，修正后的经验公式可能无法准确预测其总振动特性。修正后的经验公式虽然考虑了一些主要的结构参数和工况因素，但对于一些复杂的非线性因素，如材料的非线性特性、结构的几何非线性等，还无法完全考虑。在未来的研究中，可以进一步拓展修正后经验公式的适用范围，考虑更多的非线性因素，提高其预测精度和可靠性。五、结论与展望5.1研究成果总结本研究成功引入超单元法对船舶总振动经验公式进行修正，取得了一系列具有重要理论和实践价值的成果。通过深入剖析船舶总振动理论以及超单元法的基本原理，明确了超单元法在船舶结构分析中的独特优势。超单元法基于动态子结构法，能够将复杂的船舶结构模型分解为多个超单元和残余结构进行独立分析，大大降低了计算的复杂性，提高了计算效率。在船舶总振动特性分析中，超单元法利用三维模型进行计算，使结构分析更加接近实际情况，有效提高了预测的准确性。在对超单元法实施步骤的研究中，详细阐述了结构模型划分、超单元生成与缩聚以及超单元与残余结构组合分析的全过程。在结构模型划分阶段，依据船舶结构的几何形状、功能特点和连接方式等因素，将船舶合理划分为多个超单元和残余结构，为后续分析奠定了基础。在超单元生成与缩聚过程中，采用模态综合法和动力缩聚法等方法，生成超单元并对其自由度进行缩聚，进一步提高了计算效率。在超单元与残余结构的组合分析中，通过合理确定连接关系，将超单元与残余结构的运动方程进行组合求解，得到了整个结构的动力学响应，为船舶总