TMS BOOK · ACADEMY 讲义

流固耦合(FSI)与流致振动

大纲的完整展开版——讲师授课蓝本 / 学员自学材料

J5-04 流固耦合(FSI)与流致振动

课程代码 J5-04 · 板块 J 仿真分析 CAE / 结构/热应力/疲劳 FEA 时长 约 4.5 小时(8 讲) 前置 J5-03 振动、模态与随机振动疲劳;J2-01 CFD 基础、湍流模型与网格划分 本讲义定位 讲师授课蓝本 / 学员自学讲义,是 J5-04 大纲的完整展开版


按任务导航

你手上的事 直接看 还要配
审一份别人交来的 FSI 报告,判它能不能信 第1讲:责任链六段,逐段问「这一段错了,谁会告诉我」 第7讲的交接卡四件;第6讲看时间步是不是由分析频段倒推
给一根具体的管列出激励源清单 第2讲:七域按它在回路上的位置挂点,⛔ 不按机理回忆 第3讲三条量化路(阶次/水锤/宽带);第2讲涡脱的两步射程判据
定「这次要不要上双向耦合」 第5讲:五条路,⛔ 不是「单向/双向」两条 第5讲三条判据:质量比 ρ_l/ρ_s · 变形改不改流道 · 是不是稳定性问题
已经出了异响或焊缝裂纹,要交一条改进 第8讲:抑制四类按「改哪一端」分 第7讲四对频率对账;第8讲案例的四步砍源
写一份交给疲劳的输入规格 第7讲:应力口径 · 频段与带宽 · 是输入还是响应 · 工况时长与占比 第1讲两条口径:p̃=0-peak、f_n=湿模态

⛔ 本课不覆盖:

栏一 · 不讲哪些题目——随机振动模型本体、PSD 三区间判读与 Miner 累积 → J5-03 振动、模态与随机振动疲劳;疲劳评价方法本体与自建 S-N → J5-07 钎缝/焊缝接头的疲劳评价方法:三级应力口径、等效结构应力法与自建 S-N,许用值与安全系数 → J5-01 换热器/冷板结构强度与耐压 FEA;CFD 的湍流模型、y+ 与网格无关性 → J2-01 CFD 基础、湍流模型与网格划分(Fluent/STAR-CCM+);汽蚀会不会发生、NPSH 裕量够不够 → G5-04 泵的气蚀、空气分离与可靠性;泵内部脉动的产生机理 → G5-07 电子油泵本体零件级设计:容积式构型几何与排量、间隙链、困油与流量脉动;2 kHz 以上的噪声与声辐射 → J5-05 声学 CAE:气动噪声与结构声辐射仿真;隔振件本身的设计与选型 → I4-01 支架、悬置与减振结构设计;锤击模态与振动台的执行 → L6-04 结构动态特性试验:模态、动刚度与传递率实测;通用结构建模(单元选型、连接与约束)→ J5-06 结构 FEA 建模方法基础:单元、约束、本构与结果可信度;阀切换速率的标定 → G8-02 多通阀 + 泵 + 传感器的一体化集成设计。

栏二 · ⛔ 真正不给的数——① ISO 16750-3 的限值、试验等级、谱形与时长:本课未取得原文核对,只写去向,版次以现行目录为准,引用前须自行核原文;② 管路振动与异响的通过/不通过阈值:无统一编号的通用标准,去本企业规范与整车 NVH 目标值;③ 具体泵/压缩机的 N_pocket 与脉动幅值:向供应商索取排量结构参数与脉动测试报告,或按 G5-07 电子油泵本体零件级设计:容积式构型几何与排量、间隙链、困油与流量脉动 从泵结构反推。

引言:六段里只有一段会自己报错

项目上丢过来的从来不是一道仿真题,而是一个结果——「这根 Chiller 进口管有异响」「这道焊缝三千公里就裂了」。于是有人打开 CFD 与结构求解器,把两边接起来跑一遍,云图很漂亮,模态列表打出来,结论写「激励频率与固有频率相差较远,不共振,安全」。这份报告的每一步都跑通了,而它可以整份是错的,并且不会有任何东西提示你。本课要教的,就是在这种报告面前问对问题。

钉子①:本课的产出不是「一个振动幅值」,是两样东西。一样是一次频率带对账——激励带、管系声学模态带、结构湿模态带,三条带画到同一根频率轴上,看它们重不重叠;另一样是一个口径写死的应力幅 σ_a(宽带激励时是应力 PSD),连同一张交接卡交给 J5-03 振动、模态与随机振动疲劳 去算寿命。到此为止:本课⛔ 不产出寿命、⛔ 不产出许用值、⛔ 不判汽蚀会不会发生。读者最常越的界,是拿 FSI 算出的振动幅值(mm/s 或 g)直接当疲劳输入——中间缺的是应力提取与 J5-03 振动、模态与随机振动疲劳 的方法,那一步不在本课自动发生。

