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


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


/p>

 

 

 

Результаты вычисления граничных параметров представлены в таблице 3.1.

Таблица 3.1

1

9

2

10

3

11

4

12

5

13

6

14

7

15

8

16

 

Как видно из этого примера, не требуется переставлять строки матриц А и В

для исключения нулевых ведущих элементов, что является дополнительным преимуществом MATLAB по сравнению с непосредственным программированием метода исключения Гаусса в таких средах, как Delphi, Visual Fortran, Visual C++ и т.п.

 

3.2.2. Построение эпюр прогибов, углов поворота, поперечных сил

          и изгибающих моментов

Эпюры параметров изгиба можно строить для каждого стержня в отдельности, используя соотношения метода начальных параметров. При этом удобно выполнять построение эпюр в одном окне, для чего можно привлекать процедуру subplot. Для соблюдения требуемого масштаба можно использовать процедуру axis. Масштабная сетка эпюр появится, если задействовать процедуру grid on. Выполнение этих требований приводит к следующим операторам построения эпюр

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

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

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

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

Протокол построения эпюр удобно поместить в отдельный М-файл и по мере необходимости корректировать масштаб изображения графиков. Выражения для параметров изгиба стержня по методу начальных параметров имеют вид

              (3.12)

где элементы от внешней нагрузки запишутся следующим образом:

 (3.13)

Положительные направления внешней нагрузки и их координаты показаны на рис. 3.3

 

Рис.3.3

Для балки по рис. 3.1 протоколы построения эпюр примут вид:

Стержень 0 – 1

x = 0 : 0.001 : 2.0;

EIn

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

EIfi

= - (X(2,1) – X(4,1).*x.^2/2);

Q = X(4,1); M = X(4,1).*x;




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