已知函数 f(x)=13x3+ax2+xf\left( x\right) = \frac{1}{3}{x}^{3} + a{x}^{2} + xf(x)=31x3+ax2+x 在 (0,+∞)\left( {0, + \infty }\right)(0,+∞) 上单调递增，则实数 aaa 的取值范围为( )A. (−∞,−1]( - \infty , - 1\rbrack(−∞,−1] B. [−1,1]\left\lbrack {-1,1}\right\rbrack[−1,1] C. [1,+∞)\lbrack 1, + \infty )[1,+∞) D. [−1,+∞)\lbrack - 1, + \infty )[−1,+∞)