CN112883488A - 一种基于船波相对运动的砰击载荷预报方法 - Google Patents

一种基于船波相对运动的砰击载荷预报方法 Download PDF

Info

Publication number
CN112883488A
CN112883488A CN202110157867.8A CN202110157867A CN112883488A CN 112883488 A CN112883488 A CN 112883488A CN 202110157867 A CN202110157867 A CN 202110157867A CN 112883488 A CN112883488 A CN 112883488A
Authority
CN
China
Prior art keywords
ship
slamming
surface pressure
wave
water
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110157867.8A
Other languages
English (en)
Other versions
CN112883488B (zh
Inventor
焦甲龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN202110157867.8A priority Critical patent/CN112883488B/zh
Publication of CN112883488A publication Critical patent/CN112883488A/zh
Application granted granted Critical
Publication of CN112883488B publication Critical patent/CN112883488B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63B—SHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B71/00—Designing vessels; Predicting their performance
    • B63B71/20—Designing vessels; Predicting their performance using towing tanks or model basins for designing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Algebra (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Computing Systems (AREA)
  • Chemical & Material Sciences (AREA)
  • Ocean & Marine Engineering (AREA)
  • Mechanical Engineering (AREA)
  • Combustion & Propulsion (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
  • Force Measurement Appropriate To Specific Purposes (AREA)

Abstract

本发明公开了一种基于船波相对运动的砰击载荷预报方法，包括如下步骤：获取船型信息和计算工况信息；获取船体表面压力计算点位置处的砰击压力系数；采用理论计算或模型试验的方法，获取船体表面压力计算点的入水过程中的垂向相对位移和垂向相对速度，即获取船体表面压力计算点与波浪相对运动的时历曲线；根据船波相对运动的位移时历曲线，判断船体表面压力计算点的砰击发生时间；对于船体表面压力计算点位于水面以下的入水阶段的时间段，根据砰击压力的简化计算公式计算不同时刻的砰击压力值。该方法可基于理论计算或试验测量所得的船体与波浪的相对运动或相对速度信息，直接快速高效地计算出船体表面指定位置处的砰击压力载荷的时历曲线。

Description

一种基于船波相对运动的砰击载荷预报方法
技术领域
本发明涉及船舶砰击载荷预报技术领域，特别涉及一种基于船波相对运动的砰击载荷预报方法。
背景技术
船舶在恶劣海况中航行时会产生六自由度摇荡运动，剧烈的垂荡和纵摇运动会使船体艏艉等部位发生频繁的出水和入水现象，船体在入水过程中会与波浪发生相互接触与撞击，进而导致砰击现象的发生。巨大的砰击载荷不仅会引起全船结构的颤振响应，还可能会引起局部结构的破坏。因此，准确预报与评估船舶在波浪中的砰击压力载荷是十分重要的，也是船体结构强度设计与评估的重要前提条件。
船舶砰击压力载荷的预报方法包括理论计算、数值模拟和试验方法。为了简化船舶与波浪的相互作用问题，研究者们往往将船艏部入水砰击问题与船舶在波浪中的运动问题分离开并对结构的入水问题进行了独立研究。例如，采用Wagner理论和von Karman理论可以估算二维楔形体和平板的入水砰击载荷；采用相似解等解析方法可以计算二维楔形体和三维回转体的入水砰击载荷；采用基于边界元法的势流理论和基于有限体积法的CFD技术可以计算更加复杂几何形状的物体与流体相互作用的砰击压力；采用船体局部几何模型或简化后的二维剖面楔形体模型开展落体试验可以测量物面上不同位置处的入水砰击压力。
为了更准确地预报船舶在波浪中航行时的砰击压力载荷，研究者们将船舶运动及波浪力的求解与砰击载荷的求解进行耦合计算，即在船舶运动控制微分方程中引入了砰击力的贡献，甚至还考虑了船体结构的水弹性振动对于砰击载荷的影响。迄今，相关的势流理论、CFD技术也已取得了显著成绩，并可以准确预报船舶运动响应与局部砰击载荷。另一方面，在船模的艏部外飘区域或底部区域安装若干个压力传感器，通过开展水池模型耐波性试验或波浪载荷试验也可准确测量船体表面上的砰击压力。然而，上述理论或试验方法中的砰击压力载荷的计算或测量较为复杂、难度较高、工作量较大。例如，在势流理论中需要采用基于边界元法的Wagner模型求解船体剖面在入水过程中的砰击压力；在CFD方法中需要求解RANS方程和连续性方程得到船体表面上的入水砰击压力；在模型试验中需要在船体外表面布置多个压力传感器并进行高频采样从而测量砰击压力。
发明内容
本发明的目的在于克服现有技术的缺陷和不足，提供了一种基于船波相对运动的砰击载荷预报方法。该方法可基于理论计算或试验测量所得的船体与波浪的相对运动或相对速度信息，直接快速高效地计算出船体表面指定位置处的砰击压力载荷的时历曲线。
本发明的目的可以通过如下技术方案实现：一种基于船波相对运动的砰击载荷预报方法，包括如下步骤：
获取船型信息和计算工况信息；
获取船体表面压力计算点位置处的砰击压力系数；
采用理论计算或模型试验的方法，获取船体表面压力计算点的入水过程中的垂向相对位移和垂向相对速度，即获取船体表面压力计算点与波浪相对运动的时历曲线；
根据船波相对运动的位移时历曲线，判断船体表面压力计算点的砰击发生时间；
对于船体表面压力计算点位于水面以下的入水阶段的时间段，采用下述公式计算砰击压力值：
Figure BDA0002935132120000021
其中，P为砰击压力值，ρ为水密度，K为砰击压力系数；v为船体表面压力计算点与波浪的垂向相对运动速度，t为时间。
进一步地，砰击压力系数K根据下述公式计算得到：
Figure BDA0002935132120000022
其中，P0为基于数值仿真或模型试验所得到的船体表面压力计算点在整个入水过程中的砰击压力峰值，V为对应时刻的入水速度。
进一步地，采用三维线性势流理论对船波垂向相对位移和垂向相对速度进行计算，假设船体表面压力计算点的坐标为Q(xb,yb,zb)，则该点的垂向位移zp和垂向速度vp分别为：
Figure BDA0002935132120000031
Figure BDA0002935132120000032
其中，z、
Figure BDA0002935132120000033
θ分别为船舶重心位置处的垂荡、横摇角和纵摇角；
船体表面压力计算点Q(xb,yb,zb)在波浪传播方向轴上的位置为：
X＝(Ut+xb)cosβ-yb sinβ (5)
其中，U为船舶航速，β为浪向角；
船体表面压力计算点位置处对应的波浪的波面高程可表示为：
ζ(X,t)＝ζacos{ωet+k[(Ut+xb)cosβ-yb sinβ]} (6)
其中，k为波数；
则船体表面压力计算点位置处对应的波面的垂向速度可对公式(6)求导得到：
vξ(t)＝ωeζa sin{ωet+k[(Ut+xb)cosβ-yb sinβ]} (7)
则船体表面压力计算点Q(xb,yb,zb)与波浪的垂向相对位移和垂向相对运动速度分别为：
w(t)＝zp(t)-ζ(t) (8)
v(t)＝vp(t)-vζ(t) (9)。
进一步地，根据三维线性势流理论可知，船舶重心位置处的垂荡、横摇角和纵摇角随时间的变化可看作是简谐运动，可表达为：
Figure BDA0002935132120000034
其中，ζa为波幅，ωe为遭遇频率，Rz、
Figure BDA0002935132120000035
Rθ分别为单位波幅下的船舶垂荡、横摇角和纵摇角的振幅，εz、
Figure BDA0002935132120000036
εθ分别为单位波幅下的船舶垂荡、横摇角和纵摇角的初相位。
进一步地，Rz、
Figure BDA0002935132120000041
Rθ、εz、
Figure BDA0002935132120000042
εθ可以根据三维线性势流理论计算得到。
进一步地，采用船舶耐波性试验对船波垂向相对位移和垂向相对速度进行测量，通过在水池试验船模的艏部相应剖面安装电容式浪高仪，测量船模在波浪中航行时的船波垂向相对位移，对垂向相对位移进行时间求导得到垂向相对速度。
进一步地，根据船波相对运动的位移时历曲线，将其与船体表面压力计算点与静水面之间的距离作差，结果大于零对应的时间段即为船体表面压力计算点与水接触的持续时间。
进一步地，船体表面压力计算点与水接触的持续时间段包括艏部入水和艏部出水两个阶段，艏部出水过程中相对速度为负值，认为砰击压力值是零。
进一步地，艏部入水过程中船体表面压力计算点的下潜深度逐渐增加，对应于相对位移时历曲线的递增区间。
进一步地，艏部出水过程中船体表面压力计算点的下潜深度逐渐减小，对应于相对位移时历曲线的递减区间。
本发明与现有技术相比，具有如下优点和有益效果：本发明一种基于船波相对运动的砰击载荷预报方法可基于理论计算或试验测量所得的船体与波浪的相对运动或相对速度信息，直接快速高效地计算出船体表面指定位置处的砰击压力载荷的时历曲线。解决了现有的理论计算与模型试验中砰击压力载荷预报的直接计算的繁杂性、测量成本高等缺陷。
附图说明
图1是本发明实施例中基于船波相对运动的砰击载荷预报方法流程框图；
图2是本发明实施例中船体表面压力计算点所在横剖面的位置图；
图3是图2中压力计算点Q所在横剖面的正视图；
图4是本发明实施例中基于ANSYS/LS-DYNA软件模拟的计算点在入水过程中的压力变化图；
图5是本发明实施例中船舶在波浪中运动引起的计算点入水的示意图；
图6是本发明实施例中船波垂向相对位置随时间的变化图；
图7是本发明实施例中砰击发生与持续时间示意图；
图8是本发明实施例中砰击压力计算结果示意图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。
如图1所示，一种基于船波相对运动的砰击载荷预报方法，包括如下步骤：
步骤一：获取船型信息和计算工况信息。船型信息包括几何型线信息、船型主尺度、排水量、重量分布、重心位置、转动惯量、船体表面压力计算点位置等。计算工况信息包括海况参数(规则波的波长或周期、波高)、船舶航速、浪向角。
步骤二：获取船体表面压力计算点位置处的砰击压力系数。
船壳上不同位置处的砰击压力系数仅取决于与该位置相关的船壳几何形状信息，而与入水速度无关。因此，不同位置处的砰击压力系数可以预先通过理论计算或模型试验的方式获得砰击过程的压力信号，进而根据下述公式计算得到：
Figure BDA0002935132120000051
其中，K为所要计算的砰击压力系数，P0为基于数值仿真或模型试验所得到的船体表面压力计算点在整个入水过程中的砰击压力峰值，V为对应时刻的入水速度，ρ为水密度。
步骤三：采用理论计算或模型试验的方法，获取船体表面压力计算点的入水过程中的垂向相对位移和垂向相对速度，即获取船体表面压力计算点与波浪相对运动的时历曲线。
采用三维线性势流理论对船波垂向相对位移和垂向相对速度进行计算，假设船体表面压力计算点的坐标为Q(xb,yb,zb)，则该点的垂向位移zp和垂向速度vp分别为：
Figure BDA0002935132120000052
Figure BDA0002935132120000061
其中，z、
Figure BDA0002935132120000062
θ分别为船舶重心位置处的垂荡、横摇角和纵摇角；
船体表面压力计算点Q(xb,yb,zb)在波浪传播方向轴上的位置为：
X＝(Ut+xb)cosβ-yb sinβ (5)
其中，U为船舶航速，β为浪向角；
船体表面压力计算点位置处对应的波浪的波面高程可表示为：
ζ(X,t)＝ζacos{ωet+k[(Ut+xb)cosβ-yb sinβ]} (6)
其中，k为波数；
则船体表面压力计算点位置处对应的波面的垂向速度可对公式(6)求导得到：
vξ(t)＝ωeζa sin{ωet+k[(Ut+xb)cosβ-yb sinβ]} (7)
则船体表面压力计算点Q(xb,yb,zb)与波浪的垂向相对位移和垂向相对运动速度分别为：
w(t)＝zp(t)-ζ(t) (8)
v(t)＝vp(t)-vζ(t) (9)
根据三维线性势流理论可知，船舶重心位置处的垂荡、横摇角和纵摇角随时间的变化可看作是简谐运动，可表达为：
Figure BDA0002935132120000063
其中，ζa为波幅，ωe为遭遇频率，Rz、
Figure BDA0002935132120000064
Rθ分别为单位波幅下的船舶垂荡、横摇角和纵摇角的振幅，εz、
Figure BDA0002935132120000065
εθ分别为单位波幅下的船舶垂荡、横摇角和纵摇角的初相位。Rz、
Figure BDA0002935132120000066
Rθ、εz、
Figure BDA0002935132120000067
εθ可以根据三维线性势流理论计算得到。
也可以采用船舶耐波性试验对船波垂向相对位移和垂向相对速度进行直接测量，通过在水池试验船模的艏部相应剖面安装电容式浪高仪，测量船模在波浪中航行时的船波垂向相对位移，对垂向相对位移进行时间求导得到垂向相对速度。
步骤四：根据船波相对运动的位移时历曲线，判断船体表面压力计算点的砰击发生时间。
判断船体表面压力计算点Q(xb,yb,zb)与波面的相对位置。步骤三已经求出了船波相对运动的位移时历曲线，将其与船体表面压力计算点与静水面之间的距离作差，结果大于零对应的时间段即为船体表面压力计算点与水接触的持续时间。这是因为只有当水面高度超过计算点时才会产生砰击压力，否则计算点未与水发生接触，压力为零(将大气压力设为零点)。
船体表面压力计算点与水接触的持续时间段包括艏部入水和出水两个阶段，入水过程中船体表面压力计算点的下潜深度逐渐增加，对应于相对位移历时曲线的递增区间。出水过程中船体表面压力计算点的下潜深度逐渐减小，对应于相对位移历时曲线的递减区间，此时的相对速度为负值，因此不考虑此阶段的砰击载荷，也认为砰击压力值是零。
步骤五：根据步骤四所得到的船体表面压力计算点Q(xb,yb,zb)与波面的相对位置，判断船体表面压力计算点的砰击现象发生和持续时间，并计算不同时刻的砰击压力值。对于船体表面压力计算点位于水面以上的时间段及虽位于水面以下但处于出水状态的时间段，压力值为零。对于船体表面压力计算点位于水面以下的入水阶段的时间段，采用下述公式计算砰击压力值：
Figure BDA0002935132120000071
其中，P为砰击压力值，K为砰击压力系数(仅取决于物面形状，已通过步骤二得到)；v为船体表面点与波浪的垂向相对运动速度(已通过步骤三得到)。
实施例一
本实施例中所选取船型主要参数见表1，本实施例中针对1:40的缩尺模型的砰击压力载荷进行研究。如图2、3所示，选择船艏部外飘区域的艏柱所在横剖面上的点为例进行分析，船体表面压力计算点在实船随船坐标系中的具体位置为Q(127,1.23,10)，单位m。所选取的计算工况为：规则波迎浪工况，波长船长比λ/L＝1.0、波高H＝6.4m，实船航速V＝8kn，浪向角β＝0(迎浪)。
表1船型主要参数
内容 实船 模型
缩尺比 1:1 1:40
水线长L/m 132 3.3
型宽B/m 20.8 0.52
型深D/m 9.8 0.245
吃水T/m 6.5 0.1625
排水量Δ/t 10620 0.166
重心纵向位置(距艉柱)x<sub>g</sub>/m 64.11 1.60275
船壳上不同位置处的砰击压力系数仅取决于与该位置相关的船壳几何形状信息，而与入水速度无关。因此，本实施例中使用基于CFD原理的ANSYS/LS-DYNA软件对船体表面压力计算点所在的船体横剖面进行二维剖面入水过程的水动力仿真模拟，计算得到船体表面压力计算点入水过程中的砰击压力时历曲线。如图4、5所示，根据数值仿真所得到的船体表面压力计算点在整个入水过程中的砰击压力峰值P0为253.2kPa和对应时刻的入水速度V为12.1m/s，利用公式(2)可以计算得到船体表面压力计算点位置的砰击压力系数为K＝3.45。
基于三维线性势流理论计算的船舶在波浪中的运动响应，可以得到单位波幅入射波作用下的船舶重心位置处的垂荡、横摇角和纵摇角的幅频响应和初相位信息。进而根据公式(8)得到船波相对位置随时间的变化关系，见图6，然后利用公式(9)得到船体表面压力计算点的入水过程中的砰击速度。
判断船体表面压力计算点Q(xb,yb,zb)与波面的相对位置。步骤三已经求出了船波相对运动的位移时历曲线，将其与船体表面压力计算点与静水面之间的距离作差，结果大于零对应的时间段即为砰击压力的可能持续时间，见图7。这是因为只有当水面高度超过计算点时才会产生砰击压力，否则计算点未与水发生接触，压力为零(将大气压力设为零点)。
进一步地，船体表面压力计算点与水接触的持续时间段包括艏部入水和出水两个阶段。入水过程中砰击压力计算点的下潜深度逐渐增加，对应于相对位移历时曲线的递增区间，此时船壳与水发生相互撞击产生高压。出水过程中砰击压力计算点的下潜深度逐渐减小，对应于相对位移历时曲线的递减区间，此时的相对速度为负值，因此不考虑此阶段的砰击载荷，也认为砰击压力是零。
根据步骤四所得到的船体表面压力计算点Q(xb,yb,zb)与波面的相对位置，判断船体表面压力计算点是否与水发生接触，并计算不同时刻的砰击压力值。对于计算点位于水面以上的时间段及虽位于水面以下但处于出水过程的时间段，砰击压力值为零。对于计算点位于水面以下的入水过程时间段，采用公式(1)计算压力值，计算结果见图8。其中砰击压力系数K为仅取决于物面形状，已通过步骤二得到；船体表面点与波浪的垂向相对运动速度v已通过步骤三得到。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种基于船波相对运动的砰击载荷预报方法，其特征在于，包括如下步骤：
获取船型信息和计算工况信息；
获取船体表面压力计算点位置处的砰击压力系数；
采用理论计算或模型试验的方法，获取船体表面压力计算点的入水过程中的垂向相对位移和垂向相对速度，即获取船体表面压力计算点与波浪相对运动的时历曲线；
根据船波相对运动的位移时历曲线，判断船体表面压力计算点的砰击发生时间；
对于船体表面压力计算点位于水面以下的入水阶段的时间段，采用下述公式计算砰击压力值：
Figure FDA0002935132110000011
其中，P为砰击压力值，ρ为水密度，K为砰击压力系数；v为船体表面压力计算点与波浪的垂向相对运动速度，t为时间。
2.根据权利要求1所述的一种基于船波相对运动的砰击载荷预报方法，其特征在于，砰击压力系数K根据下述公式计算得到：
Figure FDA0002935132110000012
其中，P0为基于数值仿真或模型试验所得到的船体表面压力计算点在整个入水过程中的砰击压力峰值，V为对应时刻的入水速度。
3.根据权利要求1所述的一种基于船波相对运动的砰击载荷预报方法，其特征在于，采用三维线性势流理论对船波垂向相对位移和垂向相对速度进行计算，假设船体表面压力计算点的坐标为Q(xb,yb,zb)，则该点的垂向位移zp和垂向速度vp分别为：
Figure FDA0002935132110000015
Figure FDA0002935132110000013
其中，z、
Figure FDA0002935132110000014
θ分别为船舶重心位置处的垂荡、横摇角和纵摇角；
船体表面压力计算点Q(xb,yb,zb)在波浪传播方向轴上的位置为：
X＝(Ut+xb)cosβ-ybsinβ (5)
其中，U为船舶航速，β为浪向角；
船体表面压力计算点位置处对应的波浪的波面高程可表示为：
ζ(X,t)＝ζacos{ωet+k[(Ut+xb)cosβ-ybsinβ]} (6)
其中，k为波数；
则船体表面压力计算点位置处对应的波面的垂向速度可对公式(6)求导得到：
vξ(t)＝ωeζasin{ωet+k[(Ut+xb)cosβ-ybsinβ]} (7)
则船体表面压力计算点Q(xb,yb,zb)与波浪的垂向相对位移和垂向相对运动速度分别为：
w(t)＝zp(t)-ζ(t) (8)
v(t)＝vp(t)-vζ(t) (9)。
4.根据权利要求1所述的一种基于船波相对运动的砰击载荷预报方法，其特征在于，根据三维线性势流理论可知，船舶重心位置处的垂荡、横摇角和纵摇角随时间的变化可看作是简谐运动，可表达为：
Figure FDA0002935132110000021
其中，ζa为波幅，ωe为遭遇频率，Rz、
Figure FDA0002935132110000022
Rθ分别为单位波幅下的船舶垂荡、横摇角和纵摇角的振幅，εz、
Figure FDA0002935132110000023
εθ分别为单位波幅下的船舶垂荡、横摇角和纵摇角的初相位。
5.根据权利要求4所述的一种基于船波相对运动的砰击载荷预报方法，其特征在于，Rz、
Figure FDA0002935132110000024
Rθ、εz、
Figure FDA0002935132110000025
εθ可以根据三维线性势流理论计算得到。
6.根据权利要求1所述的一种基于船波相对运动的砰击载荷预报方法，其特征在于，采用船舶耐波性试验对船波垂向相对位移和垂向相对速度进行测量，通过在水池试验船模的艏部相应剖面安装电容式浪高仪，测量船模在波浪中航行时的船波垂向相对位移，对垂向相对位移进行时间求导得到垂向相对速度。
7.根据权利要求1所述的一种基于船波相对运动的砰击载荷预报方法，其特征在于，根据船波相对运动的位移时历曲线，将其与船体表面压力计算点与静水面之间的距离作差，结果大于零对应的时间段即为船体表面压力计算点与水接触的持续时间。
8.根据权利要求7所述的一种基于船波相对运动的砰击载荷预报方法，其特征在于，船体表面压力计算点与水接触的持续时间段包括艏部入水和艏部出水两个阶段，艏部出水过程中相对速度为负值，认为砰击压力值是零。
9.根据权利要求8所述的一种基于船波相对运动的砰击载荷预报方法，其特征在于，艏部入水过程中船体表面压力计算点的下潜深度逐渐增加，对应于相对位移时历曲线的递增区间。
10.根据权利要求8所述的一种基于船波相对运动的砰击载荷预报方法，其特征在于，艏部出水过程中船体表面压力计算点的下潜深度逐渐减小，对应于相对位移时历曲线的递减区间。
CN202110157867.8A 2021-02-05 2021-02-05 一种基于船波相对运动的砰击载荷预报方法 Active CN112883488B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110157867.8A CN112883488B (zh) 2021-02-05 2021-02-05 一种基于船波相对运动的砰击载荷预报方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110157867.8A CN112883488B (zh) 2021-02-05 2021-02-05 一种基于船波相对运动的砰击载荷预报方法

Publications (2)

Publication Number Publication Date
CN112883488A true CN112883488A (zh) 2021-06-01
CN112883488B CN112883488B (zh) 2022-09-20

Family

ID=76057315

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110157867.8A Active CN112883488B (zh) 2021-02-05 2021-02-05 一种基于船波相对运动的砰击载荷预报方法

Country Status (1)

Country Link
CN (1) CN112883488B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113705121A (zh) * 2021-08-26 2021-11-26 大连理工大学 一种基于运动参数的二维结构物入水砰击力分解和计算方法
CN113978645A (zh) * 2021-09-16 2022-01-28 中国人民解放军海军工程大学 常规水面舰艇全寿命周期的垂向组合弯矩数据处理方法及系统
CN116757123A (zh) * 2023-08-11 2023-09-15 南京航空航天大学 基于二维半理论的水面飞行器结构着水载荷预报方法
CN116911135A (zh) * 2023-07-24 2023-10-20 武汉理工大学 计及波浪下砰击载荷的非线性水弹性时域计算评估方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3671954A (en) * 1969-03-15 1972-06-20 Manchester Liners Ltd Slamming indicator system for marine devices
KR20050068077A (ko) * 2003-12-29 2005-07-05 삼성중공업 주식회사 선박의 슬래밍 충격손상 방지장치
US6996493B1 (en) * 2004-02-11 2006-02-07 The United States Of America As Represented By The Secretary Of The Navy Programmed computation of predicted pressure loads on marine vessels
KR20150062238A (ko) * 2013-11-28 2015-06-08 한국해양과학기술원 파랑중부가저항 계측 장치
CN106248343A (zh) * 2016-07-07 2016-12-21 哈尔滨工程大学 一种二维运动砰击试验装置
CN106528924A (zh) * 2016-09-28 2017-03-22 哈尔滨工程大学 一种应用于侧壁式气垫船的湿甲板砰击预报方法
US20170369127A1 (en) * 2014-12-30 2017-12-28 Centro Per Gli Studi Di Tecnica Navale - Cetena S.P.A. Structural monitoring system of the hull of a ship integrated with a navigation decision support system
CN108416133A (zh) * 2018-02-28 2018-08-17 中国船舶工业集团公司第七0八研究所 一种船舶尾部节能导管极限工况下屈服强度评估方法
CN108549616A (zh) * 2018-02-28 2018-09-18 哈尔滨工程大学 基于g-n波浪模型的船舶在恶劣海况中大幅运动的预报方法
US20190072211A1 (en) * 2016-02-03 2019-03-07 China National Offshore Oil Corporation Real-time analysis system for operation of working ship based on ship attitude measurement
US20190079112A1 (en) * 2017-09-13 2019-03-14 Hefring ehf. Methods and systems for wave slam monitoring
CN112307419A (zh) * 2020-10-19 2021-02-02 华南理工大学 一种船舶垂向运动响应的快速预报方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3671954A (en) * 1969-03-15 1972-06-20 Manchester Liners Ltd Slamming indicator system for marine devices
KR20050068077A (ko) * 2003-12-29 2005-07-05 삼성중공업 주식회사 선박의 슬래밍 충격손상 방지장치
US6996493B1 (en) * 2004-02-11 2006-02-07 The United States Of America As Represented By The Secretary Of The Navy Programmed computation of predicted pressure loads on marine vessels
KR20150062238A (ko) * 2013-11-28 2015-06-08 한국해양과학기술원 파랑중부가저항 계측 장치
US20170369127A1 (en) * 2014-12-30 2017-12-28 Centro Per Gli Studi Di Tecnica Navale - Cetena S.P.A. Structural monitoring system of the hull of a ship integrated with a navigation decision support system
US20190072211A1 (en) * 2016-02-03 2019-03-07 China National Offshore Oil Corporation Real-time analysis system for operation of working ship based on ship attitude measurement
CN106248343A (zh) * 2016-07-07 2016-12-21 哈尔滨工程大学 一种二维运动砰击试验装置
CN106528924A (zh) * 2016-09-28 2017-03-22 哈尔滨工程大学 一种应用于侧壁式气垫船的湿甲板砰击预报方法
US20190079112A1 (en) * 2017-09-13 2019-03-14 Hefring ehf. Methods and systems for wave slam monitoring
CN108416133A (zh) * 2018-02-28 2018-08-17 中国船舶工业集团公司第七0八研究所 一种船舶尾部节能导管极限工况下屈服强度评估方法
CN108549616A (zh) * 2018-02-28 2018-09-18 哈尔滨工程大学 基于g-n波浪模型的船舶在恶劣海况中大幅运动的预报方法
CN112307419A (zh) * 2020-10-19 2021-02-02 华南理工大学 一种船舶垂向运动响应的快速预报方法

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
JIALONG JIAO ET AL: "CFD Simulation of Ship Seakeeping Performance and Slamming Loads in Bi-Directional CrossWave", 《JOURNAL OF MARINE SCIENCE AND ENGINEERING》 *
WANG S ET AL: "Experimental and numerical study of the slamming load on the bow of a chemical tanker in irregular waves", 《OCEAN ENGINEERING》 *
任慧龙 等: "砰击载荷作用下船艏结构瞬态响应研究", 《中国舰船研究》 *
何广华 等: "S－175 船关键剖面的入水砰击模拟", 《哈尔滨工程大学学报》 *
司海龙 等: "船艏底部砰击压力概率预报方法研究", 《中国造船》 *
文志飞: "船舶砰击载荷的计算方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
王强: "船艏外飘结构砰击载荷计算研究", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》 *
王虎虎 等: "一种新型尾滑道式船载小艇收放系统可行性研究", 《船舶工程》 *
陈占阳 等: "外张砰击压力时间分布的模型试验", 《上海交通大学学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113705121A (zh) * 2021-08-26 2021-11-26 大连理工大学 一种基于运动参数的二维结构物入水砰击力分解和计算方法
CN113705121B (zh) * 2021-08-26 2024-03-22 大连理工大学 一种基于运动参数的二维结构物入水砰击力分解和计算方法
CN113978645A (zh) * 2021-09-16 2022-01-28 中国人民解放军海军工程大学 常规水面舰艇全寿命周期的垂向组合弯矩数据处理方法及系统
CN116911135A (zh) * 2023-07-24 2023-10-20 武汉理工大学 计及波浪下砰击载荷的非线性水弹性时域计算评估方法
CN116911135B (zh) * 2023-07-24 2024-02-13 武汉理工大学 计及波浪下砰击载荷的非线性水弹性时域计算评估方法
CN116757123A (zh) * 2023-08-11 2023-09-15 南京航空航天大学 基于二维半理论的水面飞行器结构着水载荷预报方法
CN116757123B (zh) * 2023-08-11 2023-11-07 南京航空航天大学 基于二维半理论的水面飞行器结构着水载荷预报方法

Also Published As

Publication number Publication date
CN112883488B (zh) 2022-09-20

Similar Documents

Publication Publication Date Title
CN112883488B (zh) 一种基于船波相对运动的砰击载荷预报方法
Kvalsvold et al. Hydroelastic modeling of wet deck slamming on multihull vessels
Yettou et al. Experimental study on the water impact of a symmetrical wedge
CN110096734B (zh) 一种浅水中大型船舶阻力和流场特性的分析方法及系统
CN109655229A (zh) 一种小水线面双船体波浪载荷试验模型设计及测试方法
Zong et al. Evolution of slamming load and flow field in water-entry process of trimaran ship section
Jiao et al. Experimental and numerical analysis of hull girder vibrations and bow impact of a large ship sailing in waves
Whelan Wetdeck slamming of high-speed catamarans with a centre bow
Guo et al. Numerical simulation and experimental studies on aft hull local parameterized non-geosim deformation for correcting scale effects of nominal wake field
Hossain et al. Forces, ship motions and velocity wake field for KRISO container ship model in regular head waves
Ling et al. Free surface effect on the hydrodynamics of an underwater vehicle hullform, The DARPA SUBOFF
Elhadad et al. Experimental and Cfd Resistance Validation of Naval Combatant Dtmb 5415 Model
Van’t Veer et al. Bilge keel loads and hull pressures created by bilge keels fitted to a rotating cylinder
Veen et al. A 2D Smoothed Particle Hydrodynamics theory for calculating slamming loads on ship hull sections
Zheng et al. Simulation and experimental research on hydrodynamic performance of fully attached submarine in the surface and near the surface modes
Fonfach et al. Numerical study of the hydrodynamic interaction between ships in viscous and inviscid flow
Hu et al. Investigation on the effect of bow configuration on motion and green water for tumblehome vessels based on cartesian grid finite difference method
Chen et al. Experimental study of the effect of internal waves on the rotational hydrodynamics of underwater vehicle
Brizzolara et al. Wind heeling moments on very large ships. Some insights through CFD results
Du et al. Influence of the draft to ship dynamics in the virtual tank based on openfoam
Jiao et al. Slamming and Green Water Loads on Bow-Flare Ship in Regular Head Waves Investigated by Hydroelasticity Theory and Experiment
Wang et al. Numerical simulations of KCS parametric rolling in head waves
Li et al. Influence of wave amplitude on seakeeping performance of a high-speed trimaran with T-foil in head seas
Ling et al. Numerical prediction and analysis of motion response of high speed planing craft in regular waves
Lungu CFD prediction of ship-bank interaction

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant