This should be simple fix but I am stuck. I am attempting to take the real part of an extensive complex number involving a variable. I was able to make it work with my first expression (delta1) but the same method fails with my second expression (delta2). I don't understand why it would work for one and not the other. Any help would be much appreciated!
Here's my code:
k1 = 2;
k2 = 2;
k12 = 2;
m1 = 1;
m2 = 1;
f = 1;
b = .1;
M = {{m1, 0}, {0, m2}};
K1 = {{k1 + k12, -k12}, {-k12, k2 + k12}};
a = {a1, a2};
F0 = {f, f};
B = {b, b};
Solve[Det[K1 - \[Omega]^2*M] == 0];
Sol = LinearSolve[(K1 + I*B*\[Omega] - \[Omega]^2*M), F0];
C1 = Sol[[1]];
C2 = Sol[[2]];
\[Delta]1 =
Simplify[Im[ComplexExpand[C1]], Element[\[Omega], Reals]]/
Simplify[Re[ComplexExpand[C1]], Element[\[Omega], Reals]]
\[Delta]2 =
Simplify[Im[ComplexExpand[C1]], Element[\[Omega], Reals]]/
Simplify[Re[ComplexExpand[C1]], Element[\[Omega], Reals]]
Attachments: