CN105759632B - 一种水下机器人伺服阀控制液压推进器动力学仿真系统及其仿真方法 - Google Patents

Abstract

本发明公开了一种水下机器人伺服阀控制液压推进器动力学仿真系统及其仿真方法。包括伺服阀传递函数模块、伺服阀流量计算模块、恒压变量泵模块、伺服阀控制液压推进器动力学系统模块、螺旋桨推力计算模块、螺旋桨转矩计算模块、螺旋桨推力与转矩系数计算模块和螺旋桨进速系数计算模块；本发明可较真实的模拟通过电液伺服阀来控制液压马达两端油液流量、压力、扭矩和转速的变化，进而调节螺旋桨输出转矩和推力的动态过程；可较真实的模拟螺旋桨负载连续变化条件下，液压马达的流量、压力扭矩和转速的动态响应过程。本发明对作业型水下机器人的伺服阀控制液压推进器动力学仿真系统提供了一种可行、有效的解决方法。

Description

一种水下机器人伺服阀控制液压推进器动力学仿真系统及其 仿真方法

技术领域

本发明属于无人水下机器人运动控制的技术领域，尤其涉及针对作业型遥控水下机器人(ROV—Remotely Operated Vehicle)的，一种水下机器人伺服阀控制液压推进器动力学仿真系统及其仿真方法。

背景技术

作业型ROV的最大特点是能在危险深海环境中完成高强度、大负荷和高精度作业，它是我国实施深海资源开发战略不可或缺的重大技术装备之一。开发作业型ROV仿真训练模拟器，可在虚拟环境中对作业人员进行岗位培训，并对水下作业进行预演，提高水下作业效率和保证作业安全，具有重要的工程价值和现实意义。

本发明中所涉及的作业型ROV安装8台液压螺旋桨推进器，其中水平推进器4台，垂直推进器4台。为了实现作业型ROV的六自由度运动控制，需要控制8台推进器同时动作和输出推力。由于液压推进器动力学系统的执行存在较大的惯性和滞后特性，研究伺服阀控制液压推进器的动力学仿真系统，对研究作业型ROV实际装备和仿真训练模拟器具有重要的工程价值，也成为无人水下机器人控制系统研究的关键技术之一。

由于螺旋桨的推力与其进水速度、液压马达的转速和压力具有复杂的非线性关系，螺旋桨的负载转矩与其水流进速、液压马达转速和流量也具有复杂的非线性关系，按照传递函数的方法难以建立水下机器人期望推力和螺旋桨推力之间的线性化数学模型。因此，国内相关领域专家对关于伺服阀控制液压推进器动力学仿真系统的研究成果较少。同时，国外很少有公开关于作业型ROV伺服阀控制液压推进器动力学系统仿真技术方面的相关文献与资料。

本发明提出了一种水下机器人伺服阀控制液压推进器动力学系统的仿真方法，建立了一种适合作业型水下机器人的伺服阀控制液压推进器的动力学仿真系统。本发明对液压推进器动力学仿真系统的研究提供了一种有效和实用的解决方法，对作业型ROV仿真训练模拟器和ROV实际作业装备的研制都具有重要的工程价值和指导意义。

发明内容

本发明的目的是提供一种准确直观、逻辑清晰的水下机器人伺服阀控制液压推进器动力学仿真系统。本发明的目的还包括提供一种实时性强的水下机器人伺服阀控制液压推进器动力学仿真系统的仿真方法。

一种水下机器人伺服阀控制液压推进器动力学仿真系统，其特征在于：包括伺服阀传递函数模块、伺服阀流量计算模块、恒压变量泵模块、伺服阀控制液压推进器动力学系统模块、螺旋桨推力计算模块、螺旋桨转矩计算模块、螺旋桨推力与转矩系数计算模块和螺旋桨进速系数计算模块；

将与期望控制推力成正比的控制电流i_v输入到伺服阀传递函数模块，输出一个伺服阀阀芯开口位移值x_v；将伺服阀阀芯开口位移值x_v输入到伺服阀流量计算模块，得到伺服阀流量Q_v输入给伺服阀控制液压推进器动力学系统模块；将恒压变量泵的工作压力p_r输入到伺服阀流量计算模块后输送给伺服阀控制液压推进器动力学系统模块；伺服阀控制液压推进器动力学系统模块根据接收的伺服阀流量Q_v和t-1时刻的马达油液压力p_t-1和角速度ω_t-1，得到马达的由液压力p_t、输出扭矩M_g和角速度ω_t；

