A numerical study to a generalized Korteweg-de Vries (KdV) equation is adopted to model the propagation and disintegration of large-amplitude internal solitary waves (ISWs) in the South China Sea (SCS). Based on theoretical analysis and in situ measurements, the drag coefficient of the Chezy friction is regarded as inversely proportional to the initial amplitude of an ISW, rather than a constant as assumed in the previous studies. Numerical simulations of ISWs propagating from a deep basin to a continental shelf are performed with the generalized KdV model. It is found that the depression waves are disintegrated into several solitons on the continental shelf due to the variable topography. It turns out that the amplitude of the leading ISW reaches a maximum at the shelf break, which is consistent with the field observation in the SCS. Moreover, a dimensionless parameter defining the relative importance of the variable topography and friction is presented.