Anexo S: Programa dos torques de reação com atrito 192
(* Ponto na barra de cima *)
{a7, b7} = {L1/2 + q8*Cos[t7], q8*Sin[t7]};
(* calculo dos pontos nas extremidades da barra de cima *)
{a5, b5} = LJ*{Sin[t7], -Cos[t7]} + {a7, b7};
{a3, b3} = LJ*{-Sin[t7], Cos[t7]} + {a7, b7};
(* calculo dos pontos de conexao das barras verticais *)
sol6 = Solve[{(x - L1)^2 + y^2 == LA^2, (x - a5)^2 +
(y - b5)^2 == LB^2}, {x, y}];
sol2 = Solve[{x^2 + y^2 == LD^2, (x - a3)^2 +
(y - b3)^2 == LE^2}, {x, y}];
{a61, b61} = ReplaceAll[{x, y}, sol6[[1]]];
{a62, b62} = ReplaceAll[{x, y}, sol6[[2]]];
{a21, b21} = ReplaceAll[{x, y}, sol2[[1]]];
{a22, b22} = ReplaceAll[{x, y}, sol2[[2]]];
(* achar os pontos q. dao certo *)
fsol = {};
If[((a61 - a62)*b5 - (b61 - b62)*(a5 - L1)) >= 0 &&
((a21 - a22)*b3 - (b21 - b22)*a3) >= 0,
fsol = Append[fsol, {{a61, b61}, {a22, b22}}]];
If[((a62 - a61)*b5 - (b62 - b61)*(a5 - L1)) >= 0 &&
((a21 - a22)*b3 - (b21 - b22)*a3) >= 0,
fsol = Append[fsol, {{a62, b62}, {a22, b22}}]];
If[((a61 - a62)*b5 - (b61 - b62)*(a5 - L1)) >= 0 &&
((a22 - a21)*b3 - (b22 - b21)*a3) >= 0,
fsol = Append[fsol, {{a61, b61}, {a21, b21}}]];
If[((a62 - a61)*b5 - (b62 - b61)*(a5 - L1)) >= 0 &&
((a22 - a21)*b3 - (b22 - b21)*a3) >= 0,
fsol = Append[fsol, {{a62, b62}, {a21, b21}}]];
barra = Graphics[{RGBColor[1, 0, 0],
Line[{{L1/2, 0}, {a7, b7}}]}];
For[i = 1, i <= Length[fsol], i++,
(* Limpar variaveis *)
Clear[a6, b6, a2, b2, vertices, ssol, OE, EC, F34, F54,
F17, F12, F16, M12, M16, t3, kt3, kt5, t1, hex];
{a6, b6} = fsol[[i]][[1]];
{a2, b2} = fsol[[i]][[2]];
vertices = {{0, 0}, {a2, b2}, {a3, b3}, {a5, b5},