基于解析复势的翼型、柱体及串列结构附加质量的深度剖析与应用研究

基于解析复势的翼型、柱体及串列结构附加质量的深度剖析与应用研究一、绪论1.1研究背景与意义在流体力学的研究领域中，附加质量是一个关键的概念，它描述了物体在流体中运动时，由于流体的惯性作用，对物体运动产生的一种等效质量影响。当物体在流体中加速或减速时，周围的流体也会随之产生加速度，这种流体的运动变化会对物体施加反作用力，其效果等同于在物体自身质量上增加了一部分额外的质量，这部分额外质量就是附加质量。附加质量的存在深刻影响着物体在流体中的动力学特性，它不仅改变了物体的运动方程，还对物体的稳定性、振动特性以及能量消耗等方面产生重要作用。因此，准确理解和计算附加质量对于深入研究流体与物体的相互作用至关重要。翼型和柱体是流体力学中常见的基本结构，它们广泛应用于航空航天、船舶工程、能源开发等众多领域。在航空航天领域，飞机的机翼采用特定的翼型设计，通过精确控制翼型的形状和参数，实现高效的升力产生和阻力降低，从而确保飞机能够稳定飞行并具备良好的性能。在船舶工程中，船体的形状与柱体结构密切相关，其设计直接影响船舶的航行性能，包括阻力、推进效率以及操纵稳定性等。合理的柱体设计可以减少船舶在水中航行时的阻力，提高燃油经济性，同时增强船舶在复杂海况下的稳定性和操纵性。而在能源开发领域，风力发电机的叶片同样基于翼型原理设计，以最大限度地捕获风能并将其转化为电能，其性能的优劣直接关系到能源的利用效率和发电成本。串列组合结构则是由多个翼型或柱体按照一定的排列方式组合而成的复杂结构形式。这种结构在实际工程中也有着广泛的应用，例如，在多体船设计中，多个船体以串列形式组合，通过优化各船体之间的间距和相对位置，可以有效减小兴波阻力，提高船舶的航行效率和稳定性。在航空领域，一些特殊的飞行器采用串列翼布局，这种布局方式可以改善飞行器的气动性能，提高其机动性和飞行效率。在海洋工程中，串列圆柱结构常用于海洋平台的支撑系统，其复杂的流体动力学特性对于平台的稳定性和安全性至关重要。对翼型、柱体及其串列组合结构的附加质量进行深入研究，具有重要的现实应用价值。在航空航天领域，精确计算翼型的附加质量有助于优化飞机机翼的设计，提高飞机的飞行性能和燃油效率。通过准确把握附加质量对机翼动力学特性的影响，可以更好地预测飞机在各种飞行条件下的响应，从而提高飞行安全性。在船舶工程中，了解柱体和串列组合结构的附加质量对于船舶的设计和航行性能评估至关重要。它可以帮助工程师优化船体结构，降低船舶在水中的阻力，提高推进效率，同时增强船舶在风浪中的稳定性，减少航行风险。在能源开发领域，对于风力发电机叶片和海洋能源采集装置等结构的附加质量研究，可以优化设备的设计，提高能源转换效率，降低能源开发成本。综上所述，研究翼型、柱体及其串列组合结构的附加质量，不仅能够深化我们对流体与物体相互作用的理解，为流体力学理论的发展提供重要支撑，还能为相关工程领域的设计和优化提供关键的理论依据，对于推动航空航天、船舶工程、能源开发等行业的技术进步具有重要意义。1.2国内外研究现状在流体力学领域，关于钝体、儒可夫斯基翼型、串列多圆柱以及儒可夫斯基翼型与圆柱组合的附加质量研究一直是重要的课题，国内外学者在这些方面取得了丰硕的成果。在钝体附加质量的研究方面，早期的研究主要集中在理论分析上。学者们基于势流理论，对简单形状的钝体，如圆球、圆柱等，进行了深入的理论推导，试图精确求解其附加质量。例如，对于圆球在无界理想流体中的附加质量，通过经典的势流理论，可以得到精确的解析解，其附加质量与圆球的体积和流体密度密切相关。然而，对于形状更为复杂的钝体，理论求解面临巨大挑战。随着计算机技术的飞速发展，数值计算方法逐渐成为研究复杂钝体附加质量的重要手段。CFD技术，如Fluent、CFX等软件，被广泛应用于模拟钝体周围的流场，通过数值计算得到物体所受的流体作用力，进而推算出附加质量。在一些研究中，利用Fluent软件的动网格技术，模拟钝体在流体中的变速运动，通过监测物体所受的合力，结合动力学方程，成功计算出了复杂钝体的附加质量。实验研究也是不可或缺的一部分。学者们通过在实验室中构建各种钝体模型，利用高精度的测量设备，如力传感器、粒子图像测速仪（PIV）等，测量钝体在流体中运动时所受的力和周围流场的速度分布，从而直接或间接地获取附加质量数据。尽管取得了这些进展，但目前对于钝体在复杂流场条件下，如强湍流、多相流等环境中的附加质量研究仍显不足，相关的理论模型和计算方法有待进一步完善。儒可夫斯基翼型作为一种经典的翼型，其附加质量的研究也受到了广泛关注。早期的研究主要围绕其在理想流体中的绕流问题展开，通过保角变换等数学方法，将复杂的翼型绕流问题转化为简单的圆柱绕流问题，从而求解出翼型的速度势和流函数，进而得到附加质量的解析表达式。这些解析结果为理解儒可夫斯基翼型的气动特性提供了重要的理论基础。随着研究的深入，考虑粘性影响的数值模拟方法逐渐兴起。通过求解Navier-Stokes方程，利用有限元、有限体积等数值方法，能够更准确地模拟翼型在实际粘性流体中的流动情况，得到更为精确的附加质量结果。研究还发现，翼型的几何参数，如厚度、弯度等，对附加质量有着显著的影响。增加翼型的厚度会使附加质量增大，因为厚度的增加导致翼型周围流体的扰动增强，流体惯性作用增大。然而，目前对于儒可夫斯基翼型在非定常流动，如动态失速、颤振等工况下的附加质量研究还不够深入，相关的研究成果相对较少。串列多圆柱结构在工程中有着广泛的应用，如海洋立管、换热器管束等，其附加质量的研究具有重要的工程意义。国内外学者在这方面开展了大量的研究工作。在实验研究方面，通过在风洞或水洞中进行串列圆柱的模型实验，测量不同间距、不同雷诺数下圆柱所受的力和周围流场的特性，从而分析串列圆柱的附加质量变化规律。实验结果表明，串列圆柱的附加质量不仅与圆柱的直径、间距有关，还受到雷诺数的影响。当圆柱间距较小时，上游圆柱的尾流会对下游圆柱产生显著的干扰，导致下游圆柱的附加质量发生明显变化。在数值模拟方面，采用CFD方法对串列圆柱绕流进行模拟，能够详细分析流场的细节，如涡脱落、尾流相互作用等对附加质量的影响。一些研究利用大涡模拟（LES）方法，能够更准确地捕捉到流场中的湍流结构，为串列圆柱附加质量的研究提供了更丰富的信息。然而，由于串列多圆柱绕流的复杂性，目前的研究成果在一些关键问题上尚未达成完全一致，如涡脱落的频率和模式对附加质量的具体影响机制等，仍有待进一步深入研究。对于儒可夫斯基翼型与圆柱组合的附加质量研究相对较少，但近年来也逐渐受到关注。这种组合结构在一些特殊的工程应用中具有独特的优势，如在某些新型飞行器的设计中，采用儒可夫斯基翼型与圆柱的组合布局，以实现特定的气动性能。目前的研究主要集中在数值模拟和理论分析方面。通过数值模拟方法，如边界元法、有限元法等，求解组合结构周围的流场，得到附加质量的数值结果。在理论分析方面，尝试将儒可夫斯基翼型和圆柱的绕流理论相结合，建立适用于该组合结构的附加质量计算模型。然而，由于组合结构的几何形状和流场特性更为复杂，目前的研究还处于探索阶段，相关的研究成果还不够完善，需要进一步开展深入的研究工作，以揭示其附加质量的变化规律和影响因素。国内外在钝体、儒可夫斯基翼型、串列多圆柱以及儒可夫斯基翼型与圆柱组合的附加质量研究方面取得了一定的进展，但仍存在许多有待解决的问题和研究空白。在未来的研究中，需要进一步加强理论分析、数值模拟和实验研究的有机结合，深入探究各种因素对附加质量的影响机制，为相关工程领域的设计和优化提供更加准确和可靠的理论依据。1.3研究内容与方法本文主要围绕翼型、柱体及其串列组合结构的附加质量展开深入研究，综合运用理论分析与数值计算相结合的方法，全面剖析其附加质量的特性与影响因素。在理论分析方面，以解析复势理论为核心，深入推导翼型、柱体及其串列组合结构的附加质量公式。对于圆柱，基于其在流体中运动的基本原理，利用解析复势构建其流场的数学模型，通过对速度势和流函数的精确求解，推导得出圆柱附加质量的理论表达式。对于椭圆柱，考虑其特殊的几何形状，运用保角变换等数学手段，将复杂的椭圆柱绕流问题转化为便于处理的形式，进而推导出椭圆柱的解析复势，在此基础上得出其附加质量的计算公式。对于儒可夫斯基翼型，依据儒可夫斯基变换，将翼型绕流与圆柱绕流建立联系，通过对圆柱绕流问题的深入分析，得到儒可夫斯基翼型的解析复势和附加质量公式。在推导串列双圆柱和儒可夫斯基翼型与圆柱组合结构的附加质量时，充分考虑各结构之间的相互干扰效应，运用Milne-Thomson圆定理等理论，结合叠加原理，构建组合结构的流场模型，从而推导出相应的附加质量公式。在数值计算方面，借助专业的CFD软件Fluent进行模拟分析。首先，对翼型、柱体及其串列组合结构进行精确的几何建模，确保模型能够准确反映实际结构的形状和尺寸。然后，依据实际工况设置合理的边界条件，如入口流速、出口压力、壁面条件等，同时选择合适的湍流模型，以准确模拟流体的流动特性。在计算过程中，利用Fluent的动网格技术，模拟物体在流体中的运动过程，通过监测物体所受的流体作用力，结合动力学方程，计算得到附加质量的数值结果。将数值计算结果与理论推导结果进行详细对比，分析两者之间的差异，验证理论推导的准确性，并深入探讨数值计算中可能存在的误差来源。通过理论推导与数值计算相结合的方式，本文全面研究了翼型、柱体及其串列组合结构的附加质量。分析了几何参数，如圆柱半径、椭圆柱的长轴和短轴、儒可夫斯基翼型的厚度和弯度等，以及结构参数，如串列双圆柱的间隔距离、儒可夫斯基翼型与圆柱的相对位置等，对附加质量的影响规律。同时，探讨了不同流场条件，如雷诺数、来流速度等，对附加质量的作用机制，为相关工程领域的设计和优化提供了全面而准确的理论依据和数据支持。二、附加质量理论基础2.1现有解析分析方法概述2.1.1势函数解析法势函数解析法作为附加质量计算中的经典方法，其基本原理根植于势流理论。在理想流体假设下，即流体无粘性、不可压缩且流动无旋，流场的运动可以通过一个标量函数——速度势函数来描述。对于二维流动，速度分量u和v可以表示为速度势函数\varphi对坐标的偏导数，即u=\frac{\partial\varphi}{\partialx}，v=\frac{\partial\varphi}{\partialy}。在附加质量计算中，势函数解析法的应用通常围绕着求解物体绕流问题展开。对于简单几何形状的物体，如圆柱、圆球等，通过建立合适的坐标系，并利用边界条件，可以推导出精确的速度势函数表达式。以圆柱绕流为例，在极坐标系下，均匀来流绕圆柱的速度势函数可以表示为\varphi=U_{\infty}r\cos\theta(1-\frac{a^2}{r^2})，其中U_{\infty}为来流速度，r为极径，\theta为极角，a为圆柱半径。基于此速度势函数，通过对流体动能的计算，再结合附加质量的定义，即单位加速度下流体对物体施加的力，就可以得到圆柱的附加质量。对于圆柱在二维平面运动中，沿x方向和y方向的附加质量m_{xx}和m_{yy}相等，且m_{xx}=m_{yy}=\rho\pia^2，其中\rho为流体密度。然而，势函数解析法存在明显的局限性。当面对复杂几何形状的物体时，精确求解速度势函数变得极为困难，甚至无法实现。对于具有不规则外形的翼型，由于其边界条件的复杂性，很难找到合适的数学表达式来描述速度势函数。即使对于一些相对规则但稍复杂的形状，如椭圆柱，虽然可以通过保角变换等数学手段将其绕流问题转化为相对简单的形式来求解速度势，但过程繁琐且需要深厚的数学基础。势函数解析法通常基于理想流体假设，忽略了流体的粘性效应。在实际工程应用中，粘性对物体绕流和附加质量的影响往往不可忽视，特别是在低雷诺数或边界层内的流动，这使得势函数解析法的计算结果与实际情况存在一定偏差。2.1.2流线变形解析法流线变形解析法的核心思想是通过分析物体在流体中运动时流线的变形情况来确定附加质量。当物体在流体中运动时，会引起周围流体的流动，流线会发生弯曲和变形。流线的变形程度与物体的运动状态、几何形状以及流体的性质密切相关。该方法认为，附加质量的产生源于流体为了适应物体的运动而改变自身的流动状态，这种改变所需要的能量等价于在物体上附加了一定的质量。在具体应用中，流线变形解析法通常借助于一些数学模型和物理假设来描述流线的变形。对于简单的流动情况，可以通过理论分析建立流线方程，进而分析流线的变形特征。在均匀来流绕平板的流动中，可以基于势流理论和边界层理论，推导出平板附近流线的方程，通过观察流线在平板前后的弯曲程度来估算附加质量。对于更复杂的流动，如串列多圆柱绕流，由于圆柱之间的相互干扰，流线的变形更加复杂，需要采用数值模拟或实验测量的方法来获取流线信息。然而，该方法在处理复杂流场时面临诸多难点和挑战。复杂流场中存在着各种复杂的流动现象，如边界层分离、涡旋的生成与脱落等，这些现象会导致流线的形态变得极为复杂，难以用简单的数学模型进行准确描述。在串列双圆柱绕流中，上游圆柱的尾流会对下游圆柱产生强烈的干扰，导致下游圆柱周围的流线出现剧烈的扭曲和分离，准确捕捉这些流线的变化并将其转化为附加质量的计算变得非常困难。流场的非线性特性也是一个重要问题。随着流速的增加或物体运动状态的变化，流场往往会呈现出非线性特征，此时流线的变形不再是简单的线性关系，传统的基于线性假设的流线变形解析方法难以适用。而且，该方法对实验测量和数值模拟的依赖程度较高，实验测量过程中存在测量误差，数值模拟也面临着计算精度和计算效率的问题，这都会影响到基于流线变形解析法得到的附加质量结果的准确性和可靠性。2.2基于解析复势的研究方法2.2.1动能及附加质量理论推导在势流理论的框架下，流体被假定为理想流体，即无粘性、不可压缩且流动无旋。对于二维不可压缩无旋流动，可引入复势函数W(z)=\varphi(x,y)+i\psi(x,y)，其中\varphi(x,y)为速度势函数，\psi(x,y)为流函数，z=x+iy为复变量。速度分量u和v可通过复势函数表示为u=\frac{\partial\varphi}{\partialx}=\frac{\partial\psi}{\partialy}，v=\frac{\partial\varphi}{\partialy}=-\frac{\partial\psi}{\partialx}。流体的动能T是研究附加质量的关键物理量。对于单位长度的二维流动，其动能可表示为：T=\frac{1}{2}\rho\iint_{S}(u^{2}+v^{2})dxdy其中\rho为流体密度，S为流场区域。将速度分量用复势函数表示并代入动能表达式，通过格林公式等数学变换，可将动能积分转化为沿流场边界C的线积分：T=\frac{1}{2}\rho\oint_{C}\psi\frac{\partial\varphi}{\partialn}ds其中\frac{\partial\varphi}{\partialn}为速度势函数沿边界的法向导数，ds为边界微元弧长。当物体在流体中作匀速直线运动时，设物体的速度为U，方向沿x轴正方向。此时，复势函数可表示为W(z)=Uf(z)，其中f(z)是与物体形状和位置相关的解析函数。将复势函数代入动能表达式，经过一系列的数学推导（包括对解析函数的性质运用、积分运算等），可得到动能与物体速度U的关系为T=\frac{1}{2}m^{*}U^{2}，其中m^{*}即为附加质量。以圆柱在均匀来流中的运动为例，均匀来流速度为U_{\infty}，圆柱半径为a。其复势函数为W(z)=U_{\infty}(z+\frac{a^{2}}{z})，通过上述动能与附加质量的推导过程，可求得圆柱在二维平面运动中，沿x方向和y方向的附加质量m_{xx}和m_{yy}相等，且m_{xx}=m_{yy}=\rho\pia^{2}。这一结果表明，圆柱的附加质量与流体密度和圆柱的横截面积密切相关，体现了流体惯性对圆柱运动的影响。对于其他形状的物体，如椭圆柱、儒可夫斯基翼型等，同样可基于其特定的解析复势，按照上述推导思路，得出相应的附加质量表达式。2.2.2Fluent数值计算原理Fluent软件是计算流体动力学（CFD）领域中广泛应用的一款商业软件，它基于有限体积法对控制方程进行离散求解，能够精确模拟各种复杂的流体流动现象，在附加质量计算中发挥着重要作用。在Fluent中，控制方程主要包括质量守恒方程、动量守恒方程和能量守恒方程。对于不可压缩流体，质量守恒方程（连续性方程）可表示为：\frac{\partialu_{i}}{\partialx_{i}}=0其中u_{i}为速度分量，x_{i}为坐标分量。动量守恒方程（Navier-Stokes方程）在笛卡尔坐标系下的形式为：\rho(\frac{\partialu_{i}}{\partialt}+u_{j}\frac{\partialu_{i}}{\partialx_{j}})=-\frac{\partialp}{\partialx_{i}}+\mu(\frac{\partial^{2}u_{i}}{\partialx_{j}\partialx_{j}})+F_{i}其中\rho为流体密度，p为压力，\mu为动力粘度，F_{i}为体积力分量。在处理实际工程问题时，还需根据具体情况考虑能量守恒方程以及其他相关的物理模型，如湍流模型、多相流模型等。在数值计算过程中，首先需要对计算域进行网格划分，将连续的计算域离散化为一系列的控制体积。Fluent支持多种网格类型，包括结构化网格和非结构化网格。结构化网格具有规则的拓扑结构，计算效率高，适用于几何形状简单的模型；非结构化网格则能更好地适应复杂的几何形状，通过灵活的网格布局，准确捕捉物体表面和流场中的关键特征。在对翼型、柱体及其串列组合结构进行数值模拟时，根据模型的几何复杂度，合理选择网格类型和划分策略，以确保计算精度和效率。对于简单的圆柱模型，可采用结构化的圆形网格，保证网格在圆柱表面的正交性，提高计算精度；而对于形状复杂的儒可夫斯基翼型，非结构化的三角形或四面体网格则能更好地贴合翼型表面，准确描述其几何特征。网格划分完成后，将控制方程在每个控制体积上进行积分，得到离散化的方程组。Fluent采用有限体积法对控制方程进行离散，该方法的基本思想是将控制方程在每个控制体积上进行积分，将偏微分方程转化为代数方程组。在离散过程中，通过对速度、压力等物理量在控制体积界面上的插值和近似处理，构建离散化的方程组。对于对流项的离散，可采用一阶迎风差分格式、二阶迎风差分格式或更高级的QUICK格式等，不同的格式具有不同的精度和计算稳定性；对于扩散项的离散，通常采用中心差分格式，以保证计算的准确性。离散化后的方程组通过迭代求解的方式得到数值解。Fluent提供了多种求解器，如压力-速度耦合算法（SIMPLE算法、SIMPLEC算法、PISO算法等），用于求解压力和速度的耦合关系。在求解过程中，通过不断迭代更新速度和压力场，直至满足收敛条件。收敛条件通常根据计算精度要求设置，如残差收敛标准、物理量的变化率等。当计算结果满足收敛条件时，认为得到了稳定的数值解。在附加质量计算中，Fluent利用动网格技术模拟物体在流体中的运动。动网格技术允许网格随着物体的运动而变形，从而准确捕捉物体与流体之间的相对运动。通过设置物体的运动边界条件，如速度、加速度等，模拟物体在流体中的各种运动状态。在模拟串列双圆柱的运动时，可分别设置两个圆柱的运动速度和方向，考虑它们之间的相互干扰对附加质量的影响。在计算过程中，Fluent实时监测物体所受的流体作用力，包括压力力和粘性力。根据牛顿第二定律，物体所受的合力等于其质量与加速度的乘积，通过测量物体在运动过程中所受的合力，并结合物体的加速度，利用公式m^{*}=\frac{F}{a}（其中F为物体所受的合力，a为物体的加速度），即可计算得到附加质量。2.3坐标转换与相关定理2.3.1保角变换保角变换，又称为共形映射，是复变函数论中的一个重要概念，在流体力学中有着广泛的应用。从数学定义上讲，若函数w=f(z)在区域D内解析，且对于D内的任意一点z，其导数f'(z)\neq0，则称函数w=f(z)在区域D内是保角的。这意味着在z平面上的两条曲线在交点处的夹角，经过变换w=f(z)后，在w平面上对应曲线的夹角大小和方向都保持不变。保角变换不仅保持角度不变，还具有一些其他重要性质。它能够保持曲线的光滑性，若z平面上的曲线是光滑的，那么经过保角变换后，在w平面上对应的曲线依然光滑。保角变换还具有保圆性，即z平面上的圆或直线，在w平面上对应的图形要么是圆，要么是直线。在流体力学中，保角变换主要用于简化复杂流场的几何形状，从而方便附加质量的计算。对于具有复杂边界形状的流场，直接求解其控制方程往往非常困难，甚至无法得到解析解。通过保角变换，可以将z平面上复杂的边界形状映射到w平面上相对简单的形状，如圆形、上半平面或带形域等。在研究具有不规则边界的翼型绕流问题时，可以利用保角变换将翼型的边界映射为圆形边界，这样就可以将复杂的翼型绕流问题转化为相对简单的圆柱绕流问题。对于圆柱绕流问题，基于势流理论，已经有较为成熟的解析求解方法，可以方便地得到流场的速度势、流函数以及附加质量等物理量的解析表达式。然后，再通过逆变换，将w平面上的解转换回z平面，从而得到原复杂流场的解。保角变换还可以用于处理多个物体的绕流问题，通过合理选择变换函数，将多个物体的边界映射为简单的形状，利用叠加原理求解流场，进而计算附加质量。2.3.2儒可夫斯基变换儒可夫斯基变换是一种特殊的保角变换，在翼型和柱体附加质量分析中具有重要的应用，其数学表达式为z=\zeta+\frac{b^{2}}{\zeta}，其中z=x+iy是物理平面（z平面）上的复变量，\zeta=\xi+i\eta是变换平面（\zeta平面）上的复变量，b为常数。从几何意义上看，儒可夫斯基变换具有独特的映射特性。在\zeta平面上，以原点为圆心、半径为b的圆，经过儒可夫斯基变换后，在z平面上对应的是一个平板。当\zeta平面上的圆半径R>b时，经过儒可夫斯基变换，在z平面上得到的是一个具有一定厚度和弯度的儒可夫斯基翼型。通过调整\zeta平面上圆的参数，如圆心位置、半径大小等，可以得到不同形状的儒可夫斯基翼型，这使得儒可夫斯基变换成为设计和分析翼型的有力工具。在翼型附加质量分析中，儒可夫斯基变换的应用主要基于将复杂的翼型绕流问题转化为相对简单的圆柱绕流问题。对于儒可夫斯基翼型绕流，首先通过儒可夫斯基变换将儒可夫斯基翼型的边界映射为\zeta平面上的圆形边界。在\zeta平面上，对于均匀来流绕圆柱的流动，其复势函数可以根据势流理论得到精确的表达式。对于均匀来流速度为U_{\infty}，绕半径为R的圆柱的流动，其复势函数为W(\zeta)=U_{\infty}(\zeta+\frac{R^{2}}{\zeta})。然后，利用儒可夫斯基变换的逆变换，将\zeta平面上的复势函数转换回z平面，得到儒可夫斯基翼型绕流的复势函数。基于复势函数，可以进一步计算流场的速度分布、压力分布以及附加质量等物理量。在计算附加质量时，根据前面所述的基于解析复势的附加质量推导方法，通过对复势函数的分析和运算，得到儒可夫斯基翼型的附加质量表达式。对于柱体附加质量分析，当柱体的形状可以通过儒可夫斯基变换与圆形建立联系时，同样可以利用该变换简化分析过程，得到柱体的附加质量。2.3.3圆到椭圆的坐标变换圆到椭圆的坐标变换是一种常用的坐标转换方式，在流体力学中，当需要将圆形流场转换为椭圆形流场进行附加质量分析时，该变换起着关键作用。在极坐标系下，圆的方程可以表示为r=a（a为圆的半径），而椭圆的方程可以表示为\frac{x^{2}}{A^{2}}+\frac{y^{2}}{B^{2}}=1（A和B分别为椭圆的长半轴和短半轴）。为了实现从圆到椭圆的坐标转换，引入变换x=A\frac{r}{a}\cos\theta，y=B\frac{r}{a}\sin\theta。其中r和\theta是极坐标变量，x和y是直角坐标变量。当r=a时，x=A\cos\theta，y=B\sin\theta，这正是椭圆的参数方程，表明该变换可以将半径为a的圆映射为长半轴为A、短半轴为B的椭圆。在将圆形流场转换为椭圆形流场进行附加质量分析时，具体步骤如下：首先，对于圆形流场，基于势流理论建立其复势函数。对于均匀来流绕半径为a的圆柱的流动，复势函数为W(z)=U_{\infty}(z+\frac{a^{2}}{z})（z=x+iy）。然后，将坐标变换x=A\frac{r}{a}\cos\theta，y=B\frac{r}{a}\sin\theta代入复势函数中，实现复势函数从圆形流场坐标到椭圆形流场坐标的转换。在转换后的复势函数基础上，利用基于解析复势的附加质量推导方法，通过对复势函数的求导、积分等运算，得到椭圆形流场中物体的附加质量表达式。在进行坐标变换和附加质量计算过程中，需要注意以下要点：一是要确保坐标变换的准确性，特别是在代入复势函数和进行变量替换时，要严格按照变换公式进行操作，避免出现错误。二是要考虑边界条件的转换，圆形流场和椭圆形流场的边界条件不同，在进行坐标变换后，需要根据新的边界条件对复势函数进行修正和调整，以保证计算结果的正确性。2.3.4Milne-Thomson圆定理Milne-Thomson圆定理是流体力学中用于解决圆柱绕流等问题的重要定理，在附加质量计算中有着广泛的应用。该定理的内容为：若W(z)是解析函数，且在z平面上除圆|z-z_{0}|=a外的区域内解析，那么函数W_{1}(z)=W(z)+\overline{W(\frac{a^{2}}{z-z_{0}}+z_{0})}在z平面上除圆|z-z_{0}|\leqa外的区域内解析，且在圆|z-z_{0}|=a上满足W_{1}(z)的虚部为常数，即\text{Im}[W_{1}(z)]|_{|z-z_{0}|=a}=C（C为常数）。这意味着W_{1}(z)所代表的流场在圆|z-z_{0}|=a上满足物面条件，可用于描述绕该圆柱的流动。在解决圆柱绕流问题时，Milne-Thomson圆定理的应用十分广泛。对于均匀来流绕圆柱的流动，假设均匀来流的复势函数为W_{0}(z)=U_{\infty}z（U_{\infty}为来流速度），根据Milne-Thomson圆定理，绕半径为a，圆心在z_{0}处的圆柱的复势函数可以表示为W(z)=U_{\infty}z+\overline{U_{\infty}(\frac{a^{2}}{z-z_{0}}+z_{0})}。通过对该复势函数的分析，可以得到圆柱绕流的速度分布、压力分布等信息。在计算附加质量时，利用前面推导的基于解析复势的附加质量公式，将绕圆柱的复势函数代入其中，通过对复势函数的相关运算，得到圆柱在流体中的附加质量。在实际应用中，对于多个圆柱的串列组合结构，也可以利用Milne-Thomson圆定理结合叠加原理来求解流场和附加质量。对于串列双圆柱结构，先分别考虑每个圆柱对流场的影响，利用Milne-Thomson圆定理得到每个圆柱单独存在时的复势函数，然后根据叠加原理，将两个圆柱的复势函数相加，得到串列双圆柱的复势函数，进而计算其附加质量。2.4本章小结本章深入探讨了附加质量的理论基础，系统阐述了现有解析分析方法，包括势函数解析法和流线变形解析法，详细剖析了它们的基本原理、应用方式以及各自存在的局限性。重点介绍了基于解析复势的研究方法，通过严谨的理论推导，明确了动能与附加质量之间的紧密关系，并详细阐述了Fluent数值计算的原理，为后续的研究提供了坚实的理论和方法支撑。深入研究了坐标转换与相关定理，如保角变换、儒可夫斯基变换、圆到椭圆的坐标变换以及Milne-Thomson圆定理等。保角变换通过将复杂流场边界映射为简单形状，有效简化了流场分析；儒可夫斯基变换则为翼型和柱体附加质量分析提供了关键的转换方式，实现了复杂翼型绕流问题向简单圆柱绕流问题的转化；圆到椭圆的坐标变换为椭圆形流场中物体附加质量的计算提供了有效的途径；Milne-Thomson圆定理在解决圆柱绕流及串列组合结构问题时发挥了重要作用，通过该定理可以方便地得到圆柱绕流的复势函数，进而计算附加质量。这些坐标转换与相关定理在后续对翼型、柱体及其串列组合结构附加质量的深入分析中，将起到不可或缺的作用，是构建理论模型和进行数值计算的重要工具。三、单体附加质量分析3.1圆柱附加质量分析在理想流体中，当圆柱作匀速直线运动时，其周围流场的解析复势是研究附加质量的关键。假设圆柱半径为a，以速度U沿x轴正方向在均匀来流中运动，均匀来流速度为U_{\infty}。根据势流理论，圆柱绕流的复势函数可表示为W(z)=U_{\infty}(z+\frac{a^{2}}{z})，其中z=x+iy为复变量。从理论推导附加质量的过程如下：首先，根据复势函数求速度分量。速度分量u和v可通过复势函数对z的导数得到，u-iv=\frac{dW}{dz}，对W(z)=U_{\infty}(z+\frac{a^{2}}{z})求导可得u-iv=U_{\infty}(1-\frac{a^{2}}{z^{2}})，将z=x+iy代入，经过复数运算可得到u和v在直角坐标系下的表达式。然后，计算流体的动能T。对于单位长度的二维流动，动能T=\frac{1}{2}\rho\iint_{S}(u^{2}+v^{2})dxdy，通过将速度分量表达式代入，并利用格林公式将面积分转化为沿流场边界的线积分，经过一系列复杂的数学推导（包括对复变函数积分性质的运用、边界条件的处理等），最终得到动能T=\frac{1}{2}\rho\pia^{2}U^{2}。根据附加质量的定义，当物体在流体中运动时，动能可表示为T=\frac{1}{2}m^{*}U^{2}，其中m^{*}为附加质量。对比可得圆柱在二维平面运动中，沿x方向和y方向的附加质量m_{xx}和m_{yy}相等，且m_{xx}=m_{yy}=\rho\pia^{2}。为了深入分析圆柱半径、运动速度等因素对附加质量的影响规律，通过数值模拟和理论分析相结合的方式进行研究。当圆柱半径a发生变化时，附加质量m^{*}=\rho\pia^{2}，可以明显看出附加质量与半径的平方成正比关系。随着半径的增大，附加质量会迅速增加，这是因为半径的增大导致圆柱周围扰动的流体体积增大，流体的惯性作用增强，从而使附加质量增大。在实际工程中，对于海洋平台的支撑圆柱结构，若增大圆柱半径，其在海流中的附加质量将显著增加，这对平台的稳定性和动力学响应会产生重要影响。当圆柱运动速度U改变时，从附加质量的计算公式m_{xx}=m_{yy}=\rho\pia^{2}来看，速度并不直接影响附加质量的大小。但在实际流场中，速度的变化会引起流场的变化，如雷诺数Re=\frac{\rhoUD}{\mu}（D为特征长度，对于圆柱D=2a，\mu为动力粘度）会随着速度的变化而改变。当速度增大，雷诺数增大，流场可能从层流转变为湍流，流场的复杂性增加，虽然附加质量的理论值不变，但流体对圆柱的作用力特性会发生改变，进而影响圆柱的运动状态和动力学响应。3.2椭圆柱附加质量分析3.2.1椭圆柱的解析复势在理想流体中，推导椭圆柱的解析复势是研究其附加质量的基础。通过保角变换，将z平面上的椭圆柱绕流问题转化为\zeta平面上的圆柱绕流问题。引入儒可夫斯基变换的一般形式z=\zeta+\frac{b^{2}}{\zeta}，其中z=x+iy是物理平面（z平面）上的复变量，\zeta=\xi+i\eta是变换平面（\zeta平面）上的复变量，b为常数。对于长半轴为A、短半轴为B的椭圆柱，设\zeta平面上与之对应的圆柱半径为a，圆心位于原点。根据儒可夫斯基变换的性质，当\zeta平面上的圆经过变换后，在z平面上形成椭圆柱，此时有关系A=a+\frac{b^{2}}{a}，B=a-\frac{b^{2}}{a}。联立这两个方程，求解可得a=\frac{A+B}{2}，b=\sqrt{AB}。在\zeta平面上，均匀来流绕半径为a的圆柱的复势函数为W(\zeta)=U_{\infty}(\zeta+\frac{a^{2}}{\zeta})，其中U_{\infty}为来流速度。将z=\zeta+\frac{b^{2}}{\zeta}代入上式，通过变量代换和数学运算，得到z平面上椭圆柱绕流的解析复势W(z)。在推导过程中，假设流体为理想流体，即无粘性、不可压缩且流动无旋，同时忽略了流体的重力和其他外力的影响，仅考虑流体与椭圆柱之间的惯性相互作用。通过这种简化假设，能够利用势流理论和复变函数的方法，较为简洁地推导得到椭圆柱的解析复势，为后续附加质量的计算奠定基础。3.2.2椭圆柱的附加质量计算在得到椭圆柱的解析复势后，利用前面推导的动能与附加质量关系公式来计算椭圆柱的附加质量。首先，根据解析复势W(z)求速度分量。速度分量u和v可通过复势函数对z的导数得到，即u-iv=\frac{dW}{dz}。对W(z)进行求导运算，由于W(z)是通过复杂的变换得到的，求导过程需要运用复合函数求导法则和复变函数的求导公式。经过一系列的求导和化简，得到u和v在直角坐标系下的表达式。然后，计算流体的动能T。对于单位长度的二维流动，动能T=\frac{1}{2}\rho\iint_{S}(u^{2}+v^{2})dxdy，将速度分量表达式代入动能积分式中。通过将积分区域从z平面转换到\zeta平面，利用z与\zeta之间的变换关系以及相应的雅克比行列式，将积分转化为在\zeta平面上的积分。经过复杂的积分运算（包括对复变函数积分性质的运用、积分区间的确定等），最终得到动能T的表达式。根据附加质量的定义，当物体在流体中运动时，动能可表示为T=\frac{1}{2}m^{*}U^{2}，其中m^{*}为附加质量，U为物体的运动速度。对比动能T的表达式和\frac{1}{2}m^{*}U^{2}，可以得到椭圆柱的附加质量m^{*}的表达式。对于沿x方向和y方向运动的椭圆柱，其附加质量m_{xx}和m_{yy}的表达式分别为：m_{xx}=\rho\pi\frac{(A+B)^{2}}{4}m_{yy}=\rho\pi\frac{(A-B)^{2}}{4}其中\rho为流体密度，A和B分别为椭圆柱的长半轴和短半轴。这两个表达式清晰地展示了椭圆柱的附加质量与流体密度、长半轴和短半轴之间的关系。3.2.3附加质量结果分析分析椭圆柱长轴、短轴长度变化对附加质量的影响，对于理解椭圆柱在流体中的动力学特性至关重要。当椭圆柱长轴（长半轴A）增长时，从附加质量的表达式m_{xx}=\rho\pi\frac{(A+B)^{2}}{4}和m_{yy}=\rho\pi\frac{(A-B)^{2}}{4}可以看出，m_{xx}和m_{yy}都会增大。这是因为长轴的增长使得椭圆柱周围扰动的流体体积增大，流体的惯性作用增强，从而导致附加质量增加。在实际工程中，如海洋工程中的柱状结构，若其形状接近椭圆柱且长轴增长，在海流作用下，其附加质量增大，对结构的稳定性和受力情况会产生显著影响，可能需要加强结构设计以承受更大的流体作用力。当短轴（短半轴B）增长时，m_{xx}同样增大，而m_{yy}的变化则需要具体分析。当A\gtB时，随着B的增大，m_{yy}先减小后增大，存在一个使m_{yy}最小的短轴长度。这是由于短轴的变化会改变椭圆柱的形状，进而影响流体的绕流模式和惯性作用分布。研究椭圆柱运动方向与附加质量之间的关系，发现椭圆柱在不同方向上的附加质量存在差异。沿长轴方向（x方向）运动时，附加质量为m_{xx}=\rho\pi\frac{(A+B)^{2}}{4}；沿短轴方向（y方向）运动时，附加质量为m_{yy}=\rho\pi\frac{(A-B)^{2}}{4}。由于(A+B)^{2}\gt(A-B)^{2}（A\gt0，B\gt0），所以m_{xx}\gtm_{yy}，即椭圆柱沿长轴方向运动时的附加质量大于沿短轴方向运动时的附加质量。这是因为沿长轴方向运动时，椭圆柱对流体的扰动范围更大，引起的流体惯性作用更强。这种运动方向与附加质量的关系在实际应用中具有重要意义，例如在飞行器设计中，若飞行器的某些部件形状类似椭圆柱，其运动方向的选择需要考虑附加质量的影响，以优化飞行器的动力学性能和操控稳定性。3.2.4Fluent数值计算与验证为了验证理论分析的正确性，建立椭圆柱在流体中运动的数值模型，利用Fluent软件进行数值计算。首先，在建模软件中创建精确的椭圆柱几何模型，确保长半轴A和短半轴B的尺寸准确无误。将几何模型导入Fluent软件后，进行网格划分。考虑到椭圆柱表面的曲率变化，采用非结构化网格进行划分，在椭圆柱表面和近壁区域加密网格，以提高计算精度，准确捕捉流场的变化。在圆柱表面，网格尺寸设置为较小的值，如0.01倍的短半轴长度，以确保能够准确描述边界层内的流动。设置合理的边界条件和参数，入口边界设置为速度入口，根据实际工况给定来流速度U_{\infty}；出口边界设置为压力出口，给定出口压力值；壁面边界设置为无滑移壁面，即流体在椭圆柱表面的速度为零。选择合适的湍流模型，由于椭圆柱绕流可能存在湍流现象，采用标准k-\epsilon湍流模型，该模型在处理一般的湍流流动时具有较好的准确性和稳定性。在计算过程中，利用Fluent的动网格技术，模拟椭圆柱在流体中的运动。设置椭圆柱的运动速度和方向，使其分别沿长轴方向和短轴方向运动，监测运动过程中椭圆柱所受的流体作用力。将数值计算结果与理论计算结果对比，分析两者之间的差异。在相同的工况条件下，对比沿长轴方向和短轴方向运动时的附加质量。对于沿长轴方向运动的情况，理论计算得到的附加质量为m_{xx理论}=\rho\pi\frac{(A+B)^{2}}{4}，数值计算得到的附加质量为m_{xx数值}，计算两者的相对误差\delta_{xx}=\frac{|m_{xx数值}-m_{xx理论}|}{m_{xx理论}}\times100\%。对于沿短轴方向运动的情况，同理计算相对误差\delta_{yy}。通过对比发现，在一定的误差范围内，数值计算结果与理论计算结果基本吻合。误差的来源主要包括网格划分的精度、湍流模型的近似性以及数值计算过程中的截断误差等。通过网格无关性验证和湍流模型的对比分析，可以进一步减小误差，提高数值计算的准确性，从而验证理论分析的正确性。3.3儒可夫斯基翼型附加质量分析3.3.1儒可夫斯基翼型的解析复势儒可夫斯基翼型在航空领域具有重要地位，其解析复势的推导基于儒可夫斯基变换。儒可夫斯基变换的数学表达式为z=\zeta+\frac{b^{2}}{\zeta}，其中z=x+iy是物理平面（z平面）上的复变量，\zeta=\xi+i\eta是变换平面（\zeta平面）上的复变量，b为常数。在推导儒可夫斯基翼型的解析复势时，假设\zeta平面上有一个半径为a（a>b），圆心位于原点的圆。当\zeta平面上的这个圆经过儒可夫斯基变换后，在z平面上就会形成儒可夫斯基翼型。在\zeta平面上，均匀来流绕半径为a的圆柱的复势函数为W(\zeta)=U_{\infty}(\zeta+\frac{a^{2}}{\zeta})，其中U_{\infty}为来流速度。将儒可夫斯基变换z=\zeta+\frac{b^{2}}{\zeta}代入W(\zeta)中，通过变量代换求解\zeta关于z的表达式。对z=\zeta+\frac{b^{2}}{\zeta}进行变形可得\zeta^{2}-z\zeta+b^{2}=0，根据一元二次方程求根公式\zeta=\frac{z\pm\sqrt{z^{2}-4b^{2}}}{2}。为了确定合适的根，需要根据变换的物理意义和边界条件进行判断。通常选择满足在无穷远处边界条件的根，这里选择\zeta=\frac{z+\sqrt{z^{2}-4b^{2}}}{2}（当z\to\infty时，\zeta\toz，符合无穷远处的流动特性）。将\zeta=\frac{z+\sqrt{z^{2}-4b^{2}}}{2}代入W(\zeta)=U_{\infty}(\zeta+\frac{a^{2}}{\zeta})，经过一系列的复数运算和化简（包括根式运算、通分、合并同类项等），最终得到z平面上儒可夫斯基翼型绕流的解析复势W(z)。在这个变换过程中，参数b起着关键作用，它决定了儒可夫斯基翼型的厚度和弯度等几何特征。当b的值发生变化时，\zeta平面上圆的变换关系也会改变，从而导致在z平面上生成的儒可夫斯基翼型的形状发生变化。增大b的值，会使翼型的厚度增加，同时弯度也会发生相应的改变，这是因为b的变化影响了z=\zeta+\frac{b^{2}}{\zeta}中\frac{b^{2}}{\zeta}这一项对z的贡献，进而改变了翼型的几何形状。参数a也对翼型的形状产生影响，a主要影响翼型的整体大小和一些细节特征。增大a的值，会使翼型在z平面上的尺寸相应增大，同时也会改变翼型表面的速度分布和压力分布，因为a在复势函数W(\zeta)=U_{\infty}(\zeta+\frac{a^{2}}{\zeta})中参与了速度势和流函数的计算，从而影响了整个流场的特性。3.3.2儒可夫斯基翼型的附加质量计算在得到儒可夫斯基翼型的解析复势W(z)后，根据动能与附加质量的关系公式T=\frac{1}{2}m^{*}U^{2}（其中T为流体动能，m^{*}为附加质量，U为物体运动速度）来计算附加质量。首先，根据解析复势W(z)求速度分量。速度分量u和v可通过复势函数对z的导数得到，即u-iv=\frac{dW}{dz}。由于儒可夫斯基翼型的解析复势W(z)是通过复杂的变换得到的，求导过程需要运用复合函数求导法则和复变函数的求导公式。对W(z)进行求导时，先对W(z)中的每一项分别求导，再根据复合函数求导法则进行组合。经过一系列复杂的求导和化简运算（包括对根式函数、分式函数的求导以及复数运算等），得到u和v在直角坐标系下的表达式。然后，计算流体的动能T。对于单位长度的二维流动，动能T=\frac{1}{2}\rho\iint_{S}(u^{2}+v^{2})dxdy，将速度分量表达式代入动能积分式中。通过将积分区域从z平面转换到\zeta平面，利用z与\zeta之间的变换关系以及相应的雅克比行列式，将积分转化为在\zeta平面上的积分。在积分过程中，需要根据\zeta平面上圆的边界条件确定积分区间，再运用复变函数的积分方法进行计算。经过复杂的积分运算（包括对复变函数积分性质的运用、积分区间的确定、积分变量的替换等），最终得到动能T的表达式。根据附加质量的定义，对比动能T的表达式和\frac{1}{2}m^{*}U^{2}，可以得到儒可夫斯基翼型的附加质量m^{*}的表达式。在计算过程中，考虑了流体为理想流体，即无粘性、不可压缩且流动无旋，同时忽略了流体的重力和其他外力的影响，仅考虑流体与儒可夫斯基翼型之间的惯性相互作用。假设翼型表面为物面边界，流体在翼型表面满足无穿透条件，即\frac{\partial\varphi}{\partialn}=0（\varphi为速度势函数，n为翼型表面的法向方向），这一条件在积分运算和附加质量计算中起到了重要作用，保证了计算结果的准确性。通过这样的计算过程，得到了儒可夫斯基翼型的附加质量表达式，该表达式反映了附加质量与流体密度、翼型的几何参数（如a、b等）以及来流速度之间的关系。3.3.3几何参数对附加质量的影响分析儒可夫斯基翼型的厚度、弯度等几何参数变化时，附加质量的变化规律，对于深入理解翼型在流体中的动力学特性具有重要意义。当儒可夫斯基翼型的厚度增加时，从附加质量的计算公式和流场特性分析可知，附加质量会增大。这是因为厚度的增加使得翼型周围扰动的流体体积增大，流体的惯性作用增强。在推导附加质量计算公式时，翼型的几何形状通过复势函数影响速度分布，进而影响动能计算。厚度增加，复势函数中的相关项发生变化，导致速度场中翼型周围的速度梯度增大，更多的流体被带动参与运动，从而使动能增加，根据T=\frac{1}{2}m^{*}U^{2}，附加质量m^{*}也随之增大。在实际航空应用中，对于飞机机翼，如果采用厚度较大的儒可夫斯基翼型，在飞行过程中，机翼周围的气流受到更大的扰动，附加质量增大，这会对飞机的飞行性能产生影响，如可能需要更大的动力来维持飞行速度和姿态。当翼型弯度改变时，附加质量也会相应变化。增大翼型弯度，会使翼型上下表面的压力差发生改变，进而影响流场的速度分布和动能。弯度的增加使得翼型上表面的气流速度更快，下表面的气流速度相对较慢，导致翼型周围的流场更加复杂，流体的惯性作用增强。从复势函数的角度来看，弯度的变化会改变儒可夫斯基变换中相关参数的关系，使得复势函数发生变化，从而影响速度场和动能的计算。当翼型弯度增大时，复势函数中的某些项会使翼型表面的速度分布发生改变，导致流体对翼型的作用力特性发生变化，动能增加，附加质量增大。在不同飞行攻角下，附加质量同样会发生变化。随着攻角的增大，翼型上下表面的压力差进一步增大，升力增加，同时流场中的分离现象可能加剧，导致翼型周围的流场更加不稳定。攻角的变化会改变翼型与来流的相对角度，从而影响复势函数中的速度势和流函数，使得速度场和压力场发生改变。在大攻角下，翼型表面可能出现气流分离，形成涡旋，这会导致更多的流体参与到复杂的运动中，动能增加，附加质量增大。但当攻角增大到一定程度时，可能会出现失速现象，此时附加质量的变化规律会更加复杂，需要进一步深入研究。3.3.4Fluent数值计算与分析为了进一步验证理论分析的正确性，并深入研究儒可夫斯基翼型的附加质量特性，利用Fluent软件对儒可夫斯基翼型进行数值模拟。首先，在建模软件中创建精确的儒可夫斯基翼型几何模型。根据儒可夫斯基变换的参数a和b，准确确定翼型的形状和尺寸，确保模型的准确性。将几何模型导入Fluent软件后，进行网格划分。考虑到儒可夫斯基翼型表面的曲率变化和边界层特性，采用非结构化网格进行划分，并在翼型表面和近壁区域加密网格，以提高计算精度，准确捕捉流场的变化。在翼型表面，网格尺寸设置为较小的值，如0.01倍的翼型弦长，以确保能够准确描述边界层内的流动。设置合理的边界条件和参数。入口边界设置为速度入口，根据实际工况给定来流速度U_{\infty}；出口边界设置为压力出口，给定出口压力值；壁面边界设置为无滑移壁面，即流体在翼型表面的速度为零。选择合适的湍流模型，由于儒可夫斯基翼型绕流可能存在湍流现象，采用标准k-\epsilon湍流模型，该模型在处理一般的湍流流动时具有较好的准确性和稳定性。在计算过程中，利用Fluent的动网格技术，模拟儒可夫斯基翼型在流体中的运动。设置翼型的运动速度和方向，监测运动过程中翼型所受的流体作用力。对模拟结果进行分析，重点关注压力分布、流场速度等与附加质量的关系。通过Fluent模拟得到翼型表面的压力分布云图，可以清晰地看到在不同位置处压力的变化情况。在翼型前缘，压力较高；在翼型上表面，压力逐渐降低，形成低压区；在翼型后缘，压力又逐渐升高。压力分布的这种变化与附加质量密切相关，压力差产生的合力会影响翼型的运动，进而影响附加质量。通过分析流场速度矢量图和流线图，可以了解流场中速度的分布和流动趋势。在翼型周围，速度分布不均匀，存在速度梯度较大的区域，这些区域的流体运动对附加质量有重要影响。将数值模拟得到的附加质量结果与理论计算结果进行对比，分析两者之间的差异。在相同的工况条件下，计算数值结果与理论结果的相对误差\delta=\frac{|m_{数值}-m_{理论}|}{m_{理论}}\times100\%。通过对比发现，在一定的误差范围内，数值模拟结果与理论计算结果基本吻合。误差的来源主要包括网格划分的精度、湍流模型的近似性以及数值计算过程中的截断误差等。通过网格无关性验证和湍流模型的对比分析，可以进一步减小误差，提高数值计算的准确性，从而验证理论分析的正确性，为儒可夫斯基翼型在工程实际中的应用提供更可靠的依据。3.4本章小结本章深入研究了圆柱、椭圆柱和儒可夫斯基翼型这三种单体结构在流体中的附加质量。通过基于解析复势的理论推导，得到了它们各自的附加质量表达式，明确了附加质量与结构几何参数、流体性质之间的关系。对于圆柱，附加质量与半径的平方成正比，且在二维平面运动中，沿x方向和y方向的附加质量相等，其表达式为m_{xx}=m_{yy}=\rho\pia^{2}，这表明圆柱半径的变化会显著影响附加质量的大小，半径增大，附加质量迅速增加，反映了圆柱周围扰动流体体积与附加质量的正相关关系。椭圆柱的附加质量则与长半轴A和短半轴B密切相关，沿长轴方向运动时附加质量m_{xx}=\rho\pi\frac{(A+B)^{2}}{4}，沿短轴方向运动时附加质量m_{yy}=\rho\pi\frac{(A-B)^{2}}{4}。长轴增长会使附加质量增大，短轴增长时，m_{xx}增大，m_{yy}的变化则较为复杂，先减小后增大，且椭圆柱沿长轴方向运动时的附加质量大于沿短轴方向运动时的附加质量，这体现了椭圆柱形状对附加质量的显著影响，以及不同运动方向下流体惯性作用的差异。儒可夫斯基翼型的附加质量受厚度、弯度等几何参数以及飞行攻角的影响。厚度增加，附加质量增大，这是由于厚度增大导致翼型周围扰动的流体体积增大，流体惯性作用增强；弯度改变会使翼型上下表面的压力差和流场速度分布发生变化，从而影响附加质量；随着飞行攻角的增大，附加质量也会增大，但攻角过大可能导致失速现象，使附加质量变化规律更为复杂。通过Fluent数值计算对理论分析结果进行了验证和进一步分析。在数值模拟过程中，通过合理设置边界条件、选择合适的湍流模型以及精确的网格划分，得到了与理论计算基本吻合的结果，验证了理论分析的正确性。同时，数值模拟还能够直观地展示流场的压力分布和速度分布，深入分析这些分布与附加质量之间的关系，为理解翼型和柱体在流体中的动力学特性提供了更丰富的信息。四、多体附加质量分析4.1串列双圆柱附加质量分析4.1.1串列双圆柱的解析复势在推导串列双圆柱的解析复势时，考虑两个半径均为a的圆柱，它们沿x轴方向串列分布，中心距为L。均匀来流速度为U_{\infty}，方向沿x轴正方向。首先，对于单个圆柱在均匀来流中的绕流，根据势流理论，其复势函数为W_1(z)=U_{\infty}(z+\frac{a^{2}}{z})，其中z=x+iy为复变量。当存在串列双圆柱时，利用Milne-Thomson圆定理和叠加原理来推导其解析复势。设第一个圆柱位于原点z_1=0，第二个圆柱位于z_2=L处。对于第一个圆柱，其复势为W_{1}(z)=U_{\infty}(z+\frac{a^{2}}{z})。对于第二个圆柱，由于它受到第一个圆柱的干扰，根据Milne-Thomson圆定理，其复势需要考虑第一个圆柱的影响。第二个圆柱单独存在时，在均匀来流中的复势为W_{20}(z)=U_{\infty}(z-L+\frac{a^{2}}{z-L})，但由于第一个圆柱的存在，其复势应为W_{2}(z)=U_{\infty}(z-L+\frac{a^{2}}{z-L})+\overline{U_{\infty}(\frac{a^{2}}{(z-L)}+\L)}。根据叠加原理，串列双圆柱的解析复势W(z)为两个圆柱复势之和，即W(z)=W_{1}(z)+W_{2}(z)=U_{\infty}(z+\frac{a^{2}}{z})+U_{\infty}(z-L+\frac{a^{2}}{z-L})+\overline{U_{\infty}(\frac{a^{2}}{(z-L)}+\L)}。在推导过程中，考虑了两圆柱之间的相互干扰对复势的影响。当两个圆柱距离较近时，它们之间的流体相互作用增强，这种相互作用体现在复势函数中。在W_{2}(z)的表达式中，\overline{U_{\infty}(\frac{a^{2}}{(z-L)}+\L)}这一项就是考虑了第一个圆柱对第二个圆柱的干扰。第一个圆柱的存在改变了第二个圆柱周围的流场，使得第二个圆柱周围的速度势和流函数发生变化，从而影响其复势。这种干扰在实际工程中有着重要的意义，例如在海洋平台的支撑结构中，如果采用串列双圆柱结构，圆柱之间的相互干扰会影响整个结构在海流中的动力学响应，通过准确推导解析复势，可以更好地理解和预测这种影响。4.1.2双圆柱的附加质量计算在得到串列双圆柱的解析复势W(z)后，根据动能与附加质量的关系公式T=\frac{1}{2}m^{*}U^{2}（其中T为流体动能，m^{*}为附加质量，U为物体运动速度）来计算双圆柱各自的附加质量。首先，根据解析复势W(z)求速度分量。速度分量u和v可通过复势函数对z的导数得到，即u-iv=\frac{dW}{dz}。对W(z)=U_{\infty}(z+\frac{a^{2}}{z})+U_{\infty}(z-L+\frac{a^{2}}{z-L})+\overline{U_{\infty}(\frac{a^{2}}{(z-L)}+\L)}求导，由于W(z)的表达式较为复杂，求导过程需要运用复合函数求导法则和复变函数的求导公式。经过一系列复杂的求导和化简运算（包括对分式函数、根式函数的求导以及复数运算等），得到u和v在直角坐标系下的表达式。然后，计算流体的动能T。对于单位长度的二维流动，动能T=\frac{1}{2}\rho\iint_{S}(u^{2}+v^{2})dxdy，将速度分量表达式代入动能积分式中。通过将积分区域从z平面转换到\zeta平面（若需要进行坐标变换），利用z与\zeta之间的变换关系以及相应的雅克比行列式，将积分转化为在\zeta平面上的积分。在积分过程中，需要根据边界条件确定积分区间，再运用复变函数的积分方法进行计算。经过复杂的积分运算（包括对复变函数积分性质的运用、积分区间的确定、积分变量的替换等），最终得到动能T的表达式。根据附加质量的定义，对比动能T的表达式和\frac{1}{2}m^{*}U^{2}，可以得到串列双圆柱各自的附加质量。设第一个圆柱的附加质量为m_{1}^{*}，第二个圆柱的附加质量为m_{2}^{*}，它们可以表示为一个矩阵形式：\begin{pmatrix}m_{11}^{*}&m_{12}^{*}\\m_{21}^{*}&m_{22}^{*}\end{pmatrix}其中m_{11}^{*}表示第一个圆柱自身运动时的附加质量，m_{22}^{*}表示第二个圆柱自身运动时的附加质量，m_{12}^{*}和m_{21}^{*}表示两个圆柱之间的相互附加质量，反映了两圆柱之间的相互干扰对附加质量的影响。m_{11}^{*}和m_{22}^{*}的大小与单个圆柱在均匀来流中的附加质量以及两圆柱之间的相互作用有关。当两圆柱距离较远时，相互作用较弱，m_{11}^{*}和m_{22}^{*}接近单个圆柱的附加质量；当两圆柱距离较近时，相互作用增强，m_{11}^{*}和m_{22}^{*}会发生明显变化。m_{12}^{*}和m_{21}^{*}则体现了一个圆柱的运动对另一个圆柱附加质量的影响，其大小和正负与两圆柱的相对位置、运动方向以及来流条件等因素有关。4.1.3圆柱半径和间隔距离对附加质量的影响分析两圆柱半径变化以及它们之间间隔距离改变时，附加质量的变化趋势，对于理解串列双圆柱在流体中的动力学特性具有重要意义。当两圆柱半径a增大时，从附加质量的计算公式和流场特性分析可知，附加质量会增大。对于单个圆柱，附加质量与半径的平方成正比，在串列双圆柱中，随着半径的增大，每个圆柱周围扰动的流体体积增大，流体的惯性作用增强，从而使各自的附加质量m_{11}^{*}和m_{22}^{*}增大。从相互附加质量m_{12}^{*}和m_{21}^{*}来看，半径的增大也会使两圆柱之间的相互作用增强，因为半径增大导致圆柱周围流场的影响范围扩大，两圆柱之间的流体相互干扰更加明显，使得m_{12}^{*}和m_{21}^{*}的绝对值可能增大。在实际工程中，对于海洋平台的串列圆柱支撑结构，如果增大圆柱半径，在海流作用下，圆柱的附加质量增大，对平台的稳定性和受力情况会产生显著影响，可能需要加强结构设计以承受更大的流体作用力。当圆柱间隔距离L改变时，附加质量也会发生相应变化。当间隔距离L增大时，两圆柱之间的相互作用减弱。从各自的附加质量m_{11}^{*}和m_{22}^{*}来看，它们会逐渐接近单个圆柱在均匀来流中的附加质量，因为间隔距离增大，圆柱之间的干扰减小，每个圆柱周围的流场逐渐趋近于单个圆柱绕流的情况。对于相互附加质量m_{12}^{*}和m_{21}^{*}，随着间隔距离的增大，其绝对值会减小，这表明两圆柱之间的相互影响逐渐减弱。当间隔距离减小到一定程度时，两圆柱之间的相互作用增强，m_{12}^{*}和m_{21}^{*}的绝对值增大，可能会导致两圆柱的运动状态发生明显变化，例如在换热器管束中，如果串列圆柱的间隔距离过小，流体在管束中的流动会受到强烈干扰，导致附加质量增大，影响换热器的性能。在不同工况下，如不同的来流速度、流体密度等条件下，圆柱半径和间隔距离对附加质量的影响规律基本相似，但影响的程度可能会有所不同。在高流速的来流条件下，圆柱半径和间隔距离的变化对附加质量的影响可能更加显著，因为流速增大使得流体的惯性作用增强，圆柱周围的流场更加复杂，从而放大了圆柱半径和间隔距离变化对附加质量的影响。4.1.4Fluent数值计算与验证为了验证理论分析的正确性，利用Fluent软件对串列双圆柱进行数值计算。首先，在建模软件中创建精确的串列双圆柱几何模型，确保两个圆柱的半径a和中心距L的尺寸准确无误。将几何模型导入Fluent软件后，进行网格划分。考虑到圆柱表面的曲率变化和两圆柱之间的流场特性，采用非结构化网格进行划分，并在圆柱表面和两圆柱之间的区域加密网格，以提高计算精度，准确捕捉流场的变化。在圆柱表面，网格尺寸设置为较小的值，如0.01倍的圆柱半径，在两圆柱之间的区域，根据中心距的大小合理调整网格尺寸，以确保能够准确描述两圆柱之间的流体相互作用。设置合理的边界条件和参数。入口边界设置为速度入口，根据实际工况给定来流速度U_{\infty}；出口边界设置为压力出口，给定出口压力值；壁面边界设置为无滑移壁面，即流体在圆柱表面的速度为零。选择合适的湍流模型，由于串列双圆柱绕流可能存在湍流现象，采用标准k-\epsilon湍流模型，该模型在处理一般的湍流流动时具有较好的准确性和稳定性。在计算过程中，利用Fluent的动网格技术，模拟串列双圆柱在流体中的运动。设置圆柱的运动速度和方向，监测运动过程中圆柱所受的流体作用力。将数值计算结果与理论结果对比，分析两者之间的差异。在相同的工况条件下，对比两个圆柱各自的附加质量以及相互附加质量。对于第一个圆柱的附加质量m_{1}^{*}，理论计算得到的值为m_{1理论}^{*}，数值计算得到的值为m_{1数值}^{*}，计算两者的相对误差\delta_{1}=\frac{|m_{1数值}^{*}-m_{1理论}^{*}|}{m_{1理论}^{*}}\times100\%。同理，计算第二个圆柱的附加质量相对误差\delta_{2}以及相互附加质量的相对误差\delta_{12}和\delta_{21}。通过对比发现，在一定的误差范围内，数值计算结果与理论计算结果基本吻合。误差的来源主要包括网格划分的精度、湍流模型的近似性以及数值计算过程中的截断误差等。通过网格无关性验证和湍流模型的对比分析，可以进一步减小误差，提高数值计算的准确性，从而验证理论分析的正确性，为串列双圆柱在工程实际中的应用提供更可靠的依据。4.2儒可夫斯基翼型与圆柱附加质量分析4.2.1儒可夫斯基翼型与圆柱的解析复势在推导儒可夫斯基翼型与圆柱组合结构的解析复势时，考虑一个儒可夫斯基翼型和一个圆柱的相互作用。设儒可夫斯基翼型通过儒可夫斯基变换由\zeta平面上半径为a_1（a_1>b_1），圆心位于原点的圆变换得到，其变换关系为z=\zeta+\frac{b_1^{2}}{\zeta}。圆柱半径为a_2，圆心位于z_0处。均匀来流速度为U_{\infty}，方向沿x轴正方向。对于儒可夫斯基翼型，在\zeta平面上，其绕流的复势函数为W_1(\zeta)=U_{\infty}(\zeta+\frac{a_1^{2}}{\zeta})。通过儒可夫斯基变换z=\zeta+\frac{b_1^{2}}{\zeta}，将\zeta平面上的复势函数转换到z平面，得到儒可夫斯基翼型在z平面上的复势W_1(z)。在转换过程中，需要通过求解\zeta关于z的表达式，即对z=\zeta+\frac{b_1^{2}}{\zeta}变形为\zeta^{2}-z\zeta+b_1^{2}=0，然后根据物理意义和边界条件选择合适的根\zeta=\frac{z+\sqrt{z^{2}-4b_1^{2}}}{2}（当z\to\infty时，\zeta\toz，符合无穷远处的流动特性），再代入W_1(\zeta)进行化简。对于圆柱，根据Milne-Thomson圆定理，其复势函数为W_2(z)=U_{\infty}(z-z_0+\frac{a_2^{2}}{z-z_0})+\overline{U_{\infty}(\frac{a_2^{2}}{(z-z_0)}+z_0)}。考虑两者的相互干扰，根据叠加原理，组合结构的解析复势W(z)为儒可夫斯基翼型复势W_1(z)与圆柱复势W_2(z)之和，即W(z)=W_1(z)+W_2(z)。在推导过程中，难点主要在于处理儒可夫斯基翼型和圆柱之间的相互干扰。由于两者的形状和位置关系复杂，它们之间的流场相互作用难以直接描述。为解决这个问题，通过精确的数学变换和理论分析，利用儒可夫斯基变换将儒可夫斯基翼型的复杂形状转化为便于分析的圆形，再结合Milne-Thomson圆定理考虑圆柱的影响，运用叠加原理将两者的复势相加，从而得到组合结构的解析复势。这种方法充分利用了已有的理论成果，通过合理的数学处理，有效解决了相互干扰带来的难题，为后续附加质量的计算奠定了基础。4.2.2儒可夫斯基翼型与圆柱的附加质量计算在得到儒可夫斯基翼型与圆柱组合结构的解析复势W(z)后，依据动能与附加质量的关系公式T=\frac{1}{2}m^{*}U^{2}（其中T为流体动能，m^{*}为附加质量，U为物体运动速度）来计算各自的附加质量。首先，根据解析复势W(z)求速度分量。速度分量u和v可通过复势函数对z的导数得到，即u-iv=\frac{dW}{dz}。由于W(z)是儒可夫斯基翼型复势W_1(z)与圆柱复势W_2(z)的叠加，求导过程需要分别对W_1(z)和W_2(z)求导，再运用复合函数求导法则和复变函数的求导公式进行组合。对W_1(z)求导时，因为它是通过复杂的儒可夫斯基变换得到的，需要对变换后的表达式进行细致的求导运算；对W_2(z)求导时，要考虑其包含分式函数和共轭复数的特点，运用相应的求导规则。经过一系列复杂的求导和化简运算（包括对根式函数、分式函数的求导以及复数运算等），得到u和v在直角坐标系下的表达式。然后，计算流体