固体温度场的随机边界元分析

针对遭受随机热边界条件的结构,提出并建立了一种固体温度场分析的随机边界元算法.算例表明本文建立的随机边界元法是有效的.     

适用于Kappel桨的面元法尾涡模型

Kappel桨在梢部剧烈弯折,这样的几何结构能够抑制梢涡的形成,起到梢部加载的作用。为了模拟桨叶弯折后对整个桨叶受力分布的影响,本文通过CFD计算分析几何尺寸相似的弯折三维翼和平直三维翼的剖面环量分布,并以此为基础就面元法分析Kappel桨时的尾涡进行修正。利用修正尾涡模型的面元法计算Kappel桨的敞水性能,与未修改尾涡面的面元法计算结果作比较,并就面元法计算Kappel桨的Morino库塔条件和等压库塔条件进行相关说明。     

翼/板结合部涡旋流动结构与壁面脉动压力的大涡模拟研究

壁面脉动压力是湍流非定常特性的重要表征,而且是重要的流体动力声源,所以对于脉动压力的研究已经成为流声耦合研究领域的重要课题。文章采用大涡模拟方法结合动态Smagorinsky亚格子涡模型,数值计算了两型翼/板结合部在不同流速下的湍流脉动压力,并展示了结合部马蹄涡流动结构。通过将计算结果与试验结果进行对比分析可知,数值计算方法能准确地预报湍流脉动压力,从而为理解翼/板结合部非定常流动的物理机理奠定了基础。     

大型集装箱船迎浪航行的外飘砰击压力研究

该文采用直接计算法研究一艘大型集装箱船迎浪航行时的外飘砰击压力特点。先基于三维势流理论Rankine面元法,应用DNV船级社的商业软件Wasim计算出该船在压载时不同海况下以不同航速航行时在非规则波中时域内的非线性运动响应时历,再根据Wagner型冲击理论计算得到该集装箱船艏外飘区域的外飘砰击压力。最后,依据概率随机理论,得到不同工况下的外飘砰击概率以及外飘砰击压力极值,并分析概率和极值特点。     

考虑涡泄模式的低质量比圆柱涡激运动数值模拟

本文与经典实验进行对比,对二维圆柱在不同来流速度下的涡激运动进行研究。将运动系统简化为质量(m)-弹簧(k)-阻尼(ζ)系统,分析浮式圆柱运动的控制方程并通过4阶Runge-Kutta法求解运动微分方程,借助UDF编程嵌入到Fluent求解器中进行求解,结合动网格技术实现流固耦合,对比Jauvtis和Williamson的经典实验以验证数值模拟的可靠性,再现了SS,2S,2T和2P的涡泄模式。虽得到了与Jauvtis和Williamson物理模型实验相近的一些典型结果,没有做更进一-步的研究,但通过实验数据与数值模拟的对比,能够区别研究相同质量比下其他不同参数对涡激运动特性的影响。对m^*=1的浮式圆柱与低质量比圆柱m^*-2.6进行涡激运动数值模拟研究对比,发现不同约化速度下对运动频率、涡泄模型等涡激运动特性均有不同程度的影响。     

柔性管异形表面结构对其涡激振动影响的数值研究

涡激振动（VIV）是导致海洋立管疲劳破坏的主要诱因,因此,抑制振动是立管设计需关注的重要内容。本文通过运用双向流-固耦合数值计算方法,开展翼型管、纹理管与光滑管的涡激振动响应特性的对比研究,探讨扰流翼板与表面纹理等异形表面结构对涡激振动被动抑制的有效性。计算结果表明：异形表面结构能够有效降低柔性管体的振动响应,缩短锁振区间,并且影响尾流区泄涡发放模式,减小管体流体力系数,起到明显抑制管体涡激振动的作用,且在较高的速度工况下,这种抑制效果更加显著。该研究可为采用被动式振动抑制的翼型管和纹理管等异形结构的海洋立管设计提供一定的理论参考。     

基于OpenFOAM的无转角和有转角楔形体舱段入水砰击载荷数值模拟研究

