This work was supported by Scientific Research Fund of Hunan Provincial Education Department Nos.18A325, Research and Innovation team of Hunan Institute of Science and Technology (Grant No. 2019-TD-15), NNSF of China Grant Nos. 12071413, NSF of Guangxi (2018GXNSFDA138002), the European Union's Horizon 2020 Research and Innovation Programme under the Marie Sklodowska-Curie (823731CONMECH). Also, this work was partially supported by Hunan Provincial Key Laboratory of Mathematical Modeling and Analysis in Engineering of Changsha University of Science and Technology Grant No. 2018MMAEZD05.