首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到17条相似文献,搜索用时 171 毫秒
1.
弹性力学问题的局部边界积分方程方法   总被引:21,自引:0,他引:21  
龙述尧  许敬晓 《力学学报》2000,32(5):566-578
提出了弹性力学平面问题的局部边界积分方程方法。这种方法是一种无网格方法,它采用移动最小二乘近似试函数,且只包含中心在所考虑节点的局部边界上的边界积分。它易于施加本质边界条件。所得系统矩阵是一个带状稀疏矩阵。它组合了伽辽金有限元法、整体边界元法和无单元伽辽金法的优点。该方法可以容易推广到求解非线性问题以及非均匀介质的力学问题。计算了两个弹性力学平面问题的例子,给出了位移和能量的索波列夫模,所得计算结果证明:该方法是一种具有收敛快、精度高、简便有效的通用方法。  相似文献   

2.
加权最小二乘无网格法   总被引:29,自引:0,他引:29  
张雄  胡炜  潘小飞  陆明万 《力学学报》2003,35(4):425-431
在最小二乘法和移动最小二乘近似的基础上提出了加权最小二乘无网格法.该方法除节点外又引入了一些辅助点,控制方程在所有节点和辅助点处的残差用最小二乘法予以消除,边界条件用罚函数法引入.另外对移动最小二乘近似进行了改进,并给出了最小二乘法中泛函的简化格式,因而提高了计算效率.与配点法相比,新方法精度高,稳定性好,并且系数矩阵是对称正定矩阵.与Galerkin法相比,该方法不需要进行高斯积分,因而计算量小.算例表明该方法具有效率高、精度高和稳定性好等优点,并且易于实现.  相似文献   

3.
弹性力学的一种边界无单元法   总被引:31,自引:7,他引:24  
程玉民  陈美娟 《力学学报》2003,35(2):181-186
首先对移动最小二乘副近法进行了研究,针对其容易形成病态方程的缺点,提出了以带权的正交函数作为基函数的方法-改进的移动最小二乘副近法,改进的移动最小二乘逼近法比原方法计算量小,精度高,且不会形成病态方程组,然后,将弹性力学的边界积分方程方法与改进的移动最小二乘逼近法结合,提出了弹性力学的一种边界无单元法,这种边界无单元法法是边界积分方程的无网格方法,与原有的边界积分方程的无网格方法相比,该方法直接采用节点变量的真实解为基本未知量,是边界积分方程无网格方法的直接解法,更容易引入界条件,且具有更高的精度,最后给出了弹性力学的边界无单元法的数值算例,并与原有的边界积分方程的无网格方法进行了较为详细的比较和讨论。  相似文献   

4.
用无网格局部Petrov-Galerkin法分析非线性地基梁   总被引:3,自引:1,他引:2  
龙述尧 《力学季刊》2002,23(4):547-551
利用无网格局部Petrov-Galerkin法求解了非线性地基梁。在Petrov-Galerkin方法中,采用移动最小二乘(MLS)近似函数作为场主量挠度的试函数并取移动最小二乘近似函数中的体验函数作为近似场函数的加权函数,采用罚因子法施加本质边界条件。文末给出了两个计算实例,算例的结果表明,Petrov-galerkin法不仅能成功地分析线性地基梁,而且也适用于解非线性地基梁,在分析非线性地基梁时具有收敛快,稳定性好的优点。  相似文献   

5.
基于Voronoi结构的无网格局部Petrov-Galerkin方法   总被引:26,自引:2,他引:24  
基于自然邻结点近似位移函数提出了一种用于求解弹性力学平面问题的无网格局部局部Petrov-Galerkin方法。这种方法在结构求解域Ω内任意布置离散的结点,并且利用需求结点的自然邻结点和Voronoi结构来构造整腐朽 求解的近似位移函数,对于构造好的近似位移函数,在局部Petrov-Galerkin方法建立整体求解的平控制方程,这样平衡方程的积分可在背景三角积分网格的形心上解析计算得到,而采用标准Galerkin方法的自然单元法需要三个数值积分点。该方法能够准确地施加边界条件,得到的系统矩阵是带状稀疏矩阵,对软件用户来说,这它学是一种安全的,真正的无网格方法,所得计算结果表明,该方法的计算精度与有限元四边界单元相当,但计算和形成系统平衡方程的时间比有限元法四边界单元提高了将近一倍,是一种理想的数值求解方法。  相似文献   

