二维位势边界元法高阶单元几乎奇异积分半解析算法

准确计算几乎奇异积分是边界元法难题之一。目前,对于一般的高阶单元的几乎奇异积分尚缺乏通用高效的计算方法。本文在单元局部坐标系中表征了二维高阶单元的几何特征,提出了源点相对高阶单元的接近度概念。针对二维位势边界元法的3节点二次等参单元,构造出与单元积分核具有相同几乎奇异性的近似奇异核函数。从二维位势几乎奇异积分单元积分核中扣除近似奇异核函数,把几乎奇异积分项转换为规则积分和奇异积分两部分之和,规则积分部分用常规Gauss数值积分计算,奇异积分部分由导出的解析公式计算,从而建立了二维位势问题高阶单元几乎强奇异和超奇异积分的半解析算法。算例结果表明了本文半解析算法的有效性和计算精度。     

A NEW SEMI-ANALYTIC ALGORITHM OF NEARLY SINGULAR INTEGRALS IN HIGH ORDER BOUNDARY ELEMENT ANALYSIS OF 2D ELASTICITY

针对边界元法中高阶单元中几乎奇异积分计算难题,解剖了二维边界元法高阶单元的几何特征,定义源点相对高阶单元的接近度。将高阶单元上奇异积分核函数用近似奇异函数逼近,从而分离出积分核中主导的奇异函数部分,其奇异积分核分解为规则核函 数和奇异核函数两项积分之和。规则核函数用常规高斯数值积分,再对奇异核函数积分导出解析公式,从而建立了一种新的半解析法,用于高阶边界单元上几乎强奇异和超奇异积分计算。给出3个算例,采用边界元法高阶单元的半解析法计算了弹性力学薄体结构和近边界点位移/应力,并与线性边界元正则化算法结果作了比较,结果表明提出的二次元的半解析算法更加有效。特别是分析薄体结构,采用正则化算法的线性边界元分析比有限元有显著优势,而用提出的二次边界元半解析算法分析比其线性元的有效接近度又减小了4个量级。     

三维边界元法高阶元几乎奇异积分半解析法

分析了三维边界元法高阶曲面单元几何特征,定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系,构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数,分离出积分核中主导的奇异函数部分,将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算,奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ,对ρ积分建立解析计算列式,对θ积分应用常规Gauss数值积分计算,从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例,采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度,并与线性单元正则化算法计算结果对比,结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效.      

Helmholtz边界积分方程中奇异积分间接求解方法

提出了间接求解传统Helmholtz边界积分方程CBIE的强奇异积分和自由项系数,以及Burton-Miller边界积分方程BMBIE中的超强奇异积分的特解法。对于声场的内域问题,给出了满足Helmholtz控制方程的特解,间接求出了CBIE中的强奇异积分和自由项系数。对于声场外域对应的BMBIE中的超强奇异积分,按Guiggiani方法计算其柯西主值积分需要进行泰勒级数展开的高阶近似,公式繁复,实施困难。本文给出了满足Helmholtz控制方程和Sommerfeld散射条件的特解,提出了间接求出超强奇异积分的方法。推导了轴对称结构外场问题的强奇异积分中的柯西主值积分表达式,并通过轴对称问题算例证明了本文方法的高效性。数值结果表明,对于内域问题,采用本文特解法的计算结果优于直接求解强奇异积分和自由项系数的结果,且本文的特解法可避免针对具体几何信息计算自由项系数,因而具有更好的适用性。对于外域问题,两者精度相当,但本文的特解法可避免对核函数进行高阶泰勒级数展开,更易于数值实施。     

A NEW SEMI-ANALYTIC ALGORITHM OF NEARLY SINGULAR INTEGRALS IN HIGH ORDER BOUNDARY ELEMENT ANALYSIS OF 3D POTENTIAL

分析了三维边界元法高阶曲面单元几何特征,定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系,构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数,分离出积分核中主导的奇异函数部分,将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算,奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ,对ρ积分建立解析计算列式,对θ积分应用常规Gauss数值积分计算,从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例,采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度,并与线性单元正则化算法计算结果对比,结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效.     

弹性板边界元中计算边界高阶导数的渐近收敛积分方程

本文针对板弯曲边界元方法中计算边界曲率等高阶导数项时边界积分方程中出现的高阶奇异积分项,通过对未知挠曲函数作渐近展开并加以适当摄动,获得了渐近收敛的边界积分方程。采用这一方法计算板边界上的曲率分布,获得了满意的数值结果。     

三维非规则非均匀边界元网格的简便的高精度算法

对三维直接边界元中一阶奇异积分、一阶近奇异积分以及非奇异积分进行统一处理，给出了一种提高积分计算精度的简便有效的方法，对非规则非均匀边界元网格可获得比一般方法高得多的计算精度，非常适合边界形状比较复杂的三维实际问题的边界元分析．     

二维正交异性位势问题高阶边界元几乎奇异积分半解析算法