合理预报船首局部结构砰击载荷是船舶设计和研究者重点关心的问题。文章针对无转角和有转角楔形体舱段入水砰击问题,采用OpenFOAM开源软件,开发了基于两相流求解器InterDyMFOAM数值模拟结构入水砰击的程序包,开展了预报研究工作。数值模拟了入水过程的压力和加速度时域响应历程,与入水砰击模型实验进行了对比分析,吻合很好。同时讨论了网格划分等对结果的影响和砰击的三维效应。结果表明,本文提出的入水砰击载荷预报方法计算效率高,具有推广应用到预报复杂三维船首结构砰击问题的潜力。     

基于一种固体区域迭代算法的圆柱涡激振动数值计算

利用Fluent平台的用户自定义程序（UDF）以及动网格模型,实现了圆柱运动方程的一种迭代求解算法,分别对层流、湍流状态下,弹性支承圆柱体在一定约化速度下的涡激响应进行了数值模拟,探讨了不同阻尼比对涡激响应的影响。结果表明：采用该迭代求解算法对弹性支承圆柱涡激振动的预测结果较为合理;随着阻尼比的逐渐增加,初始支振幅、升阻力系数时程曲线将由多频率拍振,最终变为单一频率主导的振动,且涡激振幅逐渐减小;除了质量-阻尼比联合参数m*ζ外,阻尼比ζ本身也应作为一个重要的涡激影响参数单独进行考量。     

不同间距比下串联圆柱涡激振动数值模拟研究

针对不同间距比下小尺寸串联双圆柱在双自由度下的涡激振动位移响应,受力特性以及运动轨迹等进行仿真研究,应用Fluent中的SST k-ωMesh技术进行网格划分,运用动网格技术实现流固耦合。研究结果表明当S/D=3时,上游、下游两圆柱的“锁定”区间均为Vr=7-9,当S/D=4与S/D=5时,两串联圆柱的“锁定”区间范围增大。同时随着串联两圆柱间距比增大,上下游圆柱之间的影响逐渐减小,导致升力系数相应减小。随着约化速度增大,上下游圆柱的阻力系数均趋于稳定。当S/D=3时串联两圆柱运动的规律性较差,随着间距比增大上游圆柱“8”字形变得更加标准,下游圆柱轨迹则呈现多样化的特点。     

近水面椭球体非线性兴波的数值模拟

本文开发了一种去奇异面元法用于计算以定常前进速度运动的近水面椭球体引起的三维非线性自由面势流.采用求解非线性深水波的Stokes理论求解该定常非线性问题;物面用NURBS(非均匀有理B样条)表达;孤立的Rankine源分布在物体内部和自由面上方.给出了波型、波阻以及收敛特性随一些参数的变化规律.波阻计算结果和他人的其它数值结果进行了比较,其吻合程度令人满意.     

螺旋桨旋涡发放数值模拟

基于LSDYNA软件和任意拉格朗日—欧拉（ALE）方法,建立螺旋桨与流体相互作用的有限元模型,数值计算螺旋桨固定桨叶在前方来流下的流体动力特性和结构力学性能,分析总结桨叶后方旋涡运动和桨叶内部应力变化的规律。结果表明,螺旋桨桨叶的不规则形状导致了桨叶后方旋涡的相互作用,桨叶单元所在位置的厚度以及该位置与约束位置的距离决定了单元应力的大小。数值模拟得出的结论对研究螺旋桨＂唱音＂现象有一定的参考价值。     

面元法预报螺旋桨表面非定常压力分布

本文建立了扰动速度势面元法预报螺旋桨表面非定常压力分布的理论和数值方法。该方法把桨叶和桨毂表面离散为若干四边形双曲面元,每个面元上布置等强度源汇和偶极子分布,螺旋桨尾涡面也离散为布置等强度偶极子的四边形双曲面元。所有的时域参数通过傅里叶级数展开转化为频域参数,使得在时间域内的求解转化为每一谐调阶上的求解。桨叶随边处通过迭代和采用广义逆矩阵方法在每一谐调阶上实现非线性等压库塔条件。桨叶表面非定常压力     