将t-1时刻的螺旋桨的进水速度V_t-1和螺旋桨的转速n_t-1输入到螺旋桨进速系数计算模块，得到t时刻螺旋桨的进速系数J_t；将t时刻螺旋桨的进速系数J_t输入到螺旋桨推力与转矩系数计算模块，得到t时刻螺旋桨的转矩系数K_M和推力系数K_T；

将t时刻螺旋桨的推力系数K_T、马达的角速度ω_t输入到螺旋桨推力计算模块得到t时刻螺旋桨推力T_t；将t时刻螺旋桨的转矩系数K_M、马达的转速ω_t输入到螺旋桨转矩计算模块，得t时刻螺旋桨的转矩M_t。

本发明一种水下机器人伺服阀控制液压推进器动力学仿真系统的仿真方法，

包括以下步骤：

步骤一：将与期望控制推力成正比的控制电流i_v输入到伺服阀传递函数模块，输出一个伺服阀阀芯开口位移值x_v；

步骤二：根据伺服阀阀芯开口位移值x_v，得到伺服阀流量Q_v：

Q_v＝K_vx_v

其中，K_v——伺服阀的流量增益系数；

步骤三：计算当前t时刻的马达负载流量Q_t：

Q_t＝K_qx_v-K_cp

其中：p——液压马达两端油液压力差；K_q——伺服阀流量增益；K_c——伺服阀流量压力增益；

步骤四：由当前t时刻的马达负载流量Q_t、t-1时刻的马达油液压力p_t-1和角速度ω_t-1，根据液压马达流量连续性方程确定当前t时刻的马达油液压力增加速度

步骤五：由马达油液压力增加速度采样时间Δt和马达油液压力p_t-1，确定当前t时刻马达油液压力p_t：

步骤六：判断马达油液压力p_t是否到达恒压变量泵的工作压力p_r，若达到，则马达油液压力p_t＝p_r，否则马达油液压力p_t＝p_t；

步骤七：根据当前t时刻马达压力p_t，计算t时刻液压马达输出扭矩M_g：

M_g＝D_Mp_t；

步骤八：由t时刻液压马达的输出扭矩M_g、t-1时刻液压马达的负载转矩M_t-1和角速度ω_t-1，得到液压马达的角加速度

其中：M——液压马达的外负载转矩，即为螺旋桨转矩；液压马达的角加速度；J_M——螺旋桨和液压马达轴的转动惯量；M_t——液压马达的负载转矩；

步骤九：由马达角加速度采样时间Δt和t-1时刻的马达角速度ω_t-1，计算t时刻马达角速度ω_t：

步骤十：根据t-1时刻的螺旋桨的进水速度V_t-1和螺旋桨的转速n_t-1，得到t时刻螺旋桨的进速系数J_t：

其中，D——螺旋桨的直径；

步骤十一：根据t时刻螺旋桨的进速系数J_t，通过螺旋桨的敞水特性曲线得到t时刻螺旋桨的转矩系数K_M和推力系数K_T；

步骤十二：根据t时刻螺旋桨的推力系数K_T、螺旋桨的角速度ω_t得到t时刻螺旋桨推力T_t：

T_t＝K_Tρn_t ²D⁴

步骤十三：根据t时刻螺旋桨的转矩系数K_M、螺旋桨的转速ω_t，得到t时刻螺旋桨的转矩M_t：

M_t＝K_Mρn_t ²D⁵；

步骤十四：按照p_t-1＝p_t、Q_t-1＝Q_t、M_t-1＝M_t、ω_t-1＝ω_t重新赋初值，返回步骤一。

一种水下机器人伺服阀控制液压推进器动力学仿真系统的仿真方法，还可以包括：

1、伺服阀传递函数模块中的伺服阀传递函数为：

其中，X_v(s)——为伺服阀阀芯位移x_v的拉氏变换；I(s)——为伺服阀输入电流i_v的拉氏变换；K_sv——伺服阀的放大系数；T_v——伺服阀的时间常数。

2、所述的液压马达流量连续性方程为：

其中：Q——液压马达负载油液流量；D_M——液压马达排量；V_t——阀腔、马达腔和连接管道的总容积；β_e——等效体积弹性模量；——液压马达的压力增量；C_tc——液压马达的总泄露系数；θ——液压马达角位移；——液压马达角速度。