钉子②:这条链有六段,而只有第 ⑤ 段会自己报错。① 激励识别与量化 · ② 流体域传递与声学放大 · ③ 耦合方式选型 · ④ 结构建模与湿模态 · ⑤ 耦合求解与映射 · ⑥ 应力提取与交接。第 ⑤ 段错了,求解器会喊出来:映射失败、不收敛、时间步不同步、界面能量不守恒。其余五段错了,一律长成同一副样子——算出来了、云图挺好看、结论是「不共振,安全」。所以判「这次 FSI 能不能信」,逐段问的是「这一段错了,谁会告诉我」,⛔ 不是问「求解跑通了没有」。

六段责任链横向排列,每段方框下方回答同一个问句「这一段错了,谁会告诉我?」:只有第 ⑤ 段耦合求解与映射写着求解器会报映射失败/不收敛/时间步不同步/界面能量不守恒,其余五段一律写「⛔ 没有人」;链左端入口框写「脉动/汽蚀已经算出来」,右端出口框写「σ_a 或应力 PSD + 四件口径」。本图为定性示意,⛔ 不得据图取任何数值。
图1 读出三件:① 六段里只有一段会自己报错 ⇒ 可信度⛔ 不能靠「求解跑通了」来证明;② 本课的入口与出口都被上游定死,⛔ 不在本课中间产生;③ 评审一次流固耦合分析,该问的是「另外五段各由谁把关」。本图为定性示意,⛔ 不得据图取任何数值。

反钉子:最容易被读反的一条是「频率差得远 ⇒ 安全」。在车用两相/液冷管路上,用来判「差得远」的那两个数被系统性地朝同一个方向算高了——于是错误不是随机的,它稳定产出假绿灯。

读者会做的那个具体动作是这样的:拿求解器直接给出的干模态 f_n,dry(比如 200 Hz),拿用液相声速 c_l 算出来的管系驻波 f_a(比如 c_l/2L_ac ≈ 500 Hz),和压缩机排气阶次 f_p(比如 100 Hz)摆在一起,得出「三个数互不接近、留有余量」⇒ 不做后续 FSI,卡箍跨距按布置方不方便定,缓冲容积省掉。(这三个数是一组教学假设值,⛔ 不对应任何平台,只用来演示这个判断长什么样。)

两个误差为什么同向:

结构侧——管内充满液体时,流体附加质量把固有频率压下来,裸管干模态是上界。薄壁液冷金属管充液后一阶弯曲频率相对干模态典型下降 10%~30%,多数在 20% 附近;卡箍若按理想固支建而不按实测刚度建,高估幅度同样可达数十个百分点。⚠ 典型值、⛔ 不是验收依据,你的平台须实测确认。→ 拿到它的路:用声固耦合或附加质量做一次干湿对比计算(第 5 讲第 ④ 条路),再用 L6-04 结构动态特性试验:模态、动刚度与传递率实测 的锤击实测标定卡箍刚度。

流体侧——EXV 之后是两相段,而均相混合物的声速比液相和气相都低(Wood 式给出一条 U 型曲线,谷底在中间空泡份额处)。可复算的一组:α=0.5、ρ_l/ρ_g≈90、c_l≈500 m/s、c_g≈145 m/s 代入 1/(ρ_m c_m²) = α/(ρ_g c_g²) + (1−α)/(ρ_l c_l²),得 c_m ≈ 30 m/s;把 α 扫到 0.9 与 0.95,c_m 也只有约 50 与 65 m/s ⇒ 常见空泡份额区间内,两相段的 c 一律落在几十 m/s 量级,比液相低一到两个数量级。⚠ 输入的四个物性须以本平台工况现算,⛔ 不许沿用别的工况;该式给的是量级判据,⛔ 不作设计输入。后果是:真实的 f_a 及其各次谐波整片落进激励带,而按 c_l 算会把它们整体推高一个量级——那一整对频率比对连算都不会算到。

两个偏差都朝「更安全」的方向走 ⇒ 余量是假的。⛔ 反过来说也要点破:这不是「双向 FSI 才能救」。用双向强耦合而仍然拿液相声速与干模态起步,上面两处错的地方一处都没动,只是把同一个错误算得更贵。

本课的处置写在口径上,第 1 讲一次钉死、全课不改:裸写 f_n 恒指湿模态(含流体附加质量与内压刚化,卡箍按有限刚度),干模态一律写 f_n,dry 并注明是上界;裸写 c 恒指管内有效声速(含管壁柔度的 Korteweg 修正),两相段另按 Wood 式分段算,⛔ 全程不取一个 c。

