驱逐舰耐波性与人因中晕动症相关性分析

熊虎，卢晓平

海军工程大学舰船工程系，湖北武汉430033

驱逐舰耐波性与人因中晕动症相关性分析

熊虎，卢晓平

海军工程大学舰船工程系，湖北武汉430033

舰船耐波性不仅是衡量舰船总体性能的重要指标，也对舰船上人员晕动症的发生有直接影响，因此，良好的耐波性设计是保障舰员舒适度和绩效的基础。采用切片法理论编写耐波性计算程序，对DTMB 5415标模进行耐波性计算，得到在迎浪中不同航速下船艏、重心处的垂向加速度响应，分析DTMB 5415标模对应的驱逐舰升沉运动时垂向加速度对人因中晕动症造成的影响，给出驱逐舰耐波性与人因中晕动症的数值关系。结果表明：船艏和重心处垂向加速度理论计算结果正确，准确度满足晕动症预报的要求；计算所得晕动症发生率规律合理，可作为驱逐舰船型晕动症发生率分析和预报，以及人因改善的资料。

耐波性；切片法；人因工程学；晕动症；驱逐舰

0 引 言

随着舰船自动化程度的不断提升，舰员编制数量持续减少，人因工程学的运用成为驱逐舰设计过程中不可分割的一部分。其重要性体现在2个方面：一是人因工程学的合理运用为减少舰员数量，提高舰船自动化程度创造了条件；二是人因工程学的充分应用能直接提高舰员绩效，从而提升舰船的安全性及综合作战能力。据统计，约80%的灾难、事故、过失和失败归因于人的因素（即“人因”）［1］。

在海军舰船的设计和使用中，与人因联系最为密切的是舰船耐波性。人因的研究与舰船耐波性研究的关系大致可以表述为：耐波性研究是人因研究的重要基础和一个方面，而人因工程研究和人因问题的解决是耐波性研究的深入和主要应用目标之一。要深入和有效地开展海军舰船的设计和使用中人因的研究，必须具备坚实的耐波性研究基础。

当前耐波性的研究方法主要有船模试验法、CFD方法和理论计算方法［2］。船模试验方法受到船模尺寸、流场扰动、测量精度、人为因素的限制，并且试验方法还需要耗费大量的人力、物力与财力，通过试验方法获取耐波性结果并不总是令人满意。近年来，CFD方法发展迅速，但由于船舶耐波性问题属于变化物面、带自由表面的非线性流动问题，船体摇摆运动、入射波运动、船舶航行和摇摆兴波等各种运动及其受力的叠加耦合极其复杂，使得CFD在解决耐波性问题时精度不够、耗时巨大，难以获取理想的结果。势流理论计算方法主要有切片法、二维半方法和三维方法［3］。三维方法由于计算复杂且精度较二维方法并没有显著的提升，没有得到广泛应用；二维半方法是一种适于计算舰船高速运动时耐波性的方法，最近几年得到了广泛应用，计算结果较为准确；切片法是一种二维方法，经过半个世纪的发展，切片理论已成为耐波性的基本理论之一，常规的切片理论适于预报细长型舰船在斜浪中六个自由度的运动［4］。无论是模型试验方法、CFD方法，还是切片法、二维半方法和三维方法，国内耐波性的研究尚未发展到考量舰船设计和使用中人因问题的新层次，极少见到舰船耐波性与舰船人因指标相关性的研究。在国外，这方面的研究早见端倪，至今已初显规模。

本文拟在耐波性切片法计算的基础上，针对大型驱逐舰标模DTMB 5415开展耐波性与人因中晕动症关系的研究，得出晕动症与垂向加速度的相关性，以及改进船型和耐波性设计对降低晕动症发生率的重要性。

1 人因中的晕动症

1.1 晕动症的概述

晕动症是由最近存储线索与实际垂向线索之间的感官冲突或不匹配造成的。按照惯例，在海洋环境中，认为晕动症最主要的组成因素是船舶垂向运动（垂荡、纵摇和横摇的垂直分量），因为船舶其他的线性运动和角运动对晕动症的影响通常较小［5］。

