CN110174239A - L型造波机的拐角处造波的方法 - Google Patents

L型造波机的拐角处造波的方法 Download PDF

Info

Publication number
CN110174239A
CN110174239A CN201910555255.7A CN201910555255A CN110174239A CN 110174239 A CN110174239 A CN 110174239A CN 201910555255 A CN201910555255 A CN 201910555255A CN 110174239 A CN110174239 A CN 110174239A
Authority
CN
China
Prior art keywords
wave
paddle
propulsion device
arc
type
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
CN201910555255.7A
Other languages
English (en)
Other versions
CN110174239B (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.)
BEIJING SINFOTEK TECHNOLOGY Co Ltd
Original Assignee
BEIJING SINFOTEK TECHNOLOGY Co Ltd
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 BEIJING SINFOTEK TECHNOLOGY Co Ltd filed Critical BEIJING SINFOTEK TECHNOLOGY Co Ltd
Priority to CN201910555255.7A priority Critical patent/CN110174239B/zh
Publication of CN110174239A publication Critical patent/CN110174239A/zh
Application granted granted Critical
Publication of CN110174239B publication Critical patent/CN110174239B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M10/00—Hydrodynamic testing; Arrangements in or on ship-testing tanks or water tunnels
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13—Differential equations

Abstract

本发明L型造波机的拐角处造波的方法包括：在由纵造波板、弧形造波板和横造波板所形成图中建立一个坐标系，在不同水深下，L型造波机在拐角处造波造出目标波浪的平均速度为V x 和V y ，在正弦波的一个循环周期内，选择至少50个时间点，分别计算出在上述至少50个时间点处的每个推进装置的Vx和Vy；然后利用差分法计算顺序为1、2……n的推进装置在上述正弦波的一个循环周期内至少50个时间点的推进位移，将n个推进装置在上述至少50个时间点处的推进位移输入给控制装置，控制装置分别对n个推进装置的推进位移进行控制，重复上述循环，造出行进方向与X轴的夹角为的0.5*d*sin(wt)目标波浪。

Description

