The load transfer method is important for the settlement prediction of axially loaded piles, but in multi-layered complex soils, it lacks analytical solutions. Traditional numerical methods such as the finite element method suffer from strong dependence on mesh generation, time-consuming iterative calculations, and high computational costs for back-analysis. This paper proposes a load transfer analysis model based on a Domain Decomposition Physics-Informed Neural Network. A multi-subnet parallel architecture is adopted to simulate multi-layered soils, solving the problem of inter-layer stress–strain discontinuity through interface coupling and gradient continuity constraints; a non-dimensionalization system and a hard constraint mechanism are introduced to enhance training efficiency and physical consistency; and a two-stage analysis framework comprising surrogate model forward analysis and field data inversion is established. Numerical experimental results indicate that the forward analysis of this model is in high agreement with FEM simulation results, and computational efficiency is improved by six orders of magnitude; based on a small amount of field static load test data, multi-layer soil parameters are accurately inverted, achieving more precise pile settlement prediction than FEM. Comparative analysis validates the effectiveness of the domain decomposition multi-subnet over a single network, demonstrating extensibility to hyperbolic and exponential multi-soil constitutive models.