The use of close-fitting PPE is essential to prevent exposure to dispersed airborne matter, including the COVID-19 virus. The current pandemic has increased pressure on the healthcare system, leading to medical professionals using high-grade PPE for prolonged times, resulting in device-insduced skin injuries. This study focuses on computationally improving the interaction between skin and PPE to reduce the likelihood of discomfort and tissue damage. A finite element model is developed to simulate the movement of PPE against the face during day-to-day tasks. Due to limited available data on skin characteristics and how these vary interpersonally between sexes, races and ages, the main objective of this study was to establish the effects and trends that mask modifications have on the resulting subsurface strain energy density distribution in the skin. These modifications include the material, geometry and interfacial properties. Overall, the results show that skin injury can be reduced by using softer mask materials, whilst friction against the skin should be minimised, e.g. through use of micro-textures, humidity control and topical creams. Furthermore, the contact area between the mask and skin should be maximised, whilst the use of soft materials with incompressible behaviours (e.g. most elastomers) should be avoided.