Постановка краевых задач для поперечных колебаний струны. Явная схема для численного решения одномерного уравнения колебаний

Страницы работы

Содержание работы

ЛАБОРАТОРНАЯ РАБОТА ПММ-2-2-10 ЧИСЛЕННОЕ РЕШЕНИЕ ОДНОМЕРНОГО УРАВНЕНИЯ КОЛЕБАНИЙ ПО ЯВНОЙ КОНЕЧНО-РАЗНОСТНОЙ СХЕМЕ

Цель работы: практическая реализация и исследование явной схемы решения одномерного уравнения колебаний.

ВОПРОСЫ ДЛЯ ИЗУЧЕНИЯ ПО КОНСПЕКТУ ЛЕКЦИЙ

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

2.  Различные формы граничных условий.

3.  Явная схема для численного решения одномерного уравнения колебаний.

4.  Учет начальных условий в явной схеме для уравнения колебаний.

5.  Учет граничных условий в явной схеме для уравнения колебаний.

6.  Последовательность вычислений по явной схеме.

СОДЕРЖАНИЕ ЗАДАНИЯ

Дана смешанная краевая задача для уравнения теплопроводности

Utt=a2Uxx+F(x,t)

в области 0≤x≤L0≤t≤Tmax с соответствующими начальными и граничными условиями (вариант задачи – см. в Приложении).

Выполнить следующие действия в среде Mathcad:

1)  Проверить точное решение задачи подстановкой в уравнение, начальные и граничные условия с применением символьных вычислений;

2)  Записать в текстовом блоке физический  смысл граничных условий;

3)  Составить функцию – программный блок Mathcad для численного решения краевой задачи по явной конечно-разностной схеме. 

Значения шага по координате h и шага по времени τ включить в общий список параметров.

Внутри программного блока вычисления начинать с расчета  количества шагов по координате N=L/h,  и по времени  M=Tmax/τ

4)  Получить численное решение краевой задачи при Tmax=1, h=0.05, a2=1построить график численного решения при t=Tmax

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

На основании условия устойчивости вычислить максимально допустимое значение шага по времени

τmax=h/(a2)1/2

и выбрать величину шага по времени τ=τmax (или несколько меньше - с учетом требования согласованности сетки (Tmax/τ – целое число !).

(Примечание: для некоторых задач устойчивое решение не реализуется при τ=τmax. В таких случаях следует уменьшить значении τ, не нарушая требование согласованности.

5)  С помощью отдельной функции – программного блока вычислить точное решение в тех же точках и в те же моменты времени и вычислить максимальную абсолютную погрешность численного решения;

6)  Построить общий график численного и точного решения при t=Tmax, убедиться в их согласованности;

7)  Выполнить те же действия при h=0.025 и определить, во сколько раз уменьшилась максимальная погрешность при уменьшении h.

8)  Используя результаты расчета при h=0.025,  выполнить анимацию графика для иллюстрации колебаний, связав номер столбца расчетной матрицы U (момент времени) с номером кадра FRAME.

ПРИЛОЖЕНИЕ. ВАРИАНТЫ ПОСТАНОВКИ КРАЕВОЙ ЗАДАЧИ (a2=1)

Граничные условия

Начальные условия

f(x,t)

Uточн(x,t)

1

Ux(0,t)=1;

Ux(1,t)=1;

U(x,0)=3cos(πx)+x

+2;

Ut(x,0)=0

2cos(πx)-2π2cos(πt)

3cos(πx)+2cos(πt) +x

2

U(0,t)=0; Ux(1,t)-U(1,t)=- π

U(x,0)=x2sin(πx) Ut(x,0)=3x

2x2-2)sin(πx)- -4πx cos(πx)

x2sin(πx)+3xt

3

Ux(0,t)= -1 Ux(1,t)=1

U(x,0)= x2-x Ut(x,0)=2π

-2π2 sin(πt)-2

x2-x+2sin(πt)

4

U(0,t)=tcos(πt) Ux(1,t)= - π

U(x,0)= xsin(πx) Ut(x,0)=1

π2xsin(πx)-

2πcos(πx)-

- π2tcos(πt) -2π sin(πt)

xsin(πx)+tcos(πt)

5

U(0,t)=0 Ux(1,t)-U(1,t)=t

U(x,0)=0 Ut(x,0)= x2-x

2(x-t)

xt2+t(x2-x)

6

U(0,t)= -4t Ux(1,t)-2U(1,t)=4t

U(x,0)=0 Ut(x,0)= -4x

2(x2- t2)

x2t2-4xt

