Автор работы: Пользователь скрыл имя, 07 Мая 2012 в 18:59, контрольная работа
МКЭ для двумерной краевой задачи для эллиптического уравнения в декартовой системе координат. Базисные функции линейные на треугольниках. Краевые условия всех типов. Коэффициент разложить по линейным базисным функциям. Матрицу СЛАУ генерировать в разреженном строчном формате. Для решения СЛАУ использовать МСГ или ЛОС с неполной факторизацией.
аналитическое решение ui |
численное решение u*i |
погрешность |
Первые краевые условия | ||
5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 |
5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000001e+000 |
0.000000000000000e+000 0.000000000000000e+000 0.000000000000000e+000 0.000000000000000e+000 0.000000000000001e+000 |
Первые и вторые краевые условия | ||
5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 |
5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000001e+000 |
0.000000000000000e+000 0.000000000000000e+000 0.000000000000000e+000 0.000000000000000e+000 0.000000000000001e+000 |
Первые и третьи краевые условия | ||
5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 |
5.000000000000000e+000 5.000000000000000e+000 5.000000000000001e+000 5.000000000000000e+000 5.000000000000001e+000 |
0.000000000000000e+000 0.000000000000000e+000 0.000000000000001e+000 0.000000000000000e+000 0.000000000000001e+000 |
Все краевые условия | ||
5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 |
5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000000e+000 5.000000000000001e+000 |
0.000000000000000e+000 0.000000000000000e+000 0.000000000000000e+000 0.000000000000000e+000 0.000000000000001e+000 |
Тест 2:
Аналитическое решение
4 конечных элементов
аналитическое решение ui |
численное решение u*i |
погрешность |
Первые краевые условия | ||
7.000000000000000e+000 1.100000000000000e+001 1.300000000000000e+001 9.000000000000000e+000 1.000000000000000e+001 |
7.000000000000000e+000 1.100000000000000e+001 1.300000000000000e+001 9.000000000000000e+000 1.000000000000000e+001 |
|
Первые и вторые краевые условия | ||
7.000000000000000e+000 1.100000000000000e+001 1.300000000000000e+001 9.000000000000000e+000 1.000000000000000e+001 |
7.000000000000000e+000 1.100000000000000e+001 1.227272727272727e+001 9.727272727272728e+000 9.999999999999998e+000 |
0.000000000000000E+00 0.000000000000000E+00 7.272727272727990E-01 7.272727272727200E-01 1.065814103640150E-14 |
Первые и третьи краевые условия | ||
7.000000000000000e+000 1.100000000000000e+001 1.300000000000000e+001 9.000000000000000e+000 1.000000000000000e+001 |
7.000000000000000e+000 1.100000000000000e+001 1.241463414634147e+001 9.585365853658535e+000 1.000000000000000e+001 |
0.000000000000000E+00 0.000000000000000E+00 5.853658536585990E-01 5.853658536585300E-01 0.000000000000000E+00 |
Все краевые условия | ||
7.000000000000000e+000 1.100000000000000e+001 1.300000000000000e+001 9.000000000000000e+000 1.000000000000000e+001 |
7.000000000000000e+000 1.100000000000000e+001 1.286391047976414e+001 9.449276333422674e+000 1.004532967032967e+001 |
0.000000000000000E+00 0.000000000000000E+00 1.360895202359010E-01 4.492763334226700E-01 4.532967032960040E-02 |
Тест 3:
Аналитическое решение
4 конечных элементов
аналитическое решение ui |
численное решение u*i |
погрешность |
Первые краевые условия | ||
2.000000000000000e+000 2.000000000000000e+000 4.000000000000000e+000 4.000000000000000e+000 3.000000000000000e+000 |
2.000000000000000e+000 2.000000000000000e+000 4.000000000000000e+000 4.000000000000000e+000 3.000000000000001e+000 |
0.000000000000000e+000 0.000000000000000e+000 0.000000000000000e+000 0.000000000000000e+000 0.000000000000001e+000 |
Первые и вторые краевые условия | ||
2.000000000000000e+000 2.000000000000000e+000 4.000000000000000e+000 4.000000000000000e+000 3.000000000000000e+000 |
2.000000000000000e+000 2.000000000000000e+000 4.000000000000001e+000 3.999999999999999e+000 3.000000000000000e+000 |
0.000000000000000E+00 0.000000000000000E+00 0.000000000000000E+00 1.021405182655140E-14 0.000000000000000E+00 |
Первые и третьи краевые условия | ||
2.000000000000000e+000 2.000000000000000e+000 4.000000000000000e+000 4.000000000000000e+000 3.000000000000000e+000 |
2.000000000000000e+000 2.813559322033900e+000 3.186440677966103e+000 4.000000000000000e+000 3.000000000000001e+000 |
0.000000000000000E+00 8.135593220339000E-01 8.135593220339000E-01 0.000000000000000E+00 0.000000000000000E+00 |
Все краевые условия | ||
2.000000000000000e+000 2.000000000000000e+000 4.000000000000000e+000 4.000000000000000e+000 3.000000000000000e+000 |
2.000000000000000e+000 2.819480226525307e+000 3.992092984601503e+000 4.062757482577259e+000 3.126547863562431e+000 |
0.000000000000000E+00 8.194802265253000E-01 7.907015398500140E-03 6.275748257724970E-02 1.265478635624300E-01 |
Тест 4:
Аналитическое решение
4 конечных элемента
аналитическое решение ui |
численное решение u*i |
погрешность |
Первые краевые условия | ||
4.000000000000000e+000 4.000000000000000e+000 1.600000000000000e+001 1.600000000000000e+001 9.000000000000000e+000
|
4.000000000000000e+000 4.000000000000000e+000 1.600000000000000e+001 1.600000000000000e+001 8.947368421052632e+000 |
0.000000000000000E+00 0.000000000000000E+00 0.000000000000000E+00 0.000000000000000E+00 5.263157894736990E-02 |
Первые и вторые краевые условия | ||
4.000000000000000e+000 4.000000000000000e+000 1.600000000000000e+001 1.600000000000000e+001 9.000000000000000e+000 |
4.000000000000000e+000 4.000000000000000e+000 1.601600000000001e+001 1.601600000000000e+001 8.952000000000002e+000 |
0.000000000000000E+00 0.000000000000000E+00 1.599999999999820E-02 1.599999999999820E-02 4.800000000000000E-02 |
Первые и третьи краевые условия | ||
4.000000000000000e+000 4.000000000000000e+000 1.600000000000000e+001 1.600000000000000e+001 9.000000000000000e+000 |
4.000000000000000e+000 7.258375680856893e+000 9.495663816450115e+000 1.600000000000000e+001 8.477558348294437e+000 |
0.000000000000000E+00 3.258375680856890E+00 6.504336183549890E+00 0.000000000000000E+00 5.224416517055700E-01 |
Все краевые условия | ||
4.000000000000000e+000 4.000000000000000e+000 1.600000000000000e+001 1.600000000000000e+001 9.000000000000000e+000 |
4.000000000000000e+000 6.472971372161895e+000 1.399297137216189e+001 1.600000000000000e+001 9.014807502467917e+000 |
0.000000000000000E+00 2.472971372161890E+00 2.007028627838200E+00 0.000000000000000E+00 1.480750246791020E-02 |