首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
    
Linearly conforming point interpolation method (LC‐PIM) is formulated for three‐dimensional elasticity problems. In this method, shape functions are generated using point interpolation method by adopting polynomial basis functions and local supporting nodes are selected based on the background cells. The shape functions so constructed have the Kronecker delta functions property and it allows straightforward imposition of point essential boundary conditions. Galerkin weak form is used for creating discretized system equations, and a nodal integration scheme with strain‐smoothing operation is used to perform the numerical integration. The present LC‐PIM can guarantee linear exactness and monotonic convergence for the numerical results. Numerical examples are used to examine the present method in terms of accuracy, convergence, and efficiency. Compared with the finite element method using linear elements, the LC‐PIM can achieve better efficiency, and higher accuracy especially for stresses. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

2.
    
A new support integration technique is proposed, which is similar to those used in truly mesh‐free methods. The contribution of this paper is to exploit the divergence‐free condition for the support integrals to construct quadrature formulas that only require three integration points per particle in two dimensions. Numerical examples show that the proposed integration method can achieve results that agree with manufactured closed‐form solutions. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

3.
    
This paper introduces a G space theory and a weakened weak form (W2) using the generalized gradient smoothing technique for a unified formulation of a wide class of compatible and incompatible methods. The W2 formulation works for both finite element method settings and mesh‐free settings, and W2 models can have special properties including softened behavior, upper bounds and ultra accuracy. Part I of this paper focuses on the theory and fundamentals for W2 formulations. A normed G space is first defined to include both continuous and discontinuous functions allowing the use of much more types of methods/techniques to create shape functions for numerical models. Important properties and a set of useful inequalities for G spaces are then proven in the theory and analyzed in detail. These properties ensure that a numerical method developed based on the W2 formulation will be spatially stable and convergent to the exact solutions, as long as the physical problem is well posed. The theory is applicable to any problems to which the standard weak formulation is applicable, and can offer numerical solutions with special properties including ‘close‐to‐exact’ stiffness, upper bounds and ultra accuracy. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

4.
    
Domain integration by Gauss quadrature in the Galerkin mesh‐free methods adds considerable complexity to solution procedures. Direct nodal integration, on the other hand, leads to a numerical instability due to under integration and vanishing derivatives of shape functions at the nodes. A strain smoothing stabilization for nodal integration is proposed to eliminate spatial instability in nodal integration. For convergence, an integration constraint (IC) is introduced as a necessary condition for a linear exactness in the mesh‐free Galerkin approximation. The gradient matrix of strain smoothing is shown to satisfy IC using a divergence theorem. No numerical control parameter is involved in the proposed strain smoothing stabilization. The numerical results show that the accuracy and convergent rates in the mesh‐free method with a direct nodal integration are improved considerably by the proposed stabilized conforming nodal integration method. It is also demonstrated that the Gauss integration method fails to meet IC in mesh‐free discretization. For this reason the proposed method provides even better accuracy than Gauss integration for Galerkin mesh‐free method as presented in several numerical examples. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

5.
    
A stabilized conforming (SC) nodal integration, which meets the integration constraint in the Galerkin mesh‐free approximation, is generalized for non‐linear problems. Using a Lagrangian discretization, the integration constraints for SC nodal integration are imposed in the undeformed configuration. This is accomplished by introducing a Lagrangian strain smoothing to the deformation gradient, and by performing a nodal integration in the undeformed configuration. The proposed method is independent to the path dependency of the materials. An assumed strain method is employed to formulate the discrete equilibrium equations, and the smoothed deformation gradient serves as the stabilization mechanism in the nodally integrated variational equation. Eigenvalue analysis demonstrated that the proposed strain smoothing provides a stabilization to the nodally integrated discrete equations. By employing Lagrangian shape functions, the computation of smoothed gradient matrix for deformation gradient is only necessary in the initial stage, and it can be stored and reused in the subsequent load steps. A significant gain in computational efficiency is achieved, as well as enhanced accuracy, in comparison with the mesh‐free solution using Gauss integration. The performance of the proposed method is shown to be quite robust in dealing with non‐uniform discretization. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

6.
    
Peridynamics is a non‐local mechanics theory that uses integral equations to include discontinuities directly in the constitutive equations. A three‐dimensional, state‐based peridynamics model has been developed previously for linearly elastic solids with a customizable Poisson's ratio. For plane stress and plane strain conditions, however, a two‐dimensional model is more efficient computationally. Here, such a two‐dimensional state‐based peridynamics model is presented. For verification, a 2D rectangular plate with a round hole in the middle is simulated under constant tensile stress. Dynamic relaxation and energy minimization methods are used to find the steady‐state solution. The model shows m‐convergence and δ‐convergence behaviors when m increases and δ decreases. Simulation results show a close quantitative matching of the displacement and stress obtained from the 2D peridynamics and a finite element model used for comparison. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

7.
    
The material point method (MPM) has demonstrated itself as a computationally effective particle method for solving solid mechanics problems involving large deformations and/or fragmentation of structures, which are sometimes problematic for finite element methods (FEMs). However, similar to most methods that employ mixed Lagrangian (particle) and Eulerian strategies, analysis of the method is not straightforward. The lack of an analysis framework for MPM, as is found in FEMs, makes it challenging to explain anomalies found in its employment and makes it difficult to propose methodology improvements with predictable outcomes. In this paper we present an analysis of the quadrature errors found in the computation of (material) internal force in MPM and use this analysis to direct proposed improvements. In particular, we demonstrate that lack of regularity in the grid functions used for representing the solution to the equations of motion can hamper spatial convergence of the method. We propose the use of a quadratic B‐spline basis for representing solutions on the grid, and we demonstrate computationally and explain theoretically why such a small change can have a significant impact on the reduction in the internal force quadrature error (and corresponding ‘grid crossing error’) often experienced when using MPM. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

8.
    
In this paper, issues regarding numerical integration of the discrete system of equations arising from natural neighbour (natural element) Galerkin methods are addressed. The sources of error in the traditional Delaunay triangle‐based numerical integration are investigated. Two alternative numerical integration schemes are analysed. First, a ‘local’ approach in which nodal shape function supports are exactly decomposed into triangles and circle segments is shown not to give accurate enough results. Second, a stabilized nodal quadrature scheme is shown to render high levels of accuracy, while resulting specially appropriate in a Natural Neighbour Galerkin approximation method. The paper is completed with several examples showing the performance of the proposed techniques. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

9.
    
This paper proposed a rotation‐free thin shell formulation with nodal integration for elastic–static, free vibration, and explicit dynamic analyses of structures using three‐node triangular cells and linear interpolation functions. The formulation is based on the classic Kirchhoff plate theory, in which only three translational displacements are treated as the filed variables. Based on each node, the integration domains are further formed, where the generalized gradient smoothing technique and Green divergence theorem that can relax the continuity requirement for trial function are used to construct the curvature filed. With the aid of strain smoothing operation and tensor transformation rule, the smoothed strains in the integration domain can be finally expressed by constants. The principle of virtual work is then used to establish the discretized system equations. The translational boundary conditions are imposed same as the practice of standard finite element method, while the rotational boundary conditions are constrained in the process of constructing the smoothed curvature filed. To test the performance of the present formulation, several numerical examples, including both benchmark problems and practical engineering cases, are studied. The results demonstrate that the present method possesses better accuracy and higher efficiency for both static and dynamic problems. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

10.
    
  相似文献   

11.
    
A new algorithm is developed to improve the accuracy and efficiency of the material point method for problems involving extremely large tensile deformations and rotations. In the proposed procedure, particle domains are convected with the material motion more accurately than in the generalized interpolation material point method. This feature is crucial to eliminate instability in extension, which is a common shortcoming of most particle methods. Also, a novel alternative set of grid basis functions is proposed for efficiently calculating nodal force and consistent mass integrals on the grid. Specifically, by taking advantage of initially parallelogram‐shaped particle domains, and treating the deformation gradient as constant over the particle domain, the convected particle domain is a reshaped parallelogram in the deformed configuration. Accordingly, an alternative grid basis function over the particle domain is constructed by a standard 4‐node finite element interpolation on the parallelogram. Effectiveness of the proposed modifications is demonstrated using several large deformation solid mechanics problems. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

12.
    
A superconvergent point interpolation method (SC‐PIM) is developed for mechanics problems by combining techniques of finite element method (FEM) and linearly conforming point interpolation method (LC‐PIM) using triangular mesh. In the SC‐PIM, point interpolation methods (PIM) are used for shape functions construction; and a strain field with a parameter α is assumed to be a linear combination of compatible stains and smoothed strains from LC‐PIM. We prove theoretically that SC‐PIM has a very nice bound property: the strain energy obtained from the SC‐PIM solution lies in between those from the compatible FEM solution and the LC‐PIM solution when the same mesh is used. We further provide a criterion for SC‐PIM to obtain upper and lower bound solutions. Intensive numerical studies are conducted to verify these theoretical results and show that (1) the upper and lower bound solutions can always be obtained using the present SC‐PIM; (2) there exists an αexact∈(0, 1) at which the SC‐PIM can produce the exact solution in the energy norm; (3) for any α∈(0, 1) the SC‐PIM solution is of superconvergence, and α=0 is an easy way to obtain a very accurate and superconvergent solution in both energy and displacement norms; (4) a procedure is devised to find a αprefer∈(0, 1) that produces a solution very close to the exact solution. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

13.
    
Stress‐point integration provides significant reductions in the computational effort of mesh‐free Galerkin methods by using fewer integration points, and thus facilitates the use of mesh‐free methods in applications where full integration would be prohibitively expensive. The influence of stress‐point integration on the convergence and stability properties of mesh‐free methods is studied. It is shown by numerical examples that for regular nodal arrangements, good rates of convergence can be achieved. For non‐uniform nodal arrangements, stress‐point integration is associated with a mild instability which is manifested by small oscillations. Addition of stabilization improves the rates of convergence significantly. The stability properties are investigated by an eigenvalue study of the Laplace operator. It is found that the eigenvalues of the stress‐point quadrature models are between those of full integration and nodal integration. Stabilized stress‐point integration is proposed in order to improve convergence and stability properties. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

