Suggest to use Newton's method as follows:
NonlinearModelFit[dataS, model, {wp, a0, a1, a2, a3, a4, a5, b0, b1, b2, b3, b4, b5, r0, r1, r2, r3, r4, r5}, ω , Method -> "Newton"]
As I did not have the background on the constraints you are using, I left them out. Simple enough to include them again.