有益效果：

本发明将螺旋桨负载转矩引入到液压马达的力矩平衡方程中，使用步步迭代的方法实时求解液压马达的流量、压力、扭矩和转速状态；同时由马达转速和螺旋桨进水速度来实时求解螺旋桨的进速系数；最后将马达转速和螺旋桨进速系数代入螺旋桨的负载转矩和推力方程中，从而实时求解螺旋桨的负载转矩和推力。本方法可以在考虑液压马达转速和螺旋桨负载转矩实时变化的情况下，能够较准确的描述电液伺服阀控制液压推进器动力学系统的动态响应过程。本发明中的水下机器人伺服阀控制液压推进器的动力学仿真系统具有准确直观和逻辑清晰的特点，具有能够较真实的模拟实际作业型ROV的电液伺服阀控制液压推进器的功能及动态响应过程的优点。

附图说明

图1伺服阀控制液压推进器动力学仿真系统的结构原理图；

图2液压推进器的流量、压力、扭矩、转速、转矩和推力的迭代求解方法。

具体实施方式

下面将结合附图对本发明作进一步详细说明。

本发明目的在于提供一种伺服阀控制液压推进器动力学系统的仿真方法，可应用于作业型ROV仿真训练模拟器的液压推进器动力学仿真系统。如图1所示：本发明包括：伺服阀传递函数模块、伺服阀流量计算模块、恒压变量泵模块、伺服阀控制液压推进器动力学系统模块、螺旋桨推力计算模块、螺旋桨转矩计算模块、螺旋桨推力与转矩系数计算模块、螺旋桨进速系数计算模块等。本发明将螺旋桨负载转矩引入到液压马达的力矩平衡方程中，使用步步迭代的方法实时求解液压马达的流量、压力、扭矩和转速状态；同时由马达转速和螺旋桨进水速度来实时求解螺旋桨的进速系数；最后将马达转速和螺旋桨进速系数代入螺旋桨的负载转矩和推力方程中，从而实时求解螺旋桨的负载转矩和推力。本方法可以在考虑液压马达转速和螺旋桨负载转矩实时变化的情况下，能够较准确的描述电液伺服阀控制液压推进器动力学系统的动态响应过程。本发明中的水下机器人伺服阀控制液压推进器的动力学仿真系统具有准确直观和逻辑清晰的特点，具有能够较真实的模拟实际作业型ROV的电液伺服阀控制液压推进器的功能及动态响应过程的优点。

本发明的目的是这样实现的，如图1和图2所示：

第一步，与期望控制推力成正比的控制电流i_v输入到“伺服阀传递函数”模块，按照下式(1)计算，输出一个伺服阀阀芯开口位移值x_v。

伺服阀传递函数可由下式来表示：

式中：X_v(s)——为伺服阀阀芯位移x_v的拉氏变换；I(s)——为伺服阀输入电流i_v的拉氏变换；K_sv——伺服阀的放大系数；T_v——伺服阀的时间常数。

第二步，x_v输入到“伺服阀流量计算”模块，按照下式(2)计算，伺服阀输出一个液压油流量Q_v。

恒压泵的工作压力为p₀，伺服阀的阀芯位移为x_v，则伺服阀的流量Q_v可按下式计算：

Q_v＝K_vx_v (2)

式中：K_v——伺服阀的流量增益系数。

第三步，将上一t-1时刻的液压马达油液压力p_t-1，根据下式(3)计算得到当前t时刻的马达负载流量Q_t。

伺服阀控制液压马达的负载流量可用下式表示：

Q_t＝K_qx_v-K_cp (3)

式中：p——液压马达两端油液压力差；x_v——伺服阀阀芯移动距离；K_q——伺服阀流量增益；K_c——伺服阀流量压力增益。

第四步，由当前t时刻的马达负载流量Q_t、t-1时刻的马达油液压力p_t-1和角速度ω_t-1，根据下式(4)计算当前t时刻的马达油液压力增加速度

液压马达流量连续性方程可以写为下式：

式中：Q——液压马达负载油液流量；D_M——液压马达排量；V_t——阀腔、马达腔和连接管道的总容积；β_e——等效体积弹性模量；——液压马达的压力增量；C_tc——液压马达的总泄露系数；θ——液压马达角位移；θ——液压马达角位移；——液压马达角速度。