1.2 晕动症发生率经验公式

2006年2～3月，国外开展了一项研究，127 m长的三体渡轮Benchijigua Express号（图1）在加那利群岛（Canary Islands）之间的公海进行了多次2 h运输作业。该船应用运动传感器和船用雷达仪器提供海况信息。近2 000名乘客填写了详细的调查表。据报道，由于波高和航向的原因，60%～90%的乘客至少有一些晕动症的症状。笔者所在课题组曾在某预研项目中对高速三体船的适航性进行过系统、深入的研究，研究结果表明，等排水量的单体船相对于三体船的适航性、横摇稳性都会差一些。因此，在相同的航态和排水量条件下，常规单体驱逐舰上舰员出现晕动病症状的人数比例不会低于国外对127 m三体船研究所得的结果。在国外127 m三体船的研究中，作为晕动症评估问卷（MSAQ）的项目，胃肠道症状的反应程度与被测试人员的航海次数成反比；在该项试验中，被试人员有着不同的航海经历，有些具有多次航海经验（多于9班次），有些具有一定的航海经验（2～8班次），也有一些缺乏航海经验的（0～1班次），不同航海频次乘员构成的被试组正可以用于分析得出晕动症反应程度与航海次数的相关性。另外，在晕动症评估问卷中，黑暗的夜晚乘船旅行者比早晨（可以见到地平线时）及白天乘船旅行者得分高［6］，即患晕动症的概率大。

图1 澳大利亚建造的三体船Fig.1 Australian-built trimaran Benchijigua Express

通过研究志愿受试者发现，大多数晕动症发生在船舶运动频率为0.125～0.25 Hz时。晕动症的发生率随加速度和遭受船舶运动时间的增加而增加。通常用于衡量晕动症的指标是晕动症发生率（Motion Sickness Incidence，MSI），这也是受试人员遭受船舶运动后的呕吐百分比［6］。

受试者坐在一个封闭、明亮的空调实验室中，经历2 h的单向正弦垂直加速度。选择2 h的时限，是因为晕动症发生率在持续一段时间后（经过一段时日，人们通常会适应船舶运动）不再增加。注意，这里不考虑纵摇和横摇。如果增加升沉分量可以预期到，则由非纯升沉加速度所导致的晕动症发生率超过预估值。利用国外该类项目实验测试数据和结果，McCauley等［7］得出了预测晕动症发生率的估算公式：

式中：MSI为晕动症发生率，按舰员遭受2 h船舶运动后的呕吐概率定义；Φ为高斯分布概率密度的累积值，即高斯概率密度函数的不定积分，其表达式为

式中，z为标准正态变量。

式（1）中，Na，Nt′为与船体在波浪中升沉运动的垂向加速度a、频率f和时间历程t相关的物理量，其估算公式为：

式中：a为垂直加速度的有义值，g；t为测试对象遭受的时间历程，min；f为船体在波浪中的升沉运动频率，Hz。当为正值时，Φ通过查阅数学手册确定，或由式（2）计算得到；当Na和为负值时，Φ由以下关系计算得到：

Colwell认为，晕动症发生率简化公式的计算结果与McCauley所得的数据相符，即平均误差不超过1%、标准偏差不超过6%［8］。

2 规则波中的标模晕动症发生率计算

2.1 规则波中纵摇与垂荡运动

由于舰船晕动症最主要的原因是船舶的垂向运动，通常迎浪时的垂向运动最为剧烈，故可以认为，若舰船在迎浪状态满足垂向运动数值范围，则其他浪向总能满足要求。一般在艏楼常有较大的加速度，因此，根据STF方法，通过编程计算标模DTMB 5415在迎浪中不同航速下艏部和重心处的垂向绝对加速度响应［9］。假设重心在自由液面处。

建立固连于船的坐标系o-xyz。其原点o在船重心G上，xGy平面平行于静止时的载重水线面，Gx轴正向指向船艏，xGz平面在船的纵中剖面上，Gz垂直向上，坐标系与船一起移动和振荡。垂荡和纵摇的耦合方程可表示为

