Degenerative disc disease (DDD) is a common condition in elderly population that can be painful and can significantly affect individual's quality of life. Diagnosis of DDD allows prompt corrective actions but it is challenging due to the absence of any symptoms at early stages. In studying disc degeneration, measurement of the range of motion (RoM) and loads acting on the spine are crucial factors. However, direct measurement of RoM involves increased instrumentation and risk. In this paper, an innovative method is proposed for calculating RoM, emphasizing repeatability and reliability by considering the posterior thickness of the spine. This is achieved by offsetting the position of markers in relation to the actual vertebral loci. Three geometrically identical finite element models of L3-L4 are developed from a CT scan with different types of elements, and thereafter, mesh element-related metrics are provided for the assessment of the quality of models. The model with the best mesh quality is used for further analysis, where RoM are within ranges as reported in literature and in vivo experiment results. Various kinds of stresses acting on individual components including facet joints are analysed for normal and abnormal loading conditions. The results showed that the stresses in abnormal load conditions for all components including cortical (76.67 MPa), cancellous (69.18 MPa), annulus (6.30 MPa) and nucleus (0.343 MPa) are significantly greater as compared to normal loads (49.96 MPa, 44.2 MPa, 4.28 MPa and 0.23 MPa respectively). However, stress levels for both conditions are within safe limits (167-215 MPa for cortical, 46 MPa for the annulus and 3 MPa for facets). The results obtained could be used as a baseline motion and stresses of healthy subjects based on their respective lifestyles, which could benefit clinicians to suggest corrective actions for those affected by DDD.