耐波性预报的数值方法
船舶耐波性能实验——阻尼系数测量
船舶耐波性能试验—阻尼系数测量试验学生姓名：学号：学院：船舶与建筑工程学院班级：指导教师：一、船模横摇试验的目的上风浪中航行最易发生横摇，而且横摇的幅度较大，不仅影响船员生活和工作的各个方面，严重的横摇还会危及船舶的安全乃至倾覆失事。
因此，在有关耐波性的研究中，首先关注的是要求设计横摇性能优良的船舶。
由于船舶在波浪中横摇运动的复杂性，理论计算尚未达到可用于实际的程度，因而模型试验是目前预报船舶横摇最可靠的方法。
本教学试验由下列两部分组成，即：1．船模在静水中的横摇衰减试验，目的是确定船的固有周期以及作用在船体上的水动力系数，如附连水惯性矩及阻尼系数等。
据此可根据线性运动方程计算船舶在风浪中的横摇频率响应曲线。
2．船模在规则波中的横摇试验，目的是确定船的横摇频率响应函数，可用于预报船舶在中等海况下的横摇统计特性，对于高海况的预报数值则偏高，这是由于非线性影响的缘故。
二．实验原理通过《船舶原理》课程的学习，我们知道船舶的横摇运动方程可以表示为：式中，表示横摇角、横摇角速度、横摇角加速度；Ixx’表示船舶在水中的横摇惯性矩，等于船舶在空气中的横摇惯性矩Ixx 与船舶在水中的横摇附加惯性矩之和；N为阻尼力矩系数；D为排水重量；h为横稳性高度；αm0为有效波倾；ω为波浪圆频率。
引入横摇衰减系数γ和横摇固有（圆）频率ωФωФ2=Dh/Ixx’横摇运动方程可以写成:静水中自由横摇考虑船舶在初始时刻浮于静水面上，并伴有一个静横倾角φ0，但不受波浪的作用，该船舶随后将作自由横摇运动，其表达式可以写成式中，无因次衰减系数μ和相位超前角β为自由横摇幅值随时间成指数规律衰减，而横摇角随时间成余弦变化规律。
余弦函数的周期为2π，当每增加2π时，横摇完成一个摇摆，对应的时间间隔为自由横摇周期TФ’，即：ωФ’ TФ’=2π或上式中的表示水阻尼对横摇周期的影响，实际上阻尼对周期的影响很小。
若不考虑水的阻尼，则=0，式（6）对应的自由横摇周期即为横摇固有周期。
FPSO耐波性运动与液舱晃荡全耦合数值模拟-欢迎光临万德成教授研究
图3
释放横荡自由度下液舱自由面情况
-4-
第六屆海峽兩岸水動力學研討會
图4
液舱壁抨击压力点布置及压力时历曲线
同时为了分析液舱的减摇效应，我们考虑船舶所受力矩与液舱壁上的压力的时历曲线， 如图 5 右侧所示。图 5 中横摇力矩数值均放大了 50 倍。可以看出，作用在舱壁法向为坐标轴 正向的力出现在力矩使得船舶向坐标轴负向旋转的时刻，因此，在一定程度上使得船舶横摇 的角加速度减小，从而减少了横摇。因为液舱中孔道的存在，液舱内部与船体表面为一个整 体，因此船体所得到的力矩考虑到了液舱内流体对舱壁的压力而产生的力矩。如图 5 左侧所 示，可以看出，零冲水率下的船舶所受的横摇力矩最大，因为没有液舱中流体的存在给予减 摇；其次是放开横荡自由度下的横摇力矩，固定横荡自由度的横摇力矩最小。这也和图 2 中横 摇运动幅值相吻合。 而由图 2 可以看出，放开横荡后船舶横摇幅值增大，因此，我们观察船舶所受横摇力矩 与液舱压力的时历曲线，如图 6 所示。其中横摇力矩数值放大 33 倍。可以看出，在固定横荡 运动时，如图 6 左所示，第二组舱壁抨击压力对横摇力矩有明显的影响作用，且出现在与横摇 力矩方向相反的地方，阻碍横摇进一步增大；而第三组舱壁抨击压力亦出现在与横摇力矩方 向相反的地方，虽然数值没有第二组抨击压力大，但依然有阻碍横摇的作用；而当放开横荡 运动时，如图 6 右所示，第二组和第三组舱壁抨击压力均出现在与横摇力矩方向相同的地方， 即与固定横荡自由度相比，加剧了船体的横摇运动。
-1-
第六屆海峽兩岸水動力學研討會
液舱晃荡在放开横荡自由度下的液舱的非线性自由面，Nasar 等人[11]采用实验方法将液舱固 定在驳船中，放开横摇垂荡和横荡自由度，观察其耦合现象。这些研究均表明，水平面运动 将会对液舱产生激励，且会使液舱产生剧烈晃荡，并改变船舶运动性能。 本文基于 Zhuang 和 Wan[6]的研究基础上，放开横荡自由度，并与之前所做三自由度运动 情况进行对比，观察横荡与横摇、垂荡和纵摇耦合情况下船舶内外流场的耦合情况。同时检 测了舱内的抨击压力，并与三自由度下的液舱抨击压力进行对比，分析舱内流体晃荡情况。 最后通过比较船舶运动与抨击压力的相位差，分析液舱晃荡对于船舶的减摇作用。
耐波性预报的数值方法
耐波性预报的数值方法随着海洋气象学的发展，耐波性预报成为海洋气象学研究的热点之一，其涵盖的范围广泛，包括海上工程、近岸海域和内陆水域等。
耐波性预报主要用于确定海浪对于结构物的影响，其精度和准确度对于海洋结构物的设计和施工来说至关重要。
本文将介绍一种常用的耐波性预报数值方法：波浪谱法。
波浪谱法是一种通过数学模型预测海浪高度的方法。
它基于统计学原理，将海浪高度看作频率的随机变量，利用矩阵特征值问题将其转化为谱表示，从而通过计算海浪频带的谱成分来预报海浪的高度。
波浪谱法可以分为线性和非线性两种方法。
线性波浪谱法适用于海浪高度比较小的情况，其核心公式为：$S(f) = \frac{4\pi^2}{g^2}f^{-5}e^{-\frac{5}{4}(\frac{f_p}{f})^4}$其中，$S(f)$是波浪谱，$f$是频率，$f_p$是波峰频率，$g$是重力加速度。
该公式基于对于波浪动能的能量守恒假设，即使在一定程度上违反了实际水文情况的非线性效应，但其计算精度较高，在近岸海域和海上工程的预报中得到了广泛应用。
非线性波浪谱法适用于海浪高度较大的情况，其核心公式包括：$S(f) = A(f) \times \delta(f-f_p) + B(f), A(f) =\frac{5}{16}\frac{H_m^2}{f_p^4}\frac{e^{\frac{-5}{4}(\frac{f_p}{f})^4}}{f^5}$其中，$H_m$是波高标准差，$\delta(f-f_p)$是一个峰值函数，$B(f)$是定常波谱。
非线性波浪谱法计算过程较为繁琐，需要考虑多种非线性效应，如波浪不对称性、波浪剧烈变化和能量流失等，但其计算精度更高，能够更准确地预报海浪高度。
波浪谱法还有一种改进方法——时域法。
时域法将波浪谱分解成一系列时间域中的波浪组件，通过求解波浪运动方程确定波浪高度随时间的变化。
时域法计算复杂度较高，但它可以考虑更多的非线性效应和海浪的时变性，能够更准确地预报大浪和风浪耦合下的海浪高度。
SHIPFLOW软件MOTIONS模块简介
SHIPFLOW 软件MOTIONS模块简介船舶的耐波性能是表征其在海上航行能否维持其正常功能的能力，也是历来船舶设计者和使用者十分关心的问题。
关注的焦点在于保障船舶、货物、全体人员最根本的安全，以及运营的效率即油耗的情况。
反映在耐波性试验中需要关心的数据主要有如下几部分：1）船体的运动及加速度；2）波浪中的阻力或者波浪增阻；3）船体表面上整体或者局部的载荷。
随着仿真技术的发展，如何通过数值手段对耐波性结果进行预报是当下的船舶设计工程师十分关注的问题。
MOTIONS模块是SHIPFLOW软件自6.0版本添加的船舶运动分析专用模块，可用于计算船舶在规则波和不规则波中的运动和附加阻力，也包括在静水中的阻力、升沉和纵摇。
高效的求解方法耐波性数据求解可以通过多种不同程度的近似方法获取，不同求解方法在复杂性、计算时间及计算精度上都有所差异。
按照求解的复杂度递增的顺序可以将这些求解方法依次排序：传统切片法一＞局部非线性切片法一＞线性边界元（3D）法一＞非稳态RANS方法一＞大涡模拟（LES）—＞直接数值模拟（DNS）。
从非稳态RANS方法开始采用的是粘流计算，由于时间成本高，往往并不能适用于在工程实践；而前面的几种势流求解方法虽然计算速度快，但精度较低。
MOTIONS模块中采用势流、时域、完全非线性的边界元方法，旨在填补传统势流方法与非稳态、粘流方法之间的空白。
与传统的势流方法相比，该方法具有更高的精度，同时也比现有的RANS方法具有更快的计算速度。
计算单个工况点，采用MOTIONS模块大概需要6~8小时（16核工作站），而采用RANS 方法（STAR-CCM+, Fine/MARINE, OpenFoam）在相同条件下则需要200~400小时。
条件 数值 航速[kn] 15.5 Fn 0.142浪向180°波长[Lambda/Lpp] 0.6, 0.7，…， 2.0 LCC2作为研究对象，对其波浪增阻进行计算，工况如下: 完善的功能MOTIONS 的计算域由一个自由液面及截断它的浮体组成，并具有如下特征: 1）假定计算域是一个更大的计算域的一部分；2）外部计算域是静水状态或是未受干扰的波浪流场状态。
舰船结构的耐波性能分析与优化设计
舰船结构的耐波性能分析与优化设计舰船结构的耐波性能是指舰船在海洋波浪环境下承受波浪载荷时的稳定性和安全性能。
耐波性能的分析与优化设计对于舰船的设计和建造具有重要意义。
本文将就舰船结构的耐波性能进行分析，并提出相应的优化设计方法。
1. 波浪载荷对舰船结构的影响波浪载荷是指充斥在舰船结构表面的波浪力，它对舰船结构的影响可表现为“浪荡力”和“迎浪力”。
浪荡力使得舰船产生横向和纵向的运动，而迎浪力则用于提供推进力。
2. 耐波性能分析方法2.1 线性理论分析线性理论是一种基本的理论模型，它假设波浪载荷与舰船的反应之间是线性关系。
通过线性理论分析，可以获得舰船在不同波浪条件下的波动响应和载荷分布。
2.2 非线性数值模拟非线性数值模拟是一种更加精确的耐波性能分析方法，它考虑了波浪载荷和舰船结构的非线性特性。
利用计算流体力学方法，可以模拟舰船在海上的真实工作环境，获得更准确的结果。
3. 耐波性能优化设计方法3.1 材料选择与强度优化在耐波性能优化设计中，选择合适的材料对于提高舰船结构的强度和刚度至关重要。
通过优化材料的特性，可以降低结构的振动幅值和应力，提高舰船的耐波性能。
3.2 结构形状优化结构形状对于舰船的耐波性能有着直接的影响。
通过优化结构的外形和布局，可以减小舰船与波浪的相互作用，降低波浪载荷和结构响应。
3.3 主动控制技术主动控制技术可以在一定程度上改善舰船的耐波性能。
通过采用合适的控制策略，如自适应控制和反馈控制，可以减小波浪力对舰船的影响，提高舰船的稳定性和操纵性能。
4. 耐波性能分析与优化设计案例以某型舰船为例，通过上述提到的耐波性能分析与优化设计方法，对舰船的结构进行分析和改进。
通过数值模拟和实际试验，获得了舰船在不同波浪条件下的耐波性能数据，并根据数据分析结果进行相应的结构优化设计。
5. 结论舰船结构的耐波性能是保证舰船在海上安全运行的重要因素。
本文通过分析耐波性能的影响因素，并介绍了线性理论分析和非线性数值模拟两种常用的分析方法。
船舶耐波性能及优化设计研究
船舶耐波性能及优化设计研究一、引言船舶耐波性能是衡量一艘船的重要指标之一，也是保障船舶海上安全的关键因素。
船舶在海上航行时，会面临各种波浪环境，船舶的耐波性能好坏决定了其航行的安全性和舒适性。
因此，研究船舶的耐波性能以及优化设计是一个具有重要意义的课题。
二、船舶耐波性能的影响因素船舶的耐波性能是由船体本身的设计和建造质量、船舶在海上的运动状态以及各种环境因素综合影响而形成的。
以下是影响船舶耐波性能的几个主要因素：1、船体结构设计船体结构设计是影响船舶耐波性能的最重要因素之一，它包括船体型面设计、船体尺寸比例、船舶结构强度等。
合理的船体结构设计有利于提高船舶的耐波性能。
2、载货量和船员配备随着船舶的载货量增大，船舶的大、小浪受力情况也会发生变化，会对船舶的耐波性能产生一定的影响。
而船员配备的多少也会影响船舶的艇身均衡状态和灵活性，从而影响船舶的耐波性能。
3、船舶在海上的运动状态船舶在海上的运动状态是受到风、浪、潮流等多种因素的综合影响而形成的，如航向、航速、波浪高度、波浪频率等。
这些因素会影响船舶的耐波性能。
4、波浪环境波浪环境是指船舶在海上遇到的波浪形态，包括波高、波浪频率、波浪周期等。
不同的波浪环境对船舶的耐波性能有不同的影响，需要对波浪环境进行全面的评估和分析。
三、优化船舶耐波性能的设计方法为了提高船舶的耐波性能，需要采取一些有效的优化设计方法。
以下是几种主要的设计方法：1、船体结构优化设计船体结构的优化设计可以通过数值模拟和实验测试两种方法来实现。
数值模拟主要利用计算机仿真技术分析和研究船体结构的力学性能，进行结构优化设计，实验测试则是通过对船舶模型进行真实的模拟试验，获取船体结构的力学特性数据。
2、锚泊安装和操纵策略优化对于大型船舶来说，锚泊安装和操纵策略的优化也是提高船舶的耐波性能的关键因素之一。
优化锚泊安装和操纵策略可以通过数值模拟和实验测试来实现。
3、减轻船舶载重量为了提高船舶的浮力和稳定性能，可以考虑减轻船舶的载重量。
基于数值方法的船舶纵向运动频域预报研究
基于数值方法的船舶纵向运动频域预报研究发表时间：2015-07-06T11:47:07.720Z 来源：《中国科技教育·理论版》2015年第5期供稿作者：潘文良刘可刁世伦[导读] 船舶在波浪中运动的耐波性研究，对船舶在海上的航行安全及武器使用，有着及其重要的作用。
潘文良刘可刁世伦 91550部队 116023 摘要船舶在海上的日常航行，不仅需要极短期的时域预报，也需要短期频域预报，而且还需要两者相结合的预报模式。
本文运用粘性流数值方法对WigleyIII型船模顶浪纵向运动进行了数值仿真,得到了船舶纵向摇荡的时历曲线,依据谱分析理论计算船舶摇荡谱,进而对船舶顶浪纵向运动进行频域预报。
预报结果同势流理论计算结果进行比对，效果良好。
关键词纵向运动频域预报数值仿真摇荡响应谱0引言船舶在波浪中运动的耐波性研究，对船舶在海上的航行安全及武器使用，有着及其重要的作用。
舰船在日常训练及战斗活动中，不仅需要极短期的时域预报，也需要短期频域预报，而且还需要两者相结合的预报模式。
舰船摇荡的频域预报，通常是基于舰船流体力学理论，采用概率统计方法给出预报。
本文主要是通过数值方法获取船舶摇荡的时历曲线，根据统计分析原理，对在长峰不规则波中顶浪纵向运动进行频域预报。
拓宽了船舶摇荡频域预报的新途径[1][2]。
1 船舶摇荡时历的获取本文通过粘性流方法对船舶纵向运动进行数值仿真计算的方法来获取船舶摇荡时历曲线，较好地克服了上述两种传统方法的缺点，同时提高了船舶摇荡时历的获取精度[3]。
1.1船模的选择本文选取WigleyIII型标准船模作为CFD数值模拟研究对象。
Wigley系列船模是用数学公式来表述型值的理想化船模。
若采用右手坐标系O，是纵向坐标，是横向坐标，向右舷为正，是垂向坐标，向下为正。
在此坐标系中，O点在船舯静水面上，无因次化的船型表达式为： (1)1.2网格划分粘性流数值水池试验中，网格划分的科学性和网格分布合理性对试验的精度和效率起着至关重要的作用。
基于三维势流理论的耐波性数值预报
基于三维势流理论的耐波性数值预报随着海洋工程事业的发展，对于浪尖高度、移动速度、滞留时间等海浪参数的准确预报愈加重要。
而海浪的复杂性，导致现阶段的描述方法难以满足实际需要。
本文将介绍一种基于三维势流理论的耐波性数值预报方法，并分析其在海洋工程中的应用。
三维势流理论是基于Bernoulli方程的理论，能够描述流体中复杂的运动，包括了流体的速度、压力和密度等参数，进而计算海浪的传播和变形。
在此理论下，我们可以通过求解海洋自由面上的流场，并根据动量守恒方程得到海洋波浪的传播方程。
本文所提供的数值预报方法，基于三维势流理论的基础上，利用有限元方法对其进行数值求解。
该方法具备一定的数值稳定性和精度，能够对复杂的海浪场进行预报，并为海洋工程提供可靠的模拟支持。
在预报过程中，我们可以根据实际需要选择不同的边界条件和初始条件，并通过预设网格使其适应不同的计算环境。
在求解过程中，我们需注意模型设置中的物理参数，包括海水密度、粘度、重力加速度等参数，以及模型中的阻力系数、回波系数等真实参数的设定，以确保预报结果的准确性和可靠性。
最终，基于三维势流理论的耐波性数值预报方法能够为海洋工程提供可靠且准确的波浪参数预报，包括浪高、周期、传播方向、波速等，进而对海洋工程的设计、施工和咨询提供重要参考依据。
预报结果的精度和可靠性，将在海岸防护、海上输油管道、港口码头、海上风电等领域的应用中得到验证。
海洋工程中需要预测的海浪参数包括浪高、周期、传播方向、波速等。
这些参数的准确预报对于海洋工程的设计、施工和咨询非常重要，因此需要通过数值模拟等方法进行计算和预测。
浪高是指波浪顶部相对于海平面的高度。
浪高受到多种因素的影响，包括风速、海洋流速度、水深等因素。
例如，在狂风暴雨的条件下，海浪往往会变得更加剧烈。
在海洋工程中，浪高的预报对于选择结构物的高度、决定摆放平台和绳索、决定船只的运行和供应等都非常重要。
周期是指波浪的时间周期，即相邻波峰和波谷之间的时间间隔。
穿浪双体船耐波性能预报方法研究
舶， 一个 船体 切 片 产生 的波 浪 主要 波及 顺 流 而 非 逆流 。所 以可 以从 船 首 开始 逐 步 解 决 顺 流 问题 。
这种理 论 已经应 用到具 有定 常 向前运 动 的尖首 型 船舶上 。类似 的方 法 由 Ｆｌｎｅ 应用 到兴 波阻 ａ ｉｓｎ ｔ 力增加 研究 。这 种理论 由 Ｙ ｕ ｇ 人在频 域 内研 ｅｎ 等
俗话 说 “ 上 无 风 三 尺 浪 ” 在 海 上 总 是 有 海 ，
浪 ， 员 长 期 在 海 上 生 活 ， 果 船 的 耐 波 性 比较 船 如 差 ， 员晕船 率 高 ， 多 重 要 工作 将 无 人 承 担 ， 船 很 直
主要 是 能正确 处 理 定 常 和 非定 常 流 的 交 互作 用 、 尾 封板 的效应 和计算 方 法 的数 值精 度 。全 三维 方 法计算会 花 费许 多时 间 。对 于瘦 长体 的高 速船 可
ｍ ｉ ｉｄｘｓｏ ａｅ ｐ ｇｑ ａｉｅ ｒ ａｃｐａ ｌ ａ ｐｏ ｃ ｄ ｓ ｎ ｅｃｐ ｄ ｅ ｅｉａ ｃ ｔ ｒ ａｎ ｎ ｅｅ ｆ ｅｋ ｅ ｉ ｕ ｌｉｓ ｅ ｃｅ ｔｂ ｔ ｒｊ ｔ ｅｉ －ｘ ｅｔ ｄ ｄｒｓ ｔ ｅａ Ｆ ｓ ｎ ｔ ａ ｅ ｅ ｇ ａ ｓｎ
特别 是军 船 ， 器 系 统 需 要 依据 船 的耐 波 性指 标 武
进行 整个 武器 系统 的设计 和选 型 。因此 ， 在设 计 阶段 如能 比较 准确地 确 定船 的耐 波性是 非 常必要
究船舶 运 动 。对 于高 速排 水 的船 舶 理论预 报 与试 验数据 一 致 。后 来 Ｚ ａ 人 提 出 了线性 时 域理 ｈｏ等
Ｓ ｕ ｙ ｏ ｅ ｉ ｔｏ ｅｈｏ ｏ ｏ ｙ ｏ ｅ ｋ ｅ ｎｇ ｆ ｒ Ｗ ＰＣ ｔ ｄ ｎ Ｐｒ ｄ ｃｉ ｎ Ｍ ｔ ｄ ｌ ｇ ｆＳ ａ ｅ ｐｉ ｏ
耐波性试验
40
4.1 试验方法及过程
本次横摇试验测量参数: 波浪周期、波高、横摇角
波长
1.56T
2
注意：测量波高（mm）、横摇角 为峰—峰值（度）
41
4.1 试验方法及过程
对波浪中试验中结果，国际水池议建议采 用下列无因次量表达： 纵摇幅值/波倾角 垂荡幅值/波幅
垂向加速度/波幅
首向相对运动/波幅
2Q/2ah-摇角/波倾角
2 T
e
2
g
V （遭遇频率）
45
4.1 试验方法及过程
典型的摇荡运动频响函数
0.25 横摇频响函数(超载尾楔板10度无翼） 0.2
0.15 零速 0.1
0.05
0 0 0.2 0.4 0.6 0.8 遭遇频率 1 1/s 1.2 1.4 1.6 1.8
51
感谢各位聆听
裴玉国 大连理工大学 船舶学院
52
19
2 试验设备
一种是固定的消波器，通常设置在造波机 对面。其作用是消除造波机传过来的波浪， 避免波浪反射回去 另一种是可动的边消波器，通常设置在水 池一侧或两侧。其作用是消除船行波，减少 等水时间。在进行波浪试验时应当将边消波 器提长到适当高度，以免对波浪形成干扰
20
2 试验设备
消波板
21
DS30型波高测量仪器
代入
v t Sr——斯特鲁哈尔数，当地惯性力与迁移惯性力的比值。
F 1 2 2 l v l
1
l l Sr vt vt
Sr 数以德国物理学家斯特劳哈尔（V.Strouhal）命名，他在研究风吹过金 属丝发出鸣叫声时创立此数。在研究不定常流动或脉动流时， Sr 数成为重
耐波性预报的数值方法
Ｋｅ ｙ ｗｏ ｒ ｄ ｓ ： ｓ ｅ ａ ｋ ｅ ｅ ｐ ｉ ｎ ｇ ； Ａ Ｑ ＷＡ； Ｓ Ｔ Ｆ ； ３ 一 Ｄ ｐ ｏ ｔ ｅ ｎ ｔ ｉ ａ ｌ ｌ ｆ ｏ ｗ ｔ ｈ ｅ ｏ ｒ ｙ
Ｎｕｍ ｅ ｒ ｉ ｃ ａ ｌ ｓ ｔ ｕ ｄｙ ｏ ｎ ｔ ｈｅ ｓ ｅ ａ ｋｅ ｅ ｐ ｉ ｎ ｇ ｐｒ ｅ ｄｉ ｃ ｔ ｉ ｏ ｎ
ＬＩ Ｙｕ － ｃ ｈｅ ｎ， ＷＵ Ｘｉ ， ＹＵＡＮ Ｘｉ ａ ｏ — ｌ ｏ ｎ ｇ， Ｈ ＡＮ Ｘｉ ａ ｏ — ｇ ｕ ａ ｎ ｇ
预报 。
关键 词 ： 耐 波性 ； Ａ Ｑ ＷＡ； Ｓ Ｔ Ｆ ； 三 维 势 流 理 论 中 图分 类 号 ： Ｕ ６ ６ １ ． ３ ２ 文献标 识码 ： Ａ
文章 编 号 ： １ ６ ７ ２—７ ６ ４ ９ （ ２ ０ １ ３ ） １ ０— ０ ０ １ ６— ０ ５ ｄ ｏ ｉ ： １ ０ ． ３ ４ ０ ４ ／ ｊ ． ｉ ｓ ｓ ｎ ． １ ６ ７ ２— ７ ６ ４ ９ ． ２ ０ １ ３ ． １ ０ ． ０ ０ ４
（ Ｎ ａ ｖ ｙ Ａ ｒ ｍ ｓ Ｃ ｏ ｍｍ ａ ｎ ｄ ｉ ｎ ｇ Ａ ｃ ａ ｄ ｅ ｍ ｙ ， Ｇ ｕ ａ ｎ ｇ ｚ ｈ ｏ ｕ ５ １ ０ ４ ３ ０ ， Ｃ ｈ ｉ ｎ ａ ） Ａ ｂ ｓ ｔ ｒ ａ ｃ ｔ ： Ｔ ｈ ｅ Ａ Ｑ ＷＡ， ａ ｃ ｏ ｍ ｍｅ ｒ ｃ ｉ ａ ｌ ｈ ｙ ｄ ｒ ｏ ｄ ｙ ｎ ａ ｍｉ ｃ ｓ ｏ ｆ ｔ ｗ ａ ｒ ｅ ， ｗ ｈ ｉ ｃ ｈ ｗ ａ ｓ ｂ ａ ｓ ｅ ｄ ｏ ｎ ３ － Ｄ ｐ ｏ ｔ ｅ ｎ ｔ ｉ ａ ｌ ｌ ｆ ｏ ｗ
Ａ Ｑ ＷＡ ｒ ｅ ａ ｓ ｏ ｎ ａ ｂ ｌ ｌ ｙ ｒ ｅ ｌ ｆ ｅ ｃ ｔ ｅ ｄ ｔ ｈ ｅ ｓ ｅ ａ ｋ ｅ ｅ ｐ ｉ ｎ ｇ ｐ ｅ ｆ ｒ ｏ ｒ ｍ ａ ｎ ｃ ｅ ｏ ｆ ｓ ｅ ｒ ｉ ｅ ｓ ６ ０ ， Ｓ Ｏ ｔ ｈ ａ ｔ ｉ ｔ ｃ ａ ｎ ｂ ｅ ａ ｐ ｐ ｌ ｉ ｅ ｄ ｔ ｏ ｗ ａ ｖ ｅ ｌ ｏ ａ ｄ
第十章 耐波性
第十章耐波性学习目标知识目标1.船舶摇摆对船舶的危害；2.船舶静水摇摆周期、波浪周期；波浪表观周期；3.船舶谐摇对船舶的危害及如何避免；4.减摇装置及减摇原理。
能力目标1.能计算船舶静水摇摆周期、波浪周期；波浪表观周期。
第一节概述船舶耐波性是船舶在波浪中运动特性的统称，它包括船舶在波浪中所产生的各种摇荡运动以及由这些运动引起的砰击、飞溅、上浪、失速、螺旋桨飞车和波浪弯矩变化等性能，直接影响船舶在风浪作用下维持其正常功能的能力，历来是船舶及其他海洋结构物的设计和使用者十分关心的问题。
研究船舶在波浪中产生的一系列运动，是学习船舶耐波性的首要目的。
了解了船舶在波浪中运动机理后，就可以探讨保证船舶在波浪中航行安全和维持其使用功能的措施。
在海上航行的船舶，像任何刚体一样，可以产生六个自由度的运动。
为了研究这些运动，通常采用以下右手坐标系(见图10-1)：图10-1 研究船舶运动的坐标系它是以船舶重心位置G为原点而固定于船体上的直角坐标系。
x 轴在中线面内，平行于基面，指向船首为正；z轴向上为正。
x、y和z 轴可近似认为是船体的三根惯性主轴。
Oxyz坐标系内船舶重心G沿坐标轴的直线运动船舶任意时刻的运动可以分解为在及船体绕三个坐标轴的转动。
在这些运动中又有单向运动和往复运动之分，因此共有12种运动形式，如图10-2所示。
习惯采用的名称见表10-1。
表 10-l 12种运动形式的习惯名称图 10-2船舶的运动形式图 10-3 船的遭遇浪向船的遭遇浪向见图10-3，横浪对横摇影响最大；顶浪顺浪产生纵摇；纵荡，垂荡。
耐波性是船舶在风浪中性能的总的反应，它包括以下主要内容。
1.船舶摇摆船舶摇摆是指在外力作用下船舶产生倾斜，当外力消除后船舶围绕原平衡位置所的往复运动。
其中运动显著而影响严重的是横摇、纵摇和垂荡。
2.砰击由于严重的纵摇和垂荡，船体与风浪之间产生猛烈的局部冲击现象称为砰击。
砰击多发生在船首部。
砰击发生时首柱底端或船底露出水面，然后在极短的时间内以较大的速度落入水中而发生猛烈的撞击。
船舶耐波性总结2讲解
船舶耐波性总结第一章耐波性概述一、海浪的描述、、。
船舶耐波性是船舶在波浪中运动特性的统称，它包括船舶在波浪中所产生的各种摇荡运动以及由这些运动引起的抨击、飞溅、上浪、失速、螺旋桨飞车和波浪弯矩变化等性能，直接影响船舶在风浪作用下维持正常功能的能力。
二、6个自由度的摇荡运动船舶任意时刻的运动可以分解为在Oxyz坐标系内船舶中心G沿三个坐标轴的直线运动及船体绕三个坐标轴的转动。
而这些运动中又有直线运动和往复运动垂荡对船舶航行影响最大，是研究船舶摇荡运动的主要内容。
船舶摇荡是指船舶在风浪作用下产生的摇荡运动，他们的共同特点是在平衡位置附近做周期性的震荡作用。
产生何种摇荡运动形式取决于船首方向与风浪船舶方向之间的夹角，称为遭遇浪向。
三、动力响应船舶耐波性是船舶在风浪中性能的总的反应，它主要包括船舶摇荡、砰击、上浪、失速、螺旋桨飞车。
剧烈的横摇、纵摇和垂荡对船舶产生一系列有害的影响，甚至引起惨重后果，主要表现在以下三个方面：1）、对适居性的影响；2）、对航行使用性的影响；3）、对安全性的影响；船舶在风浪中产生摇荡运动时，船体本身具有角加速度和线加速度，因此属于非定常运动。
第二章海浪与统计分析2-1 海浪概述风浪的三要素：风速、风时、风区长度。
风浪要素定义：表观波长、表观波幅、表观周期。
充分发展海浪条件：应有足够的风时和风区长度。
海浪分类：风浪、涌浪、近岸浪。
风浪的要素表示方法：统计分析方法。
2-2规则波的特性波面可以用简单的函数表达的波浪称为规则波。
A 0=cos kx -t ξξω（）A k ξξω为波面升高，为波幅，为波数，为波浪圆频率。
在深水条件下，波长T c λ、周期和波速之间存在以下关系 ：≈； 2=1.56T λ； c==1.25T λλ； 2=T πω； 2k=g ω 波浪中水质点的振荡，并没有使水质点向前移动，也没用质量传递。
但是水质点具有速度且有升高，因此波浪具有能量。
余弦波单位波表面积的波浪所具有的能量2A 1E=g 2ρξ2-3不规则波理论基础一、不规则波的基本概念 1、确定性关系和统计关系我们所讨论的不规则波引起的船舶摇荡运动等都是属于统计规律范畴之内的。
舰船耐波性试验规范
2 试验目的
测定舰船耐波性能；寻求实船试验与模型试验的相关性；检验理论预报的精度和可 信度；寻求减缓舰船摇荡运动的有效措施。
3 试验仪器、设备
3.1 保证船 试验时必须提供保证船，该船应能在试验所要求的浪级下安全航行并履行其工作职
责。 3.2 测试用仪表
测试用仪表应按国家计量法的规定经过计量检定合格并处在规定的有效周期内，其 量程和精度与试验检测的要求相适应。
试验前后应对所有仪器仪表检查、校验和标症。在试验现场待试期间应保持良好状 态。 3.3 测试仪表的安装
测试仪表应按照测试项目的要求及有关操作规程安装在合适的位置上，并注意防 潮。防止在试验过程中因松动、振动及外界环境等因素影响测试结果的正确性。陀螺仪 应安装在船体重心附近。
7 试验结果的分析和评定
7.1 舰体运动和波高幅值的统计值与其均方根值间的关系 舰体运动和波高幅值的统计值按式(1)计算：
X m = Kf m0 ………（１） 式中：Xm――舰体运动和波高的单幅统计值。对于波高，单位为m；对于纵、横角，单 位为（°）；对于垂向运动加速度，单位为m/s2；
基于三维势流理论的耐波性数值预报
基于三维势流理论的耐波性数值预报祝昊;高霄鹏【摘要】准确预报船体运动响应对于砰击等波浪载荷的计算以及合理结构设计具有重要意义。
船舶在大幅波浪中的运动呈明显的非线性，而现阶段耐波性预报多采用线性切片方法。
三维水动力分析软件 WASIM基于时域势流理论，采用Rankine面元法预报船舶在波浪中的运动响应，并考虑了多种非线性因素。
本文以标模 DTMB5512为对象，采用 WASIM预报其在不同航速下的耐波性，并与基于线性切片理论的计算结果和模型试验结果进行对比。
结果表明：利用 WASIM 计算得到的船体运动响应比其他方法更接近试验值，合理体现了船舶在风浪中的实际耐波性能。
因此，利用 WASIM能够较好地评估船舶在波浪中的非线性耐波特性。
%Accurate seakeeping perediction is of great significance to the calculation of slamming and other ship wave loads. Though the motion response of ship in large amplitude waves presents obvious nonlinearity, the linear strip theory is still adopted for most seakeeping predictions at the present stage. WASIM is a three dimensional hydrodynamic analysis software based on potential flow theory in time domain, it can predict ship motion response in waves with rankine panel method and takes many nonlinear factors into consideration. In this paper software WASIM was used to predict the seakeeping performance of standard model DTMB5512 at different speed, meanwhile the calculation results by WASIM were compared with those by the linear strip method and those by experiment. As a result, it indicated that the data of ship motion response calculated by WASIM are much closer to experimental value than others and reflectedthe actual seakeeping characteristics in waves. Therefore, it can be concluded that WASIM can give a better evaluation of the nonlinear seakeeping characteristics of ship in large amplitude waves.【期刊名称】《舰船科学技术》【年(卷),期】2016(038)007【总页数】5页(P16-20)【关键词】耐波性;数值预报;三维势流理论;非线性;线性切片法【作者】祝昊;高霄鹏【作者单位】海军工程大学舰船工程系，湖北武汉 430033;海军工程大学舰船工程系，湖北武汉 430033【正文语种】中文【中图分类】U661.32在大幅波浪中运动的船舶，由于受到波浪的强烈作用，会出现剧烈的摇荡，同时船舶与波浪间必然会出现猛烈的砰击现象。
基于AQWA的LNG运输船耐波性数值仿真计算
随着 能 源 需 求 的增 加 和人 们 环 保 意识 的增 强 ，天然 气
作 为 一种 清洁能 源越来 越受 到青 睐 。对 于液 化天 然气 的海上
长 距 离 运输 ， Ｌ Ｎ Ｇ 船 无 疑 是最 有 效 、最 经 济 的 选择 。耐 波 性 作 为船舶 安全 性能 的一个 重要 指标 ，在船 舶 的设计 和运 营过
范围不断扩大 ，理论 日趋成熟，尤其是随着高速化大容量计
算机 的 出现 ，三 维流 体动力 的数 值计算 得 到迅速 的发 展 。 目 前 ，用三 维水动 力理 论求解 辐射一 绕 射问题 一般 采用 自由格
１ ． 数值计 算理 论
１ ． １ 基本 条件
林 函数法和简单格林函数法（ 又称Ｒ ａ ｎ ｋ ｉ ｎ ｅ 源方法” 。该方法
泊分 析 ，得 出了满 载压港 和压 载到 港状 态 的六 自由度运 动传
递 函数 。
程 中必须 考 虑 。通过 数 值 计算 可 以得 到 Ｌ Ｎ Ｇ 船 在 波浪 中 的 运 动响应 ，可 以为其 设计 提供参 考 。
耐波 性 的 理论 研 究 经 历 了半 个 多 世 纪 ，并 且 取得 了重 大突破 。从 二维 理论发 展到 三维理 论 ，从 频域 分析 发展 到 时 域 分析 ，从 线性假 设发 展到 非线性 处理 ，船舶 耐波 性 的预报
本 文使 用 基于 三维 势流 理 论 的Ａ Ｑ ＷＡ 软 件 对
１ ６ ５ ０ ０ ｍ Ｌ Ｎ Ｇ 运 输 船 进行 耐波 性 分 析 ，计 算 其 不 同频 率 和 不
同浪 向下 的六 自由度 运动 响应 ，得 出其 传递 函数 ，为 Ｌ Ｎ Ｇ 船
的设计 和航行 提供 一定 的指导 。
耐波性预报的数值方法
关键词：
耐波性；ＡＱＷＡ；ＳＴＦ；三维势流理论
Ｕ６６１．３２
中图分类号：
文献标识码：
Ａ
文章编号：１６７２—７６４９（２０１３）１０—００１６—０５ Ｎｕｍｅｒｉｃａｌ ｓｔｕｄｙ
ＬＩ
ｄｏｉ：１０．３４０４／ｊ．ｉｓｓｎ．１６７２—７６４９．２０１３．１０．００４
ｏｎ
ｔｈｅ ｓｅａｋｅｅｐｉｎｇ ｐｒｅｄｉｃｔｉｏｎ
主值；Ｊｏ为Ｂｅｓｓｅｌ函数的第一项。 假设每个面元表面上的源强度为常数，通过边界 条件求解如下积分方程计算：
ｌ，（￡）＝∑日（∞。）ｆ。。ｅ州，
另Ｙ。。＝ｆ。。Ｈ（０９；），则Ｙｏ。的数学期望为
（１２）
盯乙＝ＩＨ（∞。）Ｉ ２盯弘２，
因此，船舶摇荡运动响应谱密度可表示为
（１３）
击盯盯逊掣ｄＳ。
Ｓ，（∞）＝Ｉ Ｈ（０９）Ｉ ２ｓｆ（∞）。
收稿日期：２０１２—１０—１５；修回日期：２０１２—１１—０２ 作者简介：李宇辰（１９８８一），男，硕士研究生，研究方向为舰艇操纵理论与方法。
万方数据
第１０期
李宇辰，等：耐波性预报的数值方法
・１７・
取得了重大突破。从二维理论发展到三维理论， 从频域分析发展到时域分析，从线性假设发展到 非线性处理，船舶耐波性的预报范围不断扩大， 理论日趋成熟，尤其是随着高速化大容量计算机 的出现，三维流体动力的数值计算得到迅速的发 展。目前，用三维水动力理论求解辐射一绕射问 题一般采用自由格林函数法和简单格林函数法 （又称Ｒａｎｋｉｎｅ源方法口。）。该方法是在物面和自 由面上都分布奇点的一种计算方法，与自由面格 林函数方法相比，该方法在分布奇点计算上较为 简单，并可将自由面非线性和定常势的影响考虑 进去。Ｒａｎｋｉｎｅ源方法最先是在研究开尔文波和兴 波阻力时提出，之后Ｎｏｂｌｅｓｓｅ和Ｔｅｌｓｔｅ‘３－５］在满足 不同边界条件的格林函数的数值问题上，提出了 大量快速的计算方法。目前，基于频域格林函数 的边界元法已成为大型近海结构物设计过程中一 个标准计算工具。本文使用基于三维势流理论的 ＡＱＷＡ软件对方形系数为０．７的系列６０船型进行 耐波性预报，并基于ＳＴＦ法的Ｓｅａｋｅｅｐｅｒ软件计算 结果和该船模型水池试验结果用作对比数据，计 算结果表明，基于ＡＱＷＡ的三维计算方法预报耐 波性精度高，适用性强。
数值预报方法
能相当准确地报出 3 天以内的高空、 地面形势，预报准确率己超过主观预报。 对气象要素（如降水、温度、风、云、 能见度等）的预报，准确率不高，且耗费 太多的计算时间，而且有的项目还难以直 接用数值预报方法作出预报
统计预报的优缺点： 优点: 不仅对气象要素预报， 而且对形 势预报也取得了一定的成效。 缺点: 缺乏物理基础， 特别是对预报量 有时间滞后的预报因子的相关 关系的优良性不容易提高 统计关系的稳定性也不够好。
。
有很多方法用来建立统计预报方程。 以回归法为例： 假定我们要做温度预报，从历史 资料中我们已经找到与温度 (T) 有关的 因子若干个，设为 x1, x2, x3, x4, x5,…..。 可以假设它们满足以下的方程：
T a1 x1 a2 x2 a3 x3 a4 x4 a5 x5
2.资料分析和初始化
对观测资料进行处理形成初值
观测资料与模式不匹配：
来源于每天不同时刻 来自不规则分布的观测站 是标准等压面上的数据， 与预报模
式的垂直分层不同 观 测数据可能有误差和与模式不 协调
资料处理包括以下工作： （1）资料检误 （2）客观分析 （3）资料初始化 （4）四维同化
（3）台风路径业务预报模式： 每天在世界时 00 点和 12 点做两次 48 小时的台风路径预报。 （4）环境紧急响应业务系统： 对事故紧急响应，预测核污染物扩 散和迁移的路径、浓度和沉降量. 及时提供给 WMO。
1996 年我国国家气象中心被世界气象组织 （WMO）确认为亚洲区域的环境紧急响应中心，
§5. 卫星、雷达探测资料的应用
一．气象卫星探测资料的应用 卫星自上而下观测地球大气，所拍 摄的云图直观、形象，加深了预报人员 对天气系统的理解，特别是在资料稀少 的高原、沙漠、海洋上，卫星资料起着 更为重要的作用。
FPSO耐波性研究
FPSO耐波性研究亓俊良;杨贵强;温宝贵【摘要】本节讨论了FPSO的耐波性(seakeeping),包括船体在波浪中产生的各种摇摆运动、拍击、上浪和波浪弯矩的变化.它关系到FPSO在风浪中的工作能力和风暴自存的安全性.历来是设计和使用人员十分关心的问题.【期刊名称】《资源节约与环保》【年(卷),期】2012(000)005【总页数】5页(P156-159,181)【关键词】耐波性【作者】亓俊良;杨贵强;温宝贵【作者单位】中海油能源发展采油服务公司,天津塘沽,300452;中海油能源发展采油服务公司,天津塘沽,300452;中海油能源发展采油服务公司,天津塘沽,300452【正文语种】中文1、目的研究FPSO耐波性的目的在于认识和了解船舶在波浪中产生各种运动的机理，以及它们的危害；探讨防止或减小运动及其危害的措施。
在确定船舶规模时，应保证整个系统（船和系泊）的纵荡和横荡周期至少大于设计波谱峰周期的3倍。
这样做将减小周期为3—25s范围内的波频能量对整个系统的响应幅。
保证纵摇和横摇固有周期大于设计波谱峰周期，对于减小FPSO的波频运动也是很值得的。
然而，在很多情况下难以做到，安装舭龙骨对增大横摇运动阻尼很有效。
耐波性优劣是衡量船舶性能好坏的重要指标之一。
对于在环境条件恶劣海域作业的钻井船和船形FPSO，耐波性是最重要的性能。
2、FPSO耐波性的主要研究内容1）船体六个自由度的运动，其中以横摇、纵摇和垂荡影响最为重要。
剧烈的摇摆运动可能引起一系列不利，甚至严重的后果。
2）拍击，严重的纵摇和垂荡引起船体与波浪之间猛烈的相对运动，即局部冲击。
拍击多发生在船首，严重的拍击使船体振颤，船首局部结构损坏。
3）甲板上浪，船舶在波浪中剧烈摇荡，海水冲上甲板。
上浪可能对人员，设备和结构形成威胁。
船舶摇摆运动主要是由波浪扰动引起的，横摇、纵摇和垂荡影响最大，是主要的研究对象。
由于波浪是随机的，船舶摇摆也是随机的。
3、研究的理论基础1）线性叠加原理195 3年，丹尼斯和皮尔逊将随机理论中的线性叠加原理引进到船舶摇荡研究中来，这样，就把十分复杂的船舶摇荡问题转化为概率论和数理统计理论上来了，使得预报船舶对不规则波的运动响应成为可能。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
因为线性系统具有叠加性，所以系统收到不同 扰动产生的总响应可以用系统受到每个扰动产生的 响应值之和来表示，这样，船体摇荡的随机函数可 表示为
式中：ｌ，＝∞２／ｇ＝ｋｔａｎｈ（ｋｄ）；ｒ＝√（戈一ｆ）２＋（ｙ一叼）２； Ｒ＝√（戈一ｆ）２＋（Ｙ一田）２＋（ｚ—ｆ）２； Ｒ’＝￣／（菇一手）２＋（Ｙ一叼）２＋（ｚ＋２ｄ＋ｆ）２；ｐｖ为积分的
收稿日期：２０１２—１０—１５；修回日期：２０１２—１１—０２ 作者简介：李宇辰（１９８８一），男，硕士研究生，研究方向为舰艇操纵理论与方法。
万方数据
第１０期
李宇辰，等：耐波性预报的数值方法
・１７・
取得了重大突破。从二维理论发展到三维理论， 从频域分析发展到时域分析，从线性假设发展到 非线性处理，船舶耐波性的预报范围不断扩大， 理论日趋成熟，尤其是随着高速化大容量计算机 的出现，三维流体动力的数值计算得到迅速的发 展。目前，用三维水动力理论求解辐射一绕射问 题一般采用自由格林函数法和简单格林函数法 （又称Ｒａｎｋｉｎｅ源方法口。）。该方法是在物面和自 由面上都分布奇点的一种计算方法，与自由面格 林函数方法相比，该方法在分布奇点计算上较为 简单，并可将自由面非线性和定常势的影响考虑 进去。Ｒａｎｋｉｎｅ源方法最先是在研究开尔文波和兴 波阻力时提出，之后Ｎｏｂｌｅｓｓｅ和Ｔｅｌｓｔｅ‘３－５］在满足 不同边界条件的格林函数的数值问题上，提出了 大量快速的计算方法。目前，基于频域格林函数 的边界元法已成为大型近海结构物设计过程中一 个标准计算工具。本文使用基于三维势流理论的 ＡＱＷＡ软件对方形系数为０．７的系列６０船型进行 耐波性预报，并基于ＳＴＦ法的Ｓｅａｋｅｅｐｅｒ软件计算 结果和该船模型水池试验结果用作对比数据，计 算结果表明，基于ＡＱＷＡ的三维计算方法预报耐 波性精度高，适用性强。
万方数据
第１０期
李宇辰，等：耐波性预报的数值方法
・１９・
四边形单元，并且注意软件本身对网格质量的 要求：“。 计算时，船体模型的总重量、重心以及对重心 的惯性矩与实船应当一致。考虑到计算的速度和精 度，计算时划分的总网格数为５０５个，其中三角形 单元为３２个，四边形单元为４７２个，质量点单元为 １个。其具体模型如图１所示。
Ｘｉａｏ・ｌｏｎｇ，ＨＡＮ Ｘｉａｏ—ｇｕａｎｇ
Ｙｕ－ｃｈｅｎ，ＷＵ
Ｘｉ，ＹＵＡＮ
（Ｎａｖｙ
Ａｂｓｔｒａｃｔ：Ｔｈｅ ｔｈｅｏｒｙ．Ｉｔ ｉｓ ｕｓｅｄ
ｉｓ
ｒａｒｅ
Ａｒｍｓ Ｃｏｍｍａｎｄｉｎｇ Ａｃａｄｅｍｙ，Ｇｕａｎｇｚｈｏｕ ５ １
０４３０，Ｃｈｉｎａ）
ｗａｓ
ＡＱＷＡ，ａ
ａ
ｃｏｍｍｅｒｃｉａｌ ｈｙｄｒｏｄｙｎａｍｉｃ ｓｏｆｔｗａｒｅ，ｗｈｉｃｈ
定到定常状态后的速度势，这部分速度势与时间 无关，把它从总速度势中分离出来其余的非定常 部分，记作西，（Ｘ，Ｙ，ｚ）。速度势的分解如下式 所示： 中（戈，Ｙ，ｚ）＝西ｓ（ｘ，Ｙ，ｚ）＋中ｒ（ｘ，Ｙ，ｚ）， 式中：多。为定常势；西，为不定常势。 而不定常势西，又可以分解为人射波势、辐射 势和绕射势３部分，这样下面的表达式成立： ｆ中ｒ（戈，Ｙ，ｚ）＝①Ｊ（戈，Ｙ，ｚ）＋中月（ｘ，Ｙ，ｚ）＋ 函。（戈，Ｙ，。）， （１）
第３５卷
不为：
Ｇ（Ⅵ，锘，毗）＝去＋专＋
２ｐ吐。面顽万（ＩＸ
２１Ｔｉ
Ｆ＋ｖ瓦）ｅ－＂砑ｚ
ｐ＋
（１０）
式中：ｆ。为波幅，将式（１０）代入式（９）可）ｃｏｓｈ弘（ｚ＋ｄ）Ｊｏ（ｐ７）ｄ
ｉｉ‰ｃ。ｓ＾尼（ｚ＋ｄ）ｃｏｓｈｋ（ｆ＋ｄ）厶（ｋｒ）。（５）
关键词：
耐波性；ＡＱＷＡ；ＳＴＦ；三维势流理论
Ｕ６６１．３２
中图分类号：
文献标识码：
Ａ
文章编号：１６７２—７６４９（２０１３）１０—００１６—０５ Ｎｕｍｅｒｉｃａｌ ｓｔｕｄｙ
ＬＩ
ｄｏｉ：１０．３４０４／ｊ．ｉｓｓｎ．１６７２—７６４９．２０１３．１０．００４
ｏｎ
ｔｈｅ ｓｅａｋｅｅｐｉｎｇ ｐｒｅｄｉｃｔｉｏｎ
１
船舶在波浪中的运动理论
基本条件 在波浪中航行时，船舶的运动响应计算广泛采
［孔告圹詈盱ｏ，
．
１．１
（３）
［趴』》＿０＇ｚ一＾，
Ｌ
用三维势流理论‘６１，其基本条件如下： １）船体周围的流体是均匀、无粘性、无旋、 不可压缩的理想流体； ２）入射波为微幅规则波，即把入射波的波幅 和船体的运动作为一阶无穷小，在讨论问题时只保 留到一阶无穷小，略去二阶或更高阶的无穷小，寻 求问题的稳态解，此时可应用叠加原理； ３）考虑绕射和辐射的影响； ４）运动已持续了相当长的时间，已达稳态， 可用频域理论进行分析。 １．２三维流场中速度势的分解 根据线性假设，船舶在波浪中运动的流场速 度势可以按叠加原理分为定常势和不定常势２部 分。定常部分不考虑波浪影响，是船在静水中稳
Ｊ中，＝Ｒｅ（妒ｏｅｌ“）， １
ｏ
（２）
ｌ中一＝Ｒｅ｛∑＂，ｅ山‘）， ”１
【痧Ｄ＝Ｒｅ（妒，ｅ１“）。 式中：中，为入射势；中。为辐射势；少。为绕射势； 入射波势妒。为已知的；∞为波浪的遭遇频率；ｒ／ｉ为 第Ｊ个运动模式的运动复幅值；９ｊ为第歹个运动模式 单位幅值的摇荡运动的速度势。 １．３速度势与流体力的求解 速度势所满足的条件如下式所示： ［￡］：Ｖ ２９，＝０，（，＝１～７）， ［｜ｓ］：（Ｖ妒，一Ｆ）・／７，＝０，（Ｊ＝１～７），
ａ
ＡＱＷＡ
ｗｉｔｈ ｔｈｏｓｅ ｂｙ Ｓｅａｋｅｅｐｅｒ ｗｈｉｃｈ
ｂａｓｅｄ
ｗａｓ
ｏｎ
ｔｈｏｓｅ ｂｙ ｅｘｐｅｒｉｍｅｎｔ．Ａｓ
ｒｅｓｕｌｔ，ｉｔ
ｉｎｄｉｃａｔｅｓ ｔｈａｔ ｔｈｅ ｄａｔａ ｏｆ ｔｈｅ ｍｏｔｉｏｎ ｒｅｓｐｏｎｓｅ ｗｈｉｃｈ
ｃａｎ
ｃａｌｃｕｌａｔｅｄ ｂｙ
ｔｏ
ＡＱＷＡ
ｒｅａｓｏｎａｂｌｌｙ ｒｅｆｌｅｃｔｅｄ ｔｈｅ ｓｅａｋｅｅｐｉｎｇ ｐｅｒｆｏｒｍａｎｃｅ ｏｆ ｓｅｒｉｅｓ ６０，ＳＯ ｔｈａｔ ｉｔ
日ｌ／３＝２．００￣／ｍｏｏ，
，∞
（１５）
（１６）
巩。２
Ｊ０
Ｓｒ（ｗ）ｄｏＪ。
式中：胃。力为船舶各运动幅值的有义值；ｍ。为响应 谱的０阶谱矩。
［一∞２（Ｍ＋蝇（∞））－ｉｏ以（ｔｏ）＋Ｋ］Ｘ（ｔｏ）＝Ｆ（ｗ）。
（７）
２
数值计算方法
下面使用ＡＱＷＡ对方形系数为０．７的系列６０
式中：Ｍ。为结构质量矩阵；Ｍ。为附加质量矩阵； Ｃ为线性阻尼矩阵；Ｋ为系统刚度矩阵；Ｆ（∞）为波 浪力。
（６）
（１４）
式中：Ｉｓ，（∞）为船舶运动的响应密度；ｌ Ｈ（０９）Ｉ ２为
响应幅值算子；Ｓ，（０９）为波浪谱密度。 通过船舶运动响应谱密度可以求得各运动响应 的有义值：
在求得速度势后，船体表面的压力由线性化的 Ｂｅｒｎｏｕｌｉ方程表示，积分后可以得到作用于船体的流 体力。 １．４运动方程的建立 在随船平动的坐标系下，利用微幅、线性化的假 设条件，船舶在波浪中的运动响应Ｘ（∞）通过下面方程 得到：
第３５卷第１０期
２０１３年１０月
舰船科学技术
ＳＨＩＰ ＳＣＩＥＮＣＥ ＡＮＤ ＴＥＣＨＮＯＬＯＧＹ
Ｖ０１．３５，Ｎｏ．１０ Ｏｃｔ，２０１３
耐波性预报的数值方法
李宇辰，吴 唏，袁小龙，韩晓光
（海军兵种指挥学院，广东广州５１０４３０） 摘
要： 基于三维势流理论的水动力商业软件ＡＱＷＡ用于多种浮式结构物在波浪上运动的分析。而在国内
妒（ｘ，Ｙ㈡＝去Ⅱ祁（Ⅵ一亭ｍｆ）邯。（４）
式中：盯为源强；Ｓ为入水结构面，（石，Ｙ，彳）为流 场中域点的坐标；（ｆ，７／，ｇ－）为Ｓ上源点的坐标；Ｇ 为格林函数，并且满足所有边界条件，Ｇ可表
万方数据
・１８・
舰船科学技术 Ｙ（ｔ）为输出函数。 如果将波浪振荡用复数形式表示为 Ｘ（ｔ）＝ｆ。ｅ”‘，
Ｖ９ｊ＿０，ｚ＝一∞，
［Ｒ］：远离物面的自由面上的波向后传播。 式中：／７，为物面上点的广义法向量； ［￡］为控制
方程；［｜ｓ］为物面条件；［Ｆ］为线性化的自由表 面条件；［Ｂ］为底部条件；［Ｒ］为远方的辐射 条件。 辐射势妒ｉ（，＝１～６）和绕射势妒，可以通过三维 源汇分布理论，采用边界元法求解，一般的速度势 可表示为
ｓｔｒｕｃｔｕｒｅｓ
ｂａｓｅｄ
ｏｎ
３－Ｄ ｐｏｔｅｎｔｉａｌ ｆｌｏｗ
ｔｏ
ａｎａｌｙｓｅ
ｗｉｄｅ ｒａｎｇｅ ｏｆ ｆｌｏａｔｉｎｇ
ｍｏｖｉｎｇ ｉｎ ｔｈｅ ｗａｖｅｓ．Ｉｎ
ｄｏｍｅｓｔｉｃ，ｈｏｗｅｖｅｒ，ｉｔ
ｐａｐｅｒ ａｉｍｓ ｔｏ
ｔｏ
ｓｅｅ
ｐｅｏｐｌｅ ｃａｌｃｕｌａｔｉｎｇ ｔｈｅ ｍｏｔｉｏｎ ｒｅｓｐｏｎｓｅ ｏｆ ｓｈｉｐ ｂｙ
大量事实证明，一艘在静水中航行性能优良的 船舶，在波浪中的航行性能并不一定优良¨ｊ。由于 船舶实际上是经常航行在波涛汹涌的海面上，因此， 研究船舶在波浪上的运动越来越引起人们重视。总 体而言，船舶耐波性的指标包括横摇、纵荡、垂荡、 首摇、横荡、纵荡，以及特殊部位的位移、运动速 度、运动加速度、甲板上浪、首底碰击、螺旋桨出 水次数等运动响应的有义幅值的指标。预报这些指
ｃａｎ
ＡＱＷＡ，ｔｈｅｒｅｆｏｒｅ，ｔｈｉｓ
ｄｉｓｃｕｓｓ
ｗｈｅｔｈｅｒ
ＡＱＷＡ
ｂｅ ｕｓｅｄ
ｔｏ
ｐｒｅｄｉｃｔ ｔｈｅ ｓｅａｋｅｅｐｉｎｇ ｐｅｒｆｏｒｍａｎｃｅ ｏｆ ｓｅｒｉｅｓ
６０（Ｃ６＝０．７）．Ａｔ
ｗａｓ
ｔｈｅ ｓａｍｅ ＳＴＦ ａｎｄ
ｔｉｍｅ，ｃｏｍｐａｒｅｄ
ｔｈｅ
ｃａｌｃｕｌａｔｉｎｇ ｒｅｓｕｌｔｓ ｂｙ
ｂｅ ａｐｐｌｉｅｄ
ｗａｖｅ ｌｏａｄ
ｐｒｅｄｉｃｔｉｏｎ． Ｋｅｙ ｗｏｒｄｓ：