利用边界元法计算无界声场中结构体声辐射

建立了无界声场中结构体声辐射的边界元数学模型 .计算中以八节点曲边四边形等参元模拟结构体表面 .利用复合亥姆霍兹积分方程解决了边界元方法在计算声学外问题时解的不唯一现象 .同时利用正则化关系式将复合亥姆霍兹积分方程中的超奇异数值积分转化为弱奇异数值积分 .最后以脉动球和振荡球声辐射为例 ,验证了数值计算结果 ,表明利用边界元方法计算声学外问题时必须考虑解的不唯一问题     

含局部预应力的圆柱壳结构声辐射特性分析

以圆柱壳结构-声辐射特性为研究对象,分析了局部预应力对结构动力学特性的影响.基于Flügge圆柱壳薄壳理论,建立了含局部预应力的圆柱壳结构应力-应变方程,应用变分原理推导了含局部预应力的圆柱壳结构动力响应方程.结合结构有限元和声边界元法,对圆柱壳结构的声辐射问题进行数值分析,对比了预应力对声辐射功率和声指向性的影响.结果表明,由于局部预应力的存在,使得结构的局部和总体刚度都发生变化,并引起结构声振特性的变化.     

基于有限元/边界元的双层周期加筋板声辐射分析

采用有限元和边界元相结合的方法研究了水中双层周期加筋板结构在简谐力作用下的声辐射特性。给出了有限元方程,详细推导了边界元方程,利用有限元和间接边界元的声振耦合方程,计算了在谐力作用下结构表面振动的辐射声场,从而为揭示双层周期加筋板结构声学特性提供了一种数值方法。该数值方法可用于任意复杂结构在流体介质中的振动和声学分析。     

基于径向积分边界单元法的内场声学分析

由于Helmholtz方程的基本解是频率的函数,传统边界单元法在处理声场特征值问题时具有天生的缺陷.采用Laplace方程基本解生成积分方程,通过径向积分法将在此过程中产生的域积分项转化为边界积分.此方法克服了传统边界单元法系数矩阵对频率的依赖,同时克服了特解积分法和双重互易法对特解的依赖,将内场声学特征值问题转化为广义特征值问题.最后通过内场声辐射分析和声学特征值分析验证了算法的有效性.     

改进计算对称壳体声辐射的边界积分方程法

用边界积分方程法计算轴对称振动表面的声辐射时,必须处理好特征频率下表面Helmholtz方程无唯一解的问题和奇点附近区域上的奇异积分问题。本文把表面Helmholtz方程与关于内点的补充方程联立组成线性方程组,用最小二乘法求解,并采用极坐标变换将奇异积分转换成普通积分,从而可以方便地计算任意形状轴对称体在各个频率下的声辐射。     

辐射-导热耦合换热边界元算法

本文提出了一种求解辐射-导热耦合换热问题的边界单元算法(BEM),该方法将两种传热方式通过辐射热源耦合起来.首先,采用BEM对辐射传热方程、辐射热源方程和含有辐射热源的热传导方程进行离散;其次,利用辐射传热方程消除辐射热源方程中的辐射热流项;然后,根据Stefan-Boltzmann定律形成含有温度四次方以及热流密度表示的非线性代数方程组.出现在所有积分方程中的域积分由径向积分法转换成边界积分,形成了对于参与性介质问题也只需在边界上划分单元的纯边界元算法.最后,用Newton-Raphson迭代法对方程组进行求解.提供的数值算例将表明本文所介绍方法的正确性与有效性.     

基于声压修正边界元方法的气动噪声数值模拟

为考虑非紧致固体边界散射作用的影响,通过声模拟与边界元相结合的声压修正方法预测流体流动引起的噪声,将所分析的问题降低一维,能够满足无限远处边界条件。频率不是很高时,四极子源辐射的声强度远小于边界散射声强度,可以忽略耗时的四极子源项体积积分计算。通过傅立叶变换在频域求解声压方程,避免了繁琐的时间延迟计算。     

离心风机蜗壳振动声辐射的定量预测

针对离心风机蜗壳振动产生的声辐射,提出了一种基于边界元分析的声辐射定量预测方法:根据蜗壳声辐射的边界元模型布置振动测点,测量其法向振动的加速度谱;提取加速度谱中峰值较高的基频和旋转频率分量,并将其转换成速度谱后加栽于数值模型的对应节点上,作为声辐射问题的速度边界条件;利用边界元法计算声场及声功率,预测蜗壳振动产生的噪声.结果表明:结构振动声比气动声小很多,相对于风机进出口的气动噪声而言,蜗壳因振动辐射的噪声可以忽略不计;当风机进出口安装消声器后,蜗壳振动的辐射噪声可能成为主要噪声源.基于边界元分析的蜗壳振动辐射噪声定量预测方法也可以应用到其他壳体振动辐射噪声的预测中.     

基于径向积分边界单元法的内场声学分析

由于Helmholtz方程的基本解是频率的函数,传统边界单元法在处理声场特征值问题时具有天生的缺陷。采用Laplace方程基本解生成积分方程,通过径向积分法将在此过程中产生的域积分项转化为边界积分。此方法克服了传统边界单元法系数矩阵对频率的依赖,同时克服了特解积分法和双重互易法对特解的依赖,将内场声学特征值问题转化为广义特征值问题。最后通过内场声辐射分析和声学特征值分析验证了算法的有效性。     

