基于Matlab的双拖轮协助船舶操纵运动仿真：模型构建与应用分析

基于Matlab的双拖轮协助船舶操纵运动仿真：模型构建与应用分析一、引言1.1研究背景与意义1.1.1研究背景随着世界经济的快速发展，全球贸易量持续增长，航运业作为国际贸易的主要运输方式，变得越来越繁荣，船舶大型化的趋势也愈发显著。大型船舶凭借其载重量大、运输成本低等优势，在航运市场中占据着举足轻重的地位。然而，船舶大型化也给航运业带来了一系列新的挑战，尤其是在港池、航道等受限水域内操纵大型船舶的难度大幅增加。在受限水域中，操船水域面积有限，且周围往往存在浅滩等威胁航行安全的水域，这使得船舶必须严格限制在可航水域航行，对船舶操纵的精准度提出了极高的要求。同时，受限水域的相对水深较浅，船舶在这种环境下的操纵性能相较于深水中会变差。大型船舶自身尺度大、惯性大，这使得其操纵性能天生就劣于小船，在受限水域中这种劣势更加明显。此外，相对于小型船舶而言，大型船舶对拖轮的作用效果反应明显变差，增加了拖轮协助操纵的难度。为了降低大型船舶在受限水域航行的难度，船舶驾驶员通常需要借助拖轮协助船舶完成进出港航道、旋回水域以及靠离泊操纵。在大型船舶靠离旋回和靠离泊的过程中，单拖轮一般难以达到理想的作用效果，也难以保证船舶作业的安全，因此，两艘甚至多艘拖轮共同协助船舶操纵成为保证船舶安全的重要手段。尽管多拖轮作用在同一船舶上时，能使船舶出现相对比较理想的作用效果，然而，拖轮与拖轮的相互协作，拖轮与船舶的相互协作成为人们关注的问题。1.1.2研究意义本研究具有重要的理论与实际意义，主要体现在以下几个方面：为船舶操纵人员提供理论依据：通过定量计算与仿真试验相结合的方式，深入研究拖轮以不同协助方式作用时船舶的一般运动规律。这些研究成果能够为引航员和驾驶人员在实际操作中提供有价值的参考依据，帮助他们更好地理解拖轮协助船舶操纵的原理和方法，从而更加准确、安全地完成船舶操纵任务。为航海模拟器的研制与开发提供理论支持：航海模拟器是操船人员进行教育培训、科学研究以及环境评价的重要工具，其核心部分是船舶运动数学模型。本研究对拖轮协助船舶操纵运动的深入分析和建模，能够为航海模拟器中船舶运动数学模型的完善提供重要的理论支持，提高模拟器对船舶操纵过程的模拟精度，使其能够更真实地反映船舶在实际操纵中的运动状态，为操船人员提供更有效的培训和研究平台。1.2国内外研究现状1.2.1国外研究进展国外对于船舶操纵性的研究起步较早。1929年，美国斯蒂文斯工学院的Taylor和Norlin研究发现船舶操纵性可以从航向稳定性和回转性两个方面来衡量，首次提出了操纵性的概念，并提出了船舶运动的数学方程，为船舶运动数学模型的发展奠定了坚实基础。此后，船舶操纵性的研究不断深入和拓展。1957年，日本学者野本谦作运用控制理论来分析船舶的运动，提出了一阶近似船舶操纵运动方程和野本试验确定K、T指数的方法，为船舶操纵运动的研究提供了新的思路和方法。20世纪70年代末，日本拖曳水池委员会成立了船舶操纵运动数学模型研讨小组（JMAM），以田辺、黒田、吉田等人为代表的学者们提出了一系列MMG模型。该模型以船、桨、舵的单独性能为基础，能够简洁地表示船-桨-舵的干涉效应，并合理地表达作用于船上的各种流体动力。它是目前国际上较为流行的一种船舶运动数学模型，并且一直在不断完善之中。利用已发表的试验数据，即便研究者不具备试验条件，采用这种研究方法也能构建出具有相当精度的数学模型，用以预报船舶的动态响应，在船舶初步设计阶段进行操纵性预报具有较高的可行性。在拖轮协助船舶操纵方面，国外也进行了大量的研究。一些研究通过物理模型试验和数值模拟，深入分析拖轮的拖力、顶推位置以及与船舶的相互作用等因素对船舶操纵性能的影响。例如，研究不同拖轮配置方案下船舶在受限水域中的转向、制动和靠泊性能，为实际操作提供了科学的指导。同时，国外还在不断研发新的拖轮技术和设备，以提高拖轮协助船舶操纵的效率和安全性，如新型的全回转拖轮，具有更好的操纵灵活性和拖力性能。1.2.2国内研究情况我国的操纵性仿真预报研究起步相对较晚。早期，周昭明等人率先开展了操纵性仿真预报相关工作，为后续的研究奠定了基础。此后，国内众多科研机构和高校，如上海交通大学、哈尔滨工程大学、中国船舶科学研究中心等，纷纷投身于船舶操纵性的研究领域，取得了一系列具有重要价值的成果。在船舶操纵运动数学模型方面，国内学者在借鉴国外先进模型的基础上，结合我国船舶的特点和实际航行环境，进行了深入的研究和改进。通过理论分析、数值计算和实验研究等多种手段，对船舶在不同工况下的水动力性能进行了全面的研究，建立了更加符合我国国情的船舶操纵运动数学模型。在拖轮协助船舶操纵的研究方面，国内学者针对大型船舶在受限水域中的操纵问题，开展了大量的研究工作。研究内容涵盖拖轮的合理配置、拖轮与船舶的协同操纵策略以及拖轮协助操纵过程中的安全性评估等多个方面。例如，通过建立拖轮-船舶耦合动力学模型，模拟分析不同拖轮协助方式下船舶的运动响应，提出了优化的拖轮协助操纵方案。此外，随着计算机技术和计算流体力学（CFD）的飞速发展，国内在船舶操纵性数值模拟方面取得了显著进展。利用CFD技术，可以对船舶周围的流场进行精确模拟，深入研究船舶在各种工况下的水动力特性，为船舶操纵性的研究提供了更加准确和高效的手段。同时，国内还开展了一系列的实船试验和模型试验，对理论研究和数值模拟的结果进行验证和完善，进一步提高了我国在船舶操纵性研究领域的水平。1.3研究内容与方法1.3.1研究内容本研究旨在深入探讨双拖轮协助船舶操纵运动，主要研究内容如下：建立船舶操纵运动数学模型：基于经典的船舶操纵理论，考虑船舶在水平面内的三自由度运动，即纵荡、横荡和首摇运动，建立准确描述船舶运动的数学模型。同时，充分考虑船体水动力、附加质量及附加惯性矩、舵力和桨力等因素对船舶运动的影响，确定模型中的各项参数。分析双拖轮协助操纵的作用机制：研究双拖轮在不同位置和拖力作用下对船舶操纵运动的影响，分析拖轮与船舶之间的相互作用力和力矩，明确双拖轮协助操纵的作用机制。通过理论推导和数值计算，得出不同协助方式下船舶的运动响应，为后续的仿真分析提供理论基础。开展基于Matlab的仿真分析：利用Matlab软件平台，搭建双拖轮协助船舶操纵运动的仿真模型。通过设置不同的仿真工况，如船舶初始状态、拖轮配置方式、外界环境条件等，对船舶操纵运动进行模拟仿真。分析仿真结果，总结双拖轮协助船舶操纵的运动规律和特点，评估不同协助方式的效果。研究结果的验证与应用：将仿真结果与实际船舶操纵数据或已有的研究成果进行对比验证，评估模型的准确性和可靠性。根据研究结果，提出双拖轮协助船舶操纵的优化策略和建议，为实际船舶操纵提供科学依据和参考。1.3.2研究方法本研究采用理论分析、数学建模、仿真实验相结合的方法，具体如下：理论分析：深入研究船舶操纵性的基本理论，包括船舶运动方程、水动力理论、拖轮拖力计算等。通过对相关理论的分析和推导，明确双拖轮协助船舶操纵运动的基本原理和影响因素，为后续的研究提供理论支持。数学建模：根据理论分析的结果，建立船舶操纵运动的数学模型以及双拖轮与船舶之间的相互作用模型。运用数学方法对模型进行求解和分析，确定模型中各项参数的取值范围和计算方法，为仿真实验提供数学基础。仿真实验：基于Matlab软件平台，利用其强大的计算和绘图功能，对建立的数学模型进行仿真实现。通过设置不同的仿真参数和工况，模拟双拖轮协助船舶操纵的各种实际情况，获取大量的仿真数据。对仿真数据进行分析和处理，总结双拖轮协助船舶操纵的运动规律和特点，评估不同协助方式的效果。二、双拖轮协助船舶操纵运动的理论基础2.1船舶操纵运动数学模型2.1.1MMG模型原理MMG模型（ManeuveringModelingGroupModel），即操纵性建模小组模型，由日本船舶操纵运动数学模型研讨小组提出，是一种被广泛应用于船舶操纵性分析和仿真预测的理论模型，能够描述和模拟船舶在水中的运动，特别是转向和操纵等动态行为。该模型以船、桨、舵的单独性能为基础，通过解析船舶与水的相互作用力，简洁地表示船-桨-舵的干涉效应，合理地表达作用于船上的各种流体动力，进而预测其运动响应。MMG模型基于多自由度的数学模型，将船舶在水平面内的运动分解为纵荡、横荡和首摇三个自由度的运动。在模型中，充分考虑了船舶运动时所受到的各种力和力矩，包括船体水动力、附加质量及附加惯性矩、舵力、桨力以及外界环境干扰力等。通过建立这些力和力矩与船舶运动状态之间的数学关系，来准确描述船舶的运动过程。例如，在船舶回转过程中，模型会综合考虑舵力产生的转船力矩、船体因回转而受到的水动力矩以及螺旋桨推力的变化等因素，从而精确地模拟船舶的回转轨迹和运动参数。该模型将船舶视为一个整体系统，同时又把船舶细分为船体、舵、螺旋桨等多个子系统，分别研究每个子系统的水动力特性以及它们之间的相互干涉关系。通过这种方式，MMG模型能够较为全面地反映船舶操纵运动的本质，为船舶操纵性的研究提供了有力的工具。在实际应用中，MMG模型可以根据船舶的设计参数、航行条件以及操纵指令等输入信息，预测船舶的运动响应，如船舶的航速、航向、位置等随时间的变化情况。这对于船舶设计、航行规划、操纵性能评估以及航海模拟器的开发等方面都具有重要的意义。2.1.2MMG模型的构成要素船体水动力计算：船体在水中运动时，会受到水的作用力，这些力包括摩擦阻力、形状阻力、兴波阻力等。船体水动力的大小和方向与船舶的运动速度、姿态以及船体形状等因素密切相关。在MMG模型中，通常采用经验公式或数值计算方法来计算船体水动力。对于摩擦阻力，可根据平板摩擦阻力公式进行估算，并结合船体的湿表面积进行修正；形状阻力则通过考虑船体的横截面积、流线型等因素来确定；兴波阻力的计算较为复杂，一般采用基于理论或实验数据的半经验公式。同时，船体水动力还会受到附加质量和附加惯性矩的影响，附加质量是指由于船舶运动导致周围水的加速而产生的等效质量，附加惯性矩则是考虑水对船舶转动的惯性影响，它们会改变船舶的运动特性，在模型中需要准确考虑。螺旋桨水动力计算：螺旋桨是船舶的推进装置，其产生的推力是船舶前进的动力来源。螺旋桨的水动力计算主要涉及推力和转矩的计算。螺旋桨的推力与螺旋桨的转速、螺距、直径以及进速等因素有关，通常通过螺旋桨的敞水性征曲线或经验公式来确定。螺旋桨的敞水性征曲线是通过模型试验得到的，它反映了螺旋桨在不同进速系数下的推力系数和转矩系数。在MMG模型中，根据船舶的运动状态确定螺旋桨的进速，进而利用敞水性征曲线计算出螺旋桨的推力和转矩。此外，螺旋桨的伴流效应也会对船舶的操纵性能产生影响，伴流是指螺旋桨旋转时带动周围水一起运动形成的水流，它会使螺旋桨的实际进速与船舶的前进速度不同，在模型中需要对伴流进行合理的模拟和修正。舵水动力计算：舵是船舶操纵的重要装置，通过改变舵角可以产生舵力和舵力矩，从而实现船舶的转向。舵水动力的计算主要考虑舵角、舵面积、舵速以及水流对舵的作用力等因素。在MMG模型中，常用的舵水动力计算方法有基于薄翼理论的方法和经验公式法。基于薄翼理论的方法将舵视为薄翼，通过求解绕翼的流场来计算舵力和舵力矩；经验公式法则是根据大量的实验数据总结得到的，具有一定的工程实用性。舵力和舵力矩不仅与舵本身的参数有关，还会受到船体运动和螺旋桨尾流的影响，例如，螺旋桨尾流会使舵表面的水流速度和压力分布发生变化，从而改变舵力的大小和方向，在模型中需要考虑这些干涉效应。相互干涉关系：在MMG模型中，船体、螺旋桨和舵之间存在着复杂的相互干涉关系。螺旋桨的推力会影响船体的前进速度和周围的流场，进而影响舵力的大小和方向；舵的偏转也会改变船体的运动姿态，反过来影响螺旋桨的工作状态和船体水动力；船体的运动状态则是这些力和力矩综合作用的结果。当船舶转向时，舵角的改变会使船体产生横荡和首摇运动，船体的这些运动又会导致螺旋桨的进速和伴流发生变化，从而影响螺旋桨的推力和转矩；同时，螺旋桨尾流的变化也会对舵力产生影响，使得船舶的转向过程更加复杂。MMG模型通过建立相应的数学方程来描述这些相互干涉关系，从而准确地模拟船舶的操纵运动。2.2拖轮协助操纵的力学原理2.2.1拖轮推力计算方法拖轮的推力是拖轮协助船舶操纵的关键参数，其大小受到多种因素的影响，计算方法也较为复杂。一般来说，拖轮推力可以通过以下几种方式进行计算：基于螺旋桨理论的计算方法：拖轮的推力主要由其螺旋桨产生，基于螺旋桨理论的计算方法是从螺旋桨的工作原理出发，通过考虑螺旋桨的几何参数（如直径、螺距、叶片数等）、运行参数（如转速、进速等）以及流体的物理性质（如水的密度、粘性等）来计算推力。根据螺旋桨的动量定理，螺旋桨推力等于单位时间内通过螺旋桨盘面的流体动量的变化。在实际计算中，通常会引入一些经验系数来修正理论计算结果，以使其更符合实际情况。某拖轮螺旋桨直径为D，转速为n，进速为V，根据Betz理论，其推力T的计算公式可以表示为：T=k_1\rho\pi(\frac{D}{2})^2n^2V+k_2\rho\pi(\frac{D}{2})^4n^4，其中\rho为水的密度，k_1和k_2为经验系数，需要通过实验或经验数据来确定。利用拖轮性能曲线的方法：拖轮在设计和制造过程中，通常会通过实验测量得到其在不同工况下的推力性能曲线。这些性能曲线反映了拖轮推力与螺旋桨转速、拖轮航速等参数之间的关系。在实际应用中，可以根据拖轮的当前运行状态，如螺旋桨转速和航速，从性能曲线中直接查取对应的推力值。这种方法简单直观，但需要准确掌握拖轮的性能曲线数据，并且性能曲线通常是在特定的实验条件下得到的，实际使用时可能需要根据具体情况进行修正。考虑环境因素的影响：拖轮在实际工作中，其推力还会受到环境因素的显著影响，如水流、风力等。水流的存在会改变拖轮的相对速度，从而影响螺旋桨的进速和推力。当拖轮逆水航行时，其实际进速减小，推力会相应增加；反之，顺水航行时，进速增大，推力会减小。风力对拖轮推力的影响主要体现在风阻力上，风阻力会消耗拖轮的部分推力，使有效推力降低。在计算拖轮推力时，需要综合考虑这些环境因素的影响。可以通过建立水流和风力的数学模型，将其与拖轮的推力计算模型相结合，得到在实际环境条件下的拖轮推力。若水流速度为V_c，方向与拖轮航行方向夹角为\theta，风速度为V_w，风阻力系数为C_w，则考虑水流和风力影响后的拖轮推力T'可以表示为：T'=T-\frac{1}{2}\rhoC_wA_wV_w^2\cos(\alpha)-\rhoA_cV_c^2\cos(\theta)，其中A_w为拖轮的迎风面积，A_c为拖轮的迎流面积，\alpha为风向与拖轮轴线的夹角。2.2.2拖轮与船舶的相互作用力分析当拖轮协助船舶操纵时，拖轮与船舶之间通过缆绳或顶推装置等连接方式产生相互作用力，这些力和力矩对船舶的运动状态有着重要的影响。连接点的力和力矩：拖轮与船舶的连接点处会受到拉力或推力以及相应的力矩作用。在顶推方式下，拖轮对船舶施加一个向前的推力，同时由于顶推点与船舶重心不在同一位置，会产生一个使船舶绕重心转动的力矩。若拖轮顶推点在船舶重心前方且偏左，当拖轮施加推力时，除了使船舶向前运动外，还会使船舶产生一个向左的首摇力矩，导致船舶向左转向。在拖带方式下，缆绳对船舶施加拉力，拉力的大小和方向取决于拖轮的推力以及缆绳的布置情况。缆绳的拉力不仅会使船舶产生平移运动，还会因为缆绳与船舶轴线的夹角而产生一个转船力矩。当缆绳与船舶轴线成一定角度时，拉力的横向分力会使船舶产生横荡运动，而纵向分力则影响船舶的前进或后退速度，同时由于力臂的存在，会产生一个使船舶绕重心转动的首摇力矩。对船舶运动的影响：拖轮与船舶之间的相互作用力会改变船舶在水平面内的三自由度运动，即纵荡、横荡和首摇运动。拖轮的推力或拉力可以直接改变船舶的纵荡速度，使船舶加速、减速或保持匀速运动。拖轮施加的横向力会导致船舶产生横荡运动，使船舶偏离原来的航线。拖轮产生的转船力矩则会引起船舶的首摇运动，改变船舶的航向。在船舶靠泊过程中，通过调整拖轮的位置和推力方向，可以使船舶产生合适的横荡和首摇运动，从而准确地靠泊到指定位置。若需要船舶向右转向，可使右舷拖轮施加一个较大的推力，产生一个向右的转船力矩，同时控制左舷拖轮的推力，使船舶在转向的同时保持一定的前进速度，避免产生过大的横荡。这些相互作用力的大小和方向需要根据船舶的操纵需求、初始状态以及外界环境条件等因素进行精确控制，以实现安全、高效的船舶操纵。2.3Matlab仿真技术在船舶操纵中的应用原理2.3.1Matlab软件简介Matlab（MatrixLaboratory）是MathWorks公司推出的一款高性能的数值计算和可视化软件，在科学计算、工程设计、数据分析等众多领域都有着广泛的应用。它以矩阵作为基本的数据单位，提供了强大的数值运算功能，涵盖线性代数、微积分、数值分析、优化计算等多个数学领域，能够轻松实现各种复杂的数学算法和模型。用户可以通过简单的矩阵运算指令，快速完成矩阵的加、减、乘、除、求逆、特征值计算等操作，大大提高了计算效率和编程便利性。Matlab还具备卓越的图形绘制和可视化能力，能够将计算结果以直观的图形方式展示出来，包括二维图形（如折线图、散点图、柱状图等）和三维图形（如三维曲面图、立体图等）。通过可视化功能，用户可以更清晰地理解数据的特征和规律，发现数据之间的内在联系，从而为数据分析和决策提供有力支持。在船舶操纵仿真中，可以利用Matlab绘制船舶的运动轨迹图、速度时间曲线、航向变化曲线等，直观地展示船舶在不同操纵条件下的运动状态。此外，Matlab拥有丰富的工具箱，这些工具箱是针对特定领域的专业函数库，涵盖了信号处理、控制系统、图像处理、通信系统、金融分析等多个领域。在船舶操纵研究中，常用的工具箱包括Simulink工具箱、控制系统工具箱、优化工具箱等。Simulink工具箱提供了基于图形化建模的仿真环境，方便用户搭建各种动态系统模型；控制系统工具箱则提供了一系列用于控制系统设计、分析和仿真的函数和工具，有助于研究船舶操纵系统的性能和稳定性；优化工具箱可用于求解各种优化问题，如船舶操纵参数的优化选择，以提高船舶的操纵性能。Matlab还具有良好的开放性和可扩展性，用户可以根据自己的需求编写自定义函数和脚本，与其他软件进行数据交互和集成，进一步拓展其应用范围。2.3.2Simulink工具箱在船舶操纵仿真中的应用Simulink是Matlab的一个重要附加产品，它是一个基于图形用户界面（GUI）的多域仿真和模型设计工具，专门用于对动态系统进行建模、仿真和分析。在船舶操纵仿真中，Simulink发挥着至关重要的作用。利用Simulink进行船舶操纵仿真的第一步是搭建船舶运动模型。用户可以在Simulink的图形界面中，从各种功能模块库中拖曳所需的模块，如积分器模块、加法器模块、乘法器模块、传递函数模块等，并按照船舶运动数学模型的结构和逻辑关系，将这些模块连接起来，构建出完整的船舶运动模型。对于MMG模型，需要将船体水动力模块、螺旋桨模块、舵模块以及拖轮作用模块等进行合理连接，准确模拟船舶在各种力和力矩作用下的运动。在搭建模型的过程中，需要根据实际情况设置各个模块的参数。对于船舶的质量、惯性矩、水动力系数、螺旋桨参数、舵参数等，都需要根据船舶的设计资料和实验数据进行准确赋值。通过设置不同的参数值，可以模拟不同类型船舶在不同工况下的操纵性能。对于大型油轮和集装箱船，由于其质量和惯性较大，在设置参数时需要充分考虑这些因素对船舶操纵的影响。完成模型搭建和参数设置后，就可以进行仿真实验。在仿真过程中，用户可以设置仿真的起始时间、终止时间、步长等参数，控制仿真的运行过程。Simulink会根据模型的结构和参数，按照设定的仿真步长，逐步计算船舶在每个时间点的运动状态，如位置、速度、航向等，并将仿真结果保存下来。用户可以通过示波器模块、XYGraph模块等可视化工具，实时观察船舶的运动轨迹和各种运动参数随时间的变化情况，也可以在仿真结束后，对保存的仿真数据进行进一步的分析和处理。Simulink还支持对模型进行优化和调试。用户可以通过调整模型的结构和参数，观察仿真结果的变化，以优化模型的性能和准确性。当模型出现错误或异常时，Simulink提供了丰富的调试工具，如断点设置、单步执行、数据监控等，帮助用户快速定位和解决问题，确保仿真结果的可靠性。三、双拖轮协助船舶操纵运动数学模型的建立3.1考虑的因素3.1.1浅水对船舶操纵性的影响当船舶从深水航行至浅水或在浅水中航行时，航道水深与船舶吃水之比小于1.5时，便会发生浅水效应，这会导致船舶操纵性能发生显著变化。浅水对船舶操纵性的影响是多方面的，其根本原因在于船体周围的流场和兴波情况发生了改变。在深水中，船舶航行时船首和船尾部分水的流动具有三维空间内流动的特点，船首处水斜向（既向两侧，又向下方）向后运动，船尾处水斜向（即由两侧向纵中剖面，又向上方）向后运动。然而，在浅水中航行时，由于船底与河床之间的空间受到限制，船首或船尾部分水的流动由三维空间内的流动变为向两侧或由两侧同时向内的二维平面式的流动。这种流场的变化导致了船体周围水压分布的改变，进而产生了一系列影响。船体底部与河床之间形成狭窄的水道，使得船底流速增大。由于水的粘性，在河床及船底均会形成边界面，过水断面进一步减小，船底流速更加增加。这种由于水深受到限制使船体与水的相对速度较深水情况有所增大的速度称为回流速度。回流速度的存在不仅导致船底流速增加、压力降低，使船体下沉、吃水增加，还使得附加质量及附加惯性矩较之深水有明显增加。而且，由于船底和河床边界层厚度均自船首向船尾逐渐增加，船尾部过水断面较船首处小，流速增加更大，压力下降更甚，故船尾下沉较船首大，产生尾倾现象。回流速度的存在使浅水船周围的流速比深水船大，且其舷侧湿面积因船体下沉而增加，所以摩擦阻力增大。同时，因水流与船体的相对速度增大，压力下降亦大，首尾压力差将增大，并且船尾与河床的间隙小，易于发生涡流，因此涡流阻力也要增大，船在浅水中航行时其粘性阻力将增加。粘性阻力的增加使得船舶航行时需要克服更大的阻力，从而导致船速下降。在浅水区域，舵力也会受到影响。浅水中，舵叶周围的水流发生变化，伴流、涡流增加，这使得舵力下降。另一方面，由于相同转速时浅水中船速下降，增大了螺旋桨的滑失，又在一定程度上提高了舵力。但总的结果是舵力下降实际上并不大。然而，在浅水中舵效却明显变差，这是因为浅水中回转阻力大大增加，旋回性指数K大大减小，而追随性指数T减小的幅度比K来得慢，舵效指数K/T因而减小。船舶的旋回性和航向稳定性在浅水中也会发生改变。浅水中航行，舵力下降不大，舵力转船力矩下降也不大，但船舶旋回阻矩及虚惯矩均有较大增加，其中旋回阻矩的增加幅度更大。从船舶旋回性和追随性指数K、T来分析，K=a/b，T=lG/b，当船舶从深水进入浅水中，K、T同时减小，所以船舶旋回性变差，而追随性、航向稳定性变好。从大型油轮在不同水深与吃水之比H/d条件下，满舵35°旋回的试验可知，水深变浅，旋回初径增大，纵距也有所增加，但纵距增加率远低于旋回初径增加率，由此可见，浅水中旋回圈各要素并非整体性增加。根据多数试验结果表明，一般认为，H/d=4左右时，旋回圈开始受到浅水的影响；当H/d=1.25时，旋回圈约比深水中增大70%左右。3.1.2风、流等环境因素的作用风、流等环境因素对船舶运动有着重要的影响，在建立双拖轮协助船舶操纵运动数学模型时，必须充分考虑这些因素。风压力的计算及影响：风对船舶的作用力主要包括风压力和风力矩。风压力的大小与风速、风向、船舶受风面积以及船舶的航向等因素密切相关。目前，计算风压力的方法有多种，常见的有经验公式法和计算流体力学（CFD）方法。经验公式法是根据大量的实船试验数据总结出来的，具有一定的工程实用性。Isherwood公式通过考虑船舶的特征参数和风速来计算风压力，其表达式为：F_w=C_w\rho_aA_wV_w^2，其中F_w为风压力，C_w为风压力系数，\rho_a为空气密度，A_w为船舶迎风面积，V_w为风速。风压力系数C_w会受到船舶的形状、大小、装载状态以及风向等因素的影响，一般需要通过实验或经验数据来确定。当船舶受到风压力作用时，会产生平移和转动运动。若风从正横方向吹来，船舶会受到一个横向的风压力，使船舶产生横荡运动；同时，由于风压力作用点与船舶重心不重合，会产生一个转船力矩，使船舶发生首摇运动，改变船舶的航向。在低速航行或者干舷较高的时候，风对船舶操纵的影响特别显著，可能导致船舶操纵困难进而产生危险。当船舶在狭窄水道中航行时，如果受到较强的横风作用，可能会使船舶偏离预定航线，甚至与岸边或其他障碍物发生碰撞。流作用力的计算及影响：流对船舶的作用力主要包括流压力和流力矩。流压力的计算与水流速度、流向、船舶与水流的相对速度以及船舶的形状等因素有关。通常可以采用势流理论或基于实验数据的经验公式来计算流压力。根据势流理论，流压力可以通过求解船舶周围的流场势函数来得到。流作用力会改变船舶的运动状态。当船舶顺流航行时，水流会推动船舶前进，使船舶的实际航速增加；当船舶逆流航行时，水流会阻碍船舶前进，使船舶的实际航速降低。水流还会对船舶产生横向作用力，当船舶与水流存在一定夹角时，水流会使船舶产生横荡运动。在弯曲的河道中，水流的流速和流向会发生变化，船舶受到的流作用力也会随之改变，这对船舶的操纵提出了更高的要求。如果船舶在进入弯道时，未能及时调整航向和速度，可能会因受到水流的作用而偏离航道，发生搁浅或碰撞事故。流作用力还会与风压力、拖轮作用力等相互叠加，共同影响船舶的运动，使得船舶操纵更加复杂。3.2模型构建3.2.1船体水动力模型的建立船体在水中运动时，会受到复杂的水动力作用，这些力对船舶的操纵性能有着关键影响。在建立船体水动力模型时，通常将船舶在水平面内的运动分解为纵荡、横荡和首摇三个自由度的运动，分别考虑这三个方向上的水动力。对于纵荡方向的水动力，主要包括摩擦阻力、形状阻力和兴波阻力等。摩擦阻力是由于船体表面与水之间的摩擦力产生的，可根据平板摩擦阻力公式进行估算，并结合船体的湿表面积进行修正。形状阻力则与船体的横截面积、流线型等因素有关，一般通过经验公式或数值计算方法来确定。兴波阻力是由于船舶航行时产生的波浪所引起的阻力，其计算较为复杂，通常采用基于理论或实验数据的半经验公式。某船的纵荡水动力F_{u}可表示为：F_{u}=-(F_{f}+F_{s}+F_{w})，其中F_{f}为摩擦阻力，F_{s}为形状阻力，F_{w}为兴波阻力。横荡方向的水动力主要由船体的横剖面形状、漂角以及水的粘性等因素决定。当船舶存在漂角时，船体侧面会受到水的作用力，产生横荡水动力。横荡水动力F_{v}的计算公式可以表示为：F_{v}=\rhoV^2LdC_{yv}(\beta)，其中\rho为水的密度，V为船舶的速度，L为船长，d为吃水，C_{yv}(\beta)为横荡水动力系数，它是漂角\beta的函数。横荡水动力系数C_{yv}(\beta)通常通过实验或数值模拟来确定，它反映了横荡水动力随漂角的变化规律。首摇方向的水动力主要来自于船体的回转运动以及各方向水动力对船舶重心产生的力矩。首摇水动力矩M_{r}可表示为：M_{r}=\rhoV^2L^3C_{nr}(\beta)，其中C_{nr}(\beta)为首摇水动力矩系数，同样是漂角\beta的函数。首摇水动力矩系数C_{nr}(\beta)也需要通过实验或数值模拟来获取，它对于研究船舶的转向性能和航向稳定性具有重要意义。在实际计算中，还需要考虑附加质量和附加惯性矩的影响。附加质量是由于船舶运动导致周围水的加速而产生的等效质量，附加惯性矩则是考虑水对船舶转动的惯性影响。它们会改变船舶的运动特性，在模型中通常通过引入相应的附加质量系数和附加惯性矩系数来进行计算。某船的附加质量在纵荡方向为m_{u}，横荡方向为m_{v}，首摇方向为I_{r}，则在建立船体水动力模型时，需要将这些附加质量和附加惯性矩考虑进去，以更准确地描述船舶的运动。3.2.2螺旋桨和舵的模型构建螺旋桨和舵是船舶操纵系统中的关键部件，它们产生的力和力矩对船舶的运动起着决定性作用，因此需要建立准确的模型来描述它们的工作特性。螺旋桨模型：螺旋桨的主要作用是产生推力，推动船舶前进。螺旋桨的推力与螺旋桨的转速、螺距、直径以及进速等因素密切相关。通常通过螺旋桨的敞水性征曲线或经验公式来计算螺旋桨的推力。螺旋桨的敞水性征曲线是通过模型试验得到的，它反映了螺旋桨在不同进速系数下的推力系数和转矩系数。进速系数J的计算公式为：J=\frac{V_A}{nD}，其中V_A为螺旋桨进速，n为螺旋桨转速，D为螺旋桨直径。根据螺旋桨的敞水性征曲线，可以得到推力系数K_T与进速系数J的关系，进而计算出螺旋桨的推力T：T=\rhon^2D^4K_T(J)。螺旋桨的伴流效应也会对船舶的操纵性能产生影响。伴流是指螺旋桨旋转时带动周围水一起运动形成的水流，它会使螺旋桨的实际进速与船舶的前进速度不同。在计算螺旋桨推力时，需要考虑伴流的影响，通常通过引入伴流分数w来修正进速。螺旋桨的实际进速V_{A}'可表示为：V_{A}'=V(1-w)，其中V为船舶的前进速度。伴流分数w与船舶的类型、船体形状、螺旋桨位置等因素有关，一般通过经验公式或实验数据来确定。舵模型：舵是船舶实现转向的重要装置，通过改变舵角可以产生舵力和舵力矩，从而实现船舶的转向。舵力和舵力矩的大小与舵角、舵面积、舵速以及水流对舵的作用力等因素有关。常用的舵水动力计算方法有基于薄翼理论的方法和经验公式法。基于薄翼理论的方法将舵视为薄翼，通过求解绕翼的流场来计算舵力和舵力矩；经验公式法则是根据大量的实验数据总结得到的，具有一定的工程实用性。某船的舵力F_R可通过经验公式计算：F_R=\frac{1}{2}\rhoV_r^2A_RC_{L\alpha}\alpha，其中\rho为水的密度，V_r为舵速，A_R为舵面积，C_{L\alpha}为升力系数，\alpha为舵角。舵速V_r与船舶的前进速度和伴流有关，通常可表示为：V_r=V(1-w_r)，其中w_r为伴流系数，与螺旋桨的伴流情况相关。舵力矩M_R则可表示为：M_R=F_Rl_R，其中l_R为舵力作用点到船舶重心的距离。螺旋桨和舵之间还存在相互干涉关系。螺旋桨的尾流会使舵表面的水流速度和压力分布发生变化，从而改变舵力的大小和方向；舵的偏转也会影响螺旋桨的工作状态，如改变螺旋桨的进速和伴流。在建立模型时，需要考虑这些干涉效应，以准确描述船舶的操纵运动。3.2.3拖轮模型的融入在双拖轮协助船舶操纵运动数学模型中，拖轮模型的融入是关键环节，它直接关系到模型对实际操纵情况的模拟准确性。拖轮通过对船舶施加推力，改变船舶的运动状态，实现协助操纵的目的。拖轮的推力计算是建立拖轮模型的基础。拖轮的推力大小受到多种因素的影响，如拖轮的主机功率、螺旋桨性能、拖轮与船舶的连接方式以及外界环境条件等。通常可以采用基于螺旋桨理论的计算方法或利用拖轮性能曲线来计算拖轮推力。基于螺旋桨理论的计算方法从螺旋桨的工作原理出发，考虑螺旋桨的几何参数、运行参数以及流体的物理性质来计算推力。某拖轮螺旋桨直径为D_t，转速为n_t，进速为V_t，根据Betz理论，其推力T_t的计算公式可以表示为：T_t=k_1\rho\pi(\frac{D_t}{2})^2n_t^2V_t+k_2\rho\pi(\frac{D_t}{2})^4n_t^4，其中\rho为水的密度，k_1和k_2为经验系数，需要通过实验或经验数据来确定。利用拖轮性能曲线的方法则是根据拖轮在设计和制造过程中通过实验测量得到的不同工况下的推力性能曲线，根据拖轮的当前运行状态，从性能曲线中直接查取对应的推力值。在实际应用中，还需要考虑环境因素对拖轮推力的影响，如水流、风力等。当拖轮协助船舶操纵时，拖轮与船舶之间通过缆绳或顶推装置等连接方式产生相互作用力。在模型中，需要考虑连接点的力和力矩对船舶运动的影响。连接点的力包括拉力或推力，力矩则是由于力的作用点与船舶重心不重合而产生的。在顶推方式下，拖轮对船舶施加一个向前的推力，同时由于顶推点与船舶重心不在同一位置，会产生一个使船舶绕重心转动的力矩。若拖轮顶推点在船舶重心前方且偏左，当拖轮施加推力时，除了使船舶向前运动外，还会使船舶产生一个向左的首摇力矩，导致船舶向左转向。在拖带方式下，缆绳对船舶施加拉力，拉力的大小和方向取决于拖轮的推力以及缆绳的布置情况。缆绳的拉力不仅会使船舶产生平移运动，还会因为缆绳与船舶轴线的夹角而产生一个转船力矩。当缆绳与船舶轴线成一定角度时，拉力的横向分力会使船舶产生横荡运动，而纵向分力则影响船舶的前进或后退速度，同时由于力臂的存在，会产生一个使船舶绕重心转动的首摇力矩。在模型中，将拖轮的推力和连接点的力和力矩作为外力引入到船舶运动方程中，与船体水动力、螺旋桨推力、舵力等一起，共同决定船舶的运动状态。某船舶在双拖轮协助下的运动方程可以表示为：\begin{cases}m(\dot{u}-vr)=F_{u}+T_{t1x}+T_{t2x}+F_{p}+F_{R}\\m(\dot{v}+ur)=F_{v}+T_{t1y}+T_{t2y}\\I_z\dot{r}=M_{r}+M_{t1}+M_{t2}+M_{p}+M_{R}\end{cases}其中，m为船舶质量，I_z为船舶绕z轴的惯性矩，u、v、r分别为船舶在纵荡、横荡和首摇方向的速度分量，F_{u}、F_{v}、M_{r}分别为船体在纵荡、横荡和首摇方向的水动力和力矩，T_{t1x}、T_{t1y}、M_{t1}和T_{t2x}、T_{t2y}、M_{t2}分别为两个拖轮在相应方向上对船舶施加的力和力矩，F_{p}、M_{p}为螺旋桨的推力和力矩，F_{R}、M_{R}为舵力和舵力矩。通过求解这个方程组，可以得到船舶在双拖轮协助下的运动状态。3.3模型验证3.3.四、基于Matlab的双拖轮协助船舶操纵运动仿真实现4.1仿真平台搭建4.1.1Matlab及Simulink环境配置在进行基于Matlab的双拖轮协助船舶操纵运动仿真之前，首先需要进行Matlab及Simulink环境的配置。Matlab软件的安装过程较为常规，从MathWorks官方网站获取安装包后，运行安装程序，按照安装向导的提示逐步完成安装。在安装过程中，用户可根据自身需求选择安装路径以及所需的组件。为确保后续仿真工作的顺利进行，建议安装完整的Matlab组件，包括基础模块以及常用的工具箱。安装完成后，启动Matlab软件，进入其工作界面。Simulink是Matlab的重要扩展模块，用于动态系统的建模、仿真和分析，在船舶操纵运动仿真中发挥着关键作用。Simulink工具箱通常在Matlab安装过程中会一并安装，如果未安装，可通过Matlab的附加功能管理器进行安装。打开Matlab软件后，点击“附加功能”选项卡，在弹出的菜单中选择“获取附加功能”，在搜索框中输入“Simulink”，找到Simulink工具箱后点击安装即可。完成Simulink工具箱的安装后，还需要进行一些参数设置，以满足船舶操纵运动仿真的需求。在Matlab命令窗口中输入“simulink”，打开Simulink库浏览器。点击“文件”菜单，选择“首选项”，在弹出的首选项对话框中，可对仿真的相关参数进行设置。对于仿真求解器的选择，应根据船舶运动模型的特点进行合理设置。若模型为连续系统，可选择变步长求解器，如ode45（Dormand-Prince），该求解器适用于大多数连续系统的仿真，具有较高的精度和稳定性；若模型包含离散部分，则可选择定步长求解器，如ode5（Dormand-Prince），并设置合适的步长。在“数据导入/导出”选项中，可设置仿真结果的保存方式和路径。为便于后续对仿真数据的分析和处理，建议将仿真结果保存为Matlab的数据文件格式（.mat），并指定一个易于查找的保存路径。还可根据需要设置其他参数，如仿真的起始时间、终止时间等，以确保仿真过程符合实际需求。4.1.2仿真模型的搭建步骤在完成Matlab及Simulink环境配置后，接下来即可在Simulink中搭建双拖轮协助船舶操纵运动的仿真模型。搭建过程主要包括船舶模型、拖轮模型以及环境因素模型的搭建。搭建船舶模型时，根据前文建立的船舶操纵运动数学模型，从Simulink库浏览器中拖曳所需的模块。积分器模块用于对船舶的运动状态进行积分，从而得到船舶的位置、速度等信息；加法器模块用于将船舶所受到的各种力和力矩进行叠加；乘法器模块用于实现力和力矩与船舶运动参数之间的乘法运算；传递函数模块则可用于描述船舶的水动力特性、螺旋桨特性以及舵的特性等。将这些模块按照船舶运动数学模型的结构进行连接，构建出完整的船舶模型。在连接过程中，需确保模块之间的信号流向正确，参数设置准确。对于拖轮模型的搭建，同样从Simulink库浏览器中选择相应的模块。拖轮的推力计算模块根据拖轮的主机功率、螺旋桨性能等参数，计算出拖轮的推力大小和方向；连接点力和力矩计算模块则根据拖轮与船舶的连接方式和位置，计算出连接点处对船舶施加的力和力矩。将拖轮模型与船舶模型进行连接，将拖轮的推力和连接点的力和力矩作为外力输入到船舶模型中，实现拖轮对船舶的协助操纵。环境因素模型主要包括风、流等因素的模拟。风模型通过风速模块和风向模块，输入风的速度和方向信息，再经过风压力计算模块，根据风压力计算公式，计算出风对船舶施加的力和力矩；流模型则通过水流速度模块和水流方向模块，输入水流的速度和方向信息，经过流作用力计算模块，计算出流对船舶的作用力和力矩。将风模型和流模型与船舶模型进行连接，将风、流的作用力和力矩作为外力输入到船舶模型中，以考虑环境因素对船舶操纵运动的影响。在搭建完各个子模型后，对整个仿真模型进行整合和调试。检查模型中各个模块的参数设置是否正确，模块之间的连接是否无误。通过设置不同的初始条件和输入信号，对模型进行初步的仿真测试，观察模型的运行情况和输出结果。若发现模型存在错误或异常，及时进行排查和修正，确保模型能够准确地模拟双拖轮协助船舶操纵运动的过程。4.2仿真参数设置4.2.1船舶参数设定船舶参数的准确设定对于仿真结果的准确性至关重要。在本次仿真中，选取一艘典型的大型商船作为研究对象，其主要参数如下：船长L=200m，船宽B=30m，吃水d=10m，船舶质量m=50000t，船舶绕z轴的惯性矩I_z=1.5\times10^8kg\cdotm^2。这些参数反映了船舶的基本物理特征，对船舶的操纵性能有着重要影响。船体水动力系数是船舶操纵运动数学模型中的关键参数，它们决定了船体在水中运动时所受到的水动力大小和方向。根据船舶的类型和相关经验数据，确定纵荡水动力系数C_{xu}=-0.1，横荡水动力系数C_{yv}=-0.5，首摇水动力矩系数C_{nr}=-0.05。这些系数的取值会受到船舶的形状、大小、装载状态以及水流条件等因素的影响，在实际应用中可能需要根据具体情况进行调整。螺旋桨参数也是影响船舶操纵性能的重要因素。螺旋桨直径D=6m，螺距比P/D=0.8，螺旋桨效率\eta=0.6。螺旋桨的推力与转速、螺距、直径以及进速等因素密切相关，通过这些参数可以计算出螺旋桨在不同工况下的推力，从而为船舶提供前进的动力。舵参数包括舵面积A_R=30m^2，舵展弦比\lambda=2，最大舵角\alpha_{max}=35^{\circ}。舵力和舵力矩的大小与舵角、舵面积、舵速以及水流对舵的作用力等因素有关，通过设置合适的舵参数，可以实现对船舶航向的有效控制。4.2.2拖轮参数设定拖轮的参数设定直接影响其协助船舶操纵的效果。在仿真中，选用两艘功率相同的拖轮，每艘拖轮的主机功率P=2000kW。主机功率决定了拖轮的动力输出能力，进而影响拖轮的推力大小。拖轮的推力与主机功率、螺旋桨性能等因素有关。根据拖轮的主机功率和螺旋桨参数，通过推力计算公式，可得拖轮的最大推力T_{max}=500kN。在实际操纵过程中，拖轮的推力可根据船舶的操纵需求进行调整。拖轮与船舶的连接位置对船舶的运动状态有着重要影响。在本次仿真中，设置一艘拖轮连接在船舶的首部，距离船舶重心x_1=-80m，另一艘拖轮连接在船舶的尾部，距离船舶重心x_2=80m。连接位置的不同会导致拖轮对船舶施加的力和力矩不同，从而影响船舶的纵荡、横荡和首摇运动。拖轮的推力方向也可以根据船舶的操纵需求进行调整。在船舶转向时，可通过调整拖轮的推力方向，使船舶产生合适的转船力矩，实现船舶的转向；在船舶靠泊时，可调整拖轮的推力方向，使船舶平稳地靠近码头。4.2.3环境参数设定风、流等环境因素对船舶操纵运动有着显著影响，因此需要合理设定环境参数。在仿真中，设定风的速度V_w=10m/s，风向与船舶航向的夹角\theta_w=45^{\circ}。风的速度和方向决定了风对船舶施加的力和力矩大小和方向，会使船舶产生横荡和首摇运动，影响船舶的航向和航速。水流速度V_c=2m/s，水流方向与船舶航向的夹角\theta_c=30^{\circ}。水流的存在会改变船舶与水的相对速度，从而影响船舶的操纵性能。顺流时，船舶的实际航速会增加；逆流时，船舶的实际航速会降低。水流还会对船舶产生横向作用力，使船舶产生横荡运动。将风、流等环境因素的参数输入到环境因素模型中，通过模型计算出风、流对船舶施加的力和力矩，再将这些力和力矩作为外力输入到船舶操纵运动模型中，以模拟实际环境条件下双拖轮协助船舶操纵的过程。在不同的仿真工况下，可以根据实际情况调整环境参数，研究环境因素对船舶操纵运动的影响规律。4.3仿真实验方案设计4.3.1不同协助方式的仿真实验为了深入研究双拖轮协助船舶操纵的效果，设计了不同协助方式的仿真实验。在实验中，主要考虑双拖轮的位置和推力组合对船舶操纵运动的影响。设置双拖轮在船舶的不同位置进行协助操纵。除了前文提到的首部和尾部连接方式外，还设置双拖轮均连接在船舶的一侧，距离船舶重心分别为x_3=-50m和x_4=50m；以及双拖轮呈对角连接，一艘连接在船舶的首部左侧，另一艘连接在船舶的尾部右侧。通过改变双拖轮的连接位置，可以分析不同位置对船舶纵荡、横荡和首摇运动的影响。在首部和尾部连接时，拖轮主要对船舶的纵荡和首摇运动产生影响，有助于船舶的加速、减速和转向；而在一侧连接时，拖轮会使船舶产生较大的横荡运动，可用于调整船舶的横向位置；对角连接则可以综合产生纵荡、横荡和首摇运动，实现更加复杂的操纵。调整双拖轮的推力组合。设置双拖轮推力相同，均为T=200kN；以及双拖轮推力不同，一艘推力为T_1=150kN，另一艘推力为T_2=250kN。通过改变推力组合，可以研究不同推力大小对船舶操纵运动的影响。当双拖轮推力相同时，船舶的运动较为平稳；而当推力不同时，船舶会产生一定的转船力矩，可用于实现船舶的转向或调整航向。在每种协助方式下，设置不同的仿真时间和初始条件，观察船舶的运动轨迹、速度、航向等参数的变化情况。通过对仿真结果的分析，总结不同协助方式下船舶操纵运动的规律和特点，为实际船舶操纵提供参考依据。4.3.2不同工况下的仿真实验除了研究不同协助方式的影响外，还设计了不同工况下的仿真实验，以全面评估双拖轮协助船舶操纵的性能。在实验中，主要设置船舶的不同初始速度和航向。设定船舶的初始速度分别为V_0=5m/s、V_0=8m/s和V_0=10m/s。初始速度的不同会影响船舶的惯性和运动状态，进而影响双拖轮协助操纵的效果。在相同的协助方式下，初始速度越大，船舶的惯性越大，双拖轮需要施加更大的力和力矩才能改变船舶的运动状态。船舶的初始航向分别为\psi_0=0^{\circ}、\psi_0=30^{\circ}和\psi_0=60^{\circ}。初始航向的不同会导致船舶与风、流的相对角度不同，从而使船舶受到的环境作用力不同。在不同的初始航向条件下，双拖轮需要根据船舶的实际情况调整推力和方向，以实现有效的协助操纵。在每种工况下，结合不同的协助方式进行仿真实验，分析船舶在不同初始条件和协助方式下的运动响应。通过对比不同工况下的仿真结果，研究初始速度和航向对双拖轮协助船舶操纵的影响，为船舶在不同工况下的安全操纵提供指导。还可以进一步考虑船舶的装载状态、外界环境的变化等因素，设计更加复杂的仿真工况，以提高仿真实验的真实性和可靠性。五、仿真结果分析与讨论5.1双拖轮不同协助方式下船舶运动特性分析5.1.1横移运动分析通过对不同协助方式下船舶横移速度和位移的仿真结果进行分析，发现双拖轮的位置和推力组合对船舶横移运动有着显著影响。当双拖轮分别连接在船舶的首部和尾部且推力方向垂直于船舶轴线时，船舶能够产生较为明显的横移运动。在这种协助方式下，首部拖轮的推力使船舶首部向一侧移动，尾部拖轮的推力则使船舶尾部向相反方向移动，从而实现船舶的横移。若首部拖轮推力为T_1=200kN，尾部拖轮推力为T_2=200kN，仿真结果显示，在初始时刻，船舶横移速度为0，随着拖轮推力的持续作用，横移速度逐渐增大，在100s时达到最大值v_y=0.5m/s，此时船舶的横移位移为y=25m。继续增加仿真时间，横移速度逐渐趋于稳定，横移位移则继续增加，在300s时，横移位移达到y=100m。当双拖轮均连接在船舶的一侧时，船舶的横移运动也较为明显，但与首部和尾部连接方式相比，横移速度和位移的变化规律有所不同。在这种情况下，双拖轮的合力使船舶向一侧加速横移，横移速度的增长较为迅速。若双拖轮推力均为T=250kN，在仿真开始后的50s内，横移速度就达到了v_y=0.6m/s，横移位移为y=15m。随着时间的推移，横移速度继续增大，在150s时达到最大值v_y=0.8m/s，此时横移位移为y=60m。之后，横移速度逐渐稳定，横移位移持续增加，在300s时，横移位移达到y=150m。对角连接方式下，船舶的横移运动相对较为复杂，横移速度和位移的变化受到双拖轮推力的大小和方向以及船舶初始状态的影响。在某些情况下，对角连接方式可以使船舶在实现横移的同时，还能产生一定的首摇运动，从而实现更加灵活的操纵。若一艘拖轮连接在船舶的首部左侧，推力为T_3=150kN，另一艘拖轮连接在船舶的尾部右侧，推力为T_4=250kN，在仿真过程中，船舶不仅会产生横移运动，还会出现首摇运动。在100s时，横移速度为v_y=0.4m/s，横移位移为y=20m，同时首摇角度为\psi=10^{\circ}。随着时间的推移，横移速度和位移继续变化，首摇角度也会相应改变，在300s时，横移位移达到y=80m，首摇角度为\psi=30^{\circ}。5.1.2转向运动分析在双拖轮协助船舶转向运动方面，转向角速度和回转半径是衡量转向性能的重要指标。当双拖轮以不同方式协助船舶转向时，这两个指标会发生明显变化。在双拖轮分别位于船舶首尾且推力方向与船舶轴线成一定夹角的协助方式下，船舶能够产生较大的转向角速度。当首部拖轮推力与船舶轴线夹角为\theta_1=30^{\circ}，推力大小为T_1=200kN，尾部拖轮推力与船舶轴线夹角为\theta_2=-30^{\circ}，推力大小为T_2=200kN时，船舶的转向角速度在仿真开始后迅速增大。在50s时，转向角速度达到最大值\omega=0.05rad/s，此时船舶的回转半径为R=100m。随着转向的进行，转向角速度逐渐减小，回转半径则逐渐增大，在150s时，转向角速度降至\omega=0.03rad/s，回转半径增大到R=150m。若双拖轮均位于船舶一侧且推力方向相反，船舶也能实现转向，但转向角速度相对较小，回转半径较大。当双拖轮推力均为T=200kN，且方向相反时，在仿真开始后的100s内，转向角速度逐渐增大，在100s时达到最大值\omega=0.03rad/s，此时回转半径为R=200m。随着时间的推移，转向角速度逐渐稳定在\omega=0.02rad/s左右，回转半径则保持在R=250m左右。不同的转向协助方式对船舶的转向性能有着显著影响。在实际船舶操纵中，需要根据具体情况选择合适的拖轮协助方式，以实现高效、安全的转向操纵。如果船舶需要快速转向，应选择能够产生较大转向角速度的协助方式；而在对转向精度要求较高的情况下，则需要综合考虑转向角速度和回转半径，选择合适的拖轮配置和推力组合。5.2不同工况对双拖轮协助操纵效果的影响5.2.1船舶初始速度的影响船舶初始速度对双拖轮协助操纵的难易程度有着重要影响。通过仿真实验发现，当船舶初始速度较低时，双拖轮能够较为容易地改变船舶的运动状态，实现精确的操纵。在初始速度V_0=5m/s的工况下，双拖轮以首部和尾部连接且推力均为T=200kN的方式协助船舶转向，船舶在100s内就能够完成90^{\circ}的转向，转向过程较为平稳，横移和纵荡的位移也在可控范围内。随着船舶初始速度的增加，双拖轮协助操纵的难度逐渐增大。当初始速度提高到V_0=10m/s时，同样的双拖轮协助方式下，船舶完成90^{\circ}转向所需的时间延长至150s，且在转向过程中，船舶的横移和纵荡位移明显增大，操纵的精确性受到一定影响。这是因为初始速度越大，船舶的惯性越大，双拖轮需要施加更大的力和力矩才能改变船舶的运动状态。在较高初始速度下，双拖轮的推力对船舶运动状态的改变相对较小，需要更长的时间和更大的推力才能实现相同的操纵效果。当初始速度达到V_0=15m/s时，即使双拖轮加大推力至T=300kN，船舶完成90^{\circ}转向仍需要200s以上，且在操纵过程中，船舶的运动响应变得较为迟缓，容易出现超调现象，增加了操纵的风险。5.2.2风、流等环境因素的影响风、流等环境因素对拖轮协助船舶操纵运动轨迹有着显著影响。在有风的情况下，船舶会受到风压力的作用，从而改变其运动轨迹。当风速为V_w=10m/s，风向与船舶航向夹角为\theta_w=45^{\circ}时，若双拖轮以常规方式协助船舶航行，船舶的实际运动轨迹会向风的来向偏移。在没有风的情况下，船舶在双拖轮协助下能够沿着预定的直线轨迹航行，但在有风的情况下，船舶的轨迹会发生弯曲，在1