The objective of this study is to integrate a continuum-based deformable tire and terrain interaction model into a general-purpose physics-based simulation environment capable of offroad vehicle mobility analysis and high-performance computing potential. Specifically, the physics-based deformable tire and terrain models which were recently proposed and validated by Yamashita, et al. will be implemented into the structure of the multi-physics simulation engine Chrono. In off-road vehicle mobility analysis, empirical and analytical models have been commonly used for vehicle-terrain interaction. While these models utilize experimental data or terramechanics theories to create quick predictive mobility models, they are unable to capture the highly nonlinear behavior of soft soil deformation, which can lead to inaccurate or unreliable results. In order to resolve these limitations, the use of physics-based numerical approaches have been proposed. These methods make use of finite element and discrete element simulations to describe the interaction between the vehicle and deformable terrain. Continuum-based finite element models transfer tire forces to the terrain and model the deformation with elasto-plastic constitutive equations. Discrete element soil uses a large number of small rigid body particles to describe the microscale behavior of granular terrain, with the deformation of the soil represented by the motion and contact of the particles. While these physics-based models offer a more accurate vehicle-terrain interaction model, the solution procedure can become complex and computationally expensive since co-simulation techniques are often used. To address these issues, the analysis of physics-based full vehicle dynamics simulations utilizing high-fidelity deformable tire and terrain models in a multi-physics engine with highperformance computing capability is desired. To this end, the deformable tire model formulated using the continuum mechanics based shear deformable laminated composite shell element proposed by Yamashita, et al. was integrated into the flexible body dynamics simulation framework of Chrono. This shell element is based on the absolute nodal coordinate formulation and is defined by the global position coordinates and the transverse gradient coordinates of its four nodes. Element lockings are eliminated with the incorporation of the enhanced assumed strain (EAS) and assumed natural strain approaches (ANS). The element formulation includes an extension to model laminated composite materials. Additionally, a locking-free 9-node brick element was integrated into the Chrono framework that makes use of the curvature coordinates at iv the center of the element. This element is formulated with the Hencky strain measure such that multiplicative finite strain plasticity theory can be used to incorporate soil plasticity models, such as the capped Drucker-Prager failure criterion. With the shear deformable laminated composite shell element and plastic soil brick element integrated into the Chrono multi-physi...