本课在链上的位置:入口接三样——G5-07 电子油泵本体零件级设计:容积式构型几何与排量、间隙链、困油与流量脉动 给泵内部脉动的量级与阶次、G5-05 电子水泵水力部件设计:叶轮/蜗壳几何、内部间隙与轴向力 给几何改动与脉动之间的因果方向、G5-04 泵的气蚀、空气分离与可靠性 给汽蚀会不会发生;出口交两样——J5-03 振动、模态与随机振动疲劳 拿走应力做寿命、J5-01 换热器/冷板结构强度与耐压 FEA 拿走结果对许用值。分析频段以 0~1 kHz(结构疲劳带)为主,异响判读延伸到 2 kHz,再往上归 J5-05 声学 CAE:气动噪声与结构声辐射仿真。

八讲怎么走:第1讲把符号与口径一次钉死,并立起六段责任链与本课边界;第2讲铺开激励源的七个域,顺带把涡激振动的射程切清楚(光滑管内流没有周期性涡脱);第3讲把激励量化成数——阶次频率、Joukowsky 压升、弯头之间的不平衡力,以及宽带与阶次两条不同的路;第4讲讲流体域怎么传递与放大——有效声速、两相声速、分段做声学、亥姆霍兹腔;第5讲讲耦合方式的五条路与湿模态的三个来源;第6讲讲耦合求解本身:界面合力守恒、时间步由频段倒推、压力提取的三个陷阱;第7讲把四对频率画到一根轴上做判读,并写出交接卡;第8讲讲抑制手段的四类、案例决策链、误区里与抑制设计直接相关的六条(全集七条见「常见误区」一节)与一次可手算的练习。

读完这一节,你今天就能用的是一个「三带对账起手式」——拿到一根具体的管,不打开任何求解器,先写三行:激励带(f_p = n·(rpm/60)·N_pocket 在工况转速范围内扫过的那一段,n 至少取到 3)· 声学带(按管段分别取 c,硬管用有效声速、两相段按 Wood 估,f_a 按 c/(2L_ac) 或 c/(4L_ac),边界条件判错一次差一倍)· 结构带(f_n 按湿模态口径估,⛔ 不用求解器直接给的干模态)。三行写完,用 ±20% 作粗筛看有没有重叠——⚠ 它是粗筛线、⛔ 不是验收线,真实带宽由 Q_amp ≈ 1/(2ζ) 决定。三行里只要有一行是用干模态或液相声速写的,这次对账的结论就不能用;下一个动作是先把那一行的口径换掉再重画,⛔ 不是去加卡箍。

第 1 讲 先把词定死:符号口径、六段责任链,与本课的入口和出口

这一讲不讲方法,讲的是后面几讲能不能被判对错的前提。

流致振动横跨四套学科——管内流动、管系声学、结构动力学、疲劳。四套各有各的符号习惯,而它们在同一根管上撞得极密:D 在流体力学里是绕流体的特征尺寸、在管路工程里是管内径;c 既可以是流体的自由声速、也可以是含管壁柔度的管内有效声速;f_n 既可以是裸管的干模态、也可以是充液后的湿模态;Q 在三套学科里分别是流量、热量和共振放大因子。这些不是「写法不统一」的排版问题——两义之间差的不是精度,是量级。

所以本讲做五件事:把符号与单位一次钉死(1.1~1.3);把本课最容易出事的那个数值口径——压力脉动幅值 p̃——单独拎出来讲透(1.4);把一次流致振动分析拆成六段并逐段问「这一段错了,谁会告诉我」(1.5);说清为什么「求解跑通了」在 FSI 里是极弱的证据(1.6);把本课的入口、出口与射程边界写成可执行的交接口径(1.7)。

读完这一讲,你手上应该多三样今天就能用的东西:一张符号口径卡、一张六段自检卡、一张交接卡。1.8 把这三张卡整理成可以直接贴进评审模板的形式。

1.1 四套学科的符号在这里撞车,本节一次钉死,全课不改

本课最贵的两个错误都是符号错,⛔ 不是方法错。下表先钉死第 1 讲就要用到的那些;随后各讲会陆续引入声学与两相流侧的符号,含义与单位同样以本表的口径为准,全课不改。

符号 本课的含义 单位 在本课⛔ 不表示
D_i 管内径 mm(结构侧)/ m(代入公式时) ⛔ 不是绕流特征尺寸,⛔ 不得代入 Strouhal 公式或涡脱频率式
d_o 绕流体或内部障碍的特征尺寸——外部绕流取扁管厚度或圆柱直径,内部涡脱取孔板孔径、温度传感器套管外径或盲支管内径 m ⛔ 不是管内径
A_i 管内流通截面积(=π·D_i²/4),用于压力波不平衡力与动量项 mm²(结构侧)/ m²(公式代入) ⛔ 不是消音腔颈口面积(那是 A_n)
A_s 管壁的结构截面积(环形截面,=π·[(D_i/2+e)²−(D_i/2)²]),只用于梁模态式里的单位长度质量 ρ_s·A_s mm²(结构侧)/ m²(公式代入) ⛔ 不是管内流通截面积(那是 A_i)、⛔ 不是颈口面积(那是 A_n)
v 管内截面平均流速 m/s ⛔ 不是障碍处的缩流速度
c 管内有效声速(含管壁柔度的 Korteweg 修正) m/s ⛔ 不是流体的自由声速(那是 c_f)
p̃ 压力脉动幅值,恒指 0-peak(零到峰) kPa ⛔ 不是峰峰值、⛔ 不是 RMS
L_span 两卡箍之间的自由跨距(结构模态用) mm —
L_ac 声学管段长度(驻波与波头穿越时间用) m —
f_n 结构固有频率,恒指湿模态(含流体附加质量与工作内压刚化) Hz ⛔ 不是干模态
f_n,dry 干模态固有频率(不含流体) Hz ⛔ 不得当成 f_n 用
σ_a 应力幅——本课交给 J5-03 振动、模态与随机振动疲劳 的那一个量 MPa ⛔ 不是许用应力、⛔ 不是寿命

