论文部分内容阅读
随着海洋运输业的不断发展,船舶与海洋工程结构物尺寸越来越大,结构越来越复杂,而且出现了很多双体以及多体问题。基于边界元方法,利用计算机对结构物进行辐射绕射计算过程中线性方程组的求解时间随着物体湿表面面元数目增加而不断增加。对于物面面元数数以万计的大型离岸结构物,传统的直接方法需要花费数天时间求解该线性方程组,这对于以快速立足的边界元方法是不能接受的。为了提高线性方程组求解效率,本文将对其结构特点进行研究并采用合适的迭代方法求解。
线性方程组系数矩阵的结构特性对于迭代求解方法的选择至关重要,本文将从条件数的角度对浮体水动力学边界元分析中产生的线性方程组系数矩阵进行研究,分析其在不同计算模型下性态。间接方法主要是指迭代方法,通过对各种迭代方法以及本领域线性方程组结构特点的研究,在理论上确定适合的迭代求解方法。编写迭代法求解线性方程组求解程序并与直接法求解结果对比验证;为了加速迭代求解方法的收敛,我们采用不完全LU分解预处理技巧降低求解过程中的迭代次数,并通过对不同形状离岸结构进行计算确定预处理技巧的有效性。对于不同结构的线性方程组,迭代方法并不一定快于直接法,通过计算并对比两种方法的浮点运算量获得包含两种求解程序的优化算法。
首先本文简要介绍浮体在波浪中运动的势流理论,然后引出论文所要研究的辐射绕射线性方程组并给出其各项含义。编写矩阵条件数计算程序分析三维浮体无航速频域水动力分析中线性方程组系数矩阵性态,结果表明该矩阵条件数随着入射波频率升高和物体湿表面网格划分数目增加而增大;通过在物体内部自由面划分网格方法消除非规则频率,可以明显缓解非规则频率处条件数突变,但是其它频率处条件数将会普遍增大。对于给定线性方程组,编写基于广义最小残量迭代求解方法程序,结果表明求解效率明显高于基于LU分解的直接求解方法;在上述迭代程序基础上,采用不完全LU分解作为预处理可以显著缓解高频阶段迭代次数升高问题,而且对于含有较大内部自由面的物体,需要加大预处理中不完全LU分解的程度才可以使迭代次数降低到合适的水平;通过对比基于LU分解的直接法和广义最小残量的迭代法浮点运算量,我们获得一个优化判断公式,对于不同的线性方程组,该公式可以在两种方法之间选取运算量最小的方法进行求解。