“…Based on numerical integration method and analytic approximate method, such as Chu Fulei et al, Li Zhenping et al [29], Li Hongkun et al [30], Yang Yongfeng et al [31], Liu Changli et al [32],Chen Guo [33] MA H, ZHAO X, TENG Y, et al [34], Xiong et al [35] used Runge-Kutta method, Newmark-β method, shooting method, incremental method and incremental harmonic balance method [6,[29][30][31][32][33][34][35]. Numerical solution method generally needs the condition to be simplified, so objectively produces the disparity with the actual state; The numerical solution method needs to assign a value to the model, for this reason that this general model becomes the specific sampling, the solution obtained is equivalent to the experiment result.…”