14.
    
The meshless method is expected to become an effective procedure for realizing a CAD/CAE seamless system for analyses ranging from modelling to computation, because time‐consuming mesh generation processes are not required. In the present study, a new meshless approach, referred to as the Node‐By‐Node Meshless method is proposed, in which only nodal data is utilized to discretize the governing equations, which are derived using either the principle of virtual work or the Galerkin method. In this method, three key methodologies are utilized: (i) nodal integration using stabilization terms, (ii) interpolation by the Moving Least‐Squares Method, and (iii) a node‐by‐node iterative solver. This paper presents the formulation of the proposed method along with numerical results obtained for two‐dimensional elastostatic and eigenvalue problems. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

15.
    
A weighted least squares finite point method for compressible flow is formulated. Starting from a global cloud of points, local clouds are constructed using a Delaunay technique with a series of tests for the quality of the resulting approximations. The approximation factors for the gradient and the Laplacian of the resulting local clouds are used to derive an edge‐based solver that works with approximate Riemann solvers. The results obtained show accuracy comparable to equivalent mesh‐based finite volume or finite element techniques, making the present finite point method competitive. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

16.
The paper presents a fully meshless procedure fo solving partial differential equations. The approach termed generically the ‘finite point method’ is based on a weighted least square interpolation of point data and point collocation for evaluating the approximation integrals. Some examples showing the accuracy of the method for solution of adjoint and non-self adjoint equations typical of convective-diffusive transport and also to the analysis of compressible fluid mechanics problem are presented.  相似文献   

17.
    
Implementation of Dirichlet boundary conditions in mesh‐free methods is problematic. In Wagner and Liu (International Journal for Numerical Methods in Engineering 2001; 50 :507), a hierarchical enrichment technique is introduced that allows a simple implementation of the Dirichlet boundary conditions. In this paper, we provide some error analysis for the hierarchical enrichment mesh‐free technique. We derive optimal order error estimates for the hierarchical enrichment mesh‐free interpolants. For one‐dimensional elliptic boundary value problems, we can directly apply the interpolation error estimates to obtain error estimates for the mesh‐free solutions. For higher‐dimensional problems, derivation of error estimates for the mesh‐free solutions depends on the availability of an inverse inequality. Numerical examples in 1D and 2D are included showing the convergence behaviour of mesh‐free interpolants and mesh‐free solutions when the hierarchical enrichment mesh‐free technique is employed. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

18.
    
In part I of this paper, we have established the G space theory and fundamentals for W2 formulation. Part II focuses on the applications of the G space theory to formulate W2 models for solid mechanics problems. We first define a bilinear form, prove some of the important properties, and prove that the W2 formulation will be spatially stable, and convergent to exact solutions. We then present examples of some of the possible W2 models including the SFEM, NS‐FEM, ES‐FEM, NS‐PIM, ES‐PIM, and CS‐PIM. We show the major properties of these models: (1) they are variationally consistent in a conventional sense, if the solution is sought in a proper H space (compatible cases); (2) They pass the standard patch test when the solution is sought in a proper G space with discontinuous functions (incompatible cases); (3) the stiffness of the discretized model is reduced compared with the finite element method (FEM) model and possibly to the exact model, allowing us to obtain upper bound solutions with respect to both the FEM and the exact solutions and (4) the W2 models are less sensitive to the quality of the mesh, and triangular meshes can be used without any accuracy problems. These properties and theories have been confirmed numerically via examples solved using a number of W2 models including compatible and incompatible cases. We shall see that the G space theory and the W2 forms can formulate a variety of stable and convergent numerical methods with the FEM as one special case. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

19.
    
A new formulation of the element‐free Galerkin (EFG) method is presented in this paper. EFG has been extensively popularized in the literature in recent years due to its flexibility and high convergence rate in solving boundary value problems. However, accurate imposition of essential boundary conditions in the EFG method often presents difficulties because the Kronecker delta property, which is satisfied by finite element shape functions, does not necessarily hold for the EFG shape function. The proposed new formulation of EFG eliminates this shortcoming through the moving kriging (MK) interpolation. Two major properties of the MK interpolation: the Kronecker delta property (?I( s J)=δIJ) and the consistency property (∑In?I( x )=1 and ∑In?I( x )xIi=xi) are proved. Some preliminary numerical results are given. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

20.
    
The purpose of this paper is to propose a new quadrature formula for integrals with nearby singularities. In the boundary element method, the integrands of nearby singular boundary integrals vary drastically with the distance between the field and the source point. Especially, field variables and their derivatives at a field point near a boundary cannot be computed accurately. In the present paper a quadrature formula for ??‐isolated singularities near the integration interval, based on Lagrange interpolatory polynomials, is obtained. The error estimation of the proposed formula is also given. Quadrature formulas for regular and singular integrals with conjugate poles are derived. Numerical examples are given and the proposed quadrature rules present the expected polynomial accuracy. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

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

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