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

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

用边界元法分析物体的应力

针对固体力学中的应力问题,应用边界元法建立基本方程,并求解此方程得出求解线性弹性物体三维应力的式子.以便工程中使用计算机时的结构应力进行数值计算.     

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

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

结构自振特性的边界单元法分析

本文阐述了用边界单元法计算结构自振频率与振型的原理与具体实施方法,给出了计算支配方程系数矩阵元素用的特殊函数的逼近公式,用求复行列式模的局部极小值方法计算各阶自振频率,用最小二乘法计算相应的各阶振型.计算实例表明,用边界单元法计算结构自振特性,前处理简单,精度较高,不失为一种有效的方法.     

弹塑性结构安定性分析数值方法及应用

针对极限与安定理论上限分析中存在的问题,建立了复杂变化载荷作用下弹塑性结构安定上限分析的有限元数学规划格式. 利用研究结构在基准载荷域各个角点处安定的办法,克服了机动定理中对时间积分的困难,提出了一种直接迭代算法求解,以克服目标函数非线性非光滑所导致的困难. 该格式同时考虑了温度对材料屈服极限的影响. 算例表明本文作者所提出的安定上限分析算法具有计算效率高、收敛性好和数值精度高等优点. 由于采用了位移模式有限元,因而具有较广的适用范围,可用于复杂承压结构的工程分析.     

结构的三维非线性界面元分析

评价了用于含软弱夹层岩体弹塑性分析时常用的Ｇｏｏｄｍａｎ单元及其变形形式的单元。结合界面元法的特点，在单元交界面上引入了能模拟材料弹性及塑性变形的弹簧及滑块元件，使界面元法很好地解决了材料的弹塑性分析而无需另设置任何形式的夹层单元。文中导出了反映层状岩体横观各向同性材料的界面应力公式和非线性界面元的计算公式，并制了相应的计算程序。数值考题表明，将界面元法用于材料的非线性分析，尤其对含有各种节理面的     

用边界单元法求解弹性结构的动力响应

用与时间无关的Ｋｅｌｖｉｎ问题的基本解,作为加权函数的边界单元法求解弹性结构的动力响应．选用一组线性无关的坐标函数来近似域内点的位移,使惯性项的域积分转化为边界积分,把复杂的结构动力响应问题转化为边界上求解二阶线性常微分方程组的问题．利用Ｈｏｕｂｏｌｔ直接积分方法对时域进行离散,由初始条件逐步求出一系列离散时刻弹性结构的动力响应．文中的算例证实了该方法的可行性与精确度     

Analysis of 3-D Frictional Contact Mechanics Problems by a Boundary Element Method ＊

The development of two boundary element algorithms for solving 3-D, frictional, and linear elastostatic contact problems is reported in this paper. The algorithms employ nonconforming discretizations for solving 3-D boundary element models, which provide much needed flexibility in the boundary element modeling for 3-D contact problems. These algorithms are implemented in a new 3-D boundary element code and verified using several examples. For the numerical examples studied, the results using the new boundary element algorithms match very well with the results using a commercial finite element code, and clearly demonstrate the feasibility of the new boundary element approach for 3-D contact analysis     

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

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

边界元中的角点处理方法

本文对边界元中的角点效应(corner effect)问题和处理方法进行了分析,提出了降阶多重节点法。理论分析与实例计算表明,这是一种合理的角点处理方法,它可提高计算精度、降低计算工作量。     

新型曲面四边形边界元精细后处理方法研究

为了精确计算三维静电场的电场强度和电位分布,提出了新型曲面四边形边界元方法．在该方法中,对模型边界面进行二阶四边形单元剖分,对二阶单元顶点上的节点号重新编号,以单元的顶点为求解点,根据二阶四边形曲面参数方程,结合面积比值法定义的曲面单元顶点的形状函数,计算曲面单元顶点的函数值．与一阶平面四边形边界元相比,新型曲面边界元法在没有增加计算节点的情况下,由于采用更接近实际边界的曲面积分,计算精度将明显提高．但由于边界面采用二阶单元粗略剖分,单元数量相对较少,剖分后的模型较粗糙．虽然顶点节点上的函数值比较精确,但只能以平面线性单元的形式显示,离实际模型边界差别较大．本文就此提出边界元精细后处理方法．在该方法中,对曲面单元两边按一定步长等分,再根据曲面的参数方程把曲面单元精细显示出来．单元上新建节点的函数值可由曲面单元顶点上的函数值和面积比值法定义的形状函数插值得到．最后形成经精细显示后的新型曲面边界元方法．算例表明,经精细显示后边界面比未处理前更接近实际边界．     

带横孔圆轴三维应力分析的边界元法

本文用弹性力学边界积分方程边界元法解带根孔圆轴的三维应力集中问题。简述了方程的离散化及有关的数值技术。以带径向圆孔的圆轴为例，采用圆柱面、双线性插值边界元计算应力集中系数。初步计算结果与现有某些设计资料中的实验曲线相比较是接近的，这表明边界元法为改进及扩充工程实用应力集中数据开辟了一个途径。     

概率边界元法结构分析

本文讨论了在常规边界元法中引入概率分析的方法;给出了能够反应边界条件随机变化的边界元基本公式以及对由此而导出的结果进行概率分析的后处理方法;本方法在数值计算方面只需对原有常规边界元程序稍加改动便可实现。     

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

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

复合材料层合板非线性稳定性分析的边界元法

应用边界元法对多向铺设的复合材料层合板结构进行非线性稳定性分析。在分析中计及了几何非线性和物理非线性的影响,导出了层合板在纵向轴压力作用下屈曲的控制方程,采用双重傅立叶级数建立了边界积分方程。通过对几种玻璃／环氧复合材料层合板铺层情况的计算得到临界力,并获得一些有意义的结论。     

大型结构静力分析的界面关连矩阵法

本文提出了大型结构静力分析的一种新方法－界面关连矩阵法，此方法的优点是不仅可以求解界面自由度相等的问题，而且也可以求解界面自由度不等的问题，适用范围大于有限元－传递矩阵法，且有一定的实际应用价值     

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

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

采用有限元-边界元耦合方法计算弹塑性应力

采用有限元－边界元耦合方法对高压三通进行弹塑性分析．在处于弹塑性状态的连 接处附近区域使用边界元方法，其它部分采用有限元方法．给出了一种界面处理方法。计 算实例表明此方法有较好的计算精度和效率．     

关于边界元法及其解的精度问题探讨

本文阐述了边界元直接法与间接法的基本原理及步骤,提出了改进边界元解在边界附近精度的方法。     

Mixed Finite Element Method and Higher-Order Local Artificial Boundary Conditions for Exterior 3-D Poisson Equation

IntroductionAnalysis of partial differential equations on anunbounded domain often requires artificialboundaries to limit the problem to a boundedcomputational domain.Such situations can arise inmany applications such as geophysical calculationsinvolving …