M [x [1], 2]: = subs (a = 0, P [y] = 0, P [z] = 1, M [x [1]]);
> M [x [3], 2]: = subs (a = 0, P [y] = 0, P [z] = 1, M [x [3]]);
> M [upr, 2]: = subs (a = 0, P [y] = 0, P [z] = 1, M [upr]);
> F [upr, 2]: = subs (a = 0, P [y] = 0, P [z] = 1, F [upr]);
> R [B, 2]: = subs (a = 0, P [y] = 0, P [z] = 1, R [B]);
В
> # Обчислимо коефіцієнти матриці податливості:
> G: = matrix (2,2, [delta [11], delta [12], delta [21], delta [22]]);
В
>
>
>
> delta [21]: = 'delta [12]';
> delta [11]: = simplify (evalf (delta [11]));
> delta [12]: = simplify (evalf (delta [12]));
> delta [21]: = simplify (evalf (delta [21]));
> delta [22]: = simplify (evalf (delta [22]));
В В В В
> 'G' = convert (G, Matrix);
В
> # Вирази для ліниейних переміщень точки М:
> 'D_M [y]' = D_M [y];
> 'D_M [z]' = D_M [z];
В В
> # Поклавши в них a = 0, отримаємо:
> 'D_M [y]' = subs (a = 0, D_M [y]);
> 'D_M [z]' = subs (a = 0, D_M [z]);
В В
> # Матриця коефіцієнтів перед P [y] і P [z]:
> G: = lhs (simplex [display] ([D_M [y] = const, D_M [z] = const], [P [y], P [z]]));
В