Метод конечных разностей решение в excel

Метод конечных разностей в расчете прогибов пластин

Опубликовано 10 Апр 2016
Рубрика: Механика | 5 комментариев

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

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

Мы не будем писать дифференциальные уравнения, а воспользуемся готовой схемой замещения и составим систему уравнений в конечных разностях при помощи простого оператора формулы, а затем легко решим эту систему, используя мощь MS Excel.

Приступим к рассмотрению практической задачи и научимся применять в реальных расчетах метод конечных разностей.

Расчет в Excel  прогибов пластины в узлах сетки.

Квадратная стальная пластина нагружена равномерно распределенной нагрузкой по всей площади. Границы пластины шарнирно оперты на основание — пластина свободно лежит, опираясь контуром.

Необходимо рассчитать прогибы в разных точках пластины.

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

Наносим на чертеж пластины сетку с шагом a/4. В точках пересечения линий сетки находятся узлы, в которых и будем определять величину прогибов w_i (i— номер узла).

Присваиваем номера узлам, начиная с «0».

Номер «0» назначаем всем узлам, расположенным на контуре. Эти узлы по условиям задачи лежат на опорной поверхности и после приложения нагрузки остаются на прежнем месте – прогиб этих узлов равен нулю w₀=0.

Из условий приложения и характера нагрузки, а также квадратной формы пластины (условия центральной симметрии) очевидно, что прогиб будет одинаковым во всех узлах «1». Поэтому этим узлам и присвоен одинаковый номер.

Руководствуясь теми же признаками симметрии, назначаем характерным внутренним узлам номер «2».

Понятно, что наибольший прогиб будет в центре пластины, назначаем этому узлу номер «3». Он один, других узлов на пластине, имеющих такой же прогиб нет.

На рисунке вы видите законтурные (зеркальные или мнимые) узлы. Они участвуют в расчете, поэтому нам нужны. Номера им присвоим исходя из простых правил:

1. Если линия контура опирается в этом месте шарнирно, то законтурному (зеркальному) узлу присваивается номер «-N» Где «N» — номер симметричного относительно линии контура внутреннего узла (w_{— N}=-w_N).

2. Если контур жестко закреплен (жесткое защемление), то соответствующему зеркальному узлу присваивается номер «N» — такой же, как и симметричному внутреннему узлу (w_N=+w_N).

Числовые исходные данные о материале и размере пластины, а также о нагрузке записываем в таблицу MS Excel в ячейки D3…D7.

Расчет:

1. В рассматриваемом варианте расчета пластина делится сеткой на 4 квадрата по горизонтали и вертикали. Определяем шаг сетки в мм

в ячейке D9: δ=a/4

2. Цилиндрическую жесткость пластины в Н*мм вычисляем по общеизвестной формуле

