水下复杂圆柱壳体结构流固耦合效应研究概述
管道流固耦合振动及声传播的研究现状及展望
第5卷第2期船舶力学V ol.5N o.2 2001年4月Journal of Ship M echanics A pr.2001文章编号:1007-7294(2001)02-0082-09管道流固耦合振动及声传播的研究现状及展望刘忠族,孙玉东,吴有生(中国船舶科学研究中心,无锡214082)摘要:本文对管道振动及声传播领域的研究现状和水平分类作了比较全面的综述,并对今后开展潜艇管道振动、声传播及噪声控制的研究进行了展望。
关键词:管道系统;流固耦合作用;管道振动;声传播中图分类号:O353.1,O422文献标识码:A1引言管道振动、噪声及其控制技术的研究有着广泛的应用背景。
对于工程上的管道系统,其动力学及声学分析非常复杂,不同的流体与结构物理模型的组合可派生出不同的动力学问题。
管道系统涉及的问题很多,如水锤,流固耦合振动及稳定性,管道声传播等。
水锤直接影响到管道的安全性,同时也会激发管道系统振动。
当管内流体速度增加到某一临界值时,系统将发生失稳,这种失稳是有害的,它可以使管道系统的薄弱环节产生破坏;另外,颤振失稳也会辐射噪声。
管道振动及噪声研究领域有许多理论问题没有解决,例如:复杂空间管系振动传递特性分析,管道系统对激励的响应及衰减特性,液压系统的动态特性分析以及管道声传播及辐射特性等。
管道振动噪声研究具有重要的学术价值和工程意义。
2管道振动及声传播的研究现状2.1管道振动及噪声源的描述管道噪声的含义是:管壁结构的振动和管内流体中声波相互作用,并沿管壁和管内流体传播与辐射的过程。
作用于管道的激励有两种:力激励与声激励。
力激励主要是指泵作用于管壁与流体的动态力,声激励是指泵的水动力噪声与阀门等处的湍流噪声。
湍流噪声是由控制阀和流量调节器产生的,在控制阀和流量调节器处,通流面积突然缩小或扩大,由此产生湍流,在下游伴随着压力脉动。
这种脉动压力在流体中向下游传递,从而形成平面波而成为湍流噪声[1]。
对湍流噪声的机理已有许多研究。
CAE-LS-DYNA水下爆炸数值仿真
多物质组分的初始化
DYNA程序提供多物质组分的初始化功能流体算法对网格质量要求很高全六面体网格！流体物质的网格离散质量直接关系到求解精度DYNA里面的多物质流体网格可保持最佳质量*INITIAL_VOLUME_FRACTION_GEOMETRY基于模型现有结构或是自定义几何构形来分配全六面体映射网格界面单元可以实现不同组分的插值分配处理避免流体传输常见的“负体积”等数值错误，保证分析精度
LS-DYNA水下爆炸数值仿真
程序发展概况
背景及地位◆ LS-DYNA系列为国际主要的显式非线性动力分析程序, 广泛用于冲击、碰撞、爆破及流固耦合研究，为北约新武器开发的重要工具。◆ LS-DYNA程序被公认为是显式有限元分析的鼻祖和理论先导，是目前所有显式求解程序（包括显式板成型程序）的基础代码。◆ 丰富的接触碰撞算法和卓越的多物质流固耦合功能使其在分析高度瞬态非线性问题时具有其它软件无法比拟的优势。◆ 军工系统爆轰分析的首选软件，事实上的行业标准。
式中：ω、A、B、R1、R2是输入参数，适用于各种凝态炸药。
水或空气的状态方程
Jones-Wilkins-Lee 状态方程
*EOS_LINEAR_POLINOMIAL*EOS_GRUNEISEN*EOS_IDEAL_GAS
水下爆炸涉及的命令
多物质及耦合关键字*section_solid_ale (选择单元算法)*control_ale（指定物质输运算法）*ale_multi-material_group_options （多物质ALE(Euler)分析时指定）*initial_void_options （单物质+空PART时指定）*contact_sliding_only_options 炸药采用Lagrange描述*constrained_lagrange_in_solid 炸药采用ALE或Euler描述，结构采用Lagrange描述*initial_detonation 单点、多点、线起爆*control_explosive_shadow 隔爆、绕爆
部分浸没圆柱壳声固耦合计算的半解析法研究
部分浸没圆柱壳声固耦合计算的半解析法研究郭文杰;李天匀;朱翔;屈凯旸【摘要】部分浸没圆柱壳-流场耦合系统的声振分析是一种典型的半空间域内声固耦合问题,其振动及声学计算目前主要依赖于数?方法求解,但无论从检验数?法还是从机理上揭示其声固耦合特性,解析或半解析方法的发展都是不可或缺的.本文提出了一种半解析方法,先将声场坐标系建立在自由液面上,采用正弦三角级数来满足自由液面上的声压释放边界条件;接着基于二维Flügge薄壳理论建立了以圆柱圆心为坐标原点的壳-液耦合系统的控制方程;然后再利用Galerkin法处理声固耦合界面的速度连续条件,推导得到声压幅?与壳体位移幅?之间的关系矩阵并求解该耦合系统的振动和水下声辐射.与有限元软件Comsol进行了耦合系统自由、受迫振动和水下辐射噪声计算结的对比分析,表明本文方法准确可靠.本文的研究为解析求解弹性结构与声场部分耦合的声振问题提供了新的思路.【期刊名称】《物理学报》【年(卷),期】2018(067)008【总页数】13页(P143-155)【关键词】圆柱壳;自由液面;Galerkin法;正弦三角级数【作者】郭文杰;李天匀;朱翔;屈凯旸【作者单位】华中科技大学船舶与海洋工程学院, 武汉 430074;船舶与海洋水动力湖北省重点实验室, 武汉 430074;华中科技大学船舶与海洋工程学院, 武汉 430074;高新船舶与深海开发装备协同创新中心, 上海 200240;船舶与海洋水动力湖北省重点实验室, 武汉 430074;华中科技大学船舶与海洋工程学院, 武汉 430074;高新船舶与深海开发装备协同创新中心, 上海 200240;船舶与海洋水动力湖北省重点实验室, 武汉 430074;上海交通大学船舶海洋与建筑工程学院, 上海 200240【正文语种】中文1 引言圆柱壳结构以其优异的几何特性、力学特性被广泛应用于土木、化工、航空航天、海洋工程等诸多领域,比如输油管道、储液罐、潜艇耐压壳等.20世纪,国内外针对圆柱壳-流场耦合系统声振问题已有大量的文献报道[1−4],也形成了较为完善的理论体系,在这些工作中,大多数的研究均针对比如全充液或者全浸没问题(无限流域)这种壳体与流场完全耦合的振动、声学问题的研究.但是对于实际的工程问题,流域往往是存在边界的,尤其是当结构靠近流体边界时,或者结构内带有自由液面时,如果忽略流体边界的影响,计算结果往往是有偏差的[5].自由液面作为一类典型的边界条件,其与圆柱壳结构的耦合动力学特性吸引了许多学者对其开展研究.Huang[6]基于镜像原理研究了平面波入射下二维圆柱壳的散射声场,Hasheminejad和Azarpeyvand[7]进一步将该方法延展到研究平面波入射下1/4空间二维圆柱的声散射问题.白振国等[8]采用镜像原理建立了有限水深环境中二维圆柱壳的振动声学物理模型,并探讨了潜深对声场分布和衰减特性的影响规律.郭文杰等[9,10]基于镜像原理和Graf加法定理建立了有限水深下有限长圆柱壳-流体耦合振动的解析模型,并基于傅里叶变换及稳相法提出了该系统远场声压的快速准确预报方法.上述研究主要针对壳体结构完全浸没在水中的情况,这一类问题中,结构完全与流场耦合,其声振耦合方程的描述相对容易.实际上,还有一类问题,就是结构与流场部分耦合的情况,比如系泊状态下的潜艇、舰船等.对于系泊状态下圆柱壳声固耦合分析的研究工作则相对较少,尤其是圆柱壳轴线平行于自由液面时的解析或半解析研究更是鲜有报道,这是因为系泊状态下流体载荷解析表达式和耦合方程更难以得到.为了解决壳体部分充液的振动问题,Amabili[11]给出了两种近似的方法,一种是用以圆柱圆心为原点构成的扇形边界近似替代自由液面的方法,但是这种方法仅适用于浸没角度较小的情况;另一种是用部分环状区域代替原始边界,但这种方法也仅适用于浸没角度小于π的工况.值得注意的是,Amabili提出的第一类方法是可以推广到壳体部分浸没问题中[12],并明确指出浸没角度的适用范围仅为−π/8—π/8;但是第二类方法无法推广到外流场.基于Amabili提出的自由液面近似处理方法,Ye等[13]研究了系泊状态下无限长圆柱壳声振特性.此外,其他一些学者也研究过部分浸没或者部分充液的壳体耦合振动问题,Selmane 和Lakis[14]将圆柱壳沿周向分割成微段,再结合波动法研究了部分充液时其流固耦合特性.Ergin和Temarel[15]基于瑞利-李兹法研究了部分充液或者部分浸没时圆柱壳固有频率.Krishna和Ganesan[16]基于多项式级数及有限元法计算了部分充液圆柱壳固有频率以及附连水质量.但是这些研究工作实际上仅仅考虑了湿表面流体的影响而并未考虑自由液面边界条件的声学效应,尚无法揭示自由液面对声振特性的影响.另外还有一些学者通过实验的手段分析了部分充液的圆管自振特性[17]或者附连水质量效应[18],但依然无法从机理上求解这类问题.需要强调的是,在壳体与流场部分耦合的问题中,半充液或者半浸没情况是个特例,因为这种模型下自由液面刚好在壳体横截面的水平坐标轴上,容易得到流体载荷的解析表达式,可采用正弦三角级数[19]来自动满足自由液面的边界条件.基于这类解析表达式,Li等[20]对半浸状态下全充液圆柱壳的声辐射特性进行分析,李天匀等[21]针对有限长半充液圆柱壳的自由振动及受迫振动特性开展研究.尽管本文工作是要解决部分浸没这类更为普遍的工况,而文献[19—21]中的研究工作仅能解决半浸没这种特殊工况.但是,文献[19—21]中关于在半浸没情况下通过采用正弦三角级数来自动满足自由液面下的声压释放条件的思路却可以借鉴到本文的研究工作中来.针对更一般的壳体与流场部分耦合,如部分浸没等问题,此时自由液面与圆柱壳的圆心不共面.这给问题的描述和求解带来了很大的挑战.为解决这个问题,本文将声压函数和壳体位移函数建立在不同的坐标系下.具体来讲,将声场坐标系的原点建立在自由液面上,并采用正弦三角级数来满足自由液面上的声压释放边界条件;将壳体运动方程建立在以圆柱圆心为坐标原点的坐标系下.由此,便可以分别得到声压函数与壳体位移函数在各自坐标系下的解析表达式,为解析求解部分浸没问题奠定了基础.然后,再利用Galerkin法处理结构与声场部分声固耦合的界面上的速度连续条件,并通过两坐标间的几何关系求得声压幅值与壳体位移幅值之间的关系矩阵,最终便可以求解该声固耦合系统方程.本文方法的提出可以便捷有效地解决系泊状态下圆柱壳的声振问题,丰富了该问题理论研究的内涵,也为解析求解浮态问题提供了新的思路.2 理论分析2.1 物理模型为了便于研究,本文假设圆柱壳轴向是无限长的,且激励力沿轴向是均匀分布的,由此本文的数学物理模型是一个典型的平面应变模型(即二维模型).二维圆柱壳厚度为h,中面半径为Rs,杨氏模量为E,泊松比为µ;密度为ρs,部分浸没于密度为ρf;声速为cf 的流体中,壳体轴系与自由液面平行.如图1所示,以壳体圆心O为坐标原点的极坐标系(r,φ)为结构坐标系,φ的取值范围为0—2π,以壳体圆心正上方与自由液面的交点Q为原点的极坐标系(R,θ)为声场坐标系,θ的取值范围为0—π,图1 物理模型坐标图Fig.1.Coordinate figure of the physical model.定义Q点到O点距离为浸没深度H,浸没角度α满足sinα=H/Rs,当自由液面在壳体圆心下方时,H取值为负值.两类坐标系与平面上任一点的夹角定义为β.当自由液面在壳体圆心上方时β取值为正值,在下方时β取值为负值.径向激励力幅值为F0,激励角度为φ0.2.2 声学边界条件本文的物理模型是典型的声固耦合模型,而这项研究工作的难点和重点之一就是得到满足所有对应的声学边界条件的声压解析表达式.首先,声压p必须满足声学Helmholtz方程:其中kf=ω/cf为声波波数,ω=2πf为角频率,f表示频率,∇2表示拉普拉斯算子.其次,声压表达式还要满足无限远处Sommerfeld辐射条件:最后,自由液面处的声压需满足声压释放条件(因为研究的频率相对较高,自由表面的重力波可以忽略不计;此外由于自由液面以上空气密度远小于水,故作为真空处理)实际上,当声场坐标系原点建立在自由液面上时,可以通过采用正弦三角级数来自动满足自由液面声压为零的边界条件[19],具体形式为其中m为正弦三角级数的序数,Pm(R)为对应的声压幅值函数.对于自由液面上的任意点,角度θ=0或π,将其代入(4)式,可以得到自由液面声压的表达式:显然此类正弦三角级数可以用来满足自由液面声压释放条件,此外采用三角函数也更有利于应用分离变量法来求解Helmholtz方程.由此可以得到声压幅值函数Pm(R)的解析表达式:其中H(1)m()为第m阶第一类汉克尔函数,Am为声压幅值.由于第一类汉克尔函数在远场自动满足Sommerfeld辐射条件,因此将(6)式代入(4)式可以得到满足以上声学边界条件的声压解析表达式:2.3 壳体运动方程得到声压的解析表达式后,接下来需要建立壳体运动方程.本节采用二维Flügge薄壳理论[22],具体的方程为(为简洁,后文略去简谐时间项exp(−iωt)):其中v和w分别是壳体中面切向和径向位移,fp表示作用在圆柱壳表面的声载荷,f0表示外激励力载荷,[L]是二维Flügge薄壳方程中的微分算子矩阵,具体如下:由于壳体和声介质是部分耦合,因此(8)式中fp应表示为分段函数的形式:假设作用于圆柱壳的激励力是一个沿轴向均匀分布的无限长径向线力,激励力位于结构坐标系的(Rs,φ0)处,则激励力载荷可以表示为其中δ()表示狄利克雷函数.对于圆柱壳结构,由于周向的周期性,其位移及载荷函数可以在周向展开为傅里叶级数的形式[8]:其中Vn和Wn分别是周向和径向位移幅值,fpn和f0n分别表示壳体表面声载荷fp以及激励力载荷f0的幅值,n是傅里叶周向展开序列数.因为(9)和(13)式是壳体表面声载荷fp的不同形式的表达式,利用正交化处理,可得fpn的表达式:由图1可知,空间中任意点的结构坐标与声学坐标有如下关系:因此对于流固耦合面上任一点,令r=Rs,则(R,θ)在声学坐标系下的坐标可以由(16)式解出:将(17)式代入(15)式即可对fnp进行求解,但是由于积分中包含汉克尔函数,无法直接进行积分计算,因此本文采用离散求和的形式来近似计算.首先将积分域均分为K 段,然后取各段中点值代入求和公式中:其中F(φ)表示需要进行积分的函数,通过大量算例表明K取100时收敛性已经很好. 同理,(10)和(14)式是外激励力f0的不同形式的表达式,利用正交化处理,可以得到f0n的表达式:然后将(11)—(14)式代入到(8)式中,并进行正交化处理,可以得到解耦后的壳体运动方程:其中矩阵[T]的元素如下:T11=Ω2−n2,T12=T21=in,T22=1+K+Kn4−2Kn2−Ω2.是无量纲频率.由(20)式可以得到仅与径向位移幅值相关的控制方程:其中In=T11/det(T),In与n有关,det(T)表示矩阵[T]的行列式.很明显,求解控制方程(21)的关键在于得到径向位移幅值Wn和壳体表面声载荷幅值fpn之间的关系.因此需要通过壳体表面流体与结构速度连续条件来求解该问题,在声固耦合交界面处.具体方程为[10]:因为很难直接求解速度连续方程,采用Galerkin法进行处理,可选择的权函数有两类,一类是壳体径向位移的周向函数,另一类是声压的周向函数:因此方程(22)可以转变为Galerkin积分的弱形式:其中N是截断项数,即需要构造2N+1个积分方程.由(24)式可以得到径向位移幅值和声压幅值之间的关系:其中[Vs]和[Vp]均是2N+1阶方阵,{Wn}和{Am}分别表示径向位移幅值向量和声压幅值向量,且{Wn}=[W−N,W−N+1,···,WN−1,WN]T,{Am}=[A1,A2,···,A2N,A2N+1]T,上标T 表示转置.根据(16)式中两类坐标系的几何关系,可以将声压沿径向的导数转换到声学坐标系下:其中夹角β=3π/2−θ−φ.由此,当权函数选择壳体径向位移的周向函数exp(inφ)时,[Vs]和[Vp]中每一个元素的具体表达式如下:其中a,b分别表示矩阵的行和列.当权函数选择声压的周向函数sin(mθ)时,[Vs]和[Vp]中每一个元素的具体表达式如下:由于(27)或(28)式中积分包含汉克尔函数,无法直接进行积分计算,因此也采用(18)式中离散求和的方法来近似计算.为了方便求解耦合方程,将(21)式改写为矩阵运算方程:其中{f0n}=F0/2π{exp(iNφ0),exp[i(N − 1)φ0],···,exp(−iNφ0)}T;[G]矩阵为对角矩阵,[G]j,j=Ij−1−N;{fpn}={fp,−N,fp,−N+1,···,fp,N}T.同时,(15)式中{fpn}也可表示为矩阵的形式:{fpn}=[Tp]{Am}. (30)其中矩阵[Tp]中每一个元素的具体表达式如下:将(25)和(30)式代入方程(29)中,可以求解出径向位移幅值{Wn}:其中矩阵[J]是2N+1阶单位矩阵,并且在得到径向位移幅值{Wn}后,根据(25)式可以求解出{Am},从而可以计算声压.前文求解受迫振动时,是已知激励力及激励频率,求响应(径向位移幅值).当求解自由振动时,并没有激励源,固有频率是要求解的未知量.由此(32)式可以表示为其中{0}表示零向量,显然(33)式是个典型的特征值问题,即求解角频率ω.具体来讲就是通过定义(33)式中系数矩阵行列式值为零来求解角频率:其中F(ω)表示系数矩阵的行列式值.因为矩阵[J],[Vp],[Tp]中均含有角频率ω,实际上矩阵的行列式值F(ω)=0是个关于角频率隐式表达的超越方程,难以直接解出ω.因此本文采用搜根的方法进行求解,通过设置合适的搜根步长,逐步增大输入的角频率,当行列式值过零点时,输出对应的角频率,这个角频率即为系统的固有角频率.由此可以将超越方程的各阶角频率从小到大逐一求出,再根据可以得到各阶固有频率.3 数值计算模型参数:半径Rs=0.18 m,厚度h=0.001 m,壳体密度ρ=7850 kg/m3,杨氏模量E=206 GPa,泊松比µ=0.3,流体密度ρf=1025 kg/m3,流体声速cf=1500 m/s.3.1 收敛性分析为了说明方法的收敛性,取无量纲浸没深度H/Rs= −0.5,H/Rs=0和H/Rs=0.5,分别计算在激励力频率200,400,800 Hz时径向均方根振速Vm随截断数N的变化规律,定义其中vn= ∂w/∂t表示径向速度;则均方根速度级 VML=20lg(Vm/V0),其中基准速度V0=10−6m/s. 激励力幅值F0=1 N,激励角度φ0= π/4. 定义复杨氏模量E′=E(1+iη),结构阻尼η=0.01.计算中权函数选择位移的周向展开函数.图2 均方根速度级VML收敛性分析Fig.2.Convergence analysis of the root mean square velocity levels.从图2可见,均方根速度级随着截断项数N的增大很快趋于稳定;并且频率越高,达到收敛时截断项数N的取值越大.由图2可知,对于频率小于800 Hz时的分析计算,N取16时已足够收敛.3.2 本文方法的适用性分析对于Amabili[12]提出的自由液面近似处理方法,浸没角度α的适用范围仅为−π/8到π/8(无量纲浸没深度H/Rs约为−0.38—0.38,无量纲深度为负值表示壳体的圆心在自由液面上方),而本文方法的主要优势就是模型中壳体部分浸没时的浸没深度范围更广.为了验证本文的方法,无量纲浸没深度取值为−0.9—0.9,分别计算两种不同权函数下首阶固有频率值,并与有限元软件Comsol仿真计算结果进行对比.其中有限元模型如图3所示,流域以声学坐标原点为中心,半径取2 m,用完美匹配层模拟无限远声学边界,匹配层厚度取0.05 m.网格包含7980个域单元和756个边界单元.图3 有限元模型Fig.3.The finite element model.表1 不同浸没深度下首阶固有频率对比(单位为Hz)Table parison of the fundamental frequency with dif f erent immerged depth(in Hz).H/Rs Present(取位移权函数求解)Present(取声压权函数求解)FEM(Comsol)−0.9 15.45 15.57 15.57−0.8 15.06 15.07 15.16−0.7 14.33 14.51 14.45−0.6 13.24 13.31 13.33−0.5 12.65 12.67 12.72−0.4 12.36 12.36 12.43−0.3 12.24 12.2412.31−0.2 12.19 12.19 12.26−0.1 11.44 11.45 11.51 0 10.72 10.72 10.77 0.1 10.08 10.08 10.13 0.2 9.52 9.53 9.57 0.3 9.06 9.06 9.10 0.4 8.67 8.67 8.71 0.5 8.35 8.35 8.39 0.6 8.10 8.10 8.14 0.7 7.91 7.91 7.95 0.8 7.76 7.74 7.80 0.9——7.67从表1可以看出,选择位移函数或者选择声压函数作为权函数时,首阶固有频率计算结果符合很好.并且无论选择哪类权函数,无量纲浸没深度H/Rs均可以从−0.9变化到0.8(当无量纲浸没深度过大或者过小时计算结果难以收敛),这表明本文的浸深适用范围非常广,方法更具有一般性.此外,从表1还可以看出,本文方法计算结果与Comsol仿真计算结果符合良好,而且不同方法下首阶固有频率随浸没深度的变化规律也是一致的.这是由于随着壳体浸没深度的增大,流固耦合面增大,附连水质量也相应增加,因此固有频率会逐渐减小.3.3 准确性验证3.3.1 自由振动的验证为说明本文方法计算自由振动问题的准确性,分别取无量纲浸没深度H/Rs=−0.7和H/Rs=0.7,计算系统前十阶固有频率,并与有限元软件Comsol仿真计算结果进行对比.定义固有频率的相对误差Error=|f1−f2|/f2×100%.表2 H/Rs=−0.7时前十阶固有频率对比(单位为Hz)Table parison of thethe first ten order natural frequencies when H/Rs=−0.7(in Hz).阶数Present(f1) FEM(f2) Error/%14.33 14.45 0.83 2 15.01 15.10 0.59 3 43.66 43.94 0.64 4 45.32 45.75 0.94 5 84.41 85.09 0.80 6 90.83 91.61 0.85 7 144.34 145.77 0.98 8 144.90 145.93 0.71 9 213.75 215.40 0.77 10 219.91 219.90 0 1表3 H/Rs=0.7时前十阶固有频率对比(单位为Hz)Table parison of the the first ten order natural frequencies when H/Rs=0.7(in Hz).阶数Present(f1) FEM(f2) Error/%7.91 7.95 0.50 2 8.86 8.90 0.45 3 25.59 25.74 0.58 4 28.26 28.42 0.56 5 55.30 55.63 0.59 6 57.74 58.10 0.62 7 96.25 96.88 0.65 8 97.65 98.26 0.62 9 147.42 148.36 0.63 10 150.77 151.78 0.67 1从表2和表3可以看出,本文方法计算得到的前10阶固有频率值与Comsol仿真计算结果符合良好,最大相对误差不超过1%,说明采用本文方法计算固有频率是准确可靠的.3.3.2 受迫振动的准确性验证在分析完自由振动的准确性之后,进一步分析受迫振动的准确性. 取无量纲浸没深度H/Rs= −0.5和H/Rs=0.5,计算在激励力频率50—500 Hz时测点径向速度(频率间隔10 Hz),定义径向速度级VL=20×log(|V|/V0),其中V为径向速度,基准速度V0=10−6m/s.激励力幅值F0=1 N,激励角度φ0=0.测点位于周向角φ=π处.从图4可以看出,本文计算结果和Comsol仿真计算结果整体符合良好,说明本文方法分析受迫振动是准确可靠的.但是当频率较高时(以图4为例,大于400 Hz),误差逐渐增大.主要的原因可能是有限元计算声固耦合问题时,随着频率增大,对网格密度的要求也提高,计算精度会降低.图4 不同方法下径向速度级对比(a)H/Rs=−0.5;(b)H/Rs=parison of the radial velocity levels with dif f erentmethods:(a)H/Rs=−0.5;(b)H/Rs=0.5.3.3.3 声场求解的准确性验证为进一步说明本文方法计算声场也是准确的,取无量纲浸没深度H/Rs=−0.4和H/Rs=0.4,计算激励频率为50 Hz时声压幅值云图,并与Comsol仿真计算结果进行对比.其中激励力幅值F0=1 N,激励角度φ0=0,云图尺寸为1.2 m×0.6 m.从图5和图6可以看出,本文方法计算得到的声压云图和有限元软件Comsol仿真计算结果符合良好,由此可以说明本文方法计算声场是准确可靠的.图5 H/Rs=−0.4时不同方法下声压幅值云图对比 (a)本文方法;(b)Comsol仿真parison of the sound pressure contour map with dif f erent methods when H/Rs=−0.4:(a)Present method;(b)Comsol.图6 H/Rs=0.4时不同方法下声压幅值云图对比 (a)本文方法;(b)Comsol仿真parison of the sound pressure contour map with dif f erent methods when H/Rs=0.4:(a)Present method;(b)Comsol.此外,值得一提的是,本文方法计算效率也非常高,以图5或图6中声压云图的计算为例,在Matlab中仅仅需要不到2 s即可计算出精确稳定的结果.4 讨论在验证了本文方法计算声固耦合系统振动及声压求解均准确可靠之后,进一步分析与讨论部分浸没壳体的模态振型、无量纲浸没深度对受迫振动的影响以及远场声压的指向性特征.4.1 模态振型分析为了更直观地揭示声固部分耦合系统自振特性,取无量纲浸没深度H/Rs=0.7和无限域时计算得到的前4阶模态振型及固有频率.对比表4和表5可以看出,部分浸没工况下振型与无限域振型有差异.以H/Rs=0.7时第一阶振型为例,振型函数约为cos2φ+0.25·cosφ−0.05·cos3φ,而对于无限域壳体,对应的振型函数为cos2φ(规则的周向波型).这是由于流体在周向分布不均,壳体与结构的部分耦合破坏了圆柱壳的周向的对称性,故规则的周向波之间会发生互耦,形成复杂的振型函数.表4 H/Rs=0.7时前4阶模态振型Table 4.Modal shapes of the first four orders when H/Rs=0.7. 1 2 3 47.91 Hz 8.86 Hz 25.59 Hz 28.26 Hz表5 无限域时前4阶模态振型Table 5.Modal shapes of the first four orders when in infinite fluid. 1 2 3 46.34 Hz 6.34 Hz 20.37Hz 20.37Hz此外,由于自由液面的存在,系统仅有惟一对称轴,对称和反对称模态固有频率也存在差异.以H/Rs=0.7算例下前两阶固有频率为例,分别为7.91 Hz和8.86 Hz,有明显差异;但是对于无限域情况,由于系统具有周向对称性,对称和反对称模态固有频率是相同的,二者存在明显区别.另外,部分浸没工况下,对称和反对称模态振型函数之间的周向波互耦程度也不尽相同.仍以H/Rs=0.7工况为例,第二阶振型函数(反对称)约为sin2φ+0.38·sinφ−0.08·sin3φ,函数各成分之间的比值由对称模态时1:0.25:−0.05变为反对称时1:0.38:−0.08.而对于无限域情况,对应的振型函数为sin2φ,这也是系泊状态下自振特性区别于无限域工况的又一特征.4.2 均方根振速频谱分析3.2节中提到,随着浸没深度增大,流固耦合面会增大,附连水质量也会相应增加.为进一步解释这个观点,本文选取了H/Rs=−0.6,−0.3,0,0.3,0.6这五个深度,对比分析各工况下均方根速度级的频谱曲线.激励力幅值F0=1 N,激励角度φ0=0,激励频率为1—100 Hz,扫频间隔1 Hz.从图7可以看出,随着浸没深度增大(液面升高),共振峰均向低频移动.这也是因为浸没深度增大使得附连水质量增大,从而增加了系统的总质量,导致共振频率降低,频谱曲线整体左移.4.3 远场声压指向性分析系泊状态下的舰船等目标的水下辐射噪声对于其隐蔽性有着重要意义.因此,本节进一步开展二维系泊圆柱壳远场声压指向性特征的研究.取无量纲浸没深度H/Rs=−0.5,0,0.5,分别计算激励频率f=100,200 Hz时远场声压级指向性图.其中激励力幅值F0=1 N,激励位置φ0=π/4.远场点取自声学坐标系,半径R=1000 m,角度θ取0—π,取值间隔θ/180.定义声压级SPL=20lg(|p|/p0),其中基准声压p0=10−6m/s.从图8可以看出,所有声压指向性曲线都在角度θ=π/2时(正下方)取最大值.而且尽管激励位置并不在壳体的垂直对称轴上,但声压级关于垂直对称轴呈对称分布.以下对这些现象进行物理解释.图7 不同浸没深度下均方根速度级频谱曲线(a)H/Rs=−0.6,−0.3,0;(b)H/Rs=0,0.3,0.6Fig.7.Spectrum curves of the root mean square velocity levels at dif f erent immersiondepths:(a)H/Rs=−0.6,−0.3,0;(b)H/Rs=0,0.3,0.6.图8 不同浸没深度下声压级指向性 (a)f=100 Hz;(b)f=200 HzFig.8.Directivity of the sound pressure levels at dif f erent immersion depths:(a)f=100Hz;(b)f=200 Hz.因为研究的场点位于远场,从几何上讲,壳体湿表面上任意一点到场点的距离可以认为是近似相同的.由此水下的辐射面可以等效为一个点源,而且这个点源距离自由液面的距离在0—2Rs之间(这个距离与速度分布有关,是未知的).基于镜像原理[8],自由液面对水下声场的作用可以通过在自由液面的另一侧构造等距的(虚源到液面距离)以及反相位的虚源来实现.则水下声场可以认为是由实源和虚源共同贡献,且虚源和实源共同作用下声压在自由液面上满足为零的条件.实际上,由于实源和虚源距离较近,在低频下,这样的物理模型构成了偶极子模型.由此,可以得到。
流固耦合系统计算方法及应用软件概述
显地分开 ， 使描述物理现象 的方程 ， 特别是本构 方程需要 针对具 体的物理 现象来建立 ， 其耦 合效应 通过 描述 问题的微 分 方程而 体现。第 二大类 问题 的特征是耦合作 用仅仅发 生在两相 交界向 上， 在方程上 耦合 是 由两 相耦 合 面的 平衡 及协调 关 系 引人 的。 这类 问题 又 可 以具 体 的 分 为 三 种 情 况 。第 一 种 足 流 体 与固 体 结
问题 。张庆华 。 不考虑液体 的波动所 产生 的动压力 ， 液体静 将
压力 作 为外 载 荷 通 过 Ａ Ｓ Ｓ模 拟 加 载 ， 用 有 限 元 模 拟 出 弹 性 ＮＹ 利 罐 体 的 变形 情 况 。王 华 ， 刚 讨 论 了含 液 容 器 流 固耦 合 动 力 曹 特性 的有 限元 方 程 ， Ａ Ｓ Ｓ软 件 对 正 方 形 充 液 容 器 进 行 了模 用 ＮＹ 态分 析 ， 出 其 理 论 模 态 特 性 ， 果 表 明 该 模 型 有 效 可 靠 。徐 得 结
有描述流体现象 的变量及 描述 固体现象 的变 量。一般 而 言 ， 具 有 以下 两 点 特 征 ： （） １ 流体域 或固体域均不 可能单独求解 ； （） ２ 无法显式地 消去描述 流体运 动的独立 变量 或描述 固体
运 动 的 独 立变 量 。
・
３・ ６
广州化 工
２１ 年 ３ ０ １ ９卷第 １ 期
１ ２ 化工机械充液系统流 固耦合计算方法的发展概述 ．
充液系统 是 由固体和液体混合组成 的具有无 限多 自由度 的力学 系统 。在工程 实际 中， 化工 机械充液 系统 中液体 晃动 及 与 结 构 的耦 合 振 动 是 一 个 非 常 常 见 而 且 重 要 的 问 题 。这 类 耦 合 问 题 由 于 涉 及到 流 体 力 学 、 体 力 学 、 力 学 、 算 力 学 以及 高 动 计 等数学等学科 ， 其机理 和本 质的探讨非 常复 杂和 困难。一般 根 据腔体旋转方式 不同 ， 液体 晃动分为 绕对称轴 旋转运 动和 三 把 维 液 体 运 动 。如 航 天 器 上 的 液 体 燃 料 贮 箱 、 体 火 箭 等 属 于 前 液 者 ， 对 于 像 油 罐 车 、 液 汽 车 、 机 油 箱 等 这 类 储 运 系 统 ， 于 而 载 飞 属
气田高压注水过程中井身应力分布特征的流-固耦合分析
孔 的应力 集 中分析 ， 弹塑 性 材 料 壳体 上 开 多 孔 的应 力集 中分 析是 十分 复 杂 的 力 学 问题 ， 目前 国 内外 研
究较 少 。
１ 地层－ 水泥环－ 套管系统物理模型
对 拟 分析 的 注 水 井 的 套 管 射孑 Ｌ 段建立地层． 水
些油气田， 套 管毁 损 的注 水 井 占总套 管 毁 损 井 的
度， 对模型施加不 同的边界应力 Ｊ 。应力边界条件
为， 模型上部施加垂直地应 力 ， 水平方 向施 加水平
最 大 和最小 地 应 力 ， 井 眼 内 部施 加 注 水 压 力 ， 岩 石 内部 施 ／ Ｊ ｕ ｆ ； Ｌ 隙 压 力 。位 移 边 界 条 件 为 模 型底 部 和
侧面施加法向位移约束 ， 整个模型保证不产生刚性
应力 边界 条件 根 据 原 始 地应 力 资 料 ， 通 过 改 变
模 型的 温度 、 材料密度 、 力学参数 、 边界位移条件 ，
渗流场存在流固耦合作用 ， 是不可忽略的 ． ４ ； ② 套 管可以视为圆柱壳体 ， 但射孔造成很复杂的应力集
中问题 ， 国 内外 大部 分 研 究 都 是 基 于 弹 性 范 围 内单
科
学
技
术
与
工程 Βιβλιοθήκη １ ３卷 件， 联 立上述 方 程 就 可 以 得 到 所 求 应 力 场 ， 以 上 求
解 通 常需要借 助泛 函变 分法得 到近 似解 。 ２ ． ３ 渗流场 数学模 型
为 了研 究 射 孔 段 地层 、 水泥环 、 套 管 作 为 一 个 整体 系统 在 注 水 压 力 变 化 的情 况 下 其 各 部 分 应 力 分 布特 征 ， 假 定 岩 石 为 多 孔介 质 ， 孔 隙 间充 满 单 相 不可压 缩饱 和流 体 ， 整个渗流过程符合 达西定律 ，
基于灰色理论的水下双层壳体振动与声辐射统计能量分析
ＭＡ ｕｈｎ ，Ａ ｎ－ｉ Ｙ ＯＸｏ －ａｇ Ｙ Ｎ ａｎ２ ＩＯＸ －ｏｇ Ｔ ＯＪ ｇｑ ｏ Ａ ｉ ｇｌｎ２ Ａ Ｇ －ａ ｉ ａ。 ｎ ｉ ， Ｎ
（ Ｎ ｖｌ ｃｄｍｙ ｆ ｒ ａ ｅ ｔ ｅ ｉｇ１０ ７ ｈｎ； ａｂｎＥ ｇｎｅｉｇＵ ｉｅｓｙ Ｈ ｒｉ １０ ０ 。 ｈｎ ） １ ａａ Ａ ａｅ ｍ ｍ ｎ， ｉｎ ０ ０ ３Ｃ ｉａ ２Ｈ ｒｉ ｎｉｅｒ ｎｖｒｔ， ａｂｎ ５０ １ Ｃ ｉａ ｏＡ Ｂｊ ｎ ｉ
进 行分析 计算 。
收稿 日 ： ０ ４ ０— ３ 期 ２ ０— ６０
作者简介 ： 旭弘（９ ７ ）男 ， 士， 缪 １６ 一 ， 博 海军装备研究 院舰船所 高级工程师。
ｈ ｏ ｓ ｌ ｎ ｔｉ ｔ ｄ ｈ ｓｈ ｓ ｒｏ ｒｉ ．ｉ ａ ｙｇａ ｏｅ ａ ｔｉ ｓｄ ｔ ｅ ｄ ａａ ｔｒ， ｔｅｃ ｎ ｕｔｉ ｈｓ ｍｅ ｏ ａ ｉｈ ｅｒｒｍａｇｎＦ ｎ ｌ，ｒｙ ｆｒｃ ｓ ｓｕ ｅ ｏ ｇ ｔｔ ｅｍｏ ｅｐ ｒｍｅｅｓ ｈ ｌ ｈ
Ａ ｓ ａ ｔ Ｅ Ｓａｓｃｌ ｎｒ ｎ ｙｉ ｉａｐｉ ｎｌ ｅ ｎ ａ ｕ ｔｔｅ ｌ ｕｄ ｓ ｉ ｃｕｌ ｂｔ ｃ： Ａ（ ｔｔａ Ｅ ｅ ｙ ａ ｓ）ｓ ｐｌ ｄｔ ａａ ｚ ａｄ ｃｌ ｌ ｅｈ ｉ ｉ ｏ ｄ ｏｐｅ ｒ Ｓ ｔｉｉ ｇ Ａ ｌ ｓ ｅ ｏ ｙ ｃ ａ ｑ － ｌ ｄ
Ｖｏ ＿ Ｎｏ２ ｌ１ ０ ．
Ａ ｒ２ ０ ｐ． ０ ６
基于灰色理论的水下双层壳体振动与 声辐射统计能量分析
缪旭 弘 ・ ，陶景桥 ， 熊亮 姚 ， 杨娜娜
（ 海军装备研究院 北京 １０ ７ ； １ ００ ３ ２哈尔滨工程大学 ， 哈尔滨 １０ ０ ） ５０ １ 摘要 ： 本文采用统计能量 分析对一双层 圆柱壳体在水下的流固耦合振动 和声辐射进行分析计算 。首先对各系统 的模态参数采 用理论方法计算 。 与实验结果对 比表明 ， 误差较大 。为此 ， 本文提 出用灰色理论预测各子系统模态
圆柱形弹性壳液耦合系统模态的有限元法分析
阵和罐体刚度矩阵， 对于不可压缩流体 ， 若不考虑 自 由面重力 波 和结构 阻尼 的影 响 ， 方程 （ ） 则 １可简 化为
（ Ｍ ＋ ） ｄ＋ Ｋ ０ ＝ （） ２
式中
代表流体对 固体的作用 ， 以固体 附加质 现
量 的形式 出现 ， 时 流固耦 合 问题 就退 化为 考虑 附 这 加质量 的固体动 力 学 问题 ， 得 进 行 有 限元 分 析 大 使
计算 了两个 同轴 圆柱 壳 液耦 合 系 统 的频 率 和振 型 ； 张 素霞 ＿等 主要 针对存 在 自由液 面 的对称 弹性 矩形 ２ ］ 薄壳 充液 系统产 生 的流 固耦 合 振 动 问题 ， 了模 态 做 实验 ， 研究 了弹性 矩 形 薄 壳 充 液系 统 的 固有频 率 和 振 型 ； 到 了液体深 度对 整个耦 合系 统频 率 的影 响。 得
图 ５ 四 阶液 固耦 合 系统 模 态
模态 参数 分析结 果如 表 １ 。
表 １ 模态参数分析结果
一
３ 结
论
通 过 有 限元 法 对壳 液耦 合 系 统 的振 动 的分 析 ， 得 到 了与实验 现 象及 分 析 相 一致 的结果 ， 进一 步 证 实了实 验结论 的 可靠 性 ， 一 方 面也 表 明 用有 限元 另 法分析 壳液耦 合 问题具 有较 高 的精 度 。
模态 分 析 是研 究 复杂 机 械振 动 的重 要方 法 ， 通 过模态 分析 能够 得 到 实 际结 构 的振 动特 征 ， 而 大 从 大简化 了系 统 的数 学运 算 ， 在工 程 实 际 中得 到 广 泛 的应用 ， 到 了工 程界 的极 大重视 。 受
对于壳液耦合系统 的模态分析问题 ， 许多专家 学 者做 了大 量 的工作 ，． ｒｃｋ ＪＨｏａｅ…等 研 究 了装 有 液 体的薄壁 圆柱壳的模态振动特性 ， 给出圆环壳 内液 体 深度 和液体 层厚 度 对 系 统频 率 的影 响 ， 理 论 上 在
基于ANSYS储罐的建模研究
基于ANSYS储罐的建模研究摘要: 本文介绍了液固耦合的基本概念和储罐的有限单元类型，并利用大型有限元软件ANSYS建立了储罐液固耦合的有限元模型，为进一步开展储罐结构的静力和动力性能分析奠定基础。
关键词: 液固耦合储罐ANSYS 有限元模型Abstract：This paper introduces the basic concept of solid-liquid coupling and finite element type of tank, solid-liquid coupling finite element model of tank is established by finite element software ANSYS, and it lay the foundation for further developing the static and dynamic performance analysis of tank structure.Key Words: solid-liquid couplingtankANSYSfinite element model1.引言石油是工业的血液，在国民生产生活中有着举足轻重的作用[1]。
作为石油生产加工运输的重要设备储罐的抗震性能的好坏，就决定了石油工业的安全。
储罐是由管壁、底板和储液三部分组成，受力性能较复杂，大型有限元软件ANSYS 用于分析这种复杂结构的静力、动力、线性、非线性等响应特征时具有强大优势，可以很好地反映这种结构在各种复杂因素作用下的力学特征。
本文利用大型有限元软件ANSYS建立了储罐液固耦合的有限元模型，详细介绍了整个建模过程，为进一步开展储罐的静力和动力性能分析奠定基础。
2.液固耦合的基本概念流固耦合力学的研究对象是固体在流场作用下的各种行为以及固体变形或运动对流场的影响。
其重要特征是在于两相介质之间的相互作用，固体在流体动载荷作用下产生变形或运动，而固体的变形或运动反过来影响流场，从而改变流体载荷的分布和大小。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
抗 表达 式 以及 舷 间流 体 、 肋板 对 内外 壳体 的作 用 实 力表 达式 ， 到 了 内外 壳体耦 合 振动控 制方 程 ， 得 利用
稳 相法 求 得远 场 声压 。刘 涛 ， 陈美 霞 ， 东平 ， 骆 姚
３ ３
连 接 的外壳 采用 纵肋 加强 的双 层加 肋舷 间液 舱充 水 的圆 柱 壳体 的 振动 和 声辐 射 特 性 ， 导 了环 肋 的 阻 推
阻 抗 的阻 性 分量 代表 由结 构辐 射 的 声 能 ， 是流 体 对
结 构 附加 阻尼 的体现 ； 声辐 射 阻抗 的抗 性部 分 是 结 构 和流 体 之 间进 行 能量 交 换 的部 分 ， 流 体 附加 质 是
合 面上 的平 衡协 调条件 这几 个耦合方 程 。
般 都 是 以真 空 中的壳 体理 论为基 础进 行建模 。按 照
壳 体弹 性体 的控制 方 程 不 同 ， 壳体 理论 可分 为 严格 弹性 理论 方法和 薄壳近 似理 论方 法 。严 格弹 性理论 按 照各 向同性 弹 性体 所 满足 的方 程 导 出 ， 主要 用 于 处理 厚壳情 况 。薄壳 理论主 要针 对厚度 半径 比很 小 的壳 体 ， 用于 薄壳 的声振特 性 分析 ， 适 常用 的壳体理
主 ， 高于 吻 合 频 率 时 ， 体 的 效 应 以辐 射 阻 尼 为 而 流 主 。对 于有 限长 圆柱 壳 体 ， 以采 用 两 端加 半 无 限 可
长 的刚性 障柱来模 拟 ，ａｄ ｎⅢ 出半 无 限长 的刚 Ｓｎｍａ 指
性障柱对 有 限长 圆柱 壳表面 声压 沿轴 向分 布有 一定 影响 ， 对模态 声辐 射 阻抗 的影 响可 以忽 略 。 而 对 于 加肋 壳 体 ， 加 强筋 的处 理有 以下 几 种方 对 法 ： 一种 是 ｓ ａ 的方 法 ， 环肋 和 纵肋 的刚度 和 第 ｍｅｒ 将 偏 心 平均 分 配 到壳 体 上 ， 体 的材 料被 当作 正交 各 壳 向异性 ， 种 方 法仅 适 合 于低 频 以及密 集 加 肋 的情 这 况； 第二种 是利 用 能量变 分原 理 ， 到 壳体 结构 的运 得 动方 程 ； 三 种是 模态 展 开法 ， 立加 强筋 的运 动 第 建 方程 ， 将位 移和 内力 利用 壳体 的模态 展 开 ， 反力 以 将
Ｋｅ ｒ ｓ： ｉｒ ｔ ｎ ａ ｄ ｗａ ｅ； ｏ ｌ ｘ ｕ ｄ ｒ ａｅ ｙ ｉ ｄ ｉ ａ ｈ ｌ ｓ ｕ ｔｒ ｆ ｉ ・ｔ ｃｕ ｅ ｉｔ ｒ ｃｉ ｎ ； ｏ ｐ ｅ ｙ ｗｏ ｄ ｖ ｂ ａ ｉ ｎ ｖ ｏ ｃ ｍｐ ｅ ｎ ｅ ｗ ｔｒｃ ｌ ｒｃ ｌ ｅｌ ｔ ｃ ｅ； ｕ ｄ ｓ ｕ ｔ ｒ ｎ ｅ ａ ｔｏ ｎ ｓ ｒ ｕ ｌ ｒ ｃ ｕ ｌｄ
１ 水 下 圆柱 壳体 结构 研 究 现 状
对 水下 壳体 的研 究有无 限长 和有 限长之分 。一
柱 壳体振 动 声辐射 的物 理本质 是结 构受外 载荷 激励
而 振动 ， 在流 体介质 中辐射 形 成声场 ， 声场 反过 来对 结构 施加 反 作用 力 ， 成 声振 耦合 系 统 。研 究 水 中 形 舱 筏 的耦 合 振 动 和 声 辐 射 问题 包 括 结构 动 力 学 方 程 ， 场 声学 方程 及 Ｓ ｍｍｅｆｌ辐 射 条件 和 流 固耦 流 ｏ ｒ ｄ ｅ
熊亮 等研 究 了有 限长加 肋双 层 圆柱壳 模态 声辐射 特 性， 壳体振 动采 用壳 体理 论方 程描 述 ， 强筋 的作 用 加 等 效 为反 力和 反力矩 ， 内外 壳之 问 的流 体 载荷 采 用 Ｈｅ ｏｚ 动 方 程 和 壳 体 表 面 的 边 界 条 件 进 行 求 ｌ ｌ波 ｍｈ
中 图 分 类 号 ： ６ ． ； Ｂ５３ Ｕ６ １ ４ Ｔ ３ ４ 文 献标 识 符 ： Ａ ＤＯ 编码 ：０３ ６／ｉ ｎ１ ０．３ ５２ １． ． ８ Ｉ １ ． ９ ．ｓ． ６１５ ．０ １ ５０ ９ ｊｓ ０ ０ ０
Ｓ ｍｍ ａｙｏ ｂ ｏ ａ ｏ ｓｉ ｏ ｅｉ ｇｏ ｍ ｐｅ ｄ ｒ ｔｒ ｕ ｒ ｎＶｉｒ — ｃ ｕ ｔ Ｍ ｄ ｌ ｆ ｃ ｎ Ｃｏ ｌｘＵｎ ｅ ｗａｅ
２ ｌ 年 １ 月 ０１ ０
噪
声
与
振
动
控
制 பைடு நூலகம்
第５ 期
文 章编 号 ：０ ６ １ ５ （０ １０ ．０ ２０ １０ ．３ ５２ １）５０ ３ ．５
水下复杂圆柱壳体结构流 固耦合效应研 究概述
侯 磊 黄修长 徐时吟 ， ，
（． １海军驻 武汉 ４８ ３ 厂军代表室，武汉 ４０６ ； ３０ ０ ２ 上海 交通 大学 机 械 系统与振 动 国家重 点实验 室 ，上 海 ２０４ ． ０２ ０）
Ａｂ ｔａ ｔ： ｄ ｔ ｉ ｄ ｌ ｅ ａｕ ｅｒ ｖ ｅ ｏ ｉｒ — ｃ ｕ ｔ ｄ ｌ ｇ ｏ ｏ ｌｘ ｕ ｄ ｒ ｔｒｃ ｌ ｄ ｉａ ｈ ｌ ｓｒ ｃｕ ｅ ｉ ｓｒ ｃ Ａ ｅａ ｌ ｉ ｒ ｔｒ ｅ ｉｗ ｎ ｖ ｂ ｏ ａ ｏ ｓｉ ｍｏ ｅｉ ｆｃ ｍｐ ｅ ｎ ｅ ｗａｅ ｙ ｉ ｒｃ ｌ ｅ ｌ ｔｕ ｔ ｒ ｅ ｔ ｃ ｎ ｎ ｓ ｓ
ｄｕ ｌ ｌ ｅ ｃｌ ｄｉ ｌ ｈｌ ｒ ｂ ｄｃｌ ｄｉａ ｓｅ ， ｎ ｅｃ ｌ ｄｉａ ｓｅ ｉ ｔｒａ ｂ ｓｓ ｆ ｏｓ ｆｕ ｄｔ ｎ ｔ） ｏ ｂｅ ａ ｒ ｙｉ ｒ ａ ｓｅｌ ｉ ｅ ｙｉ ｒ ｌ ｈ ｌ ａｄｔ ｙｉ ｒ ｌｈ ｌｗ ｔ ｉ ｅｎｌ ａｅ （ｏ ｒ，ｏｎ ａｉ ｓ ｃ ｙ ｎ ｃ ，ｂ ｎ ｃ ｌ ｈ ｎ ｃ ｌ ｈｎ ｌ ｏ ｅ
ｉ ｔ ｅｅ ｎ ｈ ｄ ａ ｄ ｂ ｕ ａ ｙ ｅｅ ｎ ｈ ｄ ｆ ｉ ｌｍｅ ｔ ｔ ｏ ｎ ｏ ｄ ｒ ｌｍｅ ｔ ｔ ｏ ｎ ｅ ｍｅ ｎ ｍｅ
水 下 圆柱 壳体 结 构 是水 下潜 艇 ， 鱼雷 等 的基 本 结 构 形式 ， 振 动 噪声 问题 受 到广 泛 关注 。水 下 圆 其
壳体理论。有限长壳体和无限长壳体的一个根本区 别 是 有 限长 壳 体按 照模 态 展 开 时， 速 分布 和声 场 振
分 布 在边 界 处 无法 一一 对 应 ， 引起 不 同模 态 振速 之
间 的相互 耦 合 ， 体现 在 辐射 阻抗 中包 含互 辐 射 阻抗
水下复 杂 圆柱壳体 结构 流 固耦合效 应研 究概 述 和 自辐 射 阻 抗 。无 限长 圆柱 壳体 都 采 用直 接 展 开 ， 将 流 体加 载 效应 考 虑 为声 辐 射 阻抗 的形式 ， 声辐 射
ａ ｄ ｃ ｍｐ ｅ ｉｒ ｔ ｎ ｉｏ ａｉ ｎ ｓ ｓ ｍ ｒ ｒ ｓ ｎ ｅ ． ｅ ｏ ｓｆ ｒｍｏ ｅｉ ｇ ｄ ｎ ｍｉ ｅ ｐ ｎ ｅ ｆｆ ｉ ・ ａ ｅ ｙ ｉ ｄ ｃ ｌ ｎ ｏ ｌ ｘ ｖ ｂ ａ ｉ ｓ ｌｔ ｙ ｔ ａ ｅ ｐ ｅ ｅ ｔｄ Ｍ ｔ ｄ ｏ ｄ ｌ ｙ ａ ｃ ｒ ｓ ｏ ｓ ｓ ｏ ｕ ｄ ｌ ｄ ｄ ｃ ｌ ｒ ａ ｏ ｏ ｅ ｈ ｎ ｌ ｏ ｎｉ ｓ ｅｌ ｉ ｏ ｆｅ ｕ ｎ ｙ ｈ ｌ ｎ ｌ ｗ ｒ ｑ ｅ ｃ ，ｍｉ — ｅ ｕ ｎ ｙ ａ ｄ ｈ ｇ ｅ ｕ ｎ ｙ ｒ ｎ ｅ ｒ ｉ ｅ ．Ｓ ｍｅ ｐ ｔｎ ｉ ｌｍｏ ｅｉ ｇ ｔ ｃ ｎ ｑ ｅ ｒ ｄ ｆ ｑ ｅ ｃ ｎ ｉｈ ｆ ｑ ｅ ｃ ａ ｇ ｓ ａｅ ｇｖ ｎ ｏ ｏｅ ｔ ｄ ｌ ｅ ｈ ｉｕ ｓ ａｅ ｒ ｒ ａ ｎ ｐｏ ｉｅ ． ｒ ｖ ｄ ｄ
进行 了讨 论 。 对 内部带有 机 械设 备 以及基座 的壳体 的研 究更 加符 合潜 艇等 结构 的 实际情 况 。研究 内部 基座 的第
一
壳 。Ｍａ ｉ Ｇ ｎ ｕ 基 于 周 向导 纳 综 合 （ Ａ ｘｔ 和 ｉｏ ｘ Ｃ Ａ）
种 是 子 结 构 方 法 ， ｕ 研 究 了一 个 内部 具 有 … Ｇ ｏｕ
摘 要： 针对水下复杂圆柱壳体结构流固耦合效应 的研究对象， 从无限长， 限长 ， 有 双层 圆柱壳体， 加筋及有内部基 座和 内部复杂隔振系统等角度加以阐述 ； 针对建模方法 ， 介绍低频， 中频段和高频段的建模方法， 并介绍若干极具潜力 的复杂壳体流 固耦合建模方法。 关键词： 振动与波； 水下复杂圆柱壳体结构 ； 流固耦合效应 ； 耦合 的有限元和边界元
子结构 方法 ， 考虑 了周 向存 在不 均匀 加 强筋 ， 隔舱壁
等横 向加 强结 构 的水 下 壳体 的 动 力学 特 性 ， 方法 该 对 水下 流体加 载 壳体和 加 强筋子 结构 的导纳可 以通 过有 限元建模 或 解析法 建模 获得 。该方法 对 于旋转 对 称 结构 非 常实 用 ， 不 太适 合 于对 纵 向加 强 结构 但 的建模 。从 力学 的角度 ， 肋骨 的作 用体 现 为 环肋 和 舱 壁与 圆柱 壳之 间的径 向力 ， 向、 向力 以及 纵 向 切 纵 弯矩 。加 强筋 导致波 数转 换效 应 ， 将非辐 射 、 亚音 速 的高 波 数 成 分 转 换 成 有 辐 射 、 音 速 的 低 波 数 成 超
阻 抗 的形 式表 示 ； 四种是 将 环肋 看 作 较 窄 的 圆柱 第
解， 然后 引入 壳体 振 动方 程 ， 用 傅 氏变 换 ， 到 声 采 得
压 ， 功率 和辐射 效 率 。这些研 究 结果 表 明， 声 双层 圆 柱 壳体 的 能量传 递 ， 在低频 以声介质 耦 合为 主 ， 频 高 段 以实肋 板 为主 。殷 学 文 ｎ利用 Ｋｅｎ ｒ 壳体 理 论 ｎ ａｄ 研 究 了通 过 周期环 板 和壳 间流 体耦 合 的双层 壳体在 集 中力作 用 下 的振 动 声辐射 特性 ， 并对 壳间流体 、 环 板 的载荷 特 征进行 了数 值分 析及 其对 声辐 射 的影 响
ｇ ｅ ，ｎ ｏｈａｐ ｃ ｆｅｅｒｈｏｊｃ ｎ ｄｌ ｇａ ｐｏｃ ｅ． ｅｉ ｎｉｆ ｉ ｌ ｎ ｎ ｎｔｙｃｌ ｄｉａ ｓｅ ， ｉ ｎ ｉ ｔ ｓｅｔｏ ｓａｃ ｂｅｔ ａｄｍｏ ｅ ｎ ｐ ｒａｈ ｓＲ ｖｅ ｏ ｉ ｔｙｌ ｇａｄｆ ｉ ｌ ｉ ｒ ｌｈ ｌ ｖ ｂ ｓ ｒ ｓ ｉ ｗｓ ｎ ｎ ｅ ｏ ｉ ｅ ｙｎ ｃ ｌ