7

U(0,t)=0

Ux(1,t)-U(1,t)=

= πt sin(πt)

U(x,0)=0 Ut(x,0)=x cos(πx)

2π(t-x) sin(πx+πt)

xtcos(πx+πt)

8

U(0,t)=0 Ux(1,t)-U(1,t)=t

U(x,0)=0 Ut(x,0)=x2

2(x-t)

xt2+x2t

9

U(0,t)=0

U(1,t)=0

U(x,0)= sin(πx) Ut(x,0)=0

0

sin(πx)cos(πt)

10

U(0,t)=t Ux(1,t)-U(1,t)=t

U(x,0)=0

Ut(x,0)= cos(πx)+πx

π2tcos(πx)-

π2xsin(πt)

xsin(πt)+

+tcos(πx)

11

U(0,t)=0 Ux(1,t)-U(1,t)=-πt

U(x,0)=x Ut(x,0)= sin(πx)

-4π2xcos(2πt)+ +π2tsin(πx)

xcos(2πt)+ +tsin(πx)

12

U(0,t)=t Ux(1,t)-U(1,t)=t

U(x,0)=x Ut(x,0)= cos(πx)

π2tcos(πx)-

π2xcos(πt)

xcos(πt)+tcos(πx)

13

U(0,t)=0; Ux(1,t)-U(1,t)=A

U(x,0)=x2sin(2x) Ut(x,0)=3x

(4x2-2)sin(2x)- -8x cos(2x)

x2sin(2x)+3xt

14

Ux(0,t)=1;

Ux(1,t)=1;

U(x,0)=3cos(πx)+x

+2;

Ut(x,0)=0

2cos(πx)-2π2cos(πt)

3cos(πx)+2cos(πt) +x

15

U(0,t)=0; Ux(1,t)-U(1,t)=- π

U(x,0)=x2sin(πx) Ut(x,0)=3x

2x2-2)sin(πx)- -4πx cos(πx)

x2sin(πx)+3xt

16

U(0,t)=tcos(πt) Ux(1,t)= - π

U(x,0)= xsin(πx) Ut(x,0)=1

π2xsin(πx)-

2πcos(πx)-

- π2tcos(πt) -2π sin(πt)

xsin(πx)+tcos(πt)

17

U(0,t)=0 Ux(1,t)-U(1,t)=t

U(x,0)=0 Ut(x,0)= x2-x

2(x-t)

xt2+t(x2-x)

18

U(0,t)= -4t Ux(1,t)-2U(1,t)=4t

U(x,0)=0 Ut(x,0)= -4x

2(x2- t2)

x2t2-4xt

19

U(0,t)=0

Ux(1,t)-U(1,t)=

= πt sin(πt)

U(x,0)=0 Ut(x,0)=x cos(πx)

2π(t-x) sin(πx+πt)

xtcos(πx+πt)

20

U(0,t)=0 Ux(1,t)-U(1,t)=t

U(x,0)=0 Ut(x,0)=x2

2(x-t)

xt2+x2t

21

U(0,t)=0

U(1,t)=0

U(x,0)= sin(πx) Ut(x,0)=0

0

sin(πx)cos(πt)

22

U(0,t)=t Ux(1,t)-U(1,t)=t

U(x,0)=0

Ut(x,0)= cos(πx)+πx

π2tcos(πx)-

π2xsin(πt)

xsin(πt)+

+tcos(πx)

23

U(0,t)=0 Ux(1,t)-U(1,t)=-πt

U(x,0)=x Ut(x,0)= sin(πx)

-4π2xcos(2πt)+ +π2tsin(πx)

xcos(2πt)+ +tsin(πx)

24

U(0,t)=t Ux(1,t)-U(1,t)=t

U(x,0)=x Ut(x,0)= cos(πx)

π2tcos(πx)-

π2xcos(πt)

xcos(πt)+tcos(πx)

25

U(0,t)=0; Ux(1,t)-U(1,t)=A

U(x,0)=x2sin(2x) Ut(x,0)=3x

(4x2-2)sin(2x)- -8x cos(2x)

x2sin(2x)+3xt

СОДЕРЖАНИЕ ОТЧЕТА

Отчет по каждой лабораторной работе должен содержать распечатку МСдокумента, в котором: 1) указаны тема работы, фамилия и номера варианта исполнителя; 2) для каждой задачи приведено общее условие и условие конкретного варианта; 3) выполнены необходимые преобразования; 4) реализованы вычисления, снабженные текстовыми пояснениями.

Похожие материалы

Информация о работе