原创力文档- 1、本文档共30页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
理想流体力学大作业
学生姓名：
学 号：
2013年10月
Hess—Smith方法计算物体附加质量
作者：
摘 要：本文运用Hess-Smith方法计算了圆球、椭球和圆柱的附加质量系数以及椭球并行的干扰效应。同时，文章分析了网格变化对计算值的影响趋势。本文使用matlab语言对圆球、椭球与圆柱的模型进行了网格有限元的划分，得到各个单元的节点坐标，然后利用Hess-Smith方法对圆球、椭球及并行椭球的附加质量系数进行计算及分析。
关键字：边界元；Hess-Smith; 附加质量系数
一、物理背景
Hess-Smith方法是一种计算任意三维物体势流的方法，该方法由美国的Hess和Smith两人于20 世纪60 年代提出。Hess-Smith方法又称为分布奇点法，作为一种边界元方法，它用许多平面四边形或三角形表面单元来表示物体表面，并在每个单元上布置强度未知的源，然后在物体表面的某些考察点上满足法向速度为零的物面边界条件，得到求单元源密度的线性代数方程组。求解方程组得到源密度分布，进而可求流场内任意点的速度、压力等物理量。
二、理论依据
2.1 分布源模型的建立
为无界流中的物体表面，来流为均匀流，在无穷远处流体的速度为：
(2.1.1)
为定常速度势，并在物体外部空间域中满足拉普拉斯方程，在物面上适合不可穿透条件，在无穷远处，应该与均匀来流的速度势相同。即
（物体外） (2.1.2)
（物面上） (2.1.3)
其中，单位法线向量指向物体内部。
在速度势中分出已知的均匀来流项，记
(2.1.4)
这里的是扰动速度势，应适合以下定解条件：
(2.1.5)
用Rpq表示点p和点q之间的距离，根据格林第三公式，当p点位于物面s外部和远方控制面c的内部之空间域时，有如下公式：
(2.1.6)
由远方边界条件可知，远方封闭控制面上的积分趋于零，从而上式化为：
(2.1.7)
又由式（2.1.5）可得：
(2.1.8)
得到混合分布模型，为了得到单一分布模型表示的扰动势，在物体内部域中构造一个合适的内部解。
于上述物体外部的点P，函数在物体内???域中没有奇点，在物体内部域中对函数和用格林第三公式，得：
(2.1.9)
将（2.1.7）与（2.1.11）相减，得：
(2.1.10)
取下式定解条件中的：
(在内部) (2.1.11)
(在) (2.1.12)
则式(2.1.10)成为：
(2.1.13)
其中，
(2.1.14)
2.2分布源密度的求解
式(2.1.13)中右端分布源的法向导数极限由两部分组成，一部分是p点附近小曲面ε的贡献，另一部分是物面其余部分的贡献。法向指向取向物体内部，小曲面ε的贡献为2πσ（p）,则有如下关系式：
(2.2.1)
再结合物面条件(2.1.5)，得到
(2.2.2)
这就是分布源密度所适合的线性积分方程。
把积分方程(2.2.2)转换成线性代数方程组，即用离散量代替连续变量。把物面分成小块，记
(2.2.3)
用平面四边形或三角形来近似代替小曲面。具体做法如下，取第小块的四个顶点坐标之算术平均值，得到中心点的坐标。计算对角线向量的向量积（指向与曲面法线指向相符合），用表示该方向上的单位向量，形成以为法线且通过中心点的平面，再把四个顶点向该平面作投影，以四个投影点为顶点组成平面四边形，用代替原来的小曲面，称为单元。
通常把小范围内的分布源密度作为常数，因此只要分割不太粗，可以认为在单元上为常数，记作，从而
(2.2.4)
因此物面上的积分可以用个平面四边形（三角形）上积分之和来近似，即
(2.2.5)
上式左端的未知量是连续型变量，而上式右端的未知量是个离散量。为了求解这个未知数，须要个
文档评论（0）