x、y代表两个广义坐标，ab是有效载荷处加速度，需要用x[t]、y[t]、x'[t]、y'[t]、x''[t]、y''[t]表示出来，但为什么ab的计算结果是二维列表呢？不应该是一维吗？我想得到的是ab随时间变化曲线。
Clear["Global`*"];
equ = {x'[t] == y[t], 2 y'[t] == -0.1 y[t] - Sin[x[t]], x[0] == 0, y[0] == 1.9};
d = {x[t], y[t]} /. First@NDSolve[equ, {x, y}, {t, 0, 100}]; (*------位移------*)
v = D[d, t]; (*------速度------*)
a = D[d, {t, 2}]; (*------加速度------*)
ab =Cos[d[[2]]]*(a[[1]] - d[[1]]*v[[2]]^2) -Sin[d[[2]]]*(2*v[[1]]*v[[2]] + d[[1]]*a[[2]]) ;
Dimensions[d]
Dimensions[v]
Dimensions[Subscript[a, b]]
Plot[ab, {t, 0, 100}]

这里是ab[t]的数学表达式

不管怎么样，你先把下标改掉，我再来看你的代码问题在哪。
精品区常见问题栏有一篇专门讲下标的。总之就是新手勿用下标。

(*x、y代表两个广义坐标，ab是有效载荷处加速度，需要用x[t]、y[t]、x'[t]、y'[t]、x''[t]、y''[t]表示出来，但为什么ab的计算结果是二维列表呢？不应该是一维吗？我想得到的是ab随时间变化曲线。*)
Clear["Global`*"];
equ = {x'[t] == y[t], 2 y'[t] == -0.1y[t] - Sin[x[t]], x[0] == 0,y[0] == 1.9};
d = {x[t], y[t]} /.First@NDSolve[equ,{x, y}, {t, 0, 100}];
(*------位移------*)
v = D[d, t];
(*------速度------*)
a = D[d, {t,2}];
(*------加速度------*)
ab = Cos[d[[2]]]*(a[[1]] - d[[1]]*v[[2]]^2)-Sin[d[[2]]]*(2*v[[1]]*v[[2]] + d[[1]]*a[[2]]);
Dimensions[d]
Dimensions[v]
Dimensions[ab]
其中：
很奇怪，为什么会有Dimensions[ab]=2？输出值不该是1吗？另外，Dimensions[ab]=2的情况下竟然画出了图，还是1条曲线。

此问题已在别处得到解决。

……这个问题可以简化如下：
为什么
aaa + bbb // Dimensions
输出
{2}
？
这个在帮助里其实写了：