配套四条命名纪律,同样在这里写死,全课不改:

  1. ⛔ 本课禁止裸写 D。管内径一律写 D_i,绕流/障碍特征尺寸一律写 d_o。理由在下一节展开。
  2. ⛔ 本课禁止裸写 L。它在本课有两义且量级不同:结构侧的卡箍跨距(毫米量级)与声学侧的管段长度(米量级)。凡写 L 必带修饰语——L_span 或 L_ac。
  3. ⛔ 本课禁止裸写 Q。共振放大因子一律写 Q_amp(定义与取值见第 2 讲),体积流量一律写 q_v。Q 这个字母在流量、热量、放大因子三处都有主,本课把它整个空出来。
  4. n 只表示谐次(阶次的整数倍,n=1,2,3…)。⛔ 它不表示转速,也⛔ 不表示 S-N 曲线指数——后者是 J5-03 振动、模态与随机振动疲劳 的 m,本课不引入。

⚠ 关于 f_n 再多说一句,因为它是本课反钉子的一半:裸写 f_n 一律是湿模态,干模态必须写全 f_n,dry 并注明它是上界。薄壁液冷金属管充液后一阶弯曲频率相对干模态典型下降 10%~30%,多数在 20% 附近;卡箍若按理想固支建而不按实测刚度建,f_n 还会被再高估数十个百分点。⚠ 典型值、⛔ 不是验收依据,它随管径壁厚比、流体密度与卡箍形式变化很大。→ 拿到它的路:用声固耦合或附加质量做一次干湿对比计算(第 5 讲给做法),再用 L6-04 结构动态特性试验:模态、动刚度与传递率实测 的锤击实测标定卡箍刚度与阻尼。⛔ 全课任何一处频率对账,两侧都不许出现 f_n,dry。

四个最常见的翻车点,一并点在这里:① 从外部资料抄公式时把它的 D 直接当成本课的 D_i;② p̃ 从供应商数据表抄来是峰峰值,而本课按 0-peak 用,结果差 2 倍(1.4 单讲);③ 把空泡份额 α 与干度 x 混用(两者靠密度比换算,混用会把两相声速算错一个量级,第 4 讲会用到);④ 把 Q 当放大因子写在流量旁边,同一页上两个 Q 指两件事。

1.2 同一根管上有两个长度尺寸,只有其中一个能进 Strouhal 公式

上一节那条「禁止裸写 D」看起来像洁癖,其实它挡的是本课最容易发生、也最难被发现的一次误用。

Strouhal 数的定义式 St = f·d_o/U 里,长度是被绕流的那个物体的尺寸——圆柱直径、扁管厚度、孔板孔径。而在管路工程里,D 几乎总是指管内径。同一个字母,在同一门课里两义。于是一个非常自然的动作出现了:拿管内径代进去,算出一个几十赫兹的「涡脱频率」,再拿它去和结构模态比。

这个数是假的,而且它比没有数更危险——一个算得出来的数,最像已经想清楚了。

管纵剖示意分左右两半:左半光滑直管段标注管内径 D_i 与管内平均流速 v,旁有红色禁止标记「D_i ⛔ 不代入 St、⛔ 不代入 f_v」;右半同一根管上串了孔板、插入式温度传感器套管、盲支管三种内部障碍,每种各标出自己的特征尺寸 d_o(孔径/套管外径/支管内径),障碍下游画出周期性涡列,旁标 U(障碍处缩流速度,⛔ 不等于 v)与 f_v = St·U/d_o。图底横带写「先判场景,再算频率」。本图为定性示意,⛔ 不得据图取值,尺寸比例不代表真实几何比例。
图2 同一根管上的两个长度尺寸各归各的公式。要读出三件:① 只有障碍尺寸 d_o 能进 Strouhal 公式,管内径 D_i 进不去;② 公式里的速度是障碍处的缩流速度 U,⛔ 不是管内平均流速 v;③ 光滑管段那一半图上没有涡列——它不产生周期性涡脱,因此那一段根本不进入频率计算。本图为定性示意,⛔ 不得据图取值,尺寸比例不代表真实几何比例。

本课因此把「这段管会不会涡脱」定成两步,顺序不能颠倒:

