The seabed liquefaction under seismic loading is one of the main factors governing the overall stability of submarine pipelines.Based on Biot’s equations of dynamic consolidation for fully-saturated media

Mohr-Coulomb model and Drucker-Prager model were adopted to analyze the effective stress

and the values and regularity of excess pore water pressure in seabed along submarine pipelines were computed with ADINA

a general-purpose FEM analysis package

to estimate and analyze the momentary liquefactions of the seabed.In the numerical analysis

viscous-elastic artificial boundary was set up in the FEM model to simulate the transmission of seismic waves

which reproduced more realistically the dynamic response of the seabed along submarine pipelines under seismic loading.