式中：m为排水质量；z为垂荡运动值；I为纵摇惯性矩；θ为纵摇角；A11为垂向附加质量；A22为纵摇附加惯性矩；A12为垂荡与纵摇耦合系数；A21为纵摇与垂荡耦合系数；B11为垂荡阻尼力系数；B22为纵摇阻尼力系数；B12为垂荡与纵摇阻尼耦合力系数；B21为纵摇与垂荡阻尼耦合力系数；C11为垂荡复原力系数；C22为纵摇复原力系数；C12为垂荡与纵摇复原力耦合系数；C21为纵摇与垂荡复原力耦合系数；Fc，Fs分别为扰动力系数的实部和虚部；Mc，Ms分别为力矩系数的实部和虚部；ωe为遭遇频率；t为遭受船舶运动的时间。

由式（6）可得：

在上述方法的计算程序中，波幅取ζa=1，计算za和θa，解出舰模升沉、纵摇的频率响应，进而得出切片的绝对垂向加速度

舰模的遭遇频率为

式中：ω为波浪圆频率；g为重力加速度；v为航速，m/s；β为浪向角。

航行状态分别取傅汝德数Fr=0.1，0.2，0.3，0.4，0.5，0.6，0.7和0.8，浪向角β=180°（迎浪），计算标模DTMB 5415船艏和重心处的垂向加速度，所得值分别如图2和图3所示。图中，H13为有义波高。

图2 船艏垂向加速度Fig.2 Vertical acceleration of DTMB 5415 at the bow

图3 重心处垂向加速度Fig.3 Vertical acceleration of DTMB 5415 at the center of gravity

国内对各类舰船垂向加速度值的理论计算、试验研究、数值模拟研究已很广泛。为验证利用以上方法计算垂向加速度值的准确性和有效性，将DTMB 5415艏部与文献［9］中给出的驱逐舰艏部垂向加速度值进行了对比，如图4所示。

图4 船艏垂向加速度对比Fig.4 The contrast of vertical acceleration at the bow

图4中对DTMB 5415在Fr=0.2，0.3，0.4和0.5与文献［9］所给驱逐舰Fr=0.263，0.351，0.409和0.468的艏部垂向加速度进行了对比。文献［9］中驱逐舰船型与DTMB 5415船型同属方艉驱逐舰船型，且船型参数与傅汝德数接近，其垂向加速度特性大致相当，具有可比性。由图4可知，两驱逐舰的垂向加速度曲线趋势较为一致，数值上大致吻合，其吻合程度已足以说明本文垂向加速度计算方法的正确性和有效性。当然，由于DTMB 5415为国外典型驱逐舰船型，文献［9］中所给驱逐舰为国内早期船型，二者的垂向加速度值存在偏差也在情理和预期之中。对于国外典型驱逐舰DTMB5415的耐波性模型试验和更详细的理论计算方法，作者正在研究之中，这些工作可进一步验证本文计算方法及结果的正确性和有效性。

2.2 晕动症与垂向加速度关系的计算与分析

根据1.2节给出的晕动症发生率计算公式及2.1节中的相关计算结果，可计算出所求驱逐舰标模DTMB 5415对应的驱逐舰在各航态下相应运动频率的晕动症发生率，这里给出一个算例。

算例：根据驱逐舰标模DTMB 5415，计算对应的驱逐舰在典型有义波高典型航速（巡航速度附近）Fr=0.28时，船艏处的晕动症发生率。

当Fr=0.28时，由图2可知船艏处的垂向加速度值为0.282g。上文已确定大多数晕动症发生在船舶运动频率为0.125～0.25 Hz之间，这里取运动频率0.15 Hz进行计算，步骤如下：

1）将数据代入式（3），可得Na1=0.673；

2）将其代入式（4），人员遭受船舶运动的时间取2 h，可得Nt1′=1.995；

3）查阅标准正态分布表，可得

Φ（Na1）=0.56，Φ（Nt1′）=0.91

4）将以上数据代入式（1），得到MSI1= 50.96，即晕动症发生率为50.96%。