第五步，由马达油液压力增加速度采样时间Δt和马达油液压力p_t-1，计算当前t时刻马达油液压力p_t。

第六步，判断p_t是否到达液压泵的工作压力p_r，若达到，则马达压力p_t＝p_r，否则马达压力p_t＝p_t。

第七步，由当前t时刻马达压力p_t，根据下式(5)计算t时刻液压马达输出扭矩M_g。

液压马达的输出力矩M_g可按下式计算：

M_g＝D_Mp (5)

第八步，由t时刻液压马达的输出扭矩M_g、t-1时刻液压马达的负载转矩M_t-1和角速度ω_t-1，由下式(6)计算液压马达的角加速度

液压马达和负载力矩的平衡方程为：

式中：M——液压马达的外负载转矩，即为螺旋桨转矩；——液压马达的角加速度；J_M——螺旋桨和液压马达轴的转动惯量；M_t——液压马达的负载转矩。

液压马达(螺旋桨)的转速n和角速度ω存在如下关系：

第九步，由马达角加速度采样时间Δt和t-1时刻的马达角速度ω_t-1，计算t时刻马达角速度ω_t。

第十步，由t时刻的ω_t和转矩系数K_M，按照下式(7)-(9)计算t时刻液压螺旋桨的负载转矩M_t。

螺旋桨转矩M_t可用下式计算：

M_t＝K_Mρn²D⁵ (8)

式中：K_M——螺旋桨转矩系数；ρ——海水密度；n——螺旋桨的转速。

可由t时刻螺旋桨的进速系数，通过螺旋桨的敞水特性曲线，得到t时刻螺旋桨的转矩系数K_M和推力系数K_T。

t时刻螺旋桨的进速系数J按下式计算：

式中：n——螺旋桨的转速；V_A——螺旋桨的进速；即为螺旋桨的纵向进水速度；D——螺旋桨的直径。

第十一步，由t时刻马达角速度ω_t和实时推力系数K_T，按照下式(7)和(10)计算t时刻螺旋桨的推力T_t。

螺旋桨推力T_t可用下式计算：

T＝K_Tρn²D⁴ (10)

式中：K_T——螺旋桨推力系数。

螺旋桨的角速度ω和角位移θ存在如下关系：

第十二步，按照p_t-1＝p_t、Q_t-1＝Q_t、M_t-1＝M_t、ω_t-1＝ω_t重新赋初值。

一种水下机器人伺服阀控制液压推进器动力学系统的仿真方法，建立了一种适合水下机器人的伺服阀控制液压推进器的动力学仿真系统，提出了一种液压推进器的流量、压力、扭矩、转速、转矩和推力的迭代求解方法。本发明可较真实的模拟通过电液伺服阀来控制液压马达两端油液流量、压力、扭矩和转速的变化，进而调节螺旋桨输出转矩和推力的动态响应过程；可较真实的模拟螺旋桨负载连续变化条件下，液压马达的流量、压力和转速的动态响应过程。本发明对水下机器人的伺服阀控制液压推进器动力学仿真系统提供了一种有效和实用的解决方法。

一种液压推进器的流量、压力、扭矩、转速、螺旋桨转矩和推力的迭代计算方法，螺旋桨推力与其进水速度、液压马达流量和转速具有复杂的非线性关系，螺旋桨负载转矩与其进水速度、液压马达压力和转速也具有复杂的非线性关系，按照传递函数的方法难以建立水下机器人期望推力和螺旋桨输出推力之间的线性关系模型。本发明将螺旋桨实时负载转矩引入到液压马达的力矩平衡方程中，使用步步迭代的方法实时求解液压马达的流量、压力、扭矩和转速，进而实时求解螺旋桨的负载转矩和推力。

一种适合水下机器人的伺服阀控制液压推进器的动力学仿真系统，建立电液伺服阀的传递函数模型，将电流控制信号转换成阀芯位移信号；建立阀芯位移与液压阀流量的传递函数模型，将阀芯位移信号转换成伺服阀流量信号；建立伺服阀流量和负载流量的关系模型，将伺服阀流量转化为负载流量；建立液压螺旋桨推进器动力学系统模型，将马达负载流量转化为液压马达的油液压力、扭矩和转速，从而使螺旋桨输出转矩和推力。

