狭长封闭声场的近似圆锥束跟踪法模拟研究

狭长空间声传播特性研究对于地铁，长廊等场所的音质设计是非常有意义的，本文提出了模拟狭长封闭声场的近似圆锥束跟踪法，利用它研究了矩形和非矩形狭长声场的空间分布特性。仿真结果表明，对于狭长声场，经典扩散场理论并不适用，而利用本文的数值模拟方法能正确地分析这类声场的特性。     

矩形表面波探头声场的高斯声束叠加法

利用矩形压电晶片和有机玻璃楔块折射可激励出超声表面波,广泛用于固体近表面缺陷检测和材料特性测量.由于描述表面波三维声场的理论方法还鲜有报道,因而主要采用简化的表面波二维声场模型来定量分析这类问题.高斯声束模型近些年被广泛应用于解决超声体波传播的各种复杂问题,然而,目前还没有将其扩展应用到超声表面波的声场的计算中.通过结合表面波格林方程和矩形换能器的高斯声束模型,推导出基于高斯声束叠加的表面波三维声场解析解.进一步,将该方法与点源叠加的数值解进行了分析比较,计算结果表明表面波声场的高斯声束叠加方法在具有较好计算精度的同时,还具有更快的计算效率.     

基于多元高斯声束模型的圆柱体超声检测声场仿真

为了提高圆柱体缺陷检测的准确性,需了解其超声声场分布。针对圆柱体实际检测的要求,本文采用超声水浸检测法。通过多元高斯声束模型分别模拟了实心圆柱体及空心圆柱体中的辐射声场,计算仿真了垂直入射与倾斜入射至圆柱体时的声场,实现了检测声场的可视化。通过计算得出仿真结构与实际理论相符,表明本论文的方法能够应用在圆柱体结构超声声场仿真中。     

束晕-混沌控制中的粒子跟踪模拟研究

运用PIC程序研究了强流离子束中粒子的横向运动，发现了不加控制时束晕粒子并非一直处于晕区和施加非线性控制后粒子的横向运动被限制在一定的范围内形成环状，以及均方根半径的变化近似成为周期运动等性质.根据观察到的均方根半径及其变化率的规律性，提出了一种新的自适应控制器.用该控制器不仅能在很短的时间内完全控制住束晕，而且不需要对增益因子进行精确计算，也能在系统参数改变的情况下取得较好的控制效果. 关键词： 束晕_混沌 粒子跟踪 数值模拟 自适应控制     

适用于任意形状封闭空间的声场预测相干模型

针对现有几何声学的方法对封闭空间内声场进行预测时在中低频段出现较大误差的问题,该文提出一种近似圆锥声束追踪法和相干反射场理论相结合的声场预测新模型。在近似圆锥声束追踪法基础上,考虑声束轴线在边界多次反射时声压和相位的改变,最后计算不同声波之间的干涉效应,建立一种适用于任意形状封闭空间的声场预测相干模型。利用该模型对某一矩形封闭空间进行声场预测,通过对边界元法、Raynoise软件相干和非相干算法的预测结果和本模型的数值模拟结果对比。结果表明,文中提出的方法和边界元法的计算结果在中低频段非常吻合,两者的计算结果平均绝对误差为1.1 d B。本模型在中低频率下与同样考虑了相位的Raynoise相干算法相比有更好的准确性,在较高频率上,本模型计算结果与Raynoise相干算法计算结果非常吻合。     

脉冲响应法测量混响时间技术的研究

本文介绍了一种根据厅堂脉冲响应特性,利用计算机数字信号处理技术,迅速,准确地测量混响时间的实现方法。     

小阻尼界面房间声传输函数和声脉冲响应函数的有限元素法

有限元法可用于以声波动方程为基础通过数值计算求解室内声场,适用于分析界面阻抗非均匀分布和复杂形状房间内声场的低频特性。本文首先介绍了小阻尼界面条件下室内声场简正方式、声衰变系数、混响时间的ＦＥＭ计算方法。在此基础上导出了房间内两点之间声传输函数和声脉冲响应函数的ＦＥＭ计算模型,并以矩形房间为例详细讨论了有关细节。本文所讨论的计算模型可以映房间内不同的声源点,接收点位置上的声压频谱特性和脉冲响应的时     

基于改进声线跟踪法的室内交通噪声动态模拟

