快速多极边界元法用于扩散光学断层成像研究

在求解扩散光学断层成像中的正向问题时, 目前普遍采用有限元法, 但是随着实际模型规模的增大, 有限元法的计算量问题日益显著, 而边界元法则由于可以降低计算维度使计算量减少而备受关注. 本文以均匀的高散射介质为模型, 研究了将快速多极边界元法用于扩散光学断层成像的正向问题. 快速多极边界元法利用核函数的多极展开, 将常规边界元法中系数矩阵和迭代矢量的乘积项等价为相应四叉树结构的一次递归, 再结合广义最小残量法进行迭代求解. 将计算结果和蒙特卡罗法的模拟结果进行了比较, 表明利用快速多极边界元法的模拟结果和蒙特卡罗法的结果有很好的一致性. 研究结果验证了快速多极边界元法可以用于扩散光学断层成像, 为其大规模和实时成像带来可观的前景. 关键词： 扩散光学断层成像 边界元法 快速多极边界元法     

新的对角形式快速多极边界元法求解声学Helmholtz方程

传统外部声学Helmholtz边界积分方程无法在个人计算机上求解大规模工程问题. 为了有效解决这个问题, 将快速多极方法引入到边界积分方程中, 加速系统矩阵方程组的迭代求解. 由于在边界积分方程中引入基本解的对角形式多极扩展, 新的快速多极边界元法的计算效率与传统边界元相比显著提高, 计算量和存储量减少到O(N)量级(N为问题的自由度数). 包括含有420000个自由度的大型潜艇模型数值算例验证了快速多极边界元法的准确性和高效性, 清楚表明新算法在求解大规模声学问题中的优势, 具有良好的工程应用前景.     

新的对角形式快速多极边界元法求解声学Helmholtz方程

传统外部声学Helmholtz边界积分方程无法在个人计算机上求解大规模工程问题. 为了有效解决这个问题, 将快速多极方法引入到边界积分方程中, 加速系统矩阵方程组的迭代求解. 由于在边界积分方程中引入基本解的对角形式多极扩展, 新的快速多极边界元法的计算效率与传统边界元相比显著提高, 计算量和存储量减少到O(N)量级(N为问题的自由度数). 包括含有420000个自由度的大型潜艇模型数值算例验证了快速多极边界元法的准确性和高效性, 清楚表明新算法在求解大规模声学问题中的优势,     

边界元奇异与近奇异数值积分方法及其应用于大规模声学问题

提出了综合处理Burton-Miller方法所导致的奇异积分与近奇异积分问题的数值求积方法,以此改进了基于常量元素的常规边界元和低频快速多极边界元方法。对于奇异积分问题,利用Hadamard有限积分方法进行解决;对于近奇异积分问题,则采用极坐标变换法和PART方法(Projection and Angular&;Radial Transformation)进行克服。与解析解和LMS Virtual.Lab商业软件的结果比较验证了方法的正确性,并对比分析了奇异积分与近奇异积分对计算精度的影响。采用低频快速多极子方法以加速常规边界元法的计算效率,计算分析了计算复杂度,并成功实现了34万自由度大规模问题的计算。结果表明,近奇异积分问题主要由超奇异核函数引起,对计算精度的影响不容忽略;快速多极边界元法的精度与常规边界元法一致,但计算复杂度要远低于后者。      

新型快速多极边界元法求解电荷任意分布的二维静电场

在初始快速多极边界元法(FMM)基础上提出一种适合位势问题的新型快速多极边界元格式,并用于求解静电场问题.新型算法引入对角化概念,减少了形成局部展开系数的时间,提高计算效率.最后给出数值算例,证明了新型算法的计算精度及处理大规模问题的速度优势.     

三维声学多层快速多极子边界元及其应用

