Abstract
A numerical method for calculating the three-dimensional processes of impact interaction of elastoplastic bodies under large displacements and deformations based on the multi mesh sharp interface method and modified Godunov scheme is presented. To integrate the equations of dynamics of an elastoplastic medium, the principle of splitting in space and in physical processes is used. The solutions of the Riemann problem for first and second order accuracy for compact stencil for an elastic medium in the case of an arbitrary stress state are obtained and presented, which are used at the “predictor” step of the Godunov scheme. A modification of the scheme is described that allows one to obtain solutions in smoothness domains with a second order of accuracy on a compact stencil for moving Eulerian-Lagrangian grids. Modification is performed by converging the areas of influence of the differential and difference problems for the Riemann’s solver. The “corrector” step remains unchanged for both the first and second order accuracy schemes. Three types of difference grids are used. The first – a moving surface grid – consists of a continuous set of triangles that limit and accompany the movement of bodies; the size and number of triangles in the process of deformation and movement of the body can change. The second – a regular fixed Eulerian grid – is limited to a surface grid; separately built for each body; integration of equations takes place on this grid; the number of cells in this grid can change as the body moves. The third grid is a set of local Eulerian-Lagrangian grids attached to each moving triangle of the surface from the side of the bodies and allowing obtain the parameters on the boundary and contact surfaces. The values of the underdetermined parameters in cell’s centers near the contact boundaries on all types of grids are interpolated. Comparison of the obtained solutions with the known solutions by the Eulerian-Lagrangian and Lagrangian methods, as well as with experimental data, shows the efficiency and sufficient accuracy of the presented three-dimensional methodology.
Publisher
Uniscience Publishers LLC
Reference29 articles.
1. Godunov S.K., Zabrodin A.V., Prokopov G.P. (1961) The difference scheme for two-dimensional unsteady problems of gas dynamics and the calculation of the flow around a departing shock wave // Journal of Computational Mathematics and Mathematical Physics, 1961, 1(6): 1020-1050.
2. Godunov S.K., Zabrodin A.V., Ivanov M.Ya. et al. (1976) Numerical solution of multidimensional problems of gas dynamics. M .: Nauka, 1976.
3. S.L. Gavrilyuk, N. Favrie, R. Saurel (20080 Modeling wave dynamics of compressible elastic materials, J. Comput. Phys. 227 (5): 2941–2969.
4. Cheban V. G., Russu I. V (1976) Numerical methods for solving problems of the dynamic theory of elasticity. Chisinau: Stinnitsa, 1976.
5. P.T. Barton, D. Drikakis, E. Romenski, V.A. Titare (2009) Exact and approximate solutions of Riemann problems in non-linear elasticity, J. Comput. Phys. 228(18): 7046– 7068.