对碰撞点的有效性判断是传统声线跟踪法的一个重要步骤,在计算复杂的室内空间问题时计算量较大,为解决这一问题采用了一种与空间剖分相结合的改进的声线跟踪法。将该方法与微观交通流仿真和车辆噪声排放模型进行结合,实现了道路交通噪声透过窗户在多连通室内空间传播的动态模拟。最后采用该方法对相同道路和交通流条件下不同建筑朝向的4种布局室内噪声进行动态模拟。分析了室内交通噪声的大小和分布与建筑物布局、窗户的朝向、窗户的形式和面积等因数之间的关系。结果表明:窗户正对道路的房间比窗户侧对道路的房间噪声高4～6 d B(A),房间内等效声级与窗户面积的对数成正比关系。     

一种结合室内空间剖分的声线跟踪法

声线跟踪法是应用于室内声场预测的几何声学模拟经典方法之一,但该方法的计算效率随着室内空间布局复杂度的增加而降低。为了高效地应用于各种室内空间布局,提出了一种结合室内空间剖分的声线跟踪法,该方法将整个室内空间剖分成有限个凸多面体子空间,以减少在声线跟踪过程中与声线求交点的墙面数量,并省略了碰撞点有效性判断,显著提高了计算效率,尤其适用于大型复杂的室内空间布局。     

螺旋型电磁超声辐射声场对缺陷检测影响的 计算机模拟*

电磁超声(EMAT)作为一种非接触超声检测技术,已得到了普遍的关注和研究。然而,EMAT目前只得到了一些有限的应用,尤其在缺陷探伤应用方面还不能替代压电超声,深入探究其原因并加以改进,对推动EMAT在缺陷检测方面的应用具有重要意义。利用有限元方法建立了螺旋型EMAT模型,分析了其辐射声场和声场特征;通过将有限元计算所得到的洛伦兹力分布简化为门函数,并结合格林函数法,得到了螺旋型EMAT辐射声场的简化解析结果。这些方法可用于计算螺旋型EMAT辐射声场指向性以及在缺陷检测中比较关心的一些特征参数,如扩散角、偏离角等参数。通过螺旋型EMAT辐射声场实验,验证了理论分析结果。通过与压电超声检测理论与实验的对比,研究认为：螺旋型EMAT辐射声场的具有横波和纵波等多种波模及多种模态转换的特点,以及相对复杂的异型声场几何分布等复杂特性,使其在缺陷检测应用方面受到局限。这些研究方法对改进EMAT设计,优化控制EMAT辐射声场,拓展EMAT在缺陷检测方面的应用具有参考价值。     

Measurement of sound insulation of acoustic louvres by an impulse method

This paper presents results of a study of the sound attenuation of acoustic louvres. At the core of the study is an alternative method of measuring sound insulation, impulse response analysis, which circumvents the limitations imposed by standard and proposed standard methods. Using the impulse method, the sound transmission coefficient is measured at different angles of incidence and the angular dependency of transmission loss obtained. In the low frequency range, the transmission is governed by a mass layer effect. The value of transmission loss is independent of angle of incidence. For the mid and high frequencies, diffraction, interference and absorption determine louvre performance and an angular dependency is observed. The transmission at the angle of incidence, corresponding to a line-of-sight through the louvre blades, is the dominant contribution to the angle average value and a single measurement at the pitch of the louvre approximates the overall transmission loss. For the case considered, the geometry of the blades has little influence on the transmission at low frequencies and the mass of the blades has little influence at higher frequencies. In a companion paper, the impulse data are used to predict the insertion loss provided by the louvre when installed in a plant room.     

An approximate treatment of reflection coefficient in the phased beam tracing method for the simulation of enclosed sound fields at medium frequencies

The phased beam tracing method (PBTM) was suggested as a medium-frequency simulation technique for the calculation of impulse response, although main assumptions of geometric acoustics still hold. The phased method needs the reflection coefficient for characterizing the acoustic property of a surface and the complex wave number for describing the propagation characteristics. In this study, two types of approximate real reflection coefficients derived from the measured absorption coefficient were tested for a practical applicability. As a test example, pressure impulse responses and energy impulse responses computed from the PBTM were compared with those from the measurement and the ordinary beam tracing method. The PBTM employing the approximate reflection coefficients greatly increased the accuracy of the prediction compared to the ordinary beam tracing method, in particular at the medium-frequency range in octave bands above the Schroeder cutoff frequency. A comparison was made between angle-dependent and angle-independent reflection coefficients in the calculation of acoustic measures. Although the angle-dependent reflection coefficient yielded best matched results with measured data, but the simple angle-independent reflection coefficient can be also used with a reasonably good precision.     

