lign="justify">> D_theta [1], D_theta [2], D_theta [3], Theta [M],
> dBy (3), dBz (3),
> dAy (1), dAz (1),
> dMy (2), dMz (2)];
В
> # Висловимо рішення через a, P [y] і P [z].
> sys_graph: =
> F [upr] =-k * D_C [y],
> M [upr] =-c * D_theta [1],
> D_M [y] = dAy (1) + dMy (2)-D_theta [1] * r [1] * sin (theta [1]) -2 * D_theta [2] * r [2] * sin (theta [2]),
> D_M [z] = dAz (1) + dMz (2) + D_theta [1] * r [1] * cos (theta [1]) +2 * D_theta [2] * r [2] * cos (theta [2]),
> D_M [y] = D_C [y] + dBy (3) + dMy (2)-D_theta [3] * r [3] * sin (theta [3] + pi) - 3 * D_theta [2] * r [2] * sin (theta [2]),
> D_M [z] = dBz (3) + dMz (2) + D_theta [3] * r [3] * cos (theta [3] + pi) +3 * D_theta [2 ] * r [2] * cos (theta [2]),
>
> # Перепроектіруем локальні переміщення:
> dBy (3) =-R [B] * r [3]/E/F [3] * cos (theta [3]),
> dBz (3) =-R [B] * r [3]/E/F [3] * sin (theta [3]),
> dAy (1) = eval (V [1], z [1] = r [1]) * cos (theta [1]-pi/2),
> dAz (1) = eval (V [1], z [1] = r [1]) * sin (theta [1]-pi/2),
> dMy (2) = eval (V [2], z [2] = 3 * r [2]) * cos (theta [2]-pi/2),
> dMz (2) = eval (V [2], z [2] = 3 * r [2]) * sin (theta [2]-pi/2),
> Theta [M] = D_theta [2]-eval (diff (V [2], z [2]), z [2] = 3 * r [2]): p>
> SLV: = solve ({sys_graph},
> convert (params, set));
В
> # Обчислимо деформації елементів конструкції при вихідних навантаженнях (a = 472, P [y] = 118, P [z] = 236):
> subs (a = 472, P [y] = 118, P [z] = 236, SLV);
В В
> # Обчислимо деформації елементів конструкції при навантаженнях (a = 0, P [y] = 118, P [z] = 236):
> subs (a = 0, P [y] = 118, P [z] = 236, SLV);
В В
> assign (SLV);
> # Розрахунок матриці податливості:
> # Обчислимо матрицю податливості через інтеграл Мора:
> # Перше нагружение (a = 0, P [y] = 1, P [z] = 0):
> M [x [1], 1]: = subs (a = 0, P [y] = 1, P [z] = 0, M [x [1]]);
> M [x [3], 1]: = subs (a = 0, P [y] = 1, P [z] = 0, M [x [3]]);
> M [upr, 1]: = subs (a = 0, P [y] = 1, P [z] = 0, M [upr]);
> F [upr, 1]: = subs (a = 0, P [y] = 1, P [z] = 0, F [upr]);
> R [B, 1]: = subs (a = 0, P [y] = 1, P [z] = 0, R [B]);
В
> # Друге нагружение (a = 0, P [y] = 0, P [z] = 1):
>...