Basilio
|
|
GreenOne
|
|
|
|
|
|
|
Рег.: 14.10.2002
|
|
Сообщений: 13750
|
|
Из: Москва
|
|
Рейтинг: 3476
|
|
УРЧП/МКЭ: Граничные условия: не Неймана и не Дирихле
19.11.2006 21:54
|
|
|
суть задачи в двух словах следующая:
решается уравнение Lu=P на области (квадрат 1х1), с краевым условием u|Г=0 (Г- периметр квадрата).
Дополнительно нужно решить задачу с "условием периодичности" - то же уравнение Lu=P, но краевое условие такое: u|Г1=u|Г3 и u|Г2=u|Г4 (Г1, Г2, Г3, Г4 - соответствующие стороны квадрата).
(с тем, чтобы сравнить два полученных решения и сделать некоторые научные выводы)
с первым случаем все просто - конструирование функционала для минимизации и соответствующей системы линейных уравнений описано в десятках книг.
Непонятно, что делать во втором случае. Подскажите, где искать.
|
|
|
ABC47
|
|
русский
|
|
|
|
|
|
|
Рег.: 26.02.2005
|
|
Сообщений: 28318
|
|
Из: Волгоградской области
|
|
Рейтинг: -676
|
|
Re: УРЧП/МКЭ: Граничные условия: не Неймана и не Дирихле
[re: Basilio]
19.11.2006 22:07
|
|
|
Решай методом Фурье 
|
Убить жида, чтобы купить пистолет,
Убить жида, чтобы ты был вооружен.
(С) ИПВ |
|
|
Uranius
|
|
member
|
|
|
|
|
|
|
Рег.: 18.08.2003
|
|
Сообщений: 102
|
|
|
|
Рейтинг: 1
|
|
Re: УРЧП/МКЭ: Граничные условия: не Неймана и не Дирихле
[re: Basilio]
20.11.2006 00:18
|
|
|
Решал и с краевыми условиями первого типа и с периодическими, разницы не вижу никакой.
Два замечания:
1) В областях типа прямоугольник всегда нужно использовать сеточные методы при известной гладкости решения, несмотря на то что МКЭ - это круто.
2) Уравнения в МКЭ не обязательно получать минимизацией функционала (при эллиптичном L), а домножаем уравнение на базисные функции и интегрируем по области (метод Галеркина), при этом используется теорема интегрирования по частям, где периодические краевые условия срабатывают также хорошо, как и нулевые условия - интегралы по границе области уничтожаются!
3) Метод Фурье работает и в периодическом случае, если известен спектр оператора L с периодическими условиями.
При способе 2 дискретные уравнения получатся скорее точно такими, как и при минимизации функционала.
|
|
|
Basilio
|
|
GreenOne
|
|
|
|
|
|
|
Рег.: 14.10.2002
|
|
Сообщений: 13750
|
|
Из: Москва
|
|
Рейтинг: 3476
|
|
Re: УРЧП/МКЭ: Граничные условия: не Неймана и не Дирихле
[re: Uranius]
20.11.2006 13:12
|
|
|
сеточные методы не очень хороши в виду наличия включений и дырок в прямоугольной области довольно сложной геометрии
про Галеркина - видел, но пролистал. поботаю, спасибо.
|
|
|
Gonobobel
|
|
|
|
|
|
|
|
|
Рег.: 20.05.2006
|
|
Сообщений: 10715
|
|
|
|
Рейтинг: 4318
|
|
Re: УРЧП/МКЭ: Граничные условия: не Неймана и не Дирихле
[re: Basilio]
20.11.2006 15:48
|
|
|
Quote:
сеточные методы не очень хороши в виду наличия включений и дырок в прямоугольной области довольно сложной геометрии
сетку можно брать кривую  правда, тогда вычислительные формулы будут сложнее. Или сделать замену координат (это не то же самое, что криволинейная сетка).
Или разбить область на несколько мелких областей более простой геометрии.
|
I have retired this character... 06.05.2010. |
|
|
Uranius
|
|
member
|
|
|
|
|
|
|
Рег.: 18.08.2003
|
|
Сообщений: 102
|
|
|
|
Рейтинг: 1
|
|
Re: УРЧП/МКЭ: Граничные условия: не Неймана и не Дирихле
[re: Basilio]
20.11.2006 16:59
|
|
|
Тогда МКЭ, Галеркин. Минимизацию используют для легкого обоснования существования решения, ее не обязательно использовать в численных методах.
Может сложность вызывает то, что, например, некоторые базисные функции будут иметь часть носителя на верхней границе, часть - на нижней? Так это обычная ситуация. Для лучшего представления можно склеить область по границам, получить тор и на нем построить конечные элементы.
|
|
|
Basilio
|
|
GreenOne
|
|
|
|
|
|
|
Рег.: 14.10.2002
|
|
Сообщений: 13750
|
|
Из: Москва
|
|
Рейтинг: 3476
|
|
Re: УРЧП/МКЭ: Граничные условия: не Неймана и не Дирихле
[re: Uranius]
20.11.2006 19:56
|
|
|
заторить - это мысль, однако, возможно, тут надо будет кое-что в мешере править
|
|
|
Uranius
|
|
member
|
|
|
|
|
|
|
Рег.: 18.08.2003
|
|
Сообщений: 102
|
|
|
|
Рейтинг: 1
|
|
Re: УРЧП/МКЭ: Граничные условия: не Неймана и не Дирихле
[re: Basilio]
20.11.2006 22:27
|
|
|
Да, сетка здесь произвольная не пойдет, нужно чтоб узлы на противоположных границах были зеркальным отражением друг друга относительно осей симметрии
|
|
Список
Плоский
Дерево
