The snake-laying method is widely employed as an effective strategy for global buckling mitigation in submarine pipelines. The uneven distribution of soil resistance along pipeline routes significantly amplifies the complexity of global buckling responses in snake-laid pipelines and challenges their control mechanisms. This study establishes a finite element computational model to investigate the effects of soil resistance dis...