Sometimes a pavement deflects only because of seasonal volume change of expansive soils in the pavement subgrade. Engineering practitioners expect an implementable and straightforward analysis method for a geosynthetic-reinforced pavement subjected to the swelling/shrinkage issue of expansive clayey subgrade, in an effort to find the bending moment, shear force and tension force distributions through the reinforced pavement, which are induced from the volumetric changes of subgrade soils. The virtual load method (VLM) was proposed in the past following the Timoshenko beam theory to analyze geosynthetic-reinforced pavement on expansive soils. In the VLM, the unknown virtual distributed load was obtained in the way by applying the inverse theory for the identification of material parameters of the pavement-foundation system. It was seen that the selection of the number of material parameters to obtain virtual load significantly affects the accuracy of the structural properties of the pavement and tensile properties of the geosynthetics if the linear least square method is used. In this paper, a unique numerical scheme was proposed in the hopes of solving the issue. After a forward problem was solved numerically, the Timoshenko beam deflection was taken as a start-up for the inverse problem to back analyze the load applied to the pavement. Solutions from forward/backward examples have explicitly shown the accuracy achieved related to the bending moment in the pavement and tension in the geosynthetic reinforcements. The proposed methodology can be applied for an in-depth understanding of the geosynthetic function for the mitigation of longitudinal cracks on pavements caused by heave/shrinkage of expansive soils.