Для построения эпюр программируем выражения

№ п/п	Граничные параметры
1	= Х(1,1) = - 22.8591 кНм
2	= Х(2,1) = 1.0959 кНм2
3	= Х(3,1) = - 59.0041 кН
4	= Х(4,1) = 9.0752 кН
5	= Х(5,1) = 10.7927 кН
6	= Х(6,1) = - 7.2730 кНм
7	= Х(7,1) = - 0.9959 кН
8	= Х(8,1) = - 9.6477 кНм
9	= Х(9,1) = 20.7927 кН
10	= Х(10,1) = - 0.9959 кН
11	= Х(11,1) = - 21.7175 кН
12	= Х(12,1) = - 0.8570 кНм2
13	= Х(13,1) = 6.1450 кНм
14	= Х(14,1) = 0.9959 кН
15	= Х(15,1) = 10.7927 кН
16	= Х(16,1) = 2.8455 кНм3
17	= Х(17,1) = - 1.7751 кНм2
18	= Х(18,1) = 13.7229 кНм
19	= Х(19,1) = - 0.9959 кН
20	= Х(20,1) = - 21.7175 кН

3.3.2. Построение эпюр напряженно-деформированного состояния
          элементов рамы
Для построения эпюр программируем выражения (3.12), (3.13)
стержень 0-1
х = 0 : 0.001 : 1.0; q = 10.0;
EIv = - (X(2,1)*x - X(4,1)*x.^ 3 / 6 + q*x.^4 / 24);
EIfi
= - (X(2,1) - X(4,1)*x.^ 2 / 2 + q*x.^3 / 6);
Q = X(4,1) - q*x; M = X(4,1)*x - q*x.^ 2 / 2;
subplot (2,2,1), plot (x, EIv); axis
([0 1 - 0.4   0.4]); grid on
subplot (2,2,2), plot (x, EIfi); axis ([0 1 - 2        2]); grid on
subplot (2,2,3), plot (x, Q); axis
([0 1 - 10    10]); grid on
subplot (2,2,4), plot (x, M); axis
([0 1 - 5      5 ]); grid on
стержень 1-2
х = 0 : 0.001 : 1.0; q = 10.0;
EIv = - (X(17,1)*x - X(8,1)*x.^ 2 / 2- X(9,1)*x.^ 3 / 6 + q*x.^4 / 24);
EIfi
= - (X(17,1) - X(8,1)*x- X(9,1)*x.^ 2 / 2 + q*x.^3 / 6);
Q = X(9,1) - q*x; M = X(8,1) + X(9,1)*x - q*x.^ 2 / 2;
subplot (2,2,1), plot (x, EIv); axis
([0 1 - 0.4   0.4]); grid on
subplot (2,2,2), plot (x, EIfi); axis ([0 1 - 2        2]); grid on
subplot (2,2,3), plot (x, Q); axis
([0 1 - 22    22]); grid on
subplot (2,2,4), plot (x, M); axis
([0 1 - 10    10]); grid on
стержень 1-3
а). интервал (0;

) м.
х = 0 : 0.001 : 0.75;

EIv = - (X(16,1) + X(17,1)*x- X(18,1)*x.^
2 / 2 - X(19,1)*x.^
3 / 6);
EIfi
= - (X(17,1) - X(18,1)*x- X(19,1)*x.^ 2 / 2);
Q = X(19,1); M = X(18,1) + X(19,1)*x;
subplot (2,2,1), plot (x, EIv); axis
([0 0.25 - 3   3]); grid on
subplot (2,2,2), plot (x, EIfi); axis ([0 0.25 - 6    6]); grid on
subplot (2,2,3), plot (x, Q); axis
([0 0.25 - 2        2]); grid on
subplot (2,2,4), plot (x, M); axis
([0 0.25 - 14   14]); grid on
б). интервал (

;

) м.
х = 0.25 : 0.001 : 1.0; m = 20.0;
EIv=-(X(16,1)+X(17,1)*x-X(18,1)*x.^ 2/2-X(19,1)*x.^ 3/6+m*(x-0.25).^ 2/2);
EIfi
= - (X(17,1) - X(18,1)*x- X(19,1)*x.^ 2 / 2+m*(x-0.25));
Q = X(19,1); M = X(18,1) + X(19,1)*x - m;
subplot (2,2,1), plot (x, EIv); axis
([0.25   1.0 - 3   3]); grid on
subplot (2,2,2), plot (x, EIfi); axis ([0.25   1.0 - 6    6]); grid on
subplot (2,2,3), plot (x, Q); axis
([0.25   1.0 - 2        2]); grid on
subplot (2,2,4), plot (x, M); axis
([0.25   1.0 - 15   15]); grid on
стержень 2-4
а). интервал (0;

) м.
х = 0 : 0.001 : 0.5;
EIv = - (X(16,1) + X(12,1)*x- X(13,1)*x.^
2 / 2 - X(14,1)*x.^
3 / 6);
EIfi
= - (X(12,1) - X(13,1)*x- X(14,1)*x.^ 2 / 2);
Q = X(14,1); M = X(13,1) + X(14,1)*x;
subplot (2,2,1), plot (x, EIv); axis
([0 0.5 - 3   3]); grid on
subplot (2,2,2), plot (x, EIfi); axis ([0 0.5 - 6    6]); grid on
subplot (2,2,3), plot (x, Q); axis
([0 0.5 - 2        2]); grid on
subplot (2,2,4), plot (x, M); axis
([0 0.5 - 15   15]); grid on
б). интервал (

;

) м.
х = 0.5 : 0.001 : 1.0; f = 60.0;
EIv=-(X(16,1)+X(12,1)*x-X(13,1)*x.^ 2/2-X(14,1)*x.^ 3/6+ f *(x-0.5).^ 3/6);
EIfi
= - (X(12,1) - X(13,1)*x- X(14,1)*x.^ 2 / 2+ f *(x-0.5).^ 2/2);
Q = X(14,1) - f; M = X(13,1) + X(14,1)*x - f *(x-0.5);
subplot (2,2,1), plot (x, EIv); axis
([0.5   1.0 - 3   3]); grid on
subplot (2,2,2), plot (x, EIfi); axis ([0.5   1.0 - 6    6]); grid on
subplot (2,2,3), plot (x, Q); axis

([0.5   1.0 - 60    60]); grid on
subplot (2,2,4), plot (x, M); axis
([0.5   1.0 - 25   25]); grid on
Эпюры прогибов EIJ(x), углов поворота EIj(x), поперечных сил Q(x) и изгибающих моментов М(x) представлены на рис. 3.15

а)	b)

EIJ(x) кНм3	EIj(x) кНм2

с)	d)

Q(x) к	М(x) кНм
Рис. 3.15

Содержание раздела

Главная сайта