6.
配点类无网格法需要计算近似函数的二阶导数,因而在移动最小二乘(MLS)近似中至少要采用二次基函数。本文利用Voronoi图对双重点移动最小二乘近似法进行了改进,建立了基于Voronoi图的双重点移动最小二乘近似(VDG),并利用加权最小二乘法离散微分方程,导出了双重点最小二乘配点无网格法(MD GLS)。该方法将求解域用节点离散,并以节点为生成点建立Voronoi图,取Voronoi多边形的顶点为辅助点。近似函数及其二阶导数的计算过程可分解为两个步骤:首先用场函数节点值拟合辅助点处近似函数的一阶导数,再以辅助点处近似函数的一阶导数值拟合节点处近似函数的二阶导数。由于在每一步中只需计算MLS形函数及其一阶导数,这种近似方法需要较少的影响点和较小的影响域。同时借助于Voronoi结构的优良几何性质,可以快速地搜索影响点。研究表明,与基于MLS的加权最小二乘无网格法(MWLS)相比,这种方法可以显著提高计算效率,并且在精度和收敛性方面也有所改善。  相似文献   

7.
加权最小二乘无网格法是一种基于节点信息的纯无网格法,该方法使用最小二乘法建立系统的变分原理,通过移动最小二乘法构造近似函数,控制方程在节点处的残量使用最小二乘法予以消除,边界条件通过罚函数法引入。本文推导了瞬态热传导问题的加权最小二乘无网格计算格式,编制了相应的计算程序,算例结果表明,该方法具有精度高、前后处理简单的优点,是一种高效的的新型无网格法。  相似文献   

8.
基于改进的移动最小二乘(MLS)二阶导数近似,建立了一种求解弹性静力问题的无网格弱-强形式结合法(MLS-MWS)。该方法采用节点离散求解域,通过MLS构造形函数,将求解域划分为边界域和内部域,并分别使用控制方程的局部弱形式和强形式来建立离散系统方程。对强形式中涉及的近似函数二阶导数计算,提出了一种将其转化为求两次一阶导数的方法,与传统方法相比,该方法计算简单、精度高。MLS-MWS法结合了弱、强形式无网格法的优点,Neumann边界条件容易满足,并且只需在边界区域进行积分。文中应用该方法分析了两个弹性力学平面问题,分析结果表明本文方法具有良好的精度和收敛性。  相似文献   

9.
利用传统有限元法求解声压分布问题常常受到污染误差和色散误差的困扰.加权最小二乘无网格法(MWLS)是一种基于移动最小二乘(MLS)近似的无网格方法,求解声腔声压分布问题具有低色散、高精度的特点.然而传统的MLS近似有时容易产生病态矩阵,利用加权正交基函数构建改进的移动最小二乘(IMLS)近似,得到的系统方程为非病态的.论文基于改进的加权最小二乘无网格法(IMWLS)求解三维声腔内部声压分布.计算得到的声压分布和声压频响曲线都与参考值十分吻合,峰值误差和污染误差都比FEM的小,计算成本相比无单元伽辽金法显著降低.计算结果表明IMWLS相比传统的FEM,能在更高的频段内达到高精度,并且相比EFGM能大幅提高计算效率.  相似文献   

10.
复变量移动最小二乘法及其应用   总被引:9,自引:2,他引:7  
提出了复变量移动最小二乘法,并详细讨论了基于正交基函数的复变量移动最小二乘 法. 然后,将复变量移动最小二乘法和弹性力学的边界无单元法结合,提出了弹性力学的复 变量边界无单元法,推导了相应的公式,并给出了数值算例. 基于正交基函数的复变量移动 最小二乘法的优点是不形成病态方程组、精度高,所形成的无网格方法计算量小. 复变量边 界无单元法是边界积分方程的无网格方法的直接列式法,容易引入边界条件,且具有更高的 精度.  相似文献   

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

12.
The Burton-Miller boundary integral formulation is solved by a complex variable boundary element-free method (CVBEFM) for the boundary-only meshless analysis of acoustic problems with arbitrary wavenumbers. To regularize both strongly singular and hypersingular integrals and to avoid the computation of the solid angle and its normal derivative, a weakly singular Burton-Miller formulation is derived by considering the normal derivative of the solid angle and adopting the singularity subtraction procedures. To facilitate the implementation of the CVBEFM and the approximation of gradients of the boundary variables, a stabilized complex variable moving least-square approximation is selected in the meshless discretization procedure. The results show the accuracy and efficiency of the present CVBEFM and reveal that the method can produce satisfactory results for all wavenumbers, even for extremely large wavenumbers such as k = 10 000.  相似文献   

