En matemáticas de ingeniería, la naturaleza no saturada del suelo tiene un impacto significativo en la solución de la capacidad de carga sísmica. Sin embargo, generalmente ha sido ignorada en la literatura publicada hasta la fecha. Basado en el enfoque cinemático del análisis límite, el presente estudio propone un método para calcular la capacidad de carga de cimentaciones superficiales ubicadas en suelos no saturados, tomando como ejemplos cuatro tipos comunes de suelos. El enfoque pseudo-dinámico modificado (MPD) se utiliza para calcular las fuerzas sísmicas que varían con el tiempo y el espacio, y se utiliza el método de suma por capas para derivar la potencia generada por las fuerzas sísmicas. En el cálculo de la disipación de energía interna, este artículo introduce el esfuerzo efectivo basado en el esfuerzo de succión para derivar la expresión de cohesión a diferentes profundidades. La fórmula analítica de la capacidad de carga se obtiene mediante el principio de trabajo virtual, y su valor se optimiza mediante el algoritmo de Programación Cuadrática Secuencial (SQP). Para verificar la validez del método propuesto, los resultados actuales se comparan con las soluciones publicadas hasta ahora y se obtiene un buen acuerdo. Finalmente, se realiza un estudio paramétrico para investigar la influencia de diferentes tipos de parámetros en la capacidad de carga.