张克鹏_基于Virtual Wind Tunnel的某重型牵引车外气动性能分析
风屏障对桥上高速列车气动性能影响的数值研究_项超群
1 1 4
中 国 铁 道 科 学 第 3 5卷
题的数值求解需要消耗大量的计算机资源 , 尤其是 非定常求解的动网格算法在每一时间步内需要进行 反复迭代 , 其计算效率很低 。 因此考虑列车速度对 气动力影响时 , 已有研究多采用静网格合成风的方 法 , 即将列车运行速度的反向矢量与侧风进行矢量 合成 , 将合成得到的矢量作为边界条件 , 而列车外 壁作为流场中的壁面其速度设置 为 0。 在 没 有 风 屏 障阻 挡 时 可 以 近 似 认 为 此 假 设 是 成 立 的 , 文 献 ] 通过动静 网 格 模 型 的 对 比 分 析 , 发 现 静 网 格 [ 3 1 算法和动网格算 法 计 算 结 果 相 差 不 大 , 文 献 [ 4] 1 通过试验也证明了此假设基本可靠 。 侧风作用下运 行列车受到的气动力可以看成由两部分组成 , 即列 车风和自然界的侧风 。 当桥上设置风屏障后 , 风屏 障对外界侧风有阻挡作用 , 而由于列车风产生于主 梁和两侧风屏 障 组 成 的 U 型 槽 内 , 因 此 风 屏 障 不 阻挡其作用 , 但采用一般的合成风计算时 , 风屏障 阻挡了合成风速矢量 , 因此在有风屏障情况下合成 风的计算方法有待进一步验证 。 风屏障对列车周围 流场分布的影响研究一般基于合成风方法 , 采用动 网格算法的研究还不多见 。 静网格合成风方法和能模拟列车真实运动的动 网格算法计算机理不同 , 为了从流场特性方面研究 二者 的 区 别 , 本 文 基 于 大 型 流 体 计 算 通 用 软 件 L U E NT 对简支箱 梁 桥 上 不 同 风 屏 障 高 度 下 列 车 F 的外流场进行数值模拟 , 得出列车气动力系数随风 屏障高度的变化规律 , 并对采用动网格和静网格算 法得到的静压分布和速度分布进行对比研究 。
天津大学流体力学实验室流体力学实验指导书(简写本)
实验 10 平板湍流边界层雷诺应力分布的测量(选作)............................................. 34 实验 11 平板湍流边界层壁面摩擦阻力的测量(选作)............................................. 37 实验 12 氢气泡流动显示实验.................................................................................... 40 实验 13 水槽中湍流边界层瞬态流动结构的 TRPIV 测量 ...................................... 45 实验 14 圆柱与非圆截面柱体尾涡结构的流动显示................................................ 48 实验 15 用热线风速仪测量圆柱尾流平均速度分布................................................ 49 实验 16 用热线风速仪测量圆柱尾流湍流度和雷诺应力分布................................ 51 实验 17 水槽中圆柱尾流的 TRPIV 测量 .................................................................. 53 实验 18 弯道压力分布实验........................................................................................ 55 实验 19 半圆形防波堤分离涡的流动显示与表面压力分布的测量........................ 58 实验 20 各种管道的流动显示.................................................................................... 62 实验 21 加齿被动控制射流的流动显示.................................................................... 69 笔 记............................................................................................................................. 73
风屏障对公铁两用桥上车辆气动特性的影响
风屏障对公铁两用桥上车辆气动特性的影响甄谢友;苏洋;向活跃;李永乐【摘要】以某公铁两用桥为研究背景,通过大比尺节段模型风洞试验,使用天平测试有无风屏障时公路和铁路车辆气动特性,采用风速仪测试了桥面的风剖面分布,研究了车道和车辆类型对公路桥面车辆气动特性的影响.结果表明:设置风屏障有效降低了公路和铁路桥面的局部风速和车辆的气动力系数;公路桥面车辆气动力系数总体上随车道距风屏障距离的增加而减小,相同风屏障对大货车气动特性的降低程度相对于小货车和客车更为明显;设置风屏障后铁路桥面迎风侧和背风侧列车阻力系数的折减率基本接近,但升力系数的折减率差异较大.【期刊名称】《铁道建筑》【年(卷),期】2018(058)005【总页数】6页(P25-30)【关键词】公铁两用桁架桥;风屏障;气动力系数;风洞试验;车辆气动力;桥面局部风场【作者】甄谢友;苏洋;向活跃;李永乐【作者单位】西南交通大学桥梁工程系,四川成都610031;西南交通大学桥梁工程系,四川成都610031;西南交通大学桥梁工程系,四川成都610031;西南交通大学桥梁工程系,四川成都610031【正文语种】中文【中图分类】U448.12+1为提高桥梁通行能力,缓解交通压力,大跨度跨江跨海大桥常采用公铁两用的结构形式。
以钢桁架为主体框架,桥面采用整体性钢桥面板结构,桥面板与下弦杆连为整体,这种主梁结构具有较好的整体性、空间利用能力和合理受力性能,成为公铁两用桥梁首选形式。
如主跨312 m的芜湖长江公铁两用斜拉桥[1]和主跨504 m 的武汉天兴洲长江公铁两用斜拉桥[2]。
另一方面,公铁两用钢桁架桥断面较为复杂,上下2个桥面相互影响,在侧风作用下当车辆通过时流场变化显著。
因此为了提高车辆的行车安全性,桥面上通常要设置风屏障。
近年来,国内外学者对风屏障的防风效果进行了广泛研究。
在国内,白璐等[3]采用CFD数值模拟方法,研究防风栅的设置位置对小轿车、大客车和大货车的气动特性的影响,提出防风栅的最佳设置方案,并推荐了设计临界风速。
高速动车组开闭罩缝隙对气动性能的影响研究
基于CFD
第8卷㊀第6期2023年11月气体物理PHYSICSOFGASESVol.8㊀No.6Nov.2023㊀㊀DOI:10.19527/j.cnki.2096 ̄1642.1088基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析李泳德ꎬ㊀郭㊀力ꎬ㊀季㊀辰(中国航天空气动力技术研究院ꎬ北京100074)CharacterizationofTransonicAerodynamicDampingofRocketsBasedonCFD/CSDCouplingLIYong ̄deꎬ㊀GUOLiꎬ㊀JIChen(ChinaAcademyofAerospaceAerodynamicsꎬBeijing100074ꎬChina)摘㊀要:随着新型大推力火箭的发展ꎬ弯曲模态频率的不断降低ꎬ以及流动分离和跨声速飞行时产生的激波震荡等因素ꎬ其在跨声速飞行过程中更容易出现非定常振动发散ꎮ文章以某带助推的运载火箭模型为研究对象ꎬ通过数值计算获取火箭强迫振动时的气动阻尼ꎬ并对影响火箭气动阻尼的因素进行了分析ꎮ包括结构节点位置㊁振动振幅大小㊁脉动压力等ꎮ研究表明:助推主要起到增大气动阻尼的作用ꎻ前节点主要影响收缩段的气动阻尼ꎻ振动振幅大小和脉动压力对气动阻尼的影响可忽略不计ꎮ关键词:气动阻尼ꎻ数值计算ꎻ跨声速ꎻ气动弹性ꎻ运载火箭㊀㊀㊀收稿日期:2023 ̄09 ̄25ꎻ修回日期:2023 ̄10 ̄23第一作者简介:李泳德(1995 ̄)㊀男ꎬ工学硕士ꎬ助理工程师ꎬ主要研究方向为气动弹性分析ꎮE ̄mail:562064169@qq.com通信作者简介:季辰(1982 ̄)㊀男ꎬ工学博士ꎬ研究员ꎬ主要研究方向为气动弹性力学ꎮE ̄mail:jichen167@hotmail.com中图分类号:V475.1㊀㊀文献标志码:AAbstract:Withthedevelopmentofnewhigh ̄thrustrocketsꎬthedecreasingfrequencyofthebendingmodesoftherocketꎬaswellasthefactorssuchasflowseparationandshockoscillationsgeneratedduringtransonicflightmakeitmorepronetonon ̄constantvibration.Inthispaperꎬalaunchvehiclemodelwithboostwastakenastheresearchobjectꎬandtheaerody ̄namicdampingoftherocketduringforcedvibrationwasobtainedthroughnumericalcalculation.Thefactorsaffectingtheaerodynamicdampingoftherocketwereanalyzedꎬincludingthepositionofstructuralnodesꎬthemagnitudeofvibrationam ̄plitudeꎬpulsatingpressureandsoon.Thestudyshowsthattheboostmainlyplaystheroleofincreasingaerodynamicdamp ̄ingandthefrontnodemainlyaffectstheaerodynamicdampingofthecontractionsection.Thevibrationamplitudesizeandthepulsatingpressurehaveanegligibleeffectontheaerodynamicdamping.Keywords:aerodynamicdampingꎻnumericalcalculationꎻtransonicꎻaeroelasticityꎻlaunchvehicle引㊀言通常情况下人们认为气动力对火箭的振动起到阻尼作用ꎬ即气动阻尼为正值ꎮ然而随着大推力火箭发展ꎬ火箭的长细比逐渐加大ꎬ导致弯曲刚度越来越小ꎬ同时为了满足有效载荷的外形要求ꎬ火箭头部整流罩尺寸不断加大ꎬ后续箱体的直径却保持不变ꎬ形成了典型的锤头体外形ꎮ国内外大量的火箭研制经验表明[1 ̄9]ꎬ对于此类锤头体外形火箭的气动设计ꎬ必须要进行动态气动载荷与动态气弹稳定性分析ꎬ否则设计的疏忽可能会导致火箭结构出现毁灭性的破坏进而导致发射失败ꎮ目前常用的衡量气弹稳定性的方法是通过风洞试验来获取气动阻尼系数ꎮ早在1963年ꎬ美国国家航空航天局Ames研究中心(NASAAmesRe ̄searchCenter)采用半刚性模型开展试验研究[10]ꎬ获取火箭头部的气动阻尼来评估其稳定性ꎬ但这只能用来模拟火箭弯曲振型前节点之前部分的结构动力学特性ꎮ直到兰利研究中心(NASALangleyResearchCenter)开发了全弹性模型气动阻尼试验气体物理2023年㊀第8卷技术ꎬ其可以模拟整体的结构动力学特性以及气动外形ꎬ并应用于多款运载火箭研制[11 ̄15]ꎮ国内ꎬ中国航天空气动力技术研究院对气动阻尼问题开展过较多的研究[16 ̄20]ꎬ从模型设计方法㊁模型制作工艺㊁试验机构设计和数据处理等诸多方面ꎬ逐步改进实现了从半刚性模型到全弹性模型的过渡ꎬ并在多个型号上得到验证ꎮ然而通过风洞试验研究气动弹性问题ꎬ技术难度大ꎬ试验成本高ꎬ同时几乎不可能开展全尺寸试验ꎮ因此通过数值计算的方法开展相关研究是另一种重要的手段ꎮ刘子强等[21]实现了通过数值计算确定气动阻尼系数的技术和方法ꎬ并与试验结果进行对比ꎬ证实了该方法的可靠性ꎮ冉景洪等[22]通过模态数据结合准定常理论的方法分析了减阻杆加后体这一弹性结构的气动阻尼ꎬ结果表明减阻杆造成的分离流会对后体的气动阻尼系数产生影响ꎮ朱剑等[23]针对新一代捆绑式运载火箭发展了非结构网格下的气动阻尼计算方法ꎬ并分析了攻角㊁Mach数等参数对气动阻尼的影响ꎮ本文在之前的计算方法[23]的基础上采用IDDES模型ꎬ考虑脉动压力的影响ꎬ通过强迫振动的方式ꎬ针对捆绑式运载火箭的某一特定模态进行数值计算仿真ꎬ研究前节点位置ꎬ振动振幅ꎬ脉动压力等参数对气动阻尼的影响规律ꎮ1㊀计算方法图1为本文所用的捆绑式运载火箭的计算模型ꎬ是典型的锤头体结构ꎮ在跨声速阶段ꎬ其头部会产生激波造成激波边界层干扰ꎬ而在锤头体外形的过渡段会出现气流分离ꎮ为探究各部分气动阻尼的变化ꎬ将整个箭体分为头部㊁过渡段㊁弹身3个部分ꎮ图1㊀表面网格及区域划分Fig.1㊀Surfacegridandregiondivision1.1㊀流场仿真模型本文分别用Reynolds平均法(Reynolds ̄averagedNavier ̄StokesꎬRANS)和改进的延迟分离涡模拟(improveddelayeddetached ̄eddysimulationꎬID ̄DES)[24 ̄25]进行计算ꎬ在RANS方程中ꎬ将变量分为平均值和波动值两部分ꎬ对于速度分量有ui=ui+uᶄi其中ꎬi=1ꎬ2ꎬ3ꎬui和uᶄi分别代表平均量和波动量ꎬ对于压强和其他标量也采用类似的形式ꎬ将这种形式代入连续性方程和动量方程中ꎬ并写成张量形式∂ρ∂t+∂∂xi(ρui)=0(1)∂∂t(ρui)+∂∂xj(ρuiuj)=∂p∂xi+∂∂xjμ∂ui∂xj+∂uj∂xi-23δij∂uk∂xkæèçöø÷éëêêùûúú+∂∂xj(-ρuᶄiuᶄj)(2)其中ꎬiꎬjꎬk可分别取1ꎬ2ꎬ3ꎻρ是密度ꎻt是时间ꎻ当i=j时δij取0ꎬ否则取1ꎮ式(1)㊁(2)是RANS方程ꎬ由方程可知RANS方法将湍流脉动对平均流动的作用模化为Reynolds应力项即-ρuᶄiuᶄjꎬ之后采用湍流模型进行封闭ꎬ本文采用的湍流模型为SSTk ̄ω模型ꎬ其输运方程为∂∂t(ρk)+∂∂xi(ρkui)=∂∂xjΓk∂k∂xjæèçöø÷+Gk-Yk∂∂t(ρω)+∂∂xi(ρωui)=∂∂xjΓω∂ω∂xjæèçöø÷+Gω-Yω其中ꎬk和ω分别代表湍流动能和湍流耗散率ꎬΓk和Γω分别代表k和ω的有效扩散系数ꎬGk和Gω分别代表k和ω的生成率ꎬYk和Yω分别代表k和ω的耗散率ꎮ因此RANS方法只能计算大尺度的平均流动ꎬ本文采用IDDES方法计算脉动压力对气动阻尼的影响ꎮIDDES方法是由分离涡模拟(detached ̄eddysimulationꎬDES)方法改进而来ꎬ其本质思想与DES方法相同ꎬ是想以网格尺度和模型中的特征尺度隐式划分RANS和大涡模拟(large ̄eddysimulationꎬLES)区域ꎬ使其既能处理RANS方法无法得到的脉动场ꎬ也能降低LES方法在模拟高Reynolds数流动时所需的计算资源ꎮ区别在于当边界层较厚或者分离区域较窄时ꎬDES方法会出现如模型应力损耗(modeledstressdepletionꎬMSD)ꎬ网格诱导分离(grid ̄inducedseparationꎬGIS)以及对数层不匹配(logarithmic ̄layermismatchꎬLLM)问题[24]ꎬ而IDDES模型通过改良计算区域划分ꎬ结合延迟分离涡模拟(delayeddetached ̄eddysimulationꎬDDES)和03第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析壁面模型大涡模拟(wall ̄modeledlarge ̄eddysimula ̄tionꎬWMLES)ꎬ定义新的长度尺度解决了这些问题ꎬ具体公式详见文献[25]ꎮ流场网格如图2㊁图3所示ꎬ边界层采用棱柱层结构ꎬ并调整第1层网格高度使得y+小于1ꎬ远场部分采用六面体结构网格ꎬ与边界层的过渡层采用非结构网格ꎮ整体网格单元数量为4.2ˑ106ꎮ图2㊀y方向截面网格示意图Fig.2㊀Schematicdiagramofcross ̄sectionalgridinthey ̄direction图3㊀x方向截面网格示意图Fig.3㊀Schematicdiagramofcross ̄sectionalgridinthex ̄direction物面边界条件为无滑移壁面条件ꎬ远场采用压力远场边界条件ꎬ湍流模型采用SSTk ̄ω模型ꎬ采用密度基求解ꎬ气体黏性采用Sutherland定律ꎬ空间离散采用2阶迎风格式ꎬ对流通量采用Roe格式ꎮ1.2 结构分析模型结构与流场耦合分析过程中ꎬ结构部分可以采用模态方法描述ꎮ结构模态可以通过有限元方法与结构模态试验方法获得ꎮ本文采用有限元分析结果获得的模态ꎬ图4所示为结构的前3阶模态ꎬ本文只分析计算结果中气动阻尼最小的第2阶模态ꎮ(a)f=1.200Hz(b)f=2.460Hz(c)f=2.957Hz图4㊀结构的前3阶模态Fig.4㊀Firstthreemodesofthestructure由于火箭结构外形简单ꎬ一般不考虑其扭转影响ꎬ因此可以将其简化为简单的梁模型ꎬ这样就可以给出其模态振动方程q㊆i+2biωiq˙i+ω2iqi=fi(3)式中ꎬqi为第i阶模态的广义位移ꎬbi为第i阶模态的结构阻尼系数ꎬωi为第i阶模态的固有频率ꎬ13气体物理2023年㊀第8卷fi为第i阶模态下质量归一化的广义气动力ꎮ若将fi按照Taylor展开并略去高阶项ꎬ可以将其转化为气动阻尼项与气动刚度项的形式ꎬ则式(3)可写为q㊆i+2(bi+Bi)ωiq˙i+(Ki+1)ω2iqi=0(4)式中ꎬBi为气动阻尼系数ꎬKi为气动刚度系数ꎬ研究表明[26]ꎬ气动刚度相对于结构刚度为小量可以忽略不计ꎬ而在计算中结构阻尼往往设置为0ꎬ因此气动阻尼可以直接反映其气弹稳定性ꎮ1.3㊀气动阻尼分析原理气动阻尼的分析可以采用强迫振动或者自由振动的方式进行ꎬ这两种方法获得的时域数据不同ꎬ提取气动阻尼的方式也不同ꎮ强迫振动方法初始演化过程较短ꎬ因此计算量较小ꎬ同时能够分析某一种振动形式的气动阻尼ꎬ明确该振动形式是收敛还是发散ꎮ分析过程中能够获得不同部位与部件的气动阻尼ꎮ但是对于多模态相互作用引起的发散(例如颤振)较难预测ꎮ自由振动方法需要一定的自由演化时间才能够对时域数据进行分析ꎬ不过自由振动方法能够获得最能够吸收能量的模态及其振动频率ꎮ对于本研究所关注的问题ꎬ气动载荷对结构振动的过程中气动阻尼的影响较大ꎬ而对气动刚度与气动质量影响较小ꎬ即结构的固有振动频率受到来流的影响较小ꎬ其稳定性问题主要由气动阻尼的正㊁负引起ꎬ所以采用强迫振动方法分析ꎮ强迫振动下结构做简谐模态振动qi(t)=Asin(ωit)式中ꎬA表示振动的振幅ꎬ将其代入计算气动力的公式中[21]并做正交积分可得Bi=ʏl0Bx(x)dx=-1MiAω2iTʏl0ʏt0+Tt0G(xꎬt)cos(ωit)dtdx(5)式中ꎬMi为第i阶模态的模态质量ꎬT为整数倍周期ꎬG为广义气动力ꎮ根据式(5)便可以得到局部或分区域的气动阻尼ꎮ1.4㊀耦合计算流程首先进行模态分析ꎬ以确定结构的模态频率与振型ꎬ用以设计强迫振动的频率和振幅ꎮ非定常流场计算前先进行定常流场计算ꎬ来加快非定常计算的演化速度并增强收敛性ꎬ结构节点位移通过径向基函数(RBF)插值方法[27]映射到气动网格节点上ꎬ来进行网格的变形ꎬ这里径向基函数选用WendlandC2ꎬ如下所示φ(x)=(1-x)4(4x+1)最后将计算出来的广义力提取出来ꎬ截取演化完毕的整数倍周期ꎬ进行气动阻尼计算ꎮ耦合计算流程图如图5所示ꎮ图5㊀耦合计算流程图Fig.5㊀Flowchartofcoupledcalculation2㊀结果分析与讨论2.1㊀流场分析结果计算的来流Mach数范围为0.7~1.2ꎮ其中中截面的压力分布如图6所示ꎮ可以看出在头部出现了膨胀波以及跨声速激波ꎬ在过渡段存在流动分离ꎬ随着Mach数的增大ꎬ头部低压区域逐渐扩张ꎬ并且能明显看到ꎬ在流动再附的位置产生了再附激波ꎮ(a)Ma=0.7023第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析(b)Ma=0.75(c)Ma=0.80(d)Ma=0.85(e)Ma=0.88(f)Ma=0.90(g)Ma=0.92(h)Ma=0.96(i)Ma=0.9833气体物理2023年㊀第8卷(j)Ma=1.00(k)Ma=1.05(l)Ma=1.10图6㊀不同Mach数下的中截面压力分布Fig.6㊀PressuredistributioninthemiddlesectionatdifferentMachnumbers2.2 气动阻尼分布通过上述流场分析ꎬ可以看出火箭不同部位流动结构并不相同ꎬ在头部与箭身上ꎬ流动主要为附着流动ꎬ而在过渡段会出现较为复杂的波系结构以及流动分离ꎮ针对不同的流动结构随流向站位x的变化ꎬ设该位置上广义力与广义位移的相位差为φ(x)ꎬ并且简谐振动没有引入其他模态的广义力ꎬ则广义力的表达式为G(xꎬt)=Fgen sin[ωt+φ(x)]+F0(6)其中ꎬFgen为广义力的振动幅度ꎬF0为广义力的常数偏移量ꎮ将式(6)代入到式(5)中得到B(x)=-FgenMAω2Tʏt0+Tt0sin[ωt+φ(x)]cos(ωt)dt其中ꎬ广义力的常数偏移量F0的积分为0ꎬ因此省略ꎮ通过将等式中的正弦函数部分进行和差化积得到B(x)=-FgenMAω2Tʏt0+Tt0sin(ωt)cos[φ(x)]cos(ωt)dt+[ʏt0+Tt0sin[φ(x)]cos(ωt)cos(ωt)dt](7)式(7)中第1部分在整个周期中的积分为0ꎬ只有第2部分保留ꎬ因此得到B(x)=-Fgensin[φ(x)]MAω2Tʏt0+Tt0cos2(ωt)dt(8)式(8)中积分部分恒为正值ꎬ决定整个气动阻尼的部分只有相位角φ(x)的正弦值sin[φ(x)]ꎬ为了能够更加直观地获得相位角与气动阻尼B之间的关系ꎬ须将符号转化为对应的正弦函数转角ꎬ根据正弦关系ꎬ此转角为πꎬ因此得到B(x)=-Fgen(x)sin[φ(x)+π]MAω2Tʏt0+Tt0cos2(ωt)dt(9)图7为气动阻尼变化曲线ꎬ可以看出随着Mach数的增大ꎬ整体气动阻尼先增大后减少ꎬ在Mach数为0.98时达到最大值ꎬ过渡段与箭体的气动阻尼变化趋势与整体基本相同ꎬ而头部区域则不同ꎬ是随着Mach数的增大一直增大ꎬ只是增长速率变缓ꎮ图7㊀有助推时气动阻尼变化曲线Fig.7㊀Aerodynamicdampingchangecurvewithboost根据式(9)ꎬ得到相位角与气动阻尼B之间的关系为:当φ(x)ɪ(-πꎬ0)时ꎬ相位角滞后ꎬ气动阻尼B为负值ꎻ当φ(x)ɪ(0ꎬπ)ꎬ相位角提前ꎬ43第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析气动阻尼B为正值ꎻ为当φ(x)=0时ꎬ无相位角差别ꎬ气动阻尼B为0ꎮ在过渡段上ꎬ复杂的波系结构以及流动分离ꎬ使得气动力与结构位移之间会出现较为明显的迟滞现象ꎬ从而导致相位角φ(x)ɪ(-πꎬ0)ꎬ由此在过渡段上产生了负的气动阻尼ꎮ计算过程中的广义力与广义位移随时间变化曲线如图8所示ꎬ可以看出所有工况计算结果都表现良好ꎬ需要注意的是在非定常计算初期ꎬ演化的不完全导致广义力存在一些突变异常的结果ꎬ计算气动阻尼时须剔除ꎬ选择后面演化完全的周期ꎮ本文计算了9个周期ꎬ剔除了第1个周期出现的错误结果ꎬ采用后8个周期进行气动阻尼分析ꎮ强迫运动振幅为芯级直径的0.5%ꎮ(a)Ma=0.70㊀㊀㊀(b)Ma=0.75(c)Ma=0.80㊀㊀㊀(d)Ma=0.85(e)Ma=0.88㊀㊀㊀(f)Ma=0.9053气体物理2023年㊀第8卷(g)Ma=0.92㊀㊀㊀(h)Ma=0.96(i)Ma=0.98㊀㊀㊀(j)Ma=1.00(k)Ma=1.05㊀㊀㊀(l)Ma=1.10图8㊀不同工况下的广义力与广义位移随时间变化曲线Fig.8㊀Timedependentcurvesofgeneralizedforceandgeneralizeddisplacementunderdifferentoperatingconditions2.3㊀气动阻尼影响因素2.3.1㊀有无助推对气动阻尼的影响捆绑式运载火箭相比于传统的运载火箭ꎬ最大的区别就是在尾部四周捆绑了助推器ꎬ使得其流场特性变得复杂ꎬ因此须分析其对气动阻尼的影响ꎮ图7㊁图9分别为有无助推时气动阻尼变化曲线ꎬ可以看出随着Mach数的增大整体气动阻尼先增大后减少ꎬ在Mach数为0.98时达到最大值ꎬ过63第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析渡段与箭体的气动阻尼变化趋势与整体基本相同ꎬ而头部区域则不同ꎬ是随着Mach数的增大一直增大ꎬ只是增长速率变缓ꎮ对比两个图可知ꎬ助推主要起增大气动阻尼的作用ꎮ还可以看出有无助推情况下头部的气动阻尼变化很小ꎬ意味着在箭体尾部施加控制很难影响到头部的气动阻尼ꎬ特别是在超声速流场中ꎮ图9㊀无助推时气动阻尼变化曲线Fig.9㊀Aerodynamicdampingchangecurvewithoutboost2.3.2㊀前节点位置影响为了考察前节点位置变化对气动阻尼的影响ꎬ在保持振动频率不变㊁头部最大振型位置与振幅不变的条件下移动前节点ꎬ变化后的振型如图10所示ꎮ(a)Frontnodeafterthetransitionregion(b)Frontnodeinthetransitionregion(c)Frontnodebeforethetransitionregion图10㊀前节点变化后的振型Fig.10㊀Vibrationmodeafterthechangeofformernode根据对计算结果的分析分别获得了不同前节点位置的整体气动阻尼对比与过渡段气动阻尼对比ꎬ如图11㊁图12所示ꎬ可以看出前节点位置的改变并没有影响整体气动阻尼随Mach数增大而增大的趋势ꎬ且前节点在过渡段上与过渡段前的整体气动阻尼相差不大ꎬ而前节点在过渡段后的整体气动阻尼要高于另两种情况ꎬ因此过渡段与头部放在同一侧有助于提高气动阻尼ꎮ过渡段的气动阻尼会随着前节点的变化发生剧烈改变ꎬ前节点在过渡段前后随Mach数增大的变化规律相反ꎬ节点前后的振动相位变化导致不同节点位置过渡段的振动相位不同ꎬ进而导致气动阻尼发生变化ꎮ图11㊀不同节点位置的整体气动阻尼Fig.11㊀Overallaerodynamicdampingatdifferentnodepositions图12㊀不同节点位置的过渡段气动阻尼Fig.12㊀Aerodynamicdampingofthetransitionregionatdifferentnodepositions2.3.3㊀强迫振动振幅大小对气动阻尼的影响为了考察强迫振动振幅大小对气动阻尼的影响ꎬ在保证流场结构不发生改变的前提下ꎬ振动振幅分别为原来的一半和两倍ꎬ根据工程经验ꎬ如果振幅超过芯级直径的5%ꎬ则须考虑流场结构改变所造成的影响ꎮ图13㊁图14分别为不同振幅下的整体与头部气动阻尼ꎮ73气体物理2023年㊀第8卷图13㊀不同振幅下整体气动阻尼Fig.13㊀Overallaerodynamicdampingatdifferentamplitudes图14㊀不同振幅下头部气动阻尼Fig.14㊀Aerodynamicdampingoftheheadregionatdifferentamplitudes可以发现改变振幅无论是对整体气动阻尼还是头部气动阻尼来说变化都很小ꎬ这意味着气动阻尼的大小主要取决于气动力与结构振动的相位差ꎬ不依赖于振动幅度的大小ꎮ2.3.4㊀脉动压力对气动阻尼的影响为了模拟出脉动压力的影响ꎬ采用IDDES方法对火箭气动阻尼进行计算ꎬ计算来流Mach数为0.92ꎬ计算过程中的广义力与广义位移如图15所示ꎬ相较于图8可以看出广义力随时间变化曲线并不光滑ꎬ脉动压力的存在导致广义力由多个频率叠加而成ꎮ由于第2阶模态的频率为2.46Hzꎬ而由分离流㊁激波振荡等引起的脉动压力频率往往远大于此频率ꎬ因此这里选择3.5Hz为分界ꎬ将高于3.5Hz的部分视为由抖振脉动压力引起的广义力ꎬ低于3.5Hz的部分视为强迫振动引起的广义力ꎬ通过低通滤波把高于3.5Hz的广义力滤掉ꎬ可以获得由强迫振动引起的广义力与广义位移变化曲线ꎬ如图16所示ꎬ通过此广义力计算的气动阻尼为2.08ɢꎮ同样地ꎬ进行高通滤波将低于3.5Hz的广义力滤掉ꎬ可以获得由抖振脉动压力引起的气动阻尼为(2.94ˑ10-3)ɢꎬ由此得到脉动压力引起的气动阻尼变化为0.14%ꎬ可以忽略不计ꎮ同时使用RANS方法计算的气动阻尼为2.07ɢꎬ与IDDES的计算结果相比误差约为(2.94ˑ10-3+2.08-2.07)/2.07ʈ0.48%ꎬ这说明针对气动阻尼的模拟ꎬ抖振引起的脉动压力对气动阻尼的计算结果影响很小ꎬ起主要作用的还是广义力的变化ꎬ该变化由强迫振动引起的结构边界变化所导致ꎮ图15㊀基于IDDES的广义力与广义位移变化曲线Fig.15㊀VariationcuresofgeneralizedforceandgeneralizeddisplacementbasedonIDDES图16㊀滤波后的广义力与广义位移变化曲线Fig.16㊀Variationcuresofgeneralizedforceandgeneralizeddisplacementvariationcurveafterfiltering3㊀结论本文通过数值计算方法研究了火箭的气动阻尼特性ꎮ根据流动特征分析与理论推导ꎬ发现火箭过渡段几何外形的收缩导致该区域出现复杂的分离与激波结构ꎬ从而造成了气动力相对于结构振动83第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析相位的滞后ꎬ导致了该区域为气动负阻尼ꎬ即气动不稳定性的主要来源ꎮ在此机理的基础上ꎬ分析了前节点位置㊁振动振幅㊁脉动压力等因素对气动阻尼的影响规律ꎮ可以得出以下结论:1)助推增加了正阻尼区域的面积ꎬ从而相对于没有助推的构型起到了增加气动阻尼的作用ꎮ2)前节点位置的改变对过渡段气动阻尼影响很大ꎬ节点前后的振动方向相反ꎬ导致节点在过渡段前后的气动阻尼变化规律也截然相反ꎬ将过渡段与头部区域放在节点的同一侧有助于增加气动阻尼ꎮ3)在不改变流场结构的前提下ꎬ改变振动的振幅ꎬ气动力也会产生相应幅度的变化ꎬ因此结构振幅对气动阻尼的影响可忽略不计ꎮ4)高频部分的广义力对气动阻尼的贡献很小ꎬ即结构振动引起的广义力变化对气动阻尼起主要作用ꎬ而脉动压力对计算气动阻尼影响不大ꎬ可忽略不计ꎮ参考文献(References)[1]㊀CoeCF.Steadyandfluctuatingpressuresattransonicspeedsontwospace ̄vehiclepayloadshape[R].NASATMX ̄503ꎬ1961.[2]ColeSRJrꎬHenningTLꎬRaineyAG.NASAspaceve ̄hicledesigncriteria[R].NASASP ̄8001(REV)ꎬ1964. [3]EricssonLEꎬRedingJP.Analysisofflowseparationeffectsonthedynamicsofalargespacebooster[J].Jour ̄nalofSpacecraftandRocketsꎬ1965ꎬ2(4):481 ̄490. [4]RedingJPꎬEricssonLE.Effectofaeroelasticconsidera ̄tionsonseasat ̄Apayloadshrouddesign[J].JournalofSpacecraftandRocketsꎬ1981ꎬ18(3):241 ̄247. [5]程镇煌.宇航飞行器跨音速气动弹性问题探讨[J].上海航天ꎬ1997(6):16 ̄21.ChengZH.Discussionontransonicpneumaticelasticityofspacecraft[J].AerospaceShanghaiꎬ1997(6):16 ̄21(inChinese).[6]倪嘉敏.我国运载火箭气动设计回顾[C].近代空气动力学研讨会论文集ꎬ2005.NiJM.ReviewofaerodynamicdesignofChinaᶄslaunchvehicle[C].Proceedingsofthemodernaerodynamicssymposiumꎬ2005(inChinese).[7]臧涛成ꎬ胡焕性.大长细比弹箭弹性效应研究综述[J].弹道学报ꎬ1999ꎬ11(3):89 ̄93ꎬ96.ZangTCꎬHuHX.Areviewofgreatslendernessratioprojectileelasticeffectresearch[J].JournalofBallisticsꎬ1999ꎬ11(3):89 ̄93ꎬ96(inChinese). [8]吴志刚ꎬ杨超.细长体弹箭的气动弹性问题与研究方法[C].第九届(2005年)全国空气弹性学术交流会论文ꎬ2005.WuZGꎬYangC.Aeroelasticityproblemsandresearchmethodsofslenderbodymissiles[C].9thNationalSym ̄posiumonAeroelasticityꎬ2005(inChinese). [9]张贺ꎬ黄晓鹏.弹性细长旋转弹箭运动稳定性问题的研究进展[C].第九届(2005年)全国空气弹性学术交流会论文ꎬ2005.ZhangHꎬHuangXP.Researchprogressonthestabilityofelasticslenderrotatingprojectiles[C].9thNationalSymposiumonAeroelasticityꎬ2005(inChinese). [10]ColeHAJr.Dynamicresponseofhammerheadlaunchvehiclestotransonicbuffeting[R].NASATND ̄1982ꎬ1963.[11]BartelsREꎬWiesemanCDꎬMineckRE.ComputationalaeroelasticanalysisoftheAreslaunchvehicleduringas ̄cent[R].AIAA2010 ̄4374ꎬ2010.[12]AzevedoJLF.Aeroelasticanalysisoflaunchvehiclesintransonicflight[J].JournalofSpacecraftandRocketsꎬ1989ꎬ26(1):14 ̄23.[13]SinclairAꎬFlowersG.Low ̄orderaeroelasticmodeloflaunch ̄vehicledynamics[R].AIAA2010 ̄7725ꎬ2010. [14]DotsonKW.Transientcouplingoflaunchvehiclebendingresponseswithaerodynamicflowstatevariations[J].JournalofSpacecraftandRocketsꎬ2001ꎬ38(1):97 ̄104.[15]ColeSRꎬHenningTL.Buffetresponseofahammerheadlaunchvehiclewind ̄tunnelmodel[J].JournalofSpacecraftandRocketsꎬ1992ꎬ29(3):379 ̄385.[16]崔尔杰.流固耦合力学研究与应用进展[C].钱学森科学贡献暨学术思想研讨会论文集ꎬ2001.CuiEJ.Researchandapplicationprogressoffluid ̄struc ̄tureinteractionmechanics[C].ProceedingsofSeminarofQianXuesenScientificContributionsandAcademicThoughtsꎬ2001(inChinese).[17]冯明溪ꎬ王志安.火箭跨音速动导数和抖振实验[J].宇航学报ꎬ1987(1):55 ̄63.FengMXꎬWangZA.Experimentsoftransonicderiva ̄tivesandbuffetingofrocket[J].JournalofAstronauticsꎬ1987(1):55 ̄63(inChinese).[18]白葵ꎬ冯明溪.弹性模型实验技术[J].流体力学实验与测量ꎬ1999ꎬ13(1):38 ̄42.BaiKꎬFengMX.Aeroelasticmodelandthebuffetex ̄perimentaltechnique[J].ExperimentsandMeasurementsinFluidMechanicsꎬ1999ꎬ13(1):38 ̄42(inChinese). [19]JiCꎬRanJHꎬLiFꎬetal.Theaerodynamicdamping93气体物理2023年㊀第8卷testofelasticlaunchvehiclemodelintransonicflow[C].Proceedingsofthe64thInternationalAstronauticalCon ̄gressꎬ2013.[20]季辰ꎬ吴彦森ꎬ何岗ꎬ等.运载火箭气动阻尼风洞试验研究[C].第十二届全国空气弹性学术交流会论文集ꎬ2011.JiCꎬWuYSꎬHeGꎬetal.Experimentalstudyonaero ̄dynamicdampingwindtunneloflaunchvehicle[C].Pro ̄ceedingsofthe12thNationalAeroelasticityConferenceꎬ2011(inChinese).[21]刘子强ꎬ白葵ꎬ毛国良ꎬ等.锤头体弹性振动跨音速气动阻尼系数的确定[J].宇航学报ꎬ2002ꎬ23(6):1 ̄7.LiuZQꎬBaiKꎬMaoGLꎬetal.Thedeterminationofaerodynamicdampingonhammerheadlaunchvehiclesattransonicspeeds[J].JournalofAstronauticsꎬ2002ꎬ23(6):1 ̄7(inChinese).[22]冉景洪ꎬ刘子强ꎬ胡静ꎬ等.减阻杆气动阻尼研究[J].力学学报ꎬ2014ꎬ46(4):636 ̄641.RanJHꎬLiuZQꎬHuJꎬetal.Researchofaero ̄dampingforbluntwithspike[J].ChineseJournalofThe ̄oreticalandAppliedMechanicsꎬ2014ꎬ46(4):636 ̄641(inChinese).[23]朱剑ꎬ冉景洪ꎬ吴彦森ꎬ等.捆绑式运载火箭的气动阻尼数值计算方法[C].第十三届全国空气弹性学术交流会论文集.哈尔滨:中国力学学会ꎬ中国空气动力学会ꎬ2013.ZhuJꎬRanJHꎬWuYSꎬetal.Numericalcalculationmethodforaerodynamicdampingofbundlelaunchvehicles[C].Proceedingsofthe13thNationalAeroelasticityCon ̄ference.Harbin:ChineseSocietyofTheoreticalandAp ̄pliedMechanicsꎬChinaAeromechanicsSocietyꎬ2013(inChinese).[24]SpalartPR.Detached ̄eddysimulation[J].AnnualReviewofFluidMechanicsꎬ2009ꎬ41:181 ̄202. [25]GritskevichMSꎬGarbarukAVꎬSchützeJꎬetal.Devel ̄opmentofDDESandIDDESformulationsforthek ̄ωshearstresstransportmodel[J].FlowꎬTurbulenceandCombustionꎬ2012ꎬ88(3):431 ̄449.[26]季辰ꎬ吴彦森ꎬ侯英昱ꎬ等.捆绑式运载火箭跨声速气动阻尼特性试验研究[J].实验流体力学ꎬ2020ꎬ34(6):24 ̄31.JiCꎬWuYSꎬHouYYꎬetal.Experimentalstudyofaerodynamicdampingcharacteristicsofalaunchvehiclewithboostersintransonicflow[J].JournalofExperimentsinFluidMechanicsꎬ2020ꎬ34(6):24 ̄31(inChinese). [27]AllenCꎬRendallTCS.UnifiedapproachtoCFD ̄CSDinterpolationandmeshmotionusingradialbasisfunctions[R].AIAA2007 ̄3804ꎬ2007.04。
基于SIMPACK的车辆系统动力学性能分析
重载、曲线半径过小等问题ꎬ分析其动力学性能ꎮ 试验中根据脱轨系数、轮重减载率、轮轴横向
力等指标参数ꎬ对列车运行的平稳性及稳定性进行评估ꎮ
关键词:SIMPACKꎻ平稳性ꎻ安全性ꎻ直线运行ꎻ曲线通过
中图分类号:U270 文献标志码:B 文章编号:1671 ̄5276(2022)05 ̄0095 ̄03
根据建模后的离线分析结果ꎬ分别计算各运行速度下
平稳性等级
良好
的W y 、W z 、N MV 、Mα cy 、Mα cz 、Sα cy 及Sα cz ꎬ并根据所得数据进
表 3 客车运行平稳性指标与等级
2
TC2
E LNI
图 6 MP1 空车工况下平稳性指标平均值( 建模)
所有情况下的轮轨力以及其位置ꎬ因此还需进行更多相关
车轮滚动圆直径 / mm
耗ꎬ提高车辆的运行安全性和舒适性ꎬ基于车辆动力学理
钢轨 / ( kg / m)
种地铁线路缓和曲线线型进行研究ꎬ分析了各线型下的轮
车体质量 / kg
的研究 [6] ꎮ 周素霞等 [7] 为了减 轻 地 铁 运 行 中 的 轮 轨 磨
0.05
信号传输采用集流环装置ꎮ
0
5$
.1
DD4
B
6D2
,
AO
DFD(
0.4
图 7 轮对的蛇行运动
0.2
试验载荷工况包括:空车、重载ꎮ
列车在行使过程中会出现蛇行运动的现象ꎮ 蛇行运
0
动为非线性自激振动ꎮ 它是由蠕滑力和轮轨几何关系引
沙尘环境下高速列车气动特性分析
沙尘环境下高速列车气动特性分析作者:文恒于梦阁盛旭高殷硕来源:《青岛大学学报(工程技术版)》2020年第04期摘要:为研究沙尘环境下高速列车明线运行时的气动特性,基于剪切应力传输模型SST k w双方程湍流模型和拉格朗日离散相模型,与无沙环境下的高速列车气动特性进行比较,计算分析不同沙粒浓度、不同车速下的高速列车气动特性。
计算结果表明:沙尘环境下,当车速一定时,列车整车气动阻力、头车气动阻力、尾车气动阻力均随沙粒浓度增加而逐渐增大,且与沙粒浓度近似呈线性关系;对于气动升力,当车速一定时,头车气动升力绝对值随沙粒浓度的增加而增大,尾车气动升力随车速的增加而降低。
该研究成果可为高速列车在沙尘环境中的运行安全提供理论参考。
关键词:高速列车; 沙粒浓度; 气动特性; 气动阻力; 气动升力中图分类号: U271.91; U270.1+1 文献标识码: A近几年,随着科学技术的不断进步,国家加快铁路网建设,驱动区域经济地协调发展。
兰新高铁穿越我国西部高寒风沙区域,由于兰新高铁线路的特殊性,列车途经的百里风区、三十里风区是内陆大风天气频发的地区之一,时常会引发大风灾害性气象[1 3] ,强风地区大多缺少植被覆盖,地表裸露,而风速又往往远远大于起沙风速,因此地表的沙粒在大風的作用下撞击动车组,使得高速列车的气动特性明显变差[4] 。
沙尘环境属于多相流中的气固两相流问题,多采用欧拉欧拉模型和欧拉拉格朗日模型进行模拟计算。
欧拉欧拉模型主要描述两相的运动,用于高浓度离散相的问题;欧拉拉格朗日模型适用于离散相的体积分数在10%~12%以下的问题[5] ,而且关注离散相的运动轨迹。
C.Paz等人[6] 采用欧拉拉格朗日方法研究了高速列车在沙尘环境中的运行安全性;熊红兵等人[7] 研究了沙尘暴环境下高速列车运行时的气动特性;李田等人[8] 采用欧拉欧拉方法研究了不同沙尘暴环境下高速列车的动力学性能;倪守隆[9] 采用欧拉欧拉方法研究确定了高速列车在沙尘暴环境下运行的安全域;高琛光[10] 采用欧拉拉格朗日方法研究了不同风速、不同风向、头尾中间车设备舱的流动特性。
中国高速列车被动安全技术研究进展及思考
第53卷第5期2022年5月中南大学学报(自然科学版)Journal of Central South University (Science and Technology)V ol.53No.5May 2022中国高速列车被动安全技术研究进展及思考丁叁叁(中车青岛四方机车车辆股份有限公司,山东青岛,266111)摘要:从中国高速列车被动安全技术发展的切实需求与主要进程入手,重点阐述中国高速列车被动安全技术及工程化的总体思路,阐述列车碰撞响应、结构设计、装置研发、高效仿真优化、系统试验验证的具体过程,归纳列车耐撞设计所突破的能量管理、材料动态模拟、高效吸能方法、纠偏/防偏复合吸能导向、刚度匹配与力级分区、大系统集成、列车级评估等技术要点,建立高速列车碰撞性能仿真与试验的综合评估体系。
从车体材料—专用吸能部件—车端大部件—车辆—列车全链条维度,展示高速列车耐撞性技术的工程化路径。
最后,面向更高安全碰撞速度、更多样化材料结构、更广域运用条件、更完善评估体系与标准等,展望中国高速列车被动安全技术的发展方向。
关键词:高速列车;碰撞安全;吸能结构;变形抑制中图分类号:U270.3文献标志码:A文章编号:1672-7207(2022)05-1547-12Advances and prospect on passive safety technology for Chinesehigh-speed trainsDING Sansan(CRRC Qingdao Sifang Co.Ltd.,Qingdao 266111,China)Abstract:Starting from the main development processes of passive safety technology of high-speed trains in China,the overall idea of high-speed train technology in China was emphasized,the specific processes of train collision mechanism,process response,structure design,efficient simulation optimization and systematic experimental verification in detail were described,and the critical problems were summarized,such as energy management,dynamic simulation of materials,development of novel efficient energy-absorbing structure,correction of strong orientation,stiffness matching,integration of big system,train-level evaluation,correlation between simulation and test,etc.The comprehensive evaluation system of high-speed train collision performance收稿日期:2021−11−30;修回日期:2021−12−28基金项目(Foundation item):国家重点研发计划项目(2016YFB1200602);战略性国际科技创新合作重点专项资助项目(2018YFE0201300)(Project(2016YFB1200602)supported by the National Key Research and Development Program;Project (2018YFE0201300)supported by the Strategic International Scientific and Technological Innovation Cooperation)DOI:10.11817/j.issn.1672-7207.2022.05.001引用格式:丁叁叁.中国高速列车被动安全技术研究进展及思考[J].中南大学学报(自然科学版),2022,53(5):1547−1558.Citation:DING Sansan.Advances and prospect on passive safety technology for Chinese high-speed trains[J].Journal of Central South University(Science and Technology),2022,53(5):1547−1558.第53卷中南大学学报(自然科学版)simulation and test was gradually bined with the significant design points of material,component,large component and vehicle-train at every level,the engineering path of high-speed train crashworthiness technology was demonstrated.Finally,oriented to higher safe collision speed,more diversified materials and structures,wider application conditions,more complete evaluation system and standards,the development directions of passive safety technology for high-speed trains in China were prospected.Key words:high-speed train;collision safety;energy absorbing structure;deformation restrain列车运行安全性是铁路运输业非常重要的指标,列车发生碰撞事故造成的后果不堪设想,近20年来国内外列车碰撞事故时有发生,导致的人员伤亡及经济损失使列车被动安全防护技术的重要性日渐凸显,列车被动安全防护的核心是通过设计和开发专用的吸能装置、保持碰撞过程位形轨迹、提升整车耐撞性水平以最大限度地吸收列车碰撞时的强大冲击动能,降低碰撞作用力,保护司乘人员安全。
基于STAR_CCM_厢式半挂列车气动性能研究
第一作者:张克鹏,男,年生,工程师,从事车辆气动与热管理工作。
1985动量方程为平均速度分量,、为坐标分量,—u x (+—(—x —x —x —u x进行整车外流场分析时,其模拟风洞要减小其阻塞效应,研究表明,只有低于的风洞试验结果,阻塞干扰产生的误差才不需要修正。
根据这一理论,计算域上部留倍车高,左右各倍车宽。
从车辆空气动力学研究经验来看,在汽车风洞试验过程中,汽车尾部有一个比较大的湍流区域,这个区域的流动非常紊乱,数值计算过程中,这个区域如果模拟得不好对结果会产生重大的影响。
所以车辆尾部一般设定倍车长的区域,鉴于厢式半挂列车总车长较长,此处选择倍车长区域,保证汽车湍流区域充分发1%6586湍流动能耗散率方程:式中,为流体密度,为因平均速度梯度引起的湍流动能;e ?=4.0G A —t (—x —x ——(—x —x —(——u x —u x 图1厢式半挂列车面网格模型图2厢式半挂列车体网格模型——u x(a)A型(c)A型+B型+C型图7厢式半挂列车不同导流装置(b)A型+B型(d)A型+B型+C型+D型原型车CFD计算结果分析该厢式半挂列车原型车模型中驾驶室与车厢距离为,没有驾驶室导流罩等辅助装置。
滞止,形成一个正压区。
1图3厢式半挂列车体网格模型(中心对称面处)图4边界条件示意图速度进口固壁无滑移壁面边界滑移壁面边界导流装置对气动阻力特性的改进通过上面分析,利用比较成熟的种导流装置,进行整车气动阻力特性的改进。
在驾驶室顶部加装顶导流罩,定义为型。
侧面加装侧导流板,定义为型。
在半挂列车车厢前端增加导流装置,定义为型;在车厢后部增加导流装置,定义为型。
厢4a.A B c.C d.D 图5原型车表面压力等值云图图6原型车速度等值线云图(中心对称面处)边界条件的设定厢式半挂列车外流场的数值模拟过程中,前端来流方向的端面为入口边界,设定入口边界为速度入口,流速,湍流强度为;出口边界为压力出口,出口压力为相对于大气v p =85km/h 1%=0Pa(驾驶室与车厢距离对气动性能影响在最佳导流装置组合基础上,研究厢式半挂列车驾驶室与车厢间距变化对气动阻力系数的影响。
基于PowerFlow的重型车外气动性能分析及优化_张克鹏
FOCU2S01技5术年聚2焦月
设计·创新
基于 P ow erF low 的重型车
外气动性能分析及优化 *
张克鹏 (陕西重型汽车有限公司)
摘要:在某重型车驾驶室开发过程中,需考虑前扰流板对风阻系数及车门把手位置除尘效果的影响。利用基于 LBM 方法 的 CFD 软件 PowerFlow,对该重型车进行外流场计算;基于 CFD 分析结果,将前扰流板翅片角度下倾 15°。改进后的前扰 流板设计使整车风阻系数较原来降低 3.13%,除尘效果改善显著。改进重型车前扰流板的设计可以有效提升整车外气动 性能,该方法为前扰流板气动性能设计提供了理论依据。 关键词:前扰流板;PowerFlow;LBM 方法;气动性能;除尘;改进
- 35 -
技术聚焦 FOCUS
0 5 10 15 20 25 30 35 速度(/ m/s)
图 5 重型车车身中心对称面速度矢量云图
图 6 示出该重型牵引车中心对称切面的总压云 图。总压为静压与动压之和,总压在一定程度上就代表 气流动能。其中总压计算公式为:
Pt=Ps+1/2ρv2 式中:Pt———总压,Pa;
产生的误差不需要修正。根据这一理论,计算域上部留 车行驶过程中,前方气流首先遇到驾驶室前部,在此处
6 倍车高,左右各 5 倍车宽,汽车尾部一般设定 8 倍车 气流滞止,阻力发展曲线中风阻系数直线上升;在前面
长的区域。鉴于重型车总车长较长,这里选择 6 倍车长 罩、保险杠及挡风玻璃的大部分区域形成较大的正压,
后的数据作为输入条件。散热器、冷凝器及中冷器也采
取同样的方法,将修正后各模块单体性能试验数据作
为输入。
2 CFD 计算结果分析
2.1 气动阻力分析 气动阻力与气流流经车身的过程有关。气流在流
基于混合遗传算法的高速公路桥梁风屏障参数优化
基于混合遗传算法的高速公路桥梁风屏障参数优化作者:罗颖廖楚峰韩艳蔡春声何旭辉李凯来源:《湖南大学学报·自然科学版》2022年第01期摘要:为了优化高速公路桥梁风屏障参数,研究了风屏障参数对于车桥系统气动特性的影响.通过风洞试验考虑不同高度和透风率的风屏障,分别获取桥梁和车辆的气动力系数,进而得到桥梁在静风稳定性检验风速下的侧向位移和车辆在设计车速行驶下的失稳临界风速.采用多目标遗传算法(NSGA-II),以桥梁侧向位移与车辆临界风速为优化目标,将风屏障高度和透风率作为变量,得到相应的 Pareto 最优解集.利用数据包络法(DEA)对 Pareto 解集中个体的相对效率值进行评估,最终得到最优风屏障参数.结果表明:透风率为30%、高度为3.2 m 的风屏障对于桥梁和车辆的综合抗风效果最佳.关键词:风屏障;多目标遗传算法;Pareto 解集;数据包络法(DEA)中图分类号:U447 文献标志码:AParameter Optimization of Wind Barrier for Highway Bridge Based on Hybrid Genetic AlgorithmLUO Ying1,LIAO Chufeng1,HAN Yan1†,CAI Chunsheng1,2,HE Xuhui3,LI Kai1(1. School of Civil Engineering,Changsha University of Science & Technology,Changsha 410114,China;2. Department of Civil and Environmental Engineering,Louisiana State University,Baton Rouge 70803,USA;3. School of Civil Engineering,Central South University,Changsha 410075,China)Abstract:The impacts of wind barrier parameters on the aerodynamic characteristics of the vehicle-bridge sys ⁃ temare investigated to optimize the wind barrier parameters of the highway bridge. Based on wind tunnel tests con⁃ sidering the wind barriers with different heights and ventilation rates, the aerodynamic coefficients for the bridge and the vehicle are obtained,respectively. Afterward, the lateral displacement of the bridge under testing wind speed of aerostatic stability and the critical wind speed of aerostatic instability for the vehicle running with the design speed are calculated, respectively. A multi-objective genetic algorithm (NSGA-II) is applied to obtain the related Pareto op⁃ timal solution set. In the algorithm, the lateral displacement of the bridge and the critical wind speed for the vehicleare considered as optimal objectives. Meanwhile,the height and the ventilation rate of the wind barrier are regarded as the variables. Taking advantage of the Data Envelopment Analysis (DEA) to evaluate the relative efficiency of the individual in the Pareto solution set, the optimal parameter of the wind barrier can be determined finally. The results show that the wind barrier with the ventilation rate of 30% and the height of 3.2 m has the best performance in consid ⁃ eration of the wind-resistant effect for both the bridge and the vehicle.Key words:wind barrier;multi-objective genetic algorithm;Pareto solution set;Data Envelopment Analysis ( DEA)近些年来,由强风作用所导致的高速公路交通安全事故频发,不仅危及人身安全,也给经济财产造成了巨大损失.为提高车辆在强风作用下的行驶安全性,目前最常用的措施是在桥梁主梁两侧设置风屏障,以改善车辆行驶风环境,减小风荷载对车辆的影响[1].目前国内外已有多位学者通过风洞试验、数值模拟与现场实测等研究方法,结合空气动力学,从车辆的流场分布、气动力系数及动力响应等方面对风屏障的防风效果开展了研究分析.葛盛昌等[2]在强风地区铁路干线两侧设置挡风墙,通过现场实测发现设置挡风墙能大大降低其内侧的风速. Coleman 等[3]通过风洞试验测试了有无风屏障情况下的车辆气动特性,结果表明风屏障能有效减小车辆气动力系数. Chen等[4]基于风洞试验和 CFD数值模拟,分别得到桥梁及车辆的气动力系数,计算并对比了有无风屏障时桥梁及车辆的动力响应,分析了风屏障对行车安全性的影响.结合风洞试验和数值模拟,何玮等[5-6]探讨了风屏障参数对车桥系统的影响,发现风屏障参数的变化对列车阻力系数的影响较为明显,且列车位于迎风侧时其气动特性对风屏障参数的变化更为敏感.以上研究仅限于风屏障参数对车桥系统气动性能的影响,缺乏对风屏障参数优化的研究.针对该问题,向活跃等[7]通过风洞试验测试不同高度风屏障下车辆的气动力系数,以车辆的风荷载突变量作为评价指标,基于数据包络法(DEA)评估了风屏障的防风效果,但只能在已有的方案中进行选择,容易陷入局部最优,且只是基于对列车影响的风屏障高度优化.为了开展对车桥系统影响的风屏障参数优化,苏洋等[8]通过风洞试验与数值模拟计算了不同高度风屏障下车辆及桥梁的气动力系数,以车辆倾覆力矩与桥梁阻力系数作为优化目标,为了避免DEA方法易陷入局部最优及多目标优化不能直接判断个体优劣的缺点,引入 NSGA-II&DEA 混合算法探讨了风屏障高度的优化,但只是风屏障的单一参数优化.向活跃等[9]基于改进的网格搜索法,以车辆气动特性为目标函数,建立了风屏障防风效果的优化模型,对风屏障的高度和透风率进行优化.然而,上述研究只将桥梁与车辆的气动力系数作为研究目标,没有直接反映风屏障参数对车桥系统响应的影响,且主要针对铁路桥梁,较少涉及到公路桥梁.基于某三塔双索面叠合梁斜拉桥,本文通过风洞试验分别测量桥梁和车辆在不同高度和透风率风屏障下的气动力系数,进而分析得到主梁跨中在静风稳定检验风速下的侧向位移和车辆在设计车速下的失稳临界风速;然后利用多目标遗传算法—— NSGA-II,对风屏障参数进行优化,计算得到在不同高度及透风率下的 Pareto解集;最后,根据数据包络法(DEA),计算每个解集中个体的相对效率值,确定风屏障的最优参数.1风洞试验1.1工程背景本文基于某三塔双索面叠合梁斜拉桥开展研究,该桥的桥跨布置为(249.5+550+550+249.5) m,具体如图1所示.主梁采用钢混叠合梁,宽度为30.2 m,边主梁高2.92 m,主梁截面如图2所示.大桥跨越宽达1600 m、深约300 m 的峡谷,桥面处高300 m,风速较大,需要在两侧设置风屏障以保障车辆的行驶安全性.1.2试验概况为了研究风屏障参数对车辆及桥梁气动力的影响,分别开展了主梁测力试验与车辆测压试验.主梁节段模型测力试验在长沙理工大学风洞实验室高速段进行,其尺寸为21 m×4 m×3 m(长×宽×高).主梁节段模型缩尺比取1∶40,如图3(a)所示.车辆测压试验是在中南大学高速铁路建造技术国家工程实验室的高速铁路风洞试验系统进行,其尺寸为15 m×3 m×3 m(长×宽×高).文献[10]表明,相比小轿车,厢式货车更容易发生失稳,因此车辆测压试验中以厢式货车为研究对象,具体尺寸见表1.主梁和车辆模型缩尺比取1∶32,如图3(b)(c)所示.基于以往研究[5-6]考慮最不利情况,仅测量车辆位于迎风侧且垂直于来流时的气动力系数.风屏障透风率依照等效透风率原则分别取10%、20%和30%,见图3(d),高度分别取1.2 m、2.0 m、3.0 m 及4.0 m.试验模型具体如图3所示.1.3试验结果分析主梁的三分力系数定义参考文献[11],为简便起见,仅计算风轴坐标系下的主梁气动力系数,公式如下:式中:α为风攻角;CD (α)、CL (α)与 CM (α)分别为风轴坐标系下不同风攻角对应的主梁阻力系数、升力系数与扭矩系数;FD (α)、FL (α)与 M(α)分别为不同风攻角下主梁的阻力、升力与扭矩;U∞为试验来流平均风速,取10 m/s;ρ为空气质量密度,取1.225 kg/m3;H、L、B 分别为节段模型高度、长度和宽度.参考相关文献[5-6],主梁阻力系数 CD 对于风屏障的参数变化较为敏感.由于篇幅所限,以下仅给出阻力系数 CD 在不同风屏障参数下随风攻角的变化情况,具体见图4.从图4可以看出,在不同风屏障参数下,主梁阻力系数 CD 随风攻角的变化趋势大体一致,均在0°风攻角附近达到最小值,并随着风攻角绝对值的增大而增加.这是由于0°风攻角下的主梁断面整体阻风面积最小,随着风攻角的变化,主梁阻风面积逐渐变大.当风屏障高度较小时,CD 对风屏障透风率的变化较为敏感,但达到一定高度后透风率的影响作用减小,且均在0°风攻角附近透风率的影响较大.这是因为0°风攻角附近风屏障透风率的变化对主梁整体阻风面积的影响较大,而在大攻角下风屏障透风率的变化引起阻风面积的变化相对主梁自身阻风面积已经较小.针对这类大跨度斜拉桥,参考文献[12]对此桥进行了三维非线性静风响应分析,得到桥梁的静风响应.通过对比-5°~5°风攻角下的跨中侧向位移,-5°风攻角下的侧向位移最大,因此后续选取该风攻角下的情况予以分析.-5°风攻角下主梁跨中侧向位移随风屏障高度及透风率变化趋势如图5所示.由图5看出,主梁的侧向位移随着风屏障透风率的减少而增加;但随着风屏障高度的增加,透风率的影响逐渐减弱,这主要是由于风屏障高度达到一定值后,风屏障透风率的变化对于桥梁气动力的影响减小,这与主梁三分力系数的变化大体一致.根据参考文献[11]计算得到此桥的静风稳定检验风速为66 m/s,将该风速下的主梁侧向位移作为反映风屏障参数对于桥梁影响的指标,位移越小则桥梁越稳定,结果见表2.从表2可看出,主梁跨中侧向位移的变化趋势与主梁阻力系数 CD 的变化趋势是一致的,随着风屏障高度的增加、透风率的减少,侧向位移增大.但是4 m 高度下并不完全符合这个趋势,最大值出现在20%透风率的情况,这可能是因为此时主梁的阻力系数相差不大,升力系数和扭矩系数对侧向位移的影响变得突出.车辆气动力系数的处理参考文献[13],定义如下:式中:A 为车辆迎风侧面积;hv 为车辆质心距车辆底面距离;CS、CL、CD、Cp、CY 和 CR 分别为车辆的侧向力系数、升力系数、阻力系数、俯仰力矩系数、偏转力矩系数和侧倾力矩系数;FS、FL、FD、MP、MY 和 MR 分别为作用于车辆质心的侧力、升力、阻力、俯仰力矩、偏转力矩和侧倾力矩.以上研究仅限于风屏障参数对车桥系统气动性能的影响,缺乏对风屏障参数优化的研究.针对该问题,向活跃等[7]通过风洞试验测试不同高度风屏障下车辆的气动力系数,以车辆的风荷载突变量作为评价指标,基于数据包络法(DEA)评估了风屏障的防风效果,但只能在已有的方案中进行选择,容易陷入局部最优,且只是基于对列车影响的风屏障高度优化.为了开展对车桥系统影响的风屏障参数优化,苏洋等[8]通过风洞试验与数值模拟计算了不同高度风屏障下车辆及桥梁的气动力系数,以车辆倾覆力矩与桥梁阻力系数作为优化目标,为了避免DEA方法易陷入局部最优及多目标优化不能直接判断个体优劣的缺点,引入 NSGA-II&DEA混合算法探讨了风屏障高度的优化,但只是风屏障的单一参数优化.向活跃等[9]基于改进的网格搜索法,以车辆气动特性为目标函数,建立了风屏障防风效果的优化模型,对风屏障的高度和透风率进行优化.然而,上述研究只将桥梁与车辆的气动力系数作为研究目标,没有直接反映风屏障参数对车桥系统响应的影响,且主要针对铁路桥梁,较少涉及到公路桥梁.基于某三塔双索面叠合梁斜拉桥,本文通过风洞试验分别测量桥梁和车辆在不同高度和透风率风屏障下的气动力系数,进而分析得到主梁跨中在静风稳定检验风速下的侧向位移和车辆在设计车速下的失稳临界风速;然后利用多目标遗传算法—— NSGA-II,对风屏障参数进行优化,计算得到在不同高度及透风率下的 Pareto解集;最后,根据数据包络法(DEA),计算每个解集中个体的相对效率值,确定风屏障的最优参数.1风洞试验1.1工程背景本文基于某三塔双索面叠合梁斜拉桥开展研究,该桥的桥跨布置为(249.5+550+550+249.5) m,具体如图1所示.主梁采用钢混叠合梁,宽度为30.2 m,边主梁高2.92 m,主梁截面如图2所示.大桥跨越宽达1600 m、深约300 m 的峡谷,桥面处高300 m,风速较大,需要在两侧设置风屏障以保障车辆的行驶安全性.1.2试验概况为了研究风屏障参数对车辆及桥梁气动力的影响,分别开展了主梁测力试验与车辆测压试验.主梁节段模型测力试验在长沙理工大学风洞实验室高速段进行,其尺寸为21 m×4 m×3 m (长×宽×高).主梁节段模型缩尺比取1∶40,如图3(a)所示.车辆测压试验是在中南大学高速铁路建造技术国家工程实验室的高速铁路风洞试验系统进行,其尺寸为15 m×3 m×3 m(长×宽×高).文献[10]表明,相比小轿车,厢式货车更容易发生失稳,因此车辆测压试验中以厢式货车为研究对象,具体尺寸见表1.主梁和车辆模型缩尺比取1∶32,如图3(b)(c)所示.基于以往研究[5-6]考虑最不利情况,仅测量车辆位于迎风侧且垂直于来流时的气动力系数.风屏障透风率依照等效透风率原则分别取10%、20%和30%,见图3(d),高度分别取1.2 m、2.0 m、3.0 m 及4.0 m.试验模型具体如图3所示.1.3试验结果分析主梁的三分力系数定义参考文献[11],为简便起见,仅计算风轴坐标系下的主梁气动力系数,公式如下:式中:α为风攻角;CD (α)、CL (α)与 CM (α)分别为风轴坐标系下不同风攻角对应的主梁阻力系数、升力系数与扭矩系数;FD (α)、FL (α)与 M(α)分别为不同风攻角下主梁的阻力、升力与扭矩;U∞为试验来流平均风速,取10 m/s;ρ为空气质量密度,取1.225 kg/m3;H、L、B 分别为节段模型高度、长度和宽度.参考相关文献[5-6],主梁阻力系数 CD 对于风屏障的参数变化较为敏感.由于篇幅所限,以下仅给出阻力系数 CD 在不同风屏障参数下随风攻角的变化情况,具体见图4.从图4可以看出,在不同风屏障参数下,主梁阻力系数 CD 随风攻角的变化趋势大体一致,均在0°风攻角附近达到最小值,并随着风攻角绝对值的增大而增加.这是由于0°风攻角下的主梁断面整体阻风面积最小,随着风攻角的变化,主梁阻风面积逐渐变大.当风屏障高度较小时,CD 对风屏障透风率的变化较为敏感,但达到一定高度后透风率的影响作用减小,且均在0°风攻角附近透风率的影响较大.這是因为0°风攻角附近风屏障透风率的变化对主梁整体阻风面积的影响较大,而在大攻角下风屏障透风率的变化引起阻风面积的变化相对主梁自身阻风面积已经较小.针对这类大跨度斜拉桥,参考文献[12]对此桥进行了三维非线性静风响应分析,得到桥梁的静风响应.通过对比-5°~5°风攻角下的跨中侧向位移,-5°风攻角下的侧向位移最大,因此后续选取该风攻角下的情况予以分析.-5°风攻角下主梁跨中侧向位移随风屏障高度及透风率变化趋势如图5所示.由图5看出,主梁的侧向位移随着风屏障透风率的减少而增加;但随着风屏障高度的增加,透风率的影响逐渐减弱,这主要是由于风屏障高度达到一定值后,风屏障透风率的变化对于桥梁气动力的影响减小,这与主梁三分力系数的变化大体一致.根据参考文献[11]计算得到此桥的静风稳定检验风速为66 m/s,将该风速下的主梁侧向位移作为反映风屏障参数对于桥梁影响的指标,位移越小则桥梁越稳定,结果见表2.从表2可看出,主梁跨中侧向位移的变化趋势与主梁阻力系数 CD 的变化趋势是一致的,随着风屏障高度的增加、透风率的减少,侧向位移增大.但是4 m 高度下并不完全符合这个趋势,最大值出现在20%透风率的情况,这可能是因为此时主梁的阻力系数相差不大,升力系数和扭矩系数对侧向位移的影响变得突出.车辆气动力系数的处理参考文献[13],定义如下:式中:A 为车辆迎风侧面积;hv 为车辆质心距车辆底面距离;CS、CL、CD、Cp、CY 和 CR 分别为车辆的侧向力系数、升力系数、阻力系数、俯仰力矩系数、偏转力矩系数和侧倾力矩系数;FS、FL、FD、MP、MY 和 MR 分别为作用于车辆质心的侧力、升力、阻力、俯仰力矩、偏转力矩和侧倾力矩.以上研究仅限于风屏障参数对车桥系统气动性能的影响,缺乏对风屏障参数优化的研究.针对该问题,向活跃等[7]通过风洞试验测试不同高度风屏障下车辆的气动力系数,以车辆的风荷载突变量作为评价指标,基于数据包络法(DEA)评估了风屏障的防风效果,但只能在已有的方案中进行选择,容易陷入局部最优,且只是基于对列车影响的风屏障高度优化.为了开展对车桥系统影响的风屏障参数优化,苏洋等[8]通过风洞试验与数值模拟计算了不同高度风屏障下车辆及桥梁的气动力系数,以车辆倾覆力矩与桥梁阻力系数作为优化目标,为了避免DEA方法易陷入局部最优及多目标优化不能直接判断个体优劣的缺点,引入 NSGA-II&DEA混合算法探讨了风屏障高度的优化,但只是风屏障的单一参数优化.向活跃等[9]基于改进的网格搜索法,以车辆气动特性为目标函数,建立了风屏障防风效果的优化模型,对风屏障的高度和透风率进行优化.然而,上述研究只将桥梁与车辆的气动力系数作为研究目标,没有直接反映风屏障参数对车桥系统响应的影响,且主要针对铁路桥梁,较少涉及到公路桥梁.基于某三塔双索面叠合梁斜拉桥,本文通过风洞试验分别测量桥梁和车辆在不同高度和透风率风屏障下的气动力系数,进而分析得到主梁跨中在静风稳定检验风速下的侧向位移和车辆在设计车速下的失稳临界风速;然后利用多目标遗传算法—— NSGA-II,对风屏障参数进行优化,计算得到在不同高度及透风率下的 Pareto解集;最后,根据数据包络法(DEA),计算每个解集中个体的相对效率值,确定风屏障的最优参数.1风洞试验1.1工程背景本文基于某三塔双索面叠合梁斜拉桥开展研究,该桥的桥跨布置为(249.5+550+550+249.5) m,具体如图1所示.主梁采用钢混叠合梁,宽度为30.2 m,边主梁高2.92 m,主梁截面如图2所示.大橋跨越宽达1600 m、深约300 m 的峡谷,桥面处高300 m,风速较大,需要在两侧设置风屏障以保障车辆的行驶安全性.1.2试验概况为了研究风屏障参数对车辆及桥梁气动力的影响,分别开展了主梁测力试验与车辆测压试验.主梁节段模型测力试验在长沙理工大学风洞实验室高速段进行,其尺寸为21 m×4 m×3 m (长×宽×高).主梁节段模型缩尺比取1∶40,如图3(a)所示.车辆测压试验是在中南大学高速铁路建造技术国家工程实验室的高速铁路风洞试验系统进行,其尺寸为15 m×3 m×3 m(长×宽×高).文献[10]表明,相比小轿车,厢式货车更容易发生失稳,因此车辆测压试验中以厢式货车为研究对象,具体尺寸见表1.主梁和车辆模型缩尺比取1∶32,如图3(b)(c)所示.基于以往研究[5-6]考虑最不利情况,仅测量车辆位于迎风侧且垂直于来流时的气动力系数.风屏障透风率依照等效透风率原则分别取10%、20%和30%,见图3(d),高度分别取1.2 m、2.0 m、3.0 m 及4.0 m.试验模型具体如图3所示.1.3试验结果分析主梁的三分力系数定义参考文献[11],为简便起见,仅计算风轴坐标系下的主梁气动力系数,公式如下:式中:α为风攻角;CD (α)、CL (α)与 CM (α)分别为风轴坐标系下不同风攻角对应的主梁阻力系数、升力系数与扭矩系数;FD (α)、FL (α)与 M(α)分别为不同风攻角下主梁的阻力、升力与扭矩;U∞为试验来流平均风速,取10 m/s;ρ为空气质量密度,取1.225 kg/m3;H、L、B 分别为节段模型高度、长度和宽度.参考相关文献[5-6],主梁阻力系数 CD 对于风屏障的参数变化较为敏感.由于篇幅所限,以下仅给出阻力系数 CD 在不同风屏障参数下随风攻角的变化情况,具体见图4.从图4可以看出,在不同风屏障参数下,主梁阻力系数 CD 随风攻角的变化趋势大体一致,均在0°风攻角附近达到最小值,并随着风攻角绝对值的增大而增加.这是由于0°风攻角下的主梁断面整体阻风面积最小,随着风攻角的变化,主梁阻风面积逐渐变大.当风屏障高度较小时,CD 对风屏障透风率的变化较为敏感,但达到一定高度后透风率的影响作用减小,且均在0°风攻角附近透风率的影响较大.这是因为0°风攻角附近风屏障透风率的变化对主梁整体阻风面积的影响较大,而在大攻角下风屏障透风率的变化引起阻风面积的变化相对主梁自身阻风面积已经较小.针对这类大跨度斜拉桥,参考文献[12]对此桥进行了三维非线性静风响应分析,得到桥梁的静风响应.通过对比-5°~5°风攻角下的跨中侧向位移,-5°风攻角下的侧向位移最大,因此后续选取该风攻角下的情况予以分析.-5°风攻角下主梁跨中侧向位移随风屏障高度及透风率变化趋势如图5所示.由图5看出,主梁的侧向位移随着风屏障透风率的减少而增加;但随着风屏障高度的增加,透风率的影响逐渐减弱,这主要是由于风屏障高度达到一定值后,风屏障透风率的变化对于桥梁气动力的影响减小,这与主梁三分力系数的变化大体一致.根据参考文献[11]计算得到此桥的静风稳定检验风速为66 m/s,将该风速下的主梁侧向位移作为反映风屏障参数对于桥梁影响的指标,位移越小则桥梁越稳定,结果见表2.从表2可看出,主梁跨中侧向位移的变化趋势与主梁阻力系数 CD 的变化趋势是一致的,随着风屏障高度的增加、透风率的减少,侧向位移增大.但是4 m 高度下并不完全符合这个趋势,最大值出现在20%透风率的情况,这可能是因为此时主梁的阻力系数相差不大,升力系数和扭矩系数对侧向位移的影响变得突出.车辆气动力系数的处理参考文献[13],定义如下:式中:A 为车辆迎风侧面积;hv 为车辆质心距车辆底面距离;CS、CL、CD、Cp、CY 和 CR 分别为车辆的侧向力系数、升力系数、阻力系数、俯仰力矩系数、偏转力矩系数和侧倾力矩系数;FS、FL、FD、MP、MY 和 MR 分别为作用于车辆质心的侧力、升力、阻力、俯仰力矩、偏转力矩和侧倾力矩.以上研究仅限于风屏障参数对车桥系统气动性能的影响,缺乏对风屏障参数优化的研究.针对该问题,向活跃等[7]通过风洞试验测试不同高度风屏障下车辆的气动力系数,以车辆的风荷载突变量作为评价指标,基于数据包络法(DEA)评估了风屏障的防风效果,但只能在已有的方案中进行选择,容易陷入局部最优,且只是基于对列车影响的风屏障高度优化.为了开展对车桥系统影响的风屏障参数优化,苏洋等[8]通过风洞试验与数值模拟计算了不同高度风屏障下车辆及桥梁的气动力系数,以车辆倾覆力矩与桥梁阻力系数作为优化目标,为了避免DEA方法易陷入局部最优及多目标优化不能直接判断个体优劣的缺点,引入 NSGA-II&DEA混合算法探讨了风屏障高度的优化,但只是风屏障的单一参数优化.向活跃等[9]基于改进的网格搜索法,以车辆气动特性为目标函数,建立了风屏障防风效果的优化模型,对风屏障的高度和透风率进行优化.然而,上述研究只将桥梁与车辆的气动力系数作为研究目标,没有直接反映风屏障参数对车桥系统响应的影响,且主要针对铁路桥梁,较少涉及到公路桥梁.基于某三塔双索面叠合梁斜拉桥,本文通过风洞试验分别测量桥梁和车辆在不同高度和透风率风屏障下的气动力系数,进而分析得到主梁跨中在静风稳定检验风速下的侧向位移和车辆在设计车速下的失稳临界风速;然后利用多目标遗传算法—— NSGA-II,对风屏障参数进行优。
基于虚拟试验场的牵引车动态载荷研究
2024年第1期27doi:10.3969/j.issn.1005-2550.2024.01.005 收稿日期:2023-10-27基于虚拟试验场的牵引车动态载荷研究王庆华1,王丽荣2,陈小华2,李蒙然1,黄刚1(1.国家汽车质量检验检测中心(襄阳),襄阳441004;2. 北京福田戴姆勒汽车有限公司,北京 101400)摘 要:基于Adams软件的虚拟试验场动态载荷分解技术在乘用车耐久性能开发领域广泛应用。
对于重卡车型,由于车辆模型复杂、参数有限且测试难度大,虚拟试验场技术的应用推广受到限制。
搭建某牵引车整车多体动力学模型及虚拟试验场仿真环境,同时采集试验场工况下的实车载荷谱数据并与虚拟试验场动力学仿真分析提取的动态载荷进行对比。
使用相对伪损伤比值、频谱分析等评估比利时、扭曲路、搓板路等典型路面工况下仿真与实测载荷谱数据的差异。
结果表明:基于虚拟试验场的动态载荷提取技术可应用于牵引车车型且可实现较高的精度,是一种获取试验场耐久工况载荷谱的有效方法。
关键词:虚拟试验场;载荷分解;路面模型;牵引车中图分类号:U467 文献标识码:A 文章编号:1005-2550(2024)01-0027-07Research on Dynamic Load of Tractor Based on VPGWANG Qing-hua1, WANG Li-rong2, CHEN Xiao-hua2, LI Meng-ran1, HUANG Gang1(1.National Automobile Quality Inspection and T est Center (Xiangyang), Xiangyang 441004,China; 2. Beijing Foton Daimler Automobile Co., Ltd, Beijing 101400, China)Abstract: The dynamic load decomposition technology of VPG based on Adams is widely applied in the field of passenger car durability performance development. For heavytruck, the application and promotion of VPG are limited due to the complexity of vehiclemodels, limited parameters, and high RLDA testing difficulty. The complete vehicle multi-body dynamics model of a tractor and virtual proving ground simulation environment arebuilt based on Adams. The real vehicle load data acquisition of the proving ground eventswas carried out and compared with the dynamic loads extracted from dynamic simulationanalysis of the virtual proving ground to verify the model accuracy and load accuracy.Relative pseudo damage ratio, RMS value ratio, and spectrum analysis were used to evaluatethe differences between simulated and measured load data under typical road conditionssuch as Belgium, twisted roads, and washboard roads. It is proved that The dynamic loadextraction technology based on virtual proving ground can be applied to tractor models andachieve high accuracy, which is an effective method for obtaining the load data of provingground durability events.Key Words: Virtual Proving Ground; Load Extraction; Road Model; Tractor随着高精度路面扫描和轮胎力学模型建模等技术快速发展,基于虚拟试验场(V i r t u a l Proving Ground)的动态载荷提取技术在车型开发早期阶段即可开展,可有效缩短开发周期和试验成本[1-4]。
高速列车在隧道内和明线上交会时气动性能对比分析
坐 舒适 性 等 受 影 响 较 大 的 气 动 性 能 没 有 进 行 详 细 研 究 。文 献 [ ] 用一 维 可 压 缩 不 等 熵 非 定 常 流体 流 动 6采
简 化 的列 车 模 型 忽 略 了转 向架 、 电弓 、 轮 、 受 车 车
钩等 细 部物 件 , 采用 头车 、 中间 车 和尾 车 编 组 , 7 . 长 64
模拟 空气 黏性 对列 车气 动 力 的作用 。为真 实模 拟 高速
列车 在 隧道 内和 明线 上 的交 会 , 文 建 立 了简 化 的某 本
理论 研 究 了列车 在 隧道 中交 会 时的 瞬变 压力 波 。由于 隧道 中流场计 算 的复 杂 性 , 内外 学 者 在 计 算 时所 采 国 取 的简化 措施 ( 一维 模 型) 法模 拟 列 车车 头形状 对 如 无 隧道 内压 缩波 以及 会 车 压 力 波 的影 响 , 黏 流体 无 法 无
试 验 研 究
文 章 编 号 :0 2 7 0 ( 0 0 1 — 0 10 10 —6 2 2 1 )20 0 —5
铁道车辆 第4 卷第1 期 21 年1 月 8 2 00 2
. ..
同 列 车 在 隧 道 内和 明 线 上 交 会 时 速 气 动 性 能 对 比分 析
‘
冯 志 鹏 , 继 业 , 卫 华 张 张
用 的 e两方 程模 型 。本文 通 过 求 解 三维 瞬 态 可 压缩 的雷诺 时均 N S方 程 和 ke两 方程 湍 流模 型获得 列 车 - -
重卡驾驶室顶盖冲压仿真分析
10.16638/ki.1671-7988.2019.14.047重卡驾驶室顶盖冲压仿真分析张凯鹏(陕西重型汽车有限公司,陕西西安710200)摘要:重卡驾驶室顶盖结构复杂且尺寸较大,成形过程中会出现起皱、开裂、回弹等工艺问题,通过对板料成形过程理论和冲压缺陷预测的理论研究,应用AutoForm软件,采用动力显式算法对零件的拉延工序进行数值模拟,同时调整工艺参数,解决冲压件工艺问题。
关键词:顶盖外板;数值模拟;冲压成形;回弹中图分类号:U467 文献标识码:A 文章编号:1671-7988(2019)14-145-04Rebound of Stamping Forming of Roof Cap of Heavy TruckZhang Kaipeng(ShaanXi Heavy-duty Automobile Co., LtD, Shaanxi Xi'an 710200)Abstract: The Roof Cap of Heavy truck has complex structure and large size, so corrugation, crack and rebound may occur in theforming process. This paper researches on the theory of molding of sheet metal and the prediction of disfigurement, and use Autoform software to simulate theforming processthat based on the dynamic display algorithm. At the same time, by adjusting the process parameters, the stamping process problem was solved.Keywords: roof outer plate; numerical simulation; stamping; springbackCLC NO.: U467 Document Code: A Article ID: 1671-7988(2019)14-145-04引言重卡驾驶室白车身的典型结构是由侧围外板、顶盖外板、后围外板等大型外表面覆盖件和内部骨架件连接而成,受制于板料成形的复杂性、整车布置、模具制造工艺和调试水平的影响,板料冲压过程中产生的起皱、变薄、开裂和回弹等缺陷的控制一直是车身制造的一个重点和难点。
CRH3动车组模型侧风下明线运行的空气动力特性分析
CRH3动车组模型侧风下明线运行的空气动力特性分析作者:魏雅洁张群刘新桥来源:《计算机辅助工程》2022年第03期摘要:为探究在强风作用下高速列车的空气动力性能,建立3车编组CRH3动车组在侧风下明线运行的三维计算模型,在FLUENT软件中对其空气流场进行数值模拟。
当列车运行速度为60 m/s时,风向角分别选取0°、±3°、±6°、±9°、±12°、±15°、±18°、±21°、±24°,计算列车表面压力分布和流场域空气流速变化情况。
结果表明:车体表面正压最大值和负压最大值与风向角呈现正相关性,压力最大值位于车头正对迎风侧部位,车体表面近壁区域空气流速也与风向角具有正相关性。
关键词:侧风; 风向角; 空气流场; 明线运行; 多工况; 空气动力系数中图分类号: U270.11; TB115.1文献标志码: BAerodynamic characteristics analysis of open-line operationof CRH3 EMU under crosswindWEI Yajie, ZHANG Qun, LIU Xinqiao(INTESIM (Dalian) Co., Ltd., Dalian 116023, Liaoning, China)Abstract: To study the aerodynamic performance of high-speed trains under strong wind, the three-dimensional calculation modelof CRH3 EMU with three cars running on open-line under crosswind is built, and its air flow field issimulated in FLUENT software. When the train speed is 60 m/s, the pressure distribution on train surface and the change of air velocity inflow field is computed by separatelyselecting the wind direction angles of 0°, ±3°, ±6°, ±9°, ±12°, ±15°, ±18°, ±21° and ±24°. The results show that there is a positive correlation between the maximum positive pressure and maximum negative pressure on the car body surface and the wind direction angle, and the maximum pressure is located at the windward side of the front of the vehicle, and the air velocity near the wall of the car body surface also has a positive correlation with the wind direction angle.Key words: crosswind; wind direction angle; airflow field; open-line operation; multi-work condition; aerodynamic coefficient0 引言随着列车运行速度的提高,列车与空气的相互作用变得十分强烈。
重型载货汽车复合空气悬架导向臂支架优化设计
重型载货汽车复合空气悬架导向臂支架优化设计马凤军;张克鹏【摘要】To eliminate failure of the support of compound air suspension guide arm of a heavy truck in road test, we construct a finite element model for support of the guide arm, and use FE method for static strength and fatigue life analysis. Optimization design on the support based on the continuum structure topology optimization technology is conducted, and strength and fatigue life analysis are taken on the optimized model. The results demonstrate that the support of optimized guide arm which is improved both in strength and fatigue life features lightweight, vehicle test proves reliability of the optimization scheme.%针对某重型载货汽车复合空气悬架导向臂支架在道路试验中失效的问题,建立该导向臂支架的有限元模型,应用有限元法对其进行静强度和疲劳寿命分析。
基于连续体结构拓扑优化技术对支架进行优化设计,并对改进后的模型进行静强度、疲劳寿命的计算。
结果表明,经过优化设计的导向臂支架强度和疲劳寿命都得到提高,质量也有所减轻,且经实车试验验证了优化方案的可靠性。
基于ANSYS的某重卡储气筒支架结构优化
基于ANSYS的某重卡储气筒支架结构优化谷同金;李本荣;顾鹏飞【摘要】在某重卡车型开发过程中出现储气筒支架开裂问题,运用ANSYS工具建立有限元模型进行强度分析,文章针对分析结果提出结构优化方案,并对优化结果进行有限元校核.【期刊名称】《汽车实用技术》【年(卷),期】2018(000)014【总页数】2页(P76-77)【关键词】重卡;支架;有限元【作者】谷同金;李本荣;顾鹏飞【作者单位】安徽江淮汽车集团股份有限公司海外汽车研究所,安徽合肥 230601;安徽江淮汽车集团股份有限公司海外汽车研究所,安徽合肥 230601;安徽江淮汽车集团股份有限公司海外汽车研究所,安徽合肥 230601【正文语种】中文【中图分类】U467.3前言在采用气压制动系统的车辆上,储气筒是必不可少的零部件,一般在布置储气筒时都是根据储气筒布置位置设计专用的支架。
支架应有足够的强度保证,一旦支架发生开裂,进行影响整车制动性能,容易引起交通事故。
本文结合产品使用环境,采用ANSYS分析发现储气筒支架强度薄弱区域,进一步优化设计,满足使用要求,提高客户满意度。
1 市场问题我司某出口南美车型在客户使用过程中出现储气筒支架开裂现象,该地区路况较差,且车辆超载情况严重。
下图为市场故障照片。
图1 市场故障件2 有限元分析2.1 建立模型根据该支架的几何模型建立有限元模型。
将支架几何模型导入 ANSYS,对支架进行网格划分。
支架几何模型及有限元模型如下图所示。
图2 储气筒支架几何模型图3 储气筒支架有限元模型2.2 模型计算该支架通过上端面两圆孔与车架横梁螺栓连接,所以设置边界条件时可将支架上端面进行全约束,同时考虑车架抖动及超载情况,施加载荷时设定3倍重力加速度。
有限元分析结果如下图所示。
图4 有限元计算结果图5 优化后支架模型从上图计算结果可以看出,支架与车架横梁连接处倒角附近应力水平较高,同时在倒角处出现应力集中现象。
应力分布与该支架市场问题状态基本一致。
2023年云南省起重机司机作业Q2证理论考试练习题
2023年云南省起重机司机作业Q2证理论考试练习题1、(推断题)起重机在安装修理或使用中操作不便,可任意转变电路。
参考答案:错误2、(推断题)起重机的主要参数包括:起重量、跨度(幅度)、起上升度、工作速度及工作级别等。
参考答案:正确3、(推断题)桥式起重机金属结构主要受力构件有主梁和端梁。
参考答案:正确4、(推断题)我国规定的平安电压的上限值为:在任何状况下,两导体之间或任一导体与地面之问均不得超过沟通有效值为50V o参考答案:正确5、(推断题)从事《特种设备平安监察条例》规定的监督检验、定期检验、型式试验以及特地为特种设备生产、使用、检验检测供应无损检测服务的特种设备检验检测机构,应当经国务院特种设备平安监督管理部门核准。
参考答案:正确6、(推断题)特种设备出厂时,依据用户要求随附设计文件、产品质量合格证明、安装及使用修理说明等文件。
参考答案:错误7、(推断题)操违反操作规程属于事物的担心全因素。
参考答案:错误8、(推断题)指挥人员发出的指挥信号必需清楚、精确。
参考答案:正确9、(推断题)特种设备生产、使用单位和特种设备检验检测机构,应当接受特种设备平安监督管理部门依法进行的特种设备平安监察。
参考答案:正确10、(推断题)电梯应当至少每30日进行一次清洁、润滑、调整和检查。
参考答案:错误11、(推断题)《特种设备平安监察条例》所称特种设备是指涉及生命平安、危急性较大的锅炉、压力容器(含气瓶)、压力管道、电梯、起重机械、客运索道、大型游乐设施和场(厂)内专用机动车辆。
参考答案:正确12、(推断题)桥式起重机当吊钩处于工作位置最低点时,钢丝绳在卷筒上的缠绕,除固定绳尾的圈数外,必需不小于2圈。
参考答案:正确13、(推断题)起重机上所用的润滑脂,必需留意与原机说明的牌号相符。
参考答案:正确14、(推断题)缓冲器在碰撞之前,机构应切断限位开关。
参考答案:正确15、(推断题)液控早向阀允许液流朝一个方向流淌,反向流淌必需满意肯定条件。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
基于Virtual Wind Tunnel的某重型牵引车外气动性能分析Aerodynamic Performance Analysis of a Heavy TractorBased on Virtual Wind Tunnel张克鹏(陕西重型汽车有限公司陕西西安710200)摘要:车辆在高速行驶过程中,大部分的动力都要用来克服空气阻力,而空气阻力与整车的风阻系数成正比。
因此,降低整车风阻系数就成了整车设计过程中考虑的主要因素。
文章利用Altair公司有限元软件HyperWorks中的Virtual Wind Tunnel模块进行某重型车整车外流场分析,相比较传统CFD分析前处理工具,Virtual Wind Tunnel 采用为汽车外流场定制的流程化建模策略,能够避免很多其他因素干扰。
根据分析结果提出增加导流罩来优化整车外气动性能,结果表明:增加导流罩后整车风阻系数较原来有显著降低,该方法为某重型车气动性能优化设计提供理论依据。
关键词: Virtual Wind Tunnel 重型牵引车气动性能流程化优化Abstract:Most of the power is used to overcome air resistance when vehicles are in the process of high speed, and the drag force is proportional to the wind drag coefficient. Therefore, to reduce the drag coefficient of the vehicle is the main factor in the design process of the vehicle. The flow field is analyzed by using Virtual Wind Tunnel (VWT) in HyperWorks suite of Altair Company. Compared with the traditional tools for CFDpre-processing, Virtual Wind Tunnel provides a highly automated modeling strategy specifically designed for vehicle external aerodynamics analyses. The results indicate that the coefficient of drag decreased significantly when adding the fairing, and thus improve the external aerodynamic performance of the vehicle. Key words:Virtual Wind Tunnel,heavy tractor,aerodynamic performance,process,optimized1概述近年来,全世界的汽车制造商都依赖于耗时的风洞试验和计算流体动力学(CFD)仿真来研究汽车的空气动力学性能[1]。
今天,快速的计算机系统和尖端的数值方法允许人们在短时间内来研究复杂的流动结构。
在汽车的研发过程中,风洞试验仍然是一个不可或缺的过程,同时CFD风洞仿真也越来越受欢迎,它的运用大大的减小了实际所需的风洞试验次数[2]。
完成汽车风洞试验仿真不是一个简单的任务。
从模型的准备,到网格划分和CFD求解器设置,再到计算和后处理,这整个工作流程是复杂且耗时的。
汽车外流场仿真的特点是几何(如发动机舱)、边界条件(如旋转轮胎和滑移地面)和复杂流场结构都比较复杂,特别是汽车的尾迹区。
然而,流程定制化的环境,加上可靠的精确性,以及可扩展和稳健性的CFD求解器,将成为风洞仿真高效和主流方式。
文章利用Altair公司高度流程定制化的CFD工具Virtual Wind Tunnel对某重型车进行整车外流场CFD模拟,并通过增加车辆辅助装置对其外气动性能进行优化,优化后该车型整车风阻系数有明显降低,为整车降油耗提供设计方向。
2 Virtual Wind Tunnel工具简介Virtual Wind Tunnel的核心技术是Altair先进的计算流体动力学(CFD)求解器——AcuSolve。
它是一款通用的、基于有限元技术的CFD求解器,因而在求解速度、可扩展性、精确性和稳健性方面有着独特的优势。
AcuSolve采用雷诺平均Navier-Stokes(RANS)和分离涡(DES)技术来模拟湍流,预测流场和流动分离。
DES 技术采用RANS方法模拟近壁边界层流动,而采用大涡模拟(LES)模拟车辆尾迹的分离流动。
采用RANS方法进行稳态分析和DES方法进行瞬态模拟均可以得到精确的外流场计算结果。
由于AcuSolve在时间推进上采用了一套非常高效和稳健的数值方法,瞬态分析的耗时较少,通常都在可接受的范围内。
这使得分析人员有充分的时间来采用更贴合实际的、更精确的瞬态仿真,而不必由于考虑到时间成本而采用稳态仿真来求解车辆的流场。
Virtual Wind Tunnel用户界面友好直观,一体化的操作环境及高度自动化的仿真过程,使得用户不容易犯错进而影响仿真结果,整个界面如图1所示。
另外的一些参数控制通过AcuSolve专用的前处理AcuConsole在后台自动完成,同时AcuConsole也是一个集成在HyperWorks中的模块。
因此,Virtual Wind Tunnel高度自动定制化的特点大大降低了工程设计人员进行整车外流场分析门槛的同时,也缩短了整个分析的时间,大大提升了工程效率。
图1 Virtual Wind Tunnel定制化界面3 模型建立及边界条件3.1 虚拟风洞的建立计算模型为某重型牵引车,采用CATIA建立其三维模型。
为了便于考察模型细节处的流场情况及保证其计算的准确性,尽量保证计算模型与实际三维模型的一致性。
保留了驾驶室上对外流场影响较大的雨刮器、后视镜、遮阳罩、挡泥板、进气格栅等在内的所有复杂部件,驾驶室本体进行一系列处理形成封闭几何体,并采用1:1实际模型进行数值计算。
采用Altair公司HyperWorks仿真平台的通用前处理工具HyperMesh进行几何清理及面网格划分,最终形成三角形网格数为2157054,节点数为1078299,重型牵引车面网格模型如图2所示。
图3为Virtual Wind Tunnel中进行虚拟风洞的定义。
进行整车外流场分析时,其模拟风洞要减小其阻塞效应,研究表明,低于1%的风洞试验结果,阻塞干扰产生的误差才不需要修正。
根据这一理论,计算域上部留6倍车高,左右各5倍车宽,从车辆空气动力学研究经验来看,汽车风洞试验过程中,汽车尾部有一个比较大的湍流区域,这个区域的流动非常紊乱,数值计算过程中,这个区域如果模拟得不好对结果会产生重大的影响。
所以车辆尾部一般设定8倍车长的区域,保证汽车湍流区域充分发展;对于汽车前部,由于设定的风洞入口的风速比较稳定,为了充分模拟汽车前部的空气流动,汽车前部区域也要留出一定的空间,前部取4倍车长[3-6]。
图2重型牵引车面网格模型图3 Virtual Wind Tunnel中虚拟风洞定义考虑到重卡的整车相对轿车来说较长,同时为了减小风洞的阻塞效应,该重型牵引车前部取L S为2倍车长,尾部流场较复杂,所以车辆尾部取5倍车长,上部留6倍车高,左右各5倍车宽,根据整车的长L T(16m)、宽(2.5m)、高(4.2m),定义风洞的尺寸;车速定义为100km/h,大约28m/s。
汽车在行驶过程中,轮胎有一定的弹性变形,实际车轮和地面的接触是一个面;为了模拟车辆行驶真实情况以及加速计算收敛,将地面和轮胎相交的部分用平板封闭,细节部分如图4所示。
图4 地面与轮胎相交部门细节图3.2 边界条件的设定在进行重型牵引车外流场的数值模拟过程中,必须对车辆的关键部件进行边界条件的设定。
车辆散热模块的散热器、中冷器等按照多孔介质来处理,具体参数根据图5、图6所示该车型散热器和中冷器的压降试验数据获得。
图5 散热器压降试验特性曲线图6 中冷器压降试验特性曲线将该型车的底盘车轮及挂车车轮都定义为旋转车轮,根据之前定义的28m/s的车速及车轮的半径,可以计算出车轮的转速为493Rpm。
除了热交换器和发动机外,都定义了边界层,层数为3层,如图7所示。
图7 定义边界层数3.3 定义加密区域及计算设定汽车尾流区是一个非常大的湍流流动区域,该区域空气流动形成漩涡,流动速度变化范围大。
在汽车前后区域、汽车尾部区域、汽车底部接地区域及汽车周围对网格进行加密处理。
这样的网格可以保证更好地对汽车周围的气流流动进行模拟,得到更加精确的数值模拟结果,这里进行3层加密并定义加密单元尺寸为5mm,如图8所示。
在提交计算之前,首先进行各项设置:包括来流速度、计算类型、计算步数、是否为移动地面、车轮是否旋转、网格细密程度、边界层数等,如图9所示。
图8 定义网格加密区域图9 定义求解设置4 原型车CFD计算结果分析本文中分析对象原型车模型中驾驶室后围与车厢距离为1200mm,没有驾驶室导流罩等辅助装置。
通过Virtual Wind Tunnel工具求得该重型牵引车风阻系数C d值为0.673。
图10为该模型表面压力等值云图。
从图10可知,当车辆行驶过程中,前方气流首先遇到驾驶室前部,在车辆驾驶室前部气流滞止,在前面罩、保险杠及挡风玻璃的大部分区域形成较大的正压,迫使气流向顶部、侧面和车辆底部分流。
由驾驶室上部分流的气流流过驾驶室前顶边缘和A柱及后视镜附近时,产生小的气流分离区域,产生负压区域,且压力梯度较大。
气流流过驾驶室后,直接冲击到车厢前部高出驾驶室顶部的前表面,使部分气流滞止,形成一个正压区。
图11为原型车在中心对称面处的速度等值线云图。
从图中可以看出驾驶室后围与车厢的间隙处气流速度较低,驾驶室后部产生了一个大尺度漩涡,气流在此处能量损失较大;流过车辆底部的气流速度较低,压强较高,尾部压降较低,从车辆底部流出到达车辆尾部的气流,在这个压强差的作用下上卷,在尾部形成较强的上卷漩涡,形成较强负压,这一负压是车辆压差阻力的重要来源,漩涡随着离开车辆的距离增大而逐渐减弱。
5 整车外气动阻力特性的改进通过上面分析,文中利用比较成熟的3种导流装置,进行整车气动阻力特性的改进。
(1)在驾驶室顶部加装顶导流罩,定义为A型。
图10 原型车表面压力等值云图图11 原型车速度等值线云图(中心对称面处)(2)侧面加装侧导流板,定义为B型。