研究具有流变特性的由一个Maxwell体和任意个Kelvin体串联组成的广义Kelvin模型，并对其本构关系进行推导。在Microsoft Visual Studio平台上采用C++语言编译模型的头文件和源文件，开发出广义Kelvin模型的动态链接库(DLL)文件，在二维颗粒流程序(PFC2D)中对DLL文件进行加载调用，实现对广义Kelvin模型的开发。对含不同数量Kelvin体的广义Kelvin模型进行应力松弛模拟试验，得到接触力随时间的变化关系并与理论解对比，验证了开发模型的准确性。在此基础上，以某露天煤矿的含粉质黏土边坡为工程实例，利用开发的广义Kelvin模型对其失稳过程进行模拟，所得结果与边坡现场特殊的帽檐形破坏模式相符，而与采用普通线弹性模型模拟得到的边坡失稳形态为“圆弧滑动”不一致，说明该模型可用于研究此类含粉质黏土边坡失稳的破坏机制，并在此基础上，提出了对该类含具有流变特性的粉质黏土边坡的处治建议，研究结论可为今后类似研究提供有益的参考。

A constitutive relationship of Kelvin model’s type is developed by sequentially combing a Maxwell component and any number of Kelvin components. Through compiling the head and source files of generalized Kelvin model by using C++ on Microsoft Visual Studio platform, the dynamic link library (DLL) file of the generalized Kelvin constitutive model is developed, and a user-defined generalized Kelvin model is implemented by calling the DLL file in two-dimensional Particle Flow Code (PFC2D). Through stress relaxation tests, in which different numbers of Kelvin models are included in the generalized Kelvin model, the relationship between contact force and time is obtained, which turns out to coincide with the analytical solution. Based on these research results, the programmed generalized Kelvin model is used to simulate the failure process of a silty clay slope of an open pit coal mine. The numerical results illustrate that the failure surface is brim-shaped, whereas the elastic model does not predict such a failure mode, and instead it yields a circular sliding surface, implying that the developed generalized Kelvin model can be used to study the failure mechanism of silty clay slope.