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


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


EIv2 = - (X(6,1)*(sinh(la*x2)+sin(la*x2))/(2*la) - …

X(7,1)*(cosh(la*x2)-cos(la*x2))/(2*la^2) - …

X(8,1)*(sinh(la*x2)-sin(la*x2))/(2*la^3));

EIv3 = - (X(10,1)*(sinh(la*x3)+sin(la*x3))/(2*la) - …

X(11,1)*(cosh(la*x3)-cos(la*x3))/(2*la^2) - …

X(12,1)*(sinh(la*x3)-sin(la*x3))/(2*la^3));

EIv4 = - (X(14,1)*(sinh(la*x4)+sin(la*x4))/(2*la) - …

X(15,1)*(cosh(la*x4)-cos(la*x4))/(2*la^2) - …

X(16,1)*(sinh(la*x4)-sin(la*x4))/(2*la^3));

subplot

(2,2,1), plot (x1, EIv1); axis ([0 4 - 0.5 0.5]); grid on

subplot

(2,2,2), plot (x2, EIv2); axis ([0 6 - 0.6 0.6]); grid on

subplot

(2,2,3), plot (x3, EIv3); axis ([0 3 - 0.5 0.5]); grid on

subplot

(2,2,4), plot (x4, EIv4); axis ([0 1 - 1    1   ]); grid on

Перед выполнением данного протокола из отдельного М-файла в окно команд помещается вектор Х -

вектор значений граничных параметров балки при wi . Выполняется это прогонкой программы, решающей систему уравнений (3.16). Далее отрабатывается протокол построения форм колебаний при собственной частоте wi . Первые 5 форм собственных колебаний балки при условии, что

, представлены на рис 3.6. Можно рекомендовать также объединить протоколы вычисления вектора Х и построения форм колебаний, но только после четкого усвоения навыка построения графиков и выбора их масштабов.

Рис. 3.6

 

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

 

3.2.5. Расчет на вынужденные колебания

Рис. 3.7

 

 

 

 

 

(3.18)

Рис. 3.8

 

 

 

 

(3.19)

Рис. 3.9

 

 

 

 

(3.20)

Рис. 3.10

 

Расчет балки на вынужденные колебания сводится к решению матричного уравнения А*Х*

= -

В, где матрица А повторяет матрицу частотного уравнения (см. п.3.2.3), а матрица В представлена в этом разделе. Программа также записывается в отдельном М-файле, куда переносится матрица А* задачи определения частот собственных колебаний и дополняется элементами вектора В. Сама программа имеет вид.




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