舰船波浪载荷响应特性及砰击载荷影响因素分析.pdf

第 20 卷 第 9 期 装 备 环 境 工 程 2023 年 9 月 EQUIPMENT ENVIRONMENTAL ENGINEERING 1 收稿日期：2023-06-01；修订日期：2023-06-30 Received：2023-06-01；Revised：2023-06-30 基金项目：国家重点研发计划（2022YFB3306200）Fund：National Key R&D Program of China(2022YFB3306200)引文格式：耿彦超,杨骏,王海洋,等.舰船波浪载荷响应特性及砰击载荷影响因素分析J.装备环境工程,2023,20(9):1-11.GENG Yan-chao,YANG Jun,WANG Hai-yang,et al.Response Characteristics of Wave-induced Loads and Affecting Factor of Slamming Loads for ShipJ.Equipment Environmental Engineering,2023,20(9):1-11.舰船波浪载荷响应特性及砰击载荷影响 因素分析 耿彦超1,2，杨骏1,2，王海洋3，卞鑫1,2，汪雪良1,2（1.中国船舶科学研究中心，江苏 无锡 214082；2.深海技术科学太湖实验室，江苏 无锡 214082；3.中国船舶及海洋工程设计研究院，上海 200011）摘要：目的目的 获取波浪载荷特征，提供准确的设计载荷输入，为舰船结构设计服务。方法方法 采用波浪载荷水池模型试验和理论预报对比的方法，给出波浪载荷传递函数响应曲线，验证理论预报结果的正确性。基于长期预报结果，分析波浪载荷沿船长的分布规律，通过改变超越概率和舰船装载状态，研究二者对波浪载荷的影响。基于模型试验结果，阐释航速引起的砰击弯矩的变化规律。之后，通过改变舷侧外伸平台的宽度、角度、高度等参数，分析参数变化对砰击载荷的影响，揭示波浪载荷的变化规律。结果结果 模型试验结果同理论预报结果吻合较好。垂向剪力、扭矩长期预报结果沿船长呈现双峰值形状，垂向弯矩和水平弯矩长期预报结果沿船长呈现抛物线形状，波浪载荷长期预报结果随超越概率呈大致线性变化关系。大波高、高航速试验工况下，船体弯矩呈现特别明显的非线性特征，舷侧外伸平台宽度越大、角度越小、高度越低，砰击载荷就越严重。结论结论 波高、航速等对砰击弯矩的影响显著，舰船在使用过程中，若要减少砰击的发生，可考虑规避大浪和主动减速。在舰船结构设计阶段，若要减小舷侧外伸平台产生的砰击载荷分量，应从平台的大小、位置、角度等参数的设计确定入手。关键词：波浪载荷；长期预报；砰击载荷；航速；外伸平台 中图分类号：U661.4 文献标识码：A 文章编号：1672-9242(2023)09-0001-11 DOI：10.7643/issn.1672-9242.2023.09.001 Response Characteristics of Wave-induced Loads and Affecting Factor of Slamming Loads for Ship GENG Yan-chao1,2,YANG Jun1,2,WANG Hai-yang3,BIAN Xin1,2,WANG Xue-liang1,2(1.China Ship Scientific Research Center,Jiangsu Wuxi 214082,China;2.Taihu Laboratory of Deepsea Technological Science,Jiangsu Wuxi 214082,China;3.Marine Design&Research Institute of China,Shanghai 200011,China)ABSTRACT:The response characteristics of ship wave-induced loads are complex.The work aims to obtain the wave-induced load characteristics and provide accurate design load,so as to serve for ship structural design.The method of comparing basin model experiments with theoretical predictions was adopted to provide the response curve of wave-induced transfer function,which verified the correctness of theoretical prediction results.Based on long-term prediction results,the distribution pattern of 专刊舰船结构安全性 2 装 备 环 境 工 程 2023 年 9 月 wave-induced loads along the ship was analyzed,and the impact of exceeding probability and ship loading status on wave-induced loads was studied.By changing the exceeding probability and the loading state of the ship,the impact of the both on wave-induced loads was studied.Based on the model test results,the law of the change in slamming bending moment caused by ship speed was obtained.The effect of parameter change on slamming load was analyzed by changing the width,angle,height and other parameters of the side extended flaring platform,to reveal the change law of wave-induced load.The model test results were in good agreement with the theoretical prediction results.The long-term prediction results of vertical shear and torque showed a bimodal M shape along the ship length,while the long-term prediction results of vertical and horizontal bend-ing moments showed a parabolic shape along the ship length.The long-term prediction results of wave-induced loads showed a roughly linear relationship with the exceeding probability.Under large wave height and high speed test conditions,the bending moment of the ship exhibited particularly significant nonlinear characteristics.The wider the width,the smaller the angle,and the lower the height of the side extended flaring platform,the more severe the slamming load.Considering that wave height,speed,and other factors have a significant impact on the slamming moment,in order to reduce the occurrence of slamming dur-ing the use of ships,it is considered to avoid large wave height and actively slow down.In addition,in the design phase,in order to reduce the slamming load component generated by the side extended flaring platform,the size,position and angle of the flar-ing platform should be determined from the design.KEY WORDS:wave-induced loads;long-term;slamming loads;speed;flaring platform 舰船在波浪中的波浪载荷响应特征比较复杂，根据船体设计第一原则，应首先确定作用在船体上的环境载荷，掌握了船体环境外载荷的特点，才能为船体结构设计提供准确输入。舰船在恶劣海况下航行时，会发生砰击现象，对大型舰船而言，船长相对较大，高强度钢大量应用，船体逐渐变软，相对刚度下降，遭遇波浪发生的砰击会引起船体颤振，显著增加船体的砰击弯矩，因此应开展船体砰击颤振研究，掌握其砰击载荷特性，为船体结构设计服务。一般而言，波浪载荷的预报值可分为线性结果和非线性结果。线性结果一般是指频域下的线性长期统计值，是基于响应的载荷传递函数，结合相应的波浪散布图，在设计概率下的长期统计结果；非线性结果一般是指包含砰击等非线性成分的预报结果。近年来，国内外一些学者采用不同的方法手段，研究了船体的波浪载荷特性，在波浪载荷理论预报发展方面，取得了一些进展。波浪载荷理论预报，从船体周围流场处理方法来讲，分为二维切片理论、二维半理论与三维理论。二维切片理论是通过将船体波浪载荷三维流动问题转化为二维剖面流动问题来解决的，是目前较为成熟的预报手段。之后，学者们提出了应用二维控制方程加三维自由面条件来求解船体响应的方法，称为二维半理论，该种理论最早由Faltinsen 等1用 于 对 船 舶 耐 波 性 能 的 预 报。Hermundstad 等2应用二维半理论对船舶在波浪中的水弹性响应进行了求解，其数值求解方法和 Faltinsen等的不同之处在于没有考虑定常速度势对非定常势的影响，以避免求解定常兴波势带来的数值问题。段文洋等3基于二维时域格林函数的数值方法，应用二维半理论预报了船体的垂向运动响应，结果同模型试验符合得很好。随着科技的进步，三维方法逐渐发展起来，逐渐成为波浪载荷预报的主流。杨鹏等4-6基于混合边界元法，结合模态叠加法建立了基于内外场匹配技术的三维时域水弹性方法，通过边界积分法和在流体内域自由面施加刚性边界，消除了不规则频率问题，并且建立了有限水深的三维时域格林函数水弹性计算方法。Bingham 等7采用频域三维脉动源、三维移动脉动源、三维脉动源与水弹性理论相结合的 3种方法，对船舶在波浪中的运动及波浪载荷进行了研究，对比了 3 种方法的数值结果，发现在低航速下，基于三维移动源的载荷预报结果与三维脉动源的预报结果差异不大，而在高航速下差异较为明显，采用水弹性方法的预报值和刚体方法的预报值在部分工况中存在差异。Kim 等8-9采用 Newton 法在时域中求解了零航速船体与流场的流固耦合问题，对船体进行了三维有限元建模，把船作为弹性体，可以准确地模拟船体结构变形。李飞等10采用非线性波浪载荷方法对超大型集装箱船在不同船体梁刚度下的垂向波浪载荷响应进行了理论计算，分析了船体梁刚度对波浪载荷的影响。计算结果显示，船体梁刚度的下降会导致作用在船体梁上的波浪载荷显著增大。以上为波浪载荷的理论发展情况。在波浪载荷模型试验研究方面，也有一些新的研究成果。耿彦超等11对一艏江海直达船的波浪载荷进行了水池模型试验研究，研究分析了该船在波浪中的水弹性响应。Jiao 等12-13依据水池模型试验和大尺度海上模型试验，研究了普通舰船和中高速复合型船在二维与三维不规则波中的波浪载荷响应。Si 等14提出通过 U 形梁横剖面的面内剪流来获得船体承受扭矩的测量方法。Grammatikopoulos 等15加工制作了一艘驳船的弹性船模，并在水池中开展了相关试验研究，给出了线性与非线性等因素的影响结果。耿彦超等16第 20 卷 第 9 期 耿彦超，等：舰船波浪载荷响应特性及砰击载荷影响因素分析 3 依据模型试验结果分析了载荷特性，认为若把船当作刚体，船体弯矩响应除了有低频遭遇频率成分外，还有一定范围的倍频成分；若把船认为是弹性体，则弯矩响应除了上述成分外，还有一阶船体梁垂向总振动频率成分，且一阶船体梁垂向总振动频率也会放大倍频成分。另外，其他一些学者也对波浪载荷水池模型试验开展了分析研究。司海龙等17对波浪载荷模型试验的现状进行了阐述。焦甲龙等18开展了大尺度的海上波浪载荷模型试验。也有其他学者19-23开展了波浪载荷试验研究，对船体砰击弯矩和波激振动引起的弯矩进行了分析，模型试验方法是相似的，给出了砰击弯矩和波激振动的变化规律。上述研究工作主要是对于波浪载荷理论预报方法和水池模型试验方法研究的介绍，而分析舰船波浪载荷的变化特性及砰击载荷影响因素的论文较少，尤其是舷侧外伸平台变化引起的砰击载荷变化的研究工作更为少见。基于上述研究工作现状，为有针对性地开展波浪载荷特性分析，本文拟采用理论预报与水池模型试验相结合的手段，获得波浪载荷传递函数响应结果，并开展波浪载荷长期预报，且定量地给出了船体波浪载荷的设计值，分析其沿船长的分布规律，研究超越概率和舰船装载状态变化对波浪载荷的影响，分析航速引起的砰击弯矩变化，分析舷侧外伸平台参数变化对砰击载荷的影响，揭示舰船波浪载荷变化规律，给出舰船在使用过程中和结构设计阶段应重点关注的问题。技术路线如图 1 所示。图 1 技术路线 Fig.1 Technology roadmap 1 舰船波浪载荷的长期预报响应 1.1 模型试验 波浪载荷试验模型与实船应满足几何相似、水动力相似、运动相似、结构刚度相似和质量相似等24。基于这 5 个方面的相似理论，用于指导模型的设计、制作、惯量调整。基于试验状态工况，开展波浪载荷的水池模型试验工作。模型的主尺度见表 1。表 1 舰船模型主尺度 Tab.1 Principal dimensions of ship model 总长/m水线长/m型宽/m 型深/m 吃水/m排水量/t5.12 4.62 0.69 0.55 0.17 0.42 模型试验采取自航方式。试验模型将采用玻璃钢材料制作而成。模型沿着主船体的长度方向分为 5段，在第 3、7、10、15 站分段，用测量钢制圆管梁相互连接成一整体，试验中模型刚度仅取决于测量梁刚度，船壳模型只起提供浮力及传递流体动力作用。波浪载荷水池试验模型如图 2 所示。图 2 某大型船水池模型试验 Fig.2 Basin model test of the large ship 1.2 理论预报方法 波浪载荷预报，常分为短期预报和长期预报 2类。短期预报的时间范围为 30 min 到数小时，在此时间内，船的装载状态、航速、航向角以及海情都可以认为是固定不变的。长期预报的时间范围是数年或整个寿命期，在此时间内，上述因素都会改变的，长期预报是由许多短期预报组成的24。短期海浪可视为均值为 0 的平稳正态随机过程，当船体对波浪的响应是线性系统时，由随机过程理论可知，此时的波浪载荷（输出）亦是均值为 0 的平稳正态随机过程，波浪 S与波浪载荷 SW之间的关系见式（1）。SW(,H1/3,TZ,+)=H2(,v,+)S(,H1/3,TZ,)(1)式中：S(,H1/3,TZ,)为波浪谱密度；SW(,H1/3,TZ,+)为波浪载荷谱密度；H(,v,+)为系统传递函数，单位规则波幅下的载荷响应幅值；为波浪圆频率；v 为航速；为航向角；H1/3为有义波高；TZ为波浪特征周期。由概率论可知，对于一个零均值的平稳正态随机过程，在窄谱假设下，其幅值服从 Rayleigh 分布，它的概率密度为：202expxxfxEE(2)1/3Z2W1/3Z022,(,)d dEm HTS HT (3)4 装 备 环 境 工 程 2023 年 9 月 由此得到波浪载荷的各种特征值。由短期预报可知，在海况(H1/3,TZ)中，船舶以航向角、航速 v 运行时，其波浪载荷幅值 X 小于某个可能值 x 的概率为：201 expxfxE (4)如果认为各个短期分布是彼此独立的，那么长期概率分布将是各短期概率分布的加权组合，即波浪载荷幅值 X 大于某一定值 x 的超越概率为：21/3Z,expijkijkQ xp Xxxp HTppvE (5)一旦船舶运行海域和概率水平确定后，即可得到对应的波浪载荷特征最大值 Xmax。此值表示船舶在波浪遭遇次数为 n 的整个使用期内，最可能出现的最大波浪载荷。本文采用三维波浪载荷预报方法，计算了 21 个船体横剖面上的垂向剪力、垂向弯矩、扭矩和水平弯矩。首先建立波浪载荷的计算模型，如图 3 所示。图 3 船体波浪载荷计算模型网格 Fig.3 Calculation model grid for wave-induced loads of the ship 1.3 波浪载荷的传递函数响应 把船舯垂向弯矩的模型试验结果同三维理论预报结果进行对比，垂向弯矩传递函数如图 4 所示，横坐标为波长船长比，纵坐标为垂向弯矩的无因次化结果。可以看出，模型试验结果同理论预报结果吻合良好，曲线变化趋势接近，二者都在波长船长比为 1 附近时达到峰值。此时，由于计算和试验波高相对较低，船体垂向弯矩响应以低频成分为主。一般来讲，传递函数响应是作波浪载荷长期预报的基础，超越概率为108下的波浪载荷长期预报值可作为该船的波浪载荷设计值。1.4 波浪载荷长期预报沿船长分布 基于波浪载荷传递函数，对该大型舰船进行波浪载荷长期预报，海浪谱选择双参数 PM 谱24，波谱采用长峰波，海浪长期资料选择 No.3425。航向角选择 0180等概率分布，超越概率选择 108，船舶航速选择零航速。给出垂向剪力、扭矩、垂向弯矩和水平弯矩沿船长方向的分布。可以发现，船体垂向剪力和扭矩沿船长变化呈现双峰 M 的形状，垂向弯矩和水平弯矩沿船长变化呈现抛物线形状，如图 58 所示。图 4 顶浪 0kn 航速舯垂向弯矩传递函数模型 试验与三维计算对比 Fig.4 Comparison between model test and 3D calculation of vertical bending moment transfer function at a speed of 0 kn for head sea 图 5 垂向剪力 Fz长期预报 Fig.5 Long-term prediction of vertical shear force Fz 图 6 扭矩 Mx长期预报 Fig.6 Long-term prediction of torque Mx 从表 2 可知，垂向剪力在第 6 站达到最大，为56.3 MN。扭矩最大值在第 13 站取到，为 765 MNm。垂向弯矩最大值在第 11 站取到，为 3710 MNm。第 20 卷 第 9 期 耿彦超，等：舰船波浪载荷响应特性及砰击载荷影响因素分析 5 图 7 垂向弯矩 My长期预报 Fig.7 Long-term prediction of VBM My 图 8 水平弯矩 Mz长期预报 Fig.8 Long-term prediction of HBM Mz 表 2 载荷长期预报（零航速，满载，超越概率 108）Tab.2 Long-term load prediction(zero speed,full load,exceeding probability 108)of loads 站号 垂向剪力Fz/(MN)扭矩Mx/(MNm)垂向弯矩My/(MNm)水平弯矩Mz/(MNm)Sta20 0.0695 0.155 5.47 2.26 Sta19 31 124 187 177 Sta18 48.3 258 670 452 Sta17 56.1 409 1 300 808 Sta16 56.2 537 1 960 1 220 Sta15 51.5 638 2 580 1 630 Sta14 45.2 718 3 100 2 000 Sta13 37.5 765 3 470 2 190 Sta12 29.3 719 3 670 2 310 Sta11 27.1 651 3 710 2 400 Sta10 33.4 625 3 560 2 430 Sta9 41.7 539 3 250 2 420 Sta8 49.5 510 2 780 2 360 Sta7 54.8 566 2 180 2 190 Sta6 56.3 573 1 580 1 890 Sta5 50.6 571 1 230 1 520 Sta4 38.1 512 976 1 130 Sta3 23.2 431 721 736 Sta2 20.3 302 480 408 Sta1 17.6 181 233 165 Sta0 1.1 6.43 25.3 3.46 水平弯矩最大值在第 10 站取到，为 2 430 MNm。1.5 不同超越概率对波浪载荷长期预报的影响 长期预报一般采用超越概率 108，为了研究超越概率对波浪载荷的影响，图 9 给出了该船满载工况、零航速、超越概率为 1010101下船舯垂向弯矩长期预报沿船长的分布结果。图 10 给出了载荷较大位置处以及舯剖面处垂向弯矩的长期预报结果随超越概率变化曲线，超越概率为 1011010。图 9 不同超越概率下垂向弯矩长期预报 Fig.9 Long-term prediction of vertical bending moment with different exceeding probabilities 图 10 船舯附近不同超越概率下垂向弯矩长期预报 Fig.10 Long-term prediction of vertical bending moments near the midship with different exceeding probabilities 表 3 给出了长期预报的具体结果，发现随着超越概率的降低，同一剖面位置处的船体梁载荷逐渐增大，且随超越概率呈对数线性增长，101超越概率下载荷长期预报结果最小，1010载荷最大。超越概率的变化基本不影响载荷沿船长方向的分布。1.6 不同装载状态对波浪载荷的影响 针对本船开展满载、超载和压载工况下载荷长期 6 装 备 环 境 工 程 2023 年 9 月 表 3 不同超越概率下剖面载荷长期预报值 Tab.3 Long-term predicted values of wave-induced loads with different exceeding probabilities Fz/MN Mx/(MNm)My/(MNm)Mz/(MNm)超越概率 Sta6 Sta10 Sta13 Sta10 Sta11 Sta10 Sta10 101 5.81 3.43 81.2 67.5 397 382 238 102 12.4 7.32 172 143 839 806 517 103 19.4 11.4 267 221 1 300 1 250 815 104 26.6 15.6 365 301 1 770 1 700 1 130 105 33.9 20 465 383 2 260 2 170 1 450 106 41.4 24.4 568 467 2 730 2 630 1 780 107 48.8 28.9 670 551 3 220 3 090 2 110 108 56.3 33.4 774 635 3 700 3 560 2 440 109 63.9 38 879 721 4 190 4 030 2 780 1010 71.5 42.6 985 806 4 680 4 510 3 120 预报，表 4 给出了垂向弯矩沿船长方向分布，研究不同装载工况下船体梁载荷特性。可以看出，满载和超载工况下，垂向弯矩沿船纵向分布基本一致，超载工况稍大，均大于压载工况。满载垂向弯矩最大值在第11 站取到，为 3 700 MNm；超载最大值在第 11 站取到，为 3 740 MNm，占满载最大值的 101.08%；压载最大值在第 11 站取到，为 3 370 MNm，占满载的90.95%。表 4 垂向弯矩My长期预报（v=0 kn，超越概率为 108）Tab.4 Long-term prediction of vertical bending moment (v=0 kn,exceeding probability 108)MNm 站号 满载 压载 超载 Sta20 5.69 5.77 5 Sta19 185 148 178 Sta18 667 347 671 Sta17 1 290 822 1 310 Sta16 1 960 1 440 1 990 Sta15 2 580 2 070 2 620 Sta14 3 090 2 620 3 150 Sta13 3 460 3 040 3 520 Sta12 3 670 3 290 3 720 Sta11 3 700 3 370 3 740 Sta10 3 560 3 260 3 600 Sta9 3 250 2 990 3 280 Sta8 2 780 2 560 2 810 Sta7 2 180 2 070 2 220 Sta6 1 580 1 720 1 660 Sta5 1 230 1 430 1 350 Sta4 972 1 160 1 090 Sta3 719 875 815 Sta2 479 567 540 Sta1 232 265 257 Sta0 25.3 25.2 30 2 舰船非线性砰击载荷特性分析 2.1 试验船模垂向自振湿频率 在开展大型舰船波浪载荷水池模型试验前，为准确测量舰船的非线性砰击载荷，应先在水池中测量船模的垂向弯曲振动曲线。选用质量远小于模型质量的锤子（约 2 kg），用锤子击打试验模型刚性部位，以保证模型不被损坏。锤击力能促使模型在水中能够自由振动，实时记录试验模型在水池中的自由振动衰减曲线。对自由振动时域衰减曲线进行频谱分析，获得船体试验船模垂向自振湿频率26。垂向弯矩振动响应频谱如图 11 所示。试验船模一阶垂向总振动衰减曲线如图 12 所示。可以看到，船体梁的前三阶垂向总振动频率分别为 6.92、12.94、21.55 Hz。根据上述船体梁振动阻尼的计算方法，得到船体梁的振动对数衰减率 为 0.052 3。图 11 满载状态船舯测点垂向振动频谱图 Fig.11 Vibration spectrum diagram of VBM in the midship under full load condition 2.2 航速引起砰击弯矩的变化 砰击载荷一般会随着航速的增加而增大，在水池模型试验过程中，开展了多个试验航速的垂向弯矩测量。从图 13 可以看到，垂向弯矩基本随航速 第 20 卷 第 9 期 耿彦超，等：舰船波浪载荷响应特性及砰击载荷影响因素分析 7 图 12 满载状态船体梁一阶总振动曲线 Fig.12 First order vibration curve of VBM in the midship under full load condition 图 13 顶浪航行时舯无因次化垂向弯矩合成 成分随航速的变化 Fig.13 Variation of the composite component of VBM in the midship with speed in head sea 的增大而增大，在低航速时，船体并未发生砰击现象。当航速到达一定程度，就会有砰击的发生，当试验船航速为 22 kn 时达到最大，为零航速预报结果的1.31.4 倍。因此，航速是引起砰击发生的重要因素。2.3 砰击等非线性因素对船体总弯矩的影响 对模型试验中规则波大波高和高航速下船舯垂向弯矩的试验结果进行分析27，可以得到砰击等非线性对总载荷的影响。在 18 kn 航速、10.4 m 波高下开展了 5 个波浪周期的试验工况，见表 5。对于低频垂向弯矩（MVw），中拱和中垂较为接近；对于包含了砰击弯矩的合成弯矩（MVc），中拱和中垂就呈现出了明显的不对称现象，严重时，中拱成分约占中垂成分的 1/2。从弯矩的时域曲线（见图 14）可明显看出，船体的弯矩呈现特别明显的非线性特征，中垂大幅增加。从船舯的弯矩频谱图（图 15）中发现，载荷能量除了有船体存在的遭遇频率、多个倍频成分外，还有一阶船体梁垂向总振动频率成分，且会局部放大倍频成分。由于存在砰击等非线性因素，在发生砰击的试验工况，低频波浪弯矩占比较小，只有合成成分的 30%40%。由于波高较大，因此海洋环境波高是引起砰击发生的另外一个重要因素。表 5 大波高、高航速时的舯垂向弯矩响应分析 Tab.5 Response analysis of midship vertical bending moment at large wave height and high speed 工况 波高/m 波长船长比(MVw中拱/MVw中垂)/%(MVc中拱/MVc中垂)/%(MVw/MVc)/%01 10.4 0.8 99.10 92.49 76.86 02 10.4 0.9 100.73 60.91 31.62 03 10.4 1 100.70 52.01 38.05 04 10.4 1.1 100.77 67.17 37.72 05 10.4 1.2 101.03 70.90 40.05 图 14 03 工况波浪、船舯弯矩时域图 Fig.14 Time domain diagram of wave and midship bending moment under condition 03 8 装 备 环 境 工 程 2023 年 9 月 图 15 03 工况波浪、船舯弯矩频谱图 Fig.15 Wave and midship bending moment spectrum diagram under condition 03 2.4 舷侧外伸平台形状变化对砰击载荷的影响 该舰船舷侧有外伸平台，其大小、位置、角度等参数变化对砰击载荷有着重要影响。为研究其参数变化对砰击载荷大小的影响，采用非线性波浪载荷程序对该船在不同舷侧外伸平台形状下垂向波浪载荷响应进行理论计算。该程序在切片理论基础上，从工程应用的角度出发，综合考虑了船体非直壁以及外飘砰击等引起的非线性波浪载荷响应。把船体作为弹性梁，预报其在高浪级下的非线性合成弯矩，并考虑了波浪诱导弯矩与砰击弯矩相互叠加的相位问题。通过改变其宽度、角度、高度等参数，计算波浪弯矩和合成弯矩的大小，揭示其对砰击 载荷的影响规律28。研究发现，随着舷侧外伸平台宽度增大、角度变小、高度降低，砰击颤振变得剧烈。2.4.1 舷侧外伸平台宽度变化 假设舷侧外伸平台宽度分别为 0.24B、0.32B 和0.41B（B 为型宽），通过改变其宽度，建立不同的船体模型进行多工况的计算分析，研究其变化对波浪弯矩和合成弯矩的影响。从表 6 可以看出，随着平台宽度的增大，船体合成弯矩增加。在周期为 10 s 附近时，船体响应最大。2.4.2 舷侧外伸平台角度变化 分别建立舷侧外伸平台砰击角为 20、25和 30等 3 个模型，进行 18.5 m 极限有义波高下的多周期 表 6 14 m 有义波高不同舷侧外伸平台宽度下载荷计算结果 Tab.6 Calculation results of load for different side extended flaring platform widths with 14 m meaning wave height MNm 0.24B 0.32B 0.41B 波浪周期/s 波浪弯矩 合成弯矩 波浪弯矩 合成弯矩 波浪弯矩 合成弯矩 7 3 204.80 4 717.32 3 378.52 4 960.97 3 464.13 5 095.57 8 3 999.18 5 953.96 4 203.96 6 179.34 4 322.99 6 396.42 8.5 4 148.38 6 176.69 4 361.23 6 471.52 4 483.34 6 659.62 9 4 320.17 6 398.91 4 538.61 6 720.74 4 666.19 6 949.28 9.5 4 224.39 6 267.77 4 450.79 6 537.72 4 567.83 6 770.30 10 4 320.40 6 509.40 4 532.62 6 788.53 4 662.59 7 078.19 10.5 4 266.22 6 237.42 4 481.49 6 483.37 4 604.51 6 724.08 11 4 155.30 6 171.01 4 357.01 6 442.14 4 483.93 6 647.06 11.5 4 024.76 5 802.50 4 220.98 6 005.33 4 340.52 6 257.38 12 3 957.01 5 732.55 4 156.38 5 927.82 4 267.88 6 157.04 12.5 3 804.64 5 321.45 4 000.06 5 568.16 4 098.75 5 727.32 13 3 686.68 5 153.74 3 871.13 5 367.64 3 972.93 5 529.71 14 3 465.02 4 595.00 3 647.07 4 820.58 3 735.43 4 945.44 第 20 卷 第 9 期 耿彦超，等：舰船波浪载荷响应特性及砰击载荷影响因素分析 9 载荷响应计算。从表 7 可以发现，随着砰击角的增大，船体载荷逐渐减小。2.4.3 舷侧外伸平台高度变化 分别建立舷侧外伸平台开始高度在 0.65D、0.77D和 0.89D（D 为型深）的 3 个模型，在极限波高 18.5 m 下，从表 8 可以发现，随着舷侧外伸平台起始高度升高，船体弯矩响应值逐渐减小。在舷侧外伸平台起始高度为 0.65D 高度时，载荷最大。因为此时舷侧外伸平台较大，会明显增加船体的非线性砰击载荷，导致船体的合成弯矩大幅增加。表 7 18.5 m 有义波高不同舷侧外伸平台角度下载荷计算结果 Tab.7 Calculation results of load for different side extended flaring platform angles with 18.5 m meaning wave height MNm 20 25 30 波浪周期/s 波浪弯矩 合成弯矩 波浪弯矩 合成弯矩 波浪弯矩 合成弯矩 9 6 123.69 9 719.28 5 978.49 9 345.31 5 637.04 9 120.16 9.5 6 048.25 9