快速多极子边界元算法可以加速矩阵和向量乘法运算, 将传统边界元算法的计算量和内存占用量分别降为O(N log²N)和O(N), 适用于大型声学模型模拟计算. 本文发展了一种基于Burton-Miller方程的三维多层声学快速多极子边界元算法. 将新的自适应树状算法应用到对角形式的快速多极子边界元算法, 并使用最新提出的解析式源点矩计算公式, 进一步提高了快速多极子边界元的计算效率. 绝对软球体在内部共振频率处的散射声场计算, 验证了所发展算法在共振频率处求解的正确性. 与Bapat所提供的程序在多脉动球体辐射声场计算精度的比较, 验证了算法及程序在大型模型声学计算中的准确性, 同时显示了其求解的高效性. 最后, 将该算法用于车内声场及水下声学探测的分析计算.     

多层快速多极子算法中的两步插值技术

多层快速多极子算法(MLFMA)在快速多极子算法(FMM)的基础上按多层聚集、层间转移和多层扩散的思路以达到优化矩阵向量积的运算的目的,其中多层聚集和多层扩散过程,随着层数递增,角谱积分采样点数逐层递增,为了快速计算角谱积分,需要采用插值技术和反插值技术以提高计算效率.应用两步插值技术替代传统的单步插值技术,大幅提高了多层快速多极子层间插值反插值操作的计算效率,对于应用普通个人计算机求解特大电大尺寸问题,具有重要意义.     

二维声学多层快速多极子边界元及其应用

与模型自由度的平方成正比的存储量和计算量,使传统边界元无法应用到大型模型的计算。为此,发展了一种二维声学多层快速多极子边界元算法。通过二维Helmholtz核函数展开理论的简要介绍,推导了源点矩计算、源点矩转移、源点矩至本地展开转移、本地展开转移公式,并详细描述了二维声学快速多极子边界元算法的具体实现步骤。使用快速傅里叶插值进行源点矩和本地展开系数的多层传递。采用对角左预处理方法,改善边界方程的条件数,减少迭代求解次数。最后通过数值算例,验证了所发展的二维声学快速多极子算法的正确性和高效性。      

大规模宽频弹性动力学分析的快速定向压缩边界元法

发展一种大规模宽频弹性动力学分析的快速定向压缩边界元法.证明弹性动力学核函数具有定向低秩特性，为采用快速定向压缩算法提供理论基础.根据S波的波数，将节点之间的相互作用划分为低频相互作用和高频相互作用，并将后者进一步划分为与多个楔形区.在楔形区上，可以采用核函数的定向低秩特性进行快速计算.低频相互作用与核无关快速多极边界元法中计算方法相同，不同方向楔形区上的变换矩阵可以采用坐标系旋转的方法进行快速计算.可对任意频率进行快速谐响应分析.数值算例表明：该方法可以将宽频弹性动力学问题计算复杂度降低到O（N log^αN）.与卷积求积边界元法相结合，也可以应用于弹性动力学瞬态分析.     

子波谱方法在三维声辐射和声散射中的应用

将边界变量用二维子波展开,获得了三维任意边界条件声辐射和声散射的边界积分方程的子波谱方法.采用以子波为权函数的Gauss积分法计算子波谱方法的系数,获得了与传统边界元法相同的计算量,克服了普通积分法计算子波系数计算量大甚至难以收敛的缺点;采用Duffy的方法解决了子波谱方法中的奇异积分,使其能够用普通的Gauss积分法计算.算例表明:子波谱方法系数矩阵压缩率超过50%以后,计算精度仍然高于传统边界元方法.     

消声器传递损失预测的边界元模态展开混合方法

将边界元法和解析方法结合形成一种混合方法用于计算消声器的传递损失,消声器被划分成若干个子结构,解析方法和边界元方法被分别用于计算规则结构和不规则结构的阻抗矩阵,不同子结构之间通过阻抗矩阵连接起来。为减少计算时间,采用一种基于模态配点法的简化方法。对单级膨胀腔、双级膨胀腔和穿孔管阻性消声器的传递损失进行了计算,混合方法计算结果与解析方法和三维数值方法计算结果吻合良好。分析了混合方法的计算效率并与传统子结构方法进行了比较,混合方法能明显节省计算时间。      

Application of dual reciprocity boundary element method for prediction of acoustic characteristics of ducts and silencers with three-dimensional potential flow

