声振耦合与流固耦合
ansys多物理场耦合技术和方法
ANSYS是一种广泛应用于工程领域的仿真软件，它提供了多物理场耦合分析的能力，用于模拟和解决多个物理现象相互作用的问题。
以下是ANSYS多物理场耦合技术和方法的一些常见应用：1. 结构-热耦合（Thermo-Structural Coupling）：这种耦合方法用于分析结构在热载荷下的变形和应力响应。
它可以考虑热传导、热辐射、温度梯度等对结构性能的影响，并通过结构和热传导方程之间的相互作用来解决这些问题。
2. 结构-电磁耦合（Electromagnetic-Structural Coupling）：这种耦合方法用于研究结构在电磁场作用下的响应。
它可以考虑电磁场的电流、磁场、电磁感应等对结构的影响，并通过结构和电磁场方程之间的相互作用来解决这些问题。
3. 流体-结构耦合（Fluid-Structure Interaction, FSI）：这种耦合方法用于模拟流体和结构之间的相互作用。
它可以考虑流体力学中的压力、速度、湍流、流体-固体界面等对结构的影响，以及结构对流体的阻力、振动等反馈作用。
4. 流体-热耦合（Fluid-Thermal Coupling）：这种耦合方法用于模拟流体和热传导之间的相互作用。
它可以考虑流体在流动过程中的热对流、辐射等对热传导的影响，以及热传导对流体温度分布的影响。
5. 电磁-热耦合（Electromagnetic-Thermal Coupling）：这种耦合方法用于模拟电磁场和热传导之间的相互作用。
它可以考虑电磁能量的吸收、热产生和热扩散等对系统温度分布的影响，以及温度对电磁特性的影响。
以上只是ANSYS多物理场耦合技术和方法的一些例子，实际中还有其他类型的耦合分析，如声-结构耦合、声-流体耦合等。
通过使用这些耦合技术和方法，工程师可以更准确地模拟和分析不同物理场之间的相互作用，从而更好地优化设计和解决实际问题。
comsol多物理场耦合
流-固耦合接口
• 在流-固耦合接口中:
• 结构外部载荷是由流体力提供的 • 移动壁条件由结构速度场来设定 • 网格位移即为结构位移
• COMSOL Multiphysics自动计算壁上的流体力
- 预定义的表达式适用于总力
• 流-固耦合接口是三个物理接口的组合:
• 结构力学 • 层流 • 移动网格
网格剖分
• 规则几何结构，使用映射和扫掠网格剖分
弹性波地层传播
声固耦合接口
• 声固耦合 • 声壳耦合 • 声压电耦合
预定义耦合 接口
• 弹性波 • 多孔介质弹性波
• 多孔介质与其中饱和液的相互作用 • 基于BBiblioteka ot理论 • 快、慢纵波 • 横波
15.98KHz
声-固耦合
声-壳耦合
案例：声-壳耦合
弹性波传播问题
参数与几何
超弹性材料模型
• 右键固体力学，添加超弹性材料模型 • 在超弹性材料模型中通过下拉菜单选择Murnaghan
• 在材料定义节点中 定义其材料属性
对称边界
• 将图中所示的三个面定义为对称边界
描述位移边界
• 在左侧端面施加描述位 移边界条件，给出其在x 方向的位移表达式。
3D 稳态 FSI
• 几何: 拉伸 1.5e-5
3D 稳态 FSI
• 入口流速: 1[m/s]
3D 稳态 FSI
• 网格剖分: 物理控制, 单元尺寸- Extra Coarse
小技巧
• 使用带有平滑过渡的阶跃函数来开始瞬态分析
• 在前期尽量避免使用瞬态分析，以保证达到一致的初始条件
• 使用描述网格位移条件来降低移动网格（ALE）中 的低质量/反转网格
复合材料螺旋桨流固耦合数值计算
华中科技大学硕士学位论文摘要作为常见的船舶推进装置，螺旋桨具有几何形状简单、安装使用方便、推进效率高等特点。
但随着船舶事业的飞速发展，船舶逐渐向大型化、高速化的特点发展，对船舶螺旋桨性能的要求越来越高。
现阶段的常规材料螺旋桨具有噪声大、不耐腐蚀、使用寿命较短等缺点，因此近年来螺旋桨材料发生了巨大的变化。
由于新型材料的物理性质发生了变化，螺旋桨的流固耦合问题引起了越来越多的关注。
本文主要对DTMB 4119螺旋桨流固耦合问题进行以下方面的研究工作：(1)针对DTMB4119螺旋桨在不同工作工况下，研究流体计算网格量对螺旋桨敞水性能的影响；采用标准k-e、RNG k-e、SST k-w湍流模型，研究分析湍流模型对螺旋桨敞水性能的影响，并与参考文献给出的试验结果进行对比分析，从结果来看采用RNG k-e湍流模型的计算结果与文献误差最小。
(2)采用三种材料属性不同的各向同性材料，对螺旋桨水动力性能与结构响应特性与材料属性之间的关系进行探究。
从仿真结果来看，材料属性对前两种材料螺旋桨影响较小，流体与固体之间的相互影响可以忽略不计；而对玻璃纤维材料螺旋桨的影响较大，纤维材料螺旋桨的敞水性能、流场特性以及桨叶的应力应变、变形都有明显的变化。
(3)计算考虑铺层方式的流固耦合特性。
结果表明，层铺方式对螺旋桨水动力性能影响较小，但对螺旋桨的应力应变影响较大。
关键词：流固耦合，复合材料螺旋桨，水动力性能，结构响应特性，铺层方式华中科技大学硕士学位论文AbstractPropeller, as a common propulsion device, has the characteristics of simple structure, convenient use and high efficiency. But with the rapid development of the shipbuilding industry, the ship develops towards the characteristics of large scale and high speed, which also requires more. At present, the conventional propeller has many disadvantages, such as large noise, no corrosion resistance, short service life and so on, so many new type propellers have emerged. And due to the change of material properties, the fluid - structure interaction of propellers has attracted more and more attention.In this paper, the following two aspects are carried out on the two ways fluid- structure interaction of propeller:(1) For the DTMB4119 propeller under different working conditions, using the standard k-e, RNG k-e and SST K-W turbulence model, the open water performance of the propeller is calculated by numerical simulation. The results are compared with the experimental results given in the reference literature. The results show that the error between the RNG k-e turbulence model and the literature are minimum.(2) Three isotropic materials of alloy steel, nickel aluminum bronze and glass fiber are used to investigate the relationship between the hydrodynamic performance of a propeller, the response characteristic of the structure and the properties of the material. From the simulation results, the material properties have little influence on the first two kinds of material propellers, and the interaction between the fluid and the solid can be ignored, but the effect on the propeller of glass fiber material is larger. The open water performance, the flow field characteristics and the stress strain and deformation of the blade have obvious changes.(3), Calculating and comparing the fluid solid coupling characteristics of Considering the effect of composite stacking mode on composite propeller. The results show that the stacking mode method has little effect on the hydrodynamic performance of propeller, but it has great influence on the stress and strain of propeller.Key words:FSI, composite propeller, open water performance, structure response stacking mode华中科技大学硕士学位论文目录摘要 (I)Abstract (II)1 绪论1.1课题研究背景及意义 (1)1.2复合材料螺旋桨的应用与研究现状 (2)1.3论文主要研究内容 (7)2 数值计算基本理论2.1引言 (8)2.2计算流体力学基本原理 (8)2.3结构计算相关理论 (11)2.5本章小结 (15)3 螺旋桨水动力性能数值方法3.1引言 (16)3.2螺旋桨水动力计算 (16)3.3网格独立性验证 (21)3.4湍流模型的选取 (24)3.5螺旋桨敞水性能分析与验证 (26)3.6本章小结 (27)华中科技大学硕士学位论文4 螺旋桨流固耦合特性模拟与分析4.1引言 (28)4.2螺旋桨结构计算方法 (28)4.3螺旋桨流固耦合计算 (29)4.4本章小结 (38)5 复合材料螺旋桨流固耦合计算5.1引言 (39)5.2基于ACP复合材料分层有限元模型 (39)5.3复合材料螺旋桨流固耦合结果分析 (43)5.4本章小结 (45)6 结论与展望6.1结论 (46)6.2展望 (47)致谢 (48)参考文献 (50)华中科技大学硕士学位论文1 绪论1.1课题研究背景及意义进入新世纪后，随着全球经济的飞速发展，越来越频繁的贸易往来使得全球交通运输业发生了巨大变化。
管道流固耦合振动及声传播的研究现状及展望
第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]。
对湍流噪声的机理已有许多研究。
(整理)耦合分析
耦合分析耦合场分析的定义耦合场分析是指在有限元分析的过程中考虑了两种或者多种工程学科（物理场）的交叉作用和相互影响（耦合）。
例如压电分析考虑了结构和电场的相互作用：它主要解决由于所施加的位移载荷引起的电压分布问题，反之亦然。
其他的耦合场分析还有热-应力耦合分析，热-电耦合分析，流体-结构耦合分析，磁-热耦合分析和磁-结构耦合分析等等。
耦合场分析的类型耦合场分析的过程取决于所需解决的问题是由哪些场的耦合作用，但是，耦合场的分析最终可归结为两种不同的方法：序贯耦合方法和直接耦合方法。
序贯耦合解法序贯耦合解法是按照顺序进行两次或更多次的相关场分析。
它是通过把第一次场分析的结果作为第二次场分析的载荷来实现两种场的耦合的。
例如序贯热-应力耦合分析是将热分析得到的节点温度作为“体力”载荷施加在后序的应力分析中来实现耦合的。
直接耦合解法直接耦合解法利用包含所有必须自由度的耦合单元类型，仅仅通过一次求解就能得出耦合场分析结果。
在这种情形下，耦合是通过计算包含所有必须项的单元矩阵或单元载荷向量来实现的。
例如利用单元SOLID5，PLANE13，或SOLID98可直接进行压电分析。
流-固耦合分析主要用于解决流体(含气体)与结构之间的相互作用效应。
包括: 流-固耦合法、水弹性流体单元法、虚质量法。
（1）流-固耦合法流-固耦合法广泛用于声学和噪音控制领域中,如发动机噪声控制、汽车车厢和飞机客舱内的声场分布控制和研究等。
分析过程中,利用直接法和模态法进行动力响应分析。
流体假设是无旋的和可压缩的, 分析的基本控制方程是三维波方程, 二种特殊的单元可被用来描述流-固耦合边界。
(噪)声学载荷由节点的压力来描述, 其可以是常量, 也可以是与频率或时间相关的函数, 还可以是声流容积、通量、流率或功率谱密度函数。
由不同的结构件产品的噪声影响结果可被分别输出。
（2）水弹性流体单元法该方法通常用来求解具有结构界面、可压缩性及重力效应的广泛流体问题。
流固耦合声子晶体管路冲击振动特性研究
流固耦合声子晶体管路冲击振动特性研究*胡兵1)2) 郁殿龙1)2)† 刘江伟1)2) 朱付磊1)2) 张振方1)2)1) (国防科技大学装备综合保障技术重点实验室, 长沙 410073)2) (国防科技大学智能科学学院, 长沙 410073)(2020 年3 月19日收到; 2020 年6 月12日收到修改稿)流固耦合管路系统广泛应用于各种装备中, 通常用来传递物质和能量或者动量. 由于流固耦合效应, 管壁在流体作用下易产生强烈的振动与噪声, 对装备安全性、隐蔽性产生严重影响, 甚至造成严重破坏. 流固耦合管路振动抑制需求迫切, 意义重大. 声子晶体可以利用其带隙特性抑制特定频率范围内弹性波的传播, 在减振降噪领域具有广泛的应用前景. 本文基于声子晶体理论, 研究了流固耦合条件下的布拉格声子晶体管路冲击振动传递特性. 将传递矩阵法和有限元法相结合, 计算了能带结构与带隙特性, 重点考虑了流固耦合效应下, 不同冲击激励条件下声子晶体管路振动特性, 分析了流固耦合对声子晶体管路振动传递特性的影响.研究结果为流固耦合条件下管路系统的振动控制提供了技术参考.关键词：流固耦合, 声子晶体, 振动带隙, 冲击振动PACS：43.40.–r, 61.50.Ah, 47.11.–j, 65.20.De DOI: 10.7498/aps.69.202004141 引 言输流管路系统通常用来传递物质和能量或者动量, 因而广泛应用于各种装备中. 流体的压力波脉动和管壁结构容易产生耦合作用, 进而诱发强烈的振动与噪声, 对装备的性能与使用产生严重影响. 因此, 研究降低输流管路振动和保证管道输送安全, 在理论和实践上有着重大研究意义[1].流固耦合动力学是一门研究固液相互作用的学科, 其主要研究内容是变形固体在流体流场作用下的力学行为与变形固体形态对流场的影响之间的相互作用[2,3]. 随着计算固体力学和计算流体力学的快速发展以及各种商用有限元软件的开发使用, 流固耦合分析和研究得到快速发展, 研究结果对工程应用和装备设计起到越来越重要的参考价值[4−6].近年来, 凝聚态物理领域中声子晶体(phononic crystals)带隙(band gap)理论的不断发展和完善为振动传播控制提供了新的技术支持[7−11]. 声子晶体是某种或多种材料组成的周期性结构或复合材料. 弹性波在声子晶体内传播时, 受内部介质周期性的作用可以产生弹性波带隙, 因而可以利用声子晶体的带隙特性有效抑制带隙频率范围内的振动与噪声传播. Chen等[12−14]研究了嵌入式内部谐振器的夹层梁结构, 该结构可改善冲击载荷下弯曲振动性能, 并完成实验验证. Pai等[15]通过将阻尼元件结合到多谐振器超材料梁中, 实现两个带隙频率区域有效地合并以形成单个宽带能量吸收区域. Chen等[16]提出了一种由多层黏弹性连续介质构成的耗散超材料的微观结构设计, 可有效地衰减瞬态冲击波. Alamri等[17]介绍了具有多个Maxwell 型谐振器的耗散弹性超材料的发展, 可应用于减轻动态载荷和爆炸波衰减. Li等[18]提出了一种新颖的超晶格桁架芯夹层结构, 可用于实现脉冲波衰减和动态载荷衰减, 具有缓和冲击能力和动能吸收能* 国家自然科学基金(批准号: 11872371)和国家自然科学基金重大项目(批准号: 11991032, 11991034)资助的课题.† 通信作者. E-mail: dianlongyu@© 2020 中国物理学会 Chinese Physical Society 力. 李奇奇等[19,20]提出了一种三谐振器超材料, 以增强冲击应力波的衰减效果, 并对该超材料的多目标优化进行了分析. 而后提出了一种用于衰减冲击应力波的新型多谐振器超材料. 以上对声子晶体冲击波衰减的研究有很多的进展, 并且由理论研究逐渐进展到实际应用中, 但研究对象均为固体结构,而流固耦合声子晶体冲击振动特性研究较少.应用声子晶体带隙特性进行流固耦合管路系统设计, 可将管路系统设计成周期性复合结构或在管路上周期性地附加局域共振结构, 进而实现抑制管路系统振动传播, 为管路的减振降噪设计提供了新的技术途径和理论基础, 目前已经得到广泛研究和关注[21−23]. Koo等[24]首先确定了带有周期性弹性支撑的输液管路的带隙分布, 并通过实验很好地证明了其理论预测. Sorokin等[25]则研究了平面弹性波在充液周期壳体中的传播特性, 并分析了充液与否对带隙特性的影响, 随后在研究中发现在流固耦合情况下充液管路系统在某些频段同样存在一些“波阻带”现象[26], 还进一步研究了周期附加惯性质量系统管路的纵向振动与弯曲振动的耦合振动带隙特性[27]. 郁殿龙等[28]通过应用布拉格散射机理和局部共振机理, 实现了周期性复合管路输液的弯曲振动带隙, 并进一步应用传递矩阵法和有限元法研究了各种条件下输液管路的波传播和衰减特性[29−33], 并且进行了振动试验验证了周期性管路结构的波衰减能力[34]. 魏振东等[35]将周期管路应用到液压系统中, 提出了一种考虑流固耦合作用的频响计算方法, 并对其在高压条件下的带隙特性进行了理论和实验研究. 刘东彦等[36]首次研究了液压油流体特性对周期管路带隙特性的影响. 沈惠杰等[37,38]重点研究了周期性输液管路的壳体的稳定性并进一步提出了一种由功能梯度材料(FGM)制成的周期壳体模型[39−41], 以消除或减轻由材料参数的几何不连续性引起的应力集中. Liang等[42]基于声子晶体管路模型, 考虑管道长度变化, 研究了部署长管输送流体的波传播和带隙特性. 以上研究为声子晶体管路输送流体的振动分析奠定了基础. 但是这些研究中对流固耦合效应考虑得较少,特别是冲击激励的振动分析不够深入.本文以流固耦合声子晶体管路为对象, 考虑不同冲击激励条件对振动传递特性的影响. 首先采用传递矩阵法对未充液和充液周期管路的振动传递特性进行数值分析, 分析其带隙特性. 并利用有限元法, 研究不同冲击激励条件下(包括管壁冲击激励、流体冲击激励等)的流固耦合声子晶体管路振动传播特性及其影响规律.2 流固耦合理论与传递矩阵法2.1 充液声子晶体管路传递矩阵法充液管路振动模式有弯曲振动、轴向振动、扭转振动, 以及它们之间的复杂耦合振动, 其中弯曲振动指是图1中y方向的振动, 即垂直于管路轴线方向上的振动; 轴向振动指图1中x方向的振动,即沿管路轴线方向的振动; 扭转振动则是指管路在绕轴线进行的扭转振动, 其一般由旋转机械的主动力矩与负荷反力矩之间失去平衡引起的. 当内部流体以较高或较慢的速度流动时, 管路可能会弯曲或强烈振动, 而且在外部激励下, 管路主要会产生弯曲振动, 轴向振动和扭转振动较小, 可以忽略不计,故管路弯曲振动是主要的振动模式[43−45]. 因此, 研究弯曲振动对管路振动控制具有重要的理论意义. 目前研究管路弯曲振动的理论研究大多数基于梁模型.一般情况下, 当管路长度与管径长之比大于10时, 可以将管路考虑为欧拉梁模型. 为了计算能带结构, 这里我们对流固耦合管路进行简化, 假设管内的液体为理想液体(各向同性、均匀、不可压缩、线性), 液体流速以恒定的速度流动, 液体中未发生空泡现象, 忽略重力的影响, 充液管路弯曲振动Euler梁方程为[46,47]:(b)图 1 布拉格声子晶体管路结构示意图 (a)无限周期单元; (b)基本周期单元Fig. 1. Schematic diagram of Bragg phononic crystal pipeline structure: (a) Infinite periodic cell; (b) Basic peri-odic cell.m f m p 式中, E 为管路材料杨氏模量, I 为管路横截面转动惯量, w 为弯曲振动位移, p 为流体压力, A 为管路内横截面积, 为液体的单位长度质量, 为管路的单位长度质量, u 和t 为液体流速和时间.考虑管路未充液时, 振动方程(1)式可以简化为w (x,t )=W e kx e i ωt 对于一简谐波[48], (1)式的解可以写成 的形式, 则管路的弯曲振动方程可以表达为ωk 1,k 2,k 3,k 4对于给定的值 , (3)式中的波数k 有两个不同的虚数根和两个互为共轭的复数根, 分别记为, 则(1)式的解可以表示为x =na 构建的布拉格声子晶体管路见图1所示, 其中a 为晶格常数, l a为管路A 的长度, 则在单元n –1和单元n 之间的界面处的位移、转角、弯矩和剪切力都是连续的, 即在 处有[49]:W =[W 1,W 2,W 3,W 4]T 式中.x =na +l a 同理, 在单元n 中的管道A 和管道B 之间的交界面处的位移、转角、弯矩和剪切力都是连续的,即在 处有:根据(5)式和(6)式, 可以求得单元n 与单元n –1之间的关系为T =H2−1K 2K 1−1H 1式中 .由于声子晶体管路在x 方向的周期性, 根据Bloch 定理可以得到:q 式中 是轴向一维Bloch 波矢.因此, 无限周期管结构特征值是行列式的根:I 4×4ωq q T 式中 为 单位矩阵. 对于给定的 值, (9)式给出 的对应值, 根据 是实部还是虚部, 相应的波分别通过管道传播或被衰减. 同理, 可以从传递矩阵 获得振动传递的频率响应曲线(frequencyresponse function, FRF)[50].2.2 流固耦合管路振动基本原理通过流体力学N-S 方程理论和固体力学小变形弹性理论联合推导出的求解流固耦合的基本方程是如今使用较为普遍的模型, 主要包括4-方程、8-方程、12-方程和14-方程模型等, 目前4-方程模型和8-方程模型在实际应用中得到较为广泛的使用. 通常以Budenkov [51]在讨论Pochhammer 方程时采用的方法为参照, 对流固耦合问题构建合适的数学模型. 将流固耦合运动分解为轴向和横向运动等, 而后分别对其讨论, 最后将两类方程进行合成求解流固管路耦合振动问题.流体流动过程中遵循质量、能量、动量三大基本物理守恒定律, 当流体中混合有其他成分时, 还要遵循组分守恒定律[4]; 固体结构部分的控制方程主要遵循牛顿第二定律[52]. 流固耦合方程遵循基本的物理守恒定律, 在流固耦合交界面上满足流体域与固体域两相之间应力、温度、位移、热流量等变量的守恒或相等[53].目前, 直接耦合式解法和分离式解法是用来解决流固耦合问题最主要的两种方法. 直接耦合算法可以同时对流固耦合控制方程进行求解, 不存在时间滞后, 理论上较为理想; 分离解法则不必对流固耦合控制方程进行求解, 只需在不同求解器或者同一求解器中按照预先设置的求解顺序分别对流体和固体控制方程进行求解, 得到固体域和流体域的计算结果, 而后通过设置好的流固耦合面进行两相数据的交互传输, 当此刻收敛达到要求时进行下一时刻的计算求解, 依次计算求解最终结果.分离解法与直接耦合求解相比, 缺点主要是具有在流固耦合面上能量不完全守恒和时间滞后问题; 其优点是可以较好地结合流体力学和固体力学的方法和程序, 对内存要求较低, 计算速度较快,目前被应用于大多数商用CAE 软件中[4]. 第4节中基于ANSYS Workbench 平台建立的双向流固耦合分析模型就是采用分离解法对流固耦合方程进行求解的.3 充液声子晶体管路能带结构与带隙特性分析基于声子晶体理论, 构建布拉格声子晶体管l a l b a =l a +l b 路, 其结构示意如图1所示, 图1 (a)是无限周期单元, 图1 (b)是基本周期单元. 布拉格声子晶体管路是两种不同管壁材料A 和B 沿x 轴交替周期排布形成的周期管路, 其中单个周期单元中管段A 的长度为 , 管段B 的长度为 , 则周期管路的晶格常数 , 管路的半径为R , 管壁厚度为d .l a l b 本文研究对象为海水管路系统, 其脉动源假设为一台六叶片的离心泵, 转速为2500 r/min, 则推导出的其叶频和次倍频分别为250 Hz 和500 Hz,由于管路内径R 和管壁厚度d 尺寸为项目应用要求, 因此通过改变晶格常数a 和管段A 长度l a 和管段B 长度l b 来改变管路带隙, 为了使计算得到的管路带隙满足叶频和次频的振动控制要求, 计算中, 管段A 和B 分别采用结构钢和环氧树脂, 其材料参数如表1所示, 取管段A 的长度 = 0.25 m,管段B 的长度为 = 0.25 m, 管路的内径R =0.01 m, 管壁厚度d = 0.001 m, 管内介质为水, 密度为1000 kg/m 3, 介质内的声速为1400 m/s. 运用2.1节中由(7)式—(9)式求解传递矩阵T , 并计算上述参数下无限周期单元的能带结构图和频率响应.表 1 管路材料参数Table 1. Pipeline material parameters.材料名称杨氏模量/GPa密度/kg·m –3泊松比结构钢20078500.3环氧树脂4.3511800.3672图2是利用传递矩阵法计算的未充液布拉格声子晶体管路弯曲振动的能带结构和频率响应曲线, 其中图2 (a)是波矢实部与频率的关系曲线,可以表现带隙的频率范围, 图2 (b)表示具有5个周期的声子晶体管路的频率响应曲线. 由图可知,在0—800 Hz 的频率范围内存在两个衰减带隙-60-50-40-30-20-100102030100200300400500600700800波矢 k /p S a -1(b)频率f /H z频率f /H z频率响应FRF/dB图 2 未充液布拉格声子晶体管路的带隙特性 (a)能带结构; (b)振动频率响应曲线Fig. 2. Band gap characteristics of the liquid-unfilled Bragg phononic crystal pipeline: (a) Band structure; (b) Flexural vibration FRF.-40-30-20-10010203040100200300400500600700800(b)频率f /H z频率f /H z波矢 k /p S a -1频率响应FRF/dB图 3 充液布拉格声子晶体管路的带隙特性 (a)能带结构; (b)振动频率响应曲线Fig. 3. Band gap characteristics of the liquid-filled Bragg phononic crystal pipeline: (a) Band structure; (b) Flexural vibration FRF.70—90 Hz 和280—690 Hz, 其中观察振动频率响应曲线可知, 有限周期声子晶体管路的振动传递损失曲线较好地对应着无限周期声子晶体管路的带隙频率范围; 第二带隙的衰减强度远大于第一带隙, 且衰减最大可以低至–60 dB.图3是充液布拉格周期管路弯曲振动的能带结构和振动频率响应曲线, 可以发现充液管路在0—800 Hz 内出现3个衰减带隙, 分别是40—65 Hz 、180—340 Hz 和485—735 Hz, 表明当管路充液时,布拉格周期管路的弯曲振动带隙会向低频移动. 同时可以看出, 充液管路的第二带隙和第三带隙分别包含叶频与次倍频, 可以对离心泵引起的管路振动起到较好地抑制作用.4 流固耦合声子晶体管路振动特性仿真分析4.1 有限元模型及算法验证v in v out FRF =20×log (v out /v in )未考虑流固耦合效应时, 基于ANSYS Work-bench 平台, 利用谐响应模块, 在激励端施加幅值为1的速度信号, 记为 , 在响应端拾取速度信号, 记为 , 则由公式 计算频率响应, 得到5个周期的声子晶体管路振动频率响应曲线, 与传递矩阵法计算的振动传递响应相对比, 得到图4. 图4是5个周期未充液和充液布拉格管路的弯曲振动传递损失曲线, 其中蓝色虚线和黑色实线分别是有限元法和传递矩阵法的计算结果. 由图可知, 传递矩阵法和有限元法的计算结果比较吻合, 这有力地证明了传递矩阵法的准确性和有效性.图5是不同频率处未充液和充液布拉格声子晶体管路的速度幅值图, 其中图5(a)表示未充液声子晶体管路在f 为250, 500, 750, 1000 Hz 处的位移幅值. 可以发现, f = 500 Hz 和f = 1000 Hz 处管路的弯曲振动在前两个周期已经得到很大的衰减, 后半段的管路几乎没有发生振动, 表明布拉格声子晶体管路对带隙内的振动可以较好地抑制,且振动衰减效果明显. f = 250 Hz 和f = 750 Hz 都是带隙外的频率点, 两个频率处的位移幅值远大于带隙内的位移幅值, 且整个管路都处于振动状态, 此时, 弯曲振动可有效传递到管末端. 图5(b)表示充液声子晶体管路在f 为125, 250, 400,600 Hz 处的位移幅值. 可以发现, f = 250 Hz 和f =600 Hz 处管路的弯曲振动在前两个周期已经得到很大的衰减, 后半段的管路几乎没有发生振动, 表明充液布拉格声子晶体管路对带隙内的振动可以较好地抑制, 且振动衰减效果明显. f = 125 Hz 和f = 400 Hz 都是带隙外的频率点, 两个频率处的位移幅值远大于带隙内的位移幅值, 且整个管路都处于振动状态, 此时, 弯曲振动可有效传递到管末端.与图4中声子晶体管路振动频率响应曲线中的带隙相对应.进一步考虑流固耦合条件下声子晶体管路振动传递特性. 基于ANSYS Workbench 平台对不同冲击激励下流固管路进行流固耦合分析, 建立了双向流固耦合分析基本流程如图6所示, 参照图1声子晶体管路参数建立5个周期的声子晶体管路模型如图7所示. 模型选为瞬态结构和流体力学模频率响应F R F /d B频率响应F R F /d B频率f /Hz频率f /Hz 图 4 未充液和充液布拉格声子晶体管路弯曲振动频率响应 (a)未充液管路; (b) 充液管路Fig. 4. Frequency response of flexural vibration of liquid-unfilled and liquid-filled Bragg phononic crystal pipeline: (a) liquid-un-filled pipe; (b) liquid-filled pipe.kg /(m ·s )块. 管路长为2.5 m, 管路内径为0.02 m, 管壁厚度为0.001 m, 弹性管壁, 忽略管道阻尼影响. 流体部分选择水, 不可压缩流体, 其密度为1000 kg/m 3,温度设置为25 ℃, 动力黏性系数选择默认值0.001003 . 由于ANSYS Workbench 计算双向流固耦合时, 流体通过流固耦合面对管路结构传递的只是湍流运动产生的流体雷诺应力、黏性应力和脉动切应力的流体作用力, 不包括流体质量, 故在仿真计算轻质流体时结果较为准确, 但当ρe 管内流体为重质流体时, 仿真结果会产生较大的误差, 故需要在管路上引入附加质量∆m = r πr p 2. 本文的流体为水, 属于重质流体, 因此在管路建模时需要设置管路材料的当量密度 , 即ρ0ρf r in r out 式中, 为管路材料密度, 为流体密度, 为管路内径, 为管路外径. 管路材料选择表1中的结构钢和环氧树脂, 则仿真分析中设置的管路材料结构钢和环氧树脂的当量密度分别为12612 kg/m 3和5941.9 kg/m 3.m /s 本例中主要设置流体和液体接触的管道内壁面为流固耦合面, 设置限制两端面X , Y 和Z 方向位移为支撑边界条件, 即为固定约束(固支结构).Fluent 界面里设置不同的入口速度, 单位为 ,出口压力设置为零; 系统耦合分析界面设置耦合时间为0.0512 s 和耦合时间步0.0001 s, 选择双向耦合; 初始条件, 视为光滑管壁.(b)(a)(1) 第一带隙外f =250 Hz(2) 第一带隙内f =500 Hz(3)第二带隙外f =750 Hz (4) 第二带隙内f =1000 Hz (1) 第一带隙外f =125 Hz(2) 第一带隙内f =250 Hz(3) 第二带隙外 f =400 Hz(4) 第二带隙外f =600 Hz图 5 未充液和充液声子晶体管路不同频率处的速度幅值 (a)未充液管路; (b) 充液管路Fig. 5. Displacement amplitude of liquid-unfilled and liquid-filled phononic crystal pipeline at different frequencies: (a) liquid-un-filled pipe; (b) liquid-filled pipe.图 6 ANSYS Workbench 系统耦合配置方式Fig. 6. Coupling configuration of ANSYS Workbench system.出水口处进水口处450.00225.00675.00900.00 (mm)图 7 ANSYS 中建立流固耦合管路模型Fig. 7. Establishment of fluid-structure interaction pipelinemodel in ANSYS.4.2 管壁冲击激励条件下流固耦合声子晶体管路振动特性分析假设距离管路左端进水口处受到冲击载荷, 冲击波的形式满足以下方程:F max t 0t d 式中 = 200 N, = 0.0005 s, = 0.0001 s.冲击时间为0.0005 s, 在管路进出水口处拾取响应信号. 分别仿真计算未充液和充液管路的弯曲振动响应, 充液管路仿真分析时采用流固耦合仿真分析, 设置流速为0.管壁冲击激励的脉冲响应如图8(a), 利用快速傅里叶变换得到其频域分布如图8(b), 可以看出,由于管壁冲击激励的作用, 模型中产生的应力波具有0—5000 Hz 的宽频带.图9是5个周期声子晶体管路未充液和充液时的冲击振动特性. 其中蓝色实线和青色虚线分别表示进水口处速度信号时域以及通过快速傅里叶变换得到的频域分布和出水口处速度信号时域及快速傅里叶变换得到的频域分布. 频域图中阴影部分表示出口处的速度峰值较入口处的速度峰值小,表明在该频率范围内声子晶体管路具有较好的衰减效果. 对比图9(a)和图9 (b), 可以发现当管路充液时, 管路振动幅值无论在时域还是频域上都有所衰减, 表明当管路充液时, 由于流体质量的影响,外部冲击引起的管路振动减弱. 综合分析, 未充液声子晶体管路在270—625 Hz 内具有较好的衰减作用, 充液声子晶体管路在175—332 Hz 和488—725 Hz 内具有较好的衰减作用, 与第3节传递矩阵计算的未充液管路带隙280—690 Hz 以及充液管路带隙180—340 Hz 和485—735 Hz 相吻力/N幅值/d B时间/10-3 s频率/Hz图 8 管壁冲击脉冲响应及通过快速傅里叶变换得到的冲击模拟频域 (a) 管壁冲击时域; (b) 管壁冲击频域Fig. 8. Pipe wall shock impulse response and shock simulation frequency domain obtained by fast Fourier transform: (a) Time do-main of wall impact; (b) Frequency domain of wall impact.时间/s时间/s频率/Hz 频率/Hz10(a)(b)速度/m m S s −1速度/m m S s −1|速度(f )||速度(f )|图 9 未充液与充液声子晶体管路冲击振动响应 (a)未充液管路; (b) 充液管路Fig. 9. Shock vibration response of liquid-unfilled and liquid-filled phononic crystal pipeline: (a) liquid-unfilled pipe; (b) liquid-filled pipe.合, 表明未充液和充液布拉格声子晶体管路对管路外部冲击具有较好的抑制作用.图10表示未充液和充液声子晶体管路不同时刻的速度幅值图, 图10(a)表示未充液声子晶体管路在t 为0.001, 0.0025, 0.005, 0.01 s 处的速度幅值, 图10(b)表示充液声子晶体管路在t 为0.001,0.0025, 0.005, 0.01 s 处的速度幅值. 可以发现, 无论声子晶体管路是否充液, 其出口处的振动较入口处滞后, 且振动幅度较小, 这是由于管路入口处的冲击响应沿着管壁传播需要一定时间. 约在0.01 s,出口处振动幅度达到最大, 与图9中时域图中出口处的振动响应曲线相一致.接下来考虑流速对管壁冲击响应的影响, 设置流速分别为0, 10 m/s. 仿真分析得到流固耦合声子晶体管路在不同流速下的冲击振动响应.图11表示流固耦合声子晶体管路在流速为0和10 m/s 下的冲击振动响应, 图12表示流固耦合声子晶体管路在流速为0和10 m/s 下的出口处的振动响应, 由图可知, 当流速较小时, 流速改变,布拉格周期管路的衰减频率范围基本不变. 虽然随着流速的增大, 管路的衰减效果减弱, 但影响不大,这是由于管路在外部冲击作用下, 管路的振动主要由外部冲击引起, 管路内流对管路振动作用较小,但由于流固耦合效应, 流体运动产生的管壁作用力分布在管路内部流固耦合面上且前端由流体引起的管壁振动会向管路末端传递, 导致管路出口段的振动响应增大, 从而衰减强度降低.4.3 流体冲击激励条件下流固耦合声子晶体管路振动特性分析假设输流管路内流体流速初始状态为零, 当开泵或者开阀速度较快时, 流体流量发生急剧变化,(1)t =0.0010 s (2)t =0.0025 s (3)t =0.0050 s (4)t =0.0100 s (a)(b)(1)t =0.0010 s(2)t =0.0025 s(3)t =0.0050 s(4)t =0.0100 s图 10 未充液和充液声子晶体管路不同时刻的速度幅值 (a)未充液管路; (b) 充液管路Fig. 10. Velocity amplitude of liquid-unfilled and liquid-filled phononic crystal pipeline at different moments: (a) liquid-unfilled pipe;(b) liquid-filled pipe.时间/s时间/s(a)(b)频率/Hz频率/Hz速度/m m S s −1速度/m m S s −1|速度(f )|10|速度(f )|图 11 流固耦合声子晶体管路在不同流速下的冲击振动响应 (a)流速为0 m/s; (b) 流速为10 m/sFig. 11. Shock vibration response of fluid-structure interaction phononic crystal pipeline at different velocities of fluid: (a) Flow ve-locity is 0 m/s; (b) Flow velocity is 10 m/s.从而使流体流速发生剧烈变化, 假设流体流速发生变化时满足以下方程:V max t 0t d 式中, = 50 m/s, = 0.0005 s, = 0.0001 s.冲击时间为0.0005 s, 在管路进出水口处拾取响应信号.分别考虑单一材料管路(结构钢管)和布拉格周期管路, 以管路出口处和入口处某一截面的弯曲振动为研究对象, 仿真分析得到管路入口和出口处的速度时间曲线和速度频率曲线.图13是冲击流体激励下结构钢管和声子晶体管弯曲振动响应, 可以发现, 无论是单一材料管路还是声子晶体管路, 在时域上, 其振动响应发生突变出现在流体流速急剧变化时, 随着流速的衰减,管壁振动响应减小, 出口处管壁振动有所延迟, 且声子晶体管路的振动传递到出口处的时间较单一材料管路长, 表明声子晶体管对振动具有一定的抑制作用. 进水口处的振动响应较入水口处的响应略大, 这是由于冲击流体激励到达出口附近需要一定的时间, 流体流动过程中与管壁发生耦合作用, 能量有所耗散, 到达出口处流速略有降低, 因此振动响应略小. 在频域上, 单一材料管路的进水口和出水口处管壁的响应基本吻合, 但在170—210 Hz 范围内存在一定衰减, 而声子晶体管路在415—510 Hz 和575—625 Hz 范围内都存在衰减, 其衰减频段与第3节中带隙相接近但不完全重合, 表明冲击流体激励时, 由于管路流固耦合效应的存在, 流体通过流固耦合面对管壁的作用力分布在整个管路内部,因此在出水口处不仅存在流体在出水口处引起的管路振动, 还存在前端由流体引起的管路振动向管路末端传递. 进一步对比图14中冲击流体激励下结构钢管与声子晶体管在进水口和出水口处振动响应, 可以发现, 进水口处, 单一材料管路在时域和频域上的振动较声子晶体管路略大, 在声子晶体管路进水口处, 由于管路材料环氧树脂材料较软,流体冲击引起的管路振动较单一结构钢管路略大;出水口处, 在时域上, 声子晶体管路出现峰值时间较单一材料管路滞后, 且峰值略小, 表明声子晶体管路可以有效抑制流体冲击引起的管路振动. 在频域上, 声子晶体管路与单一材料管路振动幅值相当, 在415—505 Hz 和585—625 Hz 范围内, 声子晶体管路振动幅值较小, 表明衰减程度明显. 进一步对比进水口处单一材料管路与声子晶体管路的振动响应, 可以发现声子晶体管路对流体冲击引起频率/Hz101010101010|速度(f )|图 12 流固耦合声子晶体管路出口处不同流速冲击振动响应Fig. 12. Shock vibration response of the outlet of fluid-structure interaction phononic crystal pipeline at different velocities of fluid.1010时间/s时间/s(b)频率/Hz 频率/Hz速度/m m S s −1速度/m m S s −1|速度(f )||速度(f )|图 13 冲击流体激励下 (a)结构钢管和(b)声子晶体管弯曲振动响应Fig. 13. Flexural vibration response of (a) structural steel pipe and (b) phononic crystal pipe under shock fluid excitation.。
流固耦合概述及应用研究进展
流固耦合概述及应用研究进展摘要流固耦合力学是流体力学与固体力学交叉而生成的一门力学分支。
顾名思义，它是研究变形固体在流场作用下的各种行为以及固体位形对流场影响这二者交互作用的一门科学。
流固耦合力学的重要特征是两相介质之间的交互作用(fluid．solid interaction)：变形固体在流体载荷作用下会产生变形或运动，而变形或运动又反过来影响流场，从而改变流体载荷的分布和大小。
总体上 ,流固耦合问题按耦合机理可分为两大类:一类的特征是流固耦合作用仅仅发生在流、固两相交界面上 ,在方程上耦合是由两相耦合面的平衡及协调关系引入的;另一类的特征是流、固两相部分或全部重叠在一起 ,耦合效用通过描述问题的微分方程来实现。
1 流固耦合概述1．1引言历史上，人们对流固耦合现象的早期认识源于飞机工程中的气动弹性问题。
Wright兄弟和其它航空先驱者都曾遇到过气动弹性问题。
直到1939年二战前夕，由于飞机工业的迅猛发展，大量出现的飞机气动弹性问题的需要，有一大批科学家和工程师投入这一问题的研究。
从而，气动弹性力学开始发展成为一门独立的力学分支。
如果将与飞机颤振密切相关的气动弹性研究作为流固耦合的第一次高潮的话，则与风激振动及化工容器密切相关的研究可作为流固耦合研究的第二次高潮。
事实上，从美国ASME应用力学部召开的历次流固耦合研讨会上可以看出，流固耦合问题涉及到很多方面。
比如：空中爆炸及响应，噪声相互作用问题，气动弹性，水弹性问题，充液结构内的爆炸分析，管道中的水锤效应，充液容器的晃动及毛细流中血细胞的变形，沉浸结构的瞬态运动，流固相互冲击，板的颤振及流体引起的振动，圆柱由于热交换引起支持附件松动的非线性流固耦合系统，声音与结构的相互作用，涡流与结构的相互作用，机械工程中的机械气动弹性问题等等。
1.2流固耦合力学定义和特点流固耦合力学是流体力学与固体力学交叉而生成的--I'l力学分支。
顾名思义，它是研究变形固体在流场作用下的各种行为以及固体位形对流场影响这二者交互作用的一门科学。
基于声固耦合算法的LNG储罐的内罐湿模态分析
基于声固耦合算法的LNG储罐的内罐湿模态分析
潘传禹;李兆慈;汪常翔;姜铖;李睿麟
【期刊名称】《当代化工》
【年(卷),期】2022(51)11
【摘要】LNG储罐内罐的振动特性是研究储罐抗震性能的基础。
以容积为
1×10~4 m~3双金属全容式储罐为例,利用ANSYS建立流固耦合模型,计算了储罐在空罐、储液静压载荷和流固耦合3种条件下的储罐模态,得到储液晃动基本周期,分析了液体对于储罐模态的影响以及采用流固耦合方法的必要性。
模型的前1 000阶模态计算结果表明,在频率较低时主要振型为局部液面的波动和储液晃动;当频率逐渐升高,体现出环向多波等诸多振型,流固耦合振动以梁式振动为主。
【总页数】6页(P2662-2667)
【作者】潘传禹;李兆慈;汪常翔;姜铖;李睿麟
【作者单位】中国石油大学(北京)/油气管道输送安全国家工程实验室/城市油气输配技术北京市重点实验室
【正文语种】中文
【中图分类】TQ972
【相关文献】
1.基于声-固耦合算法的储液容器湿模态分析
2.大型 LNG 储罐液固耦合模态分析
3.基于声-固耦合算法的在轨补加飞行器贮箱湿模态分析
4.基于声固耦合理论对捕捞航行器的湿模态分析
5.基于声固耦合算法的塔釜持液塔器湿模态分析
因版权原因，仅展示原文概要，查看原文内容请购买。
管道及管路系统流固耦合振动问题的研究动态
第14卷 第3期应用力学学报V o l.14 N o.3 1997年9月CH INESE JOURNAL OF APPL IED M ECHAN I CS Sep.1997管道及管路系统流固耦合振动问题的研究动态α李 琳 喻立凡(北京航空航天大学 北京 100083)摘 要对管道及管路系统流固耦合振动问题在近二十年来的进展作了综述。
根据问题特点,将本课题分为三个分支,即从紊流到振动噪声源的研究,流2弹耦合振动的研究和声2弹耦合振动的研究。
在分别总结这三个分支的研究成果的同时指出了尚需进一步研究的某些问题。
关键词:流固耦合;管路系统;振动噪声1 引 言管路系统流固耦合振动问题有着广阔的工程背景,它的研究成果可直接应用于水利电力、机械、化工、航空航天以及核工程等各个领域。
同时在学术领域,它也是一个十分诱人的课题, Jou rnal of F lu ids and Structu res的创始人M.P.Paidou ssis称之为A M odel D ynam ical P rob2 lem[1]。
因为,它物理模型简单,描述它的数学方程容易简化,特别是管路系统容易实现,这给理论研究与实验研究协同并进提供了极大的方便。
此外,管道虽然是最简单的流固耦合系统,但它却涉及了流固耦合力学中的大多数问题,而且由于它结构的简单性还使得学者们可以分别(或侧重)研究流体的某一特性(如可压性、粘性、流速)对系统的影响。
七十年代以来,管道及管路系统流固耦合的振动问题有了长足进展。
随着问题的深入,关于这方面的研究形成了三个分枝:——研究内流诱发管系振动及噪声幅射机理;——研究具有定常流速的不可压缩流体与管道弯曲振动的耦合以及在此流速下的管道稳定性分析;管道与不可压缩流体的耦合振动可称为液2弹耦合振动。
——研究可压流体中的声波与管道振动的耦合,这种耦合振动可称为声2弹耦合振动。
α来稿日期:1995211230本文将按此三个分支简述课题的进展状况。
声学超材料的声振耦合数值分析
声学超材料的声振耦合数值分析李辉;李旭东【摘要】Based on definite element software ABAQUS,coupling model of acoustics metamaterial with pe-riodicity is built for the purpose of studying the effect of acoustic vibration coupling on structural displace-ment of acoustics metamaterial.When acoustics metamaterial is under acoustic-vibration recombination,fi-nite element model of acoustics metamaterial is simulated in computer with both acoustic and vibration ac-tion by ABAQUS acoustics finite element method and coupled acoustic-structural analysis.The result shows that structural displacement respond of metamaterial under sound-vibration recombination is bigger than that under single sound or vibration.It gives the rules that acoustic vibration coupling affects metama-terial structural displacement,so that ABAQUS can simultaneously simulate value under acoustic vibration recombination.%为了研究声振耦合对声学超材料结构位移响应的影响,以有限元软件ABAQUS为设计平台,建立具有周期性结构的声学超材料耦合模型.当声学超材料结构经受声-振复合环境激励时,通过ABAQUS声学有限元法和声固耦合算法,对声学超材料的有限元模型在声音、振动共同作用下进行计算机仿真.结果表明:声振复合环境激励下的超材料结构位移响应比声、振单独激励下位移叠加后的结果要大.揭示出声振耦合对超材料结构位移响应的影响规律,实现了ABAQUS软件同时对声振复合环境进行数值模拟.【期刊名称】《甘肃科学学报》【年(卷),期】2017(029)006【总页数】6页(P35-40)【关键词】有限元;声学超材料;声振耦合;位移;数值模拟【作者】李辉;李旭东【作者单位】兰州理工大学省部共建有色金属先进加工与再利用国家重点实验室,甘肃兰州 730050;兰州理工大学省部共建有色金属先进加工与再利用国家重点实验室,甘肃兰州 730050【正文语种】中文【中图分类】TB39“超材料”(metamaterials)是21世纪以来出现的一类新材料，通过人工的微结构单元构成复合结构或者复合材料并呈现出天然材料所不具备的超常物理性质。
流固耦合模态分析
流固耦合模态分析流固耦合法广泛地应用于汽车声学和噪声等控制领域，对空腔结构进行流固耦合模态分析，可以了解到声腔对结构模态的影响，为研究耦合系统的声学特性提供可靠的理论参考。
在我们汽车结构振动领域应用相对较多的如油箱流固耦合、声腔耦合分析等，通过进行声腔与钣金耦合分析可以了解车身结构件的振动特性及灵敏度。
模拟流体对结构动力学特性影响的分析方法有很多种，如流固耦合法、虚拟质量法等。
本期通过一个简单的案例分享了解流固耦合的分析流程及相关特性影响。
1相关理论假设流体是均匀、无粘、无旋且可压缩的理论流体，基于小位移理论，并忽略了流固动量传递及局部压力-密度线性关系，其耦合方程为：其中：Ms、Ks ---分别为结构的质量矩阵和刚度矩阵；Mf、Kf ---分别为流体的质量矩阵和刚度矩阵；A ---流固耦合矩阵；Fs、Ff ---分别为结构载荷和声载荷；u ---结构节点位移向量；p ---流体节点压力向量。
2案例实战2.1 本次建立的案例分析模型如图1所示，模拟车身钣金，命名为Structur e，材料属性为钢材，弹性模量E=210000MPa，泊松比NU=0.3，密度RHO=7.85 e-9ton/mm3。
车身结构模型如下图1所示：图1 车身结构模型2.2 根据车身结构模型建立如图2所示的声腔模型。
声控属性为空气，命名为Cavity，设置材料密度1.2e-12ton/mm^3，声速345m/s。
图2 车身声腔模型2.3 建立车身与声腔耦合模型，将车身与声腔模型导入，需要设置以下参数。
图3 结构模态求解设置图4 耦合模态求解设置图5 耦合模态工况设置图6 流固耦合参数设置图7 输出参数设置2.4 模态求解。
为了了解声腔模态与车身结构模态的相互影响，先分别单独计算出声腔和车身结构的自由模态。
声腔的第一阶自由模态为刚体模态，声腔内各点的声压幅值相同，车身结构自由模态前6阶为刚体模态。
以下结果均不包括刚体模态结果列表。
空气滤清器的声振耦合特性分析及优化
Ｋ — — 结构 模型 的刚度矩 阵
— —
结 构模 型 的质量 矩 阵
１ 基础 理 论
１１ 耦 合 系统 方程 ． 在 结构 和 闭空 问 的耦合 问题 中， 闭空 问 中流体
收 稿 日期 ： ０ ００ ．０ 修 改 日期 ：０ ０１．５ ２ １．９１ ； ２ １．２０
作 者 简 介 ： 付 洋 （９ ０ ， ， 南 周 口市 西 华 县 人 ， 士 郭 １ ８ 一） 男 河 硕
图４空气滤清器第 ９ 阶模态
Ｆ ｇ ４ Ｔ ｅ９ ｉ ． ｈ ｍｏ ｌ ｆ ｉ ｆ ｔ ｒ ａ ｄ ｏ ａ ｉ ｅ ｒ ｌ
图 １ 空气 滤 清 器 有 限 元模 型
Ｆｉ １ ｉ ｌｒｆｎｉｅｅ ｅ ｅ ｏ ｌ ｇ． Ａ ｒｆ ｅ ｔ ｌ ｍ ｎｔｍ ｄｅ ｉ ｉ
图 ５ 空 气 滤 清 器 第 １ 阶模 态 ０
下 ， ｒ ａ ｌ 软件 的 间接边 界 元模 型 与有 限元 模 型 Ｖｉｕ ｌａ ｔ ｂ
噪 声 ［。为避 免 空气 滤清 器 壳 体损 坏及 可 能 产 生 的 ２ １
附加 噪 声 ， 要对 空气 滤 清器 进行 分 析 ， 需 并且 根据 需 要进 行优 化 。
ｌ ］ ㈩ ｃＩ ｃ 、Ｄ ８６ ＝Ｊ 【 ０ ｝
Ｔ Ｉ３ Ｒ Ａ 单元 ５ ５ 。图 １ 空气滤清 器有 限元模 型 。 ４个 为
将 空气 滤清 器 的有 限 元模 型 导 入 Ｖｒ ａｌ ｉｕｌａ ｔ ｂ软
件 中， 因发 动 机 的前 ４ 频 率在 ０ ８ ０ Ｉ左 右 ， 阶 － ０ Ｉｚ 故
设定模态求解频率范围为 ０ ８ ０ ｚ ＂ ０ 。计算到３ 阶 ￣ Ｈ ３ 停止 。其 中前 ６ 阶为 刚体模 态 ， 前 ７ Ｏ 取 —１ 阶模 态 为 研 究对象 。
风机叶片流固耦合特性分析与故障诊断
风机叶片流固耦合特性分析与故障诊断一、内容综述风机叶片作为风力发电机中的关键部件，其流固耦合问题一直以来都是研究的热点。
在风能转换过程中，风机叶片不仅要承受气动载荷，还要传递由气流引起的动态力矩，这些力的相互作用导致叶片承受弯曲、扭转、疲劳等多种交变应力，从而引发材料疲劳损坏、结构变形与破坏等问题。
随着风力发电技术的飞速发展，对风机叶片的性能要求也日益提高。
风机叶片在复杂的海洋环境条件下运行，所面临的气动载荷和流固耦合问题更加严重。
在风力发电机的设计、制造和使用维护过程中，对风机叶片的流固耦合特性的深入分析和故障诊断显得尤为重要。
关于风机叶片流固耦合特性的研究成果多集中在叶片的模态分析、振动响应、气弹稳定性等方面。
通过采用先进的计算方法，如有限元法、计算流体动力学（CFD）法等，可以对风机叶片进行精确的结构分析和气动性能预测，为叶片的设计提供理论支持。
现有的研究仍存在一定的局限性。
现有模型多为线性或非线性静态模型，难以模拟风机叶片在实际运行中的流固耦合动态过程。
对于复杂海洋环境条件下的风机叶片流固耦合特性，目前的研究尚不够充分，需要进一步深入探索。
1. 风机叶片在风力发电中的重要性在风力发电系统中，风机叶片扮演着至关重要的角色。
随着全球对可再生能源的不断追求和环境保护意识的逐渐加强，风力发电作为一种清洁、可持续的能源形式，已经成为电力行业的重要组成部分。
风机叶片作为风力发电机的核心部件之一，其性能直接影响到整个风机的运行效率、稳定性和可靠性。
风机叶片的设计直接影响风机的捕获效率。
优秀的叶片设计能够最大限度地捕捉风能，将风能转换为机械能，进而转换为电能。
这不仅提高风机的整体发电效率，还有助于降低能耗，实现节能减排的目标。
风机叶片的长度、材料和结构直接关系到风力发电机的运行寿命。
高性能的叶片能够承受极端的气候条件，如强风、暴雨和雪载等，确保风力发电机在各种环境下都能稳定运行。
叶片材料的选择也至关重要，需要考虑其强度、耐久性和耐腐蚀性等因素，以保证叶片在长期使用过程中的安全和稳定。
ANSYS结构声振耦合解决方案[精]
![ANSYS结构声振耦合解决方案[精]](https://img.taocdn.com/s3/m/ebb841c8941ea76e58fa0474.png)
ANSYS声学模型
声学流体材料 流体密度 流体中声速 边界声吸收系数
F
拉伸
u
压缩
超弹特性
结构材料
弹性材料 超弹材料 粘弹材料 弹塑性材料 其他材料
s
e.
e. 0 粘弹特性 e
ANSYS声学流体单元
Solid185+Neo-Hookean
根据应变大小和 材料数据选择适 当的超弹模型
声学超弹材料
HYPER5x单元超弹性模型
包括 HYPER56, 58, 74 和 158 仅用于模拟几乎不可压缩 M 和 86 模拟Blatz-Ko可压缩泡沫类材料
ANSYS能够定义任意球面波、柱面波和任意方向的平面波等 等
ANSYS提供声传播和声-振耦合分析功能，完整覆盖低高频 声振范围，全面解决多介质、多界面的声传播和结构振动耦 合稳态、瞬态和谐振
动力学
声压级
自由振动—模态分析
声波散射、衍射、传输、
瞬态振动
辐射、衰减等参数
谐振动
结构动态变形应力等
随机振动
声波从空气传入水中
主动声纳探测
声波在管内震荡
主动声纳探测
ANSYS声学模型
FSI—流固界面
结构单元与流体单元接触作用表面 定义流体压力与结构作用界面
Mooney-Rivlin模型
两项形式拉伸应变可达90~100%； 更多项形式可以捕捉工程应力-应变曲线的拐点 5~9项形式应变可达100~200%
Arruda-Boyce模型——8链模型
基于统计的模型，需要的实验数据很少 应变可达300%
机油泵振动噪声解决方案
机油泵振动噪声解决方案
一、机油泵振动噪声
机油泵是汽车发动机上的润滑系统，高速运转会产生较为严重的振动和噪声危害，影响车内外人们的舒适性，同时影响机油泵的使用寿命。
机油泵性能分析中主要考虑结构和流体分析；结构包括强度、密封、疲劳等；流体性能包括排量、扭矩、压力脉动、空化等。
机油泵产生振动和噪声，最主要的是在工作过程中瞬间产生的一些冲击载荷的作用，包括油压和其它一些载荷。
机油泵流致振动噪声分析中，流固耦合、声振耦合分析均较为困难，工程师以实验测试和工程经验为主。
二、机油泵流致噪声机理
1、旋转的离散噪声：
旋转噪声是由于转子或叶片的转动, 周期性地激励泵体而产生的噪声，其频率为齿轮的啮合频率(叶轮机械的叶片通过频率) ，在其基频和谐波处具有较高的噪声值：
n为叶轮/齿轮转速，r/min; z为叶轮叶片数/齿数； i=1,2,3…表示谐波序号。
离散噪声分类:
(1)齿轮泵齿与齿之间的周期性的相互作用引起的噪声；
(2)齿轮泵固定安装不稳产生的周期性的振动噪声。
流固耦合声压分布计算公式
流固耦合声压分布计算公式引言。
流固耦合声压分布计算是指在流体和固体相互作用的情况下，计算声压分布的方法。
在许多工程领域，如飞机、汽车、船舶等，都需要对流固耦合声压进行分析和计算，以评估结构的噪声性能和优化设计。
本文将介绍流固耦合声压分布的计算公式及其应用。
流固耦合声压分布计算公式。
流固耦合声压分布的计算公式涉及到流体和固体的相互作用，需要考虑流场和结构的特性。
一般来说，流固耦合声压分布可以通过以下公式进行计算：P = P0 + ρv^2/2 + ρc^2/2 ζv ξc。
其中，P表示流固耦合声压分布，P0表示静压，ρ表示流体密度，v表示流体速度，c表示结构的振动速度，ζ和ξ分别表示流体和结构的阻尼系数。
这个公式综合考虑了流体和结构的运动对声压的影响，可以较为准确地描述流固耦合声压分布的特性。
在实际工程中，可以通过数值模拟或实验测试来获取流固耦合声压分布的数据，然后利用上述公式进行分析和计算。
应用案例。
流固耦合声压分布计算公式在工程实践中有着广泛的应用。
以下将介绍一些典型的应用案例。
飞机噪声分析。
在飞机设计和运行中，噪声是一个重要的问题。
飞机发动机产生的噪声会对飞机结构和周围环境造成影响，因此需要对飞机噪声进行分析和控制。
流固耦合声压分布计算公式可以用来评估飞机发动机产生的噪声在不同飞行状态下的分布情况，有助于设计和优化飞机结构以减少噪声。
汽车噪声控制。
汽车行驶时产生的噪声也是一个重要问题，特别是在城市交通拥堵的情况下，汽车噪声会对周围的居民和环境造成影响。
流固耦合声压分布计算公式可以用来评估汽车在不同速度和路况下产生的噪声分布情况，有助于设计和优化汽车结构以减少噪声。
船舶结构噪声分析。
船舶在航行时也会产生噪声，对船员和周围海洋生物造成影响。
流固耦合声压分布计算公式可以用来评估船舶结构在不同航行状态下产生的噪声分布情况，有助于设计和优化船舶结构以减少噪声。
结论。
流固耦合声压分布计算公式是对流体和固体相互作用下声压分布进行分析和计算的重要工具。
声波与流体动力学的耦合效应
声波与流体动力学的耦合效应声波与流体动力学的耦合效应是指声波在流体介质中传播时与流体动力学过程相互作用的现象。
声波是一种机械波，它通过介质中的分子振动传播能量，而流体动力学则研究流体的运动规律和性质。
当声波与流体动力学相互影响时，会产生一系列有趣的现象和效应，这对于理解声波在流体中的传播以及流体动力学过程具有重要意义。
一、声波在流体中的传播声波是一种机械波，它需要介质来传播。
在流体中，声波通过分子间的相互作用传播能量，使介质中的分子振动，从而形成声波波动。
声波在流体中的传播速度取决于介质的性质，如密度、压缩模量等。
在流体中，声波的传播速度通常比在固体中要慢，因为流体的分子之间的相互作用力较弱，导致声波传播速度较低。
二、声波与流体动力学的耦合效应1. 声波对流体的影响声波在流体中传播时会对流体产生一定的影响。
当声波通过流体时，会使流体中的分子发生振动，从而产生压力波动。
这种压力波动会引起流体中的密度、温度等参数发生变化，进而影响流体的运动状态。
例如，声波可以在流体中产生压力梯度，从而引起流体的运动，这种现象被称为声波对流体的激励效应。
2. 流体动力学对声波的影响流体动力学过程也会对声波产生影响。
流体的运动状态会改变声波在流体中的传播速度和传播方向。
当流体发生湍流运动时，会使声波传播路径发生扭曲和折射，从而影响声波的传播特性。
此外，流体中的湍流运动还会使声波受到散射和衰减，降低声波的传播效率。
三、声波与流体动力学的耦合效应在实际应用中的意义1. 医学领域在医学领域，声波与流体动力学的耦合效应被广泛应用于超声波诊断和治疗。
超声波可以通过人体组织传播，产生回波信号，从而实现对人体内部结构的成像和检测。
流体动力学过程对超声波的传播和成像具有重要影响，因此需要深入研究声波与流体动力学的耦合效应，以提高超声波诊断的准确性和可靠性。
2. 工程领域在工程领域，声波与流体动力学的耦合效应被应用于声纳技术和声波控制。
ANSYS 声学计算算例 流固耦合
ANSYS 声学计算算例水下圆柱壳体的建模与声学分析使用有限元软件ANSYS进行计算和分析时水下环肋圆柱壳体有限元模型的建立及结构声学分析主要分为以下一些步骤：1.建立壳体的实体模型（包括有圆柱壳体的建立，给圆柱壳体加环肋）；2.圆柱壳体外部流体介质的生成；3.对圆柱壳体和流体介质进行有限元4.设置流固耦合单元，并设置外部声场边界条件；5.在求解器中进行振动模态求解和受激励的谐响应求6.求解结果进行后处理分析。
，1.建立壳体的实体模型这个步骤主要是在预处理模块（PREP7）中完成首先是根据要建立的实体模型，进行单元的选取和定义这些单元的物理属性，水下圆柱壳体半径与壳体壁厚的比超过了20，根据ANSYS中单元的使用原则可以选用Shell63号薄壳单元，这种单元的有限元计算原理在前面已经介绍；环肋选用梁单元，ANSYS 提供了多种梁单元的结构形式，其中Beam188号梁单元符合作为壳体加强筋及肋骨的使用，所以在水下圆柱壳体环肋选用T的Beam188号梁单元进行建模；而流体介质根据分析中用途的不同要定义两种，一种是流体介质中的单元Fluid30号流体介质单元，一种是流体与结构接触的流固耦合面的单元选用Fluid30号流固耦合单元，在实际建模操作中还需要定义一种用于平面声场的29号单元（在计算中未用到，但在建模中需使用）；共需要定义五种单元。
Shell63壳体单元与Beam188梁单元为同一种材料，所以物理属性相同。
而Fluid30流体介质单元与Fluid30流固耦合单元物理属性也相同，及在分析中只需要定义两个物理属性即可。
具体的使用APDL命令定义为：/prep7 ！进行预处理模块et，1，30，！定义1号单元为Fluid30 流固耦合单元et，2，29 ！定义2号单元为Fluid29平面流体单元et，3，30 ，，1 ！定义3号单元为Fluid30流体介质单元et，4，63 ！定义4号单元为Shell63壳体单元et，5，188 ！定义5号单元为Beam188梁单元r，4，0.002 ！定义4号单元的厚度为2㎝mp，dens，4，7800 ！定义4号物理属性包括有密度mp，ex，4，2.1e11 !杨氏模量、mp，nuxy，4，0.3 ！泊松比mp，sonc，1，1460 ！设置水中声速mp，dens，1，1000 ！设置流体密度sectype，1，beam，T，! 选取T型梁secoffset,，orig ！设置梁的方向secdata，0.04，0.05，0.002，0.02，0，0，0，0，0，0 、所建立的圆柱壳体的参数：圆柱长为50 ㎝，半径为25 ㎝，壳体的壁厚为2 ㎝，cyl4，0，0，0.25，，5 ！形成圆面k，9，0，0，0 ! 定义原点k，10，0，0，0.5lstr，9，10 ！通过原点作直线adrag，5，6，7，8，，，9 ！通过放样形成圆柱wpoff，0，0，0.1asel，s，，，2，5asbw，all，，，！移动工作平面与选取的侧面相切…… ！重复上面操作，形成四个环肋面wpoff，0，0，-0.4 ！工作平面回到原点位置上k，31，0.2，0，0.1 ！定义环肋的方向点lsel，s，，，20 ！选择要划分为环肋的线段latt，4，5，5，，31，40，1 ！定义线段物理属性lesize，20，，，6 ！划分数目secnum，1lmesh，20 ！划分线段将上述的操作完成以后，壳体的模型基本完成,具体结构如图示图2-6 环肋圆柱壳体模型图2.圆柱壳体外部流体圆柱壳体外部的流体介质主要通过设置好的平面流体单元沿指定的线段进行放样，形成立体的流体介质单元。
固流耦合模态
固流耦合模态
固流耦合模态是指结合了固体力学和流体力学的模态分析方法。
它可
以模拟流固耦合系统中的振动和波动现象，可以广泛应用于制造、力学、机械等领域。
固流耦合模态分析方法的基本原理是将流体、结构和流体-结构相互作用的数学模型组合起来进行计算。
这种方法可以识别流体对结构的振
动响应产生的影响。
具有以下几个特点:
1.高可靠性。
固流耦合模态分析考虑了结构和流体之间的相互作用，因此更加准确及可靠。
2.高灵敏度和精度。
固流耦合模态分析考虑了流体质量和结构刚性等变量的影响，还考虑了流固相互作用，具有高灵敏度和精度。
3.多领域交叉应用。
固流耦合模态分析方法可应用于制造、机械、航空、水文、机电等多个领域，是一种交叉应用广泛的分析方法。
特别是在
气动声、地震、风力以及流体介质中的振动分析方面得到了广泛应用。
4.使用方便。
与其他分析方法相比，固流耦合模态分析方法的使用越来越方便，其软件和算法也越加完善，有助于提升工作效率和准确性。
固流耦合模态分析的应用范围十分广泛。
在制造业领域，可以用来分
析枪口、叶轮、飞机翼等各种结构件的振动和波动问题。
在机械领域，可以用来分析钢铁设备、风机、汽车等的振动噪声问题。
在航空领域，可以用来分析飞机机翼的振动、噪声等问题。
在水文领域，可以用来
分析谷物筒仓等农业设备的振动噪声问题。
总的来说，固流耦合模态分析方法凭借其精确度高、应用广泛的优点，在工程领域中得到了广泛应用。
因此，未来固流耦合模态分析方法将
在制造、力学、机械等领域的精细化发展方面更为得到发展。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
流固耦合力学是流体力学与固体力学交叉而生成的一门力学分支，它是研究变形固体在流场作用下的各种行为以及固体位形对流场影响这二者相互作用的一门科学。
流固耦合力学的重要特征是两相介质之间的相互作用，变形固体在流体载荷作用下会产生变形或运动。
变形或运动又反过来影响流，从而改变流体载荷的分布和大小，正是这种相互作用将在不同条件下产生形形色色的流固耦合现象。
流固耦合问题可由其耦合方程定义，这组方程的定义域同时有流体域与固体域。
而未知变量含有描述流体现象的变量和含有描述固体现象的变量，一般而言具有以下两点特征：1）流体域与固体域均不可单独地求解
2）无法显式地削去描述流体运动的独立变量及描述固体现象的独立变量
从总体上来看，流固耦合问题按其耦合机理可分为两大类:
第一类问题的特征是耦合作用仅仅发生在两相交界面上，在方程上的耦合是由两相耦合面上的平衡及协调来引入的如气动弹性、水动弹性等。
第二类问题的特征是两域部分或全部重叠在一起，难以明显地分开，使描述物理现象的方程，特别是本构方程需要针对具体的物理现象来建立，其耦合效应通过描述问题的微分方程来体现。
实际上流固耦合问题是场（流场与固体变形场）间的相互作用：场间不相互重叠与渗透其耦合作用通过界面力（包括多相流的相间作用力等...）起作用，若场间相互重叠与渗透其耦合作用通过建立不同与单相介质的本构方程等微分方程来实现。
（1）试列举出至少三个经典的加权残值方法，并简述伽辽金法的基本思想。
最小二乘法、配点法、子域法、伽辽金法、矩量法；
伽辽金法采用微分方程对应的弱形式，其原理为通过选取有限多项式函数（又称基函数或形函数），将它们叠加，再要求结果在求解域内及边界上的加权积分（权函数为势函数本身）满足原方程，便可以得到一组易于求解的线性代数方程，且自然边界条件能够自动满足。
应用这种方法可以将求解微分方程问题（通过方程所对应泛函的变分原理）简化成为线性方程组的求解问题。
而一个高维（多变量）的线性方程组又可以通过线性代数方法简化，从而达到求解微分方程的目的。
（2）与有限元方法（finite element method）相比无网格方法（Mesh-less method）有何优势？
（3）有限元方法（finite element method）与瑞利-里兹法有何联系？并简述瑞利-里兹法的基本思想。
瑞利-里兹法的基本思想是把连续系统离散化为有限自由度系统，然后根据机械守恒定律进行计算。
（4）请简述用有限元方法分析梁的横向振动问题的主要流程？
（5）有限元（FEM）和无限元（IEM）的主要适用范围有什么不同？写出无限元（IEM）的主要特点。
有限元法从它最初应用的固体力学领域，已推广应用到温度场、流体场、电磁场、声场等其他连续介质领域。
在固体力学领域，有限元法不仅可用于线性静力分析，也可用于动态分析，还可用于非线性、热应
力、接触、姗变、断裂、加工模拟、碰搜模拟等特殊问题的研究
无限单元是为了克服有限单元在求解无界域问
题时的不足而提出的因而原则上说可以用有限元方法求解但涉及无界域的物理问题都特别适合于用无限元方法求解无
9。