To define a function of a parameter, you need to add an underscore to the parameter, to make it a pattern:
\[Omega] = 8.0001*Pi;
\[CapitalOmega] = 2*Pi/To;
Subscript[r, n_] = n*\[CapitalOmega]/\[Omega];
\[Xi] = 0.05;
\[Rho]st = .1;
k = Fo/\[Rho]st;
Subscript[DAF, n_] =
1/Sqrt[((1 - Subscript[r, n]^2)^2 + (2 \[Xi]*Subscript[r, n])^2)];
Subscript[\[CapitalPsi], n_] =
ArcTan[(2 \[Xi]*Subscript[r, n])/(1 - Subscript[r, n]^2)];
The plot does not work because you have not given a numeric value to the symbol To. For example, the following computes fine:
To = 1;
\[Omega] = 8.0001*Pi;
\[CapitalOmega] = 2*Pi/To;
Subscript[r, n_] = n*\[CapitalOmega]/\[Omega];
\[Xi] = 0.05;
\[Rho]st = .1;
k = Fo/\[Rho]st;
Subscript[DAF, n_] =
1/Sqrt[((1 - Subscript[r, n]^2)^2 + (2 \[Xi]*Subscript[r, n])^2)];
Subscript[\[CapitalPsi], n_] =
ArcTan[(2 \[Xi]*Subscript[r, n])/(1 - Subscript[r, n]^2)];
Subscript[F,
ss1] = (Fo/(3 k)) + (Subscript[DAF, n]*(4.5 Fo/(k*(n*Pi)^2))*
Cos[(n*\[CapitalOmega]*t) - Subscript[\[CapitalPsi], n]]) /.
n -> 1;
Subscript[F,
ss5] = (Fo/(3 k)) +
Sum[(Subscript[DAF, n]*(4.5 Fo/(k*(n*Pi)^2))*
Cos[(n*\[CapitalOmega]*t) - Subscript[\[CapitalPsi], n]]), {n,
1, 5, 1}];
Subscript[F,
ss10] = (Fo/(3 k)) +
Sum[(Subscript[DAF, n]*(4.5 Fo/(k*(n*Pi)^2))*
Cos[(n*\[CapitalOmega]*t) - Subscript[\[CapitalPsi], n]]), {n,
1, 10, 1}];
Subscript[F,
sst] = {Subscript[F, ss1], Subscript[F, ss5], Subscript[F, ss10]};
Plot[Evaluate[Subscript[F, sst]], {t, 0, 2},
AxesLabel -> {"Time(s)", "Displacement(in)"}]