Naval Architecture and Ocean Engineering
数值求解方法
由于浮筒受力和各个水动力系数变化的复杂性,无法得到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