可较真实模拟螺旋桨负载连续变化条件下，液压马达的流量、压力、扭矩和转速的动态响应过程，在输入控制电流信号不变的条件下，当螺旋桨负载减小时，螺旋桨负载扭矩降低，从而使马达转速增加，马达两端油液压力降低，最终负载流量增加，转速进一步增加；当马达负载流量增加到伺服阀的最大输入流量时，马达负载流量不再增加，马达转速也不再增加，保持恒定转速运行。当螺旋桨负载增加时，马达转速、压力和流量的变化规律和以上过程相反。

本发明属于无人水下机器人运动控制技术领域，具体涉及一种作业型水下机器人(ROV—Remotely Operated Vehicle)伺服阀控制液压推进器动力学系统的仿真方法。本发明目的在于建立了一种适合水下机器人的伺服阀控制液压推进器的动力学仿真系统，并提出了一种液压推进器的压力、扭矩、转速、螺旋桨转矩和推力的迭代求解方法。本发明包括：伺服阀传递函数模块、伺服阀流量计算模块、恒压变量泵模块、伺服阀控制液压推进器动力学系统模块、螺旋桨推力计算模块、螺旋桨转矩计算模块、螺旋桨推力与转矩系数计算模块、进速系数计算模块等。本发明可较真实的模拟通过电液伺服阀来控制液压马达两端油液流量、压力、扭矩和转速的变化，进而调节螺旋桨输出转矩和推力的动态过程；可较真实的模拟螺旋桨负载连续变化条件下，液压马达的流量、压力扭矩和转速的动态响应过程。本发明对作业型水下机器人的伺服阀控制液压推进器动力学仿真系统提供了一种可行、有效的解决方法。

一种水下机器人伺服阀控制液压推进器动力学仿真系统的仿真方法：

第一步，如图2所示，按照公式p₀＝0、Q₀＝0、M₀＝0、ω₀＝0对液压推进器动力学仿真系统进行初始化。

第二步，如图1所示，将控制电流i_v输入到“伺服阀传递函数”模块，可以计算得到t时刻伺服阀的开口位移值x_v，其中控制电流i_v和推进器的期望控制推力成正比。

第三步，如图1-2所示，将伺服阀的开口位移值x_v输入到“伺服阀流量计算”模块，可以计算得到t时刻伺服阀流量Q_v。

第四步，如图2所示，将伺服阀流量Q_v输入到“计算马达负载流量”模块，并根据t-1时刻的液压马达油液压力p_t-1，由下式(12)计算当前t时刻的马达负载流量Q_t。

Q_t＝K_qx_v-K_cp_t-1 (12)

第五步，如图2所示，将马达负载流量Q_t输入到“计算马达油液压力增加速度”模块，并根据t-1时刻的液压马达油液压力p_t-1、马达角速度ω_t-1，由下式(13)计算当前t时刻的马达负载压力增加速度

第六步，如图2所示，将马达负载压力增加速度输入到“计算马达油液压力”模块，由下式(14)计算当前t时刻的马达的油液压力p_t。

第七步，如图2所示，将当前t时刻的马达的油液压力p_t输入到“判断当前是否达到当前液压泵的压力”模块。若达到当前工作压力，则p_t＝p_r；否则p_t＝p_t。

第八步，如图2所示，t时刻的马达的油液压力p_t输入到计算马达输出扭矩模块，由下式(15)计算当前t时刻的马达的输出扭矩M_g。

M_g＝D_Mp_t (15)

第九步，如图2所示，将t时刻的马达的输出扭矩M_g、t-1时刻的螺旋桨负载转矩M_t-1和螺旋桨转动角速度ω_t-1输入到“计算马达角加速度模块”，按照下式(16)计算t时刻液压马达的角加速度

式中：

第十步，如图2所示，将t时刻的液压马达的角加速度输入到“计算马达角速度”模块，按照下式(17)计算t时刻液压马达的角速度ω_t。

第十一步，如图1所示，将t-1时刻的螺旋桨的进水速度V_t-1和螺旋桨的转速n_t-1输入到“进速系数计算”模块，按照下式(18)计算t时刻螺旋桨的进速系数J_t：

式中：

第十二步，如图1所示，将t时刻螺旋桨的进速系数J_t分别输入到“螺旋桨推力系数计算”模块和“螺旋桨转矩系数计算”模块，按照对应螺旋桨的敞水特性曲线计算t时刻螺旋桨的转矩系数K_M和推力系数K_T。