重复上述步骤，计算得到各有义波高下相应运动频率和航速的晕动症发生率。标模DTMB 5415对应的驱逐舰晕动症发生率如表1所示，相应的三维视图如图5所示。

由图5可知，在大多数航态下，DTMB 5415对应的驱逐舰晕动症发生率均较低。由于我国驱逐舰排水量较美国“阿利·伯克”级驱逐舰小，在波浪中的摇摆运动（含垂向运动加速度）会更大，初步预期国内驱逐舰晕动症发生率会比上述针对美国“阿利·伯克”级驱逐舰所得结果更高。国内驱逐舰，从舰船耐波性的角度优化船型或采取其他改善耐波性的设计和运营措施，降低驱逐舰的垂向加速度值，从而减弱晕动症发生的主要因素，减少晕动症对驱逐舰整体性能的影响十分必要，对进一步提高国内驱逐舰的作战性能和人因指标均有重要意义。

表1 DTMB 5415对应的驱逐舰晕动症发生率Table 1 DTMB 5415 corresponding destroyer motion sickness incidence

图5 DTMB 5415对应的驱逐舰晕动症发生率Fig.5 DTMB 5415 corresponding destroyer motion sickness incidence

2.3 结果验证

这里结合Piscopo等［10］给出的计算方法校验本文中晕动症发生率计算结果的准确性。

式中：m4为船舶垂向运动波谱的第4阶谱矩［11］；μMSI为系数，由劳氏船级社给出的公式计算得出。

式中，m2为船舶垂向运动波谱的第2阶谱矩［11］，其中根据有义波高H13、特征周期T、遭遇频率ωe得出。

n阶谱矩：

用上述验证方法，同样计算驱逐舰标模对应的驱逐舰在有义波高典型航速Fr=0.28时，船艏处晕动症的发生率为54.19%，比文中计算值50.96%要大3.23%，误差在合理范围之内。使用该方法，计算各航态下船艏和重心处的晕动症发病率，所得晕动症发生率如表2所示，相应的三维视图如图6所示。将采用验证方法所得表2数据与文中方法计算所得表1数据进行对比，误差均在5%左右，说明本文所使用的计算方法准确性较好。

表2 验证方法计算所得晕动症发生率Table 2 Motion sickness incidence calculated by the method in Reference[10]

图6 验证方法计算所得晕动症发生率Fig.6 Motion sickness incidence calculated by the method in Reference[10]

3 结 语

随着海军舰船的不断发展，舰船设计将不可避免地与诸如人因工程等新理念相结合，需从多学科、多层次考量提升驱逐舰的耐波性、操纵性、波浪中稳性以及舰员的绩效、居住舒适度，降低人的因素对舰船战斗性能的不利影响。

本文计算了标模DTMB 5415各航态下迎浪运动时艏部和重心处的垂向加速度值，并进一步计算了不同垂向加速度值时舰员发生晕动症的概率。当然，晕动症是一个受人的自身精神状态、光照、声音、气味、船体运动状态等诸多因素共同作用影响的病症，此处只是将晕动症的主要影响因素（垂向加速度值）纳入考虑范围，可以大致预测出驱逐舰耐波性与舰员晕动症发生率之间的数值关系。未来可以经过实船试验研究，验证理论计算的准确性，完善理论计算方法，并进一步总结舰船耐波性与晕动症的数值关系，以作为舰船设计过程中耐波性和人因工程设计的一个参考。

［1］ 刘丙泉.船舶人因事故预警管理研究［D］.哈尔滨：哈尔滨工程大学，2009. LIU B Q.Research on forewarning management for accidents due to human factors of ships［D］.Harbin：Harbin Engineering University，2009（in Chinese）.

［2］ 倪崇本.基于CFD的船舶阻力性能综合研究［D］.上海：上海交通大学，2011. NIC B.A comprehensive investigation of ship resistance prediction based on CFD theory［D］. Shanghai：Shanghai Jiao Tong University，2011（in Chinese）.

