起重机俯仰减速器声-结构灵敏度分析.docVIP

起重机俯仰减速器声-结构灵敏度分析 耿飞宇1,林腾蛟1，钟声2 （重庆大学 机械传动国家重点实验室，重庆 400044401120） 引言 灵敏度设计对声学参量的影响大小利用LMS Virtual.ab、NASTRAN等软件，以箱体焊接板的厚度材料的弹性模量为设计，以场点辐射声压为声学参量，建立灵敏度分析模型，。2 声-结构灵敏度分析理论灵敏度设计对声学参量的影响大小LMS Virtual.Lab求解俯仰减速器声-结构灵敏度问题中，关键是求出器特征值灵敏度和声传递向量（ATV）。 无阻尼声-结构耦合系统离散有限元方程可表示为 (1) (2) 式中：，，——结构的质量矩阵、声流体质量矩阵和声-结构耦合质量矩阵； ，，——结构的刚度矩阵、声流体刚度矩阵和声-结构耦合刚度矩阵； ——系统响应向量，其中为结构节点位移向量，为声压向量；——系统激振力向量。 且，质量矩阵和刚度矩阵有如下关系，。假设 ， (3) 将式(3)带入式(1)得到稳态频域方程 (4) 式中：Z——阻抗矩阵，； X——结构响应幅值向量，，其中，为结构节点响应位移幅值分量，为声学节点响应声压分量； ——激振力频率。 2.1 基于模态叠加法的结构灵敏度 对于振型相对质量矩阵M进行归一化处理的结构系统来说，节点j的位移响应为 (5) 式中：——第i阶模态的参与因子； ——第i阶模态在节点j引起的模态位移值； 将式(5)带入式(4)得到 (6) 由于K、M均为非对称矩阵，且和不关于K、M正交，故式(6)不能按照一般的模态法求解，可以定义一个左特称向量[48,49]满足 (7) 于是，有 (8) 用左乘式(7)，并联合式(8)，得第i阶模态的参与因子 (9) 这里：——第i阶固有圆频率； ——第i阶模态在激振力位置Ｌ下引起的模态位移值； 把式(9)对设计参量进行求导，得到模态参与因子对设计参量的灵敏度 (10) 式中：——特征向量灵敏度； ——特征值灵敏度。 二者总称为特征灵敏度。而结构位移响应的灵敏度可以表示为 (11) 因此，我们可以通过计算特征值灵敏度和特征向量灵敏度，方便的得出减速器箱体表面速度响应对设计参量的灵敏度 (12) 2.2 基于声传递向量的声学灵敏度 在小扰动情况下，可以认为声学方程是线性的，故声场中某点的声压与结构表面的振动量之间是线性关系，这种关系可以用声传递向量（ATV）来表示。将结构表面离散成有限个单元，场点的声压可表示为 (13) 式中：——声传递向量，其物理意义是：结构的单元或节点在特定频率下的单位速度在场点上引起的声压值。 因此，可通过声传递向量，将场点的声压与结构表面的振动速度响应之间建立联系。声传递向量与结构的几何形状、阻尼、计算频率及场点位置有关，针对特定的结构和场点，这一值为常数，于是，声压级响应的灵敏度为 (14) 2.3 声-结构灵敏度计算公式 声-结构灵敏度就是描述结构在已知的振动载荷与边界条件下，设计场点的辐射声压值对结构设计参量的导数，可用公式表示为 (15) 将式(12)、(14)代入到式(15)中，可得声-结构灵敏度的数学表达式 (16) 因此，只需求得特征值灵敏度、特征向量灵敏度和声学传递向量，就能通过式(16)得到减速器的声-结构灵敏度值。 3 俯仰减速器声-结构灵敏度分析模型 本文分别建立结构有限元模型和声学边界元模型，基于有限元法计算减速器特征值灵敏度和结构振动速度灵敏度，基于边界元法计算声传递向量（ATV）和箱体表面声压灵敏度，通过耦合结构灵敏度和声学灵敏度求解减速器的声-结构灵敏度，即各场点辐射噪声声压频响对结构设计参量的敏度大小。 3.1结构灵敏度分析模型 在UG软件中建立俯仰减速器三维实体模型，并对窥油孔、螺栓孔等细微结构予以简化或省略，将简化后的三维模型导入ANSYS中，以建立结构有限元模型。 对结构进行修改一般可选择几何参数和材料参数。因此，本文选择上、下箱体的10个焊接板厚度和箱体的弹性模量为设计参量，研究场点声压对各设计参量变化的敏度。结构有限元模型建立的具体方法是：在ANSYS中，用shell63壳单元对各焊接板