Numerical implementation of iteration processes applied in the study of nonlinear boundary value problems of the theory of elasticity and filtration.

Autor: Narbut, M.
Zdroj: Vestnik St. Petersburg University: Mathematics; Jun2008, Vol. 41 Issue 2, p161-166, 6p
Abstrakt: Nonlinear elastic problems for hardening media are solved by applying the universal iteration process proposed by A.I. Koshelev in his works on the regularity of solutions to quasilinear elliptic and parabolic systems. This requires numerically solving a linear elliptic system at each step of the iteration procedure. The method is numerically implemented in the MATLAB environment by using its PDE Toolbox. A modification of the finite-element procedure is suggested in order to solve a linear algebraic system at each iteration step. The computer model is tested on simple examples. The same nonlinear problems are also solved by the method of elastic solutions, which consists in replacing the Laplace operator in the universal iteration process by the Lamé operator of linear elasticity. As is known, this iteration process converges to a weak solution of the nonlinear problem, provided that the displacements are fixed on the boundary. The method is tested on examples with stresses on the boundary. The third part of the paper is devoted to the nonlinear filtration problem. General properties of the iteration process for nonlinear parabolic systems have been studied by A.I. Koshelev and V.M. Chistyakov. The numerical implementation is based on slightly modified PDE Toolbox procedures. The program is tested on simple examples. [ABSTRACT FROM AUTHOR]
Databáze: Complementary Index