NASA NTRS2017
Landing vehicles impart thermal and strain energy onto the landing site from the retrorocket exhaust. Depending on the design of the vehicle, the energy may be great enough to cause spallation at the landing site. This damage may be minor and repairable in the case of landing on a terrestrial landing pad. For missions to other planetary bodies, the spallation may cause the landing site to become uneven and unstable, as well as damage. Simulating this phenomenon in a laboratory or computationally would require a significant amount of time and other resources. These resources typically are not available during the design phase of a mission. This paper presents a computationally-efficient model for the temperature and stress distributions that arise during landing. These quantities can be used along with existing failure criteria, such as the Hoek-Brown criterion for geological materials, to quickly determine whether spallation will occur. The stress and temperature distributions at the landing site are inherently 3D; however, there is a plane of symmetry and in that plane the distributions are 2D. Both quantities are modeled using series solutions to their governing partial differential equations (PDEs). The stress is modeled using the Airy stress potential function and its governing PDE is the biharmonic equation. The temperature is governed by Fourier's law. The models assume that stress due to gravity can be neglected, the points in the plane do not accelerate, and that the material properties are constant.