The Dual Reciprocity Boundary Element Method （DRBEM） is applied to predict the acoustic characteristics of ducts and silencers with three-dimensional potential flow, and the basic principle and numerical procedure of the proposed method are introduced. Compared to the Conventional Boundary Element Method （CBEM）, the DRBEM takes into account the second order terms of flow Mach number in the acoustic governing equation, which is suitable for the situations with higher Mach number subsonic flow. The four-pole parameters of a duct and a varying cross-sectional area expansion chamber are predicted with the DRBEM, and the predictions are compared with the one-dimensional analytical solutions and the CBEM results. The comparisons demonstrated that the present method is valid. Transmission loss of silencers with different structures was also calculated with the DRBEM. The results showed that the influence of the three-dimensional flow on the acoustic characteristics of silencers with complex structures is not negligible.     

Propagation of uncertainty in test-analysis correlation of substructured spacecraft

Organizations, such as the Air Force and NASA make critical decisions on spacecraft performance and survivability based on the results of test-analysis correlation metrics. In order to ensure the success of a new paradigm in finite element model validation where there is no system level test, uncertainty in the substructures must be propagated into the system level correlation metrics. The objective of this work is to quantify the level of accuracy required at the substructure level to produce acceptable analytical model accuracy at the system level. In preparation for future synthesized system level uncertainty analysis, a framework is presented for propagating analytical model uncertainty from a fixed interface Craig-Bampton substructure representation into a free-free substructure. Model uncertainty is parameterized in terms of test- or truth-analysis correlation metrics that are dictated by the Air Force. A statistical model is presented for these correlation metrics such that an analyst can specify a covariance matrix for uncertainty in model correlation at the fixed substructure level, and then propagate it into correlation uncertainty at the free substructure level. Development of the forward propagation approach then allows propagation of correlation uncertainty in the reverse direction from the free substructure into the fixed interface based Craig-Bampton representation. The proposed methods are applied to a typical spacecraft representation.     

Crystal plasticity-based constitutive modelling of irradiated bcc structures

A constitutive crystal plasticity model is proposed and developed for the inelastic deformation of irradiated bcc ferritic/martensitic steels. Defects found in these irradiated materials are used as substructure variables in the model. Insights from lower length- and time-scale simulations are used to frame the kinematic and substructure evolution relations of the governing deformation mechanisms. Models for evolution of mobile and immobile dislocations, as well as interstitial loops (formed due to irradiation), are developed. A rate theory-based approach is used to model the evolution of point defects generated during irradiation. The model is used to simulate the quasi-static tensile and creep response of a martensitic steel over a range of loading histories.     

Acoustic attenuation analysis of network systems by using impedance matrix method

An impedance matrix method is proposed to predict the acoustic attenuation characteristics of network systems. The system may contain several silencer modules and each module may be composed of complex components such as multiply connected tubes, portions with any-shaped cross-section and dissipative parts. The technique of substructuring is adopted and the system is divided into several substructure modules. Three strategies: boundary element method (BEM), numerical point collocation approach and numerical mode matching technique are introduced and the impedance matrix of each module may be computed by a certain appropriate methodology according to the dimensions and geometry of the substructure. Impedance matrix synthesis is employed to obtain the resultant impedance matrix and then transmission loss may be calculated. All the calculation results are verified by experimental measurements and 3-D BEM predictions.     

Substructure methods for structural condition assessment

It is commonly known that an accurate analysis of a large structure requires an accurate analytical model. This is also true for the inverse analysis of a structural system where measured structural responses are used as input to assess the structural conditions. However, an accurate model of the structure is always not available in practice. Two substructural identification methods are presented in this paper with the structure divided into substructures and with one substructure assessed at one time. In the first method, an accurate finite element model of the whole structure is assumed known. A state space method is applied to identify the external forces acting on the structure, and a damage identification method is then applied to identify the local damages using time domain information. Iterative model updating method based on the measured acceleration in the selected substructure is employed for the assessment. The second identification method requires only the finite element model of the substructure. The interface forces and the external forces acting on the target substructure are all taken as excitations and they are identified in state space. The substructure is then assessed similar to the first method. Since the target substructure for updating consists of a much reduced number of components and the identification problem is more efficient. The validation of the proposed methods is demonstrated by a truss structure with polluted measured accelerations with promising results.