Based on the Snyder-Mitchell linear model and the cross-spectral density (CSD) function, the analytical propagation formula of twisted Gaussian Schell-model (TGSM) beams in strongly nonlocal nonlinear medium (SNNM) is derived. Then the propagation characteristics of TGSM beam are studied. It is found that the soliton radius is jointly determined by the initial power, coherence length, and twist factor; the degree of spatial coherence is adjusted by changing the twist factor without affecting the soliton intensity. In the case of non-soliton properties, there is a threshold of coherence length which makes partially coherent beams have the same evolution law as completely coherent beams. Furthermore, increasing the twist factor, decreasing the coherence length and initial power can improve the beam quality of the beam propagating in SNNM. © 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement.