本发明提供了一种考虑土体空间变异性的管道撞击损伤分析方法、设备，首先，根据海床土体尺寸和海床土体强度参数统计特征，将土体离散成一系列欧拉实体单元；通过Python脚本读取每个欧拉实体单元的形心坐标；然后生成三维土体强度平稳随机场，根据初始温度场进行叠加，进而得到三维土体强度非平稳随机场，通过子程序技术与大变形有限元数值模型进行耦合分析，得到N组管道凹痕变形值，然后进行分布拟合，得到管道损伤分析结果。本发明的有益效果是：针对大变形有限元数值模型中欧拉单元在空间保持固定不变的特点，通过设置温度场及子程序技术有效捕捉土体在计算过程中的位置，实现土体强度随空间变化的特征，分析得到的管道损伤结果与实际情况更接近。