浮体水动力边界元分析中线性方程组迭代求解方法与研究.pdfVIP

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