边界元法中高阶单元上的几乎奇异积分一直难以计算。针对正交各向异性位势问题,提出一个半解析算法准确计算了其高阶单元上的几乎奇异积分。首先将正交各向异性材料中源点到单元的距离函数在局部坐标系下渐近展开,采用级数展开式构造出与奇异积分核函数具有相同奇异性的可积近似核函数;然后利用扣除法的思想,原奇异积分核减去近似积分核后再加回,几乎奇异积分便转换为规则部分和奇异部分之和,规则积分采用Gauss数值积分计算,奇异积分由文中推导出解析公式计算。通过两个正交各向异性的热传导算例表明,本文建立的高阶单元半解析算法能准确高效地计算近边界内点位势和位势梯度。     

尖端具有线性粘着力作用的加层空间的界面裂纹对稳态弹性波的散射

本文研究一类粘着型界面裂纹的弹性波散射问题．文中利用积分变换和积分方程方法推导了确定这类问题的奇异积分方程组．采用围道积分技术和切比雪夫多项式展开技术，得到了待定系数的非线性代数方程组．最后本文给出裂纹尖端粘着区的大小和界面应力的数值结果．     

含非完整界面的功能梯度压电材料Ⅲ型裂纹问题

本文研究了含非完整界面的功能梯度压电复合材料的Ⅲ型裂纹问题．此裂纹垂直于非完整界面,采用弹簧型力电耦合界面模型模拟非完整界面．界面两侧材料的性质,如弹性模量、压电常数和介电常数均假定呈指数函数形式且沿着裂纹方向变化．运用积分变换法将裂纹面条件转换为奇异积分方程,并使用Gauss-Chebyshev方法对其进行数值求解．根据算例结果讨论了一些退化问题并分析了裂纹尖端强度因子与材料的非均匀系数和非完整界面参数的关系．     

Preface

This special issue of PARTICUOLOGY is devoted to the first UK-China Particle Technology Forum taking place in Leeds, UK, on 1-3 April 2007. The forum was initiated by a number of UK and Chinese leading academics and organised by the University of Leeds in collaboration with Chinese Society of Particuology, Particle Technology Subject Group （PTSG） of the Institution of Chemical Engineers （IChemE）, Particle Characterisation Interest Group （PCIG） of the Royal Society of Chemistry （RSC） and International Fine Particle Research Institute （IFPRI）. The forum was supported financially by the Engineering and Physics Sciences Research Council （EPSRC） of United Kingdom,     

The 7th International Conference on Fluid Mechanics May 24-27,2015,Qingdao,China

正http://www.icfm7.org First Announcement and Call for PapersThe objective of International Conference on Fluid Mechanics(ICFM)is to provide a forum for researchers to exchange new ideas and recent advances in the fields of theoretical,experimental,computational Fluid Mechanics as well as interdisciplinary subjects.It was successfully convened by the Chinese Society of Theoretical and Applied Mechanics(CSTAM)in Beijing(1987,     

INFORMATION FOR CONTRIBUTORS

Contributions： The Journal, Acta Mechanica Solida Sinica, is pleased to receive papers from engineers and scientists working in various aspects of solid mechanics. All contributions are subject to critical review prior to acceptance and publication.     

Guide for authors

正Each of the sections below provides essential information for authors.We recommend that you take the time to read them before submitting a contribution to Acta Mechanica Sinica.We hope our guide to authors may help you navigate to the appropriate section.How to prepare a submission This document provides an outline of the editorial process involved in publishing a scientific paper in Acta Mechanica Sinica.     

基于鲁棒滤波的捷联导引头视线角速度估计方法

针对捷联导引头无法直接获取视线角速度等信息的问题,研究了鲁棒滤波在大气层外飞行器捷联导引头视线角速度估计中的应用。为了建立非线性滤波估计模型,考虑目标视线角速度的慢变特性,采用一阶马尔科夫模型建立了状态方程;推导了视线角速度的解耦模型,并建立了量测方程;考虑到实际应用中存在系统噪声统计特性失准的问题,基于Huber-Based鲁棒滤波方法,设计了视线角速度滤波器,并完成了基于Huber-Based滤波方法和扩展卡尔曼滤波方法的数学仿真。仿真结果表明Huber-Based滤波方法的视线角、视线角速度及视线角加速度估计精度分别达到0.1140'、0.1423'/s、0.0203'/s2,而扩展卡尔曼滤波方法的视线角、视线角速度及视线角加速度估计精度仅分别为0.6577'、0.6415'/s、0.0979'/s~2。仿真结果证明了该方法可以有效地估计出相对视线角速度等信息,并且在非高斯噪声的条件下,依然可获得较高的估计精度,具有一定的鲁棒性。     

Guide for authors

正Each of the sections below provides essential information for authors.We recommend that you take the time to read them before submitting a contribution to Acta Mechanica Sinica.We hope our guide to authors may help you navigate to the appropriate section.How to prepare a submission This document provides an outline of the editorial process involved in publishing a scientific paper in Acta Mechanica