se\[Gamma] = \[Gamma] -> 0.9;
seg = g -> 9.8;
Subscript[se\[Rho], 1] = Subscript[\[Rho], 1] -> 1024;
sev = v -> 0.3;
Subscript[se\[Rho], e] = Subscript[\[Rho], e] -> 0.9;
seE = E -> 6000;
seH = H -> 100;
fValues = Range[0.01, 100.01, (100.01 - 0.01)/100];
s22 = ((g*\[Omega]^2)/
H - (g k/
H ((1 + (d^3 E (k)^4)/(12 (1 - v^2))) Tanh[
0.2` k] + (1 + (d^3 E (k)^4 \[Gamma])/(
12 (1 - v^2) )) Tanh[0.8` k] +
k (1 - \[Gamma]) Subscript[\[Rho], e]
Tanh[0.2` k] Tanh[
0.8` k]) (1 + \[Sqrt](1 - (4 (1 - \[Gamma]) (1 + (
E (k)^4)/(12 (1 - v^2))) Tanh[0.2` k] Tanh[
0.8` k] (1 + \[Gamma] Tanh[0.2` k] Tanh[
0.8` k] +
k*Subscript[\[Rho],
e] (Tanh[
0.2` k] + \[Gamma] Tanh[0.8` k])))/((1 + (
d^3 E (k)^4)/(12 (1 - v^2) )) Tanh[
0.2` k] + (1 + (d^3 E (k)^4 \[Gamma])/(
12 (1 - v^2))) Tanh[0.8` k] +
k (1 - \[Gamma]) Subscript[\[Rho], e]
Tanh[0.2` k] Tanh[
0.8` k])^2)))/(2 (1 + \[Gamma] Tanh[0.2` k] Tanh[
0.8` k] +
k* Subscript[\[Rho],
e] (Tanh[0.2` k] + \[Gamma] Tanh[0.8` k])))) /.
se\[Gamma] /. seg /. Subscript[se\[Rho], 1] /. sev /.
Subscript[se\[Rho], e] /. seE /. seH /. d -> 0.05;
(*Define a function to calculate k for a given fValue*)
findK[fValue_] := Module[{z22, y22}, z22 = s22 /. \[Omega] -> fValue;
y22 = FindInstance[z22 == 0, k, Reals];
If[y22 === {}, None, k /. y22[[1]]]]
(*Use Map to apply the function to each fValue*)
solutions22 = Map[findK, fValues];
(*Clear temporary variables*)
Clear[findK];
solutions22

……你先想想为什么
\[Gamma] E /. se\[Gamma] /. seE
输出的结果不对。
还不明白就把精品区常见问题合集里的内容先看看。