求解声辐射特性的Burton-Miller改进型积分方程

提出了一种求解结构声辐射问题的Burton-Miller改进型边界积分方程,利用拉普拉斯方程的特性对传统边界积分方程及其法向偏导方程进行处理,转化其中与频率相关的高阶奇异积分项和柯西型积分项分别为弱奇异积分项和不含奇异性的积分项;进一步联立求解结构内外拉普拉斯问题下的边界积分方程,将与频率无关的高阶奇异积分项和柯西型积分项转化为弱奇异积分乘积的形式,以保证计算的精度.以脉动球源和横向振动球源为例,将所得结果与传统边界积分方程相比较,表明该方法不仅可以保证全波数范围内解的唯一性,且具有很高的计算精度.     

椭球外区域问题的自然边界元法

主要讨论了椭球面外部区域Laplace方程的自然边界元法.首先引入椭球坐标,通过分离变量法导出了Poisson积分公式和自然积分算子的无穷级数显示表达式.这样原无界区域问题就归化为椭球面上边界积分方程,然后再数值求解该积分方程.给出了该积分方程的变分问题的适定性和逼近解的误差估计,且该误差估计不仅依赖于网格参数而且依赖与级数截断后的项数.数值例子说明了该方法的有效性和理论分析的正确性.     

准格林函数方法在弹性扭转问题中的应用

利用Poisson方程的基本解构造一个准格林函数,这个函数满足Poisson方程的齐次边界条件．应用格林函数将边值问题化为积分方程,并通过建立一个规范化的边界方程来表示问题的边界,以克服积分方程核的奇异性．弹性扭转问题可看成是Poisson方程的边值问题,尺一函数理论保证了对于任何复杂的区域,总可以找到一个规范化方程,从而可以将弹性扭转问题化为一个无奇异性的第二类Fredholm积分方程．数值算例表明,该方法具有较高的精度,可用于力学、物理中复杂边值问题的研究。     

高速轿车二维风沙流边界层数学模型

基于普朗特的边界层思想,以不可简化成稀相的风沙流体为研究对象,以风沙流的一般方程为基础,在基本假设的条件下,建立了高速轿车车身顶部过流表面的边界层微分方程并给出沙流在近壁表面以滑移为特征的边界条件。又考虑到固体颗粒对边界层的扰动作用,引入了扰动因子。通过对边界层厚度的积分并利用相关条件,得到了边界层动量积分方程。     

浅海声波散射边值问题的小波数值解

浅海声波散射边值问题的小波数值解张伯坚广东农工商管理学院基础部,５１０５０７,广州关键词散射,Ｈｅｌｍｈｏｌｔｚ方程,周期小波分类号（中图）Ｏ１７５．２,Ｏ２４１．８２;（１９９１ＭＲ）３４Ｂ１２,４１Ｃ１５本文把浅海声波散射边值问题归结为边界积分方...     

振动噪声和气动噪声统一分析方法的理论研究

首先对振动声学和气动声学理论进行了关联性分析,证明了振动声辐射的积分方程是FW-H方程的一种特殊形式,表明声比拟理论不仅适用于气动噪声的预测,而且同样适用于振动噪声的预测,是分析噪声产生机理的普适化理论。然后利用FW-H方程计算了球脉动和振动辐射噪声的声压级,计算结果和采用振动声学理论分析得到的结果完全一致。在此基础上,对流固耦合噪声的预测方法进行了探讨分析。结果表明,流体和弹性薄壳体耦合作用激发的噪声仅取决于流体的运动速度和流固界面的作用力,而与壳体边界的振动速度无直接相关。     

三维Helmholtz方程的边界点方法

将移动最小二乘近似和边界积分方程相结合,提出了求解三维Helmholtz方程内外边值问题的无网格边界点方法.该方法用单层位势理论将Helmholtz方程转化为间接边界积分方程,并用边界点法离散间接边界积分方程.由于边界积分方程中含有基本解的积分计算时会出现弱奇异,详细推导了弱奇异积分的计算方式.数值算例表明了间接边界点法求解三维Helmholtz方程的有效性.     

瞬态热传导问题的精细积分边界元法

对于三维瞬态热传导问题,在考虑内部热源的情况下,采用双重互易边界元法(DRBEM)结合精细积分法(PIM)进行求解。该方法根据含有内部热源的各向同性介质瞬态常系数热传导问题的控制方程,通过加权余量法推导出相应的边界积分方程,然后用双互易法(DRM)处理得到的边界积分方程,将热源项和温度关于时间导数项引起的域积分通过径向基函数(RBF)逼近后转化为边界积分。之后将边界积分方程离散,得到与时间相关的一阶常系数微分方程组,最后,在获得解析解的过程中,通过PIM处理其中的矩阵指数函数(MEF)。通过三个数值算例来验证该方法的准确性和稳定性。     

三维定常对流扩散方程的双方程方法

三维定常对流扩散方程的经典边界积分方程,其类型关于未知对流扩散势导数是第一类积分方程,关于未知对流扩散势是第二类积分方程。本文从格林公式出发,通过建立位势的单、双场守恒积分公式,推导出三维定常对流扩散方程新的边界积分方程,其类型与经典方程相反。对不同的边界采用不同的方程,由此把双方程边界元方法推广到三维空间。