［3］ 马山.基于二维半理论的垂向船舶运动和波浪载荷预报［D］.哈尔滨：哈尔滨工程大学，2002. MA S.Prediction of vertical ship motions and wave loads based on 2-1/2D theory［D］.Harbin：Harbin Engineering University，2002（in Chinese）.

［4］ 戴仰山，沈进威，宋竞正.船舶波浪载荷［M］.北京：国防工业出版社，2007：39-64. DAI Y S，SHEN J W，SONG J Z.Ship wave loads［M］.Beijing：National Defense Industry Press，2007：39-64（in Chinese）.

［5］ LUCOT J.Pharmacology of motion sickness［J］.Journal of Vestibular Research，1998，8（1）：61-66.

［6］ ROSS J M.Human factors for naval marine vehicle design and operation［M］. Farnham： Ashgate Publishing，Ltd.，2009.

［7］ MCCAULEY M E，ROYAL J W，WYLIE C D，et al. Motion sickness incidence：exploratory studies of habituation，pitch and roll，and the refinement of a mathematical model［R］.［S.l.］：Canyon Research Group Inc Goleta Ca Human Factors Research Div，1976.

［8］ GUIGNARD J C，MCCAULEY M E.Motion sickness incidence induced by complex periodic waveforms［J］. Aviation，Space，and Environmental Medicine，1982，53（6）：554-563.

［9］ 彭英声.舰船耐波性基础［M］.北京：国防工业出版社，1989. PENG Y S.Ship seakeeping basis［M］.Beijing：National Defense Industry Press，1989（in Chinese）.

［10］ PISCOPO V，SCAMARDELLA A.The overall motion sickness incidence applied to catamarans［J］. International Journal of Naval Architecture and Ocean Engineering，2015，7（4）：655-669.

［11］ 李积德.船舶耐波性［M］.哈尔滨：哈尔滨工程大学出版社，2007. LIJD.Ship seakeeping［M］.Harbin：Harbin Engineering University Press，2007（in Chinese）.

Analysis of correlation between destroyer seakeeping and motion sickness in human factors

XIONG Hu，LU Xiaoping
Department of Naval Architecture Engineering，Naval University of Engineering，Wuhan 430033，China

Seakeeping is not only an important indicator of the overall performance of a ship，it also has a direct impact on the occurrence of motion sickness among ship personnel；as such,a good seakeeping design is the basis for guaranteeing the crew's comfort and work performance.In this paper,a seakeeping calculation program using the theory of the strip method is written to calculate the seakeeping of the DTMB 5415 standard model,and the vertical acceleration response at the bow and the center of gravity in head waves and under different speeds are obtained.Next,the effects on motion sickness of the vertical accelerations of DTMB 5415 corresponding to the destroyer's heaving are analyzed,and the numerical relationship of motion sickness with the destroyer's seakeeping performance is obtained.The results show that the calculated vertical accelerations at the bow and center of gravity are correct,the accuracy of which satisfies the requirements of motion sickness prediction;and the calculated regularities of the motion sickness are reasonable,and can be applied to analyze and predict the incidence of motion sickness on a ship,as well as to improve human factors in the hull form design of destroyers.

seakeeping；strip method；human factors engineering；motion sickness；destroyer

U661.32

A

10.3969/j.issn.1673-3185.2017.01.006

2016-04-25

2016-12-28 15:57

海军工程大学社会科学基金资助项目（HGDSK2015E06）

熊虎，男，1992年生，硕士生。研究方向：舰船水动力性能研究。E-mail：398820941@qq.com卢晓平（通信作者），男，1957年生，博士，教授。研究方向：舰船水动力性能研究。E-mail：luxiaoping100@163.com

http://www.cnki.net/kcms/detail/42.1755.TJ.20161228.1557.030.html期刊网址：www.ship-research.com

熊虎，卢晓平.驱逐舰耐波性与人因中晕动症相关性分析［J］.中国舰船研究，2017，12（1）：32-37. XIONG H，LU X P.Analysis of correlation between destroyer seakeeping and motion sickness in human factors［J］.Chi⁃nese Journal of Ship Research，2017，12（1）：32-37.