L型造波机的拐角处造波的方法
技术领域
本发明涉及了一种L型造波机的拐角处造波的方法，特别是针对拐角处的处理与造波的方法。
背景技术
目前国内的L型造波机在对其L型拐角的处理上，由于结构的限制，在其拐角处一般选择为塞堵或者空缺未处理，即使L型造波机在对其L型拐角上，也只有X方向和/或Y方向的波纹板进行运动造波，这样会使造出来的波不连续，产生断裂，严重时无法形成所需要的波纹，而现有的方法无法解决拐角对波浪实验的影响，满足不了实验需求。
本申请针对目前国内外关于L型造波机在拐角处的造波方法缺失，提出一种拐角处造波方法。该方法不需要对水动力函数进行推导求解，根据造波板的微分运动方程进行离散变换得到目标区域的造波信号，填补了关于拐角处造波的理论空白。
在造波过程中，使用坐标系中造波单元X、Y任一边进行造波，拐角不会对造波产生影响，当同时使用造波单元X、Y进行造波时，如果不对造波单元X、Y的连接拐角处进行处理，则造波会出现衍射，破坏造波的区域。在处理拐角处的推波板运动时，需要考虑拐角处与造波单元X、Y的联动作用。
发明内容
为了完成本申请的发明目的，本申请采用以下技术方案：
本发明的一种L型造波机的拐角处造波的方法， L型造波机包括：纵造波板、弧形造波 板、横造波板、推进装置和控制装置，弧形造波板为连接纵造波板和横造波板的90°圆弧的 弧形造波板，它由多段抽拉式弧形造波板组成，在纵造波板和横造波板上分别布置有若干 个推进装置，在弧形造波板的两端分别布置有一个推进装置，弧形造波板的中间布置有至 少一个推进装置，控制装置分别与上述所有推进装置相连，其中：用上述L型造波机造出行 进方向与X轴的夹角为的0.5*d*sin(wt)目标波浪的方法包括以下步骤：
（A）、在由纵造波板、弧形造波板和横造波板所形成图中建立一个坐标系
以弧形造波板的圆心为坐标系的原点，以穿过该原点并平行于纵造波板的直线为坐标系的Y轴，以穿过原点并平行于横造波板的直线为坐标系的X轴，其中0.5*d*sin(wt)中d为正弦波的波高，w为正弦波的周期频率，与纵造波板相连的推进装置的顺序定义为1,在其后的推进装置4依次定义为2、3、4……n；
（B）、在不同水深下，L型造波机在拐角处造波造出目标波浪的平均速度V x 和V y ，该平均速度与水深h成负相关，与波高d成正相关，正弦波的平均速度方程为：
其中k为色散系数，其单位为1/m，k由公式求得，xi为某个在推进 装置与弧形造波板的交点在坐标系中的X轴方向上坐标的绝对值，yi为某个在推进装置与 弧形造波板的交点在坐标系中的Y轴方向上坐标的绝对值,为正弦波的行进方向与X轴的 夹角，在正弦波的一个循环周期内，选择至少50个时间点，分别计算出在上述至少50个时间 点处的每个推进装置的Vx和Vy；
(C)、结合得到的平均速度方程，得到顺序为1的推进装置在某一时刻的推进位移：，在正弦波的一个循环周期内，计算出上述至少50个时间点处顺序为1 的推进装置的位移；
（D）、利用差分法计算顺序为2……n的推进装置在上述正弦波的一个循环周期内至少 50个时间点的推进位移，差分方程如下：式中，为下一 个顺序的推进装置的推进位移，为当前顺序推进装置的位移，t为时间步长，r为弧形造 波板的原始半径；
(E)、将n个推进装置在上述至少50个时间点处的推进位移输入给控制装置，控制装置 分别对n个推进装置的推进位移进行控制，重复上述循环，造出行进方向与X轴的夹角为 的0.5*d*sin(wt)目标波浪。
本发明的一种L型造波机的拐角处造波的方法，其中：在正弦波的一个循环周期内，选择50-100个时间点。
本发明的一种L型造波机的拐角处造波的方法，其中：在正弦波的一个循环周期内，选择50个时间点，所述的50个时间点分别为0、 1/50f、2/50f、...；1/f秒，f为正弦波的周期频率，w=2πf。
本发明的一种L型造波机的拐角处造波的方法，其中：在弧形造波板的圆弧上布置有个推进装置。
本发明的一种L型造波机的拐角处造波的方法，其中：所述弧形造波板由3-10段抽拉式弧形造波板组成。
本本方法在不需要对水动力函数进行推导求解的条件下，只根据造波板的微分运动方程进行离散变换既可以得到目标区域的造波信号，填补了关于拐角处造波的理论空白，该方法充分地考虑了拐角处造波单元X、Y的联动作用，用该方法造出正弦波与理论的正弦波的差异均小于2%，简化了计算方法，使得计算更快，达到了工程和实验的要求。
附图说明
图1为本发明L型造波机的俯向示意图，在图1中正弦波与X轴的夹角为0°；
图2为本发明L型造波机的俯向示意图，在图2中正弦波与X轴的夹角为90°；
图3为本发明L型造波机的俯向示意图，在图3中正弦波与X轴的夹角为45°；
图4为图1至图3中弧型造波板的放大示意图。
在图1至图4中，标号1为纵造波板；标号2为弧形造波板；标号3为横造波板；标号4为推进装置；标号5为控制装置。
在图1至图3中，①、②、③、④、⑤分别为顺序号为1、2、3、4、5的推进装置4。
具体实施方式
如图1和图4所示，本发明的L型造波机包括：纵造波板1、弧形造波板2、横造波板3、推进装置4和控制装置5，弧形造波板2的一侧与纵造波板1铰接，弧形造波板2的另一侧与横造波板3铰接，在纵造波板1和横造波板3上分别布置有若干个推进装置4，在弧形造波板2的两侧分别布置有推进装置4，根据需要，在弧形造波板2的圆弧上布置有三个推进装置4，弧形造波板2由多段抽拉式弧形造波板组成，如图4所示，弧形造波2由三段抽拉式弧形造波板组成，在实践中，弧形造波板2可以由3-10段抽拉式弧形造波板组成。所有推进装置4分别与控制装置5相连。
用本发明的上述L型造波机造出行进方向（如图1至图3的箭头方向）与X轴的夹角 为的0.5*d*sin(wt)目标波浪的方法包括以下步骤：
（A）、如图1至图3所示，在由纵造波板1、弧形造波板2和横造波板3所形成图中建立一个坐标系
以弧形造波板2的圆心为坐标系的原点，以穿过该原点并平行于纵造波板1的直线为坐标系的Y轴，以穿过原点并平行于横造波板3的直线为坐标系的X轴，其中0.5*d*sin(wt)中d为正弦波的波高，f为正弦波的周期频率，w=2πf，与纵造波板1相连的推进装置4的顺序定义为1,在其后的推进装置4依次定义为2、3、4、5；
（B）、在不同水深下，L型造波机在拐角处造波造出目标波浪的平均速度V x 和V y ，该平均速度与水深h成负相关，与波高d成正相关，正弦波的平均速度方程为：
其中k为色散系数，其单位为1/m，k由公式（该公式引用于《随机 波浪及其工程应用》的第250页公式7.5.21）求得，xi为某个在推进装置4与弧形造波板2的 交点在坐标系中的X轴方向上坐标的绝对值，yi为某个在推进装置4与弧形造波板2的交点 在坐标系中的Y轴方向上坐标的绝对值,为正弦波的行进方向（即图1至图3中的箭头方 向）与X轴的夹角，在正弦波的一个循环周期内，选择50个时间点，分别计算出在50个时间点 处即：0、 1/50f、2/50f、...；1/f秒处的上述5个推进装置4的Vx和Vy；
(C)、结合得到的平均速度方程，得到顺序为1的推进装置4在某一时刻的推进位移：，在正弦波的一个循环周期内，计算出上述至少50个时间点处顺序为1的 推进装置4的位移；
（D）、利用差分法计算顺序为2……5的推进装置4在上述正弦波的一个循环周期内至少 50个时间点的推进位移，差分方程如下：式中，为下一个顺 序的推进装置4的推进位移，为当前顺序推进装置4的位移，t为时间步长，r为弧形造波 板2的原始半径，根据该公式可以计算出在上述正弦波的一个循环周期内50个时间点处，顺 序为2……5的推进装置4的位移为：
(E)、将5个推进装置4在上述50个时间点处的推进位移输入给控制装置5，控制装置5分 别对5个推进装置4的推进位移进行控制，重复上述循环，造出行进方向与X轴的夹角为的 0.5*d*sin(wt)目标波浪。
以上描述是对本发明的解释，不是对发明的限定，本发明所限定的范围参见权利要求，在不违背本发明的精神的情况下，本发明可以作任何形式的修改。

