The　behaviour　of　suction　anchors　subjected　to　combined　average　and　cyclic　loads　in　soft　clays　is　an　essential　consideration　in　their　design.　A　new　numerical　computation　method　was　developed　based　on　a　stiffness　degradation　model　for　soft　clay　that　was　proposed　by　performing　a　series　of　undrained　cyclic　triaxial　tests　and　embedded　in　the　ABAQUS　software　package　by　encoding　the　USDFLD　subroutine　and　defining　relevant　field　variables.　The　numerical　method　predicts　the　cyclic　behaviour　of　suction　anchors　by　coupling　the　elastic-perfectly　plastic　model　with　the　Mohr-Coulomb　yield　criterion　and　the　proposed　stiffness　degradation　model　that　can　reflect　the　stiffness　degradation　and　the　accumulation　of　plastic　deformation　of　soils　around　the　anchor　during　cyclic　loading.　The　numerical　method　was　verified　by　a　comparison　with　the　model　test　results　of　suction　anchors　subjected　to　combined　average　and　cyclic　loads　in　soft　clays.　The　proposed　method　can　predict　the　cyclic　deformation　as　well　as　the　bearing　capacity　and　capture　the　nonlinearity,　hysteresis　and　cyclic　accumulation　characteristics　of　the　load-displacement　responses　of　suction　anchors　in　soft　clays.