[go: up one dir, main page]

CN107480334A - A kind of three joint machine fish hydrodynamics method for numerical simulation - Google Patents

A kind of three joint machine fish hydrodynamics method for numerical simulation Download PDF

Info

Publication number
CN107480334A
CN107480334A CN201710563908.7A CN201710563908A CN107480334A CN 107480334 A CN107480334 A CN 107480334A CN 201710563908 A CN201710563908 A CN 201710563908A CN 107480334 A CN107480334 A CN 107480334A
Authority
CN
China
Prior art keywords
joint
fish
main driving
driving joint
vertebra
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
CN201710563908.7A
Other languages
Chinese (zh)
Other versions
CN107480334B (en
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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN201710563908.7A priority Critical patent/CN107480334B/en
Publication of CN107480334A publication Critical patent/CN107480334A/en
Application granted granted Critical
Publication of CN107480334B publication Critical patent/CN107480334B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B63SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63HMARINE PROPULSION OR STEERING
    • B63H1/00Propulsive elements directly acting on water
    • B63H1/30Propulsive elements directly acting on water of non-rotary type
    • B63H1/36Propulsive elements directly acting on water of non-rotary type swinging sideways, e.g. fishtail type
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • Ocean & Marine Engineering (AREA)
  • Farming Of Fish And Shellfish (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种三关节机器鱼流体力学数值模拟方法,本发明选定计算区域为三关节仿生机器鱼的10倍,网格划分采取结构体四面体,边界条件采用多重参考模型方法,进行数字模拟仿真分析。得出周期内,压力最大集中在头部和两侧,且当t=1/4T和t=3/4T时,即尾部进行甩尾最大幅度时,受到的压力最大,尾部水流也最大,推进力最大;依据仿真结果数据,找到了现实头部破损原因,及破损具体位置出现在嘴部前部,大约在其前部1/3处,为增加其承受力强度,在其鱼头嘴部内侧前1/3处设计三个支撑架,以改善头部破损现象。

The invention discloses a numerical simulation method for hydrodynamics of a three-joint robotic fish. The selected calculation area of the invention is 10 times that of the three-joint bionic robotic fish. The grid division adopts a structure tetrahedron, and the boundary condition adopts a multi-reference model method. Digital simulation simulation analysis. It can be concluded that within the cycle, the pressure is concentrated on the head and both sides at its maximum, and when t=1/4T and t=3/4T, that is, when the tail flicks to the maximum extent, the pressure received is the largest, and the water flow at the tail is also the largest. The force is the largest; according to the simulation result data, the cause of the actual head damage was found, and the specific location of the damage appeared in the front of the mouth, about 1/3 of the front. Three support frames are designed at the front 1/3 of the inner side to improve the phenomenon of head damage.

Description

一种三关节机器鱼流体力学数值模拟方法A numerical simulation method for hydrodynamics of a three-joint robotic fish

技术领域technical field

本发明涉及一种仿生机器鱼技术领域,特别是一种三关节机器鱼流体力学数值模拟方法。The invention relates to the technical field of bionic robotic fish, in particular to a numerical simulation method for hydrodynamics of a three-joint robotic fish.

背景技术Background technique

计算流体力学CFD就是利用现有电子计算机对实际模型进行数值模拟,用于分析流体流动和传热等物理现象,得出计算结果,再对实际模型进行指导。CFD技术从上世纪至今已经得到了快速发展,通过可以处理庞大数据量的服务器,进行分析流体中的流动特性,CFD技术使人们在计算机上通过分析来预测现实模型中的流场运动现状,得出数据,并对数据进行处理。利用流体力学理论知识及相关CFD仿真软件Fluent进行分析,可以对机器鱼的游动研究提供更深的理论研究和更快捷的工具,能够带动机器鱼的理论研究。Computational fluid dynamics (CFD) is to use the existing electronic computer to carry out numerical simulation of the actual model, to analyze the physical phenomena such as fluid flow and heat transfer, to obtain the calculation results, and then to guide the actual model. CFD technology has developed rapidly since the last century. Through the analysis of the flow characteristics in the fluid through the server that can handle a huge amount of data, CFD technology enables people to predict the current situation of the flow field movement in the real model through analysis on the computer, and obtain output data and process the data. Using the theoretical knowledge of fluid mechanics and related CFD simulation software Fluent for analysis can provide deeper theoretical research and faster tools for the swimming research of robotic fish, and can drive the theoretical research of robotic fish.

发明内容Contents of the invention

本发明的目的是提供一种三关节机器鱼流体力学数值模拟方法。The purpose of the present invention is to provide a numerical simulation method for hydrodynamics of a three-joint robotic fish.

为解决上述技术问题,本发明采用的技术方案是:In order to solve the problems of the technologies described above, the technical solution adopted in the present invention is:

一种三关节机器鱼流体力学数值模拟方法,包括以下步骤:A numerical simulation method for hydrodynamics of a three-joint robotic fish, comprising the following steps:

第一步,建立物理模型;所述三关节机器鱼包括鱼头、鱼身和鱼尾,鱼身包括固定部和摆动部,固定部内安装有用于控制机器鱼运动姿态、无线充电的控制机构,摆动部包括三个主驱动关节,分别为沿鱼身体长方向依次设置的第一主驱动关节、第二主驱动关节与第三主驱动关节,固定部一端与鱼头固定连接,另一端与第一主驱动关节固定连接,第一主驱动关节、第二主驱动关节与第三主驱动关节的长度比例为1:1:1.2;The first step is to establish a physical model; the three-joint robotic fish includes a fish head, a fish body and a fish tail, the fish body includes a fixed part and a swing part, and a control mechanism for controlling the movement posture of the robot fish and wireless charging is installed in the fixed part. The swinging part includes three main driving joints, which are the first main driving joint, the second main driving joint and the third main driving joint arranged in sequence along the body length of the fish. One end of the fixed part is fixedly connected with the fish head, and the other end is connected with the second One main driving joint is fixedly connected, the length ratio of the first main driving joint, the second main driving joint and the third main driving joint is 1:1:1.2;

所述每个主驱动关节均包括一个舵机和驱动杆,驱动杆一端与舵机的输出轴连接,另一端与相邻主驱动关节的舵机固定连接;第一主驱动关节与第二主驱动关节之间、以及第二主驱动关节与第三驱动关节之间沿鱼身体长方向设置有若干个从驱动关节,每个从驱动关节包括支撑环、脊椎骨铰链和脊椎骨铰链连接杆,支撑环上具有用于固定脊椎骨铰链的卡槽,脊椎骨铰链具有第一连接部和第二连接部,第二连接部呈U型,第二连接部卡设于支撑环的卡槽中,第一连接部通过脊椎骨铰链连接杆与前一个从驱动关节的第二连接部转动连接,第二连接部通过脊椎骨铰链连接杆与后一个从驱动关节的第一连接部转动连接,多个从驱动关节的脊椎骨铰链首尾相连形成机器鱼的脊椎。Each of the main driving joints includes a steering gear and a driving rod, one end of the driving rod is connected to the output shaft of the steering gear, and the other end is fixedly connected to the steering gear of the adjacent main driving joint; the first main driving joint is connected to the second main driving joint. Between the driving joints, and between the second main driving joint and the third driving joint, several slave driving joints are arranged along the length direction of the fish body, and each slave driving joint includes a support ring, a vertebra hinge and a vertebra hinge connecting rod, and the support ring There is a card slot for fixing the spine hinge, the spine hinge has a first connection part and a second connection part, the second connection part is U-shaped, the second connection part is clamped in the card slot of the support ring, the first connection part The second connecting part of the previous slave driving joint is rotationally connected through the vertebral hinge connecting rod, and the second connecting part is rotationally connected with the first connecting part of the latter slave driving joint through the vertebral hinge connecting rod. The vertebral hinges of multiple slave driving joints Connected end to end to form the backbone of the robotic fish.

第二步,根据第一步得到的物理模型,确定该模型对周围环境的影响,确定计算区域,进行网格划分;将第一步得到的物理模型放于一定流速的水槽中进行流体仿真分析,通过 Fluent软件仿真出具有周期内一定规律循环运动的三关节机器鱼在流动的水中的流场图和压力图,设定计算区域为三关节机器鱼大小的10倍,选用四面体网格;In the second step, according to the physical model obtained in the first step, determine the impact of the model on the surrounding environment, determine the calculation area, and perform grid division; put the physical model obtained in the first step in a water tank with a certain flow rate for fluid simulation analysis , the flow field diagram and pressure diagram of the three-joint robot fish in the flowing water are simulated by Fluent software, and the calculation area is set to 10 times the size of the three-joint robot fish, and the tetrahedron grid is selected;

第三步,设置边界条件;假设三关节机器鱼所处实验环境的水流为恒定流,计算采用MRF 法,把整个计算区域看成一个计算区域,计算三关节机器鱼周围的流体特性;The third step is to set the boundary conditions; assuming that the water flow in the experimental environment of the three-joint robot fish is a constant flow, the calculation adopts the MRF method, and the entire calculation area is regarded as a calculation area, and the fluid characteristics around the three-joint robot fish are calculated;

第四步,把划分的网格图导入计算流体力学分析软件Fluent中,进行一系列的设置求解;设定初始条件:周期T=2s,水速v=0.5m/s,对周期内关键时刻做以仿真分析,选取t=0、t=1/8T、 t=1/4T、t=5/8T、t=3/4T五个时刻的压力场图和速度流场进行分析;The fourth step is to import the divided grid diagram into the computational fluid dynamics analysis software Fluent, and carry out a series of setting solutions; set the initial conditions: period T=2s, water velocity v=0.5m/s, for the critical moment in the period For simulation analysis, select the pressure field diagram and velocity flow field at five moments t=0, t=1/8T, t=1/4T, t=5/8T, t=3/4T for analysis;

第五步,对数据进行选择分析处理,得出数字模拟分析结果。The fifth step is to select, analyze and process the data, and obtain the digital simulation analysis results.

进一步地,第五步中绘制三关节机器鱼周期内压力最大值变化和尾部水流最大值变化图,得到周期内机器鱼的最大压力点为距离鱼头嘴部前1/3处,实践也证明此处容易发生破损。Furthermore, in the fifth step, the maximum pressure change and the maximum water flow at the tail of the three-joint robot fish are drawn in the cycle, and the maximum pressure point of the robot fish in the cycle is 1/3 from the front of the fish head and mouth. Practice has also proved that Breakage can easily occur here.

进一步地,为了增加头部的抗压强度,从而减少头部的破损情况,在所述三关节机器鱼的鱼头嘴部内侧、位于嘴部前1/3处设置有三个直径为4mm的受力支撑杆。Further, in order to increase the compressive strength of the head, thereby reducing the damage of the head, on the inside of the mouth of the fish head of the three-joint robot fish, three 4mm diameter rods are arranged at the front 1/3 of the mouth. Power support rod.

本发明采用上述方案后,具有如下技术效果:After the present invention adopts the above scheme, it has the following technical effects:

本发明选定计算区域为三关节仿生机器鱼的10倍,网格划分采取结构体四面体,边界条件采用多重参考模型方法,进行数字模拟仿真分析。得出周期内,压力最大集中在头部和两侧,且当t=1/4T和t=3/4T时,即尾部进行甩尾最大幅度时,受到的压力最大,尾部水流也最大,推进力最大;依据仿真结果数据,找到了现实头部破损原因,及破损具体位置出现在嘴部前部,大约在其前部1/3处,为增加其承受力强度,在其鱼头嘴部内侧前1/3处设计三个支撑架,以改善头部破损现象。The selected calculation area of the present invention is 10 times that of the three-joint bionic robot fish, the grid division adopts the structure tetrahedron, the boundary condition adopts the method of multiple reference models, and the digital simulation simulation analysis is carried out. It can be concluded that within the cycle, the pressure is concentrated on the head and both sides at its maximum, and when t=1/4T and t=3/4T, that is, when the tail flicks to the maximum extent, the pressure received is the largest, and the water flow at the tail is also the largest. The force is the largest; according to the simulation result data, the cause of the actual head damage was found, and the specific location of the damage appeared in the front of the mouth, about 1/3 of the front. Three support frames are designed at the front 1/3 of the inner side to improve the phenomenon of head damage.

附图说明Description of drawings

图1是模拟方法的流程示意图。Figure 1 is a schematic flow chart of the simulation method.

图2a、2b分别是t=0时刻压力场图和速度流场图。Figures 2a and 2b are the pressure field diagram and the velocity flow field diagram at t=0, respectively.

图3a、3b分别是t=1/8T时刻压力场图和速度流场图。Figures 3a and 3b are pressure field diagrams and velocity flow field diagrams at t=1/8T, respectively.

图4a、4b分别是t=1/4T时刻压力场图和速度流场图。Figures 4a and 4b are pressure field diagrams and velocity flow field diagrams at t=1/4T, respectively.

图5a、5b分别是t=5/8T时刻压力场图和速度流场图。Figures 5a and 5b are pressure field diagrams and velocity flow field diagrams at t=5/8T, respectively.

图6a、6b分别是t=3/4T时刻压力场图和速度流场图。Figures 6a and 6b are pressure field diagrams and velocity flow field diagrams at t=3/4T, respectively.

图7是周期内压力最大值变化图。Fig. 7 is a diagram of the change of the maximum value of pressure within a cycle.

图8是周期内流速最大值变化图。Fig. 8 is a graph showing the variation of the maximum value of the flow velocity within a period.

图9显示了周期内压力最大点离头部顶端距离。Figure 9 shows the distance from the point of maximum pressure to the top of the head during the cycle.

图10显示了三关节机器鱼的结构。Figure 10 shows the structure of the three-joint robotic fish.

图11显示了鱼头嘴部内侧设置支撑杆的结构。Figure 11 shows the structure in which the support rods are arranged inside the mouth of the fish head.

具体实施方式detailed description

如图1所示,一种三关节机器鱼流体力学数值模拟方法,包括:As shown in Figure 1, a numerical simulation method for hydrodynamics of a three-joint robotic fish includes:

第一步,依据现实条件,建立物理模型;The first step is to establish a physical model based on realistic conditions;

如图10所示,所述三关节机器鱼包括鱼头、鱼身和鱼尾,鱼身包括固定部1和摆动部2,固定部内安装有用于控制机器鱼运动姿态、无线充电的控制机构,摆动部包括三个主驱动关节,第一主驱动关节、第二主驱动关节与第三主驱动关节沿鱼身体长方向依次设置,固定部一端与鱼头固定连接,另一端与第一主驱动关节固定连接,第一主驱动关节、第二主驱动关节与第三主驱动关节的长度比例为1:1:1.2。As shown in Figure 10, the three-joint robot fish includes a fish head, a fish body and a fish tail, and the fish body includes a fixed part 1 and a swing part 2, and a control mechanism for controlling the movement posture of the robot fish and wireless charging is installed in the fixed part. The swinging part includes three main driving joints, the first main driving joint, the second main driving joint and the third main driving joint are arranged in sequence along the length direction of the fish body, one end of the fixed part is fixedly connected with the fish head, and the other end is connected with the first main driving joint. The joints are fixedly connected, and the length ratio of the first main driving joint, the second main driving joint and the third main driving joint is 1:1:1.2.

每个主驱动关节均包括一个舵机和驱动杆,驱动杆一端与舵机的输出轴连接,另一端与相邻主驱动关节的舵机固定连接,三关节机器鱼主驱动关节加尾部整体看来是一个往复运动,其中每一运动关节也是在上一运动关节的基础上进行的往复运动,机器鱼的运动关节的整体运动轨迹是沿正弦函数变化,第二主驱动关节在第一主驱动关节的基础下做往复运动,依然遵守正弦函数变化,第一主驱动关节和第二主驱动关节在运动时存在着相位差,第三主驱动关节和第二主驱动关节一样,在第二主驱动关节的基础下做往复运动,依然遵守正弦函数变化,第二主驱动关节和第三主驱动关节在运动时存在着更大相位差。Each main driving joint includes a steering gear and a driving rod. One end of the driving rod is connected to the output shaft of the steering gear, and the other end is fixedly connected to the steering gear of the adjacent main driving joint. It is a reciprocating motion, in which each motion joint is also a reciprocating motion based on the previous motion joint. The overall motion track of the motion joints of the robot fish changes along a sinusoidal function. The second main drive joint follows the first main drive. The reciprocating motion on the basis of joints still obeys the sinusoidal function change. There is a phase difference between the first main driving joint and the second main driving joint. The third main driving joint is the same as the second main driving joint. The reciprocating motion based on the driving joint still obeys the sinusoidal function change, and there is a larger phase difference between the second main driving joint and the third main driving joint when they move.

第一主驱动关节与第二主驱动关节之间、第二主驱动关节与第三驱动关节之间、第三驱动关节与鱼尾之间沿鱼身体长方向设置有若干个从驱动关节,每个从驱动关节包括支撑环3、脊椎骨铰链4和脊椎骨铰链连接杆,支撑环上具有用于固定脊椎骨铰链的卡槽,脊椎骨铰链具有第一连接部和第二连接部,第二连接部呈U型,第二连接部卡设于支撑环的卡槽中,第一连接部通过脊椎骨铰链连接杆与前一个从驱动关节的第二连接部转动连接,第二连接部通过脊椎骨铰链连接杆与后一个从驱动关节的第一连接部转动连接,多个从驱动关节的脊椎骨铰链首尾相连形成机器鱼的脊椎;从驱动关节的运动关系是依据主驱动关节运动关节而变化,鱼头嘴部内侧、位于嘴部前1/3处设置有三个直径为4mm的受力支撑杆。Between the first main driving joint and the second main driving joint, between the second main driving joint and the third driving joint, between the third driving joint and the fish tail, several slave driving joints are arranged along the fish body length direction, each A slave drive joint includes a support ring 3, a vertebra hinge 4 and a vertebra hinge connecting rod, the support ring has a draw-in slot for fixing the vertebra hinge, the vertebra hinge has a first connecting portion and a second connecting portion, and the second connecting portion is in the shape of a U Type, the second connecting part is clamped in the card slot of the support ring, the first connecting part is rotationally connected with the second connecting part of the previous slave driving joint through the vertebral hinge connecting rod, and the second connecting part is connected with the rear joint through the vertebral hinge connecting rod. The first connecting part of a slave driving joint is rotationally connected, and the spine hinges of multiple slave driving joints are connected end to end to form the spine of the robot fish; the kinematic relationship of the slave driving joints is changed according to the movement joints of the main driving joints, the inner side of the mouth of the fish head, Located at the front 1/3 of the mouth, there are three stressed support rods with a diameter of 4mm.

第二步,根据第一步得到的物理模型,确定该模型对周围环境的影响,确定计算区域,进行网格划分;In the second step, according to the physical model obtained in the first step, determine the impact of the model on the surrounding environment, determine the calculation area, and perform grid division;

将第一步得到的物理模型放于一定流速的水槽中进行流体仿真分析,通过Fluent软件仿真出具有周期内一定规律循环运动的三关节机器鱼在流动的水中的流场图和压力图,只取关键时刻进行分析,因为计算区域范围越大,计算将越准确,计算区域范围太小将影响计算数据,综合考虑,设定计算区域为三关节机器鱼大小的10倍。Put the physical model obtained in the first step in a water tank with a certain flow rate for fluid simulation analysis, and use the Fluent software to simulate the flow field diagram and pressure diagram of the three-joint robotic fish in the flowing water with a certain regular cycle movement in the cycle. Take the critical moment for analysis, because the larger the calculation area is, the more accurate the calculation will be. If the calculation area is too small, it will affect the calculation data. Considering comprehensively, set the calculation area to be 10 times the size of the three-joint robotic fish.

四面体网格对复杂区域边界和约束情形有很强的适应能力,而且快捷,灵活,工程人员工作量小,故选用四面体网格。在网格划分时考虑到三关节机器鱼边缘附近流场计算较大,故对边界部分网格进行加密,以利于后面的计算,在对计算区域进行离散时,对不同时刻的计算区域划分方式选择同一网格划分方式,表1给出不同时刻网格划分的数目。Tetrahedral grids have strong adaptability to complex regional boundaries and constraints, and are fast and flexible, and the workload of engineers is small, so tetrahedral grids are selected. Considering that the calculation of the flow field near the edge of the three-joint robotic fish is relatively large during grid division, the grid of the boundary part is encrypted to facilitate subsequent calculations. When the calculation area is discretized, the calculation area division method at different times The same grid division method is selected, and Table 1 shows the number of grid divisions at different times.

表1不同时刻网格划分的数目Table 1 Number of grid divisions at different times

第三步,设置边界条件;假设三关节机器鱼所处实验环境的水流为恒定流,计算采用多重参考模型方法,即MRF法(moving reference frame),MRF是一种稳态算法,计算量较小, MRF优点是不需要实验数据便能实现对实验水池的整体数值模拟,适用于三关节机器鱼在水中较小的体系。模拟时,把整个计算区域看成一个计算区域,计算三关节机器鱼周围的流体特性。The third step is to set the boundary conditions; assuming that the water flow in the experimental environment of the three-joint robotic fish is a constant flow, the calculation adopts the multiple reference model method, that is, the MRF method (moving reference frame). MRF is a steady-state algorithm, and the calculation amount is relatively small The advantage of MRF is that it can realize the overall numerical simulation of the experimental tank without experimental data, and it is suitable for the small system of the three-joint robotic fish in the water. During the simulation, the entire calculation area is regarded as a calculation area, and the fluid characteristics around the three-joint robotic fish are calculated.

第四步,把划分的网格图导入计算流体力学分析软件Fluent中,进行一系列的设置求解;The fourth step is to import the divided grid diagram into the computational fluid dynamics analysis software Fluent, and perform a series of setting solutions;

设定初始条件:周期T=2s,水速v=0.5m/s。对周期内关键时刻做以仿真分析,选取了t=0、 t=1/8T、t=1/4T、t=5/8T、t=3/4T五个关键时刻的压力场图和速度流场,如图2-6显示了单位周期内关键时刻的压力场图和速度流场图。Set the initial conditions: period T = 2s, water velocity v = 0.5m/s. Perform simulation analysis on the key moments in the cycle, and select the pressure field diagram and velocity flow at five key moments t=0, t=1/8T, t=1/4T, t=5/8T, t=3/4T Field, as shown in Figure 2-6 shows the pressure field diagram and velocity flow field diagram at critical moments in the unit cycle.

三关节机器鱼在进行周期性摆动时,当水速v=0.5m/s时,由图2a可看出t=0时,最大压力点出现头部正前方,为P1=1.45×10-1Pa,由图2b可看出t=0时,最大水流出现在头部两侧,为v1=6.65×10-1m/s;When the three-joint robot fish is swinging periodically, when the water velocity v=0.5m/s, it can be seen from Figure 2a that when t=0, the maximum pressure point appears directly in front of the head, which is P 1 =1.45×10 - 1 Pa, it can be seen from Figure 2b that when t=0, the maximum water flow appears on both sides of the head, which is v 1 =6.65×10 -1 m/s;

由图3a可看出t=1/8T时,最大压力点出现头部右侧,为P2=1.49×10-1Pa,由图3b可看出t=1/8T时,最大水流出现在尾部右侧,为v2=7.46×10-1m/s;It can be seen from Figure 3a that when t=1/8T, the maximum pressure point appears on the right side of the head, which is P 2 =1.49×10 -1 Pa. It can be seen from Figure 3b that when t=1/8T, the maximum water flow appears at On the right side of the tail, v 2 =7.46×10 -1 m/s;

由图4a可看出t=1/8T时,最大压力点出现头部右侧,为P3=1.56×10-1Pa,但此时仿生机器鱼身体右侧整体都受到很大压力,由图4b可看出t=1/8T时,最大水流出现在尾部左后方,为v1=8.59×10-1m/s;It can be seen from Figure 4a that when t=1/8T, the maximum pressure point appears on the right side of the head, which is P 3 =1.56×10 -1 Pa. It can be seen from Figure 4b that when t=1/8T, the maximum water flow appears at the left rear of the tail, which is v 1 =8.59×10 -1 m/s;

由图5a可看出t=5/8T时,最大压力点出现头部左侧,为P4=1.49×10-1Pa,由图5b可看出t=5/8T时,最大水流出现在尾部左侧,为v4=8.35×10-1m/s;It can be seen from Figure 5a that when t=5/8T, the maximum pressure point appears on the left side of the head, which is P 4 =1.49×10 -1 Pa. It can be seen from Figure 5b that when t=5/8T, the maximum water flow appears at On the left side of the tail, v 4 =8.35×10 -1 m/s;

由图6a可看出t=3/4T时,最大压力点出现头部左侧,为P5=1.58×10-1Pa,由图6b可看出t=3/4T时,最大水流出现在尾部右侧,为v5=8.42×10-1m/s。因为三关节机器鱼运动是连续的,虽然上面只是对关键时刻处理,但可以将其数据按周期内关键点进行连续性处理,下面对仿真分析的数据做具体处理。It can be seen from Figure 6a that when t=3/4T, the maximum pressure point appears on the left side of the head, which is P 5 =1.58×10 -1 Pa. It can be seen from Figure 6b that when t=3/4T, the maximum water flow appears at On the right side of the tail, v 5 =8.42×10 -1 m/s. Because the movement of the three-joint robot fish is continuous, although the above is only processing the key moments, its data can be processed continuously according to the key points in the cycle. The following is the specific processing of the data of the simulation analysis.

第五步,对数据进行选择分析处理,得出数字模拟分析结果;The fifth step is to select, analyze and process the data, and obtain the digital simulation analysis results;

三关节机器鱼周期内压力最大值变化和尾部水流最大值变化分别如图7和图8所示:The maximum change of pressure and the maximum change of water flow at the tail of the three-joint robotic fish cycle are shown in Figure 7 and Figure 8, respectively:

由图7可以看出三关节机器鱼当游动周期为2s时,周期内压力最大值主要集中在P=0.140~0.160pa,其中最大值波峰出现在t=1/4T和t=3/4T时刻,最大值波谷出现在t=1/2T 时刻;由图8可以看出仿生机器鱼周围流场最大值主要集中在v=0.66~0.85m/s之间,其中最大值波峰也同样出现在t=1/4T和t=3/4T时刻,最小值波谷出现在t=1/2T时刻。It can be seen from Figure 7 that when the swimming cycle of the three-joint robotic fish is 2s, the maximum pressure in the cycle is mainly concentrated at P=0.140-0.160Pa, and the maximum peaks appear at t=1/4T and t=3/4T time, the maximum trough appears at time t=1/2T; it can be seen from Figure 8 that the maximum value of the flow field around the bionic robotic fish is mainly concentrated between v=0.66~0.85m/s, and the maximum peak also appears at At time t=1/4T and t=3/4T, the minimum value trough appears at time t=1/2T.

当t=1/4T和t=3/4T时,其三个主驱动关节摆动角度达到极限,此时仿生机器鱼在进行身体弯曲达到极限,在进行最大幅度的甩尾动作,受到的压力最大,尾部水流也最大,推进力最大;当t=1/2T,其三个主驱动关节摆动极小,接近于零,此时仿生机器鱼身体弯曲最小,受到的压力最小,尾部水流也最小。When t=1/4T and t=3/4T, the swing angles of the three main drive joints reach the limit. At this time, the bionic robot fish is bending its body to the limit, and it is performing the largest tail flick action, and the pressure it receives is the greatest. , the water flow at the tail is also the largest, and the propulsion force is the largest; when t=1/2T, the swing of the three main drive joints is very small, close to zero. At this time, the body bending of the bionic robotic fish is the smallest, the pressure it receives is the smallest, and the water flow at the tail is also the smallest.

三关节机器鱼周期内最大值主要出现头部,以鱼头顶端为0处,鱼头指向鱼尾为正方向,其压力最大点处距离鱼头顶端距离周围内变化如图9所示,由图可以看出周期内最大压力点距离头部主要集中在L=0~16.5mm处,最大时刻出现在t=3/4T处,最小时刻出现在t=1/2处。在三关节仿生机器鱼设计根据实际中金龙鱼大小,其头部长度设计为L1=125mm,嘴部长度设计为L2=50mm,可以看出L=0~16.5mm主要集中在嘴部前部,在嘴部前1/3中间,此处容易发生破损,可以在内部设计三个直径为4mm的受力支撑杆5,如图11所示,支撑杆结构以增加头部的抗压强度,从而减少头部的破损情况。The maximum value in the cycle of the three-joint robotic fish mainly occurs at the head, with the top of the fish head being 0, and the fish head pointing to the fish tail as the positive direction. It can be seen from the figure that the maximum pressure point in the cycle is mainly concentrated at L=0~16.5mm from the head, the maximum moment appears at t=3/4T, and the minimum moment appears at t=1/2. In the design of the three-joint bionic robot fish, the length of the head is designed to be L 1 =125mm and the length of the mouth is designed to be L 2 =50mm according to the actual size of the Arowana. It can be seen that L=0~16.5mm is mainly concentrated in front of the mouth In the middle of the front 1/3 of the mouth, it is easy to be damaged here, and three supporting rods 5 with a diameter of 4mm can be designed inside, as shown in Figure 11, the structure of the supporting rods can increase the compressive strength of the head , thereby reducing damage to the head.

以上述依据本发明的理想实施例为启示,通过上述的说明内容,相关工作人员完全可以在不偏离本项发明技术思想的范围内,进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容,必须要根据权利要求范围来确定其技术性范围。Inspired by the above-mentioned ideal embodiment according to the present invention, through the above-mentioned description content, relevant workers can make various changes and modifications within the scope of not departing from the technical idea of the present invention. The technical scope of the present invention is not limited to the content in the specification, but must be determined according to the scope of the claims.

Claims (5)

1.一种三关节机器鱼流体力学数值模拟方法,其特征在于:包括以下步骤:1. A three-joint robot fish hydrodynamics numerical simulation method, is characterized in that: comprises the following steps: 第一步,建立物理模型;所述三关节机器鱼包括鱼头、鱼身和鱼尾,鱼身包括固定部和摆动部,固定部内安装有用于控制机器鱼运动姿态、无线充电的控制机构,摆动部包括三个主驱动关节,分别为沿鱼身体长方向依次设置的第一主驱动关节、第二主驱动关节与第三主驱动关节,固定部一端与鱼头固定连接,另一端与第一主驱动关节固定连接,第一主驱动关节、第二主驱动关节与第三主驱动关节的长度比例为1:1:1.2;The first step is to establish a physical model; the three-joint robotic fish includes a fish head, a fish body and a fish tail, the fish body includes a fixed part and a swing part, and a control mechanism for controlling the movement posture of the robot fish and wireless charging is installed in the fixed part. The swinging part includes three main driving joints, which are the first main driving joint, the second main driving joint and the third main driving joint arranged in sequence along the body length of the fish. One end of the fixed part is fixedly connected with the fish head, and the other end is connected with the second One main driving joint is fixedly connected, the length ratio of the first main driving joint, the second main driving joint and the third main driving joint is 1:1:1.2; 第二步,根据第一步得到的物理模型,确定该模型对周围环境的影响,确定计算区域,进行网格划分;将第一步得到的物理模型放于一定流速的水槽中进行流体仿真分析,通过Fluent软件仿真出具有周期内一定规律循环运动的三关节机器鱼在流动的水中的流场图和压力图,设定计算区域为三关节机器鱼大小的10倍,选用四面体网格;In the second step, according to the physical model obtained in the first step, determine the impact of the model on the surrounding environment, determine the calculation area, and perform grid division; put the physical model obtained in the first step in a water tank with a certain flow rate for fluid simulation analysis , the flow field diagram and pressure diagram of the three-joint robot fish in the flowing water are simulated by Fluent software, and the calculation area is set to 10 times the size of the three-joint robot fish, and the tetrahedron grid is selected; 第三步,设置边界条件;假设三关节机器鱼所处实验环境的水流为恒定流,计算采用MRF法,把整个计算区域看成一个计算区域,计算三关节机器鱼周围的流体特性;The third step is to set the boundary conditions; assuming that the water flow in the experimental environment of the three-joint robotic fish is a constant flow, the calculation adopts the MRF method, and the entire calculation area is regarded as a calculation area, and the fluid characteristics around the three-joint robotic fish are calculated; 第四步,把划分的网格图导入计算流体力学分析软件Fluent中,进行一系列的设置求解;设定初始条件:周期T=2s,水速v=0.5m/s,对周期内关键时刻做以仿真分析,选取t=0、t=1/8T、t=1/4T、t=5/8T、t=3/4T五个时刻的压力场图和速度流场进行分析;The fourth step is to import the divided grid diagram into the computational fluid dynamics analysis software Fluent, and carry out a series of setting solutions; set the initial conditions: period T=2s, water velocity v=0.5m/s, for the critical moment in the period For simulation analysis, the pressure field diagram and velocity flow field at five moments t=0, t=1/8T, t=1/4T, t=5/8T, t=3/4T are selected for analysis; 第五步,对数据进行选择分析处理,得出数字模拟分析结果。The fifth step is to select, analyze and process the data, and obtain the digital simulation analysis results. 2.根据权利要求1所述的三关节机器鱼流体力学数值模拟方法,其特征在于:第五步中绘制三关节机器鱼周期内压力最大值变化和尾部水流最大值变化图,得到周期内机器鱼的最大压力点为距离鱼头嘴部前1/3处。2. The hydrodynamic numerical simulation method of the three-joint robotic fish according to claim 1, characterized in that: in the fifth step, the maximum pressure change and the maximum water flow at the tail of the three-joint robotic fish are drawn in the cycle to obtain the cycle of the machine. The maximum pressure point of the fish is 1/3 from the front of the fish head and mouth. 3.根据权利要求1所述的三关节机器鱼流体力学数值模拟方法,其特征在于:所述每个主驱动关节均包括一个舵机和驱动杆,驱动杆一端与舵机的输出轴连接,另一端与相邻主驱动关节的舵机固定连接。3. The hydrodynamic numerical simulation method for three-joint robotic fish according to claim 1, wherein each main drive joint includes a steering gear and a driving rod, and one end of the driving rod is connected to the output shaft of the steering gear, The other end is fixedly connected with the steering gear of the adjacent main drive joint. 4.根据权利要求3所述的三关节机器鱼流体力学数值模拟方法,其特征在于:第一主驱动关节与第二主驱动关节之间、以及第二主驱动关节与第三驱动关节之间沿鱼身体长方向设置有若干个从驱动关节,每个从驱动关节包括支撑环、脊椎骨铰链和脊椎骨铰链连接杆,支撑环上具有用于固定脊椎骨铰链的卡槽,脊椎骨铰链具有第一连接部和第二连接部,第二连接部呈U型,第二连接部卡设于支撑环的卡槽中,第一连接部通过脊椎骨铰链连接杆与前一个从驱动关节的第二连接部转动连接,第二连接部通过脊椎骨铰链连接杆与后一个从驱动关节的第一连接部转动连接,多个从驱动关节的脊椎骨铰链首尾相连形成机器鱼的脊椎。4. The hydrodynamic numerical simulation method for three-joint robotic fish according to claim 3, characterized in that: between the first main driving joint and the second main driving joint, and between the second main driving joint and the third driving joint Several slave drive joints are arranged along the body length of the fish, each slave drive joint includes a support ring, a vertebra hinge and a vertebra hinge connecting rod, the support ring has a draw-in slot for fixing the vertebra hinge, and the vertebra hinge has a first connecting portion And the second connecting part, the second connecting part is U-shaped, the second connecting part is clamped in the slot of the support ring, the first connecting part is rotationally connected with the second connecting part of the previous slave driving joint through the vertebra hinge connecting rod , the second connecting part is rotationally connected with the first connecting part of the latter slave driving joint through a vertebra hinge connecting rod, and a plurality of vertebrae hinges of the slave driving joint are connected end to end to form the spine of the robotic fish. 5.根据权利要求1所述的三关节机器鱼流体力学数值模拟方法,其特征在于:所述三关节机器鱼的鱼头嘴部内侧、位于嘴部前1/3处设置有三个直径为4mm的受力支撑杆。5. The numerical simulation method for hydrodynamics of a three-joint robotic fish according to claim 1, wherein three joints with a diameter of 4 mm are arranged on the inner side of the mouth of the three-joint robotic fish and located at the front 1/3 of the mouth. force support rod.
CN201710563908.7A 2017-07-12 2017-07-12 A numerical simulation method for hydrodynamics of a three-joint robotic fish Expired - Fee Related CN107480334B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710563908.7A CN107480334B (en) 2017-07-12 2017-07-12 A numerical simulation method for hydrodynamics of a three-joint robotic fish

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710563908.7A CN107480334B (en) 2017-07-12 2017-07-12 A numerical simulation method for hydrodynamics of a three-joint robotic fish

Publications (2)

Publication Number Publication Date
CN107480334A true CN107480334A (en) 2017-12-15
CN107480334B CN107480334B (en) 2019-08-20

Family

ID=60595711

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710563908.7A Expired - Fee Related CN107480334B (en) 2017-07-12 2017-07-12 A numerical simulation method for hydrodynamics of a three-joint robotic fish

Country Status (1)

Country Link
CN (1) CN107480334B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109334859A (en) * 2018-09-10 2019-02-15 中国船舶工业集团公司第七0八研究所 A kind of small waterplane area twin hull ship global vibration steel test model design method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102158259A (en) * 2010-12-06 2011-08-17 河海大学 Bluetooth communication-based robofish mobile terminal control system and method
US20150120045A1 (en) * 2013-10-24 2015-04-30 Board Of Trustees Of Michigan State University Gliding robotic fish navigation and propulsion
CN104881045A (en) * 2015-06-17 2015-09-02 中国科学院自动化研究所 Bionic robot fish three-dimensional tracking method based on embedded visual guidance
CN106741762A (en) * 2016-12-13 2017-05-31 辽宁石化职业技术学院 A kind of machine fish and its Method for Underwater Target Tracking

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102158259A (en) * 2010-12-06 2011-08-17 河海大学 Bluetooth communication-based robofish mobile terminal control system and method
US20150120045A1 (en) * 2013-10-24 2015-04-30 Board Of Trustees Of Michigan State University Gliding robotic fish navigation and propulsion
CN104881045A (en) * 2015-06-17 2015-09-02 中国科学院自动化研究所 Bionic robot fish three-dimensional tracking method based on embedded visual guidance
CN106741762A (en) * 2016-12-13 2017-05-31 辽宁石化职业技术学院 A kind of machine fish and its Method for Underwater Target Tracking

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王平 等: "仿生机器鱼运动学模型优化与实验", 《智能系统学报》 *
马尧 等: "基于CFD的仿生机器鱼数值仿真分析", 《白城师范学院学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109334859A (en) * 2018-09-10 2019-02-15 中国船舶工业集团公司第七0八研究所 A kind of small waterplane area twin hull ship global vibration steel test model design method

Also Published As

Publication number Publication date
CN107480334B (en) 2019-08-20

Similar Documents

Publication Publication Date Title
Du et al. Underwater soft robot modeling and control with differentiable simulation
US8468108B1 (en) Modeling efficiency over a range of velocities in underwater vehicles
CN110238839A (en) A control method for multi-axis hole assembly using environment prediction to optimize non-model robot
Boyer et al. Fast dynamics of an eel-like robot—comparisons with Navier–Stokes simulations
CN109325303B (en) Ship rotary motion simulation method
Benkherouf et al. Efficiency of an auto-propelled flapping airfoil
CN107480334B (en) A numerical simulation method for hydrodynamics of a three-joint robotic fish
KR101199481B1 (en) Method for generating graphic hair motion and its apparatus
Eldredge Numerical simulations of undulatory swimming at moderate Reynolds number
Liu et al. Numerical modeling and application of the effects of fish movement on flow field in recirculating aquaculture system
CN103279645A (en) Carbon nano tube molecular dynamics simulation method based on GPU parallel computation
CN118504237A (en) A method and system for numerical simulation of external flow field of bionic fish
Conti et al. A free floating manipulation strategy for Autonomous Underwater Vehicles
AONO et al. Vortical structure and aerodynamics of hawkmoth hovering
Dziedzic et al. The role of computer animation in teaching technical subjects
CN104239621A (en) Line heating automatic forming method based on curve surface unfolding
CN119578200A (en) A method, system, device and medium for predicting unsteady aerodynamic force of aircraft
CN103177162B (en) A kind of hot Forecasting Methodology of thin-wall construction dynamics based on staggered iteration coupling technique
CN105224746A (en) Based on the pulley of Adams software---the analogy method of rope type objects
Du et al. Mechanics-based analysis on flagellated robots
Deng et al. Numerical simulation of an X-wing flapping wing MAV by means of a deforming overset grid method
Bublik Fast pressure prediction along the NACA airfoil using the convolution neural network
Prabu et al. Hydrodynamic Analysis of a Flopping NACA0012 Hydrofoil and Dolphin Fish-Like Model
Liu et al. Numerical investigation on the propulsive performance of flexible flapping fins using CFD/CSD method
Haq et al. Computational fluid dynamics analysis of IPMC actuated mechanical fin through two-way fluid structure interaction approach

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190820