Claims (5)

1.一种L型造波机的拐角处造波的方法， L型造波机包括：纵造波板（1）、弧形造波板 （2）、横造波板（3）、推进装置（4）和控制装置（5），弧形造波板（2）为连接纵造波板（1）和横造 波板（3）的90°圆弧的弧形造波板，它由多段抽拉式弧形造波板组成，在纵造波板（1）和横造 波板（3）上分别布置有若干个推进装置（4），在弧形造波板（2）的两端分别布置有一个推进 装置（4），弧形造波板（2）的中间布置有至少一个推进装置（4），控制装置（5）分别与上述所 有推进装置（4）相连，其特征在于：用上述L型造波机造出行进方向与X轴的夹角为的0.5* d*sin(wt)目标波浪的方法包括以下步骤：
（A）、在由纵造波板（1）、弧形造波板（2）和横造波板（3）所形成图中建立一个坐标系
以弧形造波板（2）的圆心为坐标系的原点，以穿过该原点并平行于纵造波板（1）的直线为坐标系的Y轴，以穿过原点并平行于横造波板（3）的直线为坐标系的X轴，其中0.5*d*sin(wt)中d为正弦波的波高，w为正弦波的周期频率，与纵造波板（1）相连的推进装置（4）的顺序定义为1,在其后的推进装置（4）依次定义为2、3、4……n；
（B）、在不同水深下，L型造波机在拐角处造波造出目标波浪的平均速度V x 和V y ，该平均速度与水深h成负相关，与波高d成正相关，正弦波的平均速度方程为：
其中k为色散系数，其单位为1/m，k由公式求得，xi为某个在推进装置 （4）与弧形造波板（2）的交点在坐标系中的X轴方向上坐标的绝对值，yi为某个在推进装置 （4）与弧形造波板（2）的交点在坐标系中的Y轴方向上坐标的绝对值，为正弦波的行进方 向与X轴的夹角，在正弦波的一个循环周期内，选择至少50个时间点，分别计算出在上述至 少50个时间点处的每个推进装置（4）的Vx和Vy；
(C)、结合得到的平均速度方程，得到顺序为1的推进装置（4）在某一时刻的推进位移：，在正弦波的一个循环周期内，计算出上述至少50个时间点处顺序为1的 推进装置（4）的位移；
（D）、利用差分法计算顺序为2……n的推进装置（4）在上述正弦波的一个循环周期内至 少50个时间点的推进位移，差分方程如下：式中，为下一个 顺序的推进装置（4）的推进位移，为当前顺序推进装置（4）的位移，t为时间步长，r为弧 形造波板（2）的原始半径；
(E)、将n个推进装置（4）在上述至少50个时间点处的推进位移输入给控制装置（5），控 制装置（5）分别对n个推进装置（4）的推进位移进行控制，重复上述循环，造出行进方向与X 轴的夹角为的0.5*d*sin(wt)目标波浪。
2.如权利要求1所述的L型造波机的拐角处造波的方法，其特征在于：在正弦波的一个循环周期内，选择50-100个时间点。
3.如权利要求2所述的L型造波机的拐角处造波的方法，其特征在于：在正弦波的一个循环周期内，选择50个时间点，所述的50个时间点分别为0、 1/50f、2/50f、...；1/f秒，f为正弦波的周期频率，w=2πf。
4.如权利要求3所述L型造波机的拐角处造波的方法，其特征在于：在弧形造波板（2）的圆弧上布置有3个推进装置（4）。
5.如权利要求4所述的L型造波机的拐角处造波的方法，其特征在于：所述弧形造波板（2）由3-10段抽拉式弧形造波板组成。
CN201910555255.7A 2019-06-25 2019-06-25 L型造波机的拐角处造波的方法 Active CN110174239B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910555255.7A CN110174239B (zh) 2019-06-25 2019-06-25 L型造波机的拐角处造波的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910555255.7A CN110174239B (zh) 2019-06-25 2019-06-25 L型造波机的拐角处造波的方法