船桨干扰定常空化性能数值模拟

文章采用CFD技术计算了船桨干扰定常空化性能。首先计算了不同空化数下的NACA66(MOD)水翼和DTRC4381螺旋桨的定常空化性能,同时预报了KCS船的伴流场、阻力性能和船体表面压力分布。为了验证计算船后螺旋桨空化流的可行性,文中对KCS船和KP505螺旋桨进行了整体考虑,并计算了不同空化数下的船后螺旋桨的空化覆盖面积。计算表明,文中采用的计算方法合理,计算结果与试验结果吻合。     

平台运动影响下张力腿涡激振动特性研究

张力腿涡激振动易引起系泊系统的疲劳损伤,危害张力腿平台的安全性和可靠性。该文借鉴细长立管涡激振动的相关研究,同时考虑张力腿与立管在截面尺度、预张力大小以及边界条件存在的巨大差距,采用切片理论结合CFD数值模拟方法,对平台运动影响下张力腿涡激振动开展了研究。考虑定常张力与时变张力对涡激振动的影响,采用GAMBIT软件分区建立了多切片计算网格,将编制求解涡激振动的UDF程序嵌入Fluent软件中,采用动网格技术实现流场的更新并计算作用于张力腿上的瞬时升力和拖曳力。文中还比较了仅考虑流和考虑浪、流联合作用张力腿的时变张力影响,研究了两种典型波高共计六种工况,计算结果表明：随着流速的增加,各模态权重比例将发生跳转现象;尽管由于平台运动导致张力腿受力的随机性,但平台运动对张力腿涡激振动影响较小。     

DTMB5415驱逐舰绕流的数值模拟

本文采用RANS方程及RNG湍流模型数值模拟了带方尾及声纳罩的5415驱逐舰绕流，计算结果表明对复杂的实际船形，CFD能有效地预报船体周围流场、波形和阻力等船舶设计感兴趣的水动力性能参数，为船形分析及优化提供了有力的工具。     

Level-set法预报液舱内剧烈晃荡引起的冲击压强

基于Level-set法就部分装载液体的液舱内剧烈晃荡引起的冲击压强进行了二维数值计算.为了模拟大幅晃荡引起的冲击压强,重点对数值计算较为敏感的自山液面处的物性参数、初始化方程的求解等进行了,数值处理.通过对某一矩形液舱内的液体晃荡的计算,得到了自由表而位置和压强的时间历经曲线.计算结果同实验值的比较表明:本方法可以用于计算预报大幅晃荡引起的载荷.     

基于面元法的吊舱推进器定常性能研究

基于势流理论面元法建立了吊舱推进器定常性能的计算方法.分别建立螺旋桨和吊舱的积分方程,通过在表面上布置双曲面元将方程离散为以面元上偶极强度为未知量的矩阵.螺旋桨和吊舱之间的相互影响通过迭代计算来处理.Newton-Raphson迭代过程被用来在桨叶随边满足压力Kutta条件.为避免数值求导中的奇异性,用柳泽(Yanagizawa)方法求得物体表面的速度分布.支架作为升力体处理,并通过迭代计算更新支架的尾涡形状.计算了拖式吊舱推进器的定常水动力性能,与实验结果的比较表明,计算误差在5%以内.分析了舱体对螺旋桨的影响,舱体的伴流会引起螺旋桨的载荷增大.     

基于UMAT的冰—结构相互作用数值仿真

基于ABAQUS 用户自定义子程序UMAT,采用已有的冰的弹性失效准则,针对冰与结构的相互作用进行数值模拟。首先, 通过简单算例验证本文开发的UMAT的正确性;然后, 采用开发的UMAT建立单个冰单元受刚性板挤压下的力学模型,得到力—位移曲线;最后, 通过建立圆柱与冰相互作用的仿真分析模型,探讨冰板厚度对结构反力的影响。研究结果表明,采用用户自定义程序UMAT使得模拟更为复杂的冰力学特性成为可能,可为以后开展大规模的冰—结构相互作用研究奠定基础。