第一步,先判场景。只有两类流道允许往下算——外部绕流(换热器风侧扁管尾迹、百叶窗后缘),以及内部障碍的下游(孔板、插入式温度套管、盲支管)。光滑管内的充分发展流没有周期性涡脱,它的壁面脉动来自湍流的宽带随机,⛔ 不是窄带涡脱。

第二步,才算频率。而且特征尺寸只准取 d_o,速度只准取障碍处的缩流速度 U——障碍把流道缩窄了,那里的速度高于管内平均流速 v。

场景判定的完整分支、St 的取值区间与去向,在第 2 讲展开。这里只要求你记住一条可以今天就用的判据:看见一份 FSI 报告里出现「管内径代入 Strouhal」,不必看后面,那一段的结论作废。

⚠ 走到「不涡脱」这个结论不是结束。车用内流的真问题不在涡脱那一路,而在「脉动阶次/声学模态/结构模态」三条频率带落不落在一起——那是第 3、4 讲的事。在光滑管上算出「不涡脱」就结案,是本课要防的另一种假绿灯。

1.3 单位制、压力单位与本课的分析频段

单位制。本课两侧各用各的,接口只有一个换算,写死在这里:

压力单位。脉动一律以 kPa 报(车用量级),进结构侧之前显式除以 1000。⛔ 全课不出现「bar」——它在不同资料里既可能表示绝压也可能表示表压,而两者之间差的正好是一个大气压,在低压脉动上这是个致命的量级。见到「bar」,先回问是绝压还是表压,问不到就不用这个数。

分析频段。本课的射程是 0~1 kHz(结构疲劳带)为主,异响判读延伸到 2 kHz。⛔ 2 kHz 以上的汽蚀高频噪声、声辐射与主观评价一律不在本课,去向 J5-05 声学 CAE:气动噪声与结构声辐射仿真。

这个频段不是随便划的,它由两头卡住:

⚠ 频段一旦定死,采样率与耦合时间步就跟着定死(第 6 讲会把这条链算完),⛔ 不许中途改。中途把频段改窄是本课要防的四种「让计算变顺同时让结论变安全」的错之一,见 1.6。

1.4 压力脉动幅值 p̃:本课恒指 0-peak,口径错一次寿命差 32 倍

p̃ 是本课从上游拿到的唯一一个数值输入。它错了,后面全错——频率带对账还能靠频率本身兜住,应力幅 σ_a 则是直接按比例跟着错。而它恰恰是本课所有量里最容易口径不清的一个。

三种写法在资料里混用:

口径 与本课口径的关系 谁常用它
0-peak(零到峰) 本课口径 本课全程
peak-to-peak(峰峰值) = 2 × 0-peak(恒等定义) 泵与压缩机数据表最常给、台架报告
RMS(均方根) 正弦波下 = 0-peak/√2 ≈ 0.7071 × 0-peak(√2=1.4142,纯数学常数) 振动与声学侧的习惯

⚠ RMS 那一行只对正弦成立。对塞状流的脉冲串与汽蚀的宽带脉动,它们根本不是正弦,⛔ 不能用这个系数换算——那两类必须直接给功率谱密度(第 2 讲讲谱形,第 6 讲讲怎么进模型)。

错一次的代价,可以逐位复算:把峰峰值当成 0-peak 用 ⇒ p̃ 偏大 2 倍 ⇒ 结构响应线性 ⇒ σ_a 偏大 2 倍 ⇒ 在 J5-03 振动、模态与随机振动疲劳 的 S-N 指数 m ≈ 5 下,寿命差 2⁵ = 32 倍(纯算术)。方向反过来(把 0-peak 当峰峰值折半)同样是 32 倍,只不过这一次错在把寿命算长了——那是没人会来查的那个方向。

⇒ 本课的硬规定:引用任何外部脉动数据之前,必须先问清它的口径并写在数字旁边;⛔ 全课不得两种口径并存。问不到口径的数,按「未取得」处理,⛔ 不要挑一个自己觉得合理的解释填进去。

三个易错点:

  1. 「±50 kPa」这个写法本身有歧义——它可能是 0-peak 50 kPa(峰峰 100),也可能是有人把峰峰 50 写成了这个样子。见到带正负号的脉动值,必须回问。
  2. 台架测得的是传感器那一点的值,而结构关心的是弯头、变径、盲端那一点的值。两点之间隔着管系声学放大——同一个源,在驻波的波腹和波节上能差出很多倍(第 4 讲)。⛔ 不能把测点的幅值直接搬到关注点上用。
  3. ⛔ 不要把不同工况的幅值和频率拼成一个「最恶劣工况」。那是拼出来的、不存在的工况;它给出的余量判断既不保守也不激进,只是没有物理意义。要覆盖多工况,就逐工况各跑一遍、各自带自己的时长占比交出去(1.7 的交接四件里有这一项)。

1.5 六段责任链:只有第 ⑤ 段会自己报错