Publications (2)

Publication Number Publication Date
CN110174239A true CN110174239A (zh) 2019-08-27
CN110174239B CN110174239B (zh) 2021-05-11

Family

ID=67698905

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910555255.7A Active CN110174239B (zh) 2019-06-25 2019-06-25 L型造波机的拐角处造波的方法

Country Status (1)

Country Link
CN (1) CN110174239B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62179631A (ja) * 1986-02-03 1987-08-06 Ishikawajima Harima Heavy Ind Co Ltd 造波装置
JPH0928600A (ja) * 1995-07-25 1997-02-04 Taisei Corp 波発生装置
CN102221446A (zh) * 2011-03-09 2011-10-19 中国海洋石油总公司 一种双板式内孤立波造波装置及其造波方法
CN105758617A (zh) * 2016-03-03 2016-07-13 中山大学 一种非线性多向不规则波和内波造波系统及其控制方法
CN106053016A (zh) * 2016-07-25 2016-10-26 中国海洋大学 一种推板造随机波的方法
CN108731913A (zh) * 2018-05-31 2018-11-02 大连理工大学 一种模拟特定非线性波浪序列的方法
CN109520700A (zh) * 2018-08-31 2019-03-26 中国船舶工业集团公司第七0八研究所 一种极端波浪的生成方法
CN109695221A (zh) * 2017-10-20 2019-04-30 交通运输部天津水运工程科学研究所 L型造波机造波信号生成方法、装置和l型造波机
CN109799065A (zh) * 2019-01-14 2019-05-24 哈尔滨工业大学 一种基于二阶造波理论产生连续聚焦波的方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62179631A (ja) * 1986-02-03 1987-08-06 Ishikawajima Harima Heavy Ind Co Ltd 造波装置
JPH0928600A (ja) * 1995-07-25 1997-02-04 Taisei Corp 波発生装置
CN102221446A (zh) * 2011-03-09 2011-10-19 中国海洋石油总公司 一种双板式内孤立波造波装置及其造波方法
CN105758617A (zh) * 2016-03-03 2016-07-13 中山大学 一种非线性多向不规则波和内波造波系统及其控制方法
CN106053016A (zh) * 2016-07-25 2016-10-26 中国海洋大学 一种推板造随机波的方法
CN109695221A (zh) * 2017-10-20 2019-04-30 交通运输部天津水运工程科学研究所 L型造波机造波信号生成方法、装置和l型造波机
CN108731913A (zh) * 2018-05-31 2018-11-02 大连理工大学 一种模拟特定非线性波浪序列的方法
CN109520700A (zh) * 2018-08-31 2019-03-26 中国船舶工业集团公司第七0八研究所 一种极端波浪的生成方法
CN109799065A (zh) * 2019-01-14 2019-05-24 哈尔滨工业大学 一种基于二阶造波理论产生连续聚焦波的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
XUWANG 等: "A numerical wave-maker for internal solitary waves with timely updated mass source/sink terms", 《EUROPEAN JOURNAL OF MECHANICS / B FLUIDS》 *

