Уравнение у нас a2(∂2u/∂x2+∂2u/∂y2)=0 с граничными условиями na2(∂u/∂n)=0. Оно эллиптическое (щелкнем переключатель Elliptic).
Пусть коэффициент температуропроводности а2=225, на верхней и нижней границе области температура равна 580, на окружностях (краях отверстий) 450. Введем коэффициент температуропроводности в соответствующую строку, укажем также, что нет источников тепла, нет конвекции.
Определим граничные условия. Для этого щелкнем пункт меню Boundary-Boundary Mode.
Появится возможность выделить участок границы. Если необходимо выделить несколько участков границы, то щелчок по каждому из них осуществляется при нажатой кнопке Shift. Выделим верхнюю и нижнюю границы области.
Для определения условий идем в пункт меню Boundary-Specify Boundary Conditions. В окне Boundary Conditions выберем условия Дирихле. Введем данные. Пусть температура равна 580. Затем выделим границы отверстий и установим для них температуру 450.
После этого для левой и правой границ введем условия Неймана (Рис.1.45).
Рис.1.45. Для левой и правой границ введем условия Неймана.
|
|
Теперь надо провести триангуляцию области (покрыть ее сеткой из треугольников). Для этого нажмем пункт меню Mesh-Mesh Mode и получим крупные треугольники и деформированные отверстия (рис.1.46). Для измельчения треугольников нажмем Mesh-Refine Mesh и так пару раз, пока края отверстий не станут круглыми. Чрезмерное измельчение сетки удлиняет время расчета.
Рис.1.46. Триангуляция области.
Для получения цветного графика (Рис.1.47) надо нажать Solve-PDE Solve.
Рис.1.47. Для получения цветного графика нажать Solve-PDE Solve.
Для получения трехмерного графика (Рис.1.48) надо еще нажать Plot-Parameters и включить в открывшемся окне флажок Height (3-D Plot).