流固耦合力学的研究现状

流固耦合力学的研究现状

流固耦合力学是一个相对新的力学边界分布，它是由力学和固体力学的相交引起的。它的研究对象是固体在流场作用下的各种行为以及固体变形或运动对流场的影响。流固耦合力学的重要特征是两相介质之间的相互作用:固体在流体动载荷作用下产生变形或运动,而固体的变形或运动又反过来影响到流场,从而改变流体载荷的分布和大小。总体上,流固耦合问题按耦合机理可分为两大类:一类的特征是流固耦合作用仅仅发生在流、固两相交界面上,在方程上耦合是由两相耦合面的平衡及协调关系引入的;另一类的特征是流、固两相部分或全部重叠在一起,耦合效用通过描述问题的微分方程来实现。本文就流固耦合问题的两大分类中三种基本情况进行了讨论。1洋以及生物领域的意义流固耦合作用的研究在航空、航天、水利、建筑、石油、化工、海洋以及生物领域都有着十分重要的意义。如液体晃动对火箭飞行稳定性的影响,大型贮液管在地震激励作用下产生的流固耦合作用,液体湍振对输液管道的影响。本文就如下三个大方面进行了总结。1.1输流管道非线性分析模型的研究现状流体引起输流管道振动的研究最初来源于横跨阿拉伯输油管道振动的分析。管道在众多的工业领域中应用十分广泛,作用极其重要。但是,在管道内流体流动状态的微弱变化往往引起在工作过程中的湍振现象,诱发流体、管道之间的耦合振动,动力学行为相当复杂。这使得人们很早就开始了这方面的研究,PaidoussisMP是其中最具有代表性的。输流管道的振动问题之所以能引起学者的兴趣,除因为该问题的广泛工业背景和现实意义之外,还因为输流管道虽然是最简单的流固耦合系统,但它却涉及了流固耦合的大多数问题,并且它的物理模型简单,系统比较容易实现,因而便于理论与试验的相互协同。考虑因素侧重面的不同,输液管道非线性运动方程有几种类型,它们之间有一定的差别。它们的基本假设都是:流体无粘且不可压;管道作为梁模型来处理;管道只是在平面内振动。尽管输流管道的非线性动力问题受到50多年极为广泛的研究,但至今尚没有一个公认的模型。文建立的4个独立变量(轴向位移、横向位移、流速和压力)的全耦合模型(耦合形式包含摩擦耦合、Poisson耦合、结合部耦合以及管道轴向和横向运动的耦合)在众多的非线性分析模型中是一个较为完整的模型。m¨u+mf[˙υf(1+u′)+2υf˙u′+υ2fu″+˙ωυ′f]+˙Ρ(υf+˙u)/c2F-[(1-2υ)Ρ(1+u′)]′+4fρfρ′⋅υ2f/DΚ-gmf(1-2υ)(1+u′)ω′-EΙ(7ω″ω‴+ω′ω‴)-EAp(2u″+6u′u″+2ω′ω″)/2=0(1)m¨ω+mf[˙υf(1+ω′)+2υf˙ω′+υ2fu″+ω″υ2f]+˙Ρ(υf+˙ω)/c2F-[(1-2υ)Ρω′]′-gm+EΙω″″-EΙ(u‴′ω′+6u″ω‴+4u′ω‴′)-EAp(u″ω′+u′ω″)=0(2)˙Ρ/c2F+mf[(1-2υ)(˙u+υf)u″-˙u′+υ′f]-mf(1-2υ)(˙u′+u′˙u′+ω′˙ω′)=0(3)Ρ′+mf(¨u+˙υf)+mf¨ωω′+gmfω′+Dfρfυ2f/2=0(4)随着对输流管道问题研究的深入,各种不同的分析计算方法也相继被提出。其中有限元法(FEM)是应用最为广泛的一种数值方法。文用FEM法求得了载流在流速超过临界流速之后的振幅、主频率和振形,所得结果和试验结果吻合较好。传递矩阵法(TMM)在振动工程中一直是一种行之有效的方法。与FEM相比,其传递矩阵维数不随单元的增加而增加,恒保持为7维,这正是它的优越性所在。最早用于研究桥梁与车辆之间相互作用的结构阻尼法也被采用于解决输流管道的动态响应问题。矩阵摄动法可以用于求解输流管道系统动态方程度特征值问题。该方法仅对基本解进行修正而不需要重新求解,可以大大地节省计算时间,不需要重新编制软件,可以利用现有的有限元程序计算基本解,从而可以提高计算效率。混沌问题一直被认为强非线性运动所具有的一种力学现象,一般认为输流管道的运动是弱非线性的,因此对混沌的研究比较少。文是研究输流管道混沌现象比较早的,他们采用的分析模型是在管自由端的左、右两侧放置两块永久磁铁,以便施加一种强非线性的磁力。计算和实验都表明:当作用的磁力使输液管产生屈曲(弯曲状态)后,加大流速使其超过颤振对应的临界值时,混沌摆动就会发生。研究趋势:(1)针对输流管道系统运动方程的特殊性以及工业管路复杂性的实际问题,需要提出行之有效的计算方法和简化方法;建立模型求解以达到理论指导实践和验证结果真实性的目的;(2)外激励作用下管道的非线性瞬态响应,以及两者共同作用下引起输液管分叉的特性以及什么样形式的外激励容易诱发分叉和混沌现象,目前所知甚少,是今后发展研究的难点。1.2基于流固耦合的数值分析方法1964年,美国Alaska地震引起众多石化容器在地震载荷下惨遭破坏,使国民经济受到极大损失。由于化工工业在现代工业中的地位,促使部分科技工作者对含液容器进行研究。早期的许多俄罗斯专家在这一问题上作出了开创性研究。文中,详述了含液容器震荡问题的有关基本理论。由于计算工具的落后,当时的研究主要局限在刚性不动容器或作简单运动刚体中液体的震荡。对含液容器的动力特性进行分析,是为了了解容器内流体的晃动特性以及流体对容器固有特性的影响,以便为进一步研究含液容器的动力响应和动力稳定性奠定基础。当前国际上在有关含液容器的动力特性研究方面,研究手段多以实验为主,计算则多采用现有的一些通用有限元软件进行,但这些软件通常都是用于结构分析的,在流固耦合分析方面,尤其是当流体具有较大自由表面,且考虑可压缩性时,就很难用这些通用软件进行精确的分析。由于流固耦合系统的复杂性,其求解主要立足于数值分析。起初人们自然想到的是用位移法结构分析的通用程序来求解耦合问题,不同的是只要将流体视为剪切刚度为零的固体即可。但实际计算发现,剪切刚度为零,计算中出现零能模式,方法无法推广应用。至于流体中采用压力、固体中采用位移的混合模式没有零能模态的困难,但其有限元方程中的系数矩阵是非对称的。经典的流体中采用压力、固体中采用位移的流固耦合有限元方程:[ΜS0-QΤΜf]{¨a¨p}+[ΚS1ρfQ0Κf]{ap}={FS0}(5)式中,p为流体节点压力向量,a为结构节点位移向量,Q为流固耦合矩阵,Mf和Kf分别为流体质量阵和刚度阵,Ms和Ks分别为结构质量阵和刚度阵,FS为结构载荷向量,ρf为流体密度。文对上述方程进行对称化处理,得到:[ΜSEEΤˉΜf]{¨a¨p}+[ˉΚS00ˉΚf]{ap}={00}(6)引入的旋转周期性分析方法,可大大缩短计算规模和计算时间,符合含义同(5)式。含液容器流固耦合问题经过近几十年的发展,已经建立了几种行之有效的数值方法。Amsden在20世纪60年代发展起来的MAC(marker-and-cell)方法即标记子与单元方法,该方法具有不同寻常的灵活性,但计算存储量大,很难推广到三维问题的计算。VOF(volumeoffluid)方法改进了MAC方法的存储量大以及计算量大的缺点,特别适于三维计算,是发展较为成熟的方法。但VOF方法的每一时间步都需要计算流体体积函数F的值,其重复计算量也比较大。MAC方法和VOF方法实际上都应用了有限差分法,在处理不规则边界时存在某些缺陷。FEM法的最大优点是比较容易处理各种复杂的几何形状和边界条件,具有丰富的数学结构,可以达到非常高的计算精度。但它也存在一些困难,如计算对流算子的非对称性。边界元(BEM)方法是在综合FEM和经典的边界积分方程基础上发展起来的,它的最大优点是可将求解空间降一维。在上世纪五、六十年代宇航工程的发展过程中,液体运载器推进剂的晃动问题变得十分突出。这是因为液体晃动所产生的低频横振导致火箭飞行的不稳定性。因此对液体晃动特性的研究一直是宇航中一个极为重要的课题。为保障火箭飞行的稳定性,必须对液体的晃动加以抑制。采用的防晃措施最主要和最有效的手段是使用液体晃动阻尼器,常见液体晃动阻尼器是防晃板防晃。自20世纪60年代至今我国一直在对半圆形挡板的防晃特性进行多方面的试验研究和分析工作,获得了许多成果。研究趋势:(1)含液容器的大幅晃动的数值模拟方法研究;(2)有限元并行计算技术研究;(3)非线性耦合系统的稳定性、分岔、混沌等问题的研究;(4)防晃板防晃特性的研究。1.3地下储层流体的固结机理地下储层是固体、液体和气体的多孔介质组合。经典渗流力学已对固体孔隙中的流体流动问题进行了深入的研究,但它并没有涉及到流—固耦合问题,而在油气开采、地下水抽放等过程中,由于多孔介质流体压力的变化,一方面要引起固体特性的变化;另一方面,这些变化又反过来影响多孔介质中流体的流动和压力的分布。因此,在许多情况下必须考虑流体,包括液体(油或水)、气体(天然气、煤矿瓦斯等)与多孔介质的相互作用力的问题,即应考虑多孔介质应力场与多孔介质中渗流场之间的相互耦合作用。文是最早研究地下储层流体—固体耦合问题。它将可变形、饱和的多孔介质中流体的流动作为流动—变形的耦合问题来看待,建立了一维固结模型,获得了广泛应用。文则进一步研究了三向变形材料与孔隙压力的相互作用,并在一些假设条件的基础上,建立了比较完善的三维固结理论。自此,地下储层流固耦合理论的发展主要围绕着假设不同孔隙材料的模式而得到不同的物理方程:孔隙流体假设为不可压缩的和可压缩的等。地下储层流固耦合问题的显著特点是固体区域与流体区域互相包含、互相缠绕,难以明显地划分开,属于流固耦合问题的第二大分类,因此必须将流体相与固体相视为相互重叠在一起的连续介质,在不同相的连续介质之间可以发生相互作用。这个特点使得流固耦合问题的控制方程需要针对具体的物理现象来建立,而流固耦合作用也是通过控制方程反映出来的。即在描述流体运动的控制方程中有体现固体变形影响的项,而在描述固体运动或平衡的控制方程中有体现流体流动影响的项。储层流固耦合的(单相流体流动,在地下储层可能存在油、水、气等多相,多相问题不在此详述)岩石平衡方程:∫ΩδεΤDΤ∂ε∂tdΩ-∫ΩδεΤm∂p∂tdΩ+∫ΩδεΤDΤ∂p∂t13ΚsdΩ-∫ΩδεΤDΤ∂ε0∂tdΩ-∂f∂t=0(7)流体在可变形储层中流动的连续性方程:-ΔΤ{ΚμΔ(p+ρgh)}+[mΤ-mΤDΤ3Κs]∂ε∂t[1-φΚs+φΚf-1(3Κs)2mΤDΤm]∂p∂t=0(8)平衡方程(7)和连续性方程(8)即为弹性储层中单相流体流动的流固耦合方程,符合具体含义参见文。我国经过几十年的发展,已经在弹性储层领域取得相当成果。2计算机处理技术描述流固耦合问题的最根本的方法是数学方法。深入研究流固耦合问题就必须研究和学习掌握近代经典动力学的各种数学方法。将数学上的最新研究成果应用于流固耦合系统的研究中,从机理与本质上探讨流固耦合问题,促进流固耦合问题研究的进展。尽管流固耦合问题经过几十年的发展已经比较完善,但要真正解决好流固耦合问题必须结合流体力学与固体结构分析中的各种方法与细节。流固耦合问题是强非线性的问题,求解方程的规模加大,还必须