MATLAB в инженерных и научных расчетах


Продолжение таблицы 2.3 - часть 19


EIfi

= - (X(10,1) – X(11,1).*x – X(12,1).*x.^2/2 – f1.*(x –2).^2/2);

Q = X(12,1) + f1; M = X(11,1) + X(12,1).*x + f1.*(x –2);

subplot (2, 2, 1), plot (x, EIn); axis ( [2  3  -150  150] ); grid on

subplot (2, 2, 2), plot (x, EIfi); axis ( [2  3  -150  150] ); grid on

subplot (2, 2, 3), plot (x, Q); axis ( [2  3  -100  100] ); grid on

subplot (2, 2, 4), plot (x, M); axis ( [2  3  -100  100] ); grid on

Стержень 3 – 4

x = 0 : 0.001 : 1;

EIn

= - (X(14,1).*x – X(15,1).*x.^2/2 – X(16,1).*x.^3/6);

EIfi

= - (X(14,1) – X(15,1).*x– X(16,1).*x.^2/2);

Q = X(16,1); M = X(15,1) + X(16,1).*x;

subplot (2, 2, 1), plot (x, EIn); axis ( [0  1  -150  150] ); grid on

subplot (2, 2, 2), plot (x, EIfi); axis ( [0  1  -150  150] ); grid on

subplot (2, 2, 3), plot (x, Q); axis ( [0  1  -100  100] ); grid on

subplot (2, 2, 4), plot (x, M); axis ( [0  1  -100  100] ); grid on.

Следует помнить, что в окно команд необходимо вводить значения неизвестных начальных параметров балки, т. е. вектор Х. Это легко выполняется с помощью отдельного М-файла, на котором записана процедура решения матричного уравнения неразрезной балки. После вектора Х нужно ввести также значения внешней нагрузки и их координаты. Эпюры EIn(х), EIj(х), Q(х) и М(х) представлены на рис. 3.1.

3.2.3. Определение частот собственных колебаний

Рис. 3.4

Исходные данные:

Матрица фундаментальных функций поперечных колебаний

 

 

 

 

 

 

 

 

 

;

 

 

 

 

(3.14)

 

 

Аi =

A11

A12

-A13

-A14

l4A14

A11

-A12

-A13

-l4A13

-l4A14

A11

A12

-l4A12

-l4A13

l4A14

A11

 

 

 

 

 

 

Частоты определяются как корни трансцендентного уравнения

çА*(w) ç= 0,

где матрица А*(w) берется из задачи статики с заменой матриц фундаментальных функций изгиба на матрицы функций поперечных колебаний

 

 

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

 

<


Начало  Назад  Вперед