《不规则波中浮体二阶水动力时域数值模拟》由会员分享，可在线阅读，更多相关《不规则波中浮体二阶水动力时域数值模拟（170页珍藏版）》请在金锄头文库上搜索。
1、Classified Index：U66 11 UDC：A Di ssertation for the Degree of DEngrTlo1、l 1meQ0m aln 3lmulatlon 0t 3econd一0rderHydrodyna mic Force on Floating Bodies in Irregular WavesCandidate：Xu GangSupervi sor：ProfDuan Wenyang Academic Degree Applied for：Doctor of EngineeringSpecialty：Fluid Mechanics Date of Sub
2、mission：July，2009Date of Oral Examination：March，2010University：Harbin Engineering University，-II-01L，、，哈尔滨工程大学学位论文原创性声明本人郑重声明：本论文的所有工作，是在导师的指导下， 由作者本人独立完成的。有关观点、方法、数据和文献的引用 已在文中指出，并与参考文献相对应。除文中已注明引用的内 容外，本论文不包含任何其他个人或集体已经公开发表的作品 成果。对本文的研究做出重要贡献的个人和集体，均已在文中 以明确方式标明。本人完全意识到本声明的法律结果由本人承 担。作者(签字)：；枷一 I日
3、期：弘l-年d王月I蝈哈尔滨工程大学学位论文授权使用声明本人完全了解学校保护知识产权的有关规定，即研究生在校 攻读学位期间论文工作的知识产权属于哈尔滨工程大学。哈尔滨 工程大学有权保留并向国家有关部门或机构送交论文的复印件。 本人允许哈尔滨工程大学将论文的部分或全部内容编入有关数据 库进行检索，可采用影印、缩印或扫描等复制手段保存和汇编本 学位论文，可以公布论文的全部内容。同时本人保证毕业后结合 学位论文研究课题再撰写的论文一律注明作者第一署名单位为哈 尔滨工程大学。涉密学位论文待解密后适用本声明。本论文(口在授予学位后即可日在授予学位12个月后口解密后)由哈尔滨工程大学送交有关部门进行保存、
4、汇编等。作者(签字)：i执z，丁，1导师(签字)：f)3(缈-l日期：枷l一，年明ft日如年o：；月7咽。妙一不规则波中浮体二阶水动力时域数值模拟l确k摘要随着我国近海油气资源日趋减少，油气开发逐渐走向深海，发展新型深 海平台技术已经成为国际海洋工程界研究的热点问题，而在时域内能够有效一地模拟非线性波浪力对深海平台的影响是非常重要的。为了满足深海平台设、计中对准确预报载荷和运动的要求，本文将以非线性不规则波对浮体的影响二特性及预报方法进行研究，旨在发展不规则波中深海平台二阶水动力的时域预报技术。本文将研究一种稳定的三维时域数值模拟方法，预报和分析浮体在规则 和不规则波中所受的线性以及非线性的水
5、动力，在远方辐射面上采用多次透 射边界条件(Multitransmitting formula，MTF)，在自由面上采用积分格式的自 由面边界条件(Integral form of Free surface Boundary Condition，IFBC)。为此， 本文着重在以下几个方面开展了研究工作：1本文从时域线性辐射和绕射问题入手，对基于MTF和IFBC边界元法 的参数敏感性进行研究，其中包括附加因子、平滑函数、时间步长、水深、 人工波速以及网格质量等因素的影响；并对时域数值模型的稳定性以及MTF的外传效率进行了研究：长时历的数值模拟，结果不会出现发散现象：同时， 给出在水波领域能够满足
6、数值精度要求的人工波速取值范围；此外还使用不 同形状的人工边界对一般形状浮体的水波问题进行了研究。2将基于MTF和IFBC边界元法应用到不规则波问题的数值模拟中。对 受不规则垂荡激励运动的半球辐射问题和直立圆柱桩的不规则绕射问题进行 了研究，其中包括双色波、多色波以及基于海浪谱的不规则波问题。除在低海况外伪1，3=05m和hl3=10m)，MTF都能有效地将不规则的波动传出人工边 界，在人工边界处达到无反射的效果。由于低海况的波能谱曲线比较平坦，一绝大部分能量所占的频率范围比高海况要宽的多，仅使用MTF法还不能有效。地将所有波动传出人工边界，因此本文采用一种新的水波透射器一阻尼区耦合MTF法(
7、MTF+DZ)作为人工边界条件模拟基于海浪谱的不规则波问题，并 比较了各种阻尼区耦合MTF法的消波效率，得到了最优的MTF+DZ的耦合 形式，同时给出了有效的最短阻尼区长度，此外还使用MTF+DZ法模拟了半潜柱体和半潜式平台的不规则波绕射问题，并使用MATLAB软件对不规则 的时历信号进行快速傅立叶变换，将得到的结果与Hydrostar软件的频域结果进行了比较。3本文对基于MTF和IFBC分布源法的数值精度进行了研究。以直立圆 柱桩绕射为例，将一阶绕射势及其各阶导数与频域解析解进行了比较分析， 得到面元法向突变处分布源法求得的面元切向速度精度较低的结论，并以三 维圆球绕流为例，进一步验证了上述
8、结论。此外本文就现有方法如何提高数值结果的精度进行了研究，得到基于Galerldn技术的分布源法可以提高面元。- t；哈尔滨T程大学博十学位论文切向速度的总体精度，但它对法向突变处面元切向速度的精度提高并不明显。 ：b浦一4本文就基于MTF和IFBC边界元法对二阶水波辐射和绕射问题进行了 研究。使用辅助函数求解二阶非线性自由面条件中与一阶速度势相关的高阶 导数项(包括速度势对Z的二阶导数项、速度势对Z和t的三阶导数项)，避免 使用传统的差分方法，提高了数值精度及稳定性。5本文研究了作强迫垂荡的半潜柱体的辐射问题和直立圆柱的二阶绕射问题，由于所选人工边界条件的消波效率、自由面条件中微分处理方式以
9、 l及物面上高阶导数的求解方法不同，最终导致目前二阶水动力结果之间存在 ，较大离散度。同时，在二阶绕射问题中，本文将一阶速度势及其各阶次导数 的解析解代入二阶非齐次自由面条件，得到的二阶水动力与频域解析解吻合 的非常好，也得到了本文结果与频域半解析解相差较大的主要原因一通过分 布源法得到的控制域形状突变处的面元切向速度精度比较低。此外，本文的 时域二阶结果非常稳定。说明MTF和IFBC可以有效地用于时域二阶问题的 数值模拟，在人工边界上达到无反射的效果。6本文还讨论了二阶自由面条件和物面条件齐次与非齐次对二阶结果 的影响，得到了一些有指导意义的结论。对于辐射问题，二阶速度势产生的 力主要是由非
10、线性的物面条件引起，而非线性的自由面条件对这部分力的贡 献非常小，但它对总的二阶力影响比较大。对于绕射问题，正好与辐射问题 的结论相反，这里二阶速度势产生的力主要是由非线性的自由面条件引起， 而二阶物面条件对这部分力的贡献随着频率的增加逐渐减小，同时，此处的 二阶物面条件对总的二阶力幅值影响也比较大。7本文最后对受不规则激励运动半潜柱体的二阶辐射问题以及其在无 限水深情况下二阶不规则波绕射问题进行了研究，得到的数值结果稳定性非 常高，说明本文方法可以有效用于非线性不规则波问题的时域数值模拟。关键词：二阶水动力；不规则波；人工边界条件；积分格式；MTF；MTF+DZ； IFBC；LNG；半潜式平
11、台不规则波中浮体二阶水动力时域数值模拟ABSTRACTwith the continuous decline of inshore oil and gas resources in Chinamore and more oil and gas resources exploration activities have been moved to the deep seaUnder this circumstanceit has already become all intemational hot topic to一develop new type of deep-water platform
12、s in the field of the ocean engineeringItlpis very important to effcctively simulate nonlinear wave body interaction in timedomainIn order to meet the requirement of predicting the wave loads andmovements on deepwater platforms precisely,the aim of this thesis is to developa predictable method for t
13、ime-domain simulation of second-order hydrodynamicforce on deepwater platforms in irregular wave through investigating thecharacteristics of nonlinear and irregular wave body interactionIn the thesis a stable three-dimensional timedomain methodbased onRankine source，has been applied to simulate line
14、ar and nonlinear hydrodynamics in regular or irregular wavesA MultiTransmitting Formula method(MTF)with an artificial wave speed is employed to satisfy the radiation condition to minimizethe wave reflection on the Artificial Boundary(AB)，and a stable Integration form of Free surface Boundary Condition(IFBC)is used to update velocity potential on the free surfaceMajor works done in this thesis are as follows：Firstly,the semisphere with heave motion and the cylinder with incident wave we
