用于各向异性压电材料的边界元法

给出一种用于压电材料静态问题的简单积分方程式.将压电体的控制方程按分离出各向同性及非耦合加权的形式进行重写,使各向同性弹性和势场问题的基本解可以通过积分方程列式表达,而另外两个补充方程通过压电材料的本构方程加以表达.最终,将耦合的各向异性问题转化为求解包含4个未知矢量的矩阵方程组.通过一个压电梁的计算实例说明了此方法的适用性.     

三维接触边界元法的一种误差直接估计

将作者所在研究组提出的二维弹性力学问题边界元解误差的直接估计推广到三维问题,给出了确定与域内解连续的边界位移的一种精确有效的方法。在此基础上提出将接触体接触单元间与域内解连续的边界位移之差的某种度量作为三维弹性接触问题边界元法的一种误差直接估计,并且提出了三维弹性接触问题的一种自适应边界元法计算方案。这种方案为确定没有解析解可作比较的复杂接触问题的边界元解精度提供了可能。文中对于三维弹性接触问题,给出了一个计算误差直接估计及自适应边界元法的算例。     

Kirchhoff圆板的解析积分边界元法

利用一般弯曲薄板边界为规则曲线的特点,对工程常用的圆形弯曲薄板,采用线性单元,导出Kirchhoff圆板各辅助态的边界积分解析表达式。建立问题的边界元法系统方程。从而使薄板的边界元分析完全避免通常使用的高斯积分。明显提高计算精度．给出4个不同荷载及边界条件情况的圆板的算例,计算结果表明。对于具有规则曲线边界的问题,采用解析积分的边界元法是十分有效的。     

Fast Multipole BEM for Simulation of 2-D Solids Containing Large Numbers of Cracks

The fast multipole method was used to solve the traction boundary integral equation for 2-D crack analysis, The use of both multipole and local expansions reduces both the computational complexity and the memory requirement to O(N). The multipole expansion uses a complex Taylor series expansion to reduce the number of multipole moments, The generalized minimum residual method solver (GMRES) was selected as the iterative solver, An improved preconditioner for GMRES was developed which uses less CPU time and less memory. A new initial candidate vector for the iterative solver was developed to further improve the efficiency, The numerical examples apply the method to the analysis of cracks in infinite 2-D space with the largest model having 900 000 degrees of freedom.     

近似边界元方法及其收敛性分析

以Ｌａｐｌａｃｅ方程Ｄｉｒｉｃｈｌｅｔ问题为例,为椭圆边值问题近似边界元法的建立及其收敛性分析提供了一种框架性的工作。文中给出了近似基本解的确定方法,具近似基本解的离散边界变分方程解的存在惟一性定理以及近似解的误差估计,特别给出了近似基本解中截断数和离散网格宽度应保持的匹配关系,文末给出了数值算例。     

多域直接边界元方法及其在模拟亥姆霍兹共振器方面的应用

多域直接边界元方法被用来求解结构表面带开口时的结构声学问题．详细地介绍了该方法的理论与求解过程．在开口处可任意选取一层流体介质与结构表面一起构成多域问题的边界．针对各域分别运用直接边界元方法,各域的边界元方程通过开口边界耦合起来．开口处多域边界选取的任意性被用来验证尖棱角处法向量处理方法的可行性．模拟了两种球形亥姆霍兹共振器．计算所得的共振频率与理论近似解非常吻合．     

压电材料平面问题应力强度因子的边界元计算

从压电介质的基本方程出发,利用功的互等原理,推导了边界积分方程．利用边界元分区算法研究了压电材料无限大平面含中心裂纹的断裂计算,在裂纹尖端采用１／４面力奇异单元,分别计算了应力强度因子ＫⅠ、ＫⅡ和电强度因子Ｋｅ,并讨论了电强度因子Ｋｅ随电载荷的变化规律．计算结果和理论解相吻合．     

Shakedown Analysis of 3-D Structures Using the Boundary Element Method Based on the Static Theorem

The static shakedown theorem was reformulated for the boundary element method (BEM) rather than the finite element method with Melan‘s theorem, then used to develop a numerical solution procedure for shakedown analysis. The self-equilibrium stress field was constructed by a linear combination of several basis self-equilibrium stress fields with undetermined parameters. These basis self-equilibrium stress fields were expressed as elastic responses of the body to imposed permanent strains obtained using a 3-D BEM elastic-plastic incremental analysis. The lower bound for the shakedown load was obtained from a series of nonlinear mathematical programming problems solved using the Complex method. Numerical examples verified the precision of the present method.     

热弹性力学边界元中二次元的几乎奇异积分计算

采用正则化积分算法,计算了二维热弹性力学边界元法中近边界点的几乎奇异积分。算法采用二次元划分边界,但对与内点邻近的二次单元,几何量采用线性插值,位移、面力等物理量仍采用二次插值。对此二次非等参单元上的积分采用正则化积分公式。算例证明了该文算法的有效性和精确性。     

屏蔽带状线特性阻抗的边界元分析程序