第十三步，如图2所示，将t时刻螺旋桨的推力系数K_T、螺旋桨的角速度ω_t输入到“计算螺旋桨推力”模块，按照下式(19)-(20)，可得t时刻螺旋桨推力T_t。

T_t＝K_Tρn_t ²D⁴ (20)

第十四步，如图2所示，将t时刻螺旋桨的转矩系数K_M、螺旋桨的转速ω_t分别输入到“计算螺旋桨转矩”模块，按照下式(21)可得t时刻螺旋桨的转矩M_t。

M_t＝K_Mρn_t ²D⁵ (21)

第十五步，按照下式(22)-(25)给液压推进器动力学仿真系统重新赋初值。

p_t-1＝p_t (22)

Q_t-1＝Q_t (23)

M_t-1＝M_t (24)

ω_t-1＝ω_t (25)

在以上模块各自功能调试及联调无误后，设定一个螺旋桨进速初始值，给一个期望推力值对应的控制电流i_v。在Matlab中运行仿真程序后，在输入电流控制不变的条件下，当螺旋桨负载减小时，螺旋桨负载扭矩降低，从而使马达转速增加，马达两端油液压力降低，最终负载流量增加，转速进一步增加；当马达负载流量增加到伺服阀的最大输入流量时，马达负载流量不再增加，马达转速也不再增加，保持恒定转速运行。当螺旋桨负载增加时，马达转速、压力和流量的变化规律和以上过程相反。在Matlab中运行仿真程序后，当输入与期望推力成正比的控制电流时，螺旋桨可输出与期望推力值一致的推力。本发明可较真实模拟螺旋桨负载连续变化条件下，液压马达的转速、压力及流量的动态响应过程。本发明最终应用于作业型ROV仿真训练模拟器的设计中，可较真实的模拟在8台液压推进器的共同控制作用下，作业型ROV六自由度的运动和姿态的动态响应过程。

符号说明：

i_v——伺服阀控制电流；x_v——伺服阀开口位移值；p_r——恒压泵的工作压力；Q_v——伺服阀的流量；ω——液压马达的转动角速度；M——螺旋桨的负载转矩；T——螺旋桨的推力；K_M——螺旋桨转矩系数；K_T——螺旋桨推力系数；J——螺旋桨的进速系数；V_A——螺旋桨的水流进速。

Q_t——t时刻的负载流量；——马达油液压力增加速度；p_t——t时刻的马达负载压力；M_g——t时刻的马达输出扭矩；——液压马达的角加速度；ω_t——t时刻的液压马达角速度；T_t——t时刻的螺旋桨推力；Q_t-1——t-1时刻的负载流量；p_t-1——t-1时刻的马达负载压力；ω_t-1——t-1时刻的液压马达角速度；T_t-1——t-1时刻的螺旋桨推力；Q₀——初始时刻的负载流量；p₀——初始时刻的马达负载压力；ω₀——初始时刻的液压马达角速度；M₀——初始时刻螺旋桨的负载转矩。

Claims (4)

1.一种水下机器人伺服阀控制液压推进器动力学仿真系统，其特征在于：包括伺服阀传递函数模块、伺服阀流量计算模块、恒压变量泵模块、伺服阀控制液压推进器动力学系统模块、螺旋桨推力计算模块、螺旋桨转矩计算模块、螺旋桨推力与转矩系数计算模块和螺旋桨进速系数计算模块；

将与期望控制推力成正比的控制电流i_v输入到伺服阀传递函数模块，输出一个伺服阀阀芯开口位移值x_v；将伺服阀阀芯开口位移值x_v输入到伺服阀流量计算模块，得到伺服阀流量Q_v输入给伺服阀控制液压推进器动力学系统模块；将恒压变量泵的工作压力p_r输入到伺服阀流量计算模块后输送给伺服阀控制液压推进器动力学系统模块；伺服阀控制液压推进器动力学系统模块根据接收的伺服阀流量Q_v和t-1时刻的马达油液压力p_t-1和角速度ω_t-1，得到马达的油液压力p_t、输出扭矩M_g和角速度ω_t；

