2002 年 MSC.Software 中国用户论文集
MSC.Nastran 在轻型客车 NVH 特性研究中的应用 马天飞 林逸 张建伟 (吉林大学)
摘要 : 本文以某轻型客车为研究对象,利用 MSC.Nastran 软件建立了用于整车 NVH 特性研究的有限元模型,并介绍了 MSC.Nastran 软件在 NVH 特性研究过程中的 具体应用。
1 前言 NVH 指的是 Noise(噪声)、Vibration(振动)和 Harshness(声振粗糙度), 由于它们在车辆等机械中是同时出现且密不可分的,因此常把它们放在一起进行 研究。对于汽车来说,乘员的一切来自于车的触觉和听觉感受都属于 NVH 特性研 究的范畴,表现为乘员的乘坐舒适性,而噪声又是其研究的重点。调查表明,目 前人们对乘坐舒适性的要求明显提高,这就使以改善汽车乘坐舒适性为目的的汽 车 NHV 特性的研究变得更加重要[1]。 汽车 NVH 特性的研究方法有多种,其中有限元方法(FEM)和多体(MB)系统 动力学方法相结合的分析方法较为成熟。利用它们建立的刚体、弹形体混合模型可 用于车内低频噪声(20~200Hz)的预测。车室内部的噪声预测是建立在对车身结 构和封闭空腔之间声固耦合作用(FSI)模拟的基础之上的,其模型必须能够反 映出车身结构和车室内空气的动力学特性和两者在边界上的相互作用关系[2]。 利用 MB/FEM 仿真技术确定出作用于车身上的力,再通过 FSI 分析求出车室内 部声压级,就可以评价不同工况下车室内部的结构噪声了。图 1 表示出了汽车 NVH 特性研究的一般过程。 细化车身 FE 模型
底盘 MB 模型
车身超单元模型
模拟激励
车室空腔 FE 模型
车身结构 FE 模型
整车模型
车身连接点的力
声固耦合 FE 模型
2002 年 MSC.Software 中国用户论文集
车室声学响应
图 1 整车 NVH 特性研究过程 作为大型通用有限元软件,MSC.Nastran 在整车的 NVH 特性研究中可以起到 重要的作用,被许多汽车公司所采用。下面以某轻型客车为研究对象,具体介绍 MSC.Nastran 软件在整车 NVH 特性研究中的应用。
2 MSC.Nastran 在整车建模中的应用 在研究整车 NVH 特性的过程中,为了准确模拟车室内的噪声水平,必须确定 车室声固耦合模型在各种工况下所受到的激励力,为此必须建立准确的整车模型。 下面首先利用 MSC.Nastran 软件建立细化的白车身模型,再利用超单元技术将车 身有限元模型作为柔体转换到多体软件 MSC.ADAMS 中。 2.1 白车身有限元模型 白车身是车身的主体部分,其有限元模型占据着绝大多数的自由度,因此在 分析中具有极为重要的地位。图 2 是利用 MSC.Patran 对某轻型客车的白车身 CAD 模型划分网格,建立的有限元模型。整个模型共有 24882 个节点,29593 个单元。 建模过程中对一些不重要的特征(如:小孔等)进行适当地简化,以减少模型的 自由度,便于计算。
2002 年 MSC.Software 中国用户论文集
图 2 白车身有限元模型
建立了白车身的有限元模型之后,可以对它进行模态分析,以获得它的共振 频率,再通过修改车身结构避开激励频率,从而防止产生共振。另外,通过模态 振型可以判断出车身变形较大的部位,从而可以有的放矢地改进车身刚度,减少 振动噪声的产生和传递。MSC.Nastran 的 SOL103 求解器可以高效稳定地进行结构 模态分析。图 3 为模态分析得到的车身第一阶扭曲模态(图 a.b)和第一阶弯曲模 态(图 c.d)。 在白车身有限元模型中再加上其它车身部件:如车门、发动机罩、座椅、内饰件 等,就可以得到与实际结构尽可能相近的详细车身模型,从而更加准确地模拟实 际的车身结构[3]。
(a)
(b)
2002 年 MSC.Software 中国用户论文集
(d)
(c)
图 3 白车身的模态分析结果
2.2 车身超单元模型 在建立整车模型时,为了有效地将有限元模型与多体模型结合起来,可 以利用车身模型建立超单元(Superelement)模型,再作为多体模型中的一 个部件合并到底盘多体模型中。 有限元分析中的超单元是指子结构(实际结构的一部分)的有限元模型 , 在数学模型中表现为各种物理量(如:质量、刚度、载荷等)从内点向外点或 边界点缩减得到的边界矩阵[4]。在超单元分析中,子结构的内点将被忽略,对 外只表现为外部节点的自由度。在求解出缩减得到的剩余结构(外部节点)的 位移之后,再求解内部节点的位移。 在 MSC.Nastran 中建立柔体超单元模型时保留下的外部节点就是合并到 MSC.ADAMS 的多体模型中时的连接点。在 MSC.ADAMS/FLEX 中可以在连接点处建 立各种边界条件(如:铰链或力)。 在 MSC.ADAMS/View 中可以接受的柔体模型是以模态中性文件(MNF)的形 式出现的。在 MNF 中保留着柔体的全部模态信息,它是由 NASTRAN 的输出文件 (*.out)经过 MSC.ADAMS/FLEX 中的转换工具的转换而产生的。转换工具是由 mnfx.alt 和 msc2mnf.exe 两个程序组成的:mnfx.alt 是用于产生 MSC 中间文 件的 DMAP 转换程序,它利用 NASTRAN 中的部件模态合成的超单元技术在 NASTRAN 输 出 文 件 中 提 取 所 需 的 柔 体 信 息 , 并 写 成 二 进 制 文 件 ; 而 msc2mnf.exe 则读此文件并通过 MSC.ADAMS/FLEX 中 MNF 工具箱中的一系列程 序写出相应的 MNF。 图 4 为 利 用 该 客 车 的 车 身 模 型 转 化 得 到 的 超 单 元 模 型 形 成 的 MNF 在 MSC.ADAMS/View 中建立的柔体模型。将它连接到在 MSC.ADAMS/View 中建立的 底盘多体模型中,就可以得到整车模型。最后在 MSC.ADAMS/View 中施加模拟 激励(如:路面谱、凸块冲击载荷等),测量出这些工况下底盘与车身之间连
2002 年 MSC.Software 中国用户论文集
接点的响应力,并作为车室声学分析的激励。
图 4 转换到 MSC.ADAMS/View 中的车身超单元模型
3 MSC.Nastran 在建立车室声固耦合模型中的应用 汽车车室构成封闭空腔,形成一个声学系统,其声学特征表现为与固有 频率和振型(压力分布)相联系的声学振动模态
[5]
。汽车 NVH 特性研究的重
要内容就是要模拟预测车室内的噪声。 空腔声学模态是通过边界条件与车身结构的振动相耦合的,这种边界条 件(开放的、刚性的或弹性的等)建立了车室内空气运动和车身壁板振动运动 之间的关系,因此空腔声学特性和车身结构动力学特性共同决定了车室内部 的声压。 试验表明:车身壁板的振动会改动声学振型的频率,移动节线(声压为 零的界线),并使车室内的噪声响应发生重大变化。为了准确模拟车身板件与 声腔中空气之间的相互作用,就要把车身结构视为弹性体,在 MSC.Nastran 中建立车室空腔与车身结构的声固耦合模型,并在分析耦合系统振动模态的 基础上进一步进行声学响应分析。 在声固耦合模型的动力学方程中,由于结构和空气在界面的相互作用而 导致质量和刚度阵中引入了对角线外的耦合项,这样在对系统进行模态分析 时将出现复特征值,在应用 MSC.Nastran 进行模态分析时应采用 SOL107 求解 器。 在 MSC.Nastran 中建立声学有限元模型时,可利用 MSC.Patran 提供的 HEXA 等实体单元建立三维的声学单元,并将相应的节点数据卡的第 7 域添加1,将它改为流体单元节点。在建立声固耦合模型时需要在耦合界面上将声学 单元的节点与结构单元的节点一一对应[6]。
2002 年 MSC.Software 中国用户论文集
声学单元的理想尺寸大约是每个波长六个单元,实际上通常采用的声学 单元的长度一般为 0.1~0.2 米。第一阶声学模态一般为出现在 40~80Hz 左右的 纵向声学模态,空腔越长频率越低。在空腔中加入座椅将使它的有效长度增加, 从而降低了它的声学模态频率。 在建立声固耦合模型时,一般是先建车室空腔模型,然后再建车身结构 模型,最后再把它们耦合起来。如果建模时先从结构模型入手,将使建模工作 增加不必要的麻烦。 这时的车身模型在结构上比前面提到的详细车身要简单一些,但为了准 确模拟声学响应,应在模型中加上对声学模态影响较大的座椅、内饰件等。如 果只为了定性地研究声固耦合作用对车室空腔声学模态的影响,可以将模型 进一步地简化。图 5 所示为该客车车室简化的耦合模型,模型中共有声学单元 30 个,节点 75 个;结构单元 66 个,节点 68 个。结构单元的节点与声学单元 边界的节点全部耦合。
图 5 简化的客车声固耦合模型 图 6~9 是在 MSC.Nastran 中分析得到的车室空腔的声学模态及其耦合模型 的声学模态。其中 a、b 图为自由边界的空腔声学模态,c 图为近似的刚性壁空 腔的声学模态,d 图为 3 毫米钢板柔性壁空腔的耦合声学模态。
2002 年 MSC.Software 中国用户论文集
a. 36.663Hz b. 36.663Hz
c. 36.671Hz
d. 38.796Hz
图 6 声固耦合模型的第一阶模态 由图 6 可以看出:由于该车的车室较长,其自由边界空腔的第一阶模态频 率较低,为 36.663Hz;而对于刚性壁模型来说,由于声场的边界被约束,其 声学模态的振形变化较大,但模态频率的变化不大;在柔性壁模型中,它的 模态振形与刚性壁模型的很相似,但在车室的中部其声压分布已经出现较大 的变化。另外,柔性壁耦合模型的模态频率也发生了变化。
a. 74.672 Hz b. 74.672Hz
c. 74.669Hz
d. 75.973Hz
图 7 声固耦合模型的第二阶模态 图 7 所示的第二阶声学模态仍然是纵向模态。由于模型的几何形状比较规
2002 年 MSC.Software 中国用户论文集
则而且在纵剖面上形状没有变化,因此它的模态频率基本上是第一阶频率的 二倍,其声压分布在横向上基本没有变化。另外,与第一阶模态相比较,柔性 壁、刚性壁模型之间模态频率和模态振形的变化变小了。
a. 108.73 Hz b. 108.73Hz
c. 108.75Hz
d. 110.61Hz
图 8 声固耦合模型的第三阶模态
车室空腔的第一个横向声学模态出现在 108.73Hz,由于车室的长宽比比 较大,因此它位于第二阶纵向模态之后。由图 8 可以看出:由于模型是左右对 称的,因此它的模态振形也均匀对称。
a. 114.74Hz b. 114.74Hz
2002 年 MSC.Software 中国用户论文集
c. 114.75Hz
d. 116.32Hz
图 9 声固耦合模型的第五阶模态 在第五阶声学模态(如图 9 所示)中,车室空腔的纵向模态与横向模态相 互作用,其综合振型使得声压在车室空腔内的分布更加复杂,由于它的模态 频率也比较低,对车室内声场的计算影响比较大,因此在预测车室内的噪声 时应予以充分的重视。 在 MSC.Nastran 中建立了汽车车室的声固耦合模型并在 MSC.ADAMS 软件中 计算得到了它的激励力,就可以利用 MSC.Nastran 中的 SOL108、109 等求解器 进行车室内的声学响应分析。
4 结论 在汽车低频 NVH 特性的研究中,利用 MSC.Nastran 软件能够方便地建立 车身结构模型并顺利完成车室的声固耦合分析等工作。另外,使用它可以方便 快捷地实现有限元模型与多体模型的结合,使计算机仿真工作能够更加高效 地为汽车 NVH 特性的研究提供服务。
参考 文献: 1 林逸,马天飞,姚为民,张建伟. 汽车 NVH 特性研究综述. 汽车工程 2002.3 2 D.Griffiths. Vehicle NVH Design Guide. Ford Motor Company 1993 3 顾彦,何波. MSC.Nastran 在汽车 NVH 中的应用.2001 年 MSC.Software 中国用户 论文集 4 MSC.Nastran Superelement Analysis SEMINAR NOTES. July 1999 5 [美]M.M.凯墨尔,J.A.沃尔夫 编 陈砺志 译, 现代汽车结构分析,人民交通出版社 1987.9 6 Fluid—Structure Analysis using MSC.Nastran. July 1999