A note on the proper frequency resolution for the room transfer function in the phased beam tracing method

In the simulation of acoustic response of a room, one cannot satisfactorily estimate the transfer function or the impulse response with a sparse frequency resolution. We have studied about the frequency resolution for a satisfactory estimation of room transfer function or the reverberation time using the phased beam tracing method. In order to investigate what happens with the frequency resolution, a number of frequency resolution (0.01-1 Hz), which are usually employed in the numerical analysis, was tested for a room model. It was found that the total accumulated phase in a transfer function and the late part of an impulse response of a room are influenced by the frequency resolution. Main reason for the difference is the correct detection of non-minimum phase zeros, depending on the frequency resolution. A criterion for a proper frequency resolution was suggested by considering the modal density of three-dimensional space. When a sparse dataset was initially given, we showed that cubic spline interpolation can be used to enhance the precision of detection for non-minimum phase zeros.     

激光加速电子束放射照相的模拟研究

激光加速产生高能量电子束具有源尺寸小、准单能、脉宽窄等特征.通过蒙特卡罗程序模拟研究了高能电子束的放射照相.模拟了200 MeV准直电子束照射台阶靶、厚铁靶,11 MeV点源电子束照射惯性约束聚变模型靶,以及70 MeV点源电子束在激光等离子体磁场下的偏转.结果表明激光加速电子束在探伤厚材料内部、确认薄材料界面、测量电磁场等诊断中具有高时空分辨、灵敏等能力.     

A hybrid model for the acoustic response of a two-dimensional rough surface to an impulse incident from a refracting medium

A numerical model to calculate the impulse response of a two-dimensional, impenetrable, rough surface directly in the time domain has been recently introduced [J. Acoust. Soc. Am. (2000) 107, 27]. This model is based on wedge diffraction theory and assumes that the half-space containing the source and receiver is homogeneous. In this work, the model is extended to handle media where the index of refraction varies with the distance to the surface by merging the scattering model with a ray-based propagation model. The resulting hybrid model is tested against a Finite-Difference Time-Domain (FDTD) method for the problem of backscattering from a corrugated surface in the presence of a refractive layer. This new model can be applied, for example, to calculate acoustic reverberation from the sea surface in cases where the water mass is inhomogeneous and in the presence of a subsurface bubble layer at low frequencies where dispersion is negligible. It can also be used for atmospheric propagation problems where there is a sound speed gradient overlying rough terrain.     

Normal tracing deflectometry using a secondary light source

Scanning deflectometric profilers based on an f–gθ system are typical optical tools used to measure mirror profiles at many synchrotron facilities. Unlike these profilers, which are based on a pencil beam, here a secondary light source and a pinhole are used to construct a system that automatically selects a beam that will always pass through the pinhole and propagate along the normal direction of the measured area on the surface under test. By measuring the angle variation of the selected beam, slope variations of the surface under test can be measured. Systematic errors introduced by manufacturing defects or aberrations of an optical element, which greatly degrade the performance of traditional profilers, could be minimized by using the developed method. Simulation values of the proposed method and a conventional method are compared.     

The impulse response function of apertured and misaligned ABCD optical system in phase-space

By expanding the aperture into complex Gaussian functions, the approximate analytical impulse response function of apertured and misaligned ABCD optical system in phase-space is determined, which allows us to formulate the input-output relation completely in terms of Wigner distribution functions. Its shown that the hard-aperture placed before the optical system will change the shape of output signal Wigner distribution and the misalignment of optical system will introduce the output Wigner distribution function a coordinate shift in phase-space. For examining the analytical results some numerical simulation are done and the graphical results as well as the absolute errors between analytical results and numerical integral ones are presented.     

Two-photon polymerisation of novel shapes using a conically diffracted femtosecond laser beam

Sub-micron ring, pillar and wall structures were written by two-photon polymerisation of a sol-gel resin using a femtosecond laser beam which was shaped using internal conical diffraction. The ring structure was written using a demagnified image of the ring-shaped beam which arises in conical diffraction of a narrow light beam. The pillar and wall structures were produced by imaging the Bessel beam formed by conical diffraction in combination with a converging lens.