Computational & Technology Resources
an online resource for computational,
engineering & technology publications
PROCEEDINGS OF THE TENTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL STRUCTURES TECHNOLOGY
Edited by: B.H.V. Topping, J.M. Adam, F.J. Pallarés, R. Bru and M.L. Romero
The Finite Cell Method for Elasto-Plastic Problems
A. Abedian1, J. Parvizian2, A. Düster3, H. Khademyzadeh1 and E. Rank4
1Department of Mechanical Engineering, 2Department of Industrial Engineering,
A. Abedian, J. Parvizian, A. Düster, H. Khademyzadeh, E. Rank, "The Finite Cell Method for Elasto-Plastic Problems", in B.H.V. Topping, J.M. Adam, F.J. Pallarés, R. Bru, M.L. Romero, (Editors), "Proceedings of the Tenth International Conference on Computational Structures Technology", Civil-Comp Press, Stirlingshire, UK, Paper 310, 2010. doi:10.4203/ccp.93.310
Keywords: finite cell method, elasto-plasticity, numerical integration scheme, Voronoi polygons.
The finite cell method (FCM) [1,2] can be considered as an embedding or fictitious domain method combined with high-order finite elements . The FCM, as a high order finite element method for which only an integration mesh is necessary, is implemented in this paper to elasto-plastic problems. The Newton-Raphson iterative algorithm has been adapted by a scalar integration parameter suitable for the integration scheme. The potential of the FCM to solve materially nonlinear elasto-plastic problems is demonstrated in a simple plane strain example. In this embedding method, the mesh is not necessarily conforming to the boundaries. The boundary is extended to a simple domain which can be discretized using a simple mesh. Thus, the problem of discretization is replaced by a problem of integration. The classic Gauss quadrature integration scheme is modified to improve the convergence behavior of the method to produce more accurate results at a lower cost. In the integration procedure of the FCM, a test is normally needed to find out if the integration point is in the physical domain or in the extended domain. The efficiency of the method has been improved when replacing this test by checking if the Voronoi Polygon, associated with each integration point, is in the physical (material) domain, totally or partially. Then, the corresponding weight of each integration point is modified according to the ratio of the material area to the total area of the polygon. In a further attempt, the position of the integration point for the weak discontinuity problems is changed to the centroid of the physical part of the Voronoi polygon. These two modifications have improved the convergence behavior of the method. Converging to acceptable results, even for singular problems, when the mesh does not conform to the boundaries, and the shape functions are standard high order polynomials, is the key advantage of the finite cell method. Effort to enrich the approximation space is not necessary. This paper shows that the new integration method can provide accurate results for elasto-plastic problems too.
purchase the full-text of this paper (price £20)