汽车气动噪声仿真入门：基于STAR-CCM+与VA One的联合仿真教学指南

starccm与Vaone进行汽车气动噪声仿真教学入门对做NVH入门的汽车人或者学生来说单独摸starccm算CFD过流体单独摸Vaone算SEA算结构但两个串起来搞「真实车外风噪传到人耳的全链路小闭环」总感觉是个「玄学连接点」卡脖子——今天就用最懒小白的视角把这套流程踩一遍代码尽量写得能直接偷咳参考分析也说人话不装。第一步starccm干的脏活累活——算湍流压力脉动表面声源别听学术圈说什么「分离涡DES」「延迟大涡IDDES」入门门槛高先上稳态RANSk-ω SST先算定常流场然后用瞬态的SAS或者LES LiteSST SAS更快小白选这个没问题算短时间窗口的脉动这个时间窗口得卡SEA的最低频对应的1/10周期不对别背公式直接选「覆盖100Hz-10kHz主要人耳敏感带」就行——100Hz的周期是10ms那取20ms的瞬态时间步长0.01ms别SAS建议CFL数库朗数控制在1-10跑起来顺就行。starccm与Vaone进行汽车气动噪声仿真教学入门先上starccm的基础启动和定常初始化的简单脚本别写太复杂的GUI操作步骤文字堆脚本复制进去点运行小白只会喊666// 打开一个空的工程记得替换成你自己的几何路径哈 newSession openGeometry D:/car_nvh/simple_sedan.stl // 画个外流场小白直接用圆柱盒前伸5车长、后伸15车长、宽/高各4车长的自动包络尺寸不够后面会教怎么调 createEnclosure Enclosure1 { type BOX frontOffset 5.0 * Length[Sedan] backOffset 15.0 * Length[Sedan] leftOffset 4.0 * Width[Sedan] rightOffset 4.0 * Width[Sedan] topOffset 4.0 * Height[Sedan] bottomOffset 1.0 * Height[Sedan] } // 选物理模型小白记住这套组合不可压缩空气、SST k-ω SAS、湍流壁面函数Y卡30-100就行不用卡到1的LES壁面处理、定常瞬态切换的准备 selectPhysicsModels Continuum1 { flow SegregatedFlow time Steady material Air viscous Turbulent turbulence SSTkOmega turbulenceModelOptions { SSTkOmega { curvatureCorrection ON SAS ON } } wallTreatment AllYPlus } // 边界条件入口120km/h换算成m/s就是33.333333别写整数显得专业点不是算的准一点出口压力0Pa大气压对称面Symmetry地面Ground设成移动壁面速度和入口一样方向相反车的外表面叫CarBody设成无滑移壁面 setBoundary Inlet1 { type VelocityInlet velocityMagnitude 33.3333333333 velocityDirection Vector(1, 0, 0) turbulenceSpecification IntensityAndLengthScale turbulenceIntensity 0.5 turbulenceLengthScale 0.01 * Length[Sedan] } setBoundary Outlet1 { type PressureOutlet staticPressure 0.0 } setBoundary Symmetry1 Symmetry2 Symmetry3 { type Symmetry } setBoundary Ground1 { type Wall wallMotion MovingWall wallVelocity Vector(-33.3333333333, 0, 0) } setBoundary CarBody { type Wall wallTemperature Constant temperature 300.0 } // 自动网格小白直接用Advancing Layer MeshSurface RemesherPolyhedral Mesher车外表面和地面附近加密这里偷个懒用自带的汽车边界层加密模板 createMeshOperation MeshOperation1 { baseMesh PolyhedralMesher surfaceMesher SurfaceRemesher prismMesher AdvancingLayerMesher targets [Enclosure1, Sedan] prismMesherOptions { AdvancingLayer { boundaryLayerTemplate AutomotiveBoundaryLayer targetPrismStretching 1.2 prismLayerThickness 0.005 * Length[Sedan] firstLayerThickness 0.0005 * Length[Sedan] } } surfaceMeshSizeOptions { SurfaceRemesher { baseSize 0.2 * Length[Sedan] relativeMinimumSize 0.01 relativeMaximumSize 2.0 curvatureControl ON curvatureAngle 10.0 proximityControl ON proximitySizeRatio 0.5 } } volumeMeshSizeOptions { Polyhedral { baseSize 0.2 * Length[Sedan] relativeMinimumSize 0.05 relativeMaximumSize 1.0 } } customControls [ CustomControl1 { type SurfaceControl surfaces [CarBody] surfaceMeshSizeOptions { SurfaceRemesher { baseSize 0.01 * Length[Sedan] curvatureAngle 5.0 } } prismMesherOptions { AdvancingLayer { targetPrismStretching 1.2 prismLayerThickness 0.005 * Length[Sedan] firstLayerThickness 0.0005 * Length[Sedan] numberOfLayers 10 } } }, CustomControl2 { type VolumeControl shape Cylinder center Vector(0, 0, Height[Sedan]/2) direction Vector(1, 0, 0) radius 2.0 * Width[Sedan] length 20.0 * Length[Sedan] volumeMeshSizeOptions { Polyhedral { baseSize 0.05 * Length[Sedan] } } } ] } // 初始化定常计算到残差收敛动量、连续性1e-4k、ω 1e-5监视器比如后视镜的升力系数稳定就行 initialize Continuum1 setSolverParameters Continuum1 { segregatedFlow { maxIterations 2000 } } createMonitor Monitor1 { type ForceCoefficient surfaces [MirrorLeft, MirrorRight] direction Vector(0, 1, 0) referenceArea 0.5 * Width[Sedan] * Height[Sedan] referenceLength Length[Sedan] referenceVelocity 33.3333333333 } createPlot Plot1 { monitors [Monitor1] xAxis Iteration } run Continuum1这段代码跑通的话你会得到一个稳定的后视镜涡脱落的定常流场残差应该也下来了。接下来是瞬态SAS开算保存CarBody的脉动压力这个是starccm给Vaone的「核心原料」// 把物理模型从Steady改成Transient setPhysicsModels Continuum1 { time Transient unsteadyOptions { Transient { timeStep 0.0001 numberOfTimeSteps 200 } } } // 重置迭代步用定常结果初始化瞬态 reset Continuum1 initialize Continuum1 // 监视器别忘了改成Time updatePlot Plot1 { xAxis Time } // 最关键的一步保存CarBody的Pressure Fluctuation瞬态压力减去定常平均压力保存成Vaone能读的「Star-CCM Surface Pressure Files」格式文件名记得带面的信息窗口选最后100个时间步前100个是瞬态过渡噪声不准 createExport SurfacePressureExport1 { type StarCCMSurfacePressure surfaces [CarBody] windowStart 0.01 windowEnd 0.02 timeStep 0.0001 filePath D:/car_nvh/export/CarBody_Pressure_Fluctuation.ssp subtractMeanPressure ON useLocalReferenceFrame OFF } run Continuum1这里插一句人话减平均压力太重要了不然你传给Vaone的是带几百Pa大气压波动的垃圾噪声级会直接上天。时间步长0.0001s对应奈奎斯特频率5000Hz不对是采样频率10000Hz对应奈奎斯特5000但我们要的是10kHz所以可以把时间步长改成0.00005s窗口步数改成200对应0.01s前100步还是过渡这样采样频率20000Hz奈奎斯特10kHz刚好覆盖主要人耳敏感带——奈奎斯特是什么简单说就是「你要测10kHz的信号必须每秒至少采20000次不然会混频」小白记这个就行不用管傅里叶变换的数学。第二步Vaone干的精细活——SEA建模耦合声源算人耳声压级starccm搞完了「车身上哪里响」接下来Vaone要搞「这些响声怎么通过车身板件传到车内再传到驾驶员右耳」——小白别一开始就搞复杂的整车SEA模型先上简化的四门两盖地板顶棚前围板的10个左右腔体板件的模型车门内饰、顶棚内饰这些用声学阻尼层代替就行。首先是Vaone的建模小技巧不用脚本因为Vaone的GUI建模逻辑比starccm更顺但可以提几个快捷键或者快速操作导入几何直接导入刚才的simple_sedan.stl但别忘了「Simplify Geometry」把所有小倒角、后视镜这些不对后视镜刚才导出了表面压力所以要保留哦对刚才的surface export里要把MirrorLeft、MirrorRight单独选出来吗不用CarBody已经包含了导入Vaone之后把CarBody拆分成几个声学板块就行——前围板、地板、顶棚、左前门、右前门、左后门、右后门、后备箱盖、前引擎盖、左后视镜、右后视镜、A/B/C柱外板。快速创建腔体和板块用「Create Cavity by Enclosed Surfaces」选中左前门板、前围板左半、地板左半、顶棚左半、左A/B/C柱外板自动创建左前腔体驾驶员的位置同理创建右前、左后、右后、后备箱、发动机舱、前风挡玻璃下方的雨刮腔如果有的话小白可以省略。设置板块属性所有外板用「Steel 0.8mm」厚度自己改根据你的几何来内饰用「Acoustic Foam 25mm NRC 0.8」小白可以直接用Vaone自带的材料库不用自己算杨氏模量、损耗因子这些。设置耦合外板和腔体之间自动会有「Structure-Acoustic Coupling」不用手动加。接下来是最激动人心的一步把starccm的表面声源导入Vaone操作步骤点击「Sources」→「Import」→「Star-CCM Surface Pressure Files」选中刚才导出的CarBodyPressureFluctuation.ssp。Vaone会自动识别所有导出的面你只需要把「Source Type」改成「Turbulent Boundary Layer Pressure Fluctuation (Corcos Model)」不对不对刚才我们导出的是真实的瞬态压力脉动不是用Corcos模型算的所以应该改成「Measured/Simulated Surface Pressure」Vaone里好像叫「User-Defined Surface Pressure Source」反正就是不用默认的Corcos用我们自己sim的就行。把每个声源面和对应的Vaone板块关联起来——比如导入的CarBody1是左前门板就关联到Vaone的LeftFrontDoorPlateCarBody2是右后视镜关联到RightMirrorPlate以此类推。点击「Apply」声源就加进去了。然后是设置人耳接收点点击「Receivers」→「Create」→「Point Receiver」坐标设成驾驶员右耳的位置比如(1.0, 0.3, 1.2)单位m根据你的几何来。把接收点放在左前腔体内自动会有「Acoustic Receiver」。最后是计算看结果点击「Analysis」→「Frequency Response Analysis」频率范围选「100Hz-10kHz」步长选「1/3倍频程」NVH行业标准小白不用改。点击「Run」计算速度取决于你的电脑配置简化模型大概10分钟左右。看结果点击「Results」→「Frequency Response」选中接收点和左前腔体的「Sound Pressure Level (SPL)」你会看到一条从100Hz到10kHz的曲线大概在2000Hz-5000Hz会有一个峰值因为后视镜涡脱落的频率大概在这个范围小白可以对比一下starccm里的Monitor1升力系数的FFT应该能对应上。踩坑总结小白必看不然会哭starccm里的时间窗口和采样频率刚才说了别忘减平均压力别忘采样频率是最高频率的2倍以上别忘前一半时间步是过渡别导出进去。Vaone里的面关联一定要把导入的starccm声源面和Vaone的板块一一对应不然声源会加在错误的地方结果全错。Vaone里的材料属性别随便改材料库的参数小白先用默认的等入门了再自己调损耗因子、阻尼系数这些。简化模型的重要性一开始别搞复杂的整车模型简化到10个左右腔体板块就行不然计算速度慢到爆炸而且你不知道哪里错了。今天的入门就到这里啦这套流程虽然简单但已经覆盖了汽车气动噪声仿真的核心逻辑CFD算湍流压力脉动表面声源→SEA算结构-声学耦合传播路径→算接收点声压级结果。下次可以讲讲怎么优化这个模型比如把前风挡玻璃、车门玻璃换成透明的声学材料或者把Corcos模型和真实的瞬态压力脉动对比一下。如果有哪里看不懂的欢迎在评论区留言哦