首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 125 毫秒
1.
本文提出了求解位势问题的虚边界元法,建立了位势问题的虚边界元的离散方程式,推导了离散化求系数的积分解析式。该方法与传统边界元法相比具有不存在奇异积分和边界附近精度较高等优点,可用来计算真空静电场、稳定温度场、流体绕流、介质中的渗流等各类位势问题。大量算例均获得了满意的结果。  相似文献   

2.
求解位势问题的虚边界元法   总被引:5,自引:0,他引:5  
本文提出了求解位势问题的虚边界元法,建立了位势问题的虚边界元的离散方程式,推导了离散化求系数的积分解析式。该方法与传统边界元法相比具有不存在奇异积分和边界附近精度较高等优点,可用来计算真空静电场,稳定温度场,流体绕流,介质中的渗流等各类位势问题,大量算例均获得了满意的结果。  相似文献   

3.
提出了边界元法(BEM)的一种新的实现方法——边界面法(BFM)。在传统的边界元法中,单元不仅用来进行边界积分和函数插值,而且用来近似几何体。当离散网格较稀疏时,会引起较大几何误差,因而影响计算精度。本文基于参数曲面,将几何实体的边界曲面离散为参数空间里的曲面单元,边界积分和场变量的插值都是在曲面参数空间里进行。积分点...  相似文献   

4.
本文提出求解任意形状的薄板弯曲问题的虚边界元-最小二乘法。本法首先利用薄板弯曲平衡方程的格林函数和离开实际边界上分布的未知的横向荷载和法向弯矩函数建立满足实际边界条件的积分方程;然后采用最小二乘法和沿虚边界分段离散化的待定的分布横向荷载和法向弯矩函数得到求上述积分方程离散化数值解的线性代数方程组。导出了一系列的数值积分的公式,并求解了许多例题,数值结果说明本法完全避免了奇异积分及其复杂的处理方法和耗时的运算,而且在边界及其附近区域解的精度比普通边界元(以后简称边界元)法大大地提高了。  相似文献   

5.
采用双重互易边界元法结合精细积分法求解二维含热源的瞬态热传导问题。针对边界积分方程中热源项和温度关于时间导数项引起的域积分,采用双重互易法处理,将域积分转换为边界积分。采用边界元法将边界积分方程离散后,得到关于时间的微分方程组,并利用精细积分法处理其中的指数型矩阵;对于微分方程组中由边界条件和热源项引起的非齐次项,采用解析的方法计算。为了比较精细积分-双重互易边界元法的计算效果,同时使用有限差分法计算温度对时间的导数项。通过数值算例验证了本文方法的有效性和精确性。计算结果表明:时间步长对于精细积分-双重互易边界元法的结果影响较小,而有限差分法对时间步长比较敏感且只在时间步长选取较小时有效;当选取较大时间步长时,精细积分-双重互易边界元法依然具有良好的计算精度。  相似文献   

6.
采用超参非连续元离散三维弹性力学问题边界积分方程,借助三角极坐标变换方法处理奇异积分。将超参非连续元用于多域边界元分析,解决了自由度约束问题。提出了二次缩聚的概念,提高了多域缩聚边界元法的求解效率。通过数值算例表明了本文方法的可行性和有效性。  相似文献   

7.
依据弹性力学虚边界元法的基本思想和电磁弹性固体的基本解,提出了电磁弹性固体三维问题的虚边界元-等额配点法.该方法继承传统边界元法优点的同时,有效地避免了传统边界元法的边界积分奇异性的问题.算例表明该方法有很高的精度,是求解电磁弹性固体三维问题的一个有效的数值方法.  相似文献   

8.
本文基于有限元法、边界元法和虚拟激励法,对随机激励下结构振动声辐射问题进行研究。提出了一种计算随机激励下结构振动声辐射问题的新方法,其中,有限元法用于计算结构谐振响应,边界元法用于计算结构振动声辐射,虚拟激励法结合有限元和边界元计算随机激励下结构振动声辐射问题。 数值算例表明,本文方法在计算精度上与传统方法等价,且更具高效性。  相似文献   

9.
利用非连续元离散边界积分方程,有效地解决了“角点效应”问题,对影响非连续元精度和分析效率的几个问题从数值计算的角度进行了讨论,将非连续边界元用于自适应边界元分析,给出了自适应边界元误差指示确定的一种方法,通过对具体实例分析表明了给所方法的可行性,。》  相似文献   

10.
复杂载荷三维裂纹分析双重边界元法   总被引:11,自引:1,他引:10  
陆山  黄其青 《力学学报》2002,34(5):715-725
提出可用于高温、高转速状态下的热动力机械三维含裂构件热弹性分析方法——双重边界元法.首先建立了考虑温度及离心载荷的双重边界积分方程组,并对边界积分方程组的选取及适用范围进行了讨论。然后提出角非快调元模型离散技术。接着提出超奇异积分方程分析去除奇异性方法及数值积分技术.数值算例表明计算结果与有关权函数解十分吻合,说明了用双重边界元法计算复杂载荷条件下三维应力强度因子的有效性.还讨论了有关热应力强度因子权函数解的适用范围.  相似文献   

11.
A boundary element method is proposed for studying periodic shallow water problems. The numerical model is based on the shallow water equation. The key feature of this method is that the boundary integral equations are derived using the weighted residual method and the fundamental solutions for shallow water wave problems are obtained by solving the simultaneous singular equations. The accuracy of this method is studied for the wave reflection problem in a rectangular tank. As a result of this test, it has been shown that the number of element divisions and the distribution of nodes are significant to the accuracy. For numerical examples of external problems, the wave diffraction problems due to single cylindrical, double cylindrical and plate obstructions are analysed and compared with the exact and other numerical solutions. Relatively accurate solutions are obtained.  相似文献   

12.
A formulation is developed to impose pressure-prescribed boundary conditions in the penalty finite element method. Some numerical experiments for the Poiseuille flow problem are performed to compare it with the conventional traction-prescribed boundary condition. Also the incorrectness of the traction-free outlet boundary condition for contained-flows is studied with explanatory numerical examples. Discussion is focused on the inlet and outlet boundary conditions to simulate fully developed flows. Finally, the three-dimensional flow in a bifurcated pipe is analysed with the proposed formulation.  相似文献   

13.
本文将D.Nardini和C.A.Brebbia所提出的一种动态边界积分方程新解法应用于动态断裂力学数值计算,对数值实现问题,尤其对数值稳定性及精度问题进行了详细研究,得到了能保证数值稳定性的数值解法,给出了动态断裂力学计算实例,同己有的数值结果比较,表明本文的计算是成功的。  相似文献   

14.
YuanYiwu(袁镒吾)(ReceivedOct.2,1994;CommunicatedbyChienWeizang)INTERPOLATIONPERTURBATIONMETHODFORSOLVINGTHEBOUNDARYLAYERTYPEPROB...  相似文献   

15.
The effect of small viscosity on nearly inviscid flows of an incompressible fluid through a given domain with permeable boundary is studied. The Vishik–Lyusternik method is applied to construct a boundary layer asymptotic at the outlet in the limit of vanishing viscosity. Mathematical problems with both consistent and inconsistent initial and boundary conditions at the outlet are considered. It is shown that in the former case, the viscosity leads to a boundary layer only at the outlet. In the latter case, in the leading term of the expansion there is a boundary layer at the outlet and there is no boundary layer at the inlet, but in higher order terms another boundary layer appears at the inlet. To verify the validity of the expansion, a number of simple examples are presented. The examples demonstrate that asymptotic solutions are in quite good agreement with exact or numerical solutions.  相似文献   

16.
平均源边界节点法ASBNM是一种最近提出的边界型无网格法。该方法仅使用边界节点不涉及任何单元和积分的概念,具有方法简单和程序设计容易等特点。但是,对于依赖于边界积分方程的边界型无网格法,关键问题是如何准确高效地估计影响矩阵的对角元。本文提出直接计算影响矩阵对角元的方法,是已有ASBNM法的改进,将对角元的计算转化为一个纯几何问题,因此适用于任何二维边值问题。数值算例证明了本文方法的有效性和准确性。  相似文献   

17.
A numerical method is developed to solve the coupled unsteady laminar momentum and thermal boundary layers over a circular cylinder impulsively started from rest. The present non-iterative finite difference method, which requires relatively fewer grid points in the reversed flow region than any other method, can easily handle the separating boundary layer flows. The results indicate that the present method has accuracy comparable with the earlier methods, while consuming computer time approximately one order of magnitude less. The present numerical method allowed investigation of the effect of buoyancy parameter on the starting boundary layer. The time-dependent behaviour of the boundary layer is studied in terms of the appearance of the singularity, the distribution of skin friction and wall heat flux, and the wall position of the inflection point of the velocity profile. The transient as well as buoyancy-dependent patterns of the streamlines and isotherms are also studied.  相似文献   

18.
Summary A variant of the boundary element method, called the boundary contour method, offers a further reduction in dimensionality. Consequently, boundary contour analysis of 2-D problems does not require any numerical integration at all. In a boundary contour analysis, boundary stresses can be accurately computed using the approach proposed in Ref. [1]. However, due to singularity, this approach can be used only to calculate boundary stresses at points that do not lie at an end of a boundary element. Herein, it is shown that a technique based on the displacement/velocity shape functions can overcome this drawback. Further, the approach is much simpler to apply, requires less computational effort, and provides competitive accuracy. Numerical solutions and convergence study for some well-known problems in linear elasticity and Stokes flow are presented to show the effectiveness of the proposed approach. This research was supported in part by the 2004 Ralph E. Powe Junior Faculty Enhancement Award from Oak Ridge Associated Universities and by the University of South Alabama Research Council.  相似文献   

19.
Using Green’s functions, the extended general displacement solutions of a three-dimensional crack problem in anisotropic electro-magneto-elastic (EME) bimaterials under extended loads are analyzed by the boundary element method. Then, the crack problem is reduced to solving a set of hypersingular integral equations (HIE) coupled with boundary integral equations. The singularity of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of HIE, and the exact analytical solutions of the extended singular stresses and extended stress intensity factors (SIFs) near the crack front in anisotropic EME bimaterials are given. Also, the numerical method of the HIE for a rectangular crack subjected to extended loads is put forward with the extended crack opening dislocation approximated by the product of basic density functions and polynomials. At last, numerical solutions of the extended SIFs of some examples are obtained.  相似文献   

20.
张赞  程玉民 《力学季刊》2007,28(2):333-339
无网格方法与有限元法或边界元法耦合是无网格方法处理边界条件的方法之一,在无网格方法中研究无网格方法与有限元法或边界元法耦合的研究显得非常重要.本文在无单元Galerkin法和边界元法的基础上,基于无单元Galerkin法子域和边界元法子域的界面上位移连续和面力平衡条件,提出了一种新的无单元Galerkin法和边界元法的直接耦合方法,对弹性力学问题详细推导了在整个求解域上的耦合公式.与以往的耦合法相比,这种方法简单直观,不需要增加新的耦合区域,也不需要建立新的逼近函数来保证界面位移的连续性.算例结果表明,该方法具有较好的计算精度.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号