Numerical Solution of Coupled Thermo-Elastic-Plastic Dynamic Problems

The article considers a numerical method for solving a two-dimensional coupled dynamic thermoplastic boundary value problem based on deformation theory of plasticity. Discrete equations are compiled by the finite-difference method in the form of explicit and implicit schemes. The solution of the explicit schemes is reduced to the recurrence relations regarding the components of displacement and temperature. Implicit schemes are efficiently solved using the elimination method for systems with a three diagonal matrix along the appropriate directions. In this case, the diagonal predominance of the transition matrices ensures the convergence of implicit difference schemes. The problem of a thermoplastic rectangle clamped from all sides under the action of an internal thermal field is solved numerically. The stress-strain state of a thermoplastic rectangle and the distribution of displacement and temperature over various sections and points in time have been investigated.

Stability of discrete schemes of Biot’s poroelastic equations

Summary The efficient and accurate numerical modeling of Biot’s equations of poroelasticity requires the knowledge of the exact stability conditions for a given set of input parameters. Up to now, a numerical stability analysis of the discretized elastodynamic Biot’s equations has been performed only for a few numerical schemes. We perform the von Neumann stability analysis of the discretized Biot’s equations. We use an explicit scheme for the wave propagation and different implicit and explicit schemes for Darcy’s flux. We derive the exact stability conditions for all the considered schemes. The obtained stability conditions for the discretized Biot’s equations were verified numerically in one-, two- and three-dimensions. Additionally, we present von Neumann stability analysis of the discretized linear damped wave equation considering different implicit and explicit schemes. We provide both the Matlab and symbolic Maple routines for the full reproducibility of the presented results. The routines can be used to obtain exact stability conditions for any given set of input material and numerical parameters.