/span> , то матриця коефіцієнтів перед , і буде матрицею піддатливості системи.
Додаток
Розрахунок у математичному пакеті Maple 9.0
> restart;
> Heaviside (0): = 0:
> # Рівняння рівноваги для визначення реакцій зв'язків:
> sys_stat: =
> Y [A] + Y [O] + a * r [1] ^ 2/2 * cos (theta [1]-pi/2),
> Z [A] + Z [O] + a * r [1] ^ 2/2 * sin (theta [1]-pi/2),
> M [upr] + r [1] * R [Oy [1]] + a * r [1] ^ 3/3,
> P [y]-Y [A] + R [B] * cos (theta [3]),
> P [z]-Z [A] + R [B] * sin (theta [3]),
> -2 * r [2] * P [z] * cos (pi-theta [2]) -2 * r [2] * P [y] * sin (pi-theta [2]) + [2] * R [B] * cos (theta [2] -pi/2-theta [3]),
> R [B] = R [C],
> N [C]-R [C] * sin (theta [3]),
> F [upr]-R [B] * cos (theta [3]);
В В
> # Розподіл згинальних моментів для стрижнів 1,2:
> # Перепроектіруем сили на локальні осі координат:
> R [Oy [1]]: = Y [O] * cos (theta [1]-pi/2) + Z [O] * sin (theta [1]-pi/ 2);
> R [Ay [1]]: = Y [A] * cos (theta [1]-pi/2) + Z [A] * sin (theta [1]-pi/ 2);
> R [Ay [2]]: =-Y [A] * cos (theta [2]-pi/2)-Z [A] * sin (theta [2]-pi /2);
> P [y [2]]: = P [y] * cos (theta [2]-pi/2) + P [z] * sin (theta [2]-pi/ 2);
В В В В
> # Рівняння згинальних моментів при розгляді рівноваги лівої частини для кожного з изгибающихся стрижнів:
> M [x [1]]: = M [upr] + z [1] * R [Oy [1]] + a * z [1] ^ 3/6; p>
> M [x [3]]: = R [B] * z [2] * cos (theta [2] -pi/2-theta [3]) + (z [2 ]-r [2]) * R [Ay [2]] * Heaviside (z [2]-r [2]);
В В
> # Вихідні дані:
> pi: = Pi:
k: = 1.36 * 10 ** 5;
> # Рішення рівнянь статики:
> solve ({sys_stat}, {Y [A], Z [A], Y [O], Z [O], N [C], M [upr], F [upr ], R [B], R [C]});
> assign (%);
В В
> # Побудуємо епюри згинальних моментів:
> P [y]: = 118;
P [z]: = 236;: = 472;
В В В
> plot (M [x [1]], z [1] = 0 .. r [1], color = blue, thickness = 3, title = "Епюра згинальних моментів n для першого стрижня ");
> plot (M [x [3]], z [2] = 0 .. 3 * r [2], color = blue, thickness = 3, title = "Епюра згинальних моментів n для другого стрижня ");
В В
> # Розрахунок на міцність ст...