Решение задачи продольных колебаний бруска, описанных волновым уравнением с помощью MATLAB
17

Согласно разработанному алгоритму, на первом этапе мы задаеми описываем все необходимые для расчетов переменные.

m - число шагов по времени

n - число шагов по координате

E - модуль упругости материала

ro - плотность материала

L - длина бруска

T - %время моделирования

dx=L/n - шаг по расстоянию

dt=T/m - шаг по времени

cur=dt/dx*(E/ro)^0.5 - число Курранта

Далее задаем начальные и граничные условия для расчета модели.Задаем циклический процесс для расчета. Для организации вычислений с заданным числом повторяющихся циклов используем конструкцию имеющую следующий вид:

forvar=Выражение. Инструкция..... Инструкция end

Результатом решения является матрица значений продольного перемещения на каждом шаге по времени и по расстоянию (рис. 2.2).

Рисунок 2.2. – фрагмент матрицы решений