13.
将比例边界法与无单元伽辽金法相结合,建立了反平面断裂分析的无单元伽辽金比例边界法。这是一种边界型无网格法,在环向方向上采用无单元伽辽金法进行离散,因此计算时仅需要边界上的节点信息,不需要边界元所要求的基本解。为了便于施加本质边界条件,通过建立节点值和虚拟节点值之间的关系给出了修正的移动最小二乘形函数。在径向方向上,该方法利用解析的方法求解,因此是一种半解析的数值方法。最后,给出了数值算例,并验证了所提方法后处理简单和计算精度高的特点,适合于求解反平面断裂问题。  相似文献   

14.
A meshless approach to analysis of arbitrary Kirchhoff plates by the local boundary integral equation(LBIE) method is presented. The method combines the advantageous features of, all the three methods: the Galerkin finite element method (GFEM), the boundary element method (BEM) and the element-free Galerkin method (EFGM). It is a truly meshless method, which means that the discretization is independent of geometric subdivision into elements or cells, but is only based on a set of nodes (ordered or scattered) over a domain in question. It involves only boundary integration, however, over a local boundary centered at the node in question; It poses no difficulties in satisfying the essential boundary conditions while leading to banded and sparse system matrices using the moving least square (MLS) approximations. It is shown that high accuracy can be achieved for arbitrary geometries for clamped and simply-supported edge conditions. The method is found to be simple, efficient, and attractive. Project supported by the National Science Foundation of China (No. 19972019).  相似文献   

15.
A meshless local Petrov–Galerkin (MLPG) formulation is presented for bending problems of shear deformable shallow shells with orthotropic material properties. Shear deformation of shells described by the Reissner theory is considered. Analyses of shells under static and dynamic loads are given here. For transient elastodynamic case the Laplace-transform is used to eliminate the time dependence of the field variables. A weak formulation with a unit test function transforms the set of governing equations into local integral equations on local subdomains in the plane domain of the shell. Nodal points are randomly spread in that domain and each node is surrounded by a circular subdomain to which local integral equations are applied. The meshless approximation based on the moving least-squares (MLS) method is employed for the implementation. Unknown Laplace-transformed quantities are computed from the local boundary integral equations. The time-dependent values are obtained by the Stehfest’s inversion technique.  相似文献   

16.
In this paper, the local radial point interpolation meshless method (LRPIM) is used for the analysis of two‐dimensional potential flows, based on a local‐weighted residual method with the Heaviside step function as the weighting function over a local subdomain. Trial functions are constructed using radial basis functions. The present method is a truly meshless method based only on a number of randomly located nodes. Integration over the subdomains requires only a simple integration cell to obtain the solution. No element matrix assembly is required and no special treatment is needed to impose the essential boundary conditions. The novelty of the paper is the use of a local Heaviside weight function in the LRPIM, which does not need local domain integration and integrations only on the boundary of the local domains are needed. Effects of the sizes of local subdomain and interpolation domain on the performance of the present method are investigated. The behavior of shape parameters of multiquadrics has been systematically studied. Two numerical tests in groundwater and fluid flows are presented and compared with closed‐form solutions and finite element method. The results show that the use of a local Heaviside weight function in the LRPIM is highly accurate and possesses no numerical difficulties. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

17.
A new efficient meshless method based on the element-free Galerkin method is proposed to analyze the static deformation of thin and thick plate structures in this paper. Using the new 3D shell-like kinematics in analogy to the solid-shell concept of the finite element method, discretization is carried out by the nodes located on the upper and lower surfaces of the structures. The approximation of all unknown field variables is carried out by using the moving least squares (MLS) approximation scheme in the in-plane directions, while the linear interpolation is applied through the thickness direction. Thus, different boundary conditions are defined only using displacements and penalty method is used to enforce the essential boundary conditions. The constrained Galerkin weak form, which incorporates only displacement degrees of freedom (d.o.f.s), is derived. A modified 3D constitutive relationship is adopted in order to avoid or eliminate some self-locking effects. The numeric efficiency of the proposed meshless formulation is illustrated by the numeric examples.  相似文献   

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

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