Also Published As

Publication number Publication date
CN110174239B (zh) 2021-05-11

Similar Documents

Publication Publication Date Title
CN112904842B (zh) 一种基于代价势场的移动机器人路径规划与优化方法
CN106462968B (zh) 用于校准机动车的摄影机系统的方法和装置
CN110906934B (zh) 一种基于碰撞危险系数的无人艇避障方法和系统
CN107608347A (zh) 一种分布式控制无人艇集群围追跟踪方法
CN109597417B (zh) 一种基于避碰准则的多usv群体协同避碰规划方法
CN103777639B (zh) 运动障碍环境中uuv三维航路规划方法
CN108681787A (zh) 基于改进双向快速扩展随机树算法的无人机路径优化方法
CN104570735B (zh) 一种码垛机器人弧形运动轨迹算法
CN110267832A (zh) 车辆变动控制装置
WO2018176596A1 (zh) 基于权重改进粒子群算法的无人自行车路径规划方法
CN110703752B (zh) 免疫遗传-人工势场法的无人艇双层路径规划方法
CN115328152B (zh) 基于改进人工势场法的无人车局部实时避障路径规划方法
CN114200931B (zh) 一种基于b样条曲线优化的移动机器人路径平滑方法
Wang et al. Ship route planning based on double-cycling genetic algorithm considering ship maneuverability constraint
CN105469386B (zh) 一种确定立体相机高度与俯仰角的方法及装置
CN109931943B (zh) 无人船舶全局路径规划方法及电子设备
CN104298124A (zh) 一种多拖轮协作沉井浮运操纵仿真控制方法
CN112527018A (zh) 一种欠驱动自主水下航行器的三维镇定控制方法
CN107991872A (zh) 欠驱动auv的虚拟锚泊水平面区域镇定控制方法及实现方法
CN113492907A (zh) 一种基于多点预瞄的车道保持方法及系统
CN102339478A (zh) 一种从等高线图生成数字高程模型的方法和装置
CN110174239A (zh) L型造波机的拐角处造波的方法
CN110411459B (zh) 一种使用图像序列的空间非合作目标自旋参数的辨识方法
CN110618685A (zh) 一种无人水面艇障碍物探测误差修正及安全避碰方法
Mohamed-Seghir The branch-and-bound method and genetic algorithm in avoidance of ships collisions in fuzzy environment

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