共查询到20条相似文献,搜索用时 15 毫秒
1.
2.
本文基于积分方程法研究并建立了一种模拟横向同性介质中任意各向异性异常 体三维电磁响应的高阶广义扩展Born近似(Ho-GEBA)算法. 首先利用逐次迭代技术给出积分方程的广义级数展开解, 为保证其收敛性, 引入一种各向异性条件下满足压缩映射的迭代算子. 然后利用异常体区域分解技术, 并结合扩展Born近似原理, 得到各向异性介质三维电磁响应的Ho-GEBA解. 为提高效率, 计算过程中采用并矢Green函数的解析表达式. 最后通过数值计算实例对比验证了本文算法的有效性.
关键词:
高阶广义扩展Born近似
积分方程
电磁模拟
解析Green函数 相似文献
3.
提出求解三维静电场的三角形线性插值边界元解析积分方法.针对含1/R和1/R2的积分项,将单元形状函数分解为常数项、含x的线性项和含y的线性项,从而将边界单元积分简化为6个基本积分组合,并导出其解析计算公式,避免了因形状函数改变而导致的重复计算.该方法不仅可以准确计算远离奇异情况下的边界元积分,而且可以准确计算一阶和二阶接近奇异积分以及一阶奇异积分.计算结果表明,在接近奇异积分和奇异积分比较突出的问题中,当数值积分方法不能给出正确结果时,用同样的边界元网格,解析积分方法可以给出正确的结果,提高了三维静电场线性插值边界元法的计算精度. 相似文献
4.
朱寅弢吴海军孙瑞华蒋伟康 《声学学报》2023,(6):1218-1226
核函数中保留Dirac函数的原型,形成关于时间的卷积积分,是声学时域边界元法中一种稳定、有效的时间数值积分计算方法 (CQ-BEM)。然而,传统CQ-BEM中卷积积分系数的获取有计算量大、耗时长,且对不同单元需要重新计算的问题,极大地降低了CQ-BEM法计算时域声场的效率。针对传统CQ-BEM积分系数计算效率低的问题,本文利用多项式展开定理给出了待求函数泰勒系数的解析表达与数值计算方法,建立了不同单元间待求系数的转换理论,可以在一次循环迭代内完成不同单元的积分系数的计算,大幅降低了计算量,提高了CQ-BEM方法的声场计算效率。脉动球源数值算例结果表明,在相同要求下,本文方法计算时间较传统方法减少50%以上,相对误差小5个数量级以上,且计算时间随单元数的增长率仅为传统方法的2.34%。因此,本文提出的系数计算方法能够有效提高CQ-BEM方法的时域声场计算效率,拓展了CQ-BEM在大型机电设备时域声场模拟的计算规模。 相似文献
5.
为了充分考虑海底地形随三维空间变化的海洋环境中水平方位角耦合效应对声传播的影响,建立了一种三维柱坐标系下流体高阶抛物方程算法。该算法采用泰勒近似将二维方根算子分裂成一维方根算子,并采用分裂步进的高阶Pade近似将一维方根算子写成微分算子有理分式连乘的形式,进而应用Galerkin离散化方法来处理微分算子,最终将微分方程写成矩阵方程的形式;采用能量守恒近似来处理海底边界,以考虑复杂海底对于声传播的影响;采用交替方向隐式格式,实现了三维声场的步进计算。楔形和海底山等典型海域声场仿真计算表明,相比于已有的声场计算模型,三维柱坐标系下高阶抛物方程模型可以更加精确地计算楔形海域和海底山区域的三维声场,实现水平方位全空间声场计算。 相似文献
6.
本文针对边界元法在计算薄型结构力学、裂纹扩展等物理问题时存在的积分难题,提出一种基于自适应单元细分法的高效高精度近奇异域积分计算方法,该方法基于二叉树数据结构的单元细分技术对体单元进行自适应细分,消除单元几何形状所引起的近奇异性,能直接用于计算连续核函数的近奇异域积分。针对间断核函数的近奇异域积分,在细分单元的基础上采用腔面重建算法和投影算法,重新构建源点附近的积分子单元。数值算例表明:本方法可采用较少的积分点得到准确结果,是处理近奇异域积分的一种有效方法。 相似文献
7.
8.
提出三维导体目标与导体粗糙面复合散射的解析-数值混合迭代算法,推导出三维目标与粗糙面的耦合积分方程,以及粗糙面散射的Kirchhoff近似(KA)计算式.粗糙面的KA解析计算大大降低了粗糙面求解的复杂度,与目标矩量法的混合迭代保证了计算结果的精度,使得三维体-面目标复合散射计算变得可行.由于体-面两者的高阶耦合作用明显减小,保证了该混合迭代算法的收敛性.与镜像Green函数方法的比较表明该混合算法的有效性,并讨论了粗糙面长度选择对计算结果的影响.结合Monte-Carlo方法,数值分析了理想导体Gauss
关键词:
复合散射
Kirchhoff近似
共轭梯度法
互耦迭代 相似文献
9.
提出三维导体目标与导体粗糙面复合散射的解析-数值混合迭代算法,推导出三维目标与粗糙面的耦合积分方程,以及粗糙面散射的Kirchhoff近似(KA)计算式.粗糙面的KA解析计算大大降低了粗糙面求解的复杂度,与目标矩量法的混合迭代保证了计算结果的精度,使得三维体-面目标复合散射计算变得可行.由于体-面两者的高阶耦合作用明显减小,保证了该混合迭代算法的收敛性.与镜像Green函数方法的比较表明该混合算法的有效性,并讨论了粗糙面长度选择对计算结果的影响.结合Monte-Carlo方法,数值分析了理想导体Gauss 相似文献
10.
11.
在波叠加法中,结构外部声场是在离散边界上对Green函数进行积分并叠加得到,但数值积分的计算效率较低。而等效源法虽然提高了计算效率,但其面源简化为点源的过程中存在较大的积分近似误差。针对上述两种方法的缺陷,构造了一种波函数以替代离散单元关于Green函数积分的声场。首先,利用球坐标系下Helmholtz方程的解,推导了替代矩形单元积分的一般形式波函数及效率更高的内推波函数。其次,当离散单元为正方形时,将其近似成圆形域,进一步简化了内推波函数的表达式。最后,将所构造的波函数应用于声场计算。数值结果表明,在计算单个矩形单元外部辐射声场时,构造的波函数不仅保证了计算精度,而且相比于直接积分大幅度提高了计算效率。其中,矩形域一般形式和内推形式的波函数计算效率是直接积分的5~6倍,圆形域内推波函数计算效率达到了直接积分的12~13倍。在简支板声源和立方箱体辐射声源数值算例中,圆形域内推波函数在整个计算频段的声场计算精度均高于等效源法。 相似文献
12.
13.
提出了综合处理Burton-Miller方法所导致的奇异积分与近奇异积分问题的数值求积方法,以此改进了基于常量元素的常规边界元和低频快速多极边界元方法。对于奇异积分问题,利用Hadamard有限积分方法进行解决;对于近奇异积分问题,则采用极坐标变换法和PART方法(Projection and Angular&;Radial Transformation)进行克服。与解析解和LMS Virtual.Lab商业软件的结果比较验证了方法的正确性,并对比分析了奇异积分与近奇异积分对计算精度的影响。采用低频快速多极子方法以加速常规边界元法的计算效率,计算分析了计算复杂度,并成功实现了34万自由度大规模问题的计算。结果表明,近奇异积分问题主要由超奇异核函数引起,对计算精度的影响不容忽略;快速多极边界元法的精度与常规边界元法一致,但计算复杂度要远低于后者。 相似文献
14.
15.
提出了计算任意表面形状刚性边界目标散射的基于边界积分方程的T矩阵方法(TMM-BIE).利用Helmholtz积分方程法(HIEM)计算目标表面声场,替代扩展边界法(EBCM)计算中对目标表面声场的近似处理,解决了扩展边界法不能计算任意形状目标的散射T矩阵问题.文中计算了刚性边界的球目标、有限长圆柱目标以及非对称的三维散射体-猫眼(cat's-eye)模型的散射指向性和T矩阵.通过与解析解和HIEM结果比较,证明该方法的有效性. 相似文献
16.
提出了一种求取轴对称结构任意边界条件下声辐射特性的边界元方法。采用Burton和Miller改进型公式将高阶奇异项转化为弱奇异项之和,保证声辐射参数的唯一性,且计算简单精确。将结构表面声压与振速按照旋转轴角度进行Fourier级数展开,利用级数的正交性建立各项待定系数的求解公式;然后转化格林函数的法向偏导为切向偏导,方便直接计算各项积分,并将面积分公式表示为沿结构边界的线积分和沿旋转角度的积分;进一步采用二次等参单元离散结构边界线,建立声压与振速的关系矩阵,从而确定结构声辐射参数。以脉动球源和横向振动球源为例计算,与解析解和传统边界元法结果作对比,说明该方法的有效精确性。 相似文献
17.
特殊函数的数值计算综述 总被引:1,自引:0,他引:1
本文首先介绍了在科技工作中常遇到的各类特殊函数,以计算Bessel函数和椭圆积分为例,讨论了在计算大多数特殊函数中所遇到的有关算法的一些共性问题。如采用小宗量级数展开式,大宗量渐近式,多项式近似,递推法,迭代法计算等。最后还讨论了计算某些特殊函数的零点及其它有关数值计算问题。 相似文献
18.
19.
20.
基于共形完全匹配层(CPML)吸收边界,构造一种共形完全匹配层矢量单元按层积分算法,将多层单元积分运算叠加到一层单元中进行,用多层单元剖分,以一层单元计算矩阵元素,即保留了多层单元的几何和材料信息,又减少了单元数量和计算量.数值算例表明,这种按层积分CPML吸收效果好,计算量小,可靠,高效. 相似文献