This paper presents the mathematical model consisting of conservation and balance laws (CBL) of classical continuum mechanics (CCM) and the constitutive theories derived using entropy inequality and representation theorem for thermoviscoelastic solids (TVES) matter without memory. The CBL and the constitutive theories take into account finite deformation and finite strain deformation physics. This mathematical model is thermodynamically and mathematically consistent and is ideally suited to study nonlinear dynamics of TVES and dynamic bifurcation and is used in the work presented in this paper. The finite element formulations are constructed for obtaining the solution of the initial value problems (IVPs) described by the mathematical models. Both space-time coupled as well as space-time decoupled finite element methods are considered for obtaining solutions of the IVPs. Space-time coupled finite element formulations based on space-time residual functional (STRF) that yield space-time variationally consistent space-time integral forms are considered. This approach ensures unconditional stability of the computations during the entire evolution. In the space-time decoupled finite element method based on Galerkin method with weak form for spatial discretization, the solutions of nonlinear ODEs in time resulting from the decoupling of space and time are obtained using Newmark linear acceleration method. Newton's linear method is used to obtain converged solution for the nonlinear system of algebraic equations at each time step in the Newmark method. The different aspects of the deformation physics leading to the factors that influence nonlinear dynamic response and dynamic bifurcation are established using the proposed mathematical model, the solution method and their validity is demonstrated through model problem studies presented in this paper. Energy methods and superposition techniques in any form including those