一次流致振动分析可以拆成六段。真正要问的不是「这六段是什么」,而是每一段下面那句话:

段 这一段交出什么 这一段错了,谁会告诉我?
① 激励识别与量化 哪几个激励源、各自的频率与幅值 p̃ ⛔ 没有人
② 流体域传递与声学放大 分段的有效声速 c、管系声学模态 f_a、腔体 f_H ⛔ 没有人
③ 耦合方式选型 走哪一条路、为什么这条够 ⛔ 没有人
④ 结构建模与湿模态 f_n(湿)、振型、卡箍刚度与阻尼 ζ ⛔ 没有人
⑤ 耦合求解与网格/时间步映射 收敛的耦合解、界面载荷 求解器会:映射失败/不收敛/时间步不同步/界面能量不守恒
⑥ 应力提取与交接 σ_a 或应力 PSD + 口径四件 ⛔ 没有人

那五段的失效形态是同一个:算出来了、云图挺好看、结论是「不共振,安全」。报告照常写得出,评审会上照常有人点头。

这张表是本课的骨架,理由很直接:读者对 FSI 的既有印象几乎全落在第 ⑤ 段(「怎么把 CFD 和结构求解器接起来」)——而这门课要治的七条典型错误(全集见「常见误区」一节)里,没有一条落在第 ⑤ 段,它们分布在 ①②③④⑥ 上。你把第 ⑤ 段做到极致,也不会碰到它们中的任何一个。

⇒ 把上表右列拆成六个问句,就是一张今天就能用在你手上那份报告上的自检卡:

  1. ① 这份分析里列了几个激励源?除了泵和压缩机的阶次,阀件瞬态、两相冲击、空化、内部障碍涡脱、并联流道的自激振荡,各自判过没有?(全集在第 2 讲)
  2. ② 声速取的是哪一个?是流体自由声速,还是含管壁柔度的管内有效声速?软管段和两相段分段取了没有?
  3. ③ 为什么选这条耦合路?它排除掉的另外几条,是因为不需要,还是因为没想到?
  4. ④ 报告里的固有频率是湿模态吗?卡箍是按实测刚度建的,还是按理想固支?
  5. ⑤ 界面两侧的合力核对过没有?(这一段是唯一会自己报错的,⛔ 别浪费它)
  6. ⑥ 交出去的是应力幅还是应力 PSD?它的应力口径、对应频段、是输入还是响应、对应多长时间——四件写全了吗?

三个易错点:① 把六段当成软件操作顺序背下来,而它的用法是六个问句;② 认为「耦合收敛了」说明前四段对——耦合收敛只说明数值稳定,⛔ 与载荷对不对无关;③ 把「自检」理解成再跑一遍同一个模型——再跑一遍只会得到同一个数。

1.6 「求解跑通」是极弱的证据:四种错会让计算更顺,同时让结论更安全

上一节说五段不会报错。这一节要说的比那更麻烦:有四种错不只是「不报错」,它们会让求解变得更容易——迭代更快收敛、云图更干净、机时更短。它们与人的动机同向。

这个错 它把结果推向哪 为什么求解反而更顺
载荷加小了(峰峰/零峰折反了,或漏掉压力波扫过管段时两端的不平衡力) 应力偏小 位移小、非线性弱,迭代次数下降
频段截窄了(耦合时间步取粗) 高频分量被数值阻尼吃掉,谱看起来干净 步长大 ⇒ 步数少 ⇒ 机时短
约束加死了(卡箍按理想固支建,而不是按有限刚度) f_n 抬高、位移变小 刚度矩阵条件数变好,收敛更快
流体按不可压处理 管系声学模态整个不存在 不解声学方程,压力场瞬时传遍全域

最后一行的后果最大,值得单说:不可压假设下声速趋于无穷,管系声学模态频率跟着趋于无穷 ⇒ 「激励阶次 ↔ 声学模态」「声学模态 ↔ 结构模态」这两对频率比对连算都不会算。屏幕上不会少任何东西,报告里也不会少任何一页——少掉的是两整条判据,而它们的缺席不会以任何形式出现在输出里。

⇒ 这四种错⛔ 不能靠「注意一下」来防——它们四个都在往省事和安全两个方向同时使劲。唯一挡得住的办法是把它们做成第 ⑥ 段交接时必须逐项回答的问句(1.7 的口径四件就是干这个的)。

另外两条相关的:

  • ⛔ 用稳态 CFD 结果做流致振动分析。稳态解里没有瞬态脉动,而脉动正是这门课全部的输入。稳态压力场映射过去能得到一个静应力分布,它是对的,但它回答的是耐压问题,⛔ 不是振动问题。
  • ⛔ 把「网格加密后结果不变」当成正确性证明。那证明的是网格无关性——它属于 CFD 方法本体,归 J2-01 CFD 基础、湍流模型与网格划分(Fluent/STAR-CCM+)。网格无关性说明「这个方程被解对了」,⛔ 不说明「这个方程的载荷项是对的」。载荷加小了的模型,加密网格之后还是那个加小了的载荷。

