The evaluation of the bearing capacity of strip footings generally assumes that the soil is either dry or fully saturated, which contradicts the actual condition in nature where the soil is often in a partially saturated state. Furthermore, infiltration has a significant impact on the shear strength of the soil. Following the upper bound theory of the limit analysis, this article provides a theoretical framework for assessing the bearing capacity under transient flow with linear variation in infiltration intensity for the first time. Firstly, the closed form of suction stress under linear transient infiltration is derived using Laplace transform and introduced into the Mohr–Coulomb criterion. A discrete failure mechanism with fewer variables and higher accuracy is provided to ensure kinematic admissibility. The upper bound solution for bearing capacity is obtained by solving the power balance equation. The present results are compared with results from the published literature and the finite element, confirming the validity and superiority of the theoretical framework provided. A parametric analysis is also conducted on three hypothetical soil types (fine sand, silt, and clay), and the results show that unsaturated transient infiltration has a positive influence on increasing the foundation bearing capacity. The magnitude of the influence is comprehensively controlled by factors such as soil type, saturated hydraulic conductivity, infiltration intensity, infiltration time, and water table depth. The increase in bearing capacity due to unsaturated transient infiltration can be incorporated into Terzaghi’s equation as a separate component presented in tabular form for engineering design purposes.

In engineering mathematics, the unsaturated nature of soil has a significant impact on the seismic bearing capacity solution. However, it has generally been neglected in the published literature to date. Based on the kinematic approach of limit analysis, the present study proposes a method for calculating the bearing capacity of shallow strip footings located in unsaturated soils, taking four common types of soils as examples. The modified pseudo-dynamic (MPD) approach is used to calculate the seismic forces varying with time and space, and the layerwise summation method is used to derive the power generated by the seismic forces. In the calculation of internal energy dissipation, this paper introduces the effective stress based on the suction stress to derive the cohesion expression at different depths. The analytical formula of bearing capacity is obtained by the principle of virtual work, and its value is optimized by the Sequential Quadratic Programming (SQP) algorithm. In order to verify the validity of the proposed method, the present results are compared with the solutions published so far and a good agreement is obtained. Finally, a parametric study is performed to investigate the influence of different types of parameters on the bearing capacity.