数值求解方法

由于浮筒受力和各个水动力系数变化的复杂性，无法得到FPTO的解析解，因此采用数值方法求解。采用显示改进欧拉法对运动方程和记忆函数做离散化处理，根据所选参数与初始条件，通过 MATLAB 软件编程，得到运动方程的解，最后根据式(6)得到浮筒在不同系统压力下的平均输出功率，其中规定初始条件为： z( 0)   (0), v( 0)  0 (8)数值模拟中，系统压力选取为 0~ 5MPa 范围内以 0.25MPa 为间距的点。浮筒在每个选定的系统压力上，计算一次平均输出功率，其中计算时间间隔 设定为 80s。最后，将浮筒不同系统压力下的平均输出功率绘制成曲线图，并进行浮筒取能效率的计算，定量与定性地分析浮筒的发电性能。单次数值模拟程序流程见图5。 Fb

v  0  Fr  Mg  k1

F

kv  v