1.7 本课的入口与出口都被上游定死

流致振动的问题在项目上通常是这样落到你手里的:一个结果——异响、渗漏、焊缝开裂——被丢过来,然后链上五门课的人都会觉得该由 FSI 一次算清楚。

不划界的直接后果,是本课越界给出一个自己没有判据的寿命数或限值。所以先把两个接口写死。

中央大框写「本课:流固耦合与流致振动」,左侧三个入口箭头分别标「泵内部脉动的量级与阶次(来自 G5-07)」「几何改动与脉动之间的因果方向(来自 G5-05)」「汽蚀会不会发生、NPSH 裕量(来自 G5-04)」;右侧两个出口箭头标「σ_a 或应力 PSD + 四件口径 → J5-03」「频率带对账结论 → 设计改进」;大框外围一圈灰色小框列出十一条明确不在本课的内容与各自去向,每条各带红色禁止标记。本图为定性示意,⛔ 不得据图取值。
图3 本课在整条链上的位置:两个接口都不在本课内部产生。要读出三件:① 本课的输入不是自己产生的,三条入口各有属主课;② 本课的输出是一份带口径的应力交付物,⛔ 不是寿命、⛔ 不是许用值;③ 外圈十一条是明确不在本课的内容,遇到它们应当直接去对应的课,⛔ 不在本课里找答案。本图为定性示意,⛔ 不得据图取值。

入口:本课从「脉动/汽蚀已经算出来了」这一步接手。

入口 谁给 本课拿它干什么
泵内部脉动的产生机理、量级与阶次 G5-07 电子油泵本体零件级设计:容积式构型几何与排量、间隙链、困油与流量脉动 当作 ① 段的激励输入(p̃ 与 f_p)
几何改动与脉动之间的因果方向 G5-05 电子水泵水力部件设计:叶轮/蜗壳几何、内部间隙与轴向力 判「削源」这条抑制路有没有杠杆
汽蚀会不会发生、NPSH 裕量够不够 G5-04 泵的气蚀、空气分离与可靠性 判要不要把宽带高频那一路算进来

⛔ 本课不判汽蚀会不会发生,只做「已经发生的汽蚀激起了什么结构响应」。把汽蚀余量不足的问题当结构问题去加卡箍,是这条边界上最常见的一次误判——加了卡箍,异响可能真的小了,而气泡还在原地溃灭,叶轮该蚀还是蚀。

出口:本课交出两样,到此为止。

  1. 一次频率带对账——激励带、声学模态带、结构湿模态带三条带在同一根频率轴上重不重叠。
  2. 一个口径写死的应力量:确定性谐波激励交应力幅 σ_a;宽带随机激励(汽蚀、湍流脉动、塞状流冲击)交应力 PSD。两者的应力口径都是「关注点沿最大主应力方向的单轴分量」,与 J5-03 振动、模态与随机振动疲劳 的 σ_rms 口径对齐。⛔ 不得直接交 von Mises 而不声明等效方法。

交接必须附四件,缺一件那个数就不具备被判对错的资格:

件 要回答的问题
① 口径 是 σ_a 还是应力 PSD?应力取的是单轴分量还是等效应力?关注点在哪
② 频段与带宽 这个数对应哪一段频率?如果是窄带,中心频率与带宽各是多少
③ 是输入还是响应 这是加在结构上的载荷,还是结构算出来的响应
④ 工况时长与占比 这个工况在整个服役谱里占多长时间、多大比例

⛔ 本课到此为止:不产出寿命、不产出许用值、不产出安全系数。寿命计算归 J5-03 振动、模态与随机振动疲劳,许用值与安全系数归 J5-01 换热器/冷板结构强度与耐压 FEA,接头的应力评价口径归 J5-07 钎缝/焊缝接头的疲劳评价方法:三级应力口径、等效结构应力法与自建 S-N。

射程之外的十一件,只给去向:

你要的东西 去向
随机振动模型本体、PSD 输入与三区间判读、Miner 累积 J5-03 振动、模态与随机振动疲劳(本课只补「流体激励这一路」怎么进模型)
CFD 的湍流模型、y+ 与网格无关性 J2-01 CFD 基础、湍流模型与网格划分(Fluent/STAR-CCM+)
疲劳评价方法本体(应力口径三级、自建 S-N) J5-07 钎缝/焊缝接头的疲劳评价方法:三级应力口径、等效结构应力法与自建 S-N
热循环与应变寿命 J5-02 热应力、热变形与热疲劳仿真
载荷谱提取、雨流压缩、平均应力修正式本体 I6-01 振动、热冲击与压力脉冲耐久设计
台架执行与谱代表性 L1-04 振动、热冲击与耐久台架试验
许用值与安全系数 J5-01 换热器/冷板结构强度与耐压 FEA
隔振/减振件本身的设计与选型 I4-01 支架、悬置与减振结构设计
模态锤击与振动台标定的执行 L6-04 结构动态特性试验:模态、动刚度与传递率实测
气动噪声与结构声辐射的声学 CAE J5-05 声学 CAE:气动噪声与结构声辐射仿真;空调箱风噪与风门异响 C2-11 空调箱 NVH:风噪/风门异响/鼓风机噪声/再生噪声
阀切换速率的标定与一体化集成 G8-02 多通阀 + 泵 + 传感器的一体化集成设计;通用结构建模 J5-06 结构 FEA 建模方法基础:单元、约束、本构与结果可信度

