Contact　stiffness　and　backlash　in　the　harmonic　drive　significantly　impact　a　robot＇s　positioning　accuracy　and　vibration　characteristics.　The　height　of　the　harmonic　drive　tooth　pair　is　typically　less　than　1　mm,　making　the　measurement　and　modeling　of　backlash　and　contact　stiffness　inherently　complex.　This　paper　proposes　a　contact　stiffness　and　backlash　model　by　establishing　a　correlation　between　fractal　parameters　and　tooth　contact　load.　To　obtain　the　fractal　roughness　parameters　of　the　real　machined　tooth　surface,　a　combination　of　a　noncontact　optical　profiler　and　the　RMS　method　is　employed.　Subsequently,　the　study　explores　the　influence　of　rough　tooth　surface　and　contact　force　fractal　parameters　on　contact　stiffness　and　gear　backlash.　The　results　demonstrate　the　substantial　impact　of　surface　topography　parameters　and　contact　force　on　contact　stiffness　and　backlash.　Specifically,　an　increase　in　the　fractal　dimension　correlates　with　a　reduction　in　gear　backlash　and　contact　stiffness.　Conversely,　the　fractal　roughness　parameter　exhibits　the　opposite　effect.　Notably,　an　increase　in　contact　force　enhances　contact　stiffness.