一、引言带状传输线的理论和应用已有三十余年的历史。早期关于非屏蔽(接地板无限宽)带状线和薄导带(厚度趋于零)屏蔽带状线的分析,多采用复位函数的保角变换技术,导出其特性阻抗的解析表达式。对于厚导带的屏蔽带状线,当导带置于正常的中心位置时又称矩形同轴线,也已经借助施瓦兹变换导出了特性阻抗的严格解;但当导带偏离中心位置、或不平行于外导体边界时,只能采用各种近似方法求解。Brackelmann曾用分区本征函数展开法分析导带二维偏位的屏蔽带状线,给出有关的阻抗曲线,但读数标度过于简略。有限差分和有限元法都可用于分析这些传输线,然而近年来迅速发展的边界元法具有计算量     

基于设计灵敏度分析的圆柱壳声辐射优化

采用有限元法计算结构振动,边界元法求解结构振动引起的辐射声场问题.在此基础上提出了圆柱壳结构声辐射优化模型,重点推导了结构振动声辐射灵敏度计算公式.优化模型中,声辐射功率作为目标函数,约束函数中不仅有内部场点平均声压作为声场设计指标,而且包含了质量和自振频率的结构设计指标,圆柱壳的厚度作为设计变量.通过与用差分法计算的...     

基于边界元法的高炉炉底炉缸侵蚀模型

通过边界元法建立了高炉炉底炉缸传热数学模型．采用基尔霍夫变换把非线性问题转化为线性问题，解决了利用边界元法建立高炉炉底炉缸侵蚀模型把导热系数看成常数而造成计算精度下降的问题．求解控制高炉炉底炉缸传热过程的热传导方程，再通过正交试验的方法确定满足实测边界温度分布的侵蚀边界．该模型可在线预测高炉炉底炉缸1150℃等温线的位置和形状，以了解和分析炉底炉缸的破损情况．结果表明，监测点热电偶温度值和模型计算值吻合较好．     

二维多重高和方程的BEM（法）及其应用

对二维多重非齐次调和方程应用BEM求解时将会出现域内积分项，在假定非齐次项为m次调和的情况下，将域内积分转化为边界积分，从而给出了相应的不含域内积分项的边界和积分方程，从而可化为二次齐次问题讨论。     

二维弹性力学边界条件反识别TSVD正则化法

针对二维各向同性弹性力学Cauchy问题,文章采用线性单元对边界积分方程进行离散,再引入已知的边界条件,得到包含所有待求边界条件信息的线性病态方程组。采用截断奇异值分解正则化技术求解该病态方程组,并使用L曲线法选择最优正则化参数,即奇异值截断位置,从而得到方程组的解。通过数值算例对求得的边界条件数值解与解析解进行比较,并进行误差分析,以表明截断奇异值分解算法的有效性和稳定性。通过减少已知数据中的随机偏差和增加边界单元密度可提高求解的精确度。     

基于本征正交分解加速的二维位势问题边界元法分析

基于本征正交分解(POD)的思想,采用奇异值分解法(SVD),将物理问题的求解空间分解为几何子空间和设计子空间,通过线性组合几何与设计子空间获得随机变量响应结果.与传统加速算法不同,采用径向基函数(RBF)近似设计子空间响应,实现了系统信息的压缩表达,有效降低了计算成本.采用边界元法(BEM)求解二维位势问题,并结合S...     

静,动态基本解及奇异积分的探讨

本文首先讨论了弹性力学中动力基本解和静力基本解的关系,得到了对于空间问题令频率ω→0可以从动力基本解退化到静力基本解,对于平面问题则不能从动力基本解退化为静力基本解的结论;其次求得了边界单元法中的奇异积分的解析表达式,从而证明了不管奇点是否位于角点、奇异积分在Cauchy主值意义下总是存在的问题.文中还给出了数值验证结果.     

PZT／环氧树脂0—3型压电复合材料性能的研究

制备了以环氧树脂为基体，不同ＰＺＴ材料为压电相的两种０－３型压电复合材料。研究了极化工艺参数和陶瓷粉末对压电系数ｄ－３３的影响，结果表明，复合材料的压电系数ｄ－３３先随极化合场Ｅ和极化时间ｔ的增加而增加，之后趋于恒定。ＰＺＴ体积含量增加，介电常数ε２减小，则复合材料压电系数ｄ－３３增加。     

二维Laplace方程边界元方法的误差估计

本文利用边界元方法来解决R~2中的Laplace问题,先给出该问题相应积分方程的误差估计,然后利用此及其近似解的构造,导出解及其导数的渐近误差估计.     

含界面中心裂纹的压电材料反平面问题

研究了含界面中心裂纹的不同压电材料在反平面剪切栽荷和平面内电场作用下的反平面问题．得到了用级数表示的满足控制拉普拉斯方程和可导通边界条件的基本解及应力强度因子．最后用边界配置法求解了应力强度因子与截面几何尺寸之间的关系．结果表明，边界配置法计算简便，具有广泛的应用性．     

汽车空调出风管道气动噪声分析与控制

通过耦合CFD(Computational Fluid Dynamics)与专业声学代码SYSNOISE求解汽车空调管道气动噪声,即利用LES(Large Eddy Simulation)湍流模型对空调管道的瞬态流场进行求解获得噪声源项,然后将噪声源项作为边界条件导入SYSNOISE来计算噪声的传播.根据流场分析与声场分析结果对空调管道的结构提出了两种改型方案,并对改型前后的空调系统噪声进行了测试.测试结果表明相比原始空调系统,两种方案都能有效降低噪声且方案二效果更好,驾驶员附近的噪声最大降幅达4.5 dB.