This paper introduces a SIR model with a nonlinear incidence rate and incorporates a vaccination scenario, offering a more realistic framework for disease dynamics. Assuming a vaccination coverage of $p\%$, we analyze its influence on epidemic outcomes. The model features a disease-free equilibrium $E_0$ from which we derive the basic reproduction number $\mathcal{R}_0$, serving as the threshold for disease eradication. When $\mathcal{R}_0 > 1$, an endemic equilibrium $E_1$ emerges; conversely, $\mathcal{R}_0< 1$ guarantees the global stability of $E_0$, indicating disease elimination. A transcritical bifurcation at $\mathcal{R}_0=1$ captures the transition between disease extinction and persistence, with no evidence of Hopf bifurcations as shown by limit set analysis. Sensitivity analysis of $\mathcal{R}_0$ highlights key parameters influencing transmission, informing intervention strategies. We also develop an optimal control framework to determine the most effective vaccination coverage, providing actionable insights for public health policies. Numerical simulations validate the theoretical results, illustrating how variations in $p$ impact outbreak trajectories and underscoring the importance of sustained vaccination efforts. By integrating nonlinear transmission with vaccination dynamics, this study advances epidemic modeling and offers practical tools for disease management.
[1] I. Barkana, Barbalat’s lemma and stability- misuse of a correct mathematical result?, Mathematics in Engineering, Science and Aerospace (MESA), 7(1) (2016), 197–218.
[2] D. Bernoulli, An attempt at a new analysis of the mortality caused by smallpox and of the advantages of inoculation to prevent it, Rev. Med. Virol., 14 (2004), 275–288.
[3] W. Cui and Y. Zhao, Saddle-node bifurcation and Bogdanov-Takens bifurcation of a SIRS epidemic model with nonlinear incidence rate, J. Differential Equations, 384 (2024), 252–278.
[4] R. Din, K. Shah, M. A. Alqudah, T. Abdeljawad, and F. Jarad, Mathematical study of SIR epidemic model under convex incidence rate, AIMS Math., 5(6) (2020), 7548–7561.
[5] W. O. Kermack and A. G. McKendrick, A contribution to the mathematical theory of epidemics, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character, 115(772) (1927), 700–721.
[6] W. O. Kermack and A. G. McKendrick, Contributions to the mathematical theory of epidemics-I, Bull. Math. Biol., 53 (1991), 33–55.
[7] W. O. Kermack and A. G. McKendrick, Contributions to the mathematical theory of epidemics–II, The problem of endemicity, 1932, Bull. Math. Biol., 53 (1991), 57–87.
[8] W. O. Kermack and A. G. McKendrick, Contributions to the mathematical theory of epidemics–III, Further studies of the problem of endemicity, 1933, Bull. Math. Biol., 53 (1991), 89–118.
[9] T. Khan, R. Ullah, and G. Zaman, Hepatitis B virus transmission via epidemic model, Advances in Epidemio-logical Modeling and Control of Viruses, (2023), 29–54.
[10] J. Li, Z. Teng, G. Wang, L. Zhang, and C. Hu, Stability and bifurcation analysis of an SIR epidemic model with logistic growth and saturated treatment, Chaos, Solitons and Fractals, 99 (2017), 63–71.
[11] E. M. Lotfi, M. Maziane, K. Hattaf, and N. Yousfi, Partial differential equations of an epidemic model with spatial diffusion, International Journal of Partial Differential Equations, 2014(1) (2014), 1–6.
[12] M. Martcheva, An introduction to mathematical epidemiology, New York: Springer, 61 (2015), 9–31.
[13] J. P. Mauricio de Carvalho and A. A. Rodrigues, SIR model with vaccination: bifurcation analysis, Qual. Theory Dyn. Syst., 22(3) (2023), 32.
[14] L. Perko, Differential equations and dynamical systems. Springer Science and Business Media, 7 (2013).
[15] A. Ramponi and M. E. Tessitore, Optimal social and vaccination control in the SVIR epidemic model, Mathematics, 12(7) (2024), 933.
[16] R. E. Russell, R. A. Katz, K. L. Richgels, D. P. Walsh, and E. H. Grant, A framework for modeling emerging diseases to inform management, Emerg. Infect. Dis., 23(1) (2017), 1–6.
[17] M. Shakeri and R. Kazemi, Analysis of an SIR model with nonlinear incidence rate and the full impact of vaccination (in Persian), Journal of Advanced Mathematical Modeling, 14(2) (2024), 96–108.
[18] K. Shigemoto, Various Logistic Curves in SIS and SIR Models, arXiv preprint arXiv. (2022), 2211.04186.
[19] N. Yousfi, K. Hattaf, and A. Tridane, Modeling the adaptive immune response in HBV infection, Journal of Mathematical Biology, 63 (2011), 933–957.
[20] F. Zhang, W. Cui, Y. Dai, and Y. Zhao, Bifurcations of an SIRS epidemic model with a general saturated incidence rate, Math. Biosci. Eng., 19(11) (2022), 10710–10730.
[21] Y. Zhao, W. Li, and Y. Wang, Global stability of epidemic models under discontinuous treatment strategy, Advances in Epidemiological Modeling and Control of Viruses, (2023), 323–363.
Shakeri Barzoki, M. A. , Kazemi, R. and Asheghi, R. (2026). Optimal control of an SIR model with effect of vaccination and a generic nonlinear infection rate. Computational Methods for Differential Equations, 14(4), 1749-1762. doi: 10.22034/cmde.2025.65636.3030
MLA
Shakeri Barzoki, M. A. , , Kazemi, R. , and Asheghi, R. . "Optimal control of an SIR model with effect of vaccination and a generic nonlinear infection rate", Computational Methods for Differential Equations, 14, 4, 2026, 1749-1762. doi: 10.22034/cmde.2025.65636.3030
HARVARD
Shakeri Barzoki, M. A., Kazemi, R., Asheghi, R. (2026). 'Optimal control of an SIR model with effect of vaccination and a generic nonlinear infection rate', Computational Methods for Differential Equations, 14(4), pp. 1749-1762. doi: 10.22034/cmde.2025.65636.3030
CHICAGO
M. A. Shakeri Barzoki , R. Kazemi and R. Asheghi, "Optimal control of an SIR model with effect of vaccination and a generic nonlinear infection rate," Computational Methods for Differential Equations, 14 4 (2026): 1749-1762, doi: 10.22034/cmde.2025.65636.3030
VANCOUVER
Shakeri Barzoki, M. A., Kazemi, R., Asheghi, R. Optimal control of an SIR model with effect of vaccination and a generic nonlinear infection rate. Computational Methods for Differential Equations, 2026; 14(4): 1749-1762. doi: 10.22034/cmde.2025.65636.3030