Mapping, Shape Functions, and Numerical Integration

Author(s):  
Darius Olesch ◽  
Charlotte Kuhn ◽  
Alexander Schlüter ◽  
Ralf Müller

AbstractPhase field models for fracture are energy-based and employ a continuous field variable, the phase field, to indicate cracks. The width of the transition zone of this field variable between damaged and intact regions is controlled by a regularization parameter. Narrow transition zones are required for a good approximation of the fracture energy which involves steep gradients of the phase field. This demands a high mesh density in finite element simulations if 4-node elements with standard bilinear shape functions are used. In order to improve the quality of the results with coarser meshes, exponential shape functions derived from the analytic solution of the 1D model are introduced for the discretization of the phase field variable. Compared to the bilinear shape functions these special shape functions allow for a better approximation of the fracture field. Unfortunately, lower-order Gauss-Legendre quadrature schemes, which are sufficiently accurate for the integration of bilinear shape functions, are not sufficient for an accurate integration of the exponential shape functions. Therefore in this work, the numerical accuracy of higher-order Gauss-Legendre formulas and a double exponential formula for numerical integration is analyzed.


Biomechanics ◽  
2012 ◽  
pp. 295-312
Author(s):  
Cees Oomens ◽  
Marcel Brekelmans ◽  
Frank Baaijens

Author(s):  
Márton Petö ◽  
Fabian Duvigneau ◽  
Daniel Juhre ◽  
Sascha Eisenträger

Abstract Polygonal finite elements offer an increased freedom in terms of mesh generation at the price of more complex, often rational, shape functions. Thus, the numerical integration of rational interpolants over polygonal domains is one of the challenges that needs to be solved. If, additionally, strong discontinuities are present in the integrand, e.g., when employing fictitious domain methods, special integration procedures must be developed. Therefore, we propose to extend the conventional quadtree-decomposition-based integration approach by image compression techniques. In this context, our focus is on unfitted polygonal elements using Wachspress shape functions. In order to assess the performance of the novel integration scheme, we investigate the integration error and the compression rate being related to the reduction in integration points. To this end, the area and the stiffness matrix of a single element are computed using different formulations of the shape functions, i.e., global and local, and partitioning schemes. Finally, the performance of the proposed integration scheme is evaluated by investigating two problems of linear elasticity.


1966 ◽  
Vol 25 ◽  
pp. 227-229 ◽  
Author(s):  
D. Brouwer

The paper presents a summary of the results obtained by C. J. Cohen and E. C. Hubbard, who established by numerical integration that a resonance relation exists between the orbits of Neptune and Pluto. The problem may be explored further by approximating the motion of Pluto by that of a particle with negligible mass in the three-dimensional (circular) restricted problem. The mass of Pluto and the eccentricity of Neptune's orbit are ignored in this approximation. Significant features of the problem appear to be the presence of two critical arguments and the possibility that the orbit may be related to a periodic orbit of the third kind.


Sign in / Sign up

Export Citation Format

Share Document