Traumatic injuries lead to articular cartilage lesion formation and result in the development of osteoarthritis. Recent research suggests that the early stage of mechanical injuries involve cell death (apoptosis and necrosis) and inflammation. In this thesis, we focus on building mathematical models to investigate the biological mechanism involving chondrocyte death and inflammatory responses in the process of cartilage degeneration. Chapter 1 describes the structure of articular cartilage, the process of carti- lage degeneration, and reviews of existing mathematical models. Chapter 2 presents a delay-diffusion-reaction model of cartilage lesion formation under cyclic loading. Computational methods were used to simulate the impact of varying loading stresses and erythropoietin levels. The model is parameterized with experimental results, and is therefore clinically relevant. Due to numerical limitations using delay differential equations, a new model is presented using tools for population dynamics. Chapter 3 presents an age and space-structured model of articular cartilage lesion formation un- der a single blunt impact. Age structure is introduced to represent the time delay in cytokine synthesis and cell transition. Numerical simulations produce similar tempo- ral and spatial patterns to our experimental data. In chapter 4, we extend our model under the cyclic loading setting. Chapter 5 builds a spatio-temporal model adapted from the former models, and investigates the distribution of model parameters using experimental data and statistical methods. Chapter 6 concludes.