论文部分内容阅读
边界元方法具有高精度和降维等优点,被广泛应用于波浪与结构物相互作用的水动力问题的研究中。但在离散边界积分方程所形成的线性方程组中,系数矩阵为大型非对称的稠密矩阵,方程组的求解耗费大量计算机资源,严重限制了该方法所能处理的问题的规模。探索可用于大型复杂多浮体结构水动力分析的预修正快速傅里叶变换(precorrected Fast Fourier Transform,简称pFFT)边界元方法是本研究的主要目标。预修正快速傅里叶变换(pFFT)方法可以处理任意形式格林函数问题,该方法将面元之间的相互影响分为近场和远场两部分,远场积分的计算采用将面元投影到位于均匀网格的等效奇点来近似计算,近场积分则采用直接积分方法进行计算。事实上,所有的点状奇点都位于均匀网格上。该特性使得由点状奇点所产生的速度势可以通过快速傅里叶变换方法快速求得。结合基于广义最小残量法(GMRES)的迭代方法求解线性方程组,pFFT方法可以将边界元方法的计算复杂度由O(N3)降低至O(N log N),其中N为边界控制点个数。经过科学严谨的算例验证,本研究发展的pFFT边界单元方法完全可以解决大型复杂结构水动力问题。为了能够求解复杂三维结构的水动力问题,本文探索了一种基于有限元和边界元耦合的建模及数据提取方法。针对结构物中常见的模型几何奇异性问题,提出了双重以及多重节点重置法,有效解决了典型高阶边界元方法中加权平均所导致的几何奇异性,进一步提高了水动力问题的计算精度。对于流场问题,采用了可以直接对边界高阶奇异积分方程求解的解析方法,计算结果具有很高的精度。应用Newman的广义运动模态分解方法和拉格朗日乘子法,基于波浪辐射和绕射理论建立了多浮体结构水动力响应的运动方程,应用pFFT方法分别对自由漂浮、刚性连接或铰接约束条件下,由多个浮式、潜式极薄板组成的多浮体结构水动力问题进行了分析,得到了结构的透射系数、反射系数、附加质量、阻尼系数、波浪力、运动响应RAO以及域内速度场分布。数值结果表明该pFFT方法计算精度和计算速度很高,可以用于任意结构形式,任意约束条件下多浮体系统的水动力分析。