标准。本课挂 ISO 16750-3(道路车辆电气电子设备的环境条件与试验,机械负荷部分)。本课未取得原文核对 ⇒ 只写去向不写内容:⛔ 不给任何限值、试验等级、谱形与时长;版次以发布方现行目录为准,引用前须核。另外要如实说明一件事:管路振动与噪声的评价没有统一编号的通用标准,各主机厂与 Tier1 各自定企业规范。→ 拿限值的路:查本平台的企业规范;许用值与安全系数归 J5-01 换热器/冷板结构强度与耐压 FEA,试验执行与谱代表性归 L1-04 振动、热冲击与耐久台架试验,应力评价口径归 J5-07 钎缝/焊缝接头的疲劳评价方法:三级应力口径、等效结构应力法与自建 S-N。

⚠ 比版次更容易错的是射程:先问「这条标准管的是不是这件事」。你拿到的任务书上写的可能是「按某标准做振动试验」,而你要评的那个件未必在那套体系的射程内——承压流体件走的是耐压、压力循环那一套。

三个易错点:① 拿 FSI 算出的振动幅值(mm/s 或 g)直接当疲劳输入——中间缺了应力提取与 J5-03 振动、模态与随机振动疲劳 的方法,速度和加速度不是应力;② 反过来,指望 J5-03 振动、模态与随机振动疲劳 替本课决定「激励谱长什么样」——随机振动那门课的输入是给定的谱,而流体激励这一路的谱正是本课的产出;③ 把「去向某课」读成「那门课会给我一个数」——去向给的是方法与口径,⛔ 不是数。

1.8 本讲的交付:三张卡

把前面七节压成三张可以直接贴进评审模板的东西。

卡一 · 符号口径卡(贴在报告首页,一次说清,后面不再解释)

  • 管内径写 D_i,绕流/障碍特征尺寸写 d_o;⛔ 禁止裸写 D,D_i ⛔ 不进 Strouhal 公式。
  • 声速裸写 c = 管内有效声速(含管壁柔度修正),分段取;自由声速另写 c_f。
  • 固有频率裸写 f_n = 湿模态;干模态写全 f_n,dry 并注明「上界」,⛔ 不进任何频率对账。
  • 脉动幅值 p̃ = 0-peak,单位 kPa;外部数据必须注明原口径与折算。
  • 长度写 L_span(卡箍跨距)或 L_ac(声学管段长度),⛔ 禁止裸写 L。
  • 放大因子写 Q_amp,流量写 q_v,⛔ 禁止裸写 Q;n 只表示谐次。
  • 结构侧 mm-t-s-MPa,流体侧 SI,接口 1 MPa = 1 N/mm² = 10⁶ Pa;⛔ 全课不写 bar。
  • 分析频段 0~1 kHz(疲劳),异响判读延伸到 2 kHz;以上归 J5-05 声学 CAE:气动噪声与结构声辐射仿真。

卡二 · 六段自检卡(用来评审任何一份流固耦合报告,含你自己的)

逐段问一句「这一段错了,谁会告诉我」,答案是「⛔ 没有人」的那五段,就是你必须自己动手核的地方。六个问句见 1.5,四种「让计算更顺同时让结论更安全」的错见 1.6。

⇒ 判出来之后的下一个动作:若六段里有任何一段答不上来,先补那一段,⛔ 不要先去加密网格或换更贵的耦合方式——那两件事一件都不会碰到这五段里的任何一个问题。

卡三 · 交接卡(交给做寿命的人时随数字一起交)

一行数字,四件限定:口径(σ_a 还是 PSD、单轴还是等效、关注点在哪)· 频段与带宽 · 是输入还是响应 · 工况时长与占比。

⇒ 收到一个不带这四件的应力数时,正确的动作是把它退回去,⛔ 不是自己替它补一个假设。四件里任缺一件,这个数都无法被判对错——它不是「不够准」,是没有定义。

⚠ 最后提醒一条本讲不给的东西:三条频率带(激励带、声学模态带、结构湿模态带)具体怎么算、粗筛带宽取多少、粗筛通过之后下一步做什么,都在第 3~5 讲。这里只钉死了口径——但口径钉不死,后面每一步算出来的数都不具备可比性。

后面还有 7 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做

会员专属

后续为会员深水区内容——四库数据与深度拆解。

查看会员方案