в D10: D=(E*h³)/(12*(1- μ²)

3. Для каждого узла, в котором мы хотим найти прогиб, составляем уравнение, применив метод конечных разностей.

Для составления уравнений используем оператор формулы, изображенный выше на схеме пластины.

Мысленно накладываем оператор центральной клеткой сначала на  узел «1» и каждой периферийной клеточкой на соответствующие соседние узлы. Пишем уравнение:

20*w₁-8*(w₀+w₂+w₃+w₂)+2*(w₀+w₀+w₁+w₁)+1*(w_-1+w₀+

w₁+w₀)=k*q*δ⁴/D

Думаю, понятно, как оно получилось? Уясните этот шаг, он ключевой при составлении уравнений!

После преобразований получим:

24*w₁-16*w₂-8*w₃=k*q*δ⁴/D

Выражение в правой части уравнения справедливо для распределенной нагрузки q.

Для сосредоточенной силы F выражение имеет вид: F*δ²/D

Если в точке действуют несколько различных нагрузок, то выражения их характеризующие суммируются в правой части уравнения.

Что такое коэффициент веса k? Если бы распределенная нагрузка была приложена только к центральной части пластины, ограниченной квадратом с узлами «2» по углам, то:

для узлов «2»: k=0,25

для узлов «1»: k=0,5

для узла «3»: k=1

Так как в нашем примере распределенная нагрузка приложена ко всей поверхности пластины, то для всех внутренних узлов k=1.

Накладывая оператор формулы метода конечных разностей последовательно на все уникальные внутренние узлы, нужно составить все уравнения. Из-за равномерности сплошной нагрузки и симметрии формы в нашем примере таких узлов 3, а, следовательно, и уравнений будет тоже 3.

Запишем эти уравнения в таблицу Excel в виде матрицы, заполнив область ячеек O5…Q7.

4. Заполним область R5…R7 значениями коэффициентов k.

5. Вычислим свободные члены уравнений в мм

в ячейках S5…S7: k*q*δ⁴/D

6. Определим коэффициенты обратной матрицы. Для этого впишем формулу массива

в область O9…Q11: {=МОБР(O5:Q7)}

Чтобы правильно ввести формулу массива напишем

в O9: =МОБР(O5:Q7)

Выделим область для обратной матрицы равную по размеру исходной матрице – O9…Q11, нажмем на клавиатуре клавишу F2, затем – Ctrl+Shift+Enter. Excel мгновенно рассчитал коэффициенты обратной матрицы.

7. Для вычисления корней уравнений системы (читай: искомых прогибов) нужно перемножить обратную матрицу с матрицей свободных членов уравнений. Для этого запишем формулу массива

в область ячеек S9…S11: {=МУМНОЖ(O9:Q11;S5:S7)}

Всё, задача решена. Используя метод конечных разностей, мы определили прогибы в узлах «1», «2» и «3». Наибольший прогиб ожидаемо оказался в центре пластины в узле «3».

Ниже приведены скриншоты решения этой же задачи при уменьшенном в 2,5 раза шаге сетки δ. Excel легко решает в этом случае систему из 15-и уравнений!

Что в итоге?

Расчет прогиба в центре пластины из рассмотренного примера по «точным» аналитическим формулам  дает следующие результаты при μ=0,3:

w_ц=0,00406*q*a⁴/D=3,632 мм (Тимошенко С.П.)

w_ц=0,0443*q*a⁴/(E*h³)=3,629 мм (Вайнберг Д.В.)

Прогибы, определенные методом конечных разностей:

При шаге сетки δ=a/4:

w_ц=w₃=3,604 мм (относительная погрешность -0,8%)

При шаге сетки δ=a/10:

w_ц=w₁₅=3,630 мм (относительная погрешность 0,0%)

Как видим, метод конечных разностей даже при большом шаге сетки дает незначительную погрешность результатов!

Используя рассмотренную методику, можно выполнять расчеты прямоугольных пластин с разными способами закрепления контуров и произвольным расположением распределенных и концентрированных нагрузок. Составление уравнений даже без навыка использования оператора формулы не займет, поверьте, у вас много времени. На вариант с 15-ю узлами было затрачено не более получаса.

Развитие программы возможно для решения других похожих задач.

Ссылка на скачивание файла с программой: metod-konechnyh-raznostej (xls 164KB).

Статьи с близкой тематикой

Отзывы

Задание.
Решить
краевую задачу методом конечных
разностей, построив два приближения
(две итерации) с шагом h и с шагом h/2.

Проанализировать
полученные результаты. Варианты заданий
приведены в приложении 4.

Порядок выполнения работы

1. Постройте
   вручную
   конечноразностную аппроксимацию
   краевой задачи (конечноразностную
   СЛАУ) с шагом h,
   заданным вариантом.

2. Используя
   метод конечных разностей, сформируйте
   в Excel
   систему
   линейных алгебраических конечно-разностных
   уравнений для шага h
   разбивки отрезка [a,
   b].
   Запишите эту СЛАУ на рабочем листе
   книги Excel.
   Расчетная схема приведена на рис.4.1.

3. Полученную
   СЛАУ решите методом прогонки.

4. Проверьте
   правильность решения СЛАУ с помощью
   надстройки Excel
   Поиск решения.

5. Уменьшите
   шаг сетки в 2 раза и еще раз решите
   задачу. Результаты представьте в
   графическом виде.

6. Сравните
   полученные результаты. Сделайте вывод
   о необходимости продолжения или о
   прекращении счета.

Решение краевой задачи с использованием электронных таблиц Microsoft Excel.

Пример
4.1. Методом
конечных разностей найти решение краевой
задачи
,
y(1)=1, y^’
(2)=0,5
на отрезке
x[1,
2] с шагом
h=0,2 и с шагом h=0,1. Сравнить полученные
результаты и сделать вывод о необходимости
продолжения или о прекращении счета.

Расчетная схема
для шага h=0,2 приведена на рис.4.1.

Полученное
решение (сеточную функцию) Y{1.000,
1.245, 1.474, 1.673, 1.829, 1.930}, Х
{1; 1,2; 1,4; 1,6; 1,8;2} в столбце L
и B
можно принять за первую итерацию (первое
приближение) исходной задачи.

Рис.4.1.

Для
нахождениявторой
итерации
сделайте сетку вдвое гуще (n=10, шаг h=0,1)
и повторите приведенный выше алгоритм.

Это
можно проделать на том же или на другом
листе книги Excel.
Решение (второе приближение) приведено
на рис.4.2.

Сравните
полученные приближенные решения. Для
наглядности можно построить графики
этих двух приближений (двух сеточных
функций), рис.4.3.

Рис.4.2.

Порядок построения графиков приближенных решений краевой задачи

1. Постройте
   график решения задачи для разностной
   сетки с шагом h=0,2
   (n=5).

2. Активизируйте
   уже построенный график и выберите
   команду меню
   ДиаграммаДобавить данные

3. В
   окне Новые
   данные
   укажите данные x_i,
   y_i
   для разностной сетки с шагом h/2 (n=10).

4. В
   окне Специальная
   вставка
   установите флажки в полях:

• новые
  ряды,

• категории(значение
  оси х) в первом столбце.

Как
видно из приведенных данных, два
приближенных решения краевой задачи
(две сеточные функции) отличаются друг
от друга не более, чем на 5%. Поэтому за
приближенное решение исходной задачи
принимаем вторую итерацию, т.е.

Y
{1, 1.124,
1.246, 1.364, 1.478, 1.584, 1.683, 1.772, 1.849, 1.914, 1.964}

Рис.4.3.

Лабораторная работа №5 Тема. Численные методы оптимизации. Графический метод

Задание.
Решить задачу
линейного программирования, приведенную
в приложении 5, графическим методом и
на ЭВМ. Сравнить полученные решения.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]

• #
• #
• #
• #
• #
• #
• #
• #
• #
• #
• #

Уважаемый посетитель!

Чтобы распечатать файл, скачайте его (в формате Word).

Ссылка на скачивание — внизу страницы.

Помогаю со студенческими работами здесь

Методом конечных разностей найти решение краевой задачи
Методом конечных разностей найти решение краевой задачи
y»-7*x*y’-y=5*x^2
u(0)-u'(0)=0; u(1)=1…

Метод конечных разностей для краевой задачи
Помогите пожалуйста. Надо написать программу, которая решает краевую задачу для ОДУ 2 порядка…

Решение краевой задачи для дифференциального уравнения методом сеток
Мне в университете задали данную программу написать. но к сожалению на паре мы не успели пройти…

Решение задачи методом конечных разностей
Доброе время суток, прошу помочь решить задачку методом конечных разностей в маткад, или хотя бы…

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму: