FEM: особливість матриці жорсткості


11

(σ2(x)u(x))=f(x),0x1
u(0)=u(1)=0u(0)=u(1)=0σ(x)σ0>0Au=fA

Дотримуючись схему FEM, я звожу свою проблему до задачі оптимізації

J(u)=(Au,u)2(f,u)minu
Вводяться кінцеві елементи hk(x) як
vk(x)={1(xxkh)2,x[xk1,xk+1]0,otherwise
для будь-якого k=1,,n1 , де xk=hk , h=1n . Кінцеві елементи v0(x) і vn(x) вводяться аналогічно.

Я намагаюся знайти числовий вектор α таким, щоб u(x)=k=0nαkvk(x) вирішував задачу оптимізації. У нас

J(u)=i=0nj=0nαiαj(Avi,vj)i=0n2αi(vi,f)=αTVα2αTbminα,
де bi=(f,vi) і Vi,j=(Avi,vj) . Після диференціації відносно α я отримую
Vα=b,
але тут матриця жорсткості V є сингулярною. То що мені робити? Можливо, мені доведеться вибирати інші кінцеві елементи?

Привіт, Німза, у тебе проблема з тестом, що ти знаєш точне рішення? Якщо так, спробуйте вирішити VTVα=VTb спочатку, щоб перевірити, чи є ваша основа всередині домену, якщо все виглядає правильно, то, можливо, саме неправильно поставлений BC робить матрицю єдиною. Але БК здається мені нормальним.
Shuhao Cao

Відповіді:


13

У порядку зменшення ймовірності

  1. Неправильна основа. З вашого опису виходить, що у вас є рівно дві квадратичні функції з підтримкою кожного елемента. Цей простір не є розділом єдності і не є (безперервні перші похідні). Щоб дискретизувати свою проблему четвертого порядку безпосередньо (замість того, щоб зменшити її до системи рівнянь другого порядку, наприклад), вам знадобиться база . Зауважимо, що база повинна мати можливість точно відтворювати всі лінійні функції.C1C1C1

  2. Недостатні граничні умови. Це буде очевидно очевидно, якщо обчислити та побудувати нульовий простір.

  3. Неправильна збірка. Перевірте карту від елементів до зібраного замовлення, щоб підтвердити, що це те, що ви очікували, наприклад, що воно не змінює орієнтацію елементів.

  4. Неправильне місцеве зібрання. У 1D ви можете аналітично обчислити, як виглядає матриця жорсткості елемента (можливо, для спрощеного випадку) і перевірити, чи відтворює його код.


Дякую. 1. Я думаю, що мені знадобиться основа тому що . Тоді, якщо я розглядаю лише функції, які відповідають граничним умовам, то . C2(Au,v)=01σ2(x)u(x)v(x)dxkerA={0}
Аплікація

1
базис досить, підінтегральна не повинна бути безперервною. Зауважимо, що граничні умови для другої похідної стануть граничним інтегралом. Ви можете використовувати основу для прямої дискретизації проблеми четвертого порядку, але вам потрібно буде інтегрувати терміни стрибка, як і при розривних методах Галеркіна для систем першого та другого порядку. Це не поганий метод, але він є надмірно складним в 1D, оскільки так легко створювати бази з будь-яким порядком безперервності (наприклад, сплайни). Ця стаття є прикладом " DG". C1C0C0
Джед Браун

Гаразд. Я виправив свою основу: тепер на і . Тепер це . Але метод все ще не працює. vk(x)=cos2(π2h(xxi))[xi1,xi+1]i=1,,n1C1
Аплікація

Основа повинна бути здатна відтворювати лінійні функції, але це не може. Після виправлення цього перевірте, чи інтеграли виконуються правильно, а потім перевірте граничні умови. C1
Джед Браун

0

Очевидно, що проблема має похідне замовлення ODD. Більш конкретно для більших чисел Пеклета , матриця жорсткості може не підтримувати «тонкої» форми, що створює нулі під час складання, і, отже, отримує сингулярний або іноді дуже малий визначник, який помітний коливаннями в графіку розчину.

Вирішення подібної проблеми - серед інших методів застосування штрафу. Більш конкретно це називається методом Петрова-Галеркіна .

Вибачте за моє погане розуміння англійської мови.

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.