将t-1时刻的螺旋桨的进水速度V_t-1和螺旋桨的转速n_t-1输入到螺旋桨进速系数计算模块，得到t时刻螺旋桨的进速系数J_t；将t时刻螺旋桨的进速系数J_t输入到螺旋桨推力与转矩系数计算模块，得到t时刻螺旋桨的转矩系数K_M和推力系数K_T；

将t时刻螺旋桨的推力系数K_T、马达的角速度ω_t输入到螺旋桨推力计算模块得到t时刻螺旋桨推力T_t；将t时刻螺旋桨的转矩系数K_M、马达的转速ω_t输入到螺旋桨转矩计算模块，得t时刻螺旋桨的转矩M_t。

2.一种水下机器人伺服阀控制液压推进器动力学仿真系统的仿真方法，其特征在于：

包括以下步骤：

步骤一：将与期望控制推力成正比的控制电流i_v输入到伺服阀传递函数模块，输出一个伺服阀阀芯开口位移值x_v；

步骤二：根据伺服阀阀芯开口位移值x_v，得到伺服阀流量Q_v：

Q_v＝K_vx_v

其中，K_v——伺服阀的流量增益系数；

步骤三：计算当前t时刻的马达负载流量Q_t：

Q_t＝K_qx_v-K_cp

其中：p——液压马达两端油液压力差；K_q——伺服阀流量增益；K_c——伺服阀流量压力增益；

步骤四：由当前t时刻的马达负载流量Q_t、t-1时刻的马达油液压力p_t-1和角速度ω_t-1，根据液压马达流量连续性方程确定当前t时刻的马达油液压力增加速度

步骤五：由马达油液压力增加速度采样时间Δt和马达油液压力p_t-1，确定当前t时刻马达油液压力p_t：

步骤六：判断马达油液压力p_t是否到达恒压变量泵的工作压力p_r，若达到，则马达油液压力p_t＝p_r，否则马达油液压力p_t＝p_t；

步骤七：根据当前t时刻马达压力p_t，计算t时刻液压马达输出扭矩M_g：

M_g＝D_Mp_t；

步骤八：由t时刻液压马达的输出扭矩M_g、t-1时刻液压马达的负载转矩M_t-1和角速度ω_t-1，得到液压马达的角加速度

其中：M——液压马达的外负载转矩，即为螺旋桨转矩；液压马达的角加速度；J_M——螺旋桨和液压马达轴的转动惯量；M_t——液压马达的负载转矩；

步骤九：由马达角加速度采样时间Δt和t-1时刻的马达角速度ω_t-1，计算t时刻马达角速度ω_t：

步骤十：根据t-1时刻的螺旋桨的进水速度V_t-1和螺旋桨的转速n_t-1，得到t时刻螺旋桨的进速系数J_t：

其中，D——螺旋桨的直径；

步骤十一：根据t时刻螺旋桨的进速系数J_t，通过螺旋桨的敞水特性曲线得到t时刻螺旋桨的转矩系数K_M和推力系数K_T；

步骤十二：根据t时刻螺旋桨的推力系数K_T、螺旋桨的角速度ω_t得到t时刻螺旋桨推力T_t：

T_t＝K_Tρn_t ²D⁴

步骤十三：根据t时刻螺旋桨的转矩系数K_M、螺旋桨的转速n_t，得到t时刻螺旋桨的转矩M_t：

M_t＝K_Mρn_t ²D⁵；

步骤十四：按照p_t-1＝p_t、Q_t-1＝Q_t、M_t-1＝M_t、ω_t-1＝ω_t重新赋初值，返回步骤一。

3.根据权利要求2所述的一种水下机器人伺服阀控制液压推进器动力学仿真系统的仿真方法，其特征在于：伺服阀传递函数模块中的伺服阀传递函数为：

其中，X_v(s)——为伺服阀阀芯位移x_v的拉氏变换；I(s)——为伺服阀输入电流i_v的拉氏变换；K_sv——伺服阀的放大系数；T_v——伺服阀的时间常数。

4.根据权利要求2所述的一种水下机器人伺服阀控制液压推进器动力学仿真系统的仿真方法，其特征在于：所述的液压马达流量连续性方程为：

其中：Q——液压马达负载油液流量；D_M——液压马达排量；V_t——阀腔、马达腔和连接管道的总容积；β_e——等效体积弹性模量；——液压马达的压力增量；C_tc——液压马达的总泄露系数；θ——液压马达角位移；液压马达角速度。