This work is supported by the National Natural Science Foundation of China, (Grant Nos. 11771059, 11301551), the Natural Science Foundation of Hunan Province, China (Grant No. 2017JJ3525) and Hunan Provincial Key Laboratory of Mathematical Modeling and Analysis in Engineering, China (Grant No. 2017TP1017).