CN104504215A - 基于单位分解“有限元-无网络”单元的汽车车内声场预测方法 - Google Patents
基于单位分解“有限元-无网络”单元的汽车车内声场预测方法 Download PDFInfo
- Publication number
- CN104504215A CN104504215A CN201510005758.9A CN201510005758A CN104504215A CN 104504215 A CN104504215 A CN 104504215A CN 201510005758 A CN201510005758 A CN 201510005758A CN 104504215 A CN104504215 A CN 104504215A
- Authority
- CN
- China
- Prior art keywords
- vehicle interior
- sound field
- interior sound
- automobile
- automobile vehicle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Fittings On The Vehicle Exterior For Carrying Loads, And Devices For Holding Or Mounting Articles (AREA)
Abstract
本发明公开了一种基于单位分解有限元-无网格单元的汽车车内声场预测的方法,包括如下步骤:利用汽车车内空气声腔网格模型;利用一层相连节点关系定义汽车车内空气声腔网格模型的单元节点支撑域;在汽车车内声腔网格模型的节点支撑域内利用径向点插值法构造局部近似函数,同时在汽车车内声腔的全局网格范围内利用等参单元形函数近似;按等参单元积分方式构造汽车车内声腔模型的高斯积分点处形函数及其梯度;构造汽车车内声腔问题声学方程的伽辽金弱形式,将得到汽车车内声腔网格模型的系统方程,进而预测汽车车内声腔模态频率及振型,同时预测发动机激励下汽车车内声腔的声学响应,并对计算结果进行评价。本发明方法在汽车车内声场振动噪声预测问题中能够得到比有限元法更好的计算结果。
Description
技术领域
本发明涉及汽车车内声场计算的预测方法,具体涉及一种基于单位分解“有限元-无网格”单元的汽车车内声场计算预测方法。
背景技术
目前,汽车的车内声腔声场预测分析是车身NVH(Noise、Vibration、Harshness,噪声、振动与声振粗糙度)性能CAE(computer Aided Engineering,计算机辅助工程)分析中的一项重要工作,它对车身NVH性能预测以及指导NVH性能开发都有十分重要的意义。
目前,汽车的车内声场仿真预测方法有解析法和数值法,其中解析法一般只针对简单问题和模型;数值法是目前主流的预测方法,包括有限元法和边界元法,并在许多商业软件如Nastran、Virtual.LAB等中广泛应用。许多工程师对这些方法使用比较娴熟,在汽车车内声场性能预测与分析中广泛采用。但这类方法也存在一些问题,有限元法存在声学数值色散效应,预测结果受模型网格尺寸大小和计算频率的高低的影响较大,而边界元法的系统矩阵是非对称满阵、计算效率低下,矩阵的计算复杂,其应用前景还不如有限元法广阔。
为满足汽车的车内声场分析问题预测结果的可靠性,汽车设计工程师需要在建模过程中对网格尺寸、模型质量等进行检查,这需要耗费较多的前处理时间和人力,因此,有必要对发明一些新的预测方法以降低预测结果受模型质量的影响。另外,由于汽车设计工程师们已经习惯了运用有限元网格模型进行分 析,这就要求新的预测方法最好是基于有限元网格模型,这样能更好地保证模型的通用性。
发明内容
为了克服上述现有技术中存在的缺陷,本发明的目的是提供一种基于单位分解有限元-无网格单元的汽车车内声场预测方法,该方法在形函数构造时针对汽车车内声场模型的全局域采用六面体的等参单元形函数近似,在其模型节点的支撑域内采用径向点插值法近似,形成混合的汽车车内声场模型的有限元-径向点插值形函数,该形函数集成了RPIM和FEM形函数的优点,能降低数值色散误差,提高了汽车车内声场分析的计算精度和分析频率范围。
为了实现本发明的上述目的,本发明提供一种基于单位分解有限元-无网格单元的汽车车内声场预测方法,包括如下步骤:
2、包括如下步骤:
S1、建立汽车车内声腔网格模型,即对汽车车内声腔模型进行网格划分,网格单元采用六面体单元,每个六面体单元包含8节点;
S2、利用一层相连节点关系定义汽车车内声腔模型的单元节点支撑域;
S3、构造汽车车内声场点的近似函数;
S4、构造汽车车内声场模型六面体单元形函数的完整形式及声压梯度形式;
S5,构造基于单位分解有限元-无网格单元的汽车车内声场预测模型
首先利用标准的Galerkin方法得到汽车车内声场预测模型的伽辽金弱形式;接着对其进行离散处理,得到汽车车内声场预测模型的离散方程,从而构造基于单位分解有限元-无网格单元的汽车车内声场预测模型;
S6,利用基于单位分解有限元-无网格单元的汽车车内声场预测模型对汽 车声腔的模态频率、振型进行计算;
S7,利用基于单位分解有限元-无网格单元的汽车车内声场预测模型可对车内声场的声学响应进行仿真预测。
本发明的有益效果在于;
(1)本发明采用基于单位分解有限元-无网格单元的形函数来改善车内声场预测模型的硬度,降低其预测过程中的数值误差,获得比有限元法更好的计算结果,有限元-无网格法对汽车车内声场模型的网格质量要求比有限元法要低,分析频率带宽更宽。这一优点为本发明解决汽车的车内声场预测问题提供了技术基础。利用本发明进行分析汽车车内声场预测时,采用基于单位分解有限元-无网格单元的形函数降低了车内声腔模型的硬度,减少数值色散效应,提高了车内声场频率响应分析的精度和分析频率范围,进而获得更准确的预测结果。
(2)本发明在汽车车内声场分析问题中,能够得到更好的计算效果,更高有效的分析频带宽,而且该方法对模型的质量要求更低(如网格模型粗糙),这样能降低更多的前处理时间,在汽车产品的开发设计应用中前景广阔。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是本发明基于单位分解有限元-无网格单元的汽车车内声场预测方法的流程图;
图2是本发明一种优选实施方式中汽车车内声场六面体网格示意图;
图3是本发明一种优选实施方式中汽车车内声场模型节点支撑域的选择示意图;
图4是本发明一种优选实施方式中汽车车内声场模型示意图;
图5是本发明一种优选实施方式中汽车车内声场频率响应预测结果。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
本发明中汽车车内声场预测问题采用传统的有限元网格模型,并采用基于单位分解的有限元-无网格单元形函数对汽车车内声场进行近似处理,将混合有限元-无网格单元形函数代入至汽车车内声场波动方程的伽辽金弱形式中,获得基于单位分解的有限元-无网格单元的车内声场离散波动方程,计算汽车车内声场模态频率及其振型,并根据汽车车内声场的激励边界计算汽车车内声场的频率响应,对汽车车内声场问题进行预测。
本发明提供了一种基于单位分解有限元-无网格单元的汽车车内声场预测方法,如图1所示,下面结合图1对本发明的预测方法进行详细说明。
本发明所述的基于单位分解的有限元-无网格单元的车内声场预测方法包括如下步骤:
S1,建立汽车车内声腔网格模型。
本步骤可利用传统的前处理软件(如Hypermesh)对汽车车内声腔模型进行网格划分,网格单元采用六面体单元,每个六面体单元包含8节点。网格尺寸为100mm×100mm×100mm。如图2所示的汽车车内声场的六面体网格模型,它由420个六面体单元组成。
S2,选择汽车车内声场模型各节点的支撑域。
本发明中采用径向点插值法对汽车车内声场点进行局部近似,首先要构造模型各节点的支撑域。本发明中利用一层相连节点关系来定义汽车车内声场模型的节点的支撑域。如图3所示,为了得到点i的支撑域,首先选择其自身,另外还需要选择一层相连的节点从编号1至26。
S3,构造汽车车内声场点的近似函数。
本步骤利用有限元法和径向点插值法构造汽车车内声场点的声压近似函数。首先在车内声场全局范围内利用三维有限元等参单元形函数进行全局近似,再利用径向点插值法对S2中的汽车车内声场点支撑域进行局部近似。在进行局部近似中,优先选择耦合多项式基的径向点插值法进行近似,并通过已知的支撑域节点声压向量计算径向点插值函数的各项系数。步骤3的详细过程如下:
S3.1,构造汽车车内声场点的声压变量的局部近似函数。
利用径向点插值法对S2中的汽车车内声场点的支撑域进行声压变量的局部近似,并获得声压变量的局部近似函数。在本发明的一种优选实施中,采用耦合多项式基的径向点插值法对汽车车内声场点声压变量进行局部近似。
首先,对汽车车内声场内任意点x,近似声压ph(x)可以通过以下插值方式得到:
式中n表示汽车车内声场点x的局部支撑域所包含的节点数;pi表示其节点xi(i=1,2,...,n)的声压值;Φi为相应的形函数,它由车内声场节点的支撑域内所有节点构造。
接着,利用耦合多项式基的径向点插值法,对汽车车内声场的声压进行近 似,得到:
式中Ri(x)是汽车车内声场节点的径向基函数,n是汽车车内声场节点的径向基函数的项数,qj(x)是多项式基函数,m是多项式基函数的项数。汽车车内声场节点的径向基函数采用多二次径向基函数,可以写成:
Ri(x)=[(x-xi)2+(y-yi)2+(z-zi)2+(αcdc)2]q,
式中q和αc为形状参数,分别为1.03和0.5;dc为汽车车内声场单元的等效长度。汽车车内声场模型的线性多项式基向量可以写成如下形式:qT(x)=[1 x y z]。
通过车内声场节点的支撑域内节点声压值可计算各项系数,得到车内声场内任一点的声压变量形函数表达式:
式中为矩阵G-1中的元素,其中 其中Rq为径向基函数的动差矩阵,可以写成:
矩阵qm表示为:
是qm矩阵的转置。
至此,汽车车内声场节点i的声压近似函数由RPIM近似得到,可以写成:
式中Φi为由汽车车内声场节点i通过RPIM构成的形函数;为汽车车内声场节点i的支撑域内各节点声压值所组成的向量;n为其节点i的所有支撑节点数。值得注意的是,汽车车内声场单元内所有节点的支撑点数目并不一致,n随汽车车内声场中节点i的变化而不同。
S3.2,构造汽车车内声场点的声压向量全局近似函数。
利用传统的有限元法对汽车车内声场网格模型进行分析,在全局范围内构造六面体的等参单元形函数作为汽车车内声场网格模型的全局近似函数。
汽车车内声场内任意点的声压p(x,y,z),
p(x,y,z)=Npe,
式中N为车内声场模型的三维有限元等参单元形函数。向量pe表示相应的车内声场模型三维单元的节点声压函数pi(x,y,z)(i=1,2,...,8),可表示为:
pe={p1(x,y,z) p2(x,y,z) p3(x,y,z) p4(x,y,z)
,
p5(x,y,z) p6(x,y,z) p7(x,y,z) p8(x,y,z)}T
式中pi(x,y,z)|i=1,2...8为车内声场模型节点局部近似函数,可在S3.1中获得;
S4,构造汽车车内声场模型六面体单元形函数的完整形式及声压梯度形式。
根据S3中汽车车内声场点的声压变量的近似函数形式,结合有限元法中单元形函数的组成形式,构造基于单位分解有限元-无网格单元形函数的完整形式。根据完整的单元形函数推导其声压梯度形式。
S4.1,构造汽车车内声场模型六面体单元形函数的完整形式。
利用S3中汽车车内声场点的声压变量的近似函数形式,将六面体单元中各节点组合起来,按有限元法中单元形函数的组成方式构造基于单位分解有限元-无网格单元的完整形函数。
在汽车车内声场模型中,基于单位分解有限元-无网格单元的形函数完整形式可写成:
矩阵Φ8×M由Φi(i=1,2,3,…,8)组装得到,列数M等于车内声场单元支撑域的节点数。为方便起见,前8列分别对应单元8个节点;其它列对应于汽车车内声场对于单元的支撑域节点。
如前所述,汽车车内声场内任意点的声压p(x,y,z)即可写成:
因此,基于单位分解有限元-无网格单元的汽车车内声场预测模型的形函数矩阵可以写成:
S4.2,构造汽车车内声场模型六面体单元的声压梯度计算形式。
根据S4.1中完整形函数公式,利用梯度计算方法求解其微分形式,即得到汽车车内声场模型六面体单元的声压梯度计算形式。
对形函数矩阵求微分,得到基于单位分解有限元-无网格单元的汽车车内声场预测模型中声压梯度矩阵:
Ψ,x=[Ψ1,x Ψ2,x Ψ3,x, … ΨM,x]≡N,xΦ+NΦ,x
Ψ,y=[Ψ1,y Ψ2,y Ψ3,y, … ΨM,y]≡N,yΦ+NΦ,y
Ψ,z=[Ψ1,z Ψ2,z Ψ3,z, … ΨM,z]≡N,zΦ+NΦ,z
S5,构造基于单位分解有限元-无网格单元的汽车车内声场预测模型。
首先利用标准的Galerkin方法可得到汽车车内声场预测模型的伽辽金弱形式;接着对其进行离散处理,得到汽车车内声场预测模型的离散方程,这样就构造了基于单位分解有限元-无网格单元的汽车车内声场预测模型。
利用标准的Galerkin方法,构造汽车车内声场预测模型的伽辽金弱形式如下:
在整个汽车车内声场模型内对上述伽辽金弱形式进行离散处理,得到基于单位分解有限元-无网格单元的汽车车内声场预测模型:
(K-k2M+iρωC)p=F
式中k和ω分别表示车内声场预测问题的波数和圆频率;K为其刚度矩阵,写成:K=∫ΩBTBdΩ。其中B为车内声场预测模型的声压梯度矩阵,可以写成:
M为汽车车内声场模型的质量矩阵:Μ=∫ΩΨTΨdΩ;C为汽车车内声场的声学阻尼矩阵:
F为汽车车内声场所受的激励向量:P为车内声场边界上的节点激励向量:P={p1,p2,…,pn}T。
S6,利用基于单位分解有限元-无网格单元的汽车车内声场预测模型对汽车声腔的模态频率、振型进行计算。
根据S5中汽车车内声场预测模型中K矩阵和M矩阵,按公式(K-λiM)Φi=0,求出特征值为λi和特征向量Φi。其中λi=(ω/c)2为汽车车内声场的模态频率,而特征向量Φi即汽车车内声场的模态振型。在计算过程中需要定义汽车车内的材料属性:波速c=343m/s,密度ρf=1.225kg/m3。运用基于单位分解有限元-无网格单元的汽车车内声场预测方法计算其模态频率和模态振型,同时给出了参考结果,参考结果则由精细汽车车内声场网格模型(3264六面体)通过商业软件Nastran的数值实验得到,可以视为近似结果,如表1所示。
表1
S7,利用基于单位分解有限元-无网格单元的汽车车内声场预测模型可对车内声场的声学响应进行仿真预测。
在S6的基础上,本步骤利用基于单位分解有限元-无网格单元的汽车车内声场预测模型对其声学响应进行预测。首先设置汽车车内声场预测模型的边界条件:声阻抗系数An=0.000144m/(Pa.s),载荷激励为发动机传递汽车车内声腔的振动速度V=0.001m/s,幅值频率范围50-150Hz,如图4所示。运用基于单位分解有限元-无网格单元的汽车车内声场预测方法计算驾驶员右耳处的声压频率响应,频率范围50-150Hz,步长1Hz。上述预测结果如图5所示。
从图5和表1中可以看出:基于单位分解有限元-无网格单元的汽车车内声场预测方法在分析汽车车内声场问题时具有比有限元法更高的精度,其结果更接近参考结果值,表明基于单位分解有限元-无网格单元的汽车车内声场预测方法能获得良好的分析结果。
基于以上方法介绍和数值实验分析,证明了基于单位分解有限元-无网格单元的汽车车内声场预测方法非常适合于汽车车内声场问题的预测,其预测精度比传统方法要高;对汽车声场模型的质量要求更低。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同物限定。
Claims (4)
1.一种基于单位分解“有限元-无网格”单元的汽车车内声场预测方法,其特征在于,包括如下步骤:
S1、建立汽车车内声腔网格模型,即对汽车车内声腔模型声腔模型进行网格划分,网格单元采用六面体单元,每个六面体单元包含8节点;
S2、利用一层相连节点关系定义汽车车内声腔模型的单元节点支撑域;
S3、构造汽车车内声场点的近似函数;
S4、构造汽车车内声场模型六面体单元形函数的完整形式及声压梯度形式;
S5、构造基于单位分解“有限元-无网格”单元的汽车车内声场预测模型;
首先利用标准的Galerkin方法得到汽车车内声场预测模型的伽辽金弱形式;接着对其进行离散处理,得到汽车车内声场预测模型的离散方程,从而构造基于单位分解有限元-无网格单元的汽车车内声场预测模型;
S6、利用基于单位分解“有限元-无网格”单元的汽车车内声场预测模型对汽车声腔的模态频率、振型进行计算;
S7、利用基于单位分解“有限元-无网格”单元的汽车车内声场预测模型对车内声场的声学响应进行仿真预测。
2.根据权利要求1所述的基于单位分解“有限元-无网格”单元的汽车车内声场预测方法,其特征在于,步骤3通过以下方式实现:
S3.1、构造汽车车内声场点的声压变量的局部近似函数;
利用径向点插值法对S2中的汽车车内声场点的支撑域进行声压变量的局部近似,并获得声压变量的局部近似函数;
首先,对汽车车内声场内任意点x|(x=x,y,z),近似声压ph(x)可以通过以下插值方式得到:
式中n表示汽车车内声场点x的局部支撑域所包含的节点数;pi表示其节点xi(i=1,2,...,n)的声压值;Φi为相应的形函数,它由车内声场节点的支撑域内所有节点构造;
接着,利用耦合多项式基的径向点插值法,对汽车车内声场的声压进行近似,得到:
式中Ri(x)是汽车车内声场节点的径向基函数,n是汽车车内声场节点的径向基函数的项数,qj(x)是多项式基函数,m是多项式基函数的项数;汽车车内声场节点的径向基函数采用多二次径向基函数,可以写成:
Ri(x)=[(x-xi)2+(y-yi)2+(z-zi)2+(αcdc)2]q,
式中q和αc为形状参数,分别为1.03和0.5;dc为汽车车内声场单元的等效长度;汽车车内声场模型的线性多项式基向量可以写成如下形式:qT(x)=[1 x y z];
通过车内声场节点的支撑域内节点声压值可计算各项系数,得到车内声场内任一点的声压变量形函数Φk(x)(k=1,2,...,n)表达式:
式中为矩阵G-1中的元素,其中 其中Rq为径向基函数的动差矩阵,可以写成:
矩阵qm表示为:
是qm矩阵的转置;
至此,汽车车内声场节点i的声压近似函数由RPIM近似得到,可以写成:
式中Φi为由汽车车内声场节点i通过RPIM构成的形函数;为汽车车内声场节点i的支撑域内各节点声压值所组成的向量;n为其节点i的所有支撑节点数;
S3.2、构造汽车车内声场点的声压向量全局近似函数
利用传统的有限元法对汽车车内声场网格模型进行分析,在全局范围内构造六面体的等参单元形函数作为汽车车内声场网格模型的全局近似函数;
汽车车内声场内任意点的声压p(x,y,z),
p(x,y,z)=Npe,
式中N为车内声场模型的三维有限元等参单元形函数;向量pe表示相应的车内声场模型三维单元的节点声压函数pi(x,y,z)(i=1,2,...,8),可表示为:
pe={p1(x,y,z) p2(x,y,z) p3(x,y,z) p4(x,y,z)
,
p5(x,y,z) p6(x,y,z) p7(x,y,z) p8(x,y,z)}T
式中pi(x,y,z)|i=1,2...8为车内声场模型节点局部近似函数,可在S3.1中获得;
3.根据权利要求1所述的基于单位分解有限元-无网格单元的汽车车内声场预测方法,其特征在于,步骤4通过以下方式实现:
S4.1、构造汽车车内声场模型六面体单元形函数的完整形式
利用S3中汽车车内声场点的声压变量的近似函数形式,将六面体单元中各节点组合起来,按有限元法中单元形函数的组成方式构造基于单位分解有限元-无网格单元的完整形函数;
在汽车车内声场模型中,基于单位分解“有限元-无网格”单元的形函数完整形式可写成:
矩阵Φ8×M由Φi(i=1,2,3,…,8)组装得到,列数M等于车内声场单元支撑域的节点数;
汽车车内声场内任意点的声压p(x,y,z)即可写成:
因此,基于单位分解“有限元-无网格”单元的汽车车内声场预测模型的形函数矩阵可以写成:
S4.2、构造汽车车内声场模型六面体单元的声压梯度计算形式
根据S4.1中完整形函数公式,利用梯度计算方法求解其微分形式,即得到汽车车内声场模型六面体单元的声压梯度计算形式,对形函数矩阵求微分,得到基于单位分解“有限元-无网格”单元的汽车车内声场预测模型中声压梯度矩阵Ψ,x、Ψ,y和Ψ,z:
Ψ,x=[Ψ1,x Ψ2,x Ψ3,x, …ΨM,x]≡N,xΦ+NΦ,x
Ψ,y=[Ψ1,y Ψ2,y Ψ3,y, …ΨM,y]≡N,yΦ+NΦ,y。
Ψ,z=[Ψ1,z Ψ2,z Ψ3,z, …ΨM,z]≡N,zΦ+NΦ,z
4.根据权利要求1所述的基于单位分解“有限元-无网格”单元的汽车车内声场预测方法,其特征在于,步骤5通过以下方式实现:
利用标准的Galerkin方法,构造汽车车内声场预测模型的伽辽金弱形式如下:
在整个汽车车内声场模型内对上述伽辽金弱形式进行离散处理,得到基于单位分解有限元-无网格单元的汽车车内声场预测模型:
(K-k2M+iρωC)p=F
式中k和ω分别表示车内声场预测问题的波数和圆频率;K为其刚度矩阵,写成:K=∫ΩBTBdΩ;其中B为车内声场预测模型的声压梯度矩阵,可以写成:
M为汽车车内声场模型的质量矩阵:Μ=∫ΩΨTΨdΩ;C为汽车车内声场的声学阻尼矩阵:
F为汽车车内声场所受的激励向量:P为车内声场边界上的节点激励向量:P={p1,p2,…,pn}T。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510005758.9A CN104504215A (zh) | 2015-01-07 | 2015-01-07 | 基于单位分解“有限元-无网络”单元的汽车车内声场预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510005758.9A CN104504215A (zh) | 2015-01-07 | 2015-01-07 | 基于单位分解“有限元-无网络”单元的汽车车内声场预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104504215A true CN104504215A (zh) | 2015-04-08 |
Family
ID=52945612
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510005758.9A Pending CN104504215A (zh) | 2015-01-07 | 2015-01-07 | 基于单位分解“有限元-无网络”单元的汽车车内声场预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104504215A (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104850721A (zh) * | 2015-06-03 | 2015-08-19 | 湖南大学 | 一种基于概率和区间的混合进行外声场预测的方法和装置 |
CN105279327A (zh) * | 2015-10-16 | 2016-01-27 | 奇瑞汽车股份有限公司 | 一种汽车车内振动和噪声的计算机仿真预测方法 |
CN105701286A (zh) * | 2016-01-12 | 2016-06-22 | 北京航空航天大学 | 一种预测汽车内部模糊不确定声压的数值计算方法 |
CN105760585A (zh) * | 2016-02-03 | 2016-07-13 | 北京航空航天大学 | 一种基于高阶摄动理论的模糊声压快速求解方法 |
CN106021847A (zh) * | 2016-04-28 | 2016-10-12 | 北京航空航天大学 | 一种基于高维模型表征与顶点分析策略的空腔噪声预测方法 |
CN107045573A (zh) * | 2017-04-05 | 2017-08-15 | 宁波大学 | 一种基于模糊逐维分析理论的车内振动噪声预测方法 |
CN107860469A (zh) * | 2017-11-22 | 2018-03-30 | 重庆大学 | 一种基于正交多项式拟合的变电站噪声预测方法 |
CN108200524A (zh) * | 2016-12-08 | 2018-06-22 | 浙江吉利控股集团有限公司 | 一种开口声腔模态频率的测试方法及系统 |
CN110164466A (zh) * | 2019-04-28 | 2019-08-23 | 清华大学苏州汽车研究院(相城) | 一种应用于汽车发动机主动噪声控制的车内声场可视化方法 |
CN110287570A (zh) * | 2019-06-18 | 2019-09-27 | 中国船舶工业集团公司第七0八研究所 | 一种船舶艏侧推器激励载荷转换方法 |
CN111814262A (zh) * | 2020-09-02 | 2020-10-23 | 广州汽车集团股份有限公司 | 车辆内低频声场的仿真分析方法、平台及计算机存储介质 |
CN112632818A (zh) * | 2020-12-17 | 2021-04-09 | 重庆交通大学 | 高阶梯度光滑无网格方法、系统、介质、计算机设备 |
CN113132847A (zh) * | 2021-04-13 | 2021-07-16 | 北京安声科技有限公司 | 主动降噪耳机的降噪参数确定方法及装置、主动降噪方法 |
WO2023143110A1 (zh) * | 2022-01-29 | 2023-08-03 | 华为技术有限公司 | 一种偏微分方程求解方法及其相关设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1815229A (zh) * | 2005-02-02 | 2006-08-09 | 南京汽车集团有限公司 | 一种预测汽车蒙皮表面阻尼处理结构声辐射的方法 |
CN102999678A (zh) * | 2012-12-26 | 2013-03-27 | 北京航空航天大学 | 一种基于非线性多目标区间鲁棒优化的汽车降噪方法 |
US20140019103A1 (en) * | 2012-07-11 | 2014-01-16 | Sumitomo Rubber Industries, Ltd. | Method for estimating noise performance of rolling tire |
-
2015
- 2015-01-07 CN CN201510005758.9A patent/CN104504215A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1815229A (zh) * | 2005-02-02 | 2006-08-09 | 南京汽车集团有限公司 | 一种预测汽车蒙皮表面阻尼处理结构声辐射的方法 |
US20140019103A1 (en) * | 2012-07-11 | 2014-01-16 | Sumitomo Rubber Industries, Ltd. | Method for estimating noise performance of rolling tire |
CN102999678A (zh) * | 2012-12-26 | 2013-03-27 | 北京航空航天大学 | 一种基于非线性多目标区间鲁棒优化的汽车降噪方法 |
Non-Patent Citations (2)
Title |
---|
姚凌云: "基于分区光滑理论与无网格法的声学数值方法研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
王青春 等: "《有限元方法及Midas软件在汽车结构分析中的应用》", 31 May 2014, 机械工业出版社 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104850721A (zh) * | 2015-06-03 | 2015-08-19 | 湖南大学 | 一种基于概率和区间的混合进行外声场预测的方法和装置 |
CN105279327A (zh) * | 2015-10-16 | 2016-01-27 | 奇瑞汽车股份有限公司 | 一种汽车车内振动和噪声的计算机仿真预测方法 |
CN105701286B (zh) * | 2016-01-12 | 2019-05-03 | 北京航空航天大学 | 一种预测汽车内部模糊不确定声压的数值计算方法 |
CN105701286A (zh) * | 2016-01-12 | 2016-06-22 | 北京航空航天大学 | 一种预测汽车内部模糊不确定声压的数值计算方法 |
CN105760585A (zh) * | 2016-02-03 | 2016-07-13 | 北京航空航天大学 | 一种基于高阶摄动理论的模糊声压快速求解方法 |
CN106021847A (zh) * | 2016-04-28 | 2016-10-12 | 北京航空航天大学 | 一种基于高维模型表征与顶点分析策略的空腔噪声预测方法 |
CN106021847B (zh) * | 2016-04-28 | 2018-10-02 | 北京航空航天大学 | 一种基于高维模型表征与顶点分析策略的空腔噪声预测方法 |
CN108200524B (zh) * | 2016-12-08 | 2020-08-18 | 浙江吉利控股集团有限公司 | 一种开口声腔模态频率的测试方法及系统 |
CN108200524A (zh) * | 2016-12-08 | 2018-06-22 | 浙江吉利控股集团有限公司 | 一种开口声腔模态频率的测试方法及系统 |
CN107045573B (zh) * | 2017-04-05 | 2019-07-30 | 宁波大学 | 一种基于模糊逐维分析理论的车内振动噪声预测方法 |
CN107045573A (zh) * | 2017-04-05 | 2017-08-15 | 宁波大学 | 一种基于模糊逐维分析理论的车内振动噪声预测方法 |
CN107860469A (zh) * | 2017-11-22 | 2018-03-30 | 重庆大学 | 一种基于正交多项式拟合的变电站噪声预测方法 |
CN110164466A (zh) * | 2019-04-28 | 2019-08-23 | 清华大学苏州汽车研究院(相城) | 一种应用于汽车发动机主动噪声控制的车内声场可视化方法 |
CN110287570A (zh) * | 2019-06-18 | 2019-09-27 | 中国船舶工业集团公司第七0八研究所 | 一种船舶艏侧推器激励载荷转换方法 |
CN110287570B (zh) * | 2019-06-18 | 2022-11-11 | 中国船舶工业集团公司第七0八研究所 | 一种船舶艏侧推器激励载荷转换方法 |
CN111814262A (zh) * | 2020-09-02 | 2020-10-23 | 广州汽车集团股份有限公司 | 车辆内低频声场的仿真分析方法、平台及计算机存储介质 |
CN112632818A (zh) * | 2020-12-17 | 2021-04-09 | 重庆交通大学 | 高阶梯度光滑无网格方法、系统、介质、计算机设备 |
CN112632818B (zh) * | 2020-12-17 | 2023-03-14 | 重庆交通大学 | 高阶梯度光滑无网格方法、系统、介质、计算机设备 |
CN113132847A (zh) * | 2021-04-13 | 2021-07-16 | 北京安声科技有限公司 | 主动降噪耳机的降噪参数确定方法及装置、主动降噪方法 |
CN113132847B (zh) * | 2021-04-13 | 2024-05-10 | 北京安声科技有限公司 | 主动降噪耳机的降噪参数确定方法及装置、主动降噪方法 |
WO2023143110A1 (zh) * | 2022-01-29 | 2023-08-03 | 华为技术有限公司 | 一种偏微分方程求解方法及其相关设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104504215A (zh) | 基于单位分解“有限元-无网络”单元的汽车车内声场预测方法 | |
Upadhyay et al. | Numerical analysis perspective in structural shape optimization: A review post 2000 | |
Zheng et al. | The MLS-based numerical manifold method with applications to crack analysis | |
Liu et al. | An automatic approach for the acoustic analysis of three-dimensional bounded and unbounded domains by scaled boundary finite element method | |
Song | The scaled boundary finite element method in structural dynamics | |
Cui et al. | An accurate and efficient scheme for acoustic-structure interaction problems based on unstructured mesh | |
Meldi et al. | An arbitrary Lagrangian–Eulerian approach for the simulation of immersed moving solids with Lattice Boltzmann Method | |
He et al. | A non-contact acoustic pressure-based method for load identification in acoustic-structural interaction system with non-probabilistic uncertainty | |
CN102999678A (zh) | 一种基于非线性多目标区间鲁棒优化的汽车降噪方法 | |
CN103559366A (zh) | 基于光滑有限元边界元法的车身结构声场耦合预测方法 | |
Doke et al. | Vehicle concept modeling: A new technology for structures weight reduction | |
He et al. | Acoustic simulation using α-FEM with a general approach for reducing dispersion error | |
CN104317985A (zh) | 一种基于界带有限元和拉格朗日坐标的流体仿真方法 | |
He et al. | Numerical investigation of ES-FEM with various mass re-distribution for acoustic problems | |
He et al. | Mid-frequency acoustic analysis using edge-based smoothed tetrahedron radialpoint interpolation methods | |
Jiang et al. | Acoustic characterization of ducts lined with poroelastic materials based on wave finite element method | |
CN100555294C (zh) | 主梁断面气动自激力的全过程数值模拟方法 | |
Piggott et al. | Unstructured adaptive meshes for ocean modeling | |
Wu et al. | The development of hybrid ES-FE-SEA method for mid-frequency vibration analysis of complex built-up structure | |
Liu et al. | Wave based method for flexural vibration of thin plate with general elastically restrained edges | |
Puri et al. | A Krylov–Arnoldi reduced order modelling framework for efficient, fully coupled, structural–acoustic optimization | |
CN111177958B (zh) | 一种微穿孔板声衬声学性能数值分析方法 | |
JP7542820B2 (ja) | 発泡樹脂製吸音材の内部構造を設計する装置、並びに、発泡樹脂製吸音材の音響特性を算出する装置 | |
Gagliardini et al. | Virtual SEA-FEA-based modeling of mid-frequency structure-borne noise | |
Chazot et al. | Harmonic response computation of poroelastic multilayered structures using ZPST shell elements |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150408 |
|
RJ01 | Rejection of invention patent application after publication |