Abstract: Considering the influences of sea ice thermodynamic factors on the ice conditions,such as thickness,concentration,a sea ice thermodynamic-dynamic model is established based on the smooth particle hydrodynamics (SPH).In this sea ice model,the numerical diffusion in the traditional methods of the finite difference method (FDM) and the particle-in-cell (PIC) is avoided,while the ice edges can be modeled with high precision when the thermodynamics is introduced into the sea ice numerical model of SPH.In this paper,the basic theory of SPH and the thermodynamic factors of sea ice are presented first.The governing equations of the sea ice thermodynamics-dynamics are addressed by the using the thermodynamics to calculate the ice thickness and concentration in the SPH numerical simulation of the sea ice dynamics With this current sea ice model,the sea ice in Bohai Sea is simulated in 48 hours,and the distributions of sea ice thickness and drifting velocity in the whole ice field are obtained Moreover,the ice thickness,the temperature in vertical direction and the other thermodynamic factors are also discussed in detail.With the numerical results above,it can be concluded that this sea ice thermodynamic-dy-namic model based on the SPH method is an effective numerical model,and can be applied well in the sea ice simulation of Bohai Sea.