Home About us Contact
|
Home About us Contact | ||
|
|
||
Integration Points (integration + point)
Selected AbstractsA reduced integration solid-shell finite element based on the EAS and the ANS concept,Geometrically linear problemsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 10 2009Marco Schwarze Abstract In this paper a new reduced integration eight-node solid-shell finite element is presented. The enhanced assumed strain (EAS) concept based on the Hu,Washizu variational principle requires only one EAS degree-of-freedom to cure volumetric and Poisson thickness locking. One key point of the derivation is the Taylor expansion of the inverse Jacobian with respect to the element center, which closely approximates the element shape and allows us to implement the assumed natural strain (ANS) concept to eliminate the curvature thickness and the transverse shear locking. The second crucial point is a combined Taylor expansion of the compatible strain with respect to the center of the element and the normal through the element center leading to an efficient and locking-free hourglass stabilization without rank deficiency. Hence, the element requires only a single integration point in the shell plane and at least two integration points in thickness direction. The formulation fulfills both the membrane and the bending patch test exactly, which has, to the authors' knowledge, not yet been achieved for reduced integration eight-node solid-shell elements in the literature. Owing to the three-dimensional modeling of the structure, fully three-dimensional material models can be implemented without additional assumptions. Copyright © 2009 John Wiley & Sons, Ltd. [source] The least-squares meshfree method for elasto-plasticity and its application to metal forming analysisINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 6 2005Kie-Chan Kwon Abstract A new meshfree method for the analysis of elasto-plastic deformation is presented. The method is based on the proposed first-order least-squares formulation for elasto-plasticity and the moving least-squares approximation. The least-squares formulation for classical elasto-plasticity and its extension to an incrementally objective formulation for finite deformation are proposed. In the formulation, equilibrium equation and flow rule are enforced in least-squares sense, i.e. their squared residuals are minimized, and hardening law and loading/unloading condition are enforced pointwise at each integration point. The closest point projection method for the integration of rate-form constitutive equation is inherently involved in the formulation, and thus the radial-return mapping algorithm is not performed explicitly. The proposed formulation is a mixed-type method since the residuals are represented in a form of first-order differential system using displacement and stress components as nodal unknowns. Also the penalty schemes for the enforcement of boundary and frictional contact conditions are devised and the reshaping of nodal supports is introduced to avoid the difficulties due to the severe local deformation near contact interface. The proposed method does not employ structure of extrinsic cells for any purpose. Through some numerical examples of metal forming processes, the validity and effectiveness of the method are discussed. Copyright © 2005 John Wiley & Sons, Ltd. [source] The Myoendothelial Junction: Breaking through the Matrix?MICROCIRCULATION, Issue 4 2009KATHERINE R. HEBERLEIN ABSTRACT Within the vasculature, specialized cellular extensions from endothelium (and sometimes smooth muscle) protrude through the extracellular matrix where they interact with the opposing cell type. These structures, termed myoendothelial junctions, have been cited as a possible key element in the control of several vascular physiologies and pathologies. This review will discuss observations that have led to a focus on the myoendothelial junction as a cellular integration point in the vasculature for both homeostatic and pathological conditions and as a possible independent signaling entity. We will also highlight the need for novel approaches to studying the myoendothelial junction in order to comprehend the cellular biology associated with this structure. [source] Second order homogenization method based on higher order finite elementsPROCEEDINGS IN APPLIED MATHEMATICS & MECHANICS, Issue 1 2005A. D?ster Modeling materials with lattice-like microstructures like open-cell foams requires an extended continuum mechanical setting on the macroscopic scale, e. g. a micropolar or micromorphic theory. In order to avoid the formulation of constitutive equations a higher order numerical homogenization scheme (FE2) is proposed. Therefore, each integration point possesses its own microstructure which, in the present case, consists of beam-like elements representing the cell walls. In this paper, the microstructures are discretized by continuum-based higher order locking free finite elements with high aspect ratios, leading to a numerically efficient treatment of a local displacement-driven boundary value problem according to the macroscopic strain and curvature. The resulting stress distributions in the microstructures are homogenized to macroscopic stresses and couple stresses. The approach is demonstrated by a numerical example. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim) [source] Optimal stress recovery points for higher-order bar elements by Prathap's best-fit methodINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 8 2009S. Rajendran Abstract Barlow was the first to propose a method to predict optimal stress recovery points in finite elements (FEs). Prathap proposed an alternative method that is based on the variational principle. The optimal points predicted by Prathap, called Prathap points in this paper, have been reported in the literature for linear, quadratic and cubic elements. Prathap points turn out to be the same as Barlow points for linear and quadratic bar elements but different for cubic bar element. Nevertheless, for all the three elements, Prathap points coincide with the reduced Gaussian integration points. In this paper, an alternative implementation of Prathap's best-fit method is used to compute Prathap points for higher-order (viz., 4,10th order) bar elements. The effectiveness of Prathap points as points of accurate stress recovery is verified by actual FE analysis for typical bar problems. Copyright © 2008 John Wiley & Sons, Ltd. [source] Exact integration of polynomial,exponential products with application to wave-based numerical methodsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 3 2009G. Gabard Abstract Wave-based numerical methods often require to integrate products of polynomials and exponentials. With quadrature methods, this task can be particularly expensive at high frequencies as large numbers of integration points are required. This paper presents a set of closed-form solutions for the integrals of polynomial,exponential products in two and three dimensions. These results apply to arbitrary polygons in two dimensions, and for arbitrary polygonal surfaces or polyhedral volumes in three dimensions. Quadrature methods are therefore not required for this class of integrals that can be evaluated quickly and exactly. Copyright © 2008 John Wiley & Sons, Ltd. [source] A new support integration scheme for the weakform in mesh-free methodsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 6 2010Yan Liu Abstract 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. [source] Higher-order XFEM for curved strong and weak discontinuitiesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 5 2010Kwok Wah Cheng Abstract The extended finite element method (XFEM) enables the accurate approximation of solutions with jumps or kinks within elements. Optimal convergence rates have frequently been achieved for linear elements and piecewise planar interfaces. Higher-order convergence for arbitrary curved interfaces relies on two major issues: (i) an accurate quadrature of the Galerkin weak form for the cut elements and (ii) a careful formulation of the enrichment, which should preclude any problems in the blending elements. For (i), we employ a strategy of subdividing the elements into subcells with only one curved side. Reference elements that are higher-order on only one side are then used to map the integration points to the real element. For (ii), we find that enrichments for strong discontinuities are easily extended to higher-order accuracy. In contrast, problems in blending elements may hinder optimal convergence for weak discontinuities. Different formulations are investigated, including the corrected XFEM. Numerical results for several test cases involving strong or weak curved discontinuities are presented. Quadratic and cubic approximations are investigated. Optimal convergence rates are achieved using the standard XFEM for the case of a strong discontinuity. Close-to-optimal convergence rates for the case of a weak discontinuity are achieved using the corrected XFEM. Copyright © 2009 John Wiley & Sons, Ltd. [source] A reduced integration solid-shell finite element based on the EAS and the ANS concept,Geometrically linear problemsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 10 2009Marco Schwarze Abstract In this paper a new reduced integration eight-node solid-shell finite element is presented. The enhanced assumed strain (EAS) concept based on the Hu,Washizu variational principle requires only one EAS degree-of-freedom to cure volumetric and Poisson thickness locking. One key point of the derivation is the Taylor expansion of the inverse Jacobian with respect to the element center, which closely approximates the element shape and allows us to implement the assumed natural strain (ANS) concept to eliminate the curvature thickness and the transverse shear locking. The second crucial point is a combined Taylor expansion of the compatible strain with respect to the center of the element and the normal through the element center leading to an efficient and locking-free hourglass stabilization without rank deficiency. Hence, the element requires only a single integration point in the shell plane and at least two integration points in thickness direction. The formulation fulfills both the membrane and the bending patch test exactly, which has, to the authors' knowledge, not yet been achieved for reduced integration eight-node solid-shell elements in the literature. Owing to the three-dimensional modeling of the structure, fully three-dimensional material models can be implemented without additional assumptions. Copyright © 2009 John Wiley & Sons, Ltd. [source] Adaptive through-thickness integration for accurate springback predictionINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 5 2008I. A. Burchitz Abstract Accurate numerical prediction of springback in sheet metal forming is essential for the automotive industry. Numerous factors influence the accuracy of prediction of this complex phenomenon by using the finite element method. One of them is the numerical integration through the thickness of shell elements. It is known that the traditional numerical schemes are very inefficient in elastic,plastic analysis and even for simple problems they require up to 50 integration points for an accurate springback prediction. An adaptive through-thickness integration strategy can be a good alternative. The main characteristic feature of the strategy is that it defines abscissas and weights depending on the integrand's properties and, thus, can adapt itself to improve the accuracy of integration. A concept of an adaptive through-thickness integration strategy for shell elements is presented in this paper. Its potential is demonstrated using two examples. Calculations of a simple test,bending a beam under tension,show that for a similar set of material and process parameters the adaptive rule with seven integration points performs significantly better than the traditional trapezoidal rule with 50 points. Simulations of an unconstrained cylindrical bending problem demonstrate that the adaptive through-thickness integration strategy for shell elements can guarantee an accurate springback prediction at minimal costs. Copyright © 2007 John Wiley & Sons, Ltd. [source] Development of a genetic algorithm-based lookup table approach for efficient numerical integration in the method of finite spheres with application to the solution of thin beam and plate problemsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 12 2006Suleiman BaniHani Abstract It is observed that for the solution of thin beam and plate problems using the meshfree method of finite spheres, Gaussian and adaptive quadrature schemes are computationally inefficient. In this paper, we develop a novel technique in which the integration points and weights are generated using genetic algorithms and stored in a lookup table using normalized coordinates as part of an offline computational step. During online computations, this lookup table is used much like a table of Gaussian integration points and weights in the finite element computations. This technique offers significant reduction of computational time without sacrificing accuracy. Example problems are solved which demonstrate the effectiveness of the procedure. Copyright © 2006 John Wiley & Sons, Ltd. [source] Symmetry preserving algorithm for large displacement frictionless contact by the pre-discretization penalty methodINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 15 2004D. Gabriel Abstract A three-dimensional contact algorithm based on the pre-discretization penalty method is presented. Using the pre-discretization formulation gives rise to contact searching performed at the surface Gaussian integration points. It is shown that the proposed method is consistent with the continuum formulation of the problem and allows an easy incorporation of higher-order elements with midside nodes to the analysis. Moreover, a symmetric treatment of mutually contacting surfaces is preserved even under large displacement increments. The proposed algorithm utilizes the BFGS method modified for constrained non-linear systems. The effectiveness of quadratic isoparametric elements in contact analysis is tested in terms of numerical examples verified by analytical solutions and experimental measurements. The symmetry of the algorithm is clearly manifested in the problem of impact of two elastic cylinders. Copyright © 2004 John Wiley & Sons, Ltd. [source] A numerical integration scheme for special finite elements for the Helmholtz equationINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 4 2003Peter Bettess Abstract The theory for integrating the element matrices for rectangular, triangular and quadrilateral finite elements for the solution of the Helmholtz equation for very short waves is presented. A numerical integration scheme is developed. Samples of Maple and Fortran code for the evaluation of integration abscissę and weights are made available. The results are compared with those obtained using large numbers of Gauss,Legendre integration points for a range of testing wave problems. The results demonstrate that the method gives correct results, which gives confidence in the procedures, and show that large savings in computation time can be achieved. Copyright © 2002